%頻率估計--esprit
%----------------------------------------------
%初始條件
%----------------------------------------------
clear;?
clc;

for?L=1:20;

n=1:128;???????????????
xn=2*cos（2*pi*0.05*n）+3*cos（2*pi*0.40*n）+1.2*cos（2*pi*0.42*n）+sqrt（0.32）*randn（size（n））;
p=6;????????????????????????????%階數；
m=1.5*p;??????????????????????????
%-------------------------------------------------------------------------
%1、由xn求r（0）—r（m）
%-------------------------------------------------------------------------
N=length（xn）;
rxx1=xcorr（xn‘biased‘）;??????????%求自相關函數；
rxx=rxx1（N:m+N）;??????????????%取r（0）—r（m）；
%-------------------------------------------------------------------------
%2、構造自相關矩陣Rxx和Rxy
%-------------------------------------------------------------------------
Rxxx=toeplitz（rxx）;??????????????%m+1階自相關矩陣
Rxx=Rxxx（1:m1:m）;????????????%m階??Rxx
Rxy=Rxxx（1:m2:m+1）;??????????%m階??Rxy
%-------------------------------------------------------------------------
%3、對Rxx特征分解，求出最小特征值t
%-------------------------------------------------------------------------
ev=eig（Rxx）;????????????????????%矩陣分解；
[s?i]=min（ev）;
[V?D]=eig（Rxx）;
t=D（ii）;
%-------------------------------------------------------------------------
%?4、計算Rx=（Cxx，Cxy）=（Rxx-t*I，Rxy-t*Z）
%-------------------------------------------------------------------------
%4.1求Z
II=eye（m+1）;
I=eye（m）;
Z=II（1:m2:m+1）;
%4.2求（Cxx，Cxy）
Cxx=Rxx-t*I;
Cxy=Rxy-t*Z;
syms?x;
Rx=Cxx-x*Cxy;
%-------------------------------------------------------------------------
%5、求廣義特征根w?與頻率分量f
%-------------------------------------------------------------------------
y=det（Rx）;
yy=sym2poly（y）;?????????????%提取系數
w=roots（yy）;?

for?i=1:m;
f（Li）=angle（w（i））/（2*pi）;???????????????
end;

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1869??2011-01-05?18:05??MATLAB\f_est_esprit.m

?????文件???????2630??2011-01-05?18:05??MATLAB\f_est_music.m

?????文件????????798??2011-01-05?18:05??MATLAB\f_est_pisarenko.m

?????目錄??????????0??2010-05-17?23:50??MATLAB

-----------?---------??----------?-----??----

?????????????????5297????????????????????4