华工数字信号处理实验一.docVIP

内容(1) 代码 N1=16; n=1:200; window1=[ones(1,N1),zeros(1,200-N1)]; window2=[hamming(N1),zeros(1,200-N1)]; x=cos(pi*n/10)+sin(pi*n/6)+cos(2*pi*n/5); X1=fft(x.*window1); subplot(2,6,1); stem(abs(X1)); title(range); subplot(2,6,2); stem(angle(X1)); title(phase); X2=fft(x.*window2); subplot(2,6,3); stem(abs(X2)); title(range); subplot(2,6,4); stem(angle(X2)); title(phase); N2=64; window3=[ones(1,N2),zeros(1,200-N2)]; window4=[hamming(N2),zeros(1,200-N2)]; X3=fft(x.*window3); subplot(2,6,5); stem(abs(X3)); title(range); subplot(2,6,6); stem(angle(X3)); title(phase); X4=fft(x.*window4); subplot(2,6,7); stem(abs(X4)); title(range); subplot(2,6,8); stem(angle(X4)); title(phase); N3=128; window5=[ones(1,N3),zeros(1,200-N3)]; window6=[hamming(N3),zeros(1,200-N3)]; X5=fft(x.*window5); subplot(2,6,9); stem(abs(X5)); title(range); subplot(2,6,10); stem(angle(X5)); title(phase); X6=fft(x.*window6); subplot(2,6,11); stem(abs(X6)); title(range); subplot(2,6,12); stem(angle(X6)); title(phase); 内容（2） 1.函数最高频率f=200pi/2pi=100Hz，采样频率最低：2f=200Hz，采样最大间隔1/200=0.005s=5ms，频域采样最小点数：N=200Hz/1=200个点。 2.代码： t=0:0.005:1; x=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); subplot(2,1,1); plot(t,x); n=input(序列长度); X1=fft(x,n); subplot(2,1,1); stem(abs(X1)); title(DFT_range); xlabel(k); ylabel(range); subplot(2,1,2); stem(angle(X1)); title(DFT_phase); xlabel(k); ylabel(phase); 3.以200Hz采样后图像为： n=200 4.当n减小时，取n=150： 当n增大时，取n=250 从这两幅图可以看出，不论M大于200还是小于200，都会存在误差；当M大于200时，单位脉冲向右偏移；当M小于200时向左偏移； 内容（3） 1）取 0≤n ≤9，计算其DFT； 代码： N=10; M=100; n=0:N-1; xn=cos(0.48*pi*n)+cos(0.52*pi*n); Y=fft(xn,M); figure(1); stem(n,xn); title(x[n]); xlabel(n); ylabel(range); figure(2); k=0:M-1; stem(k,abs(Y)); title(DFT_range); xlabel(k); ylabel(range); figure(3) stem(k,angle(Y)); title(DFT_phase); xlabel(k); ylabel(phase); 2）将1）中的 x[n]补零，使 0≤n ≤99，计算其DFT； 3）增加采样的个数，即对原x[n]取0≤n ≤99，计算其DFT。 只需要就1）中N=10改成N=100即可得到想要的结果： 三个步骤得到的频谱都不相同，第二个步骤的频谱包络是第一个信号的频谱。