式中{Sk}=[S2 S4 S6 S8 ],系數(shù)矩陣[aij]在形式上與m=1時(shí)相同,但其中p=ω/3,而
式中各系數(shù)d形式上與m=1時(shí)相同,但其中p=ω/3。
n=1,m=5 時(shí),有方程組
[aij]{Sk}={bi} i,j=1,2,3,4, k=2j (4-61)
式中[aij]形式上與m=1時(shí)的(4-56)相同,但其中p=ω/5,而
式中各d系數(shù)形式上與m=1時(shí)相同,但其中p=ω/5。
n=1,m=7 時(shí),有方程組
[aij]{Sk}={bi} i,j=1,2,3,4, k=2j (4-63)
式中[aij]形式上與m=1時(shí)的(4-56)相同,但其中p=ω/7,而
式中各d系數(shù)形式上與m=1時(shí)相同,但其中p=ω/7。
通過(guò)求解以上幾組四元一次線性方程組可求得S2,S4,S6,S8。在此基礎(chǔ)上根據(jù)(4-19)、(4-20)式中ф1與θ1的關(guān)系和(4-53)、(4-54)可得:
n=1,m=1時(shí)
在以上工作的基礎(chǔ)上,便可根據(jù)系統(tǒng)周期解的一次表達(dá)式(4-38),得各種情況下的解析解:
m=1時(shí),系統(tǒng)的解析解為:
x=±(sinωt-εθ1cosωt) (4-69)
式中θ1為(4-65)式。
m=3,系統(tǒng)的解析解為:
式中θ1為(4-66)式。
m=5時(shí),系統(tǒng)的解析解為:
式中θ1為(4-67)式。
m=7時(shí),系統(tǒng)的解析解為:
式中θ1為(4-68)式。以上各式中上排符號(hào)對(duì)應(yīng)正,下排符號(hào)對(duì)應(yīng)負(fù)。由頻閃方程奇點(diǎn)的穩(wěn)定性可知,取正值時(shí)對(duì)應(yīng)的解是穩(wěn)定的。
以上各式中的值由下式計(jì)算:
代入(4-73),成為一個(gè)一元四次方程
a0R4+ a1R3+ a2R2+ a3R+a4=0 (4-75)
所以可用公式直播接求解(4-75),也可用數(shù)值求解。本文同時(shí)采用了公式法和數(shù)值法(牛頓-撒網(wǎng)格法)來(lái)求解,比較表明后者較好。
至此,全部求出了(4-2)式的近似解析解。
上一頁(yè)
下一頁(yè)