Кейбір сызықты емес дифференциалдық теңдеулер жүйесінің периодтық шешімдерін табу


Жұмыс түрі:  Материал
Тегін:  Антиплагиат
Көлемі: 15 бет
Таңдаулыға:   

Кейбір сызықты емес дифференциалдық теңдеулер жүйесінің периодтық шешімдерін табу.

1 Автономды сызықты емес жүйелердің периодты шешімдері

Егер бірінші жуықтау бойынша тұрақтылықты талдау кезінде бірінші жуықтау жүйесі тұрақты (стационар) болса, онда сызықты емес жүйенің орнықтылығы туралы мәселе

d x d t = A x + F ( t , x ) ( 1. 1 ) \frac{dx}{dt} = Ax + F(t, x) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 1)

мұндағы x = { x 1 , , x n } x = \left\{ x_{1}, \ldots, x_{n} \right\} , А матрицасының меншікті мәндерімен келесі жағдайларда шешіледі.

1. А матрицасының барлық меншікті мәндерінің теріс нақты бөліктері бар.

2. Бұл матрицаның кем дегенде бір меншікті мәнінің оң нақты бөлігі бар.

Бірінші жағдайда теңдеудің нөлдік шешімі

d x d t = A x ( 1. 2 ) \frac{dx}{dt} = Ax\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 2)

орнықты, яғни ол еркін үлкен бастапқы возмущений қатысты асимптотикалық орнықты, ал ( 1. 2 ) (1. 2) теңдеуінің нөлдік шешімі асимптотикалық орнықты болады, егер сызықты емес F ( t , x ) F(t, \ x) осы талапты қанағаттандырса:

F ( t , x ) x 0 x 0 \frac{\left\ F(t, \ x) \right\}{\left\ x \right\} \rightarrow 0\ \ \left\ x \right\ \rightarrow 0

t , t 0 < t < t, \ t_{0} < t < \infty\ бойынша равномерно . Екінші жағдайда ( 1. 1 ) \ (1. 1) және ( 1. 2 ) (1. 2) теңдеулерінің нөлдік шешімі F функцияға қатысты орнықсыз.

Алайда, егер А матрицасының кем дегенде бір меншікті мәнінің нөлдік нақты бөлігі бар болса, ал қалған ρ j \rho_{j} меншікті мәндері үшін R e ρ j > 0 {Re\ \rho}_{j} > 0 шарты орындалса, онда сызықты емес жүйенің орнықтылығы туралы мәселе, жалпы айтқанда, бірінші жуықтау жүйесімен шешілмейді. Төменде көрсетілгендей, нөлдік шешімнің маңайында ( 1. 1 ) (1. 1) теңдеуінің бастапқы возмущений қатысты периодты (орнықты немесе орнықсыз) шешімдері болуы мүмкін. Мұндай шешімдер сызықты емес жүйеде ( 1. 1 ) (1. 1) теңдеуінде сипатталған ауытқуларды анықтайды.

Осы бөлімде біз периодты шешімдерді құруға қатысты кейбір мәселелерді қарастырамыз. Қалай болғанда да, есептеудің нәтижелері екі жағдайға байланысты әр түрлі: ( 1. 1 ) (1. 1) теңдеуіндегі F ( t , x ) F(t, \ x) сызықты емес функциясы t t - дан тәуелді немесе тәуелсіз.

Бұл теориялық жағынан да, практикалық жағынан да алсақ, мұндай шешімдердің болу шарттарын алу сөзсіз қызығушылық тудырады.

1. 1 Квазисызықты автономды жүйелердің периодты шешімдері

Біз осы теңдеулермен сипатталған жүйенің әрекетін қарастырамыз:

d x i d t = F i ( x 1 , , x n , μ ) , i = 1 , 2 , , n , ( 1. 3 ) \frac{dx_{i}}{dt} = F_{i}\ \left( x_{1}, \ldots, x_{n}, \ \mu \right), \ \ \ i = 1, \ 2, \ldots\, \ n, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 3)

