8点DFT68点DFT[x4(n)]8点DFT[x5(n)]54幅度32101234ω/π16点DFT5678 616点DFT[x4(n)]16点DFT[x5(n)]54幅度32101234ω/π5678 3.对模拟周期信号进行频谱分析 Fs=64;T=1/Fs; N=16;n=0:N-1; ?T的变换区间N=16 x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); ?T的变换区间N=16? X6k16=fft(x6nT); %计算x6nT的16点DFT X6k16=fftshift(X6k16); %将零频率移到频谱中心 Tp=N*T;F=1/Tp; %将零频率移到频谱中心 k=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心）???? stem(fk,abs(X6k16),'g');box on;hold on; %绘制8点DFT的幅频特性图???????? N=32;n=0:N-1; ?T的变换区间N=16
x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); %对x6(t)32点采样 X6k32=fft(x6nT); ?T的变换区间N=16 X6k32=fftshift(X6k32); %将零频率移到频谱中心 ???????? Tp=N*T;F=1/Tp; %频率分辨率F k=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心）???????????????????????????????? stem(fk,abs(X6k32),'r');box on;hold off;%绘制8点DFT的幅频特性图?????????????? title('|DFT[x_6(nT)]|');xlabel('f(Hz)');ylabel('幅度'); legend('16点DFT[x6(nT)]','32点DFT[x6(nT)]') figure N=64;n=0:N-1; ?T的变换区间N=16 x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); %对x6(t)64点采样 X6k64=fft(x6nT); %计算x6nT的64点DFT X6k64=fftshift(X6k64); %将零频率移到频谱中心???????????????????? Tp=N*T;F=1/Tp; %频率分辨率F k=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心）?? stem(fk,abs(X6k64),'.'); box on%绘制8点DFT的幅频特性图?????? title('(6a) 64点|DFT[x_6(nT)]|');xlabel('f(Hz)');ylabel('幅度'); 运行结果： |DFT[x6(nT)]|1614121016点DFT[x6(nT)]32点DFT[x6(nT)]幅度86420-40-30-20-10f(Hz)0102030
(6a) 64点|DFT[x6(nT)]|35302520幅度151050-40-30-20-100f(Hz)10203040