[問題] 做完ODE後轉FFT已回收
各位好：
我的ODE方程式 [M]y''+[K]y=[f]*sin(10*t) 假設矩陣都是6*6
我解完ODE後,將y的第一個解y1帶入FFT轉換
可以找到頻域的解 我用的是ODE45(解ODE比較慢)
ODE23tb(速度快多了)
然後我自己把原方程式多微分一次 [M]y'''+[K]y'=[f]*10*cos(10*t)
6*6的話降階就是18*18矩陣 y1~y6解出來是位移項
我一樣取y1下去轉FFT
當然第一個峰值的頻率跟我ODE沒多微分一次,算出來是一樣的
不過幅值差很大(10^3左右),也就是我多微分一次做出來的FFT
峰值比原來的小很多！
請問微分方程,多微分一次跟原始做出來不是應該會一樣嗎＠＠？
謝謝!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.66.117
※ 編輯: inoran54787 來自: 140.115.66.117 (04/16 18:10)
※ 編輯: inoran54787 來自: 122.116.13.191 (04/18 01:49)