мұндағы μ \mu - сандық параметр. ( 1. 3 ) (1. 3) жүйесінің оң жақ бөліктері уақытқа тәуелсіз болғандықтан, оның x = x ( t , μ ) x = x(t, \ \mu) \ шешіміне x = x ( t t 1 , μ ) x = x\left( t - t_{1}, \ \mu \right) \ шешімі сәйкес келеді, мұндағы t 1 t_{1}\ - еркін момент уақыт. Сондықтан бұл жүйенің периодты шешімі (егер бар болса) t 1 t_{1} ретінде алуға болатын кем дегенде бір параметрге тәуелді. Басқа параметр, яғни μ \mu параметрі, периодты шешімдерді құру барысында пайда болатын қатарлардың жинақталуын қамтамасыз ету үшін аз деп саналады. Келесі төменде берілген шарттарды орындау арқылы осы шешімдерді құрамыз:

а) Кейбір x 1 , , x n x_{1}, \ldots, x_{n}\ айнымалылар аймағында және μ = 0 \mu = 0 нүктесінің кішкентай маңайында x 1 , , x n x_{1}, \ldots, x_{n} және μ \mu бойынша F i F_{i} функциялары аналитикалық болады;

б) μ = 0 \mu = 0 кезінде ( 1. 3 ) (1. 3) теңдеулері

d x i d t = k = 1 n a i k x k ( 1. 4 ) \frac{dx_{i}}{dt} = \sum_{k = 1}^{n}{a_{ik}x_{k}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 4)

периоды T болатын x = φ ( t ) x = \varphi(t) периодты шешімі бар сызықты жүйеге келеді.

Егер ( 1. 3 ) (1. 3) теңдеулер жүйесінің периодты шешімге бар болса, онда оның T 1 T_{1} периоды, жалпы айтқанда, T - дан өзгеше болады және оны осы түрде жазуға болады:

T 1 = T + α ( μ ) ( 1. 5 ) T_{1} = T + \alpha(\mu) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 5)

Пуанкер теоремасы 1. Егер ( 1. 4 ) (1. 4) теңдеуі а) және б) шарттарын қанағаттандырса, ал периоды T болатын μ = 0 \mu = 0 нүктесінде периодты шешімі бар болса, онда

x i = φ i ( t ) , i = 1 , 2 , , n , ( 1. 6 ) x_{i} = \varphi_{i}(t), \ \ i = 1, \ 2, \ \ldots\, \ n, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 6)

бастапқы шартқа сәйкес келетін ( 1. 3 ) (1. 3) жүйесінің шешімі

x i ( 0 ) = φ ( 0 ) + α i , i = 1 , 2 , , n , ( 1. 7 ) x_{i}(0) = \varphi(0) + \alpha_{i}, \ \ i = 1, \ 2, \ \ldots\, \ n, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 7)

және осындай түрде берілетін

x i = x i ( t , α 1 , , α n , μ ) , i = 1 , 2 , , n , ( 1. 8 ) x_{i} = x_{i}\left( t, \ \alpha_{1}, \ \ldots, \alpha_{n}, \mu \right), \ \ i = 1, \ 2, \ \ldots\, \ n, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 8)

болады:

1) α 1 , , α n \alpha_{1}, \ \ldots, \alpha_{n} және μ \mu - ге қатысты аналитикалық;

2) ( 1. 5 ) (1. 5) периодымен периодты, егер α 1 , , α n \alpha_{1}, \ \ldots, \alpha_{n} шамаларын периодтылық шартынан

кси x i ( T + α , α 1 , , α n , μ ) x i ( 0 , α 1 , , α n , μ ) = 0 , i = 1 , 2 , , n , ( 1. 9 ) \equiv x_{i}\left( T + \alpha, \ \alpha_{1}, \ldots, \alpha_{n}, \mu \right) - x_{i}\left( 0, \alpha_{1}, \ldots, \alpha_{n}, \mu \right) = 0, \ \ i = 1, \ 2, \ \ldots, \ n, \ \ \ \ \ (1. 9)

табуға болатын болса, ( б і р е у і н н ә л г е т е ң д е п а л ы п , м ы с а л ы , a n д і біреуін\ нәлге\ тең\ деп\ алып, \ \ мысалы, \ \ a_{n} - ді ), μ = 0 \mu = 0 үшін μ \mu аналитикалық функциялары нөлге айналады.

