实验五 离散傅立叶变换DFT
N=20; m=10; n=0:1:N-1; x=8* (0.4).^n;
n1=mod((n+m),N); %求余
xm=x(n1+1); % 求余后加1是因为MATLAB向量下标从1开始 subplot(2,1,1) stem(n,x); title(‘原序列’); xlabel(‘n’); ylabel(‘x(n)’); subplot(2,1,2) stem(n,xm);
title(‘圆周移位序列’); xlabel(‘n’); ylabel(‘xm(n)’);
3.构造一个圆周卷积函数,专门用来计算两个序列(假设均从0点开始)的循环卷积。
Function[y]=circonv(x1,x2) xn2=[x2(1),fliplr(x2)]; xn2(length(xn2))=[]; C=xn2;R=x2; M=toeplitz(C,R); y=x1* M; 如果
x1(n)??(n)?2?(n?1)?3?(n?2)?4?(n?3),
x2(n)?4?(n)?3?(n?1)?2?(n?2)??(n?3),利用上面的圆周卷积函数求
其4点圆周卷积,并显示其结果。再利用时域圆周卷积定理,调用前面的离散傅立叶正、反变换函数计算它们的4点圆周卷积,并显示其结果。
21
实验五 离散傅立叶变换DFT
4. 调用离散傅立叶正变换函数运行下面MATLAB程序,验证DFT的共轭对
称性。
x=[1 2 4 2 6 32 6 4 2 zeros(1,247)];
x1=[x(1) x(256:-1:2)];
xep=0.5*(x+x1); xop=0.5*(x-x1); XF=dft(x,256); XEPF=dft(xep,256); XOPF=dft(xop,256); clf k=0:255; subplot(2,2,1);
plot(k/128,real(XF));grid; ylabel(‘振幅’);
title(‘原序列DFT的实部’);
subplot(2,2,2);
plot(k/128,imag(XF));grid; ylabel(‘振幅’);
title(‘原序列DFT的虚部’); subplot(2,2,3);
plot(k/128,real(XEPF));grid; ylabel(‘振幅’);
title(‘其共轭对称序列DFT’);
subplot(2,2,4);
plot(k/128,imag(XOPF));grid; ylabel(‘振幅’);
title(‘其共轭反对称序列DFT’);
四、实验仪器设备
计算机,MATLAB软件
五、实验注意事项
课前预先阅读并理解实验程序;
22
实验五 离散傅立叶变换DFT
六、思考题
1.讨论实验程序1中的.^代表什么含义?
2.讨论实验程序2中圆周移位与线性移位的关系。指令mod的含义是什么? 3 .计算实验程序3的理论计算结果,对比两种程序运行结果,结果是否一致? 4 .讨论实验程序4说明实序列的共轭对称部分的离散傅立叶变换与原序列的傅立叶变换实部的关系,共轭反对称部分的离散傅立叶变换与原序列的傅立叶变换虚部的关系。
23
实验六 快速傅立叶变换FFT及其应用
实验六 快速傅立叶变换FFT及其应用
一、实验目的
1. 利用MATLAB的快速傅立叶变换来计算信号的离散傅立叶变换。 2. 利用MATLAB程序,理解进一步离散傅立叶变换的物理意义。 3. 利用MATLAB程序,理解快速卷积算法。
二、实验原理
在MATLAB中,使用函数fft可以很容易地计算有限长序列x(n)的离散傅立叶变换X[k]。此函数有两种形式,fft(x)计算序列x(n) 的离散傅立叶变换X(k),这里X(k)的长度与x(n)的长度相等。fft(x,L)计算序列x(n) 的L点离散傅立叶变换,其中L≥N。若L>N,在计算离散傅立叶变换之前,对x(n)尾部的L-N个值进行补零。同样,离散傅立叶变换序列X(k)的离散傅立叶逆变换x(n)用函数ifft计算,它也有两种形式。 (一)、基本序列的离散傅立叶变换计算
N点离散傅立叶变换的一种物理解释就是,X[k]是x(n)以N为周期的周期延拓序列的离散傅立叶级数系数X(k)的主值区间序列,即X(k)?X(k)RN(k)。例如序列
~~cos(cos(??n)RN(n),当N=16时,cos(n)RN(n)正好是cos(n)的一个周期,所以
8888n)RN(n)的周期延拓序列就是这种单一频率的正弦序列。而当N=8时,
??cos(n)RN(n)正好是cos(n)的半个周期,cos(n)RN(n)的周期延拓就不再是单一频888率的正弦序列,而是含有丰富的谐波成分,其离散傅立叶级数的系数与N=16时的差别很大,因此对信号进行谱分析时,一定要截取整个周期,否则得到错误的频谱。
???
(二)、验证N点DFT的物理意义
假如x(n)非周期、有限长,则傅立叶变换存在,那么对X(ej?)在N个等间隔频率
?k=2πk/N, k=0,1,…,,N-1取样,则可得X(k)。
24