共享ERA/DC算法和SSI算法程序

圣剑之灵

原帖由 幽幽深谷 于 2007-6-26 22:11 发表
function [fr,damp,v,d0]=ssi2(n1,n2,h,d);

% generate R-matrix from h
[r,c]=size(h);
r(1:c-n1)=0;
for q=1:c-n1;
    r(q)=0;
    for i=1:n1;
        r(q)=h(q+i)*h(i)+r(q);
    end;
end;

% g ...

%compuet xi1
clear yf p0 d0 q0 pr dd dr qr
pi=pref(2:r,:);
[p0 d0 q0]=svd(pi);%如果不是为了求取oi-1，这一句有作用么？
xi1=pinv(oi)*pi;

heshulin

回复 幽幽深谷 的帖子
楼主 你这个ERA算法真的可以用？  我发现你构造的汉克矩阵貌似有点问题吧？

    希望可以和您进一步交流

charlinew

我也非常想要这方面的程序,楼主人真不错!不过,有些地方不明白,

夏末

回复 15 # secondye 的帖子

[r,c]=size(h);
r(1:c-n1)=0;
for q=1:c-n1;
    r(q)=0;
    for i=1:n1;
        r(q)=h(q+i)*h(i)+r(q);
    end;
end;

请问这句是什么意思
谢谢