Бұл теореманың дәлелдеуін қарастырмай, оның қолдануын практикалық түрде қарастырайық. Ең алдымен, теореманың барлық шарттарын орындағанда жүйенің ( 1. 3 ) (1. 3) периодты шешімін ( 1. 8 ) (1. 8) кішкентай μ \mu - де жинақталатын қатар

x i ( t ) = φ i ( t ) + μ φ 1 i ( t ) + μ 2 φ 2 i ( t ) + . . . , i = 1 , 2 , , n , ( 1. 10 ) x_{i}(t) = \varphi_{i}(t) + \mu\varphi_{1i}(t) + \mu^{2}\varphi_{2i}(t) + \ . . . \ \, \ \ i = 1, \ 2, \ \ldots\, \ n, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 10)

ретінде беруге болады, мұндағы φ i ( t ) \varphi_{i}(t) - ( 1. 4 ) (1. 4) вырожденной жүйесінің периодты шешімдерін анықтайтын функциялар ( 1. 6 ) (1. 6) .

( 1. 9 ) (1. 9) периодтылықтың шартынан табу керек болатын параметрдің белгісіз функцияларының саны μ , α 1 , , α n {\mu, \alpha}_{1}, \ldots, \alpha_{n} теңдеулердің санынан 1 - ге көп. Алайда, α i \alpha_{i} сандарының бірі, мысалы, α n \alpha_{n} саны, нөлге тең деп алынады, сондықтан бұл неопределенность устраняется . Уақыттың басталу кезеңін немесе тұрақты t 1 t_{1} - ді (бұл туралы жоғарыда айтылған болатын) біз ерікті түрде таңдай аламыз. Общность решентя восстанавливается тем

Егер сызықты емес теңдеулер жүйесінің ( 1. 3 ) (1. 3) периодты шешімін ( ( 1. 5 ) (1. 5) периодымен) ( 1. 10 ) (1. 10) түрінде берсек, онда x i ( t ) x_{i}(t) функциялардың периодты болуынан:

x i ( t + T + α ) = x i ( t ) x_{i}(t + T + \alpha) = x_{i}(t)

T + α T + \alpha периодымен φ i k ( t ) \varphi_{ik}(t) функцилардың периодты болуы не следует , яғни ( 1. 10 ) (1. 10) функцияларының оң жақтары, жалпы айтқанда, осы теңдеулерді қанағаттандырмайды:

φ i ( t + T + α ) = φ i ( t ) , φ i k ( t + T + α ) = φ i k ( t ) \varphi_{i}(t + T + \alpha) = \varphi_{i}(t), \ {\ \ \ \varphi}_{ik}(t + T + \alpha) = \varphi_{ik}(t)

Бұл күтпеген факт болып көрінбеуі керек, себебі ( 1. 10 ) (1. 10) шешім μ \mu параметріне тәуелді, ал ( 1. 10 ) (1. 10) оң жағында осы параметрдің дәрежелері бойынша шешімдерді разложение берілген.

Мысалы, sin ( 1 + μ ) t \sin(1 + \mu) t функциясының периоды 2 π 1 + μ \ \frac{2\pi}{1 + \mu} . μ \mu параметрінің дәрежелері бойынша оны қатарға разложение , аламыз:

sin ( 1 + μ ) t = s i n t + μ t c o s t μ 2 t 2 2 s i n t + . . . \sin(1 + \mu) t = sint + \mu tcost - \frac{\mu^{2}t^{2}}{2}sint\ + \ . . . \

Осы теңдіктің оң жағындағы μ \mu ( μ , μ 2 , ) \left( \mu, \ \ \mu^{2}, \ \ldots \right) \ әрбір дәрежесіндегі қосындылар t возрастанием шектеусіз өседі. Осы жерден шығатын маңызды қорытынды аламыз, яғни ( 1. 3 ) (1. 3) жүйенің периодты шешімін жуықтап табу кезінде қатарлардың частичные қосындысын қолдану ұсынылмайды. Олар шешімнің периодты болуы туралы түсінік бермейді.

