ΠΡΠΎΠ΅ΠΊΡΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΡΠ»Π΅ΠΊΡΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΡΠ΅ΠΏΠΈ
Π ΠΎΠ±ΡΠ΅ΠΌ ΡΠ»ΡΡΠ°Π΅ ΡΠ΅ΡΠ΅Π½ΠΈΠ΅ Π·Π°Π΄Π°ΡΠΈ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ ΠΎΡΡΡΠ΅ΡΡΠ²Π»ΡΠ΅ΡΡΡ ΠΏΡΡΠ΅ΠΌ ΠΏΠ΅ΡΠ΅Π±ΠΎΡΠ° ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½Π½ΡΠΌ ΠΎΠ±ΡΠ°Π·ΠΎΠΌ Π²ΡΠ±ΠΈΡΠ°Π΅ΠΌΡΡ Π²Π°ΡΠΈΠ°Π½ΡΠΎΠ² Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ², ΡΡΠ°Π²Π½Π΅Π½ΠΈΡ ΠΈΡ ΠΌΠ΅ΠΆΠ΄Ρ ΡΠΎΠ±ΠΎΠΉ ΠΈ Π²ΡΠ±ΠΎΡΠ° Π½Π°ΠΈΠ»ΡΡΡΠ΅Π³ΠΎ Π²Π°ΡΠΈΠ°Π½ΡΠ°. ΠΠ»Π³ΠΎΡΠΈΡΠΌ Π²ΡΠ±ΠΎΡΠ° ΠΎΡΠ΅ΡΠ΅Π΄Π½ΠΎΠ³ΠΎ Π²Π°ΡΠΈΠ°Π½ΡΠ° Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ² Π½ΠΎΡΠΈΡ Π½Π°Π·Π²Π°Π½ΠΈΠ΅ ΡΡΡΠ°ΡΠ΅Π³ΠΈΠΈ ΠΏΠΎΠΈΡΠΊΠ°. Π¦Π΅Π»Ρ ΠΊΡΡΡΠΎΠ²ΠΎΠ³ΠΎ ΠΏΡΠΎΠ΅ΠΊΡΠΈΡΠΎΠ²Π°Π½ΠΈΡ ΡΠΎΡΡΠΎΠΈΡ Π² ΠΏΡΠΈΠΎΠ±ΡΠ΅ΡΠ΅Π½ΠΈΠΈ Π½Π°Π²ΡΠΊΠΎΠ² ΠΏΡΠΈΠΌΠ΅Π½Π΅Π½ΠΈΡ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΠΌΠΎΠ΄Π΅Π»ΠΈΡΠΎΠ²Π°Π½ΠΈΡ… Π§ΠΈΡΠ°ΡΡ Π΅ΡΡ >
ΠΡΠΎΠ΅ΠΊΡΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΡΠ»Π΅ΠΊΡΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΡΠ΅ΠΏΠΈ (ΡΠ΅ΡΠ΅ΡΠ°Ρ, ΠΊΡΡΡΠΎΠ²Π°Ρ, Π΄ΠΈΠΏΠ»ΠΎΠΌ, ΠΊΠΎΠ½ΡΡΠΎΠ»ΡΠ½Π°Ρ)
- ΠΠ²Π΅Π΄Π΅Π½ΠΈΠ΅
- 1. ΠΠ°Π΄Π°Π½ΠΈΠ΅ Π½Π° ΠΊΡΡΡΠΎΠ²ΠΎΠΉ ΠΏΡΠΎΠ΅ΠΊΡ
- 2. Π€ΠΎΡΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ
- 3. ΠΠΎΠΈΡΠΊ ΠΏΠ΅ΡΠΈΠΎΠ΄ΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ
- 4. ΠΠ°ΡΠ°ΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΠ»Π΅ΠΊΡΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΡΠ΅ΠΏΠΈ
- 5. Π Π΅Π·ΡΠ»ΡΡΠ°ΡΡ ΡΠ°Π±ΠΎΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ
- Π‘ΠΏΠΈΡΠΎΠΊ ΠΈΡΠΏΠΎΠ»ΡΠ·ΡΠ΅ΠΌΠΎΠΉ Π»ΠΈΡΠ΅ΡΠ°ΡΡΡΡ
- ΠΠΈΡΡΠΈΠ½Π³ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ
- ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΠΊΠ° ΡΠ»Π΅ΠΊΡΡΠΈΡΠ΅ΡΠΊΠΈΠΉ ΡΠ΅ΠΏΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ°
- ΠΠ²Π΅Π΄Π΅Π½ΠΈΠ΅
- Π¦Π΅Π»Ρ ΠΊΡΡΡΠΎΠ²ΠΎΠ³ΠΎ ΠΏΡΠΎΠ΅ΠΊΡΠΈΡΠΎΠ²Π°Π½ΠΈΡ ΡΠΎΡΡΠΎΠΈΡ Π² ΠΏΡΠΈΠΎΠ±ΡΠ΅ΡΠ΅Π½ΠΈΠΈ Π½Π°Π²ΡΠΊΠΎΠ² ΠΏΡΠΈΠΌΠ΅Π½Π΅Π½ΠΈΡ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΠΌΠΎΠ΄Π΅Π»ΠΈΡΠΎΠ²Π°Π½ΠΈΡ, ΠΌΠ΅ΡΠΎΠ΄ΠΎΠ² Π²ΡΡΠΈΡΠ»ΠΈΡΠ΅Π»ΡΠ½ΠΎΠΉ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΠΊΠΈ, ΡΠ΅ΠΎΡΠΈΠΈ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ ΠΈ ΡΡΠ΅Π΄ΡΡΠ² Π²ΡΡΠΈΡΠ»ΠΈΡΠ΅Π»ΡΠ½ΠΎΠΉ ΡΠ΅Ρ Π½ΠΈΠΊΠΈ ΠΏΡΠΈ Π°Π½Π°Π»ΠΈΠ·Π΅ ΠΈ ΠΏΡΠΎΠ΅ΠΊΡΠΈΡΠΎΠ²Π°Π½ΠΈΠΈ ΡΠ»Π΅ΠΊΡΡΠΈΡΠ΅ΡΠΊΠΈΡ ΡΠ΅ΠΏΠ΅ΠΉ.
- 1. ΠΠ°Π΄Π°Π½ΠΈΠ΅ Π½Π° ΠΊΡΡΡΠΎΠ²ΠΎΠΉ ΠΏΡΠΎΠ΅ΠΊΡ
- ΠΡΡ ΠΎΠ΄Π½ΡΠ΅ Π΄Π°Π½Π½ΡΠ΅
- E (t) = E0+ECcos (2Ρt/T)+ESsin (2Ρt/T)
- ΠΠ½Π°ΡΠ΅Π½ΠΈΡ Π½Π΅ Π²Π°ΡΡΠΈΡΡΠ΅ΠΌΡΡ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ²:
- Π0=110 Π, ΠΠ‘=50 Π, ΠS=50 Π, R1=R2=R3=1ΠΠΌ, R4=R5=100ΠΠΌ, Π’=0,01Ρ,
- ΠΠΈΠ°ΠΏΠ°Π·ΠΎΠ½Ρ Π²Π°ΡΡΠΈΡΡΠ΅ΠΌΡΡ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ²:
- L [10-4; 10-1](ΠΠ½), C [10-6; 10-2](Π€).
- ΠΡΠΊΠΎΠΌΠ°Ρ Π²ΡΡ ΠΎΠ΄Π½Π°Ρ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½Π°Ρ UR5
- ΠΠ΅ΡΠΎΠ΄ ΠΈΠ½ΡΠ΅Π³ΡΠΈΡΠΎΠ²Π°Π½ΠΈΡ: ΡΡΠ°ΠΏΠ΅ΡΠΈΠΉ
- ΠΠ΅ΡΠΎΠ΄ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ: ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΠ½ΠΎΠ³ΠΎ ΡΠΏΡΡΠΊΠ°
- ΠΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ ΠΡ: Π½Π΅ΠΏΠΎΡΡΠ΅Π΄ΡΡΠ²Π΅Π½Π½ΡΠΉ
- ΠΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΠ΅ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ: IL, UC.
- 2. Π€ΠΎΡΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ
- ΠΠ»Ρ ΡΠΎΡΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΡ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ Π·Π°Π΄Π°Π½Π½ΠΎΠΉ ΡΠ»Π΅ΠΊΡΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΡΠ΅ΠΏΠΈ Π²ΠΎΡΠΏΠΎΠ»ΡΠ·ΡΠ΅ΠΌΡΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ ΡΠΎΡΡΠΎΡΠ½ΠΈΠΉ.
- ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΌΠΎΠ΄Π΅Π»Ρ ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»ΡΠ΅Ρ ΡΠΎΠ±ΠΎΠΉ ΡΠΈΡΡΠ΅ΠΌΡ Π΄Π²ΡΡ Π΄ΠΈΡΡΠ΅ΡΠ΅Π½ΡΠΈΠ°Π»ΡΠ½ΡΡ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠΉ ΠΏΠ΅ΡΠ²ΠΎΠ³ΠΎ ΠΏΠΎΡΡΠ΄ΠΊΠ° Ρ ΠΏΠΎΡΡΠΎΡΠ½Π½ΡΠΌΠΈ ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΠ°ΠΌΠΈ
- Π³Π΄Π΅ — fi® — ΡΡΠ½ΠΊΡΠΈΠΈ ΠΎΡ ΡΠΎΠΏΡΠΎΡΠΈΠ²Π»Π΅Π½ΠΈΠΉ
- ΠΡΠ΅Π΄ΡΡΠ°Π²ΠΈΠΌ Π΄Π°Π½Π½ΡΡ ΡΠΈΡΡΠ΅ΠΌΡ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠΉ Π² ΠΌΠ°ΡΡΠΈΡΠ½ΠΎΠΉ ΡΠΎΡΠΌΠ΅:
- ΠΈΠ»ΠΈ, ΡΡΠΈΡΡΠ²Π°Ρ Ρ Π°ΡΠ°ΠΊΡΠ΅Ρ Π·Π°Π΄Π°ΡΡΠ΅Π³ΠΎ Π½Π°ΠΏΡΡΠΆΠ΅Π½ΠΈΡ:
- Π³Π΄Π΅ X (t) = [xt(t), x2(t)] = [iL(t), uC(t)] - Π²Π΅ΠΊΡΠΎΡ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ;
- — ΠΊΠ²Π°Π΄ΡΠ°ΡΠ½Π°Ρ ΠΌΠ°ΡΡΠΈΡΠ°
- Π²ΡΠΎΡΠΎΠ³ΠΎ ΠΏΠΎΡΡΠ΄ΠΊΠ° ΠΏΠΎΡΡΠΎΡΠ½Π½ΡΡ ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΠΎΠ² ΠΏΡΠΈ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ Π² ΡΡΠ°Π²Π½Π΅Π½ΠΈΡ ΠΌΠΎΠ΄Π΅Π»ΠΈ;
- — Π²Π΅ΠΊΡΠΎΡΡ Π²Π½Π΅ΡΠ½Π΅Π³ΠΎ Π²ΠΎΠ·Π΄Π΅ΠΉΡΡΠ²ΠΈΡ
- ΠΡΡΠ°ΠΆΠ΅Π½ΠΈΠ΅ Π΄Π»Ρ ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΡ Π·Π°Π΄Π°Π½Π½ΠΎΠΉ Π²ΡΡ ΠΎΠ΄Π½ΠΎΠΉ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΠΎΠΉ ΡΠ΅ΡΠ΅Π· ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΠ΅ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ ΠΈΠΌΠ΅Π΅Ρ Π²ΠΈΠ΄:
- y = UR5 = f1®*iL + f2®*uC
- ΠΠ°ΠΏΠΈΡΠ΅ΠΌ ΡΠΈΡΡΠ΅ΠΌΡ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠΉ ΠΎΠΏΠΈΡΡΠ²Π°Π΅ΠΌΡΡ ΡΡ Π΅ΠΌΡ
- ΠΡΡ ΠΎΠ΄Π½Π°Ρ ΡΠΈΡΡΠ΅ΠΌΠ° ΡΡΠ°Π²Π½Π΅Π½ΠΈΠΉ ΠΏΠΎ Π·Π°ΠΊΠΎΠ½Π°ΠΌ ΠΠΈΡΡ Π³ΠΎΡΠ°
- ΠΠΎΠ»ΠΈΡΠ΅ΡΡΠ²ΠΎ ΡΠ·Π»ΠΎΠ²: NΡΠ·Π» = 4
- ΠΠΎΠ»ΠΈΡΠ΅ΡΡΠ²ΠΎ Π²Π΅ΡΠ²Π΅ΠΉ: NΠ² = 6
- ΠΠΎΠ»ΠΈΡΠ΅ΡΡΠ²ΠΎ Π²Π΅ΡΠ²Π΅ΠΉ, ΡΠΎΠ΄Π΅ΡΠΆΠ°ΡΠΈΡ ΡΠΎΠ»ΡΠΊΠΎ ΠΈΠ΄Π΅Π°Π»ΡΠ½ΡΠ΅ ΠΈΡΡΠΎΡΠ½ΠΈΠΊΠΈ ΠΠΠ‘: NΠ=1,
- ΠΠΎΠ»ΠΈΡΠ΅ΡΡΠ²ΠΎ Π²Π΅ΡΠ²Π΅ΠΉ, ΡΠΎΠ΄Π΅ΡΠΆΠ°ΡΠΈΡ ΠΈΡΡΠΎΡΠ½ΠΈΠΊΠΈ ΡΠΎΠΊΠ°: NJ=0,
- ΠΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΠΌΠΈ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ ΡΠ²Π»ΡΡΡΡΡ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΠ΅ iL ΠΈ uC.
- ΠΠΎ ΠΏΠ΅ΡΠ²ΠΎΠΌΡ Π·Π°ΠΊΠΎΠ½Ρ ΠΠΈΡΡ Π³ΠΎΡΠ°
- N1ΡΡ= NΡΠ·Π» -NE -1 = 4 — 1 — 1 = 2
- ΠΠΎ Π²ΡΠΎΡΠΎΠΌΡ Π·Π°ΠΊΠΎΠ½Ρ ΠΠΈΡΡ Π³ΠΎΡΠ°
- N2ΡΡ= NΠ² -NΡΠ·Π» +1 = 6 — 4 + 1 = 3
- ΠΠΎΠΌΠΏΠΎΠ½Π΅Π½ΡΠ½ΡΠ΅ ΡΠΎΠΎΡΠ½ΠΎΡΠ΅Π½ΠΈΡ Π΄Π»Ρ ΡΠ»Π΅ΠΌΠ΅Π½ΡΠΎΠ² ΡΡ Π΅ΠΌΡ:
- ΠΠ°ΠΉΠ΄Π΅ΠΌ iR5
- iR5 = iR3 — iC
- iC = iL — iR2
- ΠΠΎΠ΄ΡΡΠ°Π²ΠΈΠΌ iC ΠΈ ΠΏΠΎΠ»ΡΡΠΈΠΌ iR5 = iR3 — iL + iR2
- Π‘ Π΄ΡΡΠ³ΠΎΠΉ ΡΡΠΎΡΠΎΠ½Ρ iR5 ΠΌΠΎΠΆΠ½ΠΎ Π½Π°ΠΉΡΠΈ ΡΠ΅ΡΠ΅Π· ΡΠΎΠΎΡΠ½ΠΎΡΠ΅Π½ΠΈΠ΅
- iR4 = iR2 + iR3
- iR5 = iR4 — iL
- UR4 = -UR3 — UR5 = -iR3R3 — iR5R5
- (*)
- Π ΡΡΠ°Π²Π½Π΅Π½ΠΈΠΈ Π½Π΅ ΠΈΠ·Π²Π΅ΡΡΠ½ΠΎ iR3 — Π½Π°ΠΉΠ΄Π΅ΠΌ Π΅Π³ΠΎ.
- iC = iR3 — iR5
- UR1 = UR2 — UR3 — UC
- ΠΡΠΊΡΠ΄Π° Π²ΡΡΠ°Π·ΠΈΠΌ iC
- =>(1)
- Π ΡΡΠΎΠΌ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠΈ Π½Π΅ΠΈΠ·Π²Π΅ΡΡΠ½ΠΎ iR2
- iR2 = iR3 — iR4
- UR4= -UR5 — UR3 = -iR5R5 — iR3R3
- (2)
- ΠΏΠΎΠ΄ΡΡΠ°Π²ΠΈΠΌ (2) Π² (1)
- (3)
- ΠΏΠΎΠ΄ΡΡΠ°Π²ΠΈΠΌ (3) Π² (*)
- ΠΠ±ΠΎΠ·Π½Π°ΡΠΈΠΌ ΡΠ΅ΡΠ΅Π· Z1
- ΠΠ±ΠΎΠ·Π½Π°ΡΠΈΠΌ ΡΠ΅ΡΠ΅Π· Z2
- ΠΠ°ΠΉΠ΄Π΅ΠΌ iC
- iC = iR3 — iR5 =
- ΠΠ°ΠΉΠ΄Π΅ΠΌ iL
- UL = E — UR4 — UR2 = E + UR3 + UR5 — UR2
- UR4 = - UR3 — UR5
- ΠΠ±ΠΎΠ·Π½Π°ΡΠΈΠΌ ΡΡΠ½ΠΊΡΠΈΡ ΠΎΡ R ΠΏΠ΅ΡΠ΅Π΄ UC ΠΊΠ°ΠΊ Π1 ΠΈ Π2, Π° ΠΏΠ΅ΡΠ΅Π΄ iL ΠΊΠ°ΠΊ Π3 ΠΈ Π4. Π’ΠΎΠ³Π΄Π°:
- ΠΠ΄Π΅
- ΠΡΠ΅Π΄ΡΡΠ°Π²ΠΈΠΌ Π΄Π°Π½Π½ΡΡ ΡΠΈΡΡΠ΅ΠΌΡ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠΉ Π² ΠΌΠ°ΡΡΠΈΡΠ½ΠΎΠΉ ΡΠΎΡΠΌΠ΅:
- ΠΡΠ΅Π΄ΡΡΠ°Π²ΠΈΠΌ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΡΡ ΠΌΠΎΠ΄Π΅Π»Ρ Π² ΠΌΠ°ΡΡΠΈΡΠ½ΠΎΠΉ ΡΠΎΡΠΌΠ΅
- ,
- ΠΈΠ»ΠΈ ΡΡΠΈΡΡΠ²Π°Ρ Ρ Π°ΡΠ°ΠΊΡΠ΅Ρ Π·Π°Π΄Π°ΡΡΠ΅Π³ΠΎ Π½Π°ΠΏΡΡΠΆΠ΅Π½ΠΈΡ:
- Π³Π΄Π΅Π²Π΅ΠΊΡΠΎΡ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ;
- — ΠΊΠ²Π°Π΄ΡΠ°ΡΠ½Π°Ρ ΠΌΠ°ΡΡΠΈΡΠ° Π²ΡΠΎΡΠΎΠ³ΠΎ ΠΏΠΎΡΡΠ΄ΠΊΠ° ΠΏΠΎΡΡΠΎΡΠ½Π½ΡΡ ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΠΎΠ² ΠΏΡΠΈ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ Π² ΡΡΠ°Π²Π½Π΅Π½ΠΈΡΡ ΠΌΠΎΠ΄Π΅Π»ΠΈ;
- — Π²Π΅ΠΊΡΠΎΡΡ Π²Π½Π΅ΡΠ½Π΅Π³ΠΎ Π²ΠΎΠ·Π΄Π΅ΠΉΡΡΠ²ΠΈΡ.
- ΠΡΡΠ°ΠΆΠ΅Π½ΠΈΠ΅ Π΄Π»Ρ ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΡ Π·Π°Π΄Π°Π½Π½ΠΎΠΉ Π²ΡΡ ΠΎΠ΄Π½ΠΎΠΉ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΠΎΠΉ ΡΠ΅ΡΠ΅Π· ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΠ΅ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ ΠΈΠΌΠ΅Π΅Ρ Π²ΠΈΠ΄:
- 3. ΠΠΎΠΈΡΠΊ ΠΏΠ΅ΡΠΈΠΎΠ΄ΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ
- ΠΡΠΈ ΡΠ΅Π°Π»ΠΈΠ·Π°ΡΠΈΠΈ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ ΠΌΠΎΠ΄Π΅Π»Π΅ΠΉ ΡΠ»Π΅ΠΊΡΡΠΈΡΠ΅ΡΠΊΠΈΡ ΡΠ΅ΠΏΠ΅ΠΉ ΠΏΡΠΎΠ±Π»Π΅ΠΌΠ° ΡΡΠ°ΡΠΈΠΎΠ½Π°ΡΠ½ΡΡ ΡΠ΅ΡΠ΅Π½ΠΈΠΉ ΡΠΈΡΡΠ΅ΠΌ Π΄ΠΈΡΡΠ΅ΡΠ΅Π½ΡΠΈΠ°Π»ΡΠ½ΡΡ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠΉ ΠΎΠ΄Π½Π° ΠΈΠ· ΡΠ΅Π½ΡΡΠ°Π»ΡΠ½ΡΡ . ΠΡΠΎ ΠΎΠ±ΡΡΡΠ½ΡΠ΅ΡΡΡ ΡΠ΅ΠΌ, ΡΡΠΎ ΡΡΠ°ΡΠΈΠΎΠ½Π°ΡΠ½ΡΠ΅ ΡΠ΅ΡΠ΅Π½ΠΈΡ ΠΎΠΏΠΈΡΡΠ²Π°ΡΡ ΡΡΡΠ°Π½ΠΎΠ²ΠΈΠ²ΡΠΈΠ΅ΡΡ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ, ΠΊΠΎΡΠΎΡΡΠ΅ Π΄Π»Ρ ΡΠ»Π΅ΠΊΡΡΠΈΡΠ΅ΡΠΊΠΈΡ ΡΠ΅ΠΏΠ΅ΠΉ ΡΠ²Π»ΡΡΡΡΡ ΠΎΠ±ΡΡΠ½ΡΠΌΠΈ ΡΠΎΡΡΠΎΡΠ½ΠΈΡΠΌΠΈ ΡΡΠ½ΠΊΡΠΈΠΎΠ½ΠΈΡΠΎΠ²Π°Π½ΠΈΡ.
- ΠΠ½Π°Π»ΠΈΡΠΈΡΠ΅ΡΠΊΠΈΠΉ ΠΌΠ΅ΡΠΎΠ΄
- ΠΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ ΠΏΠ΅ΡΠΈΠΎΠ΄ΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ Π°Π½Π°Π»ΠΈΡΠΈΡΠ΅ΡΠΊΠΈΠΌ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΠΏΡΠΎΠΈΠ·Π²ΠΎΠ΄ΠΈΡΡΡ Ρ ΠΏΠΎΠΌΠΎΡΡΡ ΠΏΡΠΎΡΠ΅Π΄ΡΡ procedure Def_An, procedure Def_An_t.
- procedure Def_An_t (t:real);
- var Y: Vector;
- begin
- Y:=Dc;
- Mult_Vect_Sc (n, Y, cos (w*t));
- Add_Vect (n, D0, Y, D);
- Y:=Ds;
- Mult_Vect_Sc (n, Y, sin (w*t));
- Add_Vect (n, D, Y, D);
- end;
- ΠΠ΅ΠΏΠΎΡΡΠ΅Π΄ΡΡΠ²Π΅Π½Π½ΡΠΉ ΠΌΠ΅ΡΠΎΠ΄
- ΠΠ»Ρ ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΡ ΠΏΠ΅ΡΠΈΠΎΠ΄ΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ Π½Π΅ΠΏΠΎΡΡΠ΅Π΄ΡΡΠ²Π΅Π½Π½ΡΠΌ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΠΏΡΠΈΠΌΠ΅Π½ΡΠ΅ΡΡΡ ΠΏΡΠΎΡΠ΅Π΄ΡΡΠ° procedure Def_Nep.
- Procedure Def_nep (var X: Vector);
- var Z1, Z2:Matrix;
- Y:Vector;
- begin
- y[1]: =0;y[2]:=0;
- for i:=1 to M do
- begin
- t:=i*h;
- Def_B_t (t);
- Rung_Kut4(A, B, h, X);
- end;
- Det_L (A, Lam, def);
- M_Exp (n, A, TT, Z1);
- Mult_Matr_Sc (n, Z1,-1);
- E_Matr (n, Z2);
- Add_Matr (n, Z1, Z2,Z2);
- Obr_Matr (Z2,Z1);
- Mult_Matr_Vect (n, Z1, Y, X);
- end;
- ΠΠ΅ΡΠΎΠ΄ ΡΡΡΠ°Π½ΠΎΠ²Π»Π΅Π½ΠΈΡ
- ΠΡΠΎΡΠ΅Π΄ΡΡΠ° procedure Ust ΡΠ΅Π°Π»ΠΈΠ·ΡΠ΅Ρ ΠΏΠΎΠΈΡΠΊ ΠΏΠ΅ΡΠΈΠΎΠ΄ΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΡΡΡΠ°Π½ΠΎΠ²Π»Π΅Π½ΠΈΡ ΠΈ ΠΏΠΎΡΡΡΠΎΠ΅Π½ΠΈΠ΅ Π²ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ Π·Π°Π²ΠΈΡΠΈΠΌΠΎΡΡΠ΅ΠΉ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ ΡΠΎΡΡΠΎΡΠ½ΠΈΡ ΠΈ Π²ΡΡ ΠΎΠ΄Π½ΠΎΠΉ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΠΎΠΉ Π² ΠΏΠ΅ΡΠ΅Ρ ΠΎΠ΄Π½ΠΎΠΌ ΡΠ΅ΠΆΠΈΠΌΠ΅.
- Procedure Ust (P:Vect);
- begin
- Def_A_B (P);
- x[1]: =0;x[2]:=0;
- j:=0;
- repeat
- j:=j+1;
- X1:=X;
- for i:=1 to M do
- begin
- t:=i*h;
- Def_B_t (t);
- Rung_Kut4(A, B, h, X);
- end;
- def:=true;
- for i:=1 to n do
- if abs (x[i]-x1[i])-delt*abs (x[i]+x1[i])/2>=0 then def:=false;
- until def;
- Init;
- kx:=(GetMaxX-X0−10)/TT/j;
- ky1:=120/50;
- ky2:=120/200;
- ky3:=-100/2;
- x[1]: =0;x[2]:=0;
- init;
- xx1:=X0;y11:=Y01;y21:=Y02;y31:=Y03;
- for k:=0 to j-1 do
- for i:=1 to M do
- begin
- t:=i*h+k*TT;
- Def_B_t (t);
- Rung_Kut4(A, B, h, X);
- xx:=X0+round (kx*t);
- y1:=Y01-round (ky1*x[1]);
- y2:=Y02-round (ky2*x[2]);
- y3:=Y03-round (ky3*U5(X));
- setcolor (Red);
- line (xx1,y11,xx, y1);
- line (xx1,y21,xx, y2);
- line (xx1,y31,xx, y3);
- xx1:=xx;y11:=y1;y21:=y2;y31:=y3;
- end;
- readln;
- CloseGraph;
- end;
4. ΠΠ°ΡΠ°ΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΠ»Π΅ΠΊΡΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΡΠ΅ΠΏΠΈ
Π¦Π΅Π»Ρ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠΈΠ½ΡΠ΅Π·Π° — ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ ΡΠΈΡΠ»ΠΎΠ²ΡΡ Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ² ΡΠ»Π΅ΠΌΠ΅Π½ΡΠΎΠ². Π‘ΠΈΠ½ΡΠ΅Π· Π½ΠΎΡΠΈΡ Π½Π°Π·Π²Π°Π½ΠΈΠ΅ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ, Π΅ΡΠ»ΠΈ ΠΎΠΏΡΠ΅Π΄Π΅Π»ΡΡΡΡΡ Π½Π°ΠΈΠ»ΡΡΡΠΈΠ΅ Π² Π·Π°Π΄Π°Π½Π½ΠΎΠΌ ΡΠΌΡΡΠ»Π΅ ΡΡΡΡΠΊΡΡΡΡ ΠΈ Π·Π½Π°ΡΠ΅Π½ΠΈΡ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ². ΠΠ°Π΄Π°ΡΡ Π²ΡΠ±ΠΎΡΠ° ΠΎΠΏΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠΉ ΡΡΡΡΠΊΡΡΡΡ Π½Π°Π·ΡΠ²Π°ΡΡ ΡΡΡΡΠΊΡΡΡΠ½ΠΎΠΉ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠ΅ΠΉ, Π° ΡΠ°ΡΡΠ΅Ρ ΠΎΠΏΡΠΈΠΌΠ°Π»ΡΠ½ΡΡ Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ² ΠΏΡΠΈ Π·Π°Π΄Π°Π½Π½ΠΎΠΉ ΡΡΡΡΠΊΡΡΡΠ΅ — ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠ΅ΠΉ.
Π ΠΎΠ±ΡΠ΅ΠΌ ΡΠ»ΡΡΠ°Π΅ ΡΠ΅ΡΠ΅Π½ΠΈΠ΅ Π·Π°Π΄Π°ΡΠΈ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ ΠΎΡΡΡΠ΅ΡΡΠ²Π»ΡΠ΅ΡΡΡ ΠΏΡΡΠ΅ΠΌ ΠΏΠ΅ΡΠ΅Π±ΠΎΡΠ° ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½Π½ΡΠΌ ΠΎΠ±ΡΠ°Π·ΠΎΠΌ Π²ΡΠ±ΠΈΡΠ°Π΅ΠΌΡΡ Π²Π°ΡΠΈΠ°Π½ΡΠΎΠ² Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ², ΡΡΠ°Π²Π½Π΅Π½ΠΈΡ ΠΈΡ ΠΌΠ΅ΠΆΠ΄Ρ ΡΠΎΠ±ΠΎΠΉ ΠΈ Π²ΡΠ±ΠΎΡΠ° Π½Π°ΠΈΠ»ΡΡΡΠ΅Π³ΠΎ Π²Π°ΡΠΈΠ°Π½ΡΠ°. ΠΠ»Π³ΠΎΡΠΈΡΠΌ Π²ΡΠ±ΠΎΡΠ° ΠΎΡΠ΅ΡΠ΅Π΄Π½ΠΎΠ³ΠΎ Π²Π°ΡΠΈΠ°Π½ΡΠ° Π·Π½Π°ΡΠ΅Π½ΠΈΠΉ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ² Π½ΠΎΡΠΈΡ Π½Π°Π·Π²Π°Π½ΠΈΠ΅ ΡΡΡΠ°ΡΠ΅Π³ΠΈΠΈ ΠΏΠΎΠΈΡΠΊΠ°.
Π§Π°ΡΠ΅ Π²ΡΠ΅Π³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΠ΅ Π·Π°Π΄Π°ΡΠΈ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ ΠΏΡΠΎΠ²ΠΎΠ΄ΠΈΡΡΡ Π² Π½Π΅ΡΠΊΠΎΠ»ΡΠΊΠΎ ΡΡΠ°ΠΏΠΎΠ². ΠΠ° ΠΏΠ΅ΡΠ²ΠΎΠΌ ΡΡΠ°ΠΏΠ΅ Π²ΡΠ΄Π΅Π»ΡΡΡ ΠΎΠ±Π»Π°ΡΡΠΈ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ², Π² ΠΊΠΎΡΠΎΡΡΡ ΡΠ΅Π»Π΅Π²Π°Ρ ΡΡΠ½ΠΊΡΠΈΡ ΡΠ²Π»ΡΠ΅ΡΡΡ ΡΠ½ΠΈΠΌΠΎΠ΄Π°Π»ΡΠ½ΠΎΠΉ, Π½Π° Π²ΡΠΎΡΠΎΠΌ — ΡΡΠΎΡΠ½ΡΡΡ ΠΏΠΎΠ»ΠΎΠΆΠ΅Π½ΠΈΡ ΡΠΎΡΠ΅ΠΊ Π»ΠΎΠΊΠ°Π»ΡΠ½ΠΎΠ³ΠΎ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ° Π² ΠΎΠ±Π»Π°ΡΡΡΡ ΡΠ½ΠΈΠΌΠΎΠ΄Π°Π»ΡΠ½ΠΎΡΡΠΈ, Π° Π·Π°ΡΠ΅ΠΌ ΡΡΠ΅Π΄ΠΈ ΡΠΎΡΠ΅ΠΊ Π»ΠΎΠΊΠ°Π»ΡΠ½ΠΎΠ³ΠΎ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ° Π²ΡΠ±ΠΈΡΠ°ΡΡ ΡΠΎΡΠΊΡ Π³Π»ΠΎΠ±Π°Π»ΡΠ½ΠΎΠ³ΠΎ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ°. ΠΡΠΈ ΡΡΠΎΠΌ Π½Π°ΠΈΠ±ΠΎΠ»Π΅Π΅ ΡΠ»ΠΎΠΆΠ½ΡΠΌ ΠΈ Π½Π΅ ΠΏΠΎΠ΄Π΄Π°ΡΡΠΈΠΌΡΡ ΡΡΡΠΎΠ³ΠΎΠΉ Π°Π»Π³ΠΎΡΠΈΡΠΌΠΈΠ·Π°ΡΠΈΠΈ ΠΈ ΡΠΎΡΠΌΠ°Π»ΠΈΠ·Π°ΡΠΈΠΈ ΡΠ²Π»ΡΠ΅ΡΡΡ ΠΏΠ΅ΡΠ²ΡΠΉ ΡΡΠ°ΠΏ. ΠΠ»Ρ ΡΠ΅Π°Π»ΠΈΠ·Π°ΡΠΈΠΈ Π²ΡΠΎΡΠΎΠ³ΠΎ ΡΡΠ°ΠΏΠ° ΠΏΡΠΈΠΌΠ΅Π½ΡΡΡ ΡΠΏΠ΅ΡΠΈΠ°Π»ΡΠ½ΡΠ΅ ΠΈΡΠ΅ΡΠ°ΡΠΈΠΎΠ½Π½ΡΠ΅ ΡΠΈΡΠ»Π΅Π½Π½ΡΠ΅ ΠΌΠ΅ΡΠΎΠ΄Ρ ΠΏΠΎΠΈΡΠΊΠ° ΡΠΊΡΡΡΠ΅ΠΌΡΠΌΠΎΠ² ΡΠ½ΠΈΠΌΠΎΠ΄Π°Π»ΡΠ½ΡΡ ΡΡΠ½ΠΊΡΠΈΠΉ.
ΠΠ°ΡΠ°ΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΠΎΡΡΡΠ΅ΡΡΠ²Π»ΡΠ΅ΡΡΡ Π² ΠΏΡΠΎΡΡΡΠ°Π½ΡΡΠ²Π΅ Π΄Π²ΡΡ Π²Π°ΡΡΠΈΡΡΠ΅ΠΌΡΡ ΠΏΠ°ΡΠ°ΠΌΠ΅ΡΡΠΎΠ² ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΠ½ΠΎΠ³ΠΎ ΡΠΏΡΡΠΊΠ° Π² ΡΠΎΠΎΡΠ²Π΅ΡΡΡΠ²ΠΈΠΈ Ρ Π°Π»Π³ΠΎΡΠΈΡΠΌΠΎΠΌ:
Procedure Koord (var Pmin: Vect; var Kpmin: real);
var i: byte;
P, PL: Vect;
def:boolean;
begin
for i:=1 to s do
p[i]: =(pl[i]+pp[i])/2;
repeat
PL:=P;
for i:=1 to s do
begin
p[i]: =Gold (P, i);
end;
def:=true;
for i:=1 to s do
if abs (p[i]-pl[i])-delt*abs (p[i]+pl[i])/2>=0 then def:=false;
until def;
for i:=1 to s do
pmin[i]: =(p[i]+pl[i])/2;
Kpmin:=Kp_nep (Pmin);
end;
5. Π Π΅Π·ΡΠ»ΡΡΠ°ΡΡ ΡΠ°Π±ΠΎΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ
Π‘ΠΏΠΈΡΠΎΠΊ ΠΈΡΠΏΠΎΠ»ΡΠ·ΡΠ΅ΠΌΠΎΠΉ Π»ΠΈΡΠ΅ΡΠ°ΡΡΡΡ
Π§Π΅ΡΠ²Π΅ΡΠ³ΠΎΠ² Π. Π. ΠΠ½ΡΠΎΡΠΌΠ°ΡΠΈΠΊΠ°. Π§Π°ΡΡΡ 3. Π£ΡΠ΅Π±Π½ΠΎΠ΅ ΠΏΠΎΡΠΎΠ±ΠΈΠ΅ Π΄Π»Ρ ΡΡΡΠ΄Π΅Π½ΡΠΎΠ² ΡΠΏΠ΅ΡΠΈΠ°Π»ΡΠ½ΠΎΡΡΠΈ «ΠΡΠΎΠΌΡΡΠ»Π΅Π½Π½Π°Ρ ΡΠ»Π΅ΠΊΡΡΠΎΠ½ΠΈΠΊΠ°» Π΄ΠΈΡΡΠ°Π½ΡΠΈΠΎΠ½Π½ΠΎΠΉ ΡΠΎΡΠΌΡ ΠΎΠ±ΡΡΠ΅Π½ΠΈΡ. — Π’ΠΎΠΌΡΠΊ: Π’Π£Π‘Π£Π 1998. — 55Ρ.
ΠΠΈΡΡΠΈΠ½Π³ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ
program kurs;
Uses Crt, Graph;
Const N=2; E0=110;Ec=50;Es=50; L=0.0001; C=0.1; R1=1; R2=1; R3=R2; R4=100;R5=R4;delt=0.01;TT=0.01;w=2*pi/TT;M=100;h=TT/M;s=2; Pl: array[1.s] of real=(1e-4,1e-6);
Pp:array[1.s] of real=(1e-1,1e-2); Type Matrix=array[1.N, 1. N] of real;
Vector=array[1.N] of real;Vect=array[1.s] of real;
Var A: Matrix;
B0,Bc, Bs, B, D0,Dc, Ds, D, Lam, X, X1:Vector;
P:Vect;
t, kpmax, Z1, Z2,K1,K2,K3,K4:real;
grdr, grmd, i, j, k: integer; X0, Y01,Y02,Y03,xx, xx1, y1,y11,y2,y21,y3,y31,y4,y41,y5,y51,y6,y61:integer; kx, ky1, ky2,ky3:real;
def:boolean;
{********************ΠΡΠΎΡΠ΅Π΄ΡΡΠ° ΡΠ»ΠΎΠΆΠ΅Π½ΠΈΡ Π΄Π²ΡΡ ΠΌΠ°ΡΡΠΈΡ A_sum=A1+A2}
Procedure Add_Matr (n:byte;A1,A2:Matrix; var A_rez:Matrix);
var i, j: byte;
begin
for i:=1 to n do
for j:=1 to n do
a_rez[i, j]: =a1[i, j]+a2[i, j]
end;
{******************ΠΡΠΎΡΠ΅Π΄ΡΡΠ° ΡΠ»ΠΎΠΆΠ΅Π½ΠΈΡ Π΄Π²ΡΡ Π²Π΅ΠΊΡΠΎΡΠΎΠ² B_sum=B1+B2} Procedure Add_Vect (n:byte;B1,B2:Vector; var B_rez:Vector);
var i: byte;
begin
for i:=1 to n do
b_rez[i]: =b1[i]+b2[i]
end;
{*******************ΠΡΠΎΡΠ΅Π΄ΡΡΠ° ΡΠΌΠ½ΠΎΠΆΠ΅Π½ΠΈΡ Π΄Π²ΡΡ ΠΌΠ°ΡΡΠΈΡ A_rez=A1*A2} Procedure Mult_Matr_Matr (n:byte;A1,A2:Matrix; var A_rez:Matrix);
var i, j, k:byte;
begin
for i:=1 to n do
for j:=1 to n do
begin
a_rez[i, j]: =0;
for k:=1 to n do a_rez[i, j]: =a_rez[i, j]+a1[i, k]*a2[k, j];
end;
end;
{*******************ΠΡΠΎΡΠ΅Π΄ΡΡΠ° ΡΠΌΠ½ΠΎΠΆΠ΅Π½ΠΈΡ ΠΌΠ°ΡΡΠΈΡΡ Π½Π° Π²Π΅ΠΊΡΠΎΡ Y=A*X} Procedure Mult_Matr_Vect (n:byte;A:Matrix; X: Vector; var Y: Vector);
var i, j: byte;
begin
for i:=1 to n do
begin
y[i]: =0;
for j:=1 to n do
y[i]: =y[i]+a[i, j]*x[j];
end;
end;
{******ΠΡΠΎΡΠ΅Π΄ΡΡΠ° ΡΠΌΠ½ΠΎΠΆΠ΅Π½ΠΈΡ ΠΌΠ°ΡΡΠΈΡΡ Π½Π° Π²Π΅ΡΠ΅ΡΡΠ²Π΅Π½Π½ΠΎΠ΅ ΡΠΈΡΠ»ΠΎ A:=A*c} Procedure Mult_Matr_Sc (n:byte; var A: Matrix; Sc: real);
var i, j: byte;
begin
for i:=1 to n do
for j:=1 to n do
a[i, j]: =Sc*a[i, j];
end;
{******ΠΡΠΎΡΠ΅Π΄ΡΡΠ° ΡΠΌΠ½ΠΎΠΆΠ΅Π½ΠΈΡ Π²Π΅ΠΊΡΠΎΡΠ° Π½Π° Π²Π΅ΡΠ΅ΡΡΠ²Π΅Π½Π½ΠΎΠ΅ ΡΠΈΡΠ»ΠΎ B:=B*c} Procedure Mult_Vect_Sc (n:byte; var B: Vector; Sc: real);
var i: byte;
begin
for i:=1 to n do b[i]: =Sc*b[i];
end;
{************* Π€ΠΎΡΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΠΌΠ°ΡΡΠΈΡ Π, B0, Bc, Bs ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ }
Procedure Def_A_B (P:Vect);
begin
Z1:=1+R1*R3/(R4*R1+R4*R3-R2*R4-R3*R2);
Z2:=R4*R1+R4*R3-R2*R4-R3*R2; K1:=r4/z2*r3/z1*r1*r3/z2-r3*r4/z2+r4/z2*r3*r5/z1-r3/z2*(r4+r3)/z1*r1/z2+(r4+r3)/z2-r5/z2*r3/(r4+r5)*1/z1; K3:=r4/z1*r1*r3/z2+r4/(r4+r5)*r5/z1+r1/z2*(r3+r4)/z1+r5/(r4+r5)*1/z1; K2:=r4/z2*r3/z1*r1/z2-r4/z2-r4/z2*r3/(r4+r5)*1/z1;
K4:=r4/z1*r1/z2-r4/(r4+r5)*1/z1;
a[1,1]: =-k3/p[1];
a[1,2]:=k1/p[1];
a[2,1]:=-k4/p[2];
a[2,2]:=k2/p[2];
b0[1]:=1/p[1];b0[2]:=0;
bc[1]:=1/p[1];bc[2]:=0;
bs[1]:=1/p[1];bs[2]:=0;
end;
{******ΠΡΠΎΡΠ΅Π΄ΡΡΠ° ΡΠΎΡΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΡ Π²Π΅ΠΊΡΠΎΡΠ° B (t)=B0+Bc*cos (wt)+Bs*sin (wt)**}
Procedure Def_B_t (t:real);
begin
b[1]: =b0[1]+bc[1]*cos (w*t)+bs[1]*sin (w*t);
b[2]:=b0[2]+bc[2]*cos (w*t)+bs[2]*sin (w*t);
end;
{****** Π€ΡΠ½ΠΊΡΠΈΡ ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΡ Π²ΡΡ ΠΎΠ΄Π½ΠΎΠΉ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΠΎΠΉ***************}
function U5(X:Vector):real;
begin
U5:=-r4/(r4+r5)*r5/z1*x[1]+r4/z2*r3/(r4+r5)*r5/z1*x[2];
end;
{********ΠΡΠΎΡΠ΅Π΄ΡΡΠ° Π²ΡΡΠΈΡΠ»Π΅Π½ΠΈΡ ΠΎΠ±ΡΠ°ΡΠ½ΠΎΠΉ ΠΌΠ°ΡΡΠΈΡΡ Π²ΡΠΎΡΠΎΠ³ΠΎ ΠΏΠΎΡΡΠ΄ΠΊΠ° }
Procedure Obr_Matr (A:Matrix; var A_obr:Matrix);
var det: real;
begin
det:=a[1,1]*a[2,2]-a[1,2]*a[2,1];
A_obr[1,1]:=a[2,2]/det;
A_obr[1,2]:=-a[1,2]/det;
A_obr[2,1]:=-a[2,1]/det;
A_obr[2,2]:=a[1,1]/det;
end;
{********************ΠΡΠΎΡΠ΅Π΄ΡΡΠ° ΡΠΎΡΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΡ Π΅Π΄ΠΈΠ½ΠΈΡΠ½ΠΎΠΉ ΠΌΠ°ΡΡΠΈΡΡ E}
Procedure E_Matr (n:byte; var E: Matrix);
var i, j: byte;
begin
for i:=1 to n do
begin
for j:=1 to n do
e[i, j]: =0;
e[i, i]:=1;
end;
end;
{ΠΡΠΎΡΠ΅Π΄ΡΡΠ° Π²ΡΡΠΈΡΠ»Π΅Π½ΠΈΡ ΡΠΎΠ±ΡΡΠ²Π΅Π½Π½ΡΡ ΡΠΈΡΠ΅Π» ΠΌΠ°ΡΡΠΈΡΡ, Π Π²ΡΠΎΡΠΎΠ³ΠΎ ΠΏΠΎΡΡΠ΄ΠΊΠ° }
Procedure Det_L (A:Matrix; var L: Vector; var comp: boolean);
var p0, p1,D:real;
begin
p1:=-(a[1,1]+a[2,2]);
p0:=a[1,1]*a[2,2]-a[1,2]*a[2,1];
D:=p1*p1−4*p0;
if D>=0 then
begin
comp:=false;
l[1]: =(-p1+sqrt (D))/2;
l[2]:=(-p1-sqrt (D))/2;
end
else begin
Comp:=true;
l[1]: =-p½;
l[2]:=sqrt (-D)/2;
end
end;
{**ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ Π²Π΅ΠΊΡΠΎΡΠΎΠ² D0, Dc, Ds Π΄Π»Ρ Π°Π½Π°Π»ΠΈΡΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΠΌΠ΅ΡΠΎΠ΄Π°***} procedure Def_An;
var Z1, Z2:Matrix;
Y1,Y2:Vector;
begin
Obr_Matr (A, Z1);
Mult_Matr_Vect (n, Z1, B0,D0);
Mult_Vect_Sc (n, D0,-1);
Mult_Matr_Matr (n, A, A, Z1);
E_Matr (n, Z2);
Mult_Matr_Sc (n, Z2,4*pi*pi/TT/TT);
Add_Matr (n, Z1, Z2,Z2);
Obr_Matr (Z2,Z1);
Mult_Matr_Vect (n, A, Bc, Y1);
Y2:=Bs;
Mult_Vect_Sc (n, Y2, w);
Add_Vect (n, Y1, Y2,Y1);
Mult_Matr_Vect (n, Z1, Y1,Dc);
Mult_Vect_Sc (n, Dc,-1);
Y1:=Bc;
Mult_Vect_Sc (n, Y1, w);
Mult_Matr_Vect (n, A, Bs, Y2);
Mult_Vect_Sc (n, Y2,-1);
Add_Vect (n, Y1, Y2,Y1);
Mult_Matr_Vect (n, Z1, Y1,Ds);
end;
{**ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ Π°Π½Π°Π»ΠΈΡΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ D9t)=D0+Dc*cos (wt)+Ds*sin (wt)***} procedure Def_An_t (t:real);
varY:Vector;
begin
Y:=Dc;
Mult_Vect_Sc (n, Y, cos (w*t));
Add_Vect (n, D0, Y, D);
Y:=Ds;
Mult_Vect_Sc (n, Y, sin (w*t));
Add_Vect (n, D, Y, D);
end;
{******************** Π§ΠΈΡΠ»Π΅Π½Π½Π°Ρ ΡΡ Π΅ΠΌΠ° ΡΡΠ°ΠΏΠ΅ΡΠΈΠΈ }
Procedure Trapec (n:integer; A: Matrix; B: Vector; h: real; var X: Vector);
var Z1, Z2,Z3:Matrix;
Y1,Y2:Vector;
begin
E_Matr (n, Z1);
Mult_Matr_Sc (n, Z1,2);
Z2:=A;
Mult_Matr_Sc (n, Z2, h);
Add_Matr (n, Z1, Z2,Z3);
Mult_Matr_Vect (n, Z3, X, Y1);
Y2:=B;
Mult_Vect_Sc (n, Y2,2*h);
Add_Vect (n, Y1, Y2,Y1);
Mult_Matr_Sc (n, Z2,-1);
Add_Matr (n, Z1, Z2,Z3);
Obr_Matr (Z3,Z1);
Mult_Matr_Vect (n, Z1, Y1,X);
end;
{*******ΠΈΠ½ΠΈΡΠΈΠ°Π»ΠΈΠ·Π°ΡΠΈΡ Π³ΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΠΆΠΈΠΌΠ°******** }
Procedure Init;
begin
grdr:=VGA;
grmd:=VGAHi; {ΠΈΠ½ΠΈΡΠΈΠ°Π»ΠΈΠ·Π°ΡΠΈΡ Π³ΡΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΠΆΠΈΠΌΠ°} InitGraph (grdr, grmd,'d:bpbgi');
X0:=15;
Y01:=150;
Y02:=300;
Y03:=450;
line (X0,5,X0,GetMaxY-5);
line (X0−5,Y01,GetMaxX-5,Y01);
line (X0−5,Y02,GetMaxX-5,Y02);
line (X0−5,Y03,GetMaxX-5,Y03);
outtextxy (X0+5,5,'x1');
outtextxy (X0+5,Y01+10,'x2');
outtextxy (X0+5,Y02+10,'U4');
outtextxy (GetMaxX-10,Y01+5,'t');
outtextxy (GetMaxX-10,Y02+5,'t');
outtextxy (GetMaxX-10,Y03+5,'t');
end;
{*******ΠΡΡΠΈΡΠ»Π΅Π½ΠΈΠ΅ Π½ΠΎΡΠΌΡ ΠΌΠ°ΡΡΠΈΡΡ************************** }
function Norm (n:byte; A: Matrix):real;
var
i, j: byte;
max, sum: real;
begin
max:=abs (a[1,1]);
for j:=1 to n do
begin
sum:=0;
for i:=1 to n do
sum:=sum+abs (a[i, j]);
if sum>max then max:=sum;
end;
Norm:=max;
end;
{****************ΠΡΠΎΡΠ΅Π΄ΡΡΠ° Π²ΡΡΠΈΡΠ»Π΅Π½ΠΈΡ ΡΠΊΡΠΏΠΎΠ½Π΅Π½ΡΠΈΠ°Π»ΡΠ½ΠΎΠΉ ΠΌΠ°ΡΡΠΈΡΡ } Procedure M_Exp (n:byte; A: Matrix; t: real; var Exp_M:Matrix);
var
i, s: byte;
Z1,Z2,Z3:Matrix;
ss:real;
begin
ss:=Norm (n, A);
s:=round (ln (10*t*ss)/ln (2));
E_Matr (n, Z1);
Z2:=Z1;Z3:=Z2;
for i:=1 to s do
begin
Mult_Matr_Matr (n, A, Z2,Z3); Mult_Matr_Sc (n, Z3, t/exp (s*ln (2))/i); Add_Matr (n, Z1, Z3,Z1);
Z2:=Z3;
end;
for i:=1 to s do
begin
Mult_Matr_Matr (n, Z1, Z1,Z2);
Z1:=Z2;
end;
Exp_M:=Z1;
end;
{*ΠΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ ΠΏΠ΅ΡΠΈΠΎΠ΄ΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ Π½Π΅ΠΏΠΎΡΡΠ΅Π΄ΡΡΠ²Π΅Π½Π½ΡΠΌ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ********} Procedure Def_nep (var X: Vector);
var Z1, Z2:Matrix;
Y:Vector;
begin
y[1]: =0;y[2]:=0;
for i:=1 to M do
begin
t:=i*h;
Def_B_t (t);
Trapec (n, A, B, h, X);
end;
Det_L (A, Lam, def);
M_Exp (n, A, TT, Z1);
Mult_Matr_Sc (n, Z1,-1);
E_Matr (n, Z2);
Add_Matr (n, Z1, Z2,Z2);
Obr_Matr (Z2,Z1);
Mult_Matr_Vect (n, Z1, Y, X);
end;
{***********ΠΠΎΡΡΡΠΎΠ΅Π½ΠΈΠ΅ Π³ΡΠ°ΡΠΈΠΊΠ° ΠΏΠ΅ΡΠΈΠΎΠ΄ΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ ****************} Procedure Per_Solve (P:Vect);
begin
Def_A_B (P);
Def_An;
Def_nep (X);
Init;
kx:=(GetMaxX-X0−10)/TT;
ky1:=120/50;
ky2:=120/200;
ky3:=-100/2;
Def_An_t (0);
xx1:=X0;
y11:=Y01-round (ky1*x[1]);
y21:=Y02-round (ky2*x[2]);
y31:=Y01-round (ky1*d[1]);
y41:=Y02-round (ky2*d[2]);
y51:=Y03-round (ky3*U5(X));
y61:=Y03-round (ky3*U5(D));
for i:=1 to M do
begin
t:=i*h;
Def_An_t (t);
Def_B_t (t);
Trapec (n, A, B, h, X);
xx:=X0+round (kx*t);
y1:=Y01-round (ky1*x[1]);
y2:=Y02-round (ky2*x[2]);
y3:=Y01-round (ky1*d[1]);
y4:=Y02-round (ky2*d[2]);
y5:=Y03-round (ky3*U5(X));
y6:=Y03-round (ky3*U5(D));
setcolor (Red);
line (xx1,y11,xx, y1);
line (xx1,y21,xx, y2);
line (xx1,y51,xx, y5);
setcolor (Green);
line (xx1,y31,xx, y3);
line (xx1,y41,xx, y4);
line (xx1,y61,xx, y6); xx1:=xx;y11:=y1;y21:=y2;y31:=y3;y41:=y4;y51:=y5;y61:=y6;
end;
readln;
CloseGraph;
end;
{****************************ΠΠ΅ΡΠΎΠ΄ ΡΡΡΠ°Π½ΠΎΠ²Π»Π΅Π½ΠΈΡ****************}
Procedure Ust (P:Vect);
begin
Def_A_B (P);
x[1]: =0;x[2]:=0;
j:=0;
repeat
j:=j+1;
X1:=X;
for i:=1 to M do
begin
t:=i*h;
Def_B_t (t);
trapec (n, A, B, h, X);
end;
def:=true;
for i:=1 to n do
if abs (x[i]-x1[i])-delt*abs (x[i]+x1[i])/2>=0 then def:=false; until def;
Init;
kx:=(GetMaxX-X0−10)/TT/j;
ky1:=120/50;
ky2:=120/200;
ky3:=-100/2;
x[1]: =0;x[2]:=0;
init;
xx1:=X0;y11:=Y01;y21:=Y02;y31:=Y03;
for k:=0 to j-1 do
for i:=1 to M do
begin
t:=i*h+k*TT;
Def_B_t (t);
trapec (n, A, B, h, X);
xx:=X0+round (kx*t);
y1:=Y01-round (ky1*x[1]);
y2:=Y02-round (ky2*x[2]);
y3:=Y03-round (ky3*U5(X));
setcolor (Red);
line (xx1,y11,xx, y1);
line (xx1,y21,xx, y2);
line (xx1,y31,xx, y3);
xx1:=xx;y11:=y1;y21:=y2;y31:=y3;
end;
readln;
CloseGraph;
end;
{*ΠΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΠ° ΠΏΡΠ»ΡΡΠ°ΡΠΈΠΉ ΠΏΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ Π½Π΅ΠΏΠΎΡΡΠ΅Π΄ΡΡΠ²Π΅Π½Π½ΡΠΌ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ} function Kp_nep (P:Vect):real;
var Usr, Umax, Umin: real;
begin
Def_A_B (P);
Def_nep (X);
Usr:=0;Umin:=U5(X);Umax:=Umin;
for i:=1 to M do
begin
t:=i*h;
Def_B_t (t);
trapec (n, A, B, h, X);
Usr:=Usr+U5(X);
if U5(X)>Umax then Umax:=U5(X);
if U5(X)
Kp_nep:=(Umax-Umin)*M/Usr/2;
end;
{***************ΠΠ΅ΡΠΎΠ΄ Π·ΠΎΠ»ΠΎΡΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ }
function Gold (P:Vect;i:byte):real;
const g=0.618 034;
var r, xl, xp: real;
X1,X2:Vect;
begin
xp:=pp[i]; xl:=pl[i];X1:=P;X2:=P;
repeat r:=g*(xp-xl);x1[i]: =xl+r;x2[i]:=xp-r;
if kp_nep (X1)
until abs (xp-xl)-delt*abs (xp+xl)/2<0;
Gold:=(xp+xl)/2;
end;
{*************** ΠΠ΅ΡΠΎΠ΄ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΠ½ΠΎΠ³ΠΎ ΡΠΏΡΡΠΊΠ° } Procedure Koord (var Pmin: Vect; var Kpmin: real);
var i: byte;
P, PL: Vect;
def:boolean;
begin
for i:=1 to s do
p[i]: =(pl[i]+pp[i])/2;
repeat
PL:=P;
for i:=1 to s do
begin
p[i]: =Gold (P, i);
end;
def:=true;
for i:=1 to s do
if abs (p[i]-pl[i])-delt*abs (p[i]+pl[i])/2>=0 then def:=false; until def;
for i:=1 to s do
pmin[i]: =(p[i]+pl[i])/2;
Kpmin:=Kp_nep (Pmin);
end;
{************************ Π£ΠΏΡΠ°Π²Π»ΡΡΡΠ°Ρ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ° *********** }
Begin
ClrScr;
Koord (P, kpmax);{ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΠ½ΠΎΠ³ΠΎ ΡΠΏΡΡΠΊΠ°} Def_A_B (P);{ΡΠΎΡΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΠΌΠ°ΡΡΠΈΡ Π, Π ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΠΈ} Det_L (A, Lam, def);{ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ ΡΠΎΠ±ΡΡΠ²Π΅Π½Π½ΡΡ ΡΠΈΡΠ΅Π»}
writeln ('ΠΠΏΡΠΈΠΌΠ°Π»ΡΠ½Π°Ρ ΡΠΎΡΠΊΠ°:'); writeln ('p[1]=', p[1],'p[2]=', p[2],'kpmax=', kpmax); writeln ('Π‘ΠΎΠ±ΡΡΠ²Π΅Π½Π½ΡΠ΅ ΡΠΈΡΠ»Π°:'); {Π²ΡΠ²ΠΎΠ΄ ΡΠΎΠ±ΡΡΠ²Π΅Π½Π½ΡΡ ΡΠΈΡΠ΅Π» Π½Π° ΡΠΊΡΠ°Π½} if not def then
begin
writeln ('l[1]=', lam[1]);writeln ('l[2]=', lam[2]); writeln ('ΠΠΎΡΡΠΎΡΠ½Π½ΡΠ΅ Π²ΡΠ΅ΠΌΠ΅Π½ΠΈ:'); writeln ('tau1=', 1/abs (lam[1]));writeln ('tau2', 1/abs (lam[2])); end
else
begin
writeln ('l[1]=', lam[1],'+j', lam[2]);writeln ('l[1]=', lam[1],'-j', lam[2]); writeln ('ΠΠΎΡΡΠΎΡΠ½Π½Π°Ρ Π²ΡΠ΅ΠΌΠ΅Π½ΠΈ: tau=', 1/abs (lam[1]));
writeln ('Π‘ΠΎΠ±ΡΡΠ²Π΅Π½Π½Π°Ρ ΡΠ°ΡΡΠΎΡΠ°: fc=', lam[2]/2/pi);
end;
readln;
Ust (P);{ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ ΠΏΠ΅ΡΠΈΠΎΠ΄ΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΡΡΡΠ°Π½ΠΎΠ²Π»Π΅Π½ΠΈΡ} Per_Solve (P);{ΠΎΠΏΡΠ΅Π΄Π΅Π»Π΅Π½ΠΈΠ΅ ΠΏΠ΅ΡΠΈΠΎΠ΄ΠΈΡΠ΅ΡΠΊΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ Π½Π΅ΠΏΠΎΡΡΠ΅Π΄ΡΡΠ²Π΅Π½Π½ΡΠΌ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ}
End.