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

luyunjun023

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

ChaChing

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

luyunjun023

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

luyunjun023

已解决！已解决！

ch_j1985

解决方法分享一下？

luyunjun023

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

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

luyunjun023

谢谢主任！

ChaChing

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

luyunjun023

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

livefairy

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

livefairy

楼主是做工作模态吧
这样来构造：
将每一采样时刻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

解决了？matlab程序会编么。指导指导啊