Периодты шешімнің практикалық жуықтауын құру үшін Ляпунов әдісін қолданамыз. Бұл алдымен басқа шамаға ауыстыру арқылы ( 1. 3 ) (1. 3) жүйенің преоброзуется к сообственному времени

τ = 1 T + α t ( 1. 11 ) \tau = \frac{1}{T + \alpha}t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 11)

Нәтижесінде жаңа автономды жүйе, периодтық шешімнің периоды алынатын уақытты өлшеу бірлігі шығады. μ \mu - ден тәуелді α \alpha шамасы осы параметрдің дәрежелері бойынша қатарға разлгагается . Сондықтан

T + α = T ( 1 + h 1 μ + h 2 μ 2 + . . . ) T + \alpha = T(1 + h_{1}\mu + h_{2}\mu^{2} + \ . . . ) .

Жүйенің меншікті колебаний көрінісі ( 1. 1 ) (1. 1) h i h_{i} коэффициенттеріне в разложении φ i k \varphi_{ik} функциялары үшін периодтылық шарттарын орындауға мүмкіндік береді.

Ұсынылған процедураның жалпы схемасы ұзақтау және қызықсыздау преобразование байланысты, сондықтан оны жалпы түрде бермейміз, салыстырмалы түрде қарапайым мысалды егжей-тегжейлі талдаумен шектелеміз.

Мысал 1. 1 Берілген теңдеудің периодты шешімін

d 2 x d t 2 + k 2 x = μ f ( x , d x d t ) ( 1. 12 ) \frac{d^{2}x}{{dt}^{2}} + k^{2}x = \mu f\left( x, \ \frac{dx}{dt} \right) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 12)

μ \mu бойынша аналитикалық, μ = 0 \mu = 0 болғанда, периодты шешімге келетін теңдеуді

d 2 x d t 2 + k 2 x = 0 ( 1. 13 ) \frac{d^{2}x}{{dt}^{2}} + k^{2}x = 0\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 13)

құрыңыз.

Меншікті уақытқа көшу үшін

τ = 2 π T + α t = k t ( 1 + h 1 μ + h 2 μ 2 + . . . ) 1 \tau = \frac{2\pi}{T + \alpha}t = kt\left( 1 + h_{1}\mu + h_{2}\mu^{2} + \ . . . \right) ^{- 1}

деп алып, мұндағы k = 2 π T k = \frac{2\pi}{T} - ( 1. 13 ) (1. 13) теңдеуінің периодты шешімінің частота. Онда

d x d t = d x d τ k ( 1 + h 1 μ + h 2 μ 2 + . . . ) 1 , \frac{dx}{dt} = \frac{dx}{d\tau}k\left( 1 + h_{1}\mu + h_{2}\mu^{2} + \ . . . \right) ^{- 1},

d 2 x d t 2 = d x d τ k 2 ( 1 + h 1 μ + h 2 μ 2 + . . . ) 2 , \frac{d^{2}x}{{dt}^{2}} = \frac{dx}{d\tau}k^{2}\left( 1 + h_{1}\mu + h_{2}\mu^{2} + \ . . . \right) ^{- 2},

және ( 1. 12 ) (1. 12) теңдеуін төмендегідей жазуға болады

d 2 x d τ 2 + x ( 1 + h 1 μ + h 2 μ 2 + ) 2 = \frac{d^{2}x}{{d\tau}^{2}} + x\left( 1 + h_{1}\mu + h_{2}\mu^{2} + \ \ldots \right) ^{2} =

= μ k 2 f ( x , k d x d τ ( 1 + h 1 μ + h 2 μ 2 + . . . ) 1 ) ( 1 + h 1 μ + h 2 μ 2 + . . . ) 2 . ( 1. 14 ) = \ \frac{\mu}{k^{2}}f\left( x, \ k\frac{dx}{d\tau}\left( 1 + h_{1}\mu + h_{2}\mu^{2} + \ \ . . . \right) ^{- 1} \right) \left( 1 + h_{1}\mu + h_{2}\mu^{2} + \ \ . . . \right) ^{2}. \ \ \ \ (1. 14)

Бұл теңдеудің шешімін мына түрде іздейміз

x ( τ ) = φ ( τ ) + μ φ 1 ( τ ) + μ 2 φ 2 ( τ ) + . . . , ( 1. 15 ) x(\tau) = \varphi(\tau) + \mu\varphi_{1}(\tau) + \mu^{2}\varphi_{2}(\tau) + \ . . . \ \, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 15)

мұндағы φ ( τ ) \varphi(\tau) - ( 1. 13 ) (1. 13) теңдеуінің шешімі, ал φ 1 ( τ ) \varphi_{1}(\tau) , . . . - периоды 2 π 2\pi болатын функциялар.

( 1. 15 ) (1. 15) функцияларын ( 1. 14 ) (1. 14) теңдеулеріне қойып, алатынымыз

d 2 φ d τ 2 + μ d 2 φ 1 d τ 2 + μ 2 d 2 φ 2 d τ 2 + . . . + ( φ ( τ ) + μ φ 1 ( τ ) + \frac{d^{2}\varphi}{{d\tau}^{2}} + \mu\frac{d^{2}\varphi_{1}}{{d\tau}^{2}} + \mu^{2}\frac{d^{2}\varphi_{2}}{{d\tau}^{2}} + \ . . . + \left( \varphi(\tau) + \mu\varphi_{1}(\tau) + \right. \

+ μ 2 φ 2 ( τ ) + . . . ) ( 1 + μ h 1 + μ 2 h 2 + . . . ) 2 = μ k 2 f {+ \left. \ \mu^{2}\varphi_{2}(\tau) + \ . . . \right) \left( 1 + {\mu h}_{1} + {\mu^{2}h}_{2} + \ \ . . . \right) }^{2} = \frac{\mu}{k^{2}}f \bullet

[ φ ( τ ) + μ φ 1 ( τ ) + μ 2 φ 2 ( τ ) + . . . , k ( d φ d τ + μ d φ 1 d τ + μ 2 d φ 2 d τ + . . . ) \bullet \left\lbrack \varphi(\tau) + \mu\varphi_{1}(\tau) + \mu^{2}\varphi_{2}(\tau) + \ . . . \, \ k\left( \frac{d\varphi}{d\tau} + \mu\frac{d\varphi_{1}}{d\tau} + \mu^{2}\frac{d\varphi^{2}}{d\tau} + . . . \right) \bullet \right. \

( 1 + h 1 μ + h 2 μ 2 + . . . ) 1 ] ( 1 + h 1 μ + h 2 μ 2 + . . . ) 2 . ( 1. 16 ) \left. \ \bullet \left( 1 + h_{1}\mu + h_{2}\mu^{2} + \ \ . . . \right) ^{- 1} \right\rbrack\left( 1 + h_{1}\mu + h_{2}\mu^{2} + \ \ . . . \right) ^{2}. \ \ \ (1. 16)

Осы теңдеудің сол және оң жақтарының μ \mu параметрінің бірдей дәрежелерінің коэффициенттерін теңестіріп, φ ( τ ) , φ 1 ( τ ) , φ 2 ( τ ) , \varphi(\tau), \ \ \varphi_{1}(\tau), \ \ \varphi_{2}(\tau), \ldots қатысты рекуррентті дифференциалды теңдеулер аламыз.

μ 2 \mu^{2} болғанда, алатынымыз

d 2 φ d τ 2 + φ = 0 ( 1. 17 ) \frac{d^{2}\varphi}{{d\tau}^{2}} + \varphi = 0\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 17)

және, демек, φ ( τ ) = A cos τ + B sin τ . \varphi(\tau) = A\cos\tau + B\sin\tau.

φ ( τ ) \varphi(\tau) ( 1. 17 ) (1. 17) теңдеуінің периодты шешімі болғандықтан, әрқашан осындай уақыт мезеті τ = t 1 \tau = t_{1} табылады

d φ d τ = 0 ( τ = t 0 ) . \frac{d\varphi}{d\tau} = 0\ \left( \tau = t_{0} \right) . \

t 1 = 0 t_{1} = 0 деп алсақ болады, сондықтан B = 0 B = 0 . Осылай φ ( τ ) \varphi(\tau) ретінде шешімді таңдаймыз

φ ( τ ) = A cos τ . \varphi(\tau) = A\cos\tau.

