資源簡(jiǎn)介
首先產(chǎn)生K階Slepian窗的正交序列。在MATLAB仿真軟件中，實(shí)現(xiàn)Multitaper算法的函數(shù)為PMTM函數(shù)。PMTM函數(shù)使用的方法是改進(jìn)的周期圖法線性和非線性結(jié)合。從內(nèi)部參數(shù)和外部參數(shù)的角度分別來說明各個(gè)參數(shù)的作用及其對(duì)頻譜估計(jì)性能的影響。調(diào)整的參數(shù)分別為：Slepian序列的時(shí)間帶寬積，頻率域點(diǎn)數(shù)，輸入數(shù)據(jù)及其長(zhǎng)度，采樣頻率等。通過繪制估計(jì)得得頻譜圖來評(píng)判譜估計(jì)的性能。
代碼片段和文件信息
?t=[-1.024:0.001:1.023];
?y=sinc（2*pi*t）.*cos（200*pi*t）;
?plot（y）;
?syms?x
?y=sinc（2*pi*x）*cos（200*pi*x）*（heaviside（x+0.5）-heaviside（x-0.5））;
?fourier（y）
1/8*i*（-exp（1/2*i*w）*exp（-i*pi^2）*exp（-1/2*i*（w-2*pi^2+200*pi））*Ei（1-1/2*i*（w-2*pi^2+200*pi））-exp（1/2*i*w）*exp（-i*pi^2）*exp（-1/2*i*（w-2*pi^2-200*pi））*Ei（1-1/2*i*（w-2*pi^2-200*pi））+exp（-1/2*i*w）*exp（i*pi^2）*exp（1/2*i*（w-2*pi^2+200*pi））*Ei（11/2*i*（w-2*pi^2+200*pi））+exp（-1/2*i*w）*exp（i*pi^2）*exp（1/2*i*（w-2*pi^2-200*pi））*Ei（11/2*i*（w-2*pi^2-200*pi））+exp（1/2*i*w）*exp（i*pi^2）*exp（-1/2*i*（w+2*pi^2-200*pi））*Ei（1-1/2*i*（w+2*pi^2-200*pi））-exp（-1/2*i*w）*exp（-i*pi^2）*exp（1/2*i*（w+2*pi^2-200*pi））*Ei（11/2*i*（w+2*pi^2-200*pi））+exp（1/2*i*w）*exp（i*pi^2）*exp（-1/2*i*（w+2*pi^2+200*pi））*Ei（1-1/2*i*（w+2*pi^2+200*pi））-exp（-1/2*i*w）*exp（-i*pi^2）*exp（1/2*i*（w+2*pi^2+200*pi））*Ei（11/2*i*（w+2*pi^2+200*pi）））/pi^2
?w=[0:1:1000];
?f=w/（2*pi）;
?plot（fans）;
----------------------------------------------------------------------------------------評(píng)論
共有 條評(píng)論
川公網(wǎng)安備 51152502000135號(hào)