原帖由 幽幽深谷 于 2007-6-26 22:11 发表
function =ssi2(n1,n2,h,d);
% generate R-matrix from h
=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,:);
=svd(pi);%如果不是为了求取oi-1,这一句有作用么?
xi1=pinv(oi)*pi; 回复 幽幽深谷 的帖子
楼主 你这个ERA算法真的可以用?我发现你构造的汉克矩阵貌似有点问题吧?
希望可以和您进一步交流 我也非常想要这方面的程序,楼主人真不错!不过,有些地方不明白, 回复 15 # secondye 的帖子
=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;
请问这句是什么意思
谢谢
页:
1
[2]