luyunjun023 发表于 2008-12-20 23:21

如何用互相关函数构造hankel矩阵

如何用互相关函数构造hankel矩阵

ChaChing 发表于 2008-12-21 12:40

怎麽许多人都喜欢一贴发各版块!? 增加积分吗?
若担心发错版块, 应该由管理阶层移动版块!?

luyunjun023 发表于 2009-1-3 23:02

抱歉！不是为了增加积分，是因为没有人回答，又比较急！

luyunjun023 发表于 2009-1-6 09:22

已解决！已解决！

ch_j1985 发表于 2009-1-6 20:16

回复 地板 luyunjun023 的帖子

解决方法分享一下？

luyunjun023 发表于 2009-1-8 15:19

yy=zeros(2*n1,n2);    %空hankel矩阵
for i=1:n2
    yy(:,i)=y(i:i+2*n1-1); %hankel矩阵第i列
end
yy; hankel_a=yy/((n2)^0.5);   %最后的hankel矩阵

y为信号数列
我是初学者，不知道对不对，见笑！我只是要构造一个hankel矩阵，不管什么方法。

请问怎么把自己的发帖删掉？我想把没有人回答的帖子删掉。

[ 本帖最后由 ChaChing 于 2009-6-15 10:58 编辑 ]

ch_j1985 发表于 2009-1-8 23:28

回复 6楼 luyunjun023 的帖子

帖子超过一定时间后自己是不能删除的，可以请版主帮忙！

luyunjun023 发表于 2009-1-9 12:10

回复 7楼 ch_j1985 的帖子

谢谢主任！

ChaChing 发表于 2009-1-9 13:36

回复 6楼 luyunjun023 的帖子

不清楚LZ说的"互相关函数"的意义!
matlab中是有现成的hankel函数, help hankel吧!
希望对LZ有用

luyunjun023 发表于 2009-1-14 18:40

help 过，没有用。多谢。已解决。
新的问题：
怎样生成干扰信号，不是零均值的白噪声。

livefairy 发表于 2009-6-14 21:09

回复 10楼 luyunjun023 的帖子

>> IRF=ifft(FRF,512);%FRF为响应信号
      >> max_irf=max(abs(IRF));
      >> noise=0.005*max_irf*randn(size(IRF));%构造噪声信号
      >> sprintf('SNR = %0.5g .',20*log10(std(IRF)/std(noise))),
      >> IRF=IRF+noise;%最后得到的信号
      >> FRF_noise=fft(IRF,512); %plot(w,20*real(log10(FRF_noise)),'b'),

livefairy 发表于 2009-6-14 21:17

楼主是做工作模态吧
这样来构造：
将每一采样时刻k的各相关函数值排列成矩阵。
file:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/T_1H5L5QX%7DB%7BCXQX$8%7D4@@7.jpg
Yk为离散时刻k的输出向量；Ymref为用作参考的输出向量
采样序号k＝1，……，q，……，p……p＋q－1，得到pxq个子矩阵
file:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/$959]THA%7DB2NLG%25$_LBRZ5Q.jpg

livefairy 发表于 2009-6-14 21:25

回复 地板 luyunjun023 的帖子

解决了？matlab程序会编么。指导指导啊
页: [1]
查看完整版本: 如何用互相关函数构造hankel矩阵