式中{Sk}=[S2 S4 S6 S8 ]，系數(shù)矩陣[aij]在形式上與m=1時相同极阴，但其中p=ω/3，而
式中各系數(shù)d形式上與m=1時相同弓洒，但其中p=ω/3。
n=1律罢，m=5 時膀值，有方程組
[aij]{Sk}={bi} i，j=1误辑，2沧踏，3，4巾钉， k=2j （4-61）
式中[aij]形式上與m=1時的（4-56）相同翘狱，但其中p=ω/5,而
式中各d系數(shù)形式上與m=1時相同，但其中p=ω/5砰苍。
n=1潦匈，m=7 時，有方程組
[aij]{Sk}={bi} i赚导，j=1茬缩，2，3吼旧，4寒屯， k=2j （4-63）
式中[aij]形式上與m=1時的（4-56）相同，但其中p=ω/7滩哥，而
式中各d系數(shù)形式上與m=1時相同，但其中p=ω/7昨镊。
通過求解以上幾組四元一次線性方程組可求得S2溃恼，S4，S6梁捉，S8追粉。在此基礎(chǔ)上根據(jù)（4-19）、（4-20）式中ф1與θ1的關(guān)系和（4-53）箍属、（4-54）可得：
n=1,m=1時
在以上工作的基礎(chǔ)上短硼，便可根據(jù)系統(tǒng)周期解的一次表達式（4-38）瓤梦，得各種情況下的解析解：
m=1時笔喉，系統(tǒng)的解析解為：
x=±（
sinωt-εθ1cosωt） （4-69）
式中θ1為（4-65）式。
m=3，系統(tǒng)的解析解為：
式中θ1為（4-66）式契耿。
m=5時，系統(tǒng)的解析解為：
式中θ1為（4-67）式光叉。
m=7時援奢，系統(tǒng)的解析解為：
式中θ1為（4-68）式。以上各式中上排符號對應(yīng)正
珍语，下排符號對應(yīng)負锤岸。由頻閃方程奇點的穩(wěn)定性可知，取正值時對應(yīng)的解是穩(wěn)定的板乙。
以上各式中的值由下式計算：
代入（4-73）是偷，成為一個一元四次方程
a0R4+ a1R3+ a2R2+ a3R+a4=0 （4-75）
所以可用公式直播接求解（4-75），也可用數(shù)值求解募逞。本文同時采用了公式法和數(shù)值法（牛頓-撒網(wǎng)格法）來求解蛋铆，比較表明后者較好。
至此凡辱，全部求出了（4-2）式的近似解析解戒职。
上一頁
下一頁
