资源描述：
《离散傅里叶变换(DFT)快速傅里叶变换(FFT)-实验四.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、课程名称：数字信号处理实验成绩：指导教师：实验报告院系：信息工程学院班级：学号：姓名：日期：2011.10.30实验四：DFS、DFT与FFT1、已知某周期序列的主值序列为x(n)=[0,1,2,3,2,1,0]，编程显示2个周期的序列波形。要求：①用傅里叶级数求信号的幅度谱和相位谱，并画出图形②求傅里叶级数逆变换的图形，并与原序列进行比较程序清单：N=7;xn=[0,1,2,3,2,1,0];xn=[xn,xn];n=0:2*N-1;k=0:2*N-1;Xk=xn*exp(-1i*2*pi/N).^(n'*k);x=(Xk*exp(1i*2*pi
2、/N).^(n'*k))/N;subplot(2,2,1);stem(n,xn);title('x(n)');axis([-1,2*N,1.1*min(xn),1.1*max(xn)]);subplot(2,2,2);stem(n,abs(x));title('IDFS
3、X(k)
4、');axis([-1,2*N,1.1*min(x),1.1*max(x)]);subplot(2,2,3),stem(k,abs(Xk));title('
5、X(k)
6、');axis([-1,2*N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]);
7、subplot(2,2,4),stem(k,angle(Xk));title('arg
8、X(k)
9、');axis([-1,2*N,1.1*min(angle(Xk)),1.1*max(angle(Xk)]);程序运行结果如下图：2、已知有限长序列x(n)=[1,0.5,0,0.5,1,1,0.5,0]，要求：①求该序列的DFT、IDFT的图形；程序清单：xn=[1,0.5,0,0.5,1,1,0.5,0];N=length(xn);n=0:N-1;k=0:N-1;Xk=xn*exp(-1i*2*pi/N).^(n'*k);x=(Xk*exp(1i*
10、2*pi/N).^(n'*k))/N;subplot(2,2,1);stem(n,xn);title('x(n)');axis([-1,N,1.1*min(xn),1.1*max(xn)]);subplot(2,2,2);stem(n,abs(x));title('IDFT
11、X(k)
12、');axis([-1,N,1.1*min(x),1.1*max(x)]);subplot(2,2,3),stem(k,abs(Xk));title('
13、X(k)
14、');axis([-1,N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]);su
15、bplot(2,2,4),stem(k,angle(Xk));title('arg
16、X(k)
17、');axis([-1,N,1.1*min(angle(Xk)),1.1*max(angle(Xk))]);程序运行结果如下图：②用FFT算法求该序列的DFT、IDFT的图形；程序清单：xn=[1,0.5,0,0.5,1,1,0.5,0];N=length(xn);subplot(2,2,1);stem(n,xn);title('x(n)');k=0:N-1;Xk=fft(xn,N);subplot(2,1,2);stem(k,abs(Xk));title
18、('Xk=DFT(xn)');xn1=ifft(Xk,N);subplot(2,2,2);stem(n,xn1);title('x(n)=IDFT(Xk)程序运行结果如下图：③假定采用频率Fs=20Hz，序列长度N分别取8、32和64，用FFT计算其幅度谱和相位谱。程序清单：Ts=0.05;C=[8,32,64];forr=0:2;N=C(r+1);n=0:N-1;xn=[1,0.5,0,0.5,1,1,0.5,0,zeros(1,N-8)];D=2*pi/(N*Ts);k=floor(-(N-1)/2:(N-1)/2);X=fftshift(ff
19、t(xn,N));subplot(3,2,2*r+1);plot(k*D,abs(X));subplot(3,2,2*r+2);stairs(k*D,angle(X));end程序运行结果如下图：3、已知一个无限长序列x(n)=0.5n（n≥0），采样周期Ts=0.2s，要求序列长度分别取8、32和64，用FFT求其频谱。程序清单：Ts=0.2;C=[8,32,64];forr=0:2;N=C(r+1);n=0:N-1;xn=0.5.^n;D=2*pi/(N*Ts);k=floor(-(N-1)/2:(N-1)/2);X=fftshift(fft(
20、xn,N));subplot(3,2,2*r+1);plot(k*D,abs(X));axis([-80,80,0,3])
