Практикум
.pdf•E‰ F )*#R o% E‰ F )*#RA!% 35 f (
" ( E 0 " Q R = X(:, E) $/ E " 3 abs(diag(R)) ( 4 3g
•N )*#R% Z ƒyJyow
(
! ;#% QR ! Q 3
! ! )
"5 "
! ! Q Q = QQ = 1
! ;#& x = Qx Q
! (x , y ) = (Qx, Qy) = (x, y)
" x, y
( f ( E‰A F )*>.5.+.#‰A A 2% ( 3 Q R
( (" 0 QR ( / X 3 3"
‚ / !% 3 Q R 0 35 QR
( / X 3 E‰ F )*#R% ( j
0 3 /" ( / X
(
V 3 ( E‰ F )* .*+#‰A A2A$% ( 3 Q R (
(" 0 ( / X 3 3" ( / X 0 j( / ( / x !% 3 Q R 0
35 QR ( / X 3 E‰A F )*#R%
( x /" ( / X
( j 0 3 /" 0 ( x
(
!
" " " #$
\\
/ 0 ( ' <
! ( 5 3 0 4 % ( 4 3
4 3 / 3 0 0 3
dy(t) t = f (t, y(t)) d
t [a, b] ( ( y(a) = y0
! ( 0 Z ∆t = h ( 5Z
3
tk = tk−1 + h = a + kh,
yk y(tk) = y(a + kh), k = 0, 1, . . . , N, fk = f (tk, yk) y (xk ).
. yk+1 − y(xk+1) ( ( ( 3( 4 3
5 0 3 - ( % 3"
5 " ( 3 y(x0), . . . , y(xk) 0 y0, . . . , yk " yk+1 − y(xk+1) ,
1 0 3 ) ( ( ( 0 3 4 3
4 " 0 ( % 4 %
4 % 3
= 1 ? 9 / 8
! ( 0 3 / ( ( 3 3 0 ( (
y (xk+1) yk+1 − yk h
S
"
yk+1 = yk + h f (tk, yk), k = 0, 1, . . . , N.
$ % 4 0 4 3 0
( ( 4 h 3 % 0 %
! 0 ( ( ( 4 3
y = 3 + t − y , t [0, 1], y(0) = 1 .
5 " 0 ( 0 0 / ( 4 3" ( % ( ( " m
0 0 (( 3 / :#+A,% ( 0
:\ @+ = : :#+A,% : + ,
0 0 (( % ( 2 m" 0 (
4 ) 0 % ( 3" 0 4 h
4 N = (b − a)/h 0 (( (
3
@5. * + @5. * ,! 7 " ,! ""!
Z #7 %
+#"% ,#"% ,! :=* "O
+# "% +# % Z ,# "% ,# % Z :#+# %A,# %%
. >
T ( " 0 3 (
-5=+#+A,%
+ +5.#E; \5.* 1.+Z=> T +Z ;A \1 +*# %A; +.- ;F%
S S
V 3 " 0 ( 4 ( 0
( f 4 0 ( Z 5 m0 0 ((
:\ @+ = , , #+% , = . # /.#+%% + .$-# +%
% ( f 0 (
-5=+#+A,A; ;A+A, #+%A; ;%
+ +5.#E; \5.* 1.+Z=> T +Z ;A \1 +*# %A; +.- ;F%
' f ( 4 3 4 3 0 3 4 ( 4
33 0 3 0 y (t)"
0 Z % " /
( 0 ( 4 ( 4 3 4 %
(
• ( / (
|
|
|
|
h |
|
|
|
|
|
|
|
|
|
|
|
yk+1 = yk + |
|
f (tk , yk) + f (tk + h, yk + h f (tk, yk)) ; |
|||||||||||||
2 |
|||||||||||||||
• ( |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
yk+1 = yk + h f tk + |
h |
+ |
h |
, yk) ; |
|||||||||||
|
, yk |
|
f (tk |
||||||||||||
2 |
2 |
||||||||||||||
• ( ^ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
yk+1 = yk + |
h |
f (tk, yk) + 3f tk + |
2 |
h, yk + |
2 |
hf (tk , yk) ; |
|||||||||
|
|
|
|
|
|
||||||||||
4 |
3 |
3 |
|||||||||||||
• ( . 0 3 Di W
yk+1 = yk + |
h |
(k1 + 4k2 + k3) , |
|
|
|
||||
6 |
|
|
|
||||||
|
|
|
|
tk + 2 , yk + |
2 |
|
|||
k1 = f (tk , yk), |
k2 |
= f |
, |
||||||
|
|
|
|
|
h |
|
k1 |
|
|
k3 = f (tk + h, yk + 2k2 − k1) .
S
! "<#" : :
\5.* ! l + -
"5 [t0, tf inal] ,!
: :
! "<#% 3 : "5 5
: :
! "<#& 0 5 "
: E :
dy |
= y − t2 + 1, t [0, 4], y(0) = 0.5 |
dt |
: h = 1 h = 0.1 0 :
: E :
: J
! "<#+ : E :
dy |
= sin(t), |
t [1, 5], |
y(1) = 0 |
dt |
: h = 1 h = 0.1 0 :
: E : : J
! "<#. : E :
dy |
= |
2y |
+ t2 exp(t); |
t [1, 3], |
y(1) = 0 |
|
|
||||
dt |
t |
: h = 1 h = 0.1 0 :
: E : : J
! "<#/ 0 :
:
S W
= 2 ( 0 !
( '
V 3 4 3 ( % / % cTV2d
( ( % /
4 3( " B?=O8D
0 4
•=>.&( 4 3 ( . R 0 3
( " ( ( 3 0 4 3 ! ( % 3% % 4
•=>. 4 3 ( . R 0 3
( ( TV2 % 3%
f ( ( 4 4 3
•=>."" ( 4 ( ( # 4$ 0 (
0 3 0 ( " ( 0 5
4 3
•=>. +7 3 ( . 4 3 ( "
0 5Z ( / 3 0 3 0 5Z ( ' ( 3 5 " f ( ( 3 f (" ( =>."(
•=>."( ( 4 ( 0 ( 0 3 c S " 0 (
5 d" 0 5Z ( / 3
0 ( " 0 ( 3 " 4 =>.&(
0 4 3
•=>. 4 ( " 0 5Z ( / 5
( . 0 3 $ 0 5
0 4 3 (
/ %
•=>. + ( 0 / 0 3/ ( % 4
0 4 " 0 5Z % (
0 ( ( % ( (
S R
•7b-&@ 3 4 3 % " 4 3 (
/ % y = f (t, y)" 0 F (y(a), y(b), p) = 0
•->.-. 3 4 3 ( 0 % f 0 %
/ % % 0 %" 3
( 3 0 % % / •Q8; sƒ" 0
4 3 J:D@A:= vAŽ8D8;@A:= k‹<:@A?;B 6??=L?C ( Z
! 4 ( 4 ( 3
y = F (t, y).
( f " 4 =>."( A=>. + 0
/ %
M(t)y = F (t, y),
M 3 ( / ( . 4 =>."( " =>. " =>. + =>. +7 ( 4 3 3
M(t, y)y = F (t, y).
/" 4 " 5 ( =>. " 0 3
( / ( " 7b-&@" ( % ( 3
M(t, y)y − F (t, y).
V 3 4 3 % ( ( 3 0
0 / 4 =>."( " =>. " =>. +" =>. +7 # 0
=>. +7 =>. 3 4 3 % /
% " =>."( r % / % /
% " =>. + r ( % / %
/ %
S
= ) ' 0 #
! 0 % / 3% 3 4 3 ( / %
0 3 5Z 3 0
•% % ( " ( / %)* * (
0 4 g
•& " 0 35Z 3 [t0, tf inal] V 3 0 3 4 ( ( ( t0, tl, . . . , tf inal c 0 0 3 ( 4 3 3d 0
tspan = [t0, tl, . . . , tf inal]g
•y0 % g
• p1, p2, . . . 0 0 ( " 0 ( / 5 F g
•T, Y ( / 4 Y " 3 ( "
Z ( / T g
( 0 5 % 0 % 0
/ 3 4 3 ( / %
• E A[F =5b.*#_lA+ - A, !% ( =5b.* 0 3 ( (3
4 3 ( / %
y = F (t, y) tspan ( 3( y0 b @F 0 •vk / c0 0 (( m d
3 ( 4 Y 5 ( "
Z ( ( / T g
•E A[F =5b.*#_lA + - A , !A =-+ = % 4 " 0
0 ( 4 " 0 ( ( " 0 3 ( ( 3(
( =-+ = " / =>. .+ T 0 ( 0 ( 5 5 0 ( 0 4.5 =5 c0 ( 5 1e −3d 0 ( % 5
0 4 7 =5 c (0 0 ( 5 1e−6dg
S [
• E A[F =5b.*#_lA+ - A, !A=-+ = A- 5A- % 4 " 0
0 ( |
4 " 0 3 0 |
0 ( |
p1, p2, . . . 9 F |
3 " 3 |
0 |
options = []" 0 ( 5 3g |
|
|
( 3 c=-+ = d ( 0 m
" ( 0 ( Z 5 ( =>. .+ - 0 ( 0 % ( %" 0 ( ( 0
( 3 % 0 ( %
•=*1n= +*=5 0 4 ( (
4 3" 0 ( 5 „of f „ 2 „on„"
|| || <= max(RelT ol||y||, AbsT ol).
( 5 4 0 5 0 0
35Z % 4 3g
• .5 =5 0 0 3
( 5 10−3 / 4 ( 4 / 0 0
(i) <= max(RelT ol |y(i|, AbsT ol(i));
•7 =5 5 3 " 0 ( 5 0 3
(0 ( 10−6 % 4 3% 3 3
3 % 35Z % 4 3" ( 0 4 3
•.: . 3 0 / (
f ( ( 5
1" ( =>.&(" 4 ' 0 ( 3 3" + - V g
•o\+-\+l@ 0 / •:\ @+ = € (
( " 4 3 3 3 /
" o\+-\+l@ 0 ( 5 / 5 =>.-5=+ !
( " 0 ( " o\+-\+l@
S X
=>.-Z =>.-Z 3 0 3 ( % %
( % % 0 g
•j+ + ( " 0 ( 5
of f g
• I ( / ( V 3 My = f (t, y) (3 0 3
( / V 3 0 ( M 0 / "
0 5Z ( / ( g
• I + +.‡.-. >. @. ( ( / ( y" ( 3 3
•none • weak • strong€ 2 none 3 M(t)y = F (t, y)
3 cweakd" 3 cstrongd ( 5 M(t, y)
y" weak 0 3 ( ( ( 4 3"
0 5Z ( 00 ( / 0 4 % g
•0 + 5 j+.- 0 ( ( 4 " 0 ( 5
4 0 3 ( 0 ( ( g
•I $j+.- ( ( 4 " 0 ( 5 % 4 3%
3 tspang
V 3 " 0 ( " % % 0 (
3 ( 0 ( =>. .+ 0 (
= 8
0
. ( ( 0 ( ( / 3 b 4
/ 3
x¨(t) = −x(t) t [0, 10π]
% x(0) = 1" x˙ (0) = 0 V 0 4 (
% ( Z *1= @ 1 :\ @+ = ,>=+ Z *1= @#+A,%
S Y