( 1. 16 ) (1. 16) теңдеуіндегі μ 1 \mu^{1} - дің коэффициенттерін теңестіріп, алатынымыз

d 2 φ 2 d τ 2 + φ 1 = 2 h 0 φ + 1 k 2 f ( φ , k d φ d τ ) \frac{d^{2}\varphi_{2}}{{d\tau}^{2}} + \varphi_{1} = - 2h_{0}\varphi + \frac{1}{k^{2}}f\left( \varphi, k\frac{d\varphi}{d\tau} \right)

және, демек,

d 2 φ 1 d τ 2 + φ 1 = 4 h 1 φ + 1 k 2 f ( A cos τ , k A sin τ ) . \frac{d^{2}\varphi_{1}}{{d\tau}^{2}} + \varphi_{1} = - 4h_{1}\varphi + \frac{1}{k^{2}}f\left( A\cos\tau, - kA\sin\tau \right) .

Бұл теңдеудің жалпы шешімін мына түрде беруге болады

φ 1 ( τ ) = A 1 cos τ + B 1 sin τ 2 h 1 A 0 τ cos s sin ( τ s ) d s + \varphi_{1}(\tau) = A_{1}\cos\tau + B_{1}\sin\tau - 2h_{1}A\int_{0}^{\tau}{\cos{s\ \sin(\tau - s) }}ds +

+ 1 k 2 0 τ f ( A cos s , k A sin s ) sin ( τ s ) d s . + \frac{1}{k^{2}}\int_{0}^{\tau}{f\left( A\cos s, - kA\sin s \right) }\sin(\tau - s) ds.

Периодтылықтың шарттары

φ 1 ( 2 π ) = φ 1 ( 0 ) , d φ 1 ( 2 π ) d τ = d φ 1 ( φ 1 ( 0 ) ) d τ \varphi_{1}(2\pi) = \varphi_{1}(0), \ \ \frac{d\varphi_{1}(2\pi) }{d\tau} = \frac{d\varphi_{1}\left( \varphi_{1}(0) \right) }{d\tau}

теңдеулерге алып келеді

2 h 1 A 0 2 π cos s sin s d s = 1 k 2 0 2 π f ( A cos s , k A sin s ) sin s d s , 2h_{1}A\int_{0}^{2\pi}{\cos s\sin sds =}\frac{1}{k^{2}}\int_{0}^{2\pi}{f\left( A\cos s, - kA\sin s \right) }\sin sds,

2 h 1 A 0 2 π cos 2 s d s = 1 k 2 0 2 π f ( A cos s , k A sin s ) cos s d s , 2h_{1}A\int_{0}^{2\pi}{\cos^{2}s\ ds =}\frac{1}{k^{2}}\int_{0}^{2\pi}{f\left( A\cos s, - kA\sin s \right) }\cos sds, \

немесе

0 2 π f ( A cos s , k A sin s ) sin s d s = 0 , \int_{0}^{2\pi}{f\left( A\cos s, - kA\sin s \right) }\sin sds = 0,

( 1. 18 ) (1. 18)

h 1 = 1 2 A π k 2 0 2 π f ( A cos s , k A sin s ) cos s d s . h_{1} = \frac{1}{2A\pi k^{2}}\int_{0}^{2\pi}{f\left( A\cos s, - kA\sin s \right) }\cos sds.

Теңдеулердің біріншісі ( 1. 18 ) (1. 18) А тұрақтысын анықтайды, осылайша ол сол

φ ( τ ) = A cos τ \varphi(\tau) = A\cos\tau

( 1. 17 ) (1. 17) теңдеуіндегі шешімдерді "таңдайды", около которых имеется периодическое решение уравнения ( 1. 12 ) (1. 12) . Теңдеулердің екіншісі ( 1. 18 ) (1. 18) дает поправку частоту колебаний системы в первом приближении. Бірінші жуықтауда

T 1 T ( 1 + μ h 1 ) T_{1} \approx T\left( 1 + \mu h_{1} \right)

деп алып және сызықты емес жүйенің ( 1. 18 ) (1. 18) тербеліс жиілігін ω \omega арқылы белгілеп, кіші μ \mu болғанда табамыз

