dsp研究性學(xué)習(xí)報(bào)告頻譜計(jì)算(編輯修改稿)
2024-09-10 21:21 本頁面
【文章內(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的哈明窗對(duì)信號(hào)截短后近似計(jì)算其頻譜。試用不同的A和B的取值，確定用哈明窗能分辯的最小的譜峰間隔中c的值?！绢}目分析】 本題討論用哈明窗計(jì)算序列頻譜時(shí)的頻率分辨率?！痉抡娼Y(jié)果】【結(jié)果分析】將實(shí)驗(yàn)結(jié)果與教材中定義的窗函數(shù)的有效寬度相比較，發(fā)表你的看法。W=(wsam*T/2+(0:L1)*wsam*T/L), wsam*T=2m==(1/2+(0:L1)/L)。m=c/N對(duì)應(yīng)畫出的頻譜中，兩個(gè)譜峰橫坐標(biāo)的間隔設(shè)為m, m=c/N.m=, N=64c=*64=近似等于2【自主學(xué)習(xí)內(nèi)容】始終令A(yù)=1，改變B的值（BA），控制變量來求得分辨情況?！鹃喿x文獻(xiàn)】【發(fā)現(xiàn)問題】 (專題研討或相關(guān)知識(shí)點(diǎn)學(xué)習(xí)中發(fā)現(xiàn)的問題)：fftshift的作用正是讓正半軸部分和負(fù)半軸部分的圖像分別關(guān)于各自的中心對(duì)稱。因?yàn)橹苯佑胒ft得出的數(shù)據(jù)與頻率不是對(duì)應(yīng)的，fftshift可以糾正過來【問題探究】在離散序列頻譜計(jì)算中為何要用窗函數(shù)？用不同的窗函數(shù)對(duì)計(jì)算結(jié)果有何影響？與矩形窗相比哈明窗有何特點(diǎn)?如何選擇窗函數(shù)？答：有一些離散序列長度無限長，計(jì)算機(jī)無法處理，所以要利用窗函數(shù)進(jìn)行截短。用不同的窗函數(shù)得到的計(jì)算結(jié)果不同，與矩形窗相比哈明窗減小了旁瓣，卻加寬了主瓣寬度。【仿真程序】A=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(
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1