|
3#
發(fā)表于 2014-6-4 14:27:56
|
只看該作者
這個(gè)是很基礎(chǔ)的微分模型,很多書里應(yīng)該都有的* ~9 E% c' s1 ?* @! W" p( V8 h
我說說我的思路:
5 o% t9 p: R5 m, t0 W由牛頓定律可得m*dv/dt=-k*s(t)
6 h1 p1 g9 @# [+ D顯然這是個(gè)二階微分方程,直接解出很困難,不過初中老師都講過,彈簧振子的位移是個(gè)周期函數(shù),
" e' H. ^, Q' d根據(jù)傅里葉大哥的說法,任何周期函數(shù)都可以轉(zhuǎn)化為三角函數(shù)的組合,所以不妨令位移方程s(t)=Asin(ωt) / X* D/ R7 Z* q; A, m' ^
則dv/dt=-ω^2*Asin(ωt)
6 M3 I) Z, ~1 p# k0 X帶入微分方程,可解的ω,至于振幅A,就是你自己設(shè)定彈簧最大位移7 _" [" M0 U9 H
3 I$ W& H9 r5 r/ h* X" p; w另外微分方程中可加入阻力,想要簡單的算,你可以設(shè)為一個(gè)常數(shù)C,愿意動(dòng)腦筋,可以試試將阻力設(shè)為與速度成正比的函數(shù) |
|