|
3#
發(fā)表于 2014-6-4 14:27:56
|
只看該作者
這個是很基礎(chǔ)的微分模型,很多書里應該都有的
& c& d' Y6 \0 {4 v' q' {6 |1 J+ y我說說我的思路:: M9 m7 B5 ~3 d! h* p8 E |; y
由牛頓定律可得m*dv/dt=-k*s(t)
0 }6 k+ I! }3 G4 x' K" x顯然這是個二階微分方程,直接解出很困難,不過初中老師都講過,彈簧振子的位移是個周期函數(shù),* `* K& N. B: Y6 L5 E8 a V
根據(jù)傅里葉大哥的說法,任何周期函數(shù)都可以轉(zhuǎn)化為三角函數(shù)的組合,所以不妨令位移方程s(t)=Asin(ωt) u O- y: p3 Z& D0 Q5 y' a6 z/ _
則dv/dt=-ω^2*Asin(ωt)6 G6 d% v" a" C v4 q1 z) y
帶入微分方程,可解的ω,至于振幅A,就是你自己設(shè)定彈簧最大位移# |3 @; B1 s. M9 o
' v& P, o3 \& b# v" L9 I- F另外微分方程中可加入阻力,想要簡單的算,你可以設(shè)為一個常數(shù)C,愿意動腦筋,可以試試將阻力設(shè)為與速度成正比的函數(shù) |
|