首先解决您进行的对比有无可比性的问题；
另外要对计算的每一步都要进行‘精密’的分析,做傅立叶变换,采用的是甚么样的定义等等；
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
再问： 我给你邮箱发了一个联系方式，还包括两段程序。谢谢
