華南理工大學(xué)數(shù)字信號處理第三四次試驗答案.doc

實驗3 FFT算法的應(yīng)用一、實驗?zāi)康?、加深對離散信號的DFT的理解；2、在MATLAB中實現(xiàn)FFT算法。二、實驗內(nèi)容1、2N點實數(shù)序列N=64。用一個64點的復(fù)數(shù)FFT程序，一次算出，并繪出 的圖形。程序代碼如下：N=64;n1=0:2:(2*N-1);n2=1:2:(2*N-1);x1=cos(2*pi*7*n1/N)+1/2*cos(2*pi*9*n1/N);x2=cos(2*pi*7*n2/N)+1/2*cos(2*pi*9*n2/N);X1=fft(x1,N);X2=fft(x2,N);X=X1+X2;Y=abs(X);stem(Y);2、已知某序列在單位圓上的N=64等分樣點的Z變換為:。用N點IFFT程序計算出和。程序代碼如下：N=64;k=0:1:N-1;X=1./(1-0.8*exp(-j*2*pi*k/N);x=ifft(X,N);stem(abs(x);三、 利用至少兩種頻率估計方法編程求給定信號的頻率（要求運算過程中使用小于等于128點的FFT） 無噪聲、20dB、15dB、10dB、5dB、0dB、-5dB 已知信號的采樣率Fs= 4000 Hz方法一：譜峰檢測法程序代碼如下：function f=wxk3load signal; N=128; Fs=4000; f=zeros(1,7); y=abs(fft(x,N);m,n=max(y);f(1)=(n-1)*Fs/N; SNR=20 15 10 5 0 -5; for i=1:6 x=awgn(x,SNR(i); y=abs(fft(x,N); m,n=max(y); f(i+1)=(n-1)*Fs/N;endplot(f);hold;plot(f,.);axis(0 8 593 595);xlabel(無噪聲- 20dB - 15dB - 10dB - 5dB - 0dB - -5dB);legend(譜峰估計法);根據(jù)譜峰估計法的計算原理可知，所給信號的頻率介于593.6到593.8之間方法二：導(dǎo)數(shù)法程序代碼如下：function f=wxk4Fs=4000;N=128;n=1; delta=6;SNR=20 15 10 5 0 -5; load signal;signal=zeros(7,length(x);signal(1,:)=x;for i=1:6 signal(i+1,:)=awgn(x,SNR(i); endfor j=1:7 F=abs(fft(signal(j,:),N); ,k0=max(F); En1=0;En2=0; for k=k0-delta:k0+delta; En1=En1+(k-1)(2*n)*F(k)2; En2=En2+F(k)2; end f(j)=Fs/N*(En1/En2)(1/(2*n);endplot(f,r);hold;plot(f,.r);axis(0 8 570 620);xlabel(無噪聲 - 20dB - 15dB - 10dB - 5dB - 0dB - -5dB);l