利用M精編B實現(xiàn)信DFT的計算精選文檔07級電信(2)班劉坤洋242、掌握利用MATLAB實現(xiàn)由DFT計算線性卷積的方法(1)fft(x)，fft(x，N)，ifft(x)，ifft(x，N)的含義及用法(2)在進行DFT時選取合適的時域樣本點數(shù)N請舉例，并編程實現(xiàn)22s源程序：>>N=30;%數(shù)據(jù)的長度>>L=512;%DFT的點數(shù)>>f1=100;f2=120;>>fs=600;%抽樣頻率>>T=1/fs;%抽樣間隔>>ws=2*pi*fs;>>t=(0:N-1)*T;>>f=cos(4*pi*f1*t)+cos(4*pi*f2*t);>>F=fftshift(fft(f,L));>>w=(-ws/2+(0:L-1)*ws/L)/(2*pi);>>hd=plot(w,abs(F));>>ylabel('幅度譜')>>xlabel('頻率/Hz')>>title('mypicture')(3)在對信號進行DFT時選擇hamming窗增加頻率分辨率請舉例，并編程實現(xiàn)22s源程序：>>N=50;%數(shù)據(jù)的長度>>L=512;%DFT的點數(shù)>>f1=100;f2=150;>>fs=600;%抽樣頻率>>T=1/fs;%抽樣間隔>>ws=2*pi*fs;>>t=(0:N-1)*T;>>f=cos(4*pi*f1*t)+*cos(4*pi*f2*t);>>wh=(hamming(N))';>>f=f.*wh;>>F=fftshift(fft(f,L));>>w=(-ws/2+(0:L-1)*ws/L)/(2*pi);>>plot(w,abs(F));>>ylabel('幅度譜')>>xlabel('頻率/Hz')>>title('mypicture')>>legend('N=50')請舉例，并編程實現(xiàn)>>k=0:N-1;>>L=0:511;>>x=cos(4*pi*k*4./N);>>X=fft(x);plot(k/N,abs(X),'o');>>holdon>>XE=fft(x,512);plot(L/512,abs(XE));>>legend('N=32')請舉例，并編程實現(xiàn)h[k]={2,2,1,1}>>x=[1211];>>h=[2211];>>%確定卷積結(jié)果的長度>>L=length(x)+length(h)-1;>>XE=fft(x,L);>.HE=fft(h,L);>>y1=ifft(XE.*HE);>>%繪制卷積結(jié)果及誤差圖>>k=0:L-1;>>subplot(1,2,1);>>stem(k,real(y1));axis([0607]);>>title('ResultoflinearConvolution');>>xlabel('Timeindexk');ylabel('Amplitude');>>y2=conv(x,h);error=y1-y2;>>subplot(1,2,2);stem(k,abs(err