资源描述：
《数字信号处理报告 序列的傅里叶变换》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、信息科学与技术学院本科三年级数字信号处理实验报告专业：自动化班级：自动化一班姓名:陈银银学号：20092484指导老师：张翠芳2011年12月21日实验一序列的傅立叶变换实验目的进一步加深理解DFS,DFT算法的原理;1.研究补零问题;2.快速傅立叶变换（FFT）的应用。实验步骤1.复习DFS和DFT的定义，性质和应用；2.熟悉MATLAB语言的命令窗口、编程窗口和图形窗口的使用；3.利用提供的程序例子编写实验用程序；4.按实验内容上机实验，并进行实验结果分析；5.写出完整的实验报告，并将程序附在后面。实验
2、内容1.周期方波序列的频谱已知试画出下面四种情况下的的幅度频谱,并分析补零后，对信号频谱的影响。实验结果：1.有限长序列x(n)的DFT已知：（1）取x(n)(n=0:10)时，画出x(n)的频谱X(k)的幅度；（2）将（1）中的x(n)以补零的方式，使x(n)加长到(n:0~100)时，画出x(n)的频谱X(k)的幅度；（3）（3）取x(n)(n:0~100)时，画出x(n)的频谱X(k)的幅度。2.利用FFT进行谱分析已知：模拟信号以t=0.01n(n=0:N-1)进行采样，求N点DFT的幅值谱。请分别
3、画出N=45;N=50;N=55;N=60时的幅值曲线。实验结果：X（n）补零后：以t=0.01n(n=0:N-1)进行采样，N点DFT的幅值谱：实验代码：%内容1%L=5;N=20;L=5;N=20;n=1:N;xn=[ones(1,L),zeros(1,N-L)];Xk=dfs(xn,N);magXk=abs([Xk(N/2+1:N)Xk(1:N/2+1)]);k=[-N/2:N/2];figure(1)subplot(2,1,1);stem(n,xn);xlabel('n');ylabel('xtid
4、e(n)');title('DFSofSQ.wave:L=5,N=20');subplot(2,1,2);stem(k,magXk);axis([-N/2,N/2,0,16]);xlabel('k');ylabel('Xtide(k)');%L=5;N=40;L=5;N=40;n=1:N;xn=[ones(1,L),zeros(1,N-L)];Xk=dfs(xn,N);magXk=abs([Xk(N/2+1:N)Xk(1:N/2+1)]);k=[-N/2:N/2];figure(2)subplot(2,1,
5、1);stem(n,xn);xlabel('n');ylabel('xtide(n)');title('DFSofSQ.wave:L=5,N=40');subplot(2,1,2);stem(k,magXk);axis([-N/2,N/2,0,16]);xlabel('k');ylabel('Xtide(k)');%L=5;N=60;L=5;N=60;n=1:N;xn=[ones(1,L),zeros(1,N-L)];Xk=dfs(xn,N);magXk=abs([Xk(N/2+1:N)Xk(1:N/2+1
6、)]);k=[-N/2:N/2];figure(3)subplot(2,1,1);stem(n,xn);xlabel('n');ylabel('xtide(n)');title('DFSofSQ.wave:L=5,N=60');subplot(2,1,2);stem(k,magXk);axis([-N/2,N/2,0,16]);xlabel('k');ylabel('Xtide(k)');%L=5;N=20;L=7;N=60;n=1:N;xn=[ones(1,L),zeros(1,N-L)];Xk=dfs(
7、xn,N);magXk=abs([Xk(N/2+1:N)Xk(1:N/2+1)]);k=[-N/2:N/2];figure(4)subplot(2,1,1);stem(n,xn);xlabel('n');ylabel('xtide(n)');title('DFSofSQ.wave:L=7,N=60');subplot(2,1,2);stem(k,magXk);axis([-N/2,N/2,0,16]);xlabel('k');ylabel('Xtide(k)');%内容2xn=cos(0.48*pi*n)+
8、cos(0.52*pi*n);%²»²¹ÁãM=10;N=10;n=1:M;n1=[0:1:N-1];y1=[xn(1:1:M),zeros(1,N-M)];figure(1)subplot(2,1,1);stem(n1,y1);xlabel('n');title('signalx(n),0<=n<=100');axis([0,N,-2.5,2.5]);Y1=fft(y1);magY1=abs(Y1(1:1:N/
