【正文】 樣頻率及能夠分辨兩個譜峰所需要的點(diǎn)數(shù)，繪制頻譜。（3）對模擬周期信號進(jìn)行譜分析選擇采樣頻率fs=64Hz，變換區(qū)間N=16，32，64，在三種情況進(jìn)行譜分析。分別打印其幅頻特性，并進(jìn)行分析和討論。（1）內(nèi)容1參考程序N=8。fs=16。n=0:N1。w=fs*n/N。x=2*sin(4*pi*n/fs)+5*cos(8*pi*n/fs)。X=fft(x,N)。figure(1)stem(w(1:N/2),abs(X(1:N/2)),39。k39。,39。filled39。)。grid。xlabel(39。Hz39。)。ylabel(39。|X(f)|39。)。title(39。N=839。)。% N取6或32時進(jìn)行對比討論（2）內(nèi)容2參考程序N=21。L=256。f1=120。f2=140。fs=400。T=1/fs。ws=2*pi*fs。n=0:N1。x=cos(2*pi*f1*n*T)+cos(2*pi*f2*n*T)。X=fftshift(fft(x,L))。w=(ws/2+(0:L1)*ws/L)/(2*pi)。figure(2)subplot(2,1,1)。plot(w,abs(X))。xlabel(39。Hz39。)。ylabel(39。幅度譜39。)。axis([200,200,0,15])。title(39。N=2139。)。N=11。n=0:N1。x=cos(2*pi*f1*n*T)+cos(2*pi*f2*n*T)。X=fftshift(fft(x,L))。subplot(2,1,2)。plot(w,abs(X))。xlabel(39。Hz39。)。ylabel(39。幅度譜39。)。axis([200,200,0,15])。title(39。N=1139。)。（3）內(nèi)容3參考程序figure(3)fs=64。T=1/fs。N=16。n=0:N1。 %FFT的變換區(qū)間N=16xnT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T)。 %對x(t)16點(diǎn)采樣Xk16=fft(xnT)。 %計算xnT的16點(diǎn)DFTXk16=fftshift(Xk16)。 %將零頻率移到頻譜中心 Tp=N*T。F=1/Tp。 %頻率分辨率Fk=N/2:N/21。fk=k*F。 %產(chǎn)生16點(diǎn)DFT對應(yīng)的采樣點(diǎn)頻率（以零頻率為中心）subplot(3,1,1)。stem(fk,abs(Xk16),39。.39。)。box on %繪制8點(diǎn)DFT的幅頻特性圖title(39。16點(diǎn)|DFT[x(nT)]|39。)。xlabel(39。f(Hz)39。)。ylabel(39。幅度39。)。axis([N*F/21,N*F/21,0,*max(abs(Xk16))])N=32。n=0:N1。 %FFT的變換區(qū)間N=16xnT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T)。 %對x(t)32點(diǎn)采樣Xk32=fft(xnT)。 %計算xnT的32點(diǎn)DFTXk32=fftshift(Xk32)。 %將零頻率移到頻譜中心 Tp=N*T。F=1/Tp。 %頻率分辨率Fk=N/2:N/21。fk=k*F。 %產(chǎn)生16點(diǎn)DFT對應(yīng)的采樣點(diǎn)頻率（以零頻率為中心）subplot(3,1,2)。stem(fk,abs(Xk32),39。.39。)。box on %繪制8點(diǎn)DFT的幅頻特性圖title(39。32點(diǎn)|DFT[x(nT)]|39。)。xlabel(39。f(Hz)39。)。ylabel(39。幅度39。)。axis([N*F/21,N*F/21,0,*max(abs(Xk32))])N=64。n=0:N1。 %FFT的變換區(qū)間N=16xnT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T)。 %對x(t)64點(diǎn)采樣Xk64=fft(xnT)。 %計算xnT的64點(diǎn)DFTXk64=fftshift(Xk64)。 %將零頻率移到頻譜中心 Tp=N*T。F=1/Tp。 %頻率分辨率Fk=N/2:N/21。fk=k*F。 %產(chǎn)生16點(diǎn)DFT對應(yīng)的采樣點(diǎn)頻率（以零頻率為中心）subplot(3,1,3)。stem(fk,abs(Xk64),39。.39。)。 box on%繪制8點(diǎn)DFT的幅頻特性圖title(39。 64點(diǎn)|DFT[x (nT)]|39。)。xlabel(39。f(Hz)39。)。ylabel(39。幅度39。)。axis([N*F/21,N*F/21,0,*max(abs(Xk64))])內(nèi)容圖 實驗內(nèi)容1采樣點(diǎn)數(shù)N=8時的頻譜圖 實驗內(nèi)容1采樣點(diǎn)數(shù)N=6時的頻譜圖 實驗內(nèi)容1采樣點(diǎn)數(shù)N=32時的頻譜 實驗內(nèi)容2頻譜圖（1）實驗內(nèi)容1，根據(jù)給定模擬=2Hz， =4Hz,所以滿足采樣定理，采取樣頻率=16Hz,F=42=2Hz,根據(jù)N2/F,所以N8，當(dāng)N=6時，不滿足采樣點(diǎn)數(shù)，所以頻譜在時域混疊，而當(dāng)N=32時，滿足采樣點(diǎn)數(shù),頻譜中能夠看到有兩個頻率成分，在=2Hz, =4Hz處有譜線，。(2)實驗內(nèi)容2, =120Hz, =140Hz,取采樣頻率=400Hz,F=140120=20Hz,所以N=20，顯然，當(dāng)N=21時，能把和 兩個譜峰分辨出來，當(dāng)N=11時，兩個譜峰分辨不出來。 實驗內(nèi)容3頻譜圖（3）實驗內(nèi)容3，對莫逆周期信號進(jìn)行譜分析有3個頻率成分，=4Hz, =8Hz，=10Hz,，采樣頻率=64Hz。變換區(qū)間N=16時，觀察時間=16T=,不是的整數(shù)倍周期，所以所得頻譜不正確，變換區(qū)間N=32，64時，觀察時間=,1s,是的整數(shù)周期，所以所得頻譜正確，圖中3根譜線正好位于4Hz,8Hz,10Hz處。變換區(qū)間N=64時，頻譜幅度是變換區(qū)間N=32時的2倍，這種結(jié)果正好驗證了用DFT對中期序列譜分析的理論。（1）完成各個實驗任務(wù)和要求，附上程序清單和有關(guān)曲線。（2）簡要回答思考問題