【文章內(nèi)容簡介】
。stem(2*k/N,Xm)。xlabel(39。Normalized frequency of 32 DFT39。)。ylabel(39。Magnitude39。)。grid on。[a1 b1]=max(Xm)。disp(39。peak value 39。)。disp(a139。)。disp(39。frequency of peak value39。)。disp(k(b1)39。)。(2)N=32k=0:N1。xk=sin(*pi*k)。L=input(39。the length of DFT:L=39。)。Xm=abs(fft(xk,L))。m=0:L1。stem(2*m/L,Xm)。xlabel([39。Normalized frequency of L=39。,num2str(L)])。hold on。ylabel(39。Magnitude39。)。grid on。[a1 b1]=max(Xm)。disp(39。peak value 39。)。disp(a139。)。disp(39。frequency of peak value39。)。disp(m(b1)39。)。M23 已知一離散序列為 x [k]=AcosW0k+Bcos ( (W0+DW)k)。用長度N=64的哈明窗對信號截短后近似計算其頻譜。試用不同的A和B的取值,確定用哈明窗能分辯的最小的譜峰間隔中c的值?!绢}目分析】 本題討論用哈明窗計算序列頻譜時的頻率分辨率?!痉抡娼Y(jié)果】【結(jié)果分析】將實驗結(jié)果與教材中定義的窗函數(shù)的有效寬度相比較,發(fā)表你的看法。W=(wsam*T/2+(0:L1)*wsam*T/L), wsam*T=2m==(1/2+(0:L1)/L)。m=c/N對應畫出的頻譜中,兩個譜峰橫坐標的間隔設(shè)為m, m=c/N.m=, N=64c=*64=近似等于2【自主學習內(nèi)容】始終令A=1,改變B的值(BA),控制變量來求得分辨情況?!鹃喿x文獻】【發(fā)現(xiàn)問題】 (專題研討或相關(guān)知識點學習中發(fā)現(xiàn)的問題):fftshift的作用正是讓正半軸部分和負半軸部分的圖像分別關(guān)于各自的中心對稱。因為直接用fft得出的數(shù)據(jù)與頻率不是對應的,fftshift可以糾正過來【問題探究】在離散序列頻譜計算中為何要用窗函數(shù)?用不同的窗函數(shù)對計算結(jié)果有何影響?與矩形窗相比哈明窗有何特點?如何選擇窗函數(shù)?答:有一些離散序列長度無限長,計算機無法處理,所以要利用窗函數(shù)進行截短。用不同的窗函數(shù)得到的計算結(jié)果不同,與矩形窗相比哈明窗減小了旁瓣,卻加寬了主瓣寬度?!痉抡娉绦颉緼=input(39。A=39。)。B=input(39。B=39。)。x=A*cos(100*pi*k)+B*cos((100*pi+*pi)*k)。N=64。L=512。wn=(hamming(N))39。xn=x.*wn。Xn=fftshift(fft(