ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠ΅ ΠΈ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΡΠ΅ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ
Π‘ ΠΏΠΎΠΌΠΎΡΡΡ ΡΠΎΠ·Π΄Π°Π½Π½ΠΎΠ³ΠΎ ΠΈΠ½ΡΠ΅ΡΡΠ΅ΠΉΡΠ° Π±ΡΠ»ΠΈ ΠΏΠΎΠ»ΡΡΠ΅Π½Ρ Π΄ΠΈΠ½Π°ΠΌΠΈΡΠ΅ΡΠΊΠΈΠ΅ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊΠΈ ΡΡΡΡ ΡΠΈΠΏΠΎΠ² ΡΡΠ΄ΠΎΠ². ΠΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ΅ (Π² ΠΏΡΠΎΡΠ΅Π½ΡΠ°Ρ ΠΎΡ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ) Π΄ΠΎΠΏΡΡΡΠΈΠΌΠΎΠ΅ ΠΈΠ·ΠΌΠ΅Π½Π΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ Π·Π° Π²ΡΠ΅ΠΌΡ Πt = 0.2Ρ: ΠΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ΅ (Π² ΠΏΡΠΎΡΠ΅Π½ΡΠ°Ρ ΠΎΡ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ) Π΄ΠΎΠΏΡΡΡΠΈΠΌΠΎΠ΅ ΠΈΠ·ΠΌΠ΅Π½Π΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ Π·Π° Π²ΡΠ΅ΠΌΡ Πt = 1 Ρ: ΠΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ΅ (Π² ΠΏΡΠΎΡΠ΅Π½ΡΠ°Ρ ΠΎΡ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ) Π΄ΠΎΠΏΡΡΡΠΈΠΌΠΎΠ΅ ΠΈΠ·ΠΌΠ΅Π½Π΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ Π·Π° Π²ΡΠ΅ΠΌΡ Πt = 1 Ρ: ΠΡΠ΅ΠΌΡ Π½Π°Π±ΠΎΡΠ°… Π§ΠΈΡΠ°ΡΡ Π΅ΡΡ >
ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠ΅ ΠΈ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΡΠ΅ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ (ΡΠ΅ΡΠ΅ΡΠ°Ρ, ΠΊΡΡΡΠΎΠ²Π°Ρ, Π΄ΠΈΠΏΠ»ΠΎΠΌ, ΠΊΠΎΠ½ΡΡΠΎΠ»ΡΠ½Π°Ρ)
- Π¦Π΅Π»Ρ ΡΠ°Π±ΠΎΡΡ ΠΈ ΠΈΡΡ ΠΎΠ΄Π½ΡΠ΅ Π΄Π°Π½Π½ΡΠ΅
- 1. ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ
- 2. ΠΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½Π°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ Π² ΡΡΠ΅Π΄Π΅ Matlab
- 3. ΠΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠΉ ΠΈΠ½ΡΠ΅ΡΡΠ΅ΠΉΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ (GUI)
- 4. Π Π΅Π·ΡΠ»ΡΡΠ°ΡΡ ΠΈΡΡΠ»Π΅Π΄ΠΎΠ²Π°Π½ΠΈΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ
ΠΡΠ²ΠΎΠ΄Ρ
Π¦Π΅Π»Ρ ΡΠ°Π±ΠΎΡΡ ΠΈ ΠΈΡΡ ΠΎΠ΄Π½ΡΠ΅ Π΄Π°Π½Π½ΡΠ΅
Π¦Π΅Π»Ρ ΠΊΡΡΡΠΎΠ²ΠΎΠΉ ΡΠ°Π±ΠΎΡΡ: ΡΠ°Π·ΡΠ°Π±ΠΎΡΠΊΠ° ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ, Π½Π°ΠΏΠΈΡΠ°Π½ΠΈΠ΅ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ ΠΈ ΠΈΡΡΠ»Π΅Π΄ΠΎΠ²Π°Π½ΠΈΠ΅ Ρ Π΅Π΅ ΠΏΠΎΠΌΠΎΡΡΡ Π΄ΠΈΠ½Π°ΠΌΠΈΡΠ΅ΡΠΊΠΈΡ ΡΠ²ΠΎΠΉΡΡΠ² ΠΌΠΎΠ΄Π΅Π»ΠΈΡΡΠ΅ΠΌΡΡ ΠΎΠ±ΡΠ΅ΠΊΡΠΎΠ² Π½Π° ΠΠ.
ΠΡΡ ΠΎΠ΄Π½ΡΠ΅ Π΄Π°Π½Π½ΡΠ΅
ΠΡΠ½ΠΎΠ²Π½ΡΠ΅ Π’Π’Π Π½Π°Π΄Π²ΠΎΠ΄Π½ΡΡ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ
ΠΡΡ ΠΎΠ΄Π½ΡΠ΅ Π΄Π°Π½Π½ΡΠ΅ ΠΈ ΠΎΠ³ΡΠ°Π½ΠΈΡΠ΅Π½ΠΈΡ, Π½Π°Π»ΠΎΠΆΠ΅Π½Π½ΡΠ΅ Π½Π° ΡΠΈΡΡΠ΅ΠΌΡ, ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»Π΅Π½Ρ Π² ΡΠ°Π±Π»ΠΈΡΠ΅ 1.
Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΠ΅ ΠΊΠΎΡΠ°Π±Π»Ρ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΡΠΉ
Π’Π°Π±Π»ΠΈΡΠ° 1
ΠΡΠ½ΠΎΠ²Π½ΡΠ΅ Π’Π’Π Π½Π°Π΄Π²ΠΎΠ΄Π½ΠΎΠ³ΠΎ ΠΊΠΎΡΠ°Π±Π»Ρ | |||||||
β Π²Π°ΡΠΈΠ°Π½ΡΠ° | Π’ΠΈΠΏ ΠΊΠΎΡΠ°Π±Π»Ρ | W, Ρ | N, Π».Ρ. | V, ΡΠ·Π»Ρ/VΠΊ, V1max, V2max, ΡΠ·Π»Ρ | |||
ΠΠ΅Π³ΠΊΠΈΠΉ ΠΊΡΠ΅ΠΉΡΠ΅Ρ «ΠΡΠ°ΡΠ½ΡΠΉ ΠΊΠ°Π²ΠΊΠ°Π· | |||||||
ΠΡΠ½ΠΎΠ²Π½ΡΠ΅ Π’Π’Π ΠΊΠ°ΡΠ΅ΡΠ° | |||||||
ΠΠ°ΡΠ΅Ρ «Π’ΡΠΈΡΠΎΠ½» | 2.3 | ||||||
ΠΡΠ½ΠΎΠ²Π½ΡΠ΅ Π’Π’Π ΠΊΠΎΡΠ°Π±Π»Ρ Π½Π° ΠΏΠΎΠ΄Π²ΠΎΠ΄Π½ΡΡ ΠΊΡΡΠ»ΡΡΡ | |||||||
— ΠΠΠ | 9.4 | ||||||
ΠΠ΅ΠΎΠ±Ρ ΠΎΠ΄ΠΈΠΌΠΎ ΡΡΠ΅ΡΡΡ:
1 ΡΠ·Π΅Π» = 1 ΠΌΠΈΠ»Ρ/ΡΠ°Ρ = 1852 ΠΌ / 3600 Ρ = 0.51 ΠΌ/Ρ
1 Π».Ρ. = 735.5 ΠΡ; 1 ΠΡ = 1 ΠΠΌ/Ρ
ΠΠ³ΡΠ°Π½ΠΈΡΠ΅Π½ΠΈΡ ΡΠ°Π·ΡΠ°Π±Π°ΡΡΠ²Π°Π΅ΠΌΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ:
1. ΠΡΠ΅Π΄ΠΏΠΎΠ»Π°Π³Π°Π΅ΡΡΡ, ΡΡΠΎ ΠΊΠΎΡΠ°Π±Π»Ρ Π΄Π²ΠΈΠΆΠ΅ΡΡΡ Π½Π° ΠΏΠΎΡΡΠΎΡΠ½Π½ΠΎΠΌ ΠΊΡΡΡΠ΅.
2. ΠΠ΅ ΡΡΠΈΡΡΠ²Π°Π΅ΡΡΡ Π²ΠΎΠ»Π½Π΅Π½ΠΈΠ΅ ΠΌΠΎΡΡ, Π³ΠΈΠ΄ΡΠΎΠ΄ΠΈΠ½Π°ΠΌΠΈΡΠ΅ΡΠΊΠΈΠ΅ ΠΎΡΠΎΠ±Π΅Π½Π½ΠΎΡΡΠΈ ΠΊΠΎΡΠΏΡΡΠ°, ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΠΎΠ΅ Π²ΠΎΠ·Π΄Π΅ΠΉΡΡΠ²ΠΈΠ΅ Π²Π΅ΡΡΠ° ΠΈ Ρ. ΠΏ.
1. ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ
ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ Π½Π°Π΄Π²ΠΎΠ΄Π½ΠΎΠ³ΠΎ Π²ΠΎΠ΄ΠΎΠΈΠ·ΠΌΠ΅ΡΠ°ΡΡΠ΅Π³ΠΎ ΠΊΠΎΡΠ°Π±Π»Ρ:
Π³Π΄Π΅ xi+1 — Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ Π½Π° ΡΠ»Π΅Π΄ΡΡΡΠ΅ΠΌ ΡΠ°Π³Π΅; xi — Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ Π½Π° ΡΠ΅ΠΊΡΡΠ΅ΠΌ ΡΠ°Π³Π΅; Πxi — ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΠ΅ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ Π½Π° ΠΏΡΠΎΡΠ»ΠΎΠΌ ΡΠ°Π³Π΅; Pi — ΠΎΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ Π² ΠΏΡΠΎΡΠ΅Π½ΡΠ°Ρ ΠΎΡ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ; Fmax — ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ; Πt — ΡΠ°Π³ ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΡ Π²ΡΠ΅ΠΌΠ΅Π½ΠΈ; A — ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½Ρ ΠΏΡΠΎΠΏΠΎΡΡΠΈΠΎΠ½Π°Π»ΡΠ½ΠΎΡΡΠΈ; Πxi+1 — ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΠ΅ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ Π½Π° ΡΠ΅ΠΊΡΡΠ΅ΠΌ ΡΠ°Π³Π΅; V2max — ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΊΠΎΡΠΎΡΡΠΈ ΠΊΠΎΡΠ°Π±Π»Ρ.
ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ Π½Π°Π΄Π²ΠΎΠ΄Π½ΠΎΠ³ΠΎ ΠΊΠΎΡΠ°Π±Π»Ρ «ΠΡΠ°ΡΠ½ΡΠΉ ΠΠ°Π²ΠΊΠ°Π·»
1. Π¨Π°Π³ ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΡ Π²ΡΠ΅ΠΌΠ΅Π½ΠΈ (ΠΏΡΠΈΠ½ΠΈΠΌΠ°Π΅ΡΡΡ) Πt = 1 Ρ.
2. ΠΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ΅ (Π² ΠΏΡΠΎΡΠ΅Π½ΡΠ°Ρ ΠΎΡ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ) Π΄ΠΎΠΏΡΡΡΠΈΠΌΠΎΠ΅ ΠΈΠ·ΠΌΠ΅Π½Π΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ Π·Π° Π²ΡΠ΅ΠΌΡ Πt = 1 Ρ:
Fmax = Nmax / Vmax=(55 000*735.5)/(29*0.51)=2 735 100 H
ΠFmax =0.1Fmax (Ρ.ΠΊ. ΠΊΠΎΡΠ°Π±Π»Ρ ΠΈΠΌΠ΅Π΅Ρ Π²ΠΎΠ΄ΠΎΠΈΠ·ΠΌΠ΅ΡΠ΅Π½ΠΈΠ΅ Π±ΠΎΠ»Π΅Π΅ 10 000Ρ) ΠPi Πt (ΠFmax/Fmax) 100 50%.
3. ΠΠΎΡΡΡΠΈΡΠΈΠ΅Π½Ρ ΠΏΡΠΎΠΏΠΎΡΡΠΈΠΎΠ½Π°Π»ΡΠ½ΠΎΡΡΠΈ ΡΠΈΠ»Ρ ΡΠΎΠΏΡΠΎΡΠΈΠ²Π»Π΅Π½ΠΈΡ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ:
=12 504
4. ΠΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΠ° ΠΏΡΠΎΠΉΠ΄Π΅Π½Π½ΠΎΠ³ΠΎ ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΡ:
xi+1 = xi + (xi-xi-1) + (Pi 27 351 0012/100 — 12504(xi — xi-1) xi — xi-1) /9 030 000.
5. ΠΠ΅Π»ΠΈΡΠΈΠ½Π° ΡΠ΅ΠΊΡΡΠ΅ΠΉ ΡΠΊΠΎΡΠΎΡΡΠΈ:
Vi+1 =—Πxi+1 / Πt = (xi+1 -xi) / 1
ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠ°ΡΠ΅ΡΠ° «Π’ΡΠΈΡΠΎΠ½»
1. Π¨Π°Π³ ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΡ Π²ΡΠ΅ΠΌΠ΅Π½ΠΈ (ΠΏΡΠΈΠ½ΠΈΠΌΠ°Π΅ΡΡΡ) Πt = 1 Ρ.
2. ΠΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ΅ (Π² ΠΏΡΠΎΡΠ΅Π½ΡΠ°Ρ ΠΎΡ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ) Π΄ΠΎΠΏΡΡΡΠΈΠΌΠΎΠ΅ ΠΈΠ·ΠΌΠ΅Π½Π΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ Π·Π° Π²ΡΠ΅ΠΌΡ Πt = 1 Ρ:
Fmax = Nmax / Vmax=(150*735.5)/(29*0.51)=7459 Π ΠFmax =0,2Fmax (Ρ.ΠΊ. ΠΊΠΎΡΠ°Π±Π»Ρ ΠΈΠΌΠ΅Π΅Ρ Π²ΠΎΠ΄ΠΎΠΈΠ·ΠΌΠ΅ΡΠ΅Π½ΠΈΠ΅ ΠΌΠ΅Π½Π΅Π΅ 10 000Ρ) ΠPi Πt (ΠFmax/Fmax) 100 10%.
3. ΠΠΎΡΡΡΠΈΡΠΈΠ΅Π½Ρ ΠΏΡΠΎΠΏΠΎΡΡΠΈΠΎΠ½Π°Π»ΡΠ½ΠΎΡΡΠΈ ΡΠΈΠ»Ρ ΡΠΎΠΏΡΠΎΡΠΈΠ²Π»Π΅Π½ΠΈΡ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ:
=34.1
4. ΠΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΠ° ΠΏΡΠΎΠΉΠ΄Π΅Π½Π½ΠΎΠ³ΠΎ ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΡ:
xi+1 = xi + (xi-xi-1)+(Pi 74 591/100−34.1(xi — xi-1) xi — xi-1) /2300.
5. ΠΠ΅Π»ΠΈΡΠΈΠ½Π° ΡΠ΅ΠΊΡΡΠ΅ΠΉ ΡΠΊΠΎΡΠΎΡΡΠΈ:
Vi+1 =—Πxi+1 / Πt = (xi+1 -xi) / 1
ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠΎΡΠ°Π±Π»Ρ Π½Π° ΠΏΠΎΠ΄Π²ΠΎΠ΄Π½ΡΡ ΠΊΡΡΠ»ΡΡΡ
Π³Π΄Π΅ xi+1 — Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ Π½Π° ΡΠ»Π΅Π΄ΡΡΡΠ΅ΠΌ ΡΠ°Π³Π΅; xi — Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ Π½Π° ΡΠ΅ΠΊΡΡΠ΅ΠΌ ΡΠ°Π³Π΅; Πxi — ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΠ΅ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ Π½Π° ΠΏΡΠΎΡΠ»ΠΎΠΌ ΡΠ°Π³Π΅; Pi — ΠΎΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ Π² ΠΏΡΠΎΡΠ΅Π½ΡΠ°Ρ ΠΎΡ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ; Fmax — ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ; Πt — ΡΠ°Π³ ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΡ Π²ΡΠ΅ΠΌΠ΅Π½ΠΈ; A — ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½Ρ ΠΏΡΠΎΠΏΠΎΡΡΠΈΠΎΠ½Π°Π»ΡΠ½ΠΎΡΡΠΈ; Πxi+1 — ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΠ΅ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ Π½Π° ΡΠ΅ΠΊΡΡΠ΅ΠΌ ΡΠ°Π³Π΅; V21max — ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΊΠΎΡΠΎΡΡΠΈ ΠΊΠΎΡΠ°Π±Π»Ρ ΠΏΡΠΈ Π²ΠΎΠ΄ΠΎΠΈΠ·ΠΌΠ΅ΡΠ°ΡΡΠ΅ΠΌ ΡΠ΅ΠΆΠΈΠΌΠ΅; V22max — ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΊΠΎΡΠΎΡΡΠΈ ΠΊΠΎΡΠ°Π±Π»Ρ ΠΏΡΠΈ ΡΠ΅ΠΆΠΈΠΌΠ΅ Π³Π»ΠΈΡΡΠΈΡΠΎΠ²Π°Π½ΠΈΡ ΠΈ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ Π½Π° ΠΊΡΡΠ»ΡΡΡ .
ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠΎΡΠ°Π±Π»Ρ Π½Π° ΠΏΠΎΠ΄Π²ΠΎΠ΄Π½ΡΡ ΠΊΡΡΠ»ΡΡΡ
1. Π¨Π°Π³ ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΡ Π²ΡΠ΅ΠΌΠ΅Π½ΠΈ (ΠΏΡΠΈΠ½ΠΈΠΌΠ°Π΅ΡΡΡ) Πt = 0.2Ρ.
2. ΠΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ΅ (Π² ΠΏΡΠΎΡΠ΅Π½ΡΠ°Ρ ΠΎΡ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ) Π΄ΠΎΠΏΡΡΡΠΈΠΌΠΎΠ΅ ΠΈΠ·ΠΌΠ΅Π½Π΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ Π·Π° Π²ΡΠ΅ΠΌΡ Πt = 0.2Ρ:
Fmax = Nmax / Vmax=(150*735,5)/(32*0,51)=18 027 H
ΠFmax =0,2Fmax (Ρ.ΠΊ. ΠΊΠΎΡΠ°Π±Π»Ρ ΠΈΠΌΠ΅Π΅Ρ Π²ΠΎΠ΄ΠΎΠΈΠ·ΠΌΠ΅ΡΠ΅Π½ΠΈΠ΅ ΠΌΠ΅Π½Π΅Π΅ 10 000Ρ) ΠPi Πt (ΠFmax/Fmax)100 4%.
3. ΠΠΎΡΡΡΠΈΡΠΈΠ΅Π½Ρ ΠΏΡΠΎΠΏΠΎΡΡΠΈΠΎΠ½Π°Π»ΡΠ½ΠΎΡΡΠΈ ΡΠΈΠ»Ρ ΡΠΎΠΏΡΠΎΡΠΈΠ²Π»Π΅Π½ΠΈΡ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ:
— ΠΏΡΠΈ VΠΊ:
— ΠΏΡΠΈ VVΠΊ: .
4. ΠΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΠ° ΠΏΡΠΎΠΉΠ΄Π΅Π½Π½ΠΎΠ³ΠΎ ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΡ:
ΠΡΠΈ VΠΊ:
xi+1 = xi + (xi — xi-1) + (Pi18 027 — 410(xi — xi-1) xi — xi-1) /9400.
ΠΡΠΈ VVΠΊ:
xi+1 = xi + (xi — xi-1) + (Pi18 027 — 68(xi — xi-1) xi — xi-1) /9400.
5. ΠΠ΅Π»ΠΈΡΠΈΠ½Π° ΡΠ΅ΠΊΡΡΠ΅ΠΉ ΡΠΊΠΎΡΠΎΡΡΠΈ:
Vi+1 = Πxi+1 / Π t = (xi+1 -xi) / 1.
2. ΠΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½Π°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ Π² ΡΡΠ΅Π΄Π΅ Matlab
function varargout = Kurs_gui (varargin)
% KURS_GUI M-file for Kurs_gui.fig
% KURS_GUI, by itself, creates a new KURS_GUI or raises the existing
% singleton*.
%
% H = KURS_GUI returns the handle to a new KURS_GUI or the handle to
% the existing singleton*.
%
% KURS_GUI ('CALLBACK', hObject, eventData, handles,…) calls the local
% function named CALLBACK in KURS_GUI.M with the given input arguments.
%
% KURS_GUI ('Property','Value',…) creates a new KURS_GUI or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before Kurs_gui_OpeningFunction gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to Kurs_gui_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE’s Tools menu. Choose «GUI allows only one
% instance to run (singleton)" .
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help Kurs_gui
% Last Modified by GUIDE v2.5 12-Dec-2011 16:21:28
% Begin initialization code — DO NOT EDIT
gui_Singleton = 1;
gui_State = struct ('gui_Name', mfilename, …
'gui_Singleton', gui_Singleton, …
'gui_OpeningFcn', @Kurs_gui_OpeningFcn, …
'gui_OutputFcn', @Kurs_gui_OutputFcn, …
'gui_LayoutFcn', [], …
'gui_Callback', []);
if nargin & isstr (varargin{1})
gui_State.gui_Callback = str2func (varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn (gui_State, varargin{:});
else
gui_mainfcn (gui_State, varargin{:});
end
% End initialization code — DO NOT EDIT
% —- Executes just before Kurs_gui is made visible.
function Kurs_gui_OpeningFcn (hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved — to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to Kurs_gui (see VARARGIN)
% Choose default command line output for Kurs_gui
handles.output = hObject;
% Update handles structure
guidata (hObject, handles);
% UIWAIT makes Kurs_gui wait for user response (see UIRESUME)
% uiwait (handles.figure1);
set (handles.W,'String', 0);
set (handles.N,'String', 0);
set (handles.V1,'String', 0);
set (handles.V1max,'String', 0);
set (handles.V2max,'String', 0);
% —- Outputs from this function are returned to the command line.
function varargout = Kurs_gui_OutputFcn (hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved — to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles. output;
% —- Executes on button press in ship1.
function ship1_Callback (hObject, eventdata, handles)
% hObject handle to ship1 (see GCBO)
% eventdata reserved — to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
if get (hObject,'Value')==1
set (handles.text1,'Visible','On');
set (handles.text2,'Visible','On');
set (handles.text3,'Visible','On');
set (handles.ship2,'Value', 0);
set (handles.ship3,'Value', 0);
set (handles.W,'Visible','On');
set (handles.N,'Visible','On');
set (handles.V1,'Visible','On');
set (handles.text6,'Visible','Off');
set (handles.text7,'Visible','Off');
set (handles.V1max,'Visible','Off');
set (handles.V2max,'Visible','Off');
set (handles.W,'String', 9030);
set (handles.N,'String', 55 000);
set (handles.V1,'String', 29);
else
set (handles.text1,'Visible','Off');
set (handles.text2,'Visible','Off');
set (handles.text3,'Visible','Off');
set (handles.W,'Visible','Off');
set (handles.N,'Visible','Off');
set (handles.V1,'Visible','Off');
end
% Hint: get (hObject,'Value') returns toggle state of ship1
% —- Executes on button press in ship2.
function ship2_Callback (hObject, eventdata, handles)
% hObject handle to ship2 (see GCBO)
% eventdata reserved — to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
if get (hObject,'Value')==1
set (handles.text1,'Visible','On');
set (handles.text2,'Visible','On');
set (handles.text3,'Visible','On');
set (handles.ship1,'Value', 0);
set (handles.ship3,'Value', 0);
set (handles.W,'Visible','On');
set (handles.N,'Visible','On');
set (handles.V1,'Visible','On');
set (handles.text6,'Visible','Off');
set (handles.text7,'Visible','Off');
set (handles.V1max,'Visible','Off');
set (handles.V2max,'Visible','Off');
set (handles.W,'String', 2.3);
set (handles.N,'String', 150);
set (handles.V1,'String', 29);
else
set (handles.text1,'Visible','Off');
set (handles.text2,'Visible','Off');
set (handles.text3,'Visible','Off');
set (handles.W,'Visible','Off');
set (handles.N,'Visible','Off');
set (handles.V1,'Visible','Off');
end
% Hint: get (hObject,'Value') returns toggle state of ship2
% —- Executes on button press in ship3.
function ship3_Callback (hObject, eventdata, handles)
% hObject handle to ship3 (see GCBO)
% eventdata reserved — to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
if get (hObject,'Value')==1
set (handles.text1,'Visible','On');
set (handles.text2,'Visible','On');
set (handles.text3,'Visible','On');
set (handles.text6,'Visible','On');
set (handles.text7,'Visible','On');
set (handles.ship1,'Value', 0);
set (handles.ship2,'Value', 0);
set (handles.W,'Visible','On');
set (handles.N,'Visible','On');
set (handles.V1,'Visible','On');
set (handles.V1max,'Visible','On');
set (handles.V2max,'Visible','On');
set (handles.W,'String', 9.4);
set (handles.N,'String', 400);
set (handles.V1,'String', 11);
set (handles.V1max,'String', 13);
set (handles.V2max,'String', 32);
else
set (handles.text1,'Visible','Off');
set (handles.text2,'Visible','Off');
set (handles.text3,'Visible','Off');
set (handles.text6,'Visible','Off');
set (handles.text7,'Visible','Off');
set (handles.W,'Visible','Off');
set (handles.N,'Visible','Off');
set (handles.V1,'Visible','Off');
set (handles.V1max,'Visible','Off');
set (handles.V2max,'Visible','Off');
end
% Hint: get (hObject,'Value') returns toggle state of ship3
% —- Executes on button press in close.
function close_Callback (hObject, eventdata, handles)
% hObject handle to close (see GCBO)
% eventdata reserved — to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
close
% —- Executes on button press in change.
function change_Callback (hObject, eventdata, handles)
% hObject handle to change (see GCBO)
% eventdata reserved — to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
if get (hObject,'Value')==1
set (handles.W,'Style','edit');
set (handles.N,'Style','edit');
set (handles.V1,'Style','edit');
set (handles.V1max,'Style','edit');
set (handles.V2max,'Style','edit');
end
% —- Executes on button press in ok.
function ok_Callback (hObject, eventdata, handles)
% hObject handle to ok (see GCBO)
% eventdata reserved — to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
if (get (handles.ship1,'Value')==1 || get (handles.ship2,'Value')==1)
KursoviK
elseif get (handles.ship3,'Value')==1
KursoviKPK
end
set (handles.text15,'Visible','On');
set (handles.text16,'Visible','On');
set (handles.text17,'Visible','On');
function tab_Callback (hObject, eventdata, handles)
% hObject handle to tab (see GCBO)
% eventdata reserved — to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
if (get (handles.ship1,'Value')==1 || get (handles.ship2,'Value')==1)
KursoviK
elseif get (handles.ship3,'Value')==1
KursoviKPK
end
B=[]; j=1;
if get (handles.ship1,'Value')==1
for i=1:5:length (tt),
B (j:)=[tt (i), XX (i), vv (i), pp (i)];
j=j+1;
end
end
i (get (handles.ship2,'Value')==1 || get (handles.ship3,'Value')==1)
for i=1:1:length (tt),
B (j:)=[tt (i), XX (i), vv (i), pp (i)];
j=j+1;
end
endsprintf ('%s',' ΠΡΠ΅ΠΌΡ ΠΡΡΡ Π‘ΠΊΠΎΡΠΎΡΡΡ Π’ΡΠ³Π°')
disp (B)
Π’Π΅ΠΊΡΡ ΡΠΊΡΠΈΠΏΡΠ° KursoviK. m
clc
W_h = handles. W; % Π²ΡΠ±ΠΎΡ ΠΎΠ±ΡΠ΅ΠΊΡΠ° Π°1 ΠΈΠ· ΡΠΎΠ΄ΠΈΡΠ΅Π»ΡΡΠΊΠΎΠ³ΠΎ ΠΎΠ±ΡΠ΅ΠΊΡΠ°
W = str2double (get (W_h,'String')); % ΠΈΠ·Π²Π»Π΅ΡΠ΅Π½ΠΈΠ΅ ΡΠΈΡΠ»Π΅Π½Π½ΠΎΠ³ΠΎ Π·Π½Π°ΡΠ΅Π½ΠΈΡ ΠΈΠ· ΠΎΠ±ΡΠ΅ΠΊΡΠ° Π°1
N_h = handles. N;
N = str2double (get (N_h,'String'));
V_h = handles. V1;
V = str2double (get (V_h,'String'));
% W=9030; N=55 000; V=29;
% 1 ΡΠ°Π³ ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΡ
dt=1;
% ΠΏΠ΅ΡΠ΅Π²ΠΎΠ΄ Π² ΡΠΈΡΡ Π‘ΠΈ
W1=W*1000; N1=N*735.5; V1=V*0.51;
Fmax=N1/V1;
if get (handles.ship1,'Value')==1
dFmax=0.1*Fmax;
end
if get (handles.ship2,'Value')==1
dFmax=0.2*Fmax;
end
% 2
dP=dt*(dFmax/Fmax)*100;
% 3
A=Fmax/(V12);
% ΠΡΠΎΠΉΠ΄Π΅Π½Π½ΠΎΠ΅ ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΠ΅
xi=0; x_pre=0; P=0; t=0;
global tt; global XX; global vv; global pp;
XX=[]; tt=[]; vv=[]; pp=[];
while P<100% ΡΠ°Π·Π³ΠΎΠ½ + Π²ΡΡ ΠΎΠ΄ Π½Π° ΠΌΠ°ΠΊΡ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ
X=xi+(xi-x_pre)+(P*Fmax*(dt2)/100-A*(xi-x_pre)*abs (xi-x_pre))/W1;
v=(X-xi)/dt;
vv=[vv, v];
pp=[pp, P];
x_pre=xi; xi=X;
P=P+dP;
tt=[tt, t];
XX=[XX, X];
t=t+dt;
end;
while v<0.98*V1, % ΡΠ°Π·Π³ΠΎΠ½ Π½Π° ΠΌΠ°ΠΊΡ ΡΡΠ³Π΅ Π΄ΠΎ ΠΌΠ°ΠΊΡ ΡΠΊΠΎΡΠΎΡΡΠΈ
X=xi+(xi-x_pre)+(P*Fmax*(dt2)/100-A*(xi-x_pre)*abs (xi-x_pre))/W1;
v=(X-xi)/dt;
vv=[vv, v];
pp=[pp, P];
x_pre=xi; xi=X;
tt=[tt, t];
XX=[XX, X];
t=t+dt;
end;
while P>-100, % ΡΠΎΡΠΌΠΎΠΆΠ΅Π½ΠΈΠ΅, Π²ΡΡ ΠΎΠ΄ Π½Π° ΠΎΠ±ΡΠ°ΡΠ½ΡΡ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ
X=xi+(xi-x_pre)+(P*Fmax*(dt2)/100-A*(xi-x_pre)*abs (xi-x_pre))/W1;
v=(X-xi)/dt;
vv=[vv, v];
x_pre=xi; xi=X;
pp=[pp, P];
P=P-dP;
tt=[tt, t];
XX=[XX, X];
t=t+dt;
end;
while v>0.02, % ΡΠΎΡΠΌΠΎΠΆΠ΅Π½ΠΈΠ΅, Π΄ΠΎ Π½ΡΠ»Ρ
X=xi+(xi-x_pre)+(P*Fmax*(dt2)/100-A*(xi-x_pre)*abs (xi-x_pre))/W1;
v=(X-xi)/dt;
vv=[vv, v];
pp=[pp, P];
x_pre=xi; xi=X;
tt=[tt, t];
XX=[XX, X];
t=t+dt;
end;
if get (handles.ok,'Value')==1
set (handles.XX_plot,'Visible','On');
axes (handles.XX_plot);
cla;
plot (tt, XX), grid;
set (handles.vv_plot,'Visible','On');
axes (handles.vv_plot);
cla;
plot (tt, vv), grid
set (handles.pp_plot,'Visible','On');
axes (handles.pp_plot);
cla;
plot (tt, pp), grid
end
Π’Π΅ΠΊΡΡ ΡΠΊΡΠΈΠΏΡΠ° KursoviKPK. m
clc
% W=9.4; N=400; V_k=11; V_1max=13; V_2max=32;
W_h = handles. W; % Π²ΡΠ±ΠΎΡ ΠΎΠ±ΡΠ΅ΠΊΡΠ° Π°1 ΠΈΠ· ΡΠΎΠ΄ΠΈΡΠ΅Π»ΡΡΠΊΠΎΠ³ΠΎ ΠΎΠ±ΡΠ΅ΠΊΡΠ°
W = str2double (get (W_h,'String')); % ΠΈΠ·Π²Π»Π΅ΡΠ΅Π½ΠΈΠ΅ ΡΠΈΡΠ»Π΅Π½Π½ΠΎΠ³ΠΎ Π·Π½Π°ΡΠ΅Π½ΠΈΡ ΠΈΠ· ΠΎΠ±ΡΠ΅ΠΊΡΠ° Π°1
N_h = handles. N;
N = str2double (get (N_h,'String'));
Vk_h = handles. V1;
V_k = str2double (get (Vk_h,'String'));
V1max_h = handles. V1max;
V_1max = str2double (get (V1max_h,'String'));
V2max_h = handles. V2max;
V_2max= str2double (get (V2max_h,'String'));
% 1 ΡΠ°Π³ ΠΏΡΠΈΡΠ°ΡΠ΅Π½ΠΈΡ
dt=1;
% ΠΏΠ΅ΡΠ΅Π²ΠΎΠ΄ Π² ΡΠΈΡΡ Π‘ΠΈ
W1=W*1000; N1=N*735.5; Vk=V_k*0.51; V1max=V_1max*0.51; V2max=V_2max*0.51;
Fmax=N1/V2max; dFmax=0.2*Fmax;
% 2
dP=dt*(dFmax/Fmax)*100;
% 3
A1=Fmax/(V1max2); A2=Fmax/(V2max2);
% ΠΡΠΎΠΉΠ΄Π΅Π½Π½ΠΎΠ΅ ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΠ΅
xi=0; x_pre=0; P=0; t=0; v=0;
global tt; global XX; global vv; global pp;
XX=[]; tt=[]; vv=[]; pp=[]; AA=[];
while P<100% ΡΠ°Π·Π³ΠΎΠ½ + Π²ΡΡ ΠΎΠ΄ Π½Π° ΠΌΠ°ΠΊΡ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ
if (v
A=A1;
end
if (v>=V1max)
A=A2;
end
if (Vk<=v & v
A=A1-(v-Vk)*(A1-A2)/(V1max-Vk);
end
AA=[AA, A];
X=xi+(xi-x_pre)+(P*Fmax*(dt2)/100-A*(xi-x_pre)*abs (xi-x_pre))/W1;
v=(X-xi)/dt;
vv=[vv, v];
pp=[pp, P];
x_pre=xi; xi=X;
P=P+dP;
tt=[tt, t];
XX=[XX, X];
t=t+dt;
end;
while v<0.98*V2max, % ΡΠ°Π·Π³ΠΎΠ½ Π½Π° ΠΌΠ°ΠΊΡ ΡΡΠ³Π΅ Π΄ΠΎ ΠΌΠ°ΠΊΡ ΡΠΊΠΎΡΠΎΡΡΠΈ
if (v
A=A1;
end
if (v>=V1max)
A=A2;
end
if (Vk<=v & v
A=A1-(v-Vk)*(A1-A2)/(V1max-Vk);
end
AA=[AA, A];
X=xi+(xi-x_pre)+(P*Fmax*(dt2)/100-A*(xi-x_pre)*abs (xi-x_pre))/W1;
v=(X-xi)/dt;
vv=[vv, v];
pp=[pp, P];
x_pre=xi; xi=X;
tt=[tt, t];
XX=[XX, X];
t=t+dt;
end;
while P>-100,
if (v
A=A1;
end
if (v>=V1max)
A=A2;
end
if (Vk<=v & v
A=A1-(v-Vk)*(A1-A2)/(V1max-Vk);
end
AA=[AA, A];
X=xi+(xi-x_pre)+(P*Fmax*(dt2)/100-A*(xi-x_pre)*abs (xi-x_pre))/W1;
v=(X-xi)/dt;
vv=[vv, v];
x_pre=xi; xi=X;
pp=[pp, P];
P=P-dP;
tt=[tt, t];
XX=[XX, X];
t=t+dt;
end;
while v>0.02,
if (v
A=A1;
end
if (v>=V1max)
A=A2;
end
if (Vk<=v & v
A=A1-(v-Vk)*(A1-A2)/(V1max-Vk);
end
X=xi+(xi-x_pre)+(P*Fmax*(dt2)/100-A*(xi-x_pre)*abs (xi-x_pre))/W1;
v=(X-xi)/dt;
vv=[vv, v];
pp=[pp, P];
x_pre=xi; xi=X;
tt=[tt, t];
XX=[XX, X];
t=t+dt;
end;
if get (handles.ok,'Value')==1
set (handles.XX_plot,'Visible','On');
axes (handles.XX_plot);
cla;
plot (tt, XX), grid;
set (handles.vv_plot,'Visible','On');
axes (handles.vv_plot);
cla;
plot (tt, vv), grid
set (handles.pp_plot,'Visible','On');
axes (handles.pp_plot);
cla;
plot (tt, pp), grid
end
3. ΠΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠΉ ΠΈΠ½ΡΠ΅ΡΡΠ΅ΠΉΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ (GUI)
ΠΠ»Ρ ΡΠ΄ΠΎΠ±Π½ΠΎΠ³ΠΎ Π²ΡΠ²ΠΎΠ΄Π° ΡΠ°ΡΡΡΠΈΡΠ°Π½Π½ΠΎΠΉ ΠΈΠ½ΡΠΎΡΠΌΠ°ΡΠΈΠΈ Π² ΡΡΠ΅Π΄Π΅ MATLAB Π±ΡΠ» ΡΠ°Π·ΡΠ°Π±ΠΎΡΠ°Π½ Π³ΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠΉ ΠΈΠ½ΡΠ΅ΡΡΠ΅ΠΉΡ, ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»Π΅Π½Π½ΡΠΉ Π½Π° ΡΠΈΡ. 1.
Π ΠΈΡ. 1 ΠΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠΉ ΠΈΠ½ΡΠ΅ΡΡΠ΅ΠΉΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ
ΠΠ»Ρ ΡΠΎΠ·Π΄Π°Π½ΠΈΡ ΠΈΠ½ΡΠ΅ΡΡΠ΅ΠΉΡΠ° Π±ΡΠ»ΠΈ ΠΈΡΠΏΠΎΠ»ΡΠ·ΠΎΠ²Π°Π½Ρ ΡΠ»Π΅Π΄ΡΡΡΠΈΠ΅ ΠΊΠΎΠΌΠΏΠΎΠ½Π΅Π½ΡΡ:
1. Radio Button — ΠΏΠ΅ΡΠ΅ΠΊΠ»ΡΡΠ°ΡΠ΅Π»Ρ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ.
2. AxesΠΏΠΎΠ»Ρ Π΄Π»Ρ ΠΏΠΎΡΡΡΠΎΠ΅Π½ΠΈΡ Π³ΡΠ°ΡΠΈΠΊΠΎΠ² ΠΏΡΠΎΠΉΠ΄Π΅Π½Π½ΠΎΠ³ΠΎ ΠΏΡΡΠΈ, ΡΠΊΠΎΡΠΎΡΡΠΈ ΠΈ ΡΠΈΠ»Ρ ΡΡΠ³ΠΈ.
3. Push Button — ΠΊΠ½ΠΎΠΏΠΊΠΈ. Change parameters — Π²ΠΎΠ·ΠΌΠΎΠΆΠ½ΠΎΡΡΡ ΠΈΠ·ΠΌΠ΅Π½ΡΡΡ Π·Π½Π°ΡΠ΅Π½ΠΈΡ ΠΎΡΠ½ΠΎΠ²Π½ΡΡ Π’Π’Π ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ. ΠΠ — ΠΊΠ½ΠΎΠΏΠΊΠ°, ΠΏΡΠΈ Π½Π°ΠΆΠ°ΡΠΈΠΈ Π½Π° ΠΊΠΎΡΠΎΡΡΡ ΠΏΡΠΎΠΈΡΡ ΠΎΠ΄ΡΡ ΠΎΡΠ½ΠΎΠ²Π½ΡΠ΅ ΡΠ°ΡΡΠ΅ΡΡ ΠΈ Π²ΡΠ²ΠΎΠ΄ Π³ΡΠ°ΡΠΈΠΊΠΎΠ². Close — ΠΊΠ½ΠΎΠΏΠΊΠ° Π·Π°ΠΊΡΡΡΠΈΡ ΠΎΠΊΠ½Π°.
4. Static Text — ΡΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠΉ ΡΠ΅ΠΊΡΡ. Π Π΄Π°Π½Π½ΠΎΠΉ ΡΠ°Π±ΠΎΡΠ΅ Π±ΡΠ»Π° ΡΠ΅Π°Π»ΠΈΠ·ΠΎΠ²Π°Π½Π° Π²ΠΎΠ·ΠΌΠΎΠΆΠ½ΠΎΡΡΡ ΠΈΠ·ΠΌΠ΅Π½Π΅Π½ΠΈΡ Π²ΠΈΠ΄ΠΈΠΌΠΎΡΡΠΈ ΡΠ΅Ρ ΠΈΠ»ΠΈ ΠΈΠ½ΡΡ ΡΠ°ΡΡΠ΅ΠΉ ΡΠ΅ΡΡΠ°.
5. Edit TextΠΈΠ·ΠΌΠ΅Π½ΡΠ΅ΠΌΡΠΉ ΠΏΠΎΠ»ΡΠ·ΠΎΠ²Π°ΡΠ΅Π»Π΅ΠΌ ΡΠ΅ΠΊΡΡ. ΠΠΎΠ·ΠΌΠΎΠΆΠ½ΠΎΡΡΡ ΡΠ΅Π΄Π°ΠΊΡΠΈΡΠΎΠ²Π°Π½ΠΈΡ ΡΠ΅ΠΊΡΡΠ° ΠΏΠΎΡΠ²Π»ΡΠ΅ΡΡΡ ΠΏΡΠΈ Π½Π°ΠΆΠ°ΡΠΈΠΈ ΠΊΠ½ΠΎΠΏΠΊΠΈ Change parameters.
6. Menu — ΠΌΠ΅Π½Ρ, ΠΊ ΠΊΠΎΡΠΎΡΠΎΠΌ Π±ΡΠ»Π° ΡΠ΅Π°Π»ΠΈΠ·ΠΎΠ²Π°Π½Π° Π²ΠΎΠ·ΠΌΠΎΠΆΠ½ΠΎΡΡΡ Π²ΡΠ²Π΅ΡΡΠΈ Π² ΠΊΠΎΠΌΠ°Π½Π΄Π½ΠΎΠ΅ ΠΎΠΊΠ½ΠΎ ΡΠ°Π±Π»ΠΈΡΡ Ρ ΡΠ΅Π·ΡΠ»ΡΡΠ°ΡΠ°ΠΌΠΈ ΡΠ°Π±ΠΎΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ.
4. Π Π΅Π·ΡΠ»ΡΡΠ°ΡΡ ΠΈΡΡΠ»Π΅Π΄ΠΎΠ²Π°Π½ΠΈΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ
Π‘ ΠΏΠΎΠΌΠΎΡΡΡ ΡΠΎΠ·Π΄Π°Π½Π½ΠΎΠ³ΠΎ ΠΈΠ½ΡΠ΅ΡΡΠ΅ΠΉΡΠ° Π±ΡΠ»ΠΈ ΠΏΠΎΠ»ΡΡΠ΅Π½Ρ Π΄ΠΈΠ½Π°ΠΌΠΈΡΠ΅ΡΠΊΠΈΠ΅ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊΠΈ ΡΡΡΡ ΡΠΈΠΏΠΎΠ² ΡΡΠ΄ΠΎΠ².
Π ΠΈΡ. 2 Π Π΅Π·ΡΠ»ΡΡΠ°Ρ ΡΠ°Π±ΠΎΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ, Π³ΡΠ°ΡΠΈΠΊΠΈ Π΄Π»Ρ Π»Π΅Π³ΠΊΠΎΠ³ΠΎ ΠΊΡΠ΅ΠΉΡΠ΅ΡΠ° «ΠΡΠ°ΡΠ½ΡΠΉ ΠΠ°Π²ΠΊΠ°Π·»
Π’Π°Π±Π»ΠΈΡΠ° 2
Π’Π°Π±Π»ΠΈΡΠ° Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΠΎΡΠ½ΠΎΠ²Π½ΡΡ ΠΏΠΎΠΊΠ°Π·Π°ΡΠ΅Π»Π΅ΠΉ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ:
Π‘ ΠΏΠΎΠΌΠΎΡΡΡ ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»Π΅Π½Π½ΡΡ ΡΠ΅Π·ΡΠ»ΡΡΠ°ΡΠΎΠ² ΡΠ°Π±ΠΎΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ Π±ΡΠ»ΠΎ ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΎ:
— Π²ΡΠ΅ΠΌΡ Π½Π°Π±ΠΎΡΠ° ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ ΡΠΊΠΎΡΠΎΡΡΠΈ 14,5ΠΌ/Ρ — 115 Ρ; ΠΏΡΠΈ ΡΡΠΎΠΌ ΠΊΠΎΡΠ°Π±Π»Ρ ΠΏΡΠΎΡ ΠΎΠ΄ΠΈΡ ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΠ΅ 1009,5 ΠΌ;
— Π²ΡΠ΅ΠΌΡ ΡΠΎΡΠΌΠΎΠΆΠ΅Π½ΠΈΡ — 50 Ρ, Π½Π° ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΠ΅ — 402 ΠΌ;
— ΠΎΠ±ΡΠ΅Π΅ Π²ΡΠ΅ΠΌΡ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠ°ΡΠ΅ΡΠ° ΡΠΎΡΡΠ°Π²ΠΈΠ»ΠΎ — 165 Ρ;
— ΠΎΠ±ΡΠ°Ρ ΠΏΡΠΎΠΉΠ΄Π΅Π½Π½Π°Ρ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΌ Π΄ΠΈΡΡΠ°Π½ΡΠΈΡ — 1555.8 ΠΌ.
Π ΠΈΡ. 3 Π Π΅Π·ΡΠ»ΡΡΠ°Ρ ΡΠ°Π±ΠΎΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ, Π³ΡΠ°ΡΠΈΠΊΠΈ Π΄Π»Ρ ΠΊΠ°ΡΠ΅ΡΠ° «Π’ΡΠΈΡΠΎΠ½»
Π’Π°Π±Π»ΠΈΡΠ° 3
Π’Π°Π±Π»ΠΈΡΠ° Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΠΎΡΠ½ΠΎΠ²Π½ΡΡ ΠΏΠΎΠΊΠ°Π·Π°ΡΠ΅Π»Π΅ΠΉ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ:
Π‘ ΠΏΠΎΠΌΠΎΡΡΡ ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»Π΅Π½Π½ΡΡ ΡΠ΅Π·ΡΠ»ΡΡΠ°ΡΠΎΠ² ΡΠ°Π±ΠΎΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ Π±ΡΠ»ΠΎ ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΎ:
— Π²ΡΠ΅ΠΌΡ Π½Π°Π±ΠΎΡΠ° ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ ΡΠΊΠΎΡΠΎΡΡΠΈ 14,6ΠΌ/Ρ — 12 Ρ; ΠΏΡΠΈ ΡΡΠΎΠΌ ΠΊΠΎΡΠ°Π±Π»Ρ ΠΏΡΠΎΡ ΠΎΠ΄ΠΈΡ ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΠ΅ 115.7 ΠΌ;
— Π²ΡΠ΅ΠΌΡ ΡΠΎΡΠΌΠΎΠΆΠ΅Π½ΠΈΡ — 9 Ρ, Π½Π° ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΠ΅ — 72.8 ΠΌ;
— ΠΎΠ±ΡΠ΅Π΅ Π²ΡΠ΅ΠΌΡ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠ°ΡΠ΅ΡΠ° ΡΠΎΡΡΠ°Π²ΠΈΠ»ΠΎ — 21 Ρ;
— ΠΎΠ±ΡΠ°Ρ ΠΏΡΠΎΠΉΠ΄Π΅Π½Π½Π°Ρ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΌ Π΄ΠΈΡΡΠ°Π½ΡΠΈΡ — 118.48 ΠΌ.
Π ΠΈΡ. 4 Π Π΅Π·ΡΠ»ΡΡΠ°Ρ ΡΠ°Π±ΠΎΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ, Π³ΡΠ°ΡΠΈΠΊΠΈ Π΄Π»Ρ ΠΊΠΎΡΠ°Π±Π»Ρ Π½Π° ΠΏΠΎΠ΄Π²ΠΎΠ΄Π½ΡΡ ΠΊΡΡΠ»ΡΡΡ
Π’Π°Π±Π»ΠΈΡΠ° 4
Π’Π°Π±Π»ΠΈΡΠ° Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΠΎΡΠ½ΠΎΠ²Π½ΡΡ ΠΏΠΎΠΊΠ°Π·Π°ΡΠ΅Π»Π΅ΠΉ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ:
Π‘ ΠΏΠΎΠΌΠΎΡΡΡ ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»Π΅Π½Π½ΡΡ ΡΠ΅Π·ΡΠ»ΡΡΠ°ΡΠΎΠ² ΡΠ°Π±ΠΎΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ Π±ΡΠ»ΠΎ ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΎ:
— Π²ΡΠ΅ΠΌΡ Π½Π°Π±ΠΎΡΠ° ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ ΡΠΊΠΎΡΠΎΡΡΠΈ 16.09ΠΌ/Ρ — 23 Ρ; ΠΏΡΠΈ ΡΡΠΎΠΌ ΠΊΠΎΡΠ°Π±Π»Ρ ΠΏΡΠΎΡ ΠΎΠ΄ΠΈΡ ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΠ΅ 243.5 ΠΌ;
— Π²ΡΠ΅ΠΌΡ ΡΠΎΡΠΌΠΎΠΆΠ΅Π½ΠΈΡ — 12 Ρ, Π½Π° ΡΠ°ΡΡΡΠΎΡΠ½ΠΈΠ΅ — 102.8 ΠΌ;
— ΠΎΠ±ΡΠ΅Π΅ Π²ΡΠ΅ΠΌΡ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΠΊΠ°ΡΠ΅ΡΠ° ΡΠΎΡΡΠ°Π²ΠΈΠ»ΠΎ — 35 Ρ;
— ΠΎΠ±ΡΠ°Ρ ΠΏΡΠΎΠΉΠ΄Π΅Π½Π½Π°Ρ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΌ Π΄ΠΈΡΡΠ°Π½ΡΠΈΡ — 346 ΠΌ.
ΠΡΠ²ΠΎΠ΄Ρ
ΠΡΠΈ Π²ΡΠΏΠΎΠ»Π½Π΅Π½ΠΈΠΈ Π΄Π°Π½Π½ΠΎΠΉ ΠΊΡΡΡΠΎΠ²ΠΎΠΉ ΡΠ°Π±ΠΎΡΡ Π±ΡΠ»ΠΈ ΡΠ°Π·ΡΠ°Π±ΠΎΡΠ°Π½Ρ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΠ΅ ΠΈ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΡΠ΅ ΠΌΠΎΠ΄Π΅Π»ΠΈ ΠΎΠΏΠΈΡΠ°Π½ΠΈΡ Π΄Π²ΠΈΠΆΠ΅Π½ΠΈΡ ΡΡΡΡ Π²ΠΈΠ΄ΠΎΠ² ΡΡΠ΄ΠΎΠ². Π ΡΠΎΠΎΡΠ²Π΅ΡΡΡΠ²ΠΈΠΈ Ρ ΠΈΡΡ ΠΎΠ΄Π½ΡΠΌΠΈ Π΄Π°Π½Π½ΡΠΌΠΈ ΠΏΡΠΈ ΠΏΠΎΠΌΠΎΡΠΈ ΡΠ°Π·ΡΠ°Π±ΠΎΡΠ°Π½Π½ΠΎΠ³ΠΎ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΠΎΠ³ΠΎ ΠΈΠ½ΡΠ΅ΡΡΠ΅ΠΉΡΠ° Π² MATLAB Π±ΡΠ»ΠΈ ΠΈΡΡΠ»Π΅Π΄ΠΎΠ²Π°Π½Ρ Π΄ΠΈΠ½Π°ΠΌΠΈΡΠ΅ΡΠΊΠΈΠ΅ ΡΠ²ΠΎΠΉΡΡΠ²Π° ΠΌΠΎΠ΄Π΅Π»ΠΈΡΡΠ΅ΠΌΡΡ ΠΊΠΎΡΠ°Π±Π»Π΅ΠΉ ΠΈ ΠΏΠΎΠ»ΡΡΠ΅Π½Ρ ΠΈΡ Π΄ΠΈΠ½Π°ΠΌΠΈΡΠ΅ΡΠΊΠΈΠ΅ Ρ Π°ΡΠ°ΠΊΡΠ΅ΡΠΈΡΡΠΈΠΊΠΈ.