ΠΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΡΠ½ΠΊΡΠΈΠΈ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΡΠΎΠΏΡΡΠΆΠ΅Π½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ²
Π‘ΡΡΡ ΠΌΠ΅ΡΠΎΠ΄Π° Π·ΠΎΠ»ΠΎΡΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ ΡΠΎΡΡΠΎΠΈΡ Π² ΠΈΡΠΏΠΎΠ»ΡΠ·ΠΎΠ²Π°Π½ΠΈΠΈ ΠΎΠ΄Π½ΠΎΠΈΠΌΡΠ½Π½ΠΎΠΉ ΠΏΡΠΎΠΏΠΎΡΡΠΈΠΈ. Π’ΠΎΡΠΊΠ° Π½Π° ΠΎΡΡΠ΅Π·ΠΊΠ΅ ΡΠ²Π»ΡΠ΅ΡΡΡ Π΅Π³ΠΎ «Π·ΠΎΠ»ΠΎΡΡΠΌ ΡΠ΅ΡΠ΅Π½ΠΈΠ΅ΠΌ», Π΅ΡΠ»ΠΈ ΠΎΠ½Π° Π΄Π΅Π»ΠΈΡ ΠΎΡΡΠ΅Π·ΠΎΠΊ ΡΠ°ΠΊ, ΡΡΠΎ Π΅Π³ΠΎ Π±ΠΎΠ»ΡΡΠ°Ρ ΡΠ°ΡΡΡ ΡΠΎΠΎΡΠ½ΠΎΡΠΈΡΡΡ Ρ Π½ΠΈΠΌ ΡΠ°ΠΌΠΈΠΌ ΡΠ°ΠΊ ΠΆΠ΅, ΠΊΠ°ΠΊ ΠΈ ΠΌΠ΅Π½ΡΡΠ°Ρ ΡΠ°ΡΡΡ Ρ Π±ΠΎΠ»ΡΡΠ΅ΠΉ. ΠΠ° ΠΎΡΡΠ΅Π·ΠΊΠ΅ ΡΡΡΠ΅ΡΡΠ²ΡΠ΅Ρ Π΄Π²Π΅ ΡΠ°ΠΊΠΈΡ ΡΠΎΡΠΊΠΈ. ΠΠ±ΠΎΠ·Π½Π°ΡΠΈΠΌ ΠΈΡ Π± (1) ΠΈ Π± (2). ΠΡΡΡΡ a, b — Π»Π΅Π²Π°Ρ ΠΈ ΠΏΡΠ°Π²Π°Ρ Π³ΡΠ°Π½ΠΈΡΡ ΠΎΡΡΠ΅Π·ΠΊΠ° ΡΠΎΠΎΡΠ²Π΅ΡΡΡΠ²Π΅Π½Π½ΠΎ, ΡΠΎΠ³Π΄Π° ΠΈΡ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ ΠΌΠΎΠΆΠ½ΠΎ… Π§ΠΈΡΠ°ΡΡ Π΅ΡΡ >
ΠΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΡΠ½ΠΊΡΠΈΠΈ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΡΠΎΠΏΡΡΠΆΠ΅Π½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ² (ΡΠ΅ΡΠ΅ΡΠ°Ρ, ΠΊΡΡΡΠΎΠ²Π°Ρ, Π΄ΠΈΠΏΠ»ΠΎΠΌ, ΠΊΠΎΠ½ΡΡΠΎΠ»ΡΠ½Π°Ρ)
ΠΠΈΠ½ΠΈΡΡΠ΅ΡΡΡΠ²ΠΎ ΠΎΠ±ΡΠ°Π·ΠΎΠ²Π°Π½ΠΈΡ ΠΈ Π½Π°ΡΠΊΠΈ Π ΠΎΡΡΠΈΠΉΡΠΊΠΎΠΉ Π€Π΅Π΄Π΅ΡΠ°ΡΠΈΠΈ ΠΠΠ‘Π£ΠΠΠ Π‘Π’ΠΠΠΠΠΠ ΠΠΠ ΠΠΠΠΠΠ’ΠΠΠ¬ΠΠΠ Π£Π§Π ΠΠΠΠΠΠΠ ΠΠ«Π‘Π¨ΠΠΠ ΠΠ ΠΠ€ΠΠ‘Π‘ΠΠΠΠΠΠ¬ΠΠΠΠ ΠΠΠ ΠΠΠΠΠΠΠΠ―
«ΠΠ ΠΠΠΠ£Π ΠΠ‘ΠΠΠ ΠΠΠ‘Π£ΠΠΠ Π‘Π’ΠΠΠΠΠ«Π Π£ΠΠΠΠΠ Π‘ΠΠ’ΠΠ’»
Π€Π°ΠΊΡΠ»ΡΡΠ΅Ρ ΡΠΊΠΎΠ½ΠΎΠΌΠΈΠΊΠΈ ΠΈ ΡΠΏΡΠ°Π²Π»Π΅Π½ΠΈΡ ΠΠ°ΡΠ΅Π΄ΡΠ° ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠ² ΠΈ ΠΌΠΎΠ΄Π΅Π»Π΅ΠΉ Π² ΡΠΊΠΎΠ½ΠΎΠΌΠΈΠΊΠ΅ ΠΠ’Π§ΠΠ’ ΠΏΠΎ ΡΡΠ΅Π±Π½ΠΎΠΉ-Π²ΡΡΠΈΡΠ»ΠΈΡΠ΅Π»ΡΠ½ΠΎΠΉ ΠΏΡΠ°ΠΊΡΠΈΠΊΠ΅ Π½Π° Π±Π°Π·Π΅ ΠΊΠ°ΡΠ΅Π΄ΡΡ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΈΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠ² ΠΈ ΠΌΠΎΠ΄Π΅Π»Π΅ΠΉ Π² ΡΠΊΠΎΠ½ΠΎΠΌΠΈΠΊΠ΅ ΠΠΠ£ ΠΠΠ£ 80 116.65.9011.06 Π Π ΡΠΊΠΎΠ²ΠΎΠ΄ΠΈΡΠ΅Π»ΠΈ ΠΎΡ ΠΊΠ°ΡΠ΅Π΄ΡΡ:
ΠΏΡΠΎΡΠ΅ΡΡΠΎΡ, Π·Π°Π²Π΅Π΄ΡΡΡΠΈΠΉ ΠΊΠ°ΡΠ΅Π΄ΡΡ ΠΠΠΈΠΠ Π. Π. Π Π΅Π½Π½Π΅Ρ Π°ΡΡΠΈΡΡΠ΅Π½Ρ ΠΊΠ°ΡΠ΅Π΄ΡΡ ΠΠΠΈΠΠ Π . Π. Π¨Π°ΡΡ ΠΌΠ΅ΡΠΎΠ²Π° Π°ΡΡΠΈΡΡΠ΅Π½Ρ ΠΊΠ°ΡΠ΅Π΄ΡΡ ΠΠΠΈΠΠ Π’. Π. ΠΠ΅Π»Π΅Π½ΠΈΠ½Π° ΠΡΠΏΠΎΠ»Π½ΠΈΡΠ΅Π»ΠΈ:
Π‘ΡΡΠ΄Π΅Π½Ρ Π³ΡΡΠΏΠΏΡ 09 ΠΠΠ Π. Π. ΠΠ²Π΄ΠΎΠΊΠΈΠΌΠΎΠ² ΠΡΠ΅Π½Π±ΡΡΠ³ 2011
1. ΠΠ°Π΄Π°Π½ΠΈΠ΅ Π½Π° ΠΏΡΠ°ΠΊΡΠΈΠΊΡ
2. ΠΠ΅Π·ΡΡΠ»ΠΎΠ²Π½Π°Ρ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΡΠ½ΠΊΡΠΈΠΈ ΠΌΠ½ΠΎΠ³ΠΈΡ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ
2.1 ΠΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΡΠ½ΠΊΡΠΈΠΈ Π²Π΄ΠΎΠ»Ρ Π·Π°Π΄Π°Π½Π½ΠΎΠ³ΠΎ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΡ
2.2 ΠΡΠ°Π΄ΠΈΠ΅Π½ΡΠ½ΡΠΉ ΠΌΠ΅ΡΠΎΠ΄ Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ°
2.3 ΠΠ΅ΡΠΎΠ΄ ΡΠΎΠΏΡΡΠΆΠ΅Π½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ²
2.4 ΠΠ»ΠΎΠΊ-ΡΡ Π΅ΠΌΠ° Π°Π»Π³ΠΎΡΠΈΡΠΌΠ° ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ
2.5 ΠΠΈΡΡΠΈΠ½Π³ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ «ΠΠ΅Π·ΡΡΠ»ΠΎΠ²Π½Π°Ρ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΡΠ½ΠΊΡΠΈΠΈ»
1. ΠΠ°Π΄Π°Π½ΠΈΠ΅ Π½Π° ΠΏΡΠ°ΠΊΡΠΈΠΊΡ Π Π°Π·ΡΠ°Π±ΠΎΡΠ°ΡΡ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΠΎΠ΅ ΠΎΠ±Π΅ΡΠΏΠ΅ΡΠ΅Π½ΠΈΠ΅, ΡΠ΅Π°Π»ΠΈΠ·ΡΡΡΠ΅Π΅ Π½Π°Ρ ΠΎΠΆΠ΄Π΅Π½ΠΈΠ΅ ΠΌΠΈΠ½ΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ³ΠΎ Π·Π½Π°ΡΠ΅Π½ΠΈΡ Π·Π°Π΄Π°Π½Π½ΠΎΠΉ ΡΡΠ½ΠΊΡΠΈΠΈ ΠΌΠ½ΠΎΠ³ΠΈΡ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ ΠΈ Π΅Ρ ΡΠΎΡΠΊΡ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ° ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΡΠΎΠΏΡΡΠΆΡΠ½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ².
2. ΠΠ΅Π·ΡΡΠ»ΠΎΠ²Π½Π°Ρ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΡΠ½ΠΊΡΠΈΠΈ ΠΌΠ½ΠΎΠ³ΠΈΡ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ ΠΠ΅ΡΠΎΠ΄ ΡΠΎΠΏΡΡΠΆΡΠ½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ² ΡΠ²Π»ΡΠ΅ΡΡΡ Π±ΡΡΡΡΠΎΡΡ ΠΎΠ΄ΡΡΠΈΠΌΡΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ, Π½ΠΎ ΠΎΠ½ ΡΡ ΠΎΠ΄ΠΈΡΡΡ ΡΠΎΠ»ΡΠΊΠΎ ΠΈΠ· Π½Π°ΡΠ°Π»ΡΠ½ΡΡ ΡΠΎΡΠ΅ΠΊ, Π΄ΠΎΡΡΠ°ΡΠΎΡΠ½ΠΎ Π±Π»ΠΈΠ·ΠΊΠΈΡ ΠΊ ΡΠΎΡΠΊΠ΅ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ°.
ΠΠ°ΡΠ°Π»ΡΠ½ΠΎΠ΅ ΠΏΡΠΈΠ±Π»ΠΈΠΆΠ΅Π½ΠΈΠ΅ ΡΠ»Π΅Π΄ΡΠ΅Ρ Π½Π°ΠΉΡΠΈ Ρ ΠΏΠΎΠΌΠΎΡΡΡ ΠΏΡΠΎΡΠ΅Π΄ΡΡΡ, ΡΠ΅Π°Π»ΠΈΠ·ΡΡΡΠ΅ΠΉ ΠΌΠ΅ΡΠΎΠ΄ Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ° (ΠΈΡΠΏΠΎΠ»ΡΠ·ΡΡΡΠ΅ΠΉ ΠΏΡΠΎΡΠ΅Π΄ΡΡΡ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ Π²Π΄ΠΎΠ»Ρ Π·Π°Π΄Π°Π½Π½ΠΎΠ³ΠΎ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΡ). ΠΠ°ΡΠ΅ΠΌ ΡΠ»Π΅Π΄ΡΠ΅Ρ ΡΡΠΎΡΠ½ΠΈΡΡ ΡΠΎΡΠΊΡ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ° ΠΏΡΠΈ ΠΏΠΎΠΌΠΎΡΠΈ ΠΏΡΠΎΡΠ΅Π΄ΡΡΡ, ΡΠ΅Π°Π»ΠΈΠ·ΡΡΡΠ΅ΠΉ ΠΌΠ΅ΡΠΎΠ΄ ΡΠΎΠΏΡΡΠΆΡΠ½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ² (ΠΈΡΠΏΠΎΠ»ΡΠ·ΡΡΡΠ΅ΠΉ ΠΏΡΠΎΡΠ΅Π΄ΡΡΡ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ Π²Π΄ΠΎΠ»Ρ Π·Π°Π΄Π°Π½Π½ΠΎΠ³ΠΎ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΡ). ΠΡΠΈ ΡΡΠΎΠΌ Π½Π΅ΠΎΠ±Ρ ΠΎΠ΄ΠΈΠΌΠΎ ΡΡΠ΅ΡΡΡ Π²ΠΎΠ·ΠΌΠΎΠΆΠ½ΠΎΡΡΡ ΡΠΎΠ³ΠΎ, ΡΡΠΎ Π·Π½Π°ΡΠ΅Π½ΠΈΡ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ ΠΈ ΡΠ°ΠΌΠΎΠΉ ΡΡΠ½ΠΊΡΠΈΠΈ ΠΌΠΎΠ³ΡΡ ΠΈΠΌΠ΅ΡΡ ΡΠ°Π·Π½ΡΠ΅ ΠΌΠ°ΡΡΡΠ°Π±Ρ (ΠΈΡΠΏΠΎΠ»ΡΠ·ΠΎΠ²Π°ΡΡ ΠΌΠ°ΡΡΡΠ°Π±ΠΈΡΡΡΡΠΈΠ΅ ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΡ ΠΈ ΠΎΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΡΠ΅ ΠΊΡΠΈΡΠ΅ΡΠΈΠΈ ΠΎΡΡΠ°Π½ΠΎΠ²Π°).
2.1 ΠΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΡΠ½ΠΊΡΠΈΠΈ Π²Π΄ΠΎΠ»Ρ Π·Π°Π΄Π°Π½Π½ΠΎΠ³ΠΎ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΡ Π§ΠΈΡΠ»Π΅Π½Π½ΠΎΠ΅ ΡΠ΅ΡΠ΅Π½ΠΈΠ΅ Π·Π°Π΄Π°ΡΠΈ ΠΎΠ΄Π½ΠΎΠΌΠ΅ΡΠ½ΠΎΠΉ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ:
ΠΠΌΠ΅Π΅ΠΌ Π½Π΅ΠΊΠΎΡΠΎΡΡΡ Π½Π°ΡΠ°Π»ΡΠ½ΡΡ ΡΠΎΡΠΊΡ ΠΈ Π²Π΅ΠΊΡΠΎΡ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΡ. Π’ΡΠ΅Π±ΡΠ΅ΡΡΡ Ρ Π·Π°Π΄Π°Π½Π½ΠΎΠΉ ΡΠΎΡΠ½ΠΎΡΡΡΡ Π΅ Π½Π°ΠΉΡΠΈ ΡΠΎΡΠΊΡ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ° ΡΡΠ½ΠΊΡΠΈΠΈ Π² Π΄Π°Π½Π½ΠΎΠΌ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΠΈ ΠΈ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·ΠΈΡΡΡΡΠ΅Π΅ ΡΡΠ½ΠΊΡΠΈΡ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ Π±.
Π‘ ΠΏΠΎΠΌΠΎΡΡΡ ΠΌΠ΅ΡΠΎΠ΄Π° ΡΠ΄Π²ΠΎΠ΅Π½ΠΈΡ ΡΠ°Π³Π° Π»ΠΎΠΊΠ°Π»ΠΈΠ·ΡΠ΅ΠΌ ΠΈΠ½ΡΠ΅ΡΠ²Π°Π», Π½Π° ΠΊΠΎΡΠΎΡΠΎΠΌ Π½Π°Ρ ΠΎΠ΄ΠΈΡΡΡ Π±min. Π‘ΡΡΡ ΠΌΠ΅ΡΠΎΠ΄Π° ΡΠ΄Π²ΠΎΠ΅Π½ΠΈΡ ΡΠ°Π³Π° ΡΠΎΡΡΠΎΠΈΡ Π² ΡΠΎΠΌ, ΡΡΠΎ ΠΌΡ Π²ΡΠ±ΠΈΡΠ°Π΅ΠΌ Π½Π΅ΠΊΠΎΡΠΎΡΠΎΠ΅ Π½Π°ΡΠ°Π»ΡΠ½ΠΎΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ Π± (Π² Π΄Π°Π½Π½ΠΎΠΌ ΡΠ»ΡΡΠ°Π΅ Π²ΡΠ±Π΅ΡΠ΅ΠΌ Π½ΡΠ»Π΅Π²ΠΎΠ΅) ΠΈ Π½Π΅ΠΊΠΎΡΠΎΡΡΡ Π²Π΅Π»ΠΈΡΠΈΠ½Ρ ΡΠ°Π³Π° h. Π‘Π»Π΅Π΄ΡΡΡΡΡ ΡΠΎΡΠΊΡ ΠΈΡΠ΅ΡΠ°ΡΠΈΠΎΠ½Π½ΠΎΠΉ ΠΏΠΎΡΠ»Π΅Π΄ΠΎΠ²Π°ΡΠ΅Π»ΡΠ½ΠΎΡΡΠΈ Π±i Π½Π°Ρ ΠΎΠ΄ΠΈΠΌ ΠΏΠΎ ΡΠΎΡΠΌΡΠ»Π΅ Π±i=Π±i-1+2i-1*h. ΠΠ°ΡΠ΅ΠΌ Π²ΡΡΠΈΡΠ»ΡΠ΅ΠΌ Ρ i=Ρ 0+Π±ih. ΠΡΠ΅ΡΠ°ΡΠΈΠΈ ΠΏΡΠΎΠ΄ΠΎΠ»ΠΆΠ°Π΅ΠΌ Π΄ΠΎ ΡΠ΅Ρ ΠΏΠΎΡ, ΠΏΠΎΠΊΠ° Π½Π΅ Π²ΡΠΏΠΎΠ»Π½ΠΈΡΡΡ ΡΡΠ»ΠΎΠ²ΠΈΠ΅ f (Ρ i)>f (Ρ i-1). ΠΡΠΈ Π²ΡΠΏΠΎΠ»Π½Π΅Π½ΠΈΠΈ ΡΡΠΎΠ³ΠΎ ΡΡΠ»ΠΎΠ²ΠΈΡ ΠΏΠΎΠ»Π°Π³Π°Π΅ΠΌ Π±min=[Π±i-2; Π±i].
ΠΠ°ΡΠ΅ΠΌ ΡΡΠΎΡΠ½ΠΈΠΌ Π³ΡΠ°Π½ΠΈΡΡ ΠΈΠ½ΡΠ΅ΡΠ²Π°Π»Π° (ΠΈ, ΡΠΎΠΎΡΠ²Π΅ΡΡΡΠ²Π΅Π½Π½ΠΎ, ΠΏΡΠΈΠ±Π»ΠΈΠΆΡΠ½Π½ΠΎΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΠ΅ ΡΠΎΡΠΊΠΈ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ°) Ρ ΠΏΠΎΠΌΠΎΡΡΡ ΠΎΠ΄Π½ΠΎΠ³ΠΎ ΠΈΠ· Π±ΡΡΡΡΠΎΡΡ ΠΎΠ΄ΡΡΠΈΡ ΡΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠ². ΠΠ°ΠΏΡΠΈΠΌΠ΅Ρ, Ρ ΠΏΠΎΠΌΠΎΡΡΡ ΠΌΠ΅ΡΠΎΠ΄Π° Π·ΠΎΠ»ΠΎΡΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ.
Π‘ΡΡΡ ΠΌΠ΅ΡΠΎΠ΄Π° Π·ΠΎΠ»ΠΎΡΠΎΠ³ΠΎ ΡΠ΅ΡΠ΅Π½ΠΈΡ ΡΠΎΡΡΠΎΠΈΡ Π² ΠΈΡΠΏΠΎΠ»ΡΠ·ΠΎΠ²Π°Π½ΠΈΠΈ ΠΎΠ΄Π½ΠΎΠΈΠΌΡΠ½Π½ΠΎΠΉ ΠΏΡΠΎΠΏΠΎΡΡΠΈΠΈ. Π’ΠΎΡΠΊΠ° Π½Π° ΠΎΡΡΠ΅Π·ΠΊΠ΅ ΡΠ²Π»ΡΠ΅ΡΡΡ Π΅Π³ΠΎ «Π·ΠΎΠ»ΠΎΡΡΠΌ ΡΠ΅ΡΠ΅Π½ΠΈΠ΅ΠΌ», Π΅ΡΠ»ΠΈ ΠΎΠ½Π° Π΄Π΅Π»ΠΈΡ ΠΎΡΡΠ΅Π·ΠΎΠΊ ΡΠ°ΠΊ, ΡΡΠΎ Π΅Π³ΠΎ Π±ΠΎΠ»ΡΡΠ°Ρ ΡΠ°ΡΡΡ ΡΠΎΠΎΡΠ½ΠΎΡΠΈΡΡΡ Ρ Π½ΠΈΠΌ ΡΠ°ΠΌΠΈΠΌ ΡΠ°ΠΊ ΠΆΠ΅, ΠΊΠ°ΠΊ ΠΈ ΠΌΠ΅Π½ΡΡΠ°Ρ ΡΠ°ΡΡΡ Ρ Π±ΠΎΠ»ΡΡΠ΅ΠΉ. ΠΠ° ΠΎΡΡΠ΅Π·ΠΊΠ΅ ΡΡΡΠ΅ΡΡΠ²ΡΠ΅Ρ Π΄Π²Π΅ ΡΠ°ΠΊΠΈΡ ΡΠΎΡΠΊΠΈ. ΠΠ±ΠΎΠ·Π½Π°ΡΠΈΠΌ ΠΈΡ Π± (1) ΠΈ Π± (2). ΠΡΡΡΡ a, b — Π»Π΅Π²Π°Ρ ΠΈ ΠΏΡΠ°Π²Π°Ρ Π³ΡΠ°Π½ΠΈΡΡ ΠΎΡΡΠ΅Π·ΠΊΠ° ΡΠΎΠΎΡΠ²Π΅ΡΡΡΠ²Π΅Π½Π½ΠΎ, ΡΠΎΠ³Π΄Π° ΠΈΡ ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ ΠΌΠΎΠΆΠ½ΠΎ Π²ΡΡΠ°Π·ΠΈΡΡ, ΠΈΡΠΏΠΎΠ»ΡΠ·ΡΡ ΡΠΎΡΠΌΡΠ»Ρ. ΠΠ°Ρ ΠΎΠ΄Ρ Ρ (1)=Ρ 0+Π± (1)h ΠΈ Ρ (2)=Ρ 0+Π± (2)h, ΡΡΠ°Π²Π½ΠΈΠ²Π°Π΅ΠΌ f (Ρ (1)) ΠΈ f (Ρ (2)). ΠΡΠ»ΠΈ f (Ρ (1))>f (Ρ (2)), ΡΠΎ. ΠΡΠ»ΠΈ f (Ρ (1))
2.2 ΠΡΠ°Π΄ΠΈΠ΅Π½ΡΠ½ΡΠΉ ΠΌΠ΅ΡΠΎΠ΄ Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ° ΠΠΌΠ΅Π΅ΠΌ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΎΠ½Π½ΡΡ Π·Π°Π΄Π°ΡΡ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ ΡΡΠ½ΠΊΡΠΈΠΈ:, Π³Π΄Π΅. ΠΠ°ΠΆΠ΄ΠΎΠ΅ ΡΠ»Π΅Π΄ΡΡΡΠ΅Π΅ ΠΏΡΠΈΠ±Π»ΠΈΠΆΠ΅Π½ΠΈΠ΅ Π½Π°Ρ ΠΎΠ΄ΠΈΠΌ ΠΏΠΎ ΡΠΎΡΠΌΡΠ»Π΅. ΠΡΠ»ΠΈ Π²Π΅ΠΊΡΠΎΡ hk ΡΠ°Π²Π΅Π½ Π°Π½ΡΠΈΠ³ΡΠ°Π΄ΠΈΠ΅Π½ΡΡ ΡΡΠ½ΠΊΡΠΈΠΈ Π² ΡΠΎΡΠΊΠ΅ Ρ k, Π° Π±k Π²ΡΠ±ΠΈΡΠ°Π΅ΡΡΡ ΠΈΠ· ΡΡΠ»ΠΎΠ²ΠΈΡ, ΡΠΎ ΠΌΡ ΠΈΠΌΠ΅Π΅ΠΌ ΠΌΠ΅ΡΠΎΠ΄ Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ°.
ΠΠ΅ΡΠΎΠ΄ Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ° ΡΠ²Π»ΡΠ΅ΡΡΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠΌ ΠΏΠ΅ΡΠ²ΠΎΠ³ΠΎ ΠΏΠΎΡΡΠ΄ΠΊΠ°, ΡΠ°ΠΊ ΠΊΠ°ΠΊ ΠΈΡΠΏΠΎΠ»ΡΠ·ΡΠ΅Ρ ΠΈΠ½ΡΠΎΡΠΌΠ°ΡΠΈΡ ΠΎ Π·Π½Π°ΡΠ΅Π½ΠΈΠΈ ΠΏΠ΅ΡΠ²ΡΡ ΠΏΡΠΎΠΈΠ·Π²ΠΎΠ΄Π½ΡΡ ΡΡΠ½ΠΊΡΠΈΠΈ. Π Π³ΡΡΠΏΠΏΠ΅ ΡΠ°ΠΊΠΈΡ ΠΌΠ΅ΡΠΎΠ΄ΠΎΠ² ΠΎΡΠ½ΠΎΡΡΡΡΡ ΡΠ°Π·Π»ΠΈΡΠ½ΡΠ΅ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠ½ΡΠ΅ ΠΌΠ΅ΡΠΎΠ΄Ρ.
ΠΡΠ°Π΄ΠΈΠ΅Π½ΡΠ½ΡΠ΅ ΠΌΠ΅ΡΠΎΠ΄Ρ ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»ΡΡΡ ΡΠΎΠ±ΠΎΠΉ ΠΌΠ΅ΡΠΎΠ΄Ρ ΡΠΏΡΡΠΊΠ°, Π² ΠΊΠΎΡΠΎΡΡΡ Π² ΠΊΠ°ΡΠ΅ΡΡΠ²Π΅ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΠΉ ΡΠ±ΡΠ²Π°Π½ΠΈΡ Π²ΡΠ±ΠΈΡΠ°Π΅ΡΡΡ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΠ΅ Π°Π½ΡΠΈΠ³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠ° ΡΡΠ½ΠΊΡΠΈΠΈ:-grad
Π ΠΊΠ°ΡΠ΅ΡΡΠ²Π΅ Π°Π±ΡΠΎΠ»ΡΡΠ½ΠΎΠ³ΠΎ ΠΊΡΠΈΡΠ΅ΡΠΈΡ ΠΎΡΡΠ°Π½ΠΎΠ²Π° ΠΎΠ±ΡΡΠ½ΠΎ Π²ΡΠ±ΠΈΡΠ°ΡΡ Π»ΠΈΠ±ΠΎ Π½Π΅ΡΠ°Π²Π΅Π½ΡΡΠ²ΠΎ (ΠΊΡΠΈΡΠ΅ΡΠΈΠΉ ΡΠΈΠ»ΡΠ½ΠΎΠΉ ΡΡ ΠΎΠ΄ΠΈΠΌΠΎΡΡΠΈ) ΠΈΠ»ΠΈ Β¦gradΒ¦. Π ΠΊΠ°ΡΠ΅ΡΡΠ²Π΅ ΠΎΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ³ΠΎ ΠΊΡΠΈΡΠ΅ΡΠΈΡ ΠΎΡΡΠ°Π½ΠΎΠ²Π° ΠΎΠ±ΡΡΠ½ΠΎ Π²ΡΠ±ΠΈΡΠ°ΡΡ Π»ΠΈΠ±ΠΎ, Π»ΠΈΠ±ΠΎ .
Π‘ ΠΏΠΎΠΌΠΎΡΡΡ ΠΌΠ΅ΡΠΎΠ΄Π° ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ ΡΡΠ½ΠΊΡΠΈΠΈ Π²Π΄ΠΎΠ»Ρ Π·Π°Π΄Π°Π½Π½ΠΎΠ³ΠΎ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΡ ΠΎΡΡΡΠ΅ΡΡΠ²Π»ΡΠ΅ΡΡΡ ΠΏΠΎΠΈΡΠΊ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ° ΠΈΠ· Π·Π°Π΄Π°Π½Π½ΠΎΠΉ ΡΠΎΡΠΊΠΈ. Π ΠΌΠ΅ΡΠΎΠ΄Π΅ Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ° Π² ΠΊΠ°ΡΠ΅ΡΡΠ²Π΅ Π²Π΅ΠΊΡΠΎΡΠ° Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΡ ΠΈΡΠΏΠΎΠ»ΡΠ·ΡΡΡ Π²Π΅ΠΊΡΠΎΡ Π°Π½ΡΠΈΠ³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠ°.
ΠΠ΅ΡΠΎΠ΄ Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ° ΡΡΠ΅Π±ΡΠ΅Ρ ΡΠ΅ΡΠ΅Π½ΠΈΡ Π½Π° ΠΊΠ°ΠΆΠ΄ΠΎΠΌ ΡΠ°Π³Π΅ Π·Π°Π΄Π°ΡΠΈ ΠΎΠ΄Π½ΠΎΠΌΠ΅ΡΠ½ΠΎΠΉ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ, Ρ. Π΅. ΠΈΠ· ΡΠΎΡΠΊΠΈ Π±ΡΠ΄Π΅ΠΌ Π΄Π²ΠΈΠ³Π°ΡΡΡΡ Π² Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΠΈ Π°Π½ΡΠΈΠ³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠ° Π΄ΠΎ ΡΠ΅Ρ ΠΏΠΎΡ, ΠΏΠΎΠΊΠ° Π½Π΅ Π΄ΠΎΡΡΠΈΠ³Π½Π΅ΠΌ ΠΌΠΈΠ½ΠΈΠΌΡΠΌΠ° ΡΡΠ½ΠΊΡΠΈΠΈ Π½Π° ΡΡΠΎΠΌ Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΠΈ.
Π‘Ρ ΠΎΠ΄ΠΈΠΌΠΎΡΡΡ ΠΌΠ΅ΡΠΎΠ΄Π° Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ° Π·Π°Π²ΠΈΡΠΈΡ ΠΎΡ ΠΎΡΠ½ΠΎΡΠ΅Π½ΠΈΡ ΠΌΠ°ΠΊΡΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ³ΠΎ ΠΈ ΠΌΠΈΠ½ΠΈΠΌΠ°Π»ΡΠ½ΠΎΠ³ΠΎ ΡΠΎΠ±ΡΡΠ²Π΅Π½Π½ΡΡ ΡΠΈΡΠ΅Π» Π² ΠΎΠΊΡΠ΅ΡΡΠ½ΠΎΡΡΠΈ ΡΠΎΡΠΊΠΈ ΡΠΊΡΡΡΠ΅ΠΌΡΠΌΠ°. Π§Π΅ΠΌ Π±ΠΎΠ»ΡΡΠ΅ ΡΡΠΎ ΠΎΡΠ½ΠΎΡΠ΅Π½ΠΈΠ΅, ΡΠ΅ΠΌ Ρ ΡΠΆΠ΅ ΡΡ ΠΎΠ΄ΠΈΠΌΠΎΡΡΡ ΠΌΠ΅ΡΠΎΠ΄Π°.
ΠΡΠΈ Π½Π°Ρ ΠΎΠΆΠ΄Π΅Π½ΠΈΠΈ ΡΠ΅ΡΠ΅Π½ΠΈΡ Π·Π°Π΄Π°ΡΠΈ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠ½ΡΠΌΠΈ ΠΌΠ΅ΡΠΎΠ΄Π°ΠΌΠΈ ΠΈΡΠ΅ΡΠ°ΡΠΈΠΎΠ½Π½ΡΠΉ ΠΏΡΠΎΡΠ΅ΡΡ ΠΏΡΠΎΠ΄ΠΎΠ»ΠΆΠ°Π΅ΡΡΡ Π΄ΠΎ ΡΠ΅Ρ ΠΏΠΎΡ, ΠΏΠΎΠΊΠ° Π½Π΅ Π²ΡΠΏΠΎΠ»Π½ΠΈΡΡΡ Π²ΡΠ±ΡΠ°Π½Π½ΡΠΉ ΠΊΡΠΈΡΠ΅ΡΠΈΠΉ ΠΎΡΡΠ°Π½ΠΎΠ²Π°.
ΠΠ΅ΠΎΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΈΠ»Π»ΡΡΡΡΠ°ΡΠΈΡ ΠΌΠ΅ΡΠΎΠ΄Π° ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»Π΅Π½Π° Π½Π° ΡΠΈΡΡΠ½ΠΊΠ΅ 1.
Π ΠΈΡΡΠ½ΠΎΠΊ 1 — Π³Π΅ΠΎΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΈΠ»Π»ΡΡΡΡΠ°ΡΠΈΡ ΠΌΠ΅ΡΠΎΠ΄Π° Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ°
2.3 ΠΠ΅ΡΠΎΠ΄ ΡΠΎΠΏΡΡΠΆΠ΅Π½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ² ΠΠ΅ΡΠΎΠ΄ ΡΠΎΠΏΡΡΠΆΡΠ½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ² — ΠΈΡΠ΅ΡΠ°ΡΠΈΠΎΠ½Π½ΡΠΉ ΠΌΠ΅ΡΠΎΠ΄ Π΄Π»Ρ Π±Π΅Π·ΡΡΠ»ΠΎΠ²Π½ΠΎΠΉ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ Π² ΠΌΠ½ΠΎΠ³ΠΎΠΌΠ΅ΡΠ½ΠΎΠΌ ΠΏΡΠΎΡΡΡΠ°Π½ΡΡΠ²Π΅. ΠΡΠ½ΠΎΠ²Π½ΡΠΌ Π΄ΠΎΡΡΠΎΠΈΠ½ΡΡΠ²ΠΎΠΌ ΠΌΠ΅ΡΠΎΠ΄Π° ΡΠ²Π»ΡΠ΅ΡΡΡ ΡΠΎ, ΡΡΠΎ ΠΎΠ½ ΡΠ΅ΡΠ°Π΅Ρ ΠΊΠ²Π°Π΄ΡΠ°ΡΠΈΡΠ½ΡΡ Π·Π°Π΄Π°ΡΡ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ Π·Π° ΠΊΠΎΠ½Π΅ΡΠ½ΠΎΠ΅ ΡΠΈΡΠ»ΠΎ ΡΠ°Π³ΠΎΠ². ΠΠ΄Π½Π°ΠΊΠΎ, ΡΡΠΎΡ ΡΠ°Π±ΠΎΡΠΎΡΠΏΠΎΡΠΎΠ±Π½ΠΎΡΡΡ ΡΡΠΎΠ³ΠΎ ΠΌΠ΅ΡΠΎΠ΄Π° Π·Π°Π²ΠΈΡΠΈΡ ΠΎΡ Π½Π°ΡΠ°Π»ΡΠ½ΠΎΠ³ΠΎ ΠΏΡΠΈΠ±Π»ΠΈΠΆΠ΅Π½ΠΈΡ. ΠΠΎΡΡΠΎΠΌΡ ΠΏΠ΅ΡΠ΅Π΄ ΠΏΡΠΈΠΌΠ΅Π½Π΅Π½ΠΈΠ΅ΠΌ ΠΌΠ΅ΡΠΎΠ΄Π° ΡΠΎΠΏΡΡΠΆΡΠ½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ² Π½Π°Ρ ΠΎΠ΄ΡΡ ΡΠΎΡΠΊΡ Π½Π°ΡΠ°Π»ΡΠ½ΠΎΠ³ΠΎ ΠΏΡΠΈΠ±Π»ΠΈΠΆΠ΅Π½ΠΈΡ Ρ ΠΏΠΎΠΌΠΎΡΡΡ, Π½Π°ΠΏΡΠΈΠΌΠ΅Ρ, ΠΌΠ΅ΡΠΎΠ΄Π° Π½Π°ΠΈΡΠΊΠΎΡΠ΅ΠΉΡΠ΅Π³ΠΎ ΡΠΏΡΡΠΊΠ°.
Π Π°ΡΡΠΌΠΎΡΡΠΈΠΌ ΡΠ»Π΅Π΄ΡΡΡΡΡ ΡΠ°ΡΡΠ½ΡΡ Π·Π°Π΄Π°ΡΡ ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ:
ΠΠ΄Π΅ΡΡ — ΡΠΈΠΌΠΌΠ΅ΡΡΠΈΡΠ½Π°Ρ ΠΏΠΎΠ»ΠΎΠΆΠΈΡΠ΅Π»ΡΠ½ΠΎ ΠΎΠΏΡΠ΅Π΄Π΅Π»ΡΠ½Π½Π°Ρ ΠΌΠ°ΡΡΠΈΡΠ° ΡΠ°Π·ΠΌΠ΅ΡΠ°. Π’Π°ΠΊΠ°Ρ Π·Π°Π΄Π°ΡΠ° ΠΎΠΏΡΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ Π½Π°Π·ΡΠ²Π°Π΅ΡΡΡ ΠΊΠ²Π°Π΄ΡΠ°ΡΠΈΡΠ½ΠΎΠΉ. ΠΠΌΠ΅Π΅ΠΌ ΠΏΠΎΡΠ»Π΅Π΄ΠΎΠ²Π°ΡΠ΅Π»ΡΠ½ΠΎΡΡΡ ΡΠΎΡΠ΅ΠΊ Ρ i=Ρ i-1+Π±di-1, Π³Π΄Π΅ Π± Π²ΡΠ±ΠΈΡΠ°Π΅ΡΡΡ ΠΈΠ· ΡΡΠ»ΠΎΠ²ΠΈΡ, Π° Π½Π°ΠΏΡΠ°Π²Π»Π΅Π½ΠΈΠ΅ d — ΠΏΠΎ ΠΏΡΠ°Π²ΠΈΠ»Ρ, Π³Π΄Π΅. ΠΠ»Ρ ΠΊΠ²Π°Π΄ΡΠ°ΡΠΈΡΠ½ΠΎΠΉ ΡΡΠ½ΠΊΡΠΈΠΈ ΠΌΠ΅ΡΠΎΠ΄ ΡΠΎΠΏΡΡΠΆΡΠ½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ² ΡΡ ΠΎΠ΄ΠΈΡΡΡ Π·Π° n ΡΠ°Π³ΠΎΠ², Π³Π΄Π΅ n — ΡΠΈΡΠ»ΠΎ ΠΏΠ΅ΡΠ΅ΠΌΠ΅Π½Π½ΡΡ Π² ΡΡΠ½ΠΊΡΠΈΠΈ.
Π ΠΎΠ±ΡΠ΅ΠΌ ΡΠ»ΡΡΠ°Π΅ (Π½Π΅ΠΊΠ²Π°Π΄ΡΠ°ΡΠΈΡΠ½Π°Ρ ΡΡΠ½ΠΊΡΠΈΡ) Π΄Π»Ρ Π·Π°Π΄Π°ΡΠΈ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ ΠΈΠΌΠ΅Π΅ΠΌ, Π΅ΡΠ»ΠΈ ΠΈ, Π΅ΡΠ»ΠΈ
Π ΠΊΠ°ΡΠ΅ΡΡΠ²Π΅ Π°Π±ΡΠΎΠ»ΡΡΠ½ΠΎΠ³ΠΎ ΠΊΡΠΈΡΠ΅ΡΠΈΡ ΠΎΡΡΠ°Π½ΠΎΠ²Π° ΠΎΠ±ΡΡΠ½ΠΎ Π²ΡΠ±ΠΈΡΠ°ΡΡ Π»ΠΈΠ±ΠΎ Π½Π΅ΡΠ°Π²Π΅Π½ΡΡΠ²ΠΎ (ΠΊΡΠΈΡΠ΅ΡΠΈΠΉ ΡΠΈΠ»ΡΠ½ΠΎΠΉ ΡΡ ΠΎΠ΄ΠΈΠΌΠΎΡΡΠΈ) ΠΈΠ»ΠΈ Β¦gradΒ¦. Π ΠΊΠ°ΡΠ΅ΡΡΠ²Π΅ ΠΎΡΠ½ΠΎΡΠΈΡΠ΅Π»ΡΠ½ΠΎΠ³ΠΎ ΠΊΡΠΈΡΠ΅ΡΠΈΡ ΠΎΡΡΠ°Π½ΠΎΠ²Π° ΠΎΠ±ΡΡΠ½ΠΎ Π²ΡΠ±ΠΈΡΠ°ΡΡ Π»ΠΈΠ±ΠΎ, Π»ΠΈΠ±ΠΎ .
Π ΠΈΡΡΠ½ΠΎΠΊ 2 — ΠΠ΅ΠΎΠΌΠ΅ΡΡΠΈΡΠ΅ΡΠΊΠ°Ρ ΠΈΠ»Π»ΡΡΡΡΠ°ΡΠΈΡ ΠΌΠ΅ΡΠΎΠ΄Π° ΡΠΎΠΏΡΡΠΆΠ΅Π½Π½ΡΡ Π³ΡΠ°Π΄ΠΈΠ΅Π½ΡΠΎΠ²
2.4 ΠΠ»ΠΎΠΊ-ΡΡ Π΅ΠΌΠ° Π°Π»Π³ΠΎΡΠΈΡΠΌΠ° ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΠΈ
2.5 ΠΠΈΡΡΠΈΠ½Π³ ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΡ «ΠΠ΅Π·ΡΡΠ»ΠΎΠ²Π½Π°Ρ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΡ ΡΡΠ½ΠΊΡΠΈΠΈ»
Π°Π»Π³ΠΎΡΠΈΡΠΌ ΡΠΎΠΏΡΡΠΆΠ΅Π½Π½ΡΠΉ Π³ΡΠ°Π΄ΠΈΠ΅Π½Ρ ΠΌΠΈΠ½ΠΈΠΌΠΈΠ·Π°ΡΠΈΡ
unit Unit1;
interface
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, StdCtrls, XPMan;
type
TForm1 = class (TForm)
Button1: TButton;
Label1: TLabel;
Edit2: TEdit;
XPManifest1: TXPManifest;
Label2: TLabel;
Edit3: TEdit;
Button4: TButton;
Button5: TButton;
Label3: TLabel;
Edit4: TEdit;
edt1: TEdit;
lst1: TListBox;
lbl1: TLabel;
edt2: TEdit;
lst2: TListBox;
btn1: TButton;
btn2: TButton;
lbl2: TLabel;
lst3: TListBox;
btn3: TButton;
btn4: TButton;
edt3: TEdit;
Label4: TLabel;
procedure Button4Click (Sender: TObject);
procedure Button5Click (Sender: TObject);
procedure Button1Click (Sender: TObject);
procedure btn1Click (Sender: TObject);
procedure btn3Click (Sender: TObject);
procedure btn2Click (Sender: TObject);
procedure btn4Click (Sender: TObject);
private
{ Private declarations }
public
{ Public declarations }
end;
const nvar=5;
var x, x0, x00,xleft, d: array [0.(nvar-1)] of real;
eps, h, lam0: real;
n:integer;
macht, tipich: array [0.nvar-1] of real; //ΠΠΎΡΡ-ΡΡ ΠΌΠ°ΡΡΡΠ°Π±ΠΈΡΠΎΠ²Π°Π½ΠΈΡ ΠΈ ΡΠΈΠΏΠΈΡΠ½ΡΠ΅ Π·Π½Π°ΡΠ΅Π½ΠΈΡ Ρ ΠΈ Ρ-ΡΠΈΠΈ
Form1: TForm1;
implementation
{$R *.dfm}
function myfunc (x: array of real):real;
begin
myfunc:=sqr (x[0])+2*x[0]+2*sqr (x[1])+4*x[1]+5;
end;
procedure FindAlpha (x0,d:array of Real; h, eps: real; n: integer; var lam0: real);
var lam: array [1.3] of real;
i:integer;
Procedure newlam;
begin
lam[1]: =lam[2];
lam[2]:=lam[3];
lam[3]:=lam[3]+h;
h:=2*h;
end;
Procedure newx;
var i: integer;
begin
newlam;
for i:=0 to (n-1) do
x[i]: =x0[i]+d[i]*lam[3];
end;
procedure FindInterval;
var i: Integer;
begin
if myfunc (x)>myfunc (x0) then begin
h:=-h;
newx;
end;
while myfunc (x)
for i:=0 to (n-1) do begin
xleft[i]: =x0[i];
x0[i]:=x[i];
end;
newx;
end;
end;
procedure MakeItAccurate;
var i: integer;
lam1,lam2:real; //ΠΠΎΠ»ΠΎΡΠΎΠ΅ ΡΠ΅ΡΠ΅Π½ΠΈΠ΅
x1,x2: array [0.(nvar-1)] of real;
chosen: byte; //1 — ΠΈΠ½ΡΠ΅ΡΠ²Π°Π» [a, lam2], 2 — ΠΈΠ½ΡΠ΅ΡΠ²Π°Π» [lam1, b], 1 — ΠΈΠ½ΡΠ΅ΡΠ²Π°Π» [lam1, lam2]
procedure choose (lam1,lam2:real); //Π²ΡΠ±ΠΈΡΠ°Π΅ΠΌ Π½ΠΎΠ²ΡΠΉ ΠΎΡΡΠ΅Π·ΠΎΠΊ
var i: integer;
u, v: real; //Π§ΡΠΎΠ±Ρ ΠΏΠΎΡΡΠΎΡΠ½Π½ΠΎ ΡΡΠ½ΠΊΡΠΈΡ Π½Π΅ Π²ΡΡΠΈΡΠ»ΡΡΡ
begin
for i:=0 to (n-1) do begin
x1[i]: =x0[i]+lam1*d[i];
x2[i]:=x0[i]+lam2*d[i];
end;
u:=myfunc (x1);
v:=myfunc (x2);
if u
if u>v then chosen:=2;
if u=v then chosen:=3;
end;
procedure zollam; //Π‘ΡΠΈΡΠ°Π΅ΠΌ Π½ΠΎΠ²ΡΠ΅ Π»ΡΠΌΠ±Π΄Π°
begin
if chosen=1 then begin
lam[3]: =lam2;
lam2:=lam1;
lam1:=lam[3]-0.618*(lam[3]-lam[1]);
end;
if chosen=2 then begin
lam[1]: =lam1;
lam1:=lam2;
lam2:=lam[1]+0.618*(lam[3]-lam[1]);
end;
if chosen=3 then begin
lam[1]: =lam1;
lam[3]:=lam2;
lam1:=lam[3]-0.618*(lam[3]-lam[1]);
lam2:=lam[1]+0.618*(lam[3]-lam[1]);
end;
end;
function makestop: boolean;
var norm: real;
i:integer;
function relxnorm (x1,x2:array of real):real;
var relx: array [0.nvar-1] of real;
k:integer;
norme:real;
begin
for k:=0 to (n-1) do begin
relx[k]: =abs (x2[k]-x1[k]);
if abs (x2[k])>tipich[k] then relx[k]: =relx[k]/abs (x2[k])
else relx[k]: =relx[k]/tipich[k];
end;
norme:=0;
for k:=0 to (n-1) do
norme:=norme+sqr (relx[k]);
norme:=sqrt (norme);
relxnorm:=norme;
end;
begin
choose (lam1,lam2);
for i:=0 to (n-1) do begin
xleft[i]: =x0[i]+d[i]*lam[1];
x[i]:=x0[i]+d[i]*lam[3];
end;
if relxnorm (xleft, x)
else makestop:=false;
end;
begin
for i:=0 to (n-1) do x0[i]: =x00[i+1];
lam1:=lam[3]-0.618*(lam[3]-lam[1]);
lam2:=lam[1]+0.618*(lam[3]-lam[1]);
while not makestop do zollam;
end;
begin
lam[1]: =0;
lam[2]:=0;
lam[3]:=0;
newx;
FindInterval;
MakeItAccurate;
for i:=1 to n do
x[i]: =(xleft[i]+x[i])/2;
lam0:=(lam[3]+lam[1])/2;
end;
Procedure spusken (x0:array of real; eps, h: real; n: integer; var x1: array of real);
var xo, xn, agrad: array [0.(nvar-1)] of Real; //xold, xnew, antigrad.
alpha, gradnorm: Real;
j:integer;
procedure grads (x:array of Real; var agrad: array of Real);
var i: integer;
g:real;
gradtole: array [0.(nvar-1)] of Real;
begin
agrad[0]: =-(2*x[0]+2);
agrad[1]:=-(4*x[1]+4);
g:=myfunc (x);
if abs (g)>tipich[n] then g:=abs (g)
else g:=tipich[n];
for i:=0 to (n-1) do begin
gradtole[i]: =abs (-agrad[i]);
if abs (x[i])>tipich[i] then gradtole[i]: =gradtole[i]*abs (x[i])
else gradtole[i]: =gradtole[i]*tipich[i];
gradtole[i]:=gradtole[i]/g;
end;
gradnorm:=gradtole[0];
for i:=0 to (n-1) do
if gradnorm
end;
procedure xnew (alpha:real; x0: array of Real; var x: array of Real);
var i: integer;
begin
for i:=0 to (n-1) do
x[i]: =x0[i]+alpha*agrad[i];
end;
function stop (norm, eps: Real):Boolean;
begin
if norm
else stop:=false;
end;
begin
for j:=0 to (n-1) do xo[j]: =x0[j];
repeat
grads (xo, agrad);
FindAlpha (xo, agrad, h, eps, n, alpha);
xnew (alpha, xo, xn);
for j:=0 to (n-1) do xo[j]: =xn[j];
until stop (gradnorm, eps);
for j:=0 to (n-1) do x1[j]: =xn[j];
end;
procedure soprgrad (x0:array of Real; eps, h: real; n: Integer; var x: array of real);
var xold, xnew, nd, antigrad, predd: array [0.(nvar-1)] of real; //predd — ΠΏΡΠ΅Π΄ΡΠ΄ΡΡΠ΅Π΅ d, newd — Π½ΠΎΠ²ΠΎΠ΅ d
beta, teknorm, prednorm, alpha: Real; //ΡΠ΅ΠΊΡΡΠ°Ρ ΠΈ ΠΏΡΠ΅Π΄ΡΠ΄ΡΡΠ°Ρ Π½ΠΎΡΠΌΡ
k, counter: integer;
procedure grads (x:array of Real; var agrad: array of Real; gradnorm: real);
var g: real;
i:integer;
gradtole: array [0.(nvar-1)] of Real;
begin
agrad[0]: =-(2*x[0]+2);
agrad[1]:=-(4*x[1]+4);
g:=myfunc (x);
if abs (g)>tipich[n] then g:=abs (g)
else g:=tipich[n];
for i:=0 to (n-1) do begin
gradtole[i]: =abs (-agrad[i]);
if abs (x[i])>tipich[i] then gradtole[i]: =gradtole[i]*abs (x[i])
else gradtole[i]: =gradtole[i]*tipich[i];
gradtole[i]:=gradtole[i]/g;
end;
gradnorm:=gradtole[0];
for i:=1 to (n-1) do
if gradnorm
end;
function newbeta (ti:Integer; norm, prednorm: real):Real; //ti — ΡΠ΅ΠΊΡΡΠ°Ρ ΠΈΡΠ΅ΡΠ°ΡΠΈΡ
begin
if ((ti mod n)=0) then newbeta:=0
else newbeta:=Sqr (norm/prednorm);
end;
procedure newd (agrad, d: array of Real; beta: Real; var nd: array of real);
var i: integer;
begin
for i:=0 to (n-1) do
nd[i]: =agrad[i]+beta*d[i];
end;
function theend (norm, eps: real):boolean;
begin
if norm<(eps/1000) then theend:=true
else theend:=false;
end;
procedure findx (alpha:real; d, x0: array of Real; var x: array of Real);
var i: integer;
begin
for i:=0 to (n-1) do
x[i]: =x0[i]+alpha*d[i];
end;
begin
for k:=0 to (n-1) do xold[k]: =x0[k];
spusken (xold, eps, h, n, xnew);
grads (xnew, antigrad, prednorm);
for k:=0 to (n-1) do xold[k]: =xnew[k];
counter:=0;
repeat
beta:=newbeta (counter, teknorm, prednorm);
newd (antigrad, predd, beta, nd);
FindAlpha (xold, nd, h, eps, n, alpha);
findx (alpha, nd, xold, xnew);
for k:=0 to (n-1) do xold[k]: =xnew[k];
prednorm:=teknorm;
counter:=counter+1;
until theend (teknorm, eps);
for k:=0 to (n-1) do x[k]: =xnew[k];
end;
procedure TForm1. Button4Click (Sender: TObject);
var g: real;
begin
Lst1.AddItem (Edit3.Text, Edit3);
try
g:=StrToFloat (Edit3.Text);
except
ShowMessage ('Enter real number!');
Lst1.Items.Delete (Lst1.Items.Count-1);
end;
end;
procedure TForm1. Button5Click (Sender: TObject);
begin
Lst1.Items.Delete (Lst1.ItemIndex);
end;
procedure TForm1. Button1Click (Sender: TObject);
var i: integer;
g:real;
begin
try
eps:=StrToFloat (Edit4.Text);
h:=StrToFloat (Edt1.Text);
except
ShowMessage ('Enter real number!');
Edit4.Text:='';
edt1.Text:='';
end;
n:=Lst1.Items.Count;
for i:=0 to (n-1) do begin
x0[i]: =StrToFloat (Lst1.Items[i]);
macht[i]:=StrToFloat (Lst2.Items[i]);
tipich[i]:=StrToFloat (Lst3.Items[i]);
end;
g:=myfunc (tipich);
tipich[n]:=g;
soprgrad (x0,eps, h, n, x);
Edit2.Text:='';
for i:=0 to (n-1) do
Edit2.Text:=Edit2.Text+FloatToStr (x[i])+'; ';
end;
procedure TForm1. btn1Click (Sender: TObject);
var g: real;
begin
Lst2.AddItem (Edt2.Text, Edt2);
try
g:=StrToFloat (Edt2.Text);
except
ShowMessage ('Enter real number!');
Lst2.Items.Delete (Lst2.Items.Count-1);
end;
end;
procedure TForm1. btn3Click (Sender: TObject);
var g: real;
begin
Lst3.AddItem (Edt3.Text, Edit3);
try
g:=StrToFloat (Edt3.Text);
except
ShowMessage ('Enter real number!');
Lst3.Items.Delete (Lst1.Items.Count-1);
end;
end;
procedure TForm1. btn2Click (Sender: TObject);
begin
Lst2.Items.Delete (Lst2.ItemIndex);
end;
procedure TForm1. btn4Click (Sender: TObject);
begin
Lst3.Items.Delete (Lst3.ItemIndex);
end;
end.