关于多自由度脉冲响应的程序问题
用下面程序段求脉冲响应,其中H11,H12,H21,H22是频响函数的分量subplot(221)
f=ifft(H11);
plot(w,f);
ylabel('IRF');
grid on
subplot(222)
f=ifft(H12);
plot(w,f);
ylabel('IRF');
grid on
subplot(223)
f=ifft(H21);
plot(w,f);
xlabel('Time(sec)');
ylabel('IRF');
grid on
subplot(224)
f=ifft(H22);
plot(w,f);
xlabel('Time(sec)');
ylabel('IRF');
grid on
结果如附件,好像不对呀!
那出问题呀!!请指教!
[ 本帖最后由 eight 于 2007-11-18 11:28 编辑 ] 这好像不是脉冲响应的曲线。 ifft是傅立叶反变换呀
结果确实有错误,不知哪出毛病了 1. ifft之后只取实部
2. 增加H11中的阻尼
3. 只留ifft之后一半
页:
[1]