MCP рдХреЗ рдорд╛рдзреНрдпрдо рд╕реЗ рдХрдиреЗрдХреНрдЯ рдХрд░реЗрдВ тЖТ

рдЧрдгрдирд╛ рджрд░реНрдЬ рдХрд░реЗрдВ

рд╕реВрддреНрд░ (рдлреЙрд░реНрдореВрд▓рд╛)

рд╡рд┐рдЬреНрдЮрд╛рдкрди

рдкрд░рд┐рдгрд╛рдо

y(xn) тАФ RK4 approximation
0.7615941538
рдХреНрд▓рд╛рд╕рд┐рдХрд▓ рдЪреМрдереЗ-рдХреНрд░рдо рд░реБрдВрдЧреЗ-рдХреБрдЯреНрдЯрд╛
рдЪрд░рдг-рдЖрдХрд╛рд░ h 0.02
рдЙрдкрд╡рд┐рднрд╛рдЬрди n 50
F(x,y) рдореВрд▓реНрдпрд╛рдВрдХрди 200
i / x y (рд╕рдиреНрдирд┐рдХрдЯ)
0 0
0.02 0.0199973337
0.04 0.0399786803
0.06 0.0599281034
0.08 0.079829769
0.1 0.0996679945
0.12 0.1194272984
0.14 0.1390924477
0.16 0.1586485041
0.18 0.1780808679
0.2 0.1973753199
0.22 0.2165180612
0.24 0.2354957492
0.26 0.2542955322
0.28 0.2729050801
0.3 0.291312612
0.32 0.3095069207
0.34 0.3274773943
0.36 0.3452140335
0.38 0.3627074669
0.4 0.3799489616
0.42 0.3969304313
0.44 0.4136444414
0.46 0.4300842106
0.48 0.4462436094
0.5 0.4621171563
0.52 0.4777000112
0.54 0.4929879656
0.56 0.5079774318
0.58 0.5226654285
0.6 0.5370495658
0.62 0.5511280273
0.64 0.5648995515
0.66 0.5783634117
0.68 0.591519394
0.7 0.6043677756
0.72 0.6169093013
0.74 0.6291451598
0.76 0.6410769595
0.78 0.6527067042
0.8 0.6640367685
0.82 0.675069873
0.84 0.6858090604
0.86 0.6962576708
0.88 0.7064193184
0.9 0.7162978682
0.92 0.7258974128
0.94 0.7352222508
0.96 0.7442768652
0.98 0.7530659027
1 0.7615941538

рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ

рдпрд╣ рдЯреВрд▓ \(y' = F(x, y)\) рд░реВрдк рдХреЗ рдХрд┐рд╕реА рдкреНрд░рдердо-рдХреНрд░рдо рд╕рд╛рдзрд╛рд░рдг рдЕрд╡рдХрд▓ рд╕рдореАрдХрд░рдг рдХреЛ рджреА рдЧрдИ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рд╢рд░реНрдд \(y(x_0) = y_0\) рдХреЗ рд╕рд╛рде, \(x_0\) рд╕реЗ \(x_n\) рддрдХ рдХреЗ рдЕрдВрддрд░рд╛рд▓ рдкрд░ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд░реВрдк рд╕реЗ рд╣рд▓ рдХрд░рддрд╛ рд╣реИред рдЗрд╕рдореЗрдВ рдХреНрд▓рд╛рд╕рд┐рдХрд▓ рдЪреМрдереЗ-рдХреНрд░рдо рд░реБрдВрдЧреЗ-рдХреБрдЯреНрдЯрд╛ рд╡рд┐рдзрд┐ (RK4) рдХрд╛ рдЙрдкрдпреЛрдЧ рд╣реЛрддрд╛ рд╣реИ, рдЬреЛ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╡рд┐рд╢реНрд▓реЗрд╖рдг рдореЗрдВ рд╕рдмрд╕реЗ рдЕрдзрд┐рдХ рдкреНрд░рдЪрд▓рд┐рдд рдФрд░ рднрд░реЛрд╕реЗрдордВрдж рдПрдХрд▓-рдЪрд░рдг рдЗрдВрдЯреАрдЧреНрд░реЗрдЯрд░реЛрдВ рдореЗрдВ рд╕реЗ рдПрдХ рд╣реИред рдкрд░рд┐рдгрд╛рдо рдХреЗ рд░реВрдк рдореЗрдВ рдЖрдкрдХреЛ \((x_i, y_i)\) рдмрд┐рдВрджреБрдУрдВ рдХреА рдПрдХ рд╕рд╛рд░рдгреА рдорд┐рд▓рддреА рд╣реИ рдЬреЛ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╣рд▓ рдХрд╛ рд╕рдиреНрдирд┐рдХрдЯ рдорд╛рди рджреЗрддреА рд╣реИ, рд╕рд╛рде рд╣реА рдЕрдВрддрд┐рдо рдорд╛рди \(y(x_n)\) рднреАред рдпрд╣ рд╡рд┐рд╢реБрджреНрдз рдЧрдгрд┐рддреАрдп рдЯреВрд▓ рд╣реИ, рдЗрд╕рдХрд╛ рдХрд┐рд╕реА рджреЗрд╢ рдпрд╛ рдЗрдХрд╛рдИ рд╕реЗ рдХреЛрдИ рд╕рдВрдмрдВрдз рдирд╣реАрдВ рд╣реИред

рдЗрд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ

рджрд╛рд╣рд┐рдиреА рдУрд░ рдХрд╛ рдкрдж \(F(x,y)\) рдХреЛ x рдФрд░ y рдореЗрдВ рдПрдХ рд╡реНрдпрдВрдЬрдХ рдХреЗ рд░реВрдк рдореЗрдВ рдбрд╛рд▓реЗрдВ (рдЙрджрд╛рд╣рд░рдг рдХреЗ рд▓рд┐рдП 1-y^2, x+y, x*y, рдпрд╛ sin(x)+y)ред рд╕рдорд░реНрдерд┐рдд рдСрдкрд░реЗрдЯрд░ рд╣реИрдВ + - * / ^ рдФрд░ рдлрд▓рди рдЬреИрд╕реЗ sin, cos, tan, exp, log, ln, sqrt, abs, tanh, рд╕рд╛рде рд╣реА рд╕реНрдерд┐рд░рд╛рдВрдХ pi рдФрд░ eред рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди \(x_0\), рдкреНрд░рд╛рд░рдВрднрд┐рдХ \(y_0\), рдЕрдВрддрд┐рдо \(x_n\) рд╕реЗрдЯ рдХрд░реЗрдВ, рдФрд░ рдЪреБрдиреЗрдВ рдХрд┐ рдХрд┐рддрдиреЗ рд╕рдорд╛рди рдЙрдкрд╡рд┐рднрд╛рдЬрди \(n\) рд▓реЗрдиреЗ рд╣реИрдВред рдЬрд┐рддрдиреЗ рдЕрдзрд┐рдХ рдЙрдкрд╡рд┐рднрд╛рдЬрди рд╣реЛрдВрдЧреЗ, рд╕рдЯреАрдХрддрд╛ рдЙрддрдиреА рд╣реА рдЕрдзрд┐рдХ рд╣реЛрдЧреА, рдХреНрдпреЛрдВрдХрд┐ RK4 рдХреА рд╡реИрд╢реНрд╡рд┐рдХ рддреНрд░реБрдЯрд┐ \(h^4\) рдХреЗ рдЕрдиреБрдкрд╛рдд рдореЗрдВ рдШрдЯрддреА рд╣реИред

рд╕реВрддреНрд░ рдХреА рд╡реНрдпрд╛рдЦреНрдпрд╛

рдЕрдВрддрд░рд╛рд▓ рдХреЛ \(n\) рд╕рдорд╛рди рдЪрд░рдгреЛрдВ рдореЗрдВ рдмрд╛рдБрдЯрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЬрд┐рдирдореЗрдВ рд╕реЗ рдкреНрд░рддреНрдпреЗрдХ рдХреА рдЪреМрдбрд╝рд╛рдИ \(h = (x_n - x_0)/n\) рд╣реЛрддреА рд╣реИред рд╣рд░ рдЪрд░рдг рдкрд░ RK4 рдврд╛рд▓ (slope) рдХрд╛ рдЪрд╛рд░ рдмрд╛рд░ рдорд╛рди рдирд┐рдХрд╛рд▓рддрд╛ рд╣реИ: рдПрдХ рдмрд╛рд░ рд╢реБрд░реБрдЖрдд рдореЗрдВ (\(k_1\)), рджреЛ рдмрд╛рд░ рдордзреНрдп-рдмрд┐рдВрджреБ рдкрд░ (\(k_2, k_3\)), рдФрд░ рдПрдХ рдмрд╛рд░ рдЕрдВрдд рдореЗрдВ (\(k_4\))ред рдЕрдЧрд▓рд╛ рдорд╛рди рдПрдХ рднрд╛рд░рд┐рдд рдФрд╕рдд рд╣реЛрддрд╛ рд╣реИ, $$y_{n+1} = y_n + \frac{h}{6}\left(k_1 + 2k_2 + 2k_3 + k_4\right)$$ рдЗрд╕рд╕реЗ рдЪреМрдереЗ рдХреНрд░рдо рддрдХ рдХреА рддреНрд░реБрдЯрд┐ рдХреЗ рдкрдж рдХрдЯ рдЬрд╛рддреЗ рд╣реИрдВ, рдЬрд┐рд╕рд╕реЗ рд╕реНрдерд╛рдиреАрдп рдХрдЯреМрддреА рддреНрд░реБрдЯрд┐ \(O(h^5)\) рдФрд░ рд╡реИрд╢реНрд╡рд┐рдХ рддреНрд░реБрдЯрд┐ \(O(h^4)\) рдорд┐рд▓рддреА рд╣реИред

рдПрдХ RK4 рдЪрд░рдг рдХрд╛ рдЖрд░реЗрдЦ рдЬреЛ рд╕реНрдЯреЗрдк рд╕рд╛рдЗрдЬрд╝ h рдкрд░ рдЪрд╛рд░ рдврд▓рд╛рди рдЧрдгрдирд╛рдПрдБ k1 рд╕реЗ k4 рджрд┐рдЦрд╛рддрд╛ рд╣реИ
RK4 рд╣рд░ рдЪрд░рдг рдореЗрдВ рдЪрд╛рд░ рдврд▓рд╛рди рдЕрдиреБрдорд╛рдиреЛрдВ (k1тАУk4) рдХреЛ рдорд┐рд▓рд╛рдХрд░ \((x_n, y_n)\) рд╕реЗ \((x_{n+1}, y_{n+1})\) рддрдХ рдмрдврд╝рддрд╛ рд╣реИред

рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг

\(x_0 = 0\), \(y_0 = 0\), \(x_n = 1\) рдФрд░ \(n = 10\) (\(h = 0.1\)) рдХреЗ рд╕рд╛рде \(y' = 1 - y^2\) рдХреЛ рд╣рд▓ рдХрд░реЗрдВред рд╕рдЯреАрдХ рд╣рд▓ \(y = \tanh(x)\) рд╣реИред рдкрд╣рд▓рд╛ RK4 рдЪрд░рдг \(y_1 = 0.0996679\) рджреЗрддрд╛ рд╣реИ, рдЬреЛ \(\tanh(0.1) = 0.0996680\) рд╕реЗ рдореЗрд▓ рдЦрд╛рддрд╛ рд╣реИред рд╕рднреА рджрд╕ рдЪрд░рдгреЛрдВ рдХреЗ рдмрд╛рдж \(y(1) = 0.7615942\) рдорд┐рд▓рддрд╛ рд╣реИ, рдЬреЛ \(\tanh(1) = 0.7615942\) рд╕реЗ рд╕рд╛рдд рдЕрдВрдХреЛрдВ рддрдХ рдореЗрд▓ рдЦрд╛рддрд╛ рд╣реИред

рд░реЗрдЦрд╛ рдЪрд╛рд░реНрдЯ рдЬреЛ рд╕рдЯреАрдХ рд╣рд▓ рд╡рдХреНрд░ рдХреА рддреБрд▓рдирд╛ рдЙрд╕рдХреЗ рдкрд╛рд╕ рд╕реНрдерд┐рдд RK4 рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рдмрд┐рдВрджреБрдУрдВ рд╕реЗ рдХрд░рддрд╛ рд╣реИ
RK4 рдХреЗ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рдмрд┐рдВрджреБ (рдбреЙрдЯреНрд╕) рдкреВрд░реЗ рдЕрдВрддрд░рд╛рд▓ рдореЗрдВ рд╕рдЯреАрдХ рд╣рд▓ рд╡рдХреНрд░ рдХрд╛ рдмрд╛рд░реАрдХреА рд╕реЗ рдЕрдиреБрд╕рд░рдг рдХрд░рддреЗ рд╣реИрдВред

рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рдкреНрд░рд╢реНрди

RK4, рдСрдпрд▓рд░ рд╡рд┐рдзрд┐ рд╕реЗ рдмреЗрд╣рддрд░ рдХреНрдпреЛрдВ рд╣реИ? рдСрдпрд▓рд░ рд╡рд┐рдзрд┐ рд╣рд░ рдЪрд░рдг рдореЗрдВ рдХреЗрд╡рд▓ рдПрдХ рдврд╛рд▓ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреА рд╣реИ (рддреНрд░реБрдЯрд┐ \(O(h)\))ред RK4 рдЪрд╛рд░ рдврд╛рд▓реЛрдВ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рдХреЗ рдЙрдирдХрд╛ рдФрд╕рдд рдирд┐рдХрд╛рд▓рддреА рд╣реИ, рдФрд░ рдЙрд╕реА рдЪрд░рдг-рдЖрдХрд╛рд░ рдкрд░ \(O(h^4)\) рд╕рдЯреАрдХрддрд╛ рд╣рд╛рд╕рд┐рд▓ рдХрд░рддреА рд╣реИ тАФ рдЗрд╕рд▓рд┐рдП рдХрд┐рд╕реА рдирд┐рд░реНрдзрд╛рд░рд┐рдд рд╕рдЯреАрдХрддрд╛ рдХреЗ рд▓рд┐рдП рдХрд╣реАрдВ рдХрдо рдЪрд░рдг рдЪрд╛рд╣рд┐рдПред

рдореБрдЭреЗ рдХрд┐рддрдиреЗ рдЪрд░рдг рдЪреБрдирдиреЗ рдЪрд╛рд╣рд┐рдП? 50 рд╕реЗ рд╢реБрд░реВ рдХрд░реЗрдВред рдпрджрд┐ рд╣рд▓ рд╕рд╣рдЬ (smooth) рд╣реИ, рддреЛ рдЖрдорддреМрд░ рдкрд░ рдЗрддрдирд╛ рд╣реА рдкрд░реНрдпрд╛рдкреНрдд рд╣реЛрддрд╛ рд╣реИ; рддреЗрдЬрд╝реА рд╕реЗ рдмрджрд▓рдиреЗ рд╡рд╛рд▓реА рдпрд╛ рд▓рдЧрднрдЧ рд╕реНрдЯрд┐рдлрд╝ (near-stiff) рд╕рдорд╕реНрдпрд╛рдУрдВ рдХреЗ рд▓рд┐рдП рдЗрд╕реЗ 100, 200 рдпрд╛ 500 рддрдХ рдмрдврд╝рд╛рдПрдБред

рдЕрдЧрд░ рдореБрдЭреЗ NaN рдпрд╛ Infinity рдорд┐рд▓реЗ рддреЛ рдХреНрдпрд╛ рдХрд░реЗрдВ? рд╕рдВрднрд╡ рд╣реИ рдХрд┐ рд╣рд▓ рдЕрдкрд╕рд░рд┐рдд (diverge) рд╣реЛ рдЧрдпрд╛ рд╣реЛ, рдпрд╛ \(F(x,y)\) рдиреЗ рдХрд┐рд╕реА рдЕрдорд╛рдиреНрдп рд╕рдВрдХреНрд░рд┐рдпрд╛ рдХрд╛ рд╕рд╛рдордирд╛ рдХрд┐рдпрд╛ рд╣реЛ (рдЬреИрд╕реЗ рдЛрдгрд╛рддреНрдордХ рд╕рдВрдЦреНрдпрд╛ рдХрд╛ log рдпрд╛ рд╢реВрдиреНрдп рд╕реЗ рднрд╛рдЧ)ред рд╡реНрдпрдВрдЬрдХ рдХреА рдЬрд╛рдБрдЪ рдХрд░реЗрдВ рдФрд░ рдЫреЛрдЯрд╛ рдЕрдВрддрд░рд╛рд▓ рдпрд╛ рдЕрдзрд┐рдХ рдЪрд░рдг рдЖрдЬрд╝рдорд╛рдПрдБред

рдЕрдВрддрд┐рдо рдЕрдкрдбреЗрдЯ: