Π Π°ΡΡΠ΅Ρ ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊ ΠΎΠ±ΡΠ΅ΠΊΡΠ° Ρ ΠΏΠΎΠΌΠΎΡΡΡ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄ΠΈΠ½Π°ΠΌΠΈΠΊΠΈ
Π‘ΡΠ°Π²Π½ΠΈΡΡ ΡΠ΅Π·ΡΠ»ΡΡΠ°ΡΡ ΡΠ°ΡΡΠ΅ΡΠ° ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊ Ρ ΠΈΡΠΏΠΎΠ»ΡΠ·ΠΎΠ²Π°Π½ΠΈΠ΅ΠΌ ΠΌΠΎΠ΄Π΅Π»ΠΈ ΡΡΠ°ΡΠΈΠΊΠΈ (H1 — Π»Π°Π±. ΡΠ°Π±.1) ΠΈ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄ΠΈΠ½Π°ΠΌΠΈΠΊΠΈ (H3 — Π»Π°Π±. ΡΠ°Π±.3). ΠΎΠ±ΡΠ΅ΠΊΡ ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠΉ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊΠ° ΠΌΠΎΠ΄Π΅Π»Ρ ΠΠ±ΡΠ΅ΠΊΡ ΠΌΠΎΠ΄Π΅Π»ΠΈΡΠΎΠ²Π°Π½ΠΈΡ. Π Π°ΡΡΡΠΈΡΠ°ΡΡ ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠ΅ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊΠΈ ΠΎΠ±ΡΠ΅ΠΊΡΠ° ΠΏΠΎ ΠΈΡΡ ΠΎΠ΄Π½ΡΠΌ Π΄Π°Π½Π½ΡΠΌ, ΠΏΠΎΠ»ΡΡΠ΅Π½Π½ΡΠΌ Π΄Π»Ρ Π»Π°Π±. ΡΠ°Π±.1. Π‘ΠΎΡΡΠ°Π²ΠΈΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ ΡΠ°ΡΡΠ΅ΡΠ° ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊ ΠΎΠ±ΡΠ΅ΠΊΡΠ° Ρ ΠΏΠΎΠΌΠΎΡΡΡ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄ΠΈΠ½Π°ΠΌΠΈΠΊΠΈ. Π‘ΠΎΡΡΠ°Π²ΠΈΡΡ… Π§ΠΈΡΠ°ΡΡ Π΅ΡΡ >
Π Π°ΡΡΠ΅Ρ ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊ ΠΎΠ±ΡΠ΅ΠΊΡΠ° Ρ ΠΏΠΎΠΌΠΎΡΡΡ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄ΠΈΠ½Π°ΠΌΠΈΠΊΠΈ (ΡΠ΅ΡΠ΅ΡΠ°Ρ, ΠΊΡΡΡΠΎΠ²Π°Ρ, Π΄ΠΈΠΏΠ»ΠΎΠΌ, ΠΊΠΎΠ½ΡΡΠΎΠ»ΡΠ½Π°Ρ)
ΠΠ°ΠΏΠΈΡΠ°ΡΡ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠ΅ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄ΠΈΠ½Π°ΠΌΠΈΠΊΠΈ Π³ΠΈΠ΄ΡΠ°Π²Π»ΠΈΡΠ΅ΡΠΊΠΎΠΉ Π΅ΠΌΠΊΠΎΡΡΠΈ.
2. Π‘ΠΎΡΡΠ°Π²ΠΈΡΡ Π±Π»ΠΎΠΊ ΡΡ Π΅ΠΌΡ ΠΎΠ±ΡΠ΅ΠΊΡΠ°, ΡΠΊΠ°Π·Π°Π² ΡΠΎΡΡΠ°Π²Π»ΡΡΡΠΈΠ΅ Π²Π΅ΠΊΡΠΎΡΠ° Π²Ρ ΠΎΠ΄Π° Π₯ ΠΈ Π²Π΅ΠΊΡΠΎΡΠ° Π²ΡΡ ΠΎΠ΄Π° Y.
Π‘ΠΎΡΡΠ°Π²ΠΈΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ ΡΠ°ΡΡΠ΅ΡΠ° ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊ ΠΎΠ±ΡΠ΅ΠΊΡΠ° Ρ ΠΏΠΎΠΌΠΎΡΡΡ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄ΠΈΠ½Π°ΠΌΠΈΠΊΠΈ.
Π Π°ΡΡΡΠΈΡΠ°ΡΡ ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠ΅ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊΠΈ ΠΎΠ±ΡΠ΅ΠΊΡΠ° ΠΏΠΎ ΠΈΡΡ ΠΎΠ΄Π½ΡΠΌ Π΄Π°Π½Π½ΡΠΌ, ΠΏΠΎΠ»ΡΡΠ΅Π½Π½ΡΠΌ Π΄Π»Ρ Π»Π°Π±. ΡΠ°Π±.1.
ΠΠΎΡΡΡΠΎΠΈΡΡ Π³ΡΠ°ΡΠΈΠΊΠΈ ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ ΠΈ Π΄ΠΈΠ½Π°ΠΌΠΈΡΠ΅ΡΠΊΠΈΡ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊ.
Π‘ΡΠ°Π²Π½ΠΈΡΡ ΡΠ΅Π·ΡΠ»ΡΡΠ°ΡΡ ΡΠ°ΡΡΠ΅ΡΠ° ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊ Ρ ΠΈΡΠΏΠΎΠ»ΡΠ·ΠΎΠ²Π°Π½ΠΈΠ΅ΠΌ ΠΌΠΎΠ΄Π΅Π»ΠΈ ΡΡΠ°ΡΠΈΠΊΠΈ (H1 — Π»Π°Π±. ΡΠ°Π±.1) ΠΈ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄ΠΈΠ½Π°ΠΌΠΈΠΊΠΈ (H3 — Π»Π°Π±. ΡΠ°Π±.3). ΠΎΠ±ΡΠ΅ΠΊΡ ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠΉ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊΠ° ΠΌΠΎΠ΄Π΅Π»Ρ ΠΠ±ΡΠ΅ΠΊΡ ΠΌΠΎΠ΄Π΅Π»ΠΈΡΠΎΠ²Π°Π½ΠΈΡ.
ΠΡΡ ΠΎΠ΄Π½ΡΠ΅ Π΄Π°Π½Π½ΡΠ΅:
N. | G1. | G2. | P2. | P3. | K. | H. |
2.6. | 2.4. | 1.4. | 3.1. | 7.71 068. | 2.2. |
ΠΠΈΡΡΠΈΠ½Π³ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ.
Π€Π°ΠΉΠ»-ΡΡΠ΅Π½Π°ΡΠΈΠΉ: clc. clear. close all. global G1 G2 k P2 P3. G1=2.6; G2=2.4; P2=1.4; P3=3.1; k=7.71 068; nach=[0 2.2]; H0=0; n=10;del=0.1*k;k1=k/2; for i=1:n. if i==1. k=k1; [t, H]=ode45(@fun3,nach, H0); h (1)=real (H (end));K (1)=k; figure (1). plot (t, H). title ('H=f (t)'). xlabel ('t');ylabel ('H'); grid on; else. k=k1+del*i; [t, H]=ode45(@fun3,nach, H (end)); h (i)=real (H (end)); K (i)=k; % figure (i). % plot (t, H). end. end. figure (2). plot (t, H). title ('H=f (t)'). xlabel ('t');ylabel ('H'); grid on. figure (3). plot (K, h). title ('H=f (k)'). xlabel ('k');ylabel ('H'); grid on. k=7.71 068; G1=G½;h=[]; for i=1:n. if i==1. G1=G1; [t, H]=ode45(@fun3,nach, H0); h (1)=real (H (end));G3(1)=G1; else. G1=G1+del*i; [t, H]=ode45(@fun3,nach, H (end)); h (i)=real (H (end)); G3(i)=G1; end. end. figure (4). plot (G3,h). title ('H=f (G1)'). xlabel ('G1');ylabel ('H'); grid on. hold on. G1=2.6; G2=G2/2;G=[]; h=[]; for i=1:n. if i==1. G2=G2; [t, H]=ode45(@fun3,nach, H0); h (1)=real (H (end));G (1)=G2; else. G2=G2+del*i; [t, H]=ode45(@fun3,nach, H (end)); h (i)=real (H (end)); G (i)=G2; end. end. figure (5). plot (G, h). title ('H=f (G2)'). xlabel ('G2');ylabel ('H'); grid on. |
Π€Π°ΠΉΠ»-ΡΡΠ½ΠΊΡΠΈΡ: function dydt = fun3(t, H). global k. global G2. global G1. global P2. global P3. G3=k*sqrt (H+P2-P3); dydt=[G1+G2-G3]; end. |
ΠΡΠ°ΡΠΈΠΊΠΈ Π΄ΠΈΠ½Π°ΠΌΠΈΡΠ΅ΡΠΊΠΈΡ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊ.
ΠΡΠ°ΡΠΈΠΊΠΈ ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊ.