声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: 幽幽深谷

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

[复制链接]
发表于 2010-7-15 15:44 | 显示全部楼层
本帖最后由 wdhd 于 2016-8-15 10:06 编辑
原帖由 幽幽深谷 于 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;
回复 支持 反对
分享到:

使用道具 举报

发表于 2010-8-22 14:13 | 显示全部楼层
回复 幽幽深谷 的帖子
楼主 你这个ERA算法真的可以用?  我发现你构造的汉克矩阵貌似有点问题吧?

    希望可以和您进一步交流
发表于 2010-10-16 08:17 | 显示全部楼层
我也非常想要这方面的程序,楼主人真不错!不过,有些地方不明白,
发表于 2012-4-10 21:22 | 显示全部楼层
回复 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;

请问这句是什么意思
谢谢
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 19:24 , Processed in 0.059885 second(s), 16 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表