Question

求助，脉冲响应和频响函数的问题

为什么单自由度脉冲响应函数的傅里叶变换与频响函数，用matlab作图不一样呢，理论上应该一样啊

Answer 1

首先解决您进行的对比有无可比性的问题；
另外要对计算的每一步都要进行‘精密’的分析，做傅立叶变换，采用的是甚么样的定义等等；
matlab傅立叶变换怎么定义的、怎么进行的，这些东西都搞定了才能进行对比。
您说的频响函数是幅频特性、相频特性、复频特性还是实部、虚部？用电脑做DFT时，好多取采样间隔取1，而不采用真实的值，用户也不做修正，画出的图与理论曲线当然不一致。。。。。

追问

您好，我所说的脉冲响应函数就是h(t)=1/(m*wd)*exp(-ζ*wn*t)*sin(wd*t), wd为有阻尼结构的自振频率，wn为无阻尼自振频率，h(t)经过matlab中的fft变换的结果与频响函数H（w）=1/(k-m*w^2+i*c*w)进行的比较。matlab绘制H（w）的时候是画的abs(H(w))。

追答

这回说的很明白！
应当根据：H（w）=1/(k-m*w^2+i*c*w)，再求出|H（w）|之后与matlab绘制H（w）的时候是画的abs(H(w))去比较（这是幅频特性）。
？比较之后差别是甚么？

追问

我有点愚钝啊，不知道您说的什么意思呢。我的目的是想验证一下脉冲响应函数fft变换后的结果应该和H（w）的结果是一致的，如果一致，我就达到目的了。

追答

如果你理论推导：H（w）=1/(k-m*w^2+i*c*w)没问题，再把它变成幅频特性：|H(W)|,与matlab的幅频曲线比较，若不一致，那么matlab对h(t)的fft对不对？fft时参数设置对不对？。。。告诉你一招：给出8个数（自己给出），手工算出它的FFT值，再用matlab的FFT算这8个数的，两个能对上说明你算对了。这是基本功。真想搞懂每一步都得校核。我还是想问你：比较之后差别是甚么？

追问

这个图应该是在结构的自振频率处达到一样的峰值，但是现在频率对应上了，峰值不一样，差2个数量级呢。类似于这样的图

追答

这就涉及到matlab的FFT的具体参数设置问题，用的是DFT，Δt设置的是多少？设没设置？
设置的对不对？的吃透matlab关于FFT的原理，必须和手算的理论结果一致。问一下Δt设置的是多少？

追问

Δt设置了，但是对不对还真不知道，冒昧的问一下，有没有别的联系方式啊，我这个分快不够了啊，呵呵

追答

king_zyx@163.com

追问

我给你邮箱发了一个联系方式，还包括两段程序。谢谢