快速傅立叶变换和其应用实验.docVIP

附录：实验1 快速傅里叶变换及其应用 1、程序一：n=0:1:15;%p=8不变，q变化（2,4,8）;p=8;q=2; %p=8;q=2; xa1=exp(-((n-p).^2)/q);subplot(5,2,1);plot(n,xa1,-*);xlabel(t/T);ylabel(xa(n));title(p=8 q=2)xk1=abs(fft(xa1));subplot(5,2,2);stem(n,xk1)xlabel(k);ylabel(Xa(k));title(p=8 q=2)p=8;q=4; %p=8;q=4; xa1=exp(-((n-p).^2)/q);subplot(5,2,3);plot(n,xa1,-*);xlabel(t/T);ylabel(xa(n));title(p=8 q=4)xk1=abs(fft(xa1));subplot(5,2,4);stem(n,xk1)xlabel(k);ylabel(Xa(k));title(p=8 q=4)p=8;q=8; %p=8;q=8; xa1=exp(-((n-p).^2)/q);subplot(5,2,5);plot(n,xa1,-*);xlabel(t/T);ylabel(xa(n));xk1=abs(fft(xa1));title(p=8 q=8)subplot(5,2,6);stem(n,xk1)xlabel(k);ylabel(Xa(k));title(p=8 q=8)%q=8不变,p变化(8,13,14);p=8;q=8; %p=8;q=8; xa1=exp(-((n-p).^2)/q);subplot(5,2,5);plot(n,xa1,-*);xlabel(t/T);ylabel(xa(n));xk1=abs(fft(xa1));title(p=8 q=8)subplot(5,2,6);stem(n,xk1)xlabel(k);ylabel(Xa(k));title(p=8 q=8)p=13;q=8; %p=13;q=8; xa1=exp(-((n-p).^2)/q);subplot(5,2,7);plot(n,xa1,-*);xlabel(t/T);ylabel(xa(n));xk1=abs(fft(xa1));title(p=13 q=8)subplot(5,2,8);stem(n,xk1)xlabel(k);ylabel(Xa(k));title(p=13 q=8)p=14;q=8; %p=14;q=8; xa1=exp(-((n-p).^2)/q);subplot(5,2,9);plot(n,xa1,-*);xlabel(t/T);ylabel(xa(n));title(p=14 q=8)xk1=abs(fft(xa1));subplot(5,2,10);stem(n,xk1)xlabel(k);ylabel(Xa(k));title(p=14 q=8)图1-1 （a）高斯序列的时域和频域特性（p=8，改变q）图1-1（b） 高斯序列的时域和频域特性（q =8，改变p）2、程序二：n1=0:1:15;xb1=exp(-0.1*n1).*sin(2*pi*0.0625*n1);subplot(3,2,1);plot(n1,xb1,-*);xlabel(n);ylabel(x(n));title(f=0.0625);xk1=abs(fft(xb1));subplot(3,2,2);stem(n1,xk1)xlabel(k);ylabel(X(k));title(f=0.0625);n2=0:1:15;xb2=exp(-0.1*n2).*sin(2*pi*0.4375*n2);subplot(3,2,3);plot(n2,xb2,-*);xlabel(n);ylabel(x(n));title(f=0.4375);xk2=abs(fft(xb2));subplot(3,2,4);stem(n2,xk2)xlabel(k);ylabel(X(k));title(f=0.4375);n3=0:1:15;xb3=exp(-0.1*n3).*sin(2*pi*0.5625*n3);subplot(3,2,5);plot(n3,xb3,-*);xlabel(n);ylabel(x(n));title(f=0.5625);xk3=abs(fft(xb3));subplot(3,2,6);stem(n3,xk3)xlabel(k);ylabel(X(k));title(f=0.5625)图1-2 衰减正弦序列的时域波形和幅频特性3、程序三：%N=8程序：n1=0:1:7;xc1=[0 1 2 3 4 3 2 1];subplo