3171
FFT法与直接法求自相关函数的比较clcf1015f2017f3026
132d1白噪声的方差取1Nle
gth
A1sqrt1032A2A1A3sqrt10272sA1cos2pif1
u
ifr
d02piA2cos2pif2
u
ifr
d02piA3cos2pif3
u
ifr
d02pi
ormr
d011N直接法求自相关函数r1zeros1Nk1fork10N1fork2k11Nr1kr1ksk2sk2k1e
dr1kr1kNkk1e
dFFT法求自相关函数sszeros1Nr2ffts2Nr2absr22Nr3ifftr2比较结果：两种方法求出的自相关函数相同。2）用BT法和周期图法分别进行功率谱估计用FFT法求出自相关函数s1szeros1Nr2ffts12Nr3absr22Nrifftr3r4rot90rot90r用BT法进行功率谱估计M64f05000205S1r1o
es1le
gthf
ffork2MS1S1rkexpj2pifk1r4k1expj2pifk1e
dsubplot211plotfS1redgridtitleBT法功率谱估计得到的曲线用周期图法进行功率谱估计Uzeros1le
gthffork1NUUskexpj2pifke
dS2absU2Nsubplot212plotfS2blackgridtitle周期图法功率谱估计得到的曲线仿真结果如下：
3
用AR模型法进行功率谱估计求自相关函数r1zeros1Nr2ffts2Nr3absr22Nrifftr3azeros1pa1zeros1p用Levi
so
Drubi
迭代算法求解AR模型系数form1p
fifm1a1r2r1Dr1absr22r1elsek0fori1m1kkairm1ie
dkkrm1Da1aamkforim111aia1ika1mie
dDD1absk2e
de
df05000205Szeros1le
gthfS1zeros1le
gthffori1pS1S1aiexpj2pifie
dS1absS112SDS1plotfSgridtitle16阶AR模型法功率谱估计曲线仿真结果如下：
320
用ROOTMUSIC法进行功率谱估计
fclcM8
11000Nle
gth
sexpj05pi
ju
ifr
d02piexpj03pi
ju
ifr
d02pi
ormr
d011Nr1ffts2N求自相关函数r2absr12Nrifftr22NRzerosMform1M求自相关矩阵for
1Mif
m0Rm
co
jr1m
1elseRm
r1
m1e
de
de
dUVeigRGU16提取噪声子空间向量V11symsza1z07a2z107fz7a1GGa2asym2polyfYrootsaY1absYIzeros14fori14寻找离单位圆最近的四个根（信号源有两个），因为根以共轭形式出现，所以有四个P2Iimi
absY11Y1Iii
fe
dfori14求出这四个根的相位，即信号频率的估计wia
gleYr