svd_rank.m

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

%利用SVD—TLS算法确定AR阶数2p,给定阈值为0.9995
function [P,V,S]=svd_rank(Re)               %SVD-TLSd定阶函数
Pe=20;
[U,S,V]=svd(Re);                   %矩阵奇异值分解
for i=1:(Pe+1)                     %h=min(M,Pe+1)=Pe+1
    delta(i)=S(i,i);
end;
sum=0;
    for i=1:(Pe+1)                
        sum=(delta(i))^2+sum;
    end;
v0=sum;                           %所有奇异值平方和
vv=0;
sum=0;
j=0;
while vv<0.9995
    j=j+1;
    sum=(delta(j))^2+sum;
    vv=sqrt(sum/v0);                     %归一化比值
end;
P=j;                                 %定出阶数P