快速傅立叶变换及其应用实验.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)