怎么就是求不出来，大家帮忙看看错在哪？

Juneliu

输入和输出均为平稳随机信号，信号见附件，我要想求的是结构的脉冲响应函数。因脉冲响应函数是传递函数的逆傅里叶变换，因此，我先通过信号的自谱和互谱密度求出传递函数（因是随机信号，非能量信号，不能直接做FFT），然后对其逆傅里叶变换来得到脉冲响应函数。程序如下：

clear
load sig.txt;
x=sig(:,1); %激励信号
y=sig(:,2); %响应信号（应变）
dt=0.05;%采样间隔0.05s
fs=1/0.05;%采样频率
nfft=1024;

Pxx=PSD(x,nfft,fs,[]);  %计算自功率谱
Pxy=CSD(x,y,nfft,fs,[]);  %计算互功率谱

Hw=Pxy./Pxx; ％传递函数

mc=ifft(Hw);％脉冲函数

plot(real(mc))


为了研究主要方法是否正确性，上面的过程对横坐标的问题这里先不考虑了，照上画出的图如下，很明显，这和脉冲响应函数相去甚远啊。我还试过其他的一些方法，好像都和这个结果差不多，不知道问题在哪，请高手指点一二啊，谢谢

[ 本帖最后由 zhangnan3509 于 2009-6-9 20:58 编辑 ]

jidianwangliang

输入端的响应和输出端的响应的互相关就可以代替这两个点间的脉冲响应函数
另外，用随机减量法也可以求出近似的脉冲响应函数