利用多窗法进行谱估计的仿真算法

首先产生K阶Slepian窗的正交序列。在MATLAB仿真软件中，实现Multitaper算法的函数为PMTM函数。PMTM函数使用的方法是改进的周期图法线性和非线性结合。从内部参数和外部参数的角度分别来说明各个参数的作用及其对频谱估计性能的影响。调整的参数分别为：Slepian序列的时间带宽积，频率域点数，输入数据及其长度，采样频率等。通过绘制估计得得频谱图来评判谱估计的性能。

 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）;

----------------------------------------------------------------------------------------