svd_tls.m

巧妙的运用matlab控件自编程序！很有意思哦！在信号分析领域用途广泛！

%利用SVD-TLS方法估计频率及参数
function [a,f]=svd_tls(P,Re)
Sp=zeros(P+1,P+1);
[P,V,S]=svd_rank(Re);    
for j=1:1:P
         for i=1:1:21-P
             Vj=V(:,j);
             Sp=Sp+S(j,j)^2*(Vj([i:1:i+P]))*(Vj([i:1:i+P]))';
         end
a=inv(Sp);
a=a(:,1)/a(1,1);                                      %AR参数
f=unique(abs(angle(roots(a))/(2*pi)));                %计算谐波频率
f(f>0.25)=0.5-f(f>0.25);
f(f==0)=[];
end