频响函数IFFT得脉冲响应函数

zengwj0423 · 发表于 2012-12-5 11:05

您需要登录才可以下载或查看，没有账号？我要加入

x

请问各位，我由试验测得的结构的频现函数，然后想得到脉冲响应函数，经过傅里叶逆变换可以获得，但是，IFFT后的时域脉冲响应函数为什么在一个采样周期内尾部值会突然变得很大，百思不得其解，望大神给予指点。
捕获.JPG

紫竹阁主 · 发表于 2012-12-5 11:15

问下，怎么通过傅里叶逆变换得到脉冲响应函数？

zengwj0423 · 发表于 2012-12-5 11:48

紫竹阁主发表于 2012-12-5 11:15

登录/注册后可看大图

问下，怎么通过傅里叶逆变换得到脉冲响应函数？

首先需要对频响函数扩展成共轭对称形式，再由IFFT计算得到

咯咯 · 发表于 2012-12-18 10:32

zengwj0423 发表于 2012-12-5 11:48

登录/注册后可看大图

首先需要对频响函数扩展成共轭对称形式，再由IFFT计算得到

楼主这问题解决了吗？我也遇到这种问题诶！

Jillian · 发表于 2012-12-18 10:46

fx=fft([x,zeros(1,length(y)-1)]);
fy=fft([y,zeros(1,length(x)-1)]);
g=ifft(fx.*fy)

是这样实现的吗？

zengwj0423 · 发表于 2012-12-18 17:46

Jillian 发表于 2012-12-18 10:46

登录/注册后可看大图

fx=fft([x,zeros(1,length(y)-1)]);
fy=fft([y,zeros(1,length(x)-1)]);
g=ifft(fx.*fy)

不是的，我是将频响函数扩展为伦奎斯特频率2倍的带宽，两侧数值关于这个频率共轭对称，然后再进行ifft变换

咯咯 · 发表于 2012-12-19 19:41

zengwj0423 发表于 2012-12-18 17:46

登录/注册后可看大图

不是的，我是将频响函数扩展为伦奎斯特频率2倍的带宽，两侧数值关于这个频率共轭对称，然后再进行ifft变换 ...

楼主你好！你那尾部畸变问题解决了吗？我也遇到类似的问题，不知道怎么处理

zengwj0423 · 发表于 2012-12-20 19:44

咯咯发表于 2012-12-19 19:41

登录/注册后可看大图

楼主你好！你那尾部畸变问题解决了吗？我也遇到类似的问题，不知道怎么处理

还没解决，应该和FFT算法无关，我认为是频响本身测试的问题，存在噪声干扰，或者结构非线性因素，对测试的比较理想的频响函数，IFFT得到的脉冲响应函数没有出现这种现象

咯咯 · 发表于 2012-12-20 22:13

zengwj0423 发表于 2012-12-20 19:44

登录/注册后可看大图

还没解决，应该和FFT算法无关，我认为是频响本身测试的问题，存在噪声干扰，或者结构非线性因素，对测试的 ...

我看文献上上说是吉布斯现象造成的，在函数的不连续处，就会出现畸变。离散的FFT变换相当于是加了个矩形窗截断了什么的，我也不太懂

zengwj0423 · 发表于 2012-12-21 09:59

咯咯发表于 2012-12-20 22:13

登录/注册后可看大图

我看文献上上说是吉布斯现象造成的，在函数的不连续处，就会出现畸变。离散的FFT变换相当于是加了个矩形 ...

你说的文献是哪篇，我试过了，加了汉宁窗处理后效果还是差，

咯咯 · 发表于 2012-12-21 11:24

zengwj0423 发表于 2012-12-21 09:59

登录/注册后可看大图

你说的文献是哪篇，我试过了，加了汉宁窗处理后效果还是差，

你看下这篇《有限区间信号边界效应问题的研究》。加窗以后，频域值反变换到时域不是已经不是原来的时域波形了吗，那有什么用啊？

smallqiao · 发表于 2013-4-7 09:38

这个问题解决了么？

zengwj0423 · 发表于 2013-4-7 13:03

smallqiao 发表于 2013-4-7 09:38

登录/注册后可看大图

这个问题解决了么？

应该是频响测得差的原因

smallqiao · 发表于 2013-4-11 20:04

zengwj0423 发表于 2013-4-7 13:03

登录/注册后可看大图

应该是频响测得差的原因

不仅仅这个问题，直接ifft就会有这种现象，我看到的文献也是这样的。自己做的仿真FRF也是这样的，应该不影响分析。频响函数大多在频域内计算

魔女 · 发表于 2013-7-31 15:03

zengwj0423 发表于 2012-12-18 17:46

登录/注册后可看大图

不是的，我是将频响函数扩展为伦奎斯特频率2倍的带宽，两侧数值关于这个频率共轭对称，然后再进行ifft变换 ...

楼主你好，本人刚学不怎么会编程，能不能给段程序，万分感谢！

账号		自动登录	找回密码
密码			我要加入

[FFT] 频响函数IFFT得脉冲响应函数