ω k 1 + h 1 μ k ( 1 h 1 μ ) \omega \approx \frac{k}{1 + h_{1}\mu} \approx k\left( 1 - h_{1}\mu \right)

немесе

ω = k μ 2 π A k 0 2 π f ( A cos s , k A sin s ) cos s d s . \omega = k - \frac{\mu}{2\pi Ak}\int_{0}^{2\pi}{f\left( A\cos s, - kA\sin s \right) }\cos sds.

Сондай дәрежемен бұл формуланың орнына келесі формуланы алуға болады:

ω 2 = k 2 μ π A 0 2 π f ( A cos s , k A sin s ) cos s d s . \omega^{2} = k^{2} - \frac{\mu}{\pi A}\int_{0}^{2\pi}{f\left( A\cos s, - kA\sin s \right) }\cos sds.

Демек, ( 1. 12 ) (1. 12) теңдеуінің периодты шешімінің бірінші жуықтауын осындай түрде аламыз

x ( τ , μ ) = A cos τ + μ { A 1 cos τ + B 1 sin τ 2 h 1 A 0 τ cos s sin ( t s ) d s + x(\tau, \mu) = A\cos\tau + \mu\left\{ A_{1}\cos\tau + B_{1}\sin\tau - 2h_{1}A\int_{0}^{\tau}{\cos s\sin{(t - s) ds +}} \right. \

+ 1 k 2 0 τ f ( A cos s , k A sin s ) d s } , \left. \ + \frac{1}{k^{2}}\int_{0}^{\tau}{f\left( A\cos s, - kA\sin s \right) ds} \right\},

мұндағы A 1 A_{1} және B 1 B_{1} - еркін тұрақтылар, ал h 1 h_{1} және A A ( 1. 18 ) (1. 18) формулаларымен анықталады. Собственное время τ \tau t - мен τ = k ( 4 h 1 μ ) μ t \tau = k\left( 4 - h_{1}\mu \right) \mu t формуламен байланысты.

1. 2 Крылов әдісі

Ляпунов әдісі арқылы сызықты емес жүйенің периодты шешімінің барлық жуықтауларын табуға болады. Ол үшін ( 1. 16 ) (1. 16) теңдеуіндегі μ \mu - дің келесі бірдей дәрежелеріндегі коэффициенттерді теңестіру қажет, кезек бойынша екінші, үшінші және тағы сол сияқты жуықтауларды аламыз. Бірақ сонымен қатар қажетті есептеулердің көлемі артады. Бұл жағдай әдістің практикалық қолдануының мүмкіндіктерін шектейді. Оны қажет болған кезде тек бірінші немесе екінші жуықтауды құрғанда қолдану ұсынылады. Жуықтаудың жоғарғы ретті дәрежесін алу үшін екінші, үшінші және келесі жуықтауларды құруды айтарлықтай жеңілдететін Крылов әдісін қолдану ұсынылады. Біз оның мазмұнын қарапайым мысалда қарастырамыз.

Мысал 1. 2

Процесс Пуанкер теоремасының барлық шарттарын қанағаттандыратын теңдеумен сипатталсын

d 2 x d t 2 + k 2 x + μ x 3 = 0 . ( 1. 19 ) \frac{d^{2}x}{dt^{2}} + k^{2}x + \mu x^{3} = 0. \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 19)

Егер формальді құрылған разложении

x ( t ) = φ ( t ) + μ φ 1 ( t ) + μ 2 φ 2 ( t ) + . . . ( 1. 20 ) x(t) = \varphi(t) + \mu\varphi_{1}(t) + \mu^{2}\varphi_{2}(t) + \ . . . \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1. 20)

барлық функциялар φ ( τ ) , φ 1 ( τ ) , φ 2 ( τ ) , \varphi(\tau), \ \ \varphi_{1}(\tau), \ \ \varphi_{2}(\tau), \ldots периоды бірдей периодты функциялар анықталса, онда ( 1. 20 ) (1. 20) \ қатар ( 1. 19 ) (1. 19) теңдеуінің периодты шешімі болады. Периодтылықтың шарттарын орындау үшін Крылов тиімді әдіс ұсынды, яғни қатардың құрылуымен бірге μ \mu кіші параметрі бойынша жиіліктің квадратының разложение құрылады

p 2 = k 2 + h 1 μ + k 2 μ 2 + . . . p^{2} = k^{2} + h_{1}\mu + k_{2}\mu^{2} + . . . \

( 1. 19 ) (1. 19) теңдеуінің нақты μ \mu - дің екінші дәрежесіне дейін периодты шешімін алу үшін, табатынымыз

x ( t ) = φ ( t ) + μ φ 1 ( t ) + μ 2 φ 2 ( t ) , x(t) = \varphi(t) + \mu\varphi_{1}(t) + \mu^{2}\varphi_{2}(t),

( 1. 21 ) (1. 21)

p 2 = k 2 + h 1 μ + k 2 μ 2 . p^{2} = k^{2} + h_{1}\mu + k_{2}\mu^{2}.

Осы өрнектерді ( 1. 10 ) (1. 10) теңдеуіне қоямыз

d 2 φ ( t ) d t 2 + μ d 2 φ 1 ( t ) d t 2 + μ 2 d 2 φ 2 ( t ) d t 2 + ( p 2 h 1 μ h 2 μ ) [ φ ( t ) + \frac{d^{2}\varphi(t) }{{dt}^{2}} + \mu\frac{d^{2}\varphi_{1}(t) }{{dt}^{2}} + \mu^{2}\frac{d^{2}\varphi_{2}(t) }{{dt}^{2}} + \left( p^{2} - h_{1}\mu - h_{2}\mu \right) \left\lbrack \varphi(t) + \right. \

+ μ φ 1 ( t ) + φ 3 ( t ) ] + μ [ φ ( t ) + μ φ 1 ( t ) + μ 2 φ 4 ( t ) ] 3 = 0 \left. \ + \mu\varphi_{1}(t) + \varphi_{3}(t) \right\rbrack + \mu\left\lbrack \varphi(t) + \mu\varphi_{1}(t) + \mu^{2}\varphi_{4}(t) \right\rbrack^{3} = 0

μ \mu - дің бірдей дәрежелерінде коэффициенттерін нөлге теңестіріп, φ ( t ) , φ 1 ( t ) \varphi(t), \ \ \varphi_{1}(t) және φ 2 ( t ) \varphi_{2}(t) функцияларын табу үшін келесі теңдеулерді аламыз:

... жалғасы

Сіз бұл жұмысты біздің қосымшамыз арқылы толығымен тегін көре аласыз.
Ұқсас жұмыстар
Жәй дифференциалдық теңдеулер жүйесі үшін екі нүктелі импульсті түрткілі шеттік есептің шешімін табудың алгоритмін құру
Цилиндрлік функцияларды контурлық интегралдармен өрнектеу
ДИФФЕРЕНЦИАЛДЫҚ ТЕҢДЕУЛЕРДІ ОҚЫТУ ӘДІСТЕМЕСІ
Екінші ретті сызықтық дифференциалдық теңдеулер теориясы
N-ші ретті, коэффициенттері айнымалы біртекті сызықтық дифференциалдық теңдеулерді шешу әдістері
Дифференциалдық теңдеулер туралы жалпы түсінік
Автотербелмелі жүйелер кластерінің сигнал өндіру режимдері және оларға шуыл мен флуктуациялардың әсерін тәжірибе жүзінде зерттеу
Канал арқылы таралатын жалғыз толқын және оның параметрлері
n-ші ретті, коэффициенттері айнымалы біртекті сызықтық дифференциалдық теңдеулерді жалпыланған Абель формуласын пайдаланып шешу
Теңдеудің сандық шешімі
Пәндер



Реферат Курстық жұмыс Диплом Материал Диссертация Практика Презентация Сабақ жоспары Мақал-мәтелдер 1‑10 бет 11‑20 бет 21‑30 бет 31‑60 бет 61+ бет Негізгі Бет саны Қосымша Іздеу Ештеңе табылмады :( Соңғы қаралған жұмыстар Қаралған жұмыстар табылмады Тапсырыс Антиплагиат Қаралған жұмыстар kz