§ 9-2 预估 — 校正
法
一、改进的欧拉公式
( 1 )
)(
),(
0
??
?
?
?
?
???
yay
btaytf
dt
dy
已知初值问题
??? 2
1
))(,()()( 21
t
t
dttytftyty
等价于
1111
1211
) ) ](,())(,([
2
)()(
,,,
????
??
????
????
iiiiiii
iiii
Rtytftytfhtyty
tthtttt 得且令
) ) ](,())(,([
2
)()( 111 ??? ??? iiiiii tytftytf
h
tyty
h 充分小时当
( 4)
)],(),([
2
111
0
??
?
?
?
???
?
??? iiiiii ytfytf
h
yy
y ?
并满足
的近似值与作为别用所谓改进欧拉法是指分
),,,2,1(
)()(,11
ni
tytyyy iiii
??
??
],[ )(
12
)]()([
2
)(
))](,())(,([
2
))(,(
1
3
1
111i
1
1
?
?
???
??????
??????
???
?
?
?
?
iiii
ii
t
t
iiii
t
t
tty
h
tyty
h
dtty
tytftytf
h
dttytfR
i
i
i
i
??
其中
.,2,1,0,)(,
)
2
1(12
,1
2
],,[)(
,},|),{(
),(),( 1
m ax
1
3
1
3
??????
??
?
?
?
?
?????????
?
?
??
?
iyM
f
h
Mh
f
h
baCty
ybtaytD
ytfytf
i
tt
y
i
y
y
iii
?
?
?
其中
局部截断误差满足
则改进欧拉法的时且当
解内连续
在凸区域及若函数定理
二、预估 — 校正法
)6( ),(
,)4(
)0(
1 iiii ythfyy ???
取的多次迭代解法考虑方程
.)4(),1,0(
}{)7()6(,1
2
1
),(),(
},|),{(
),( 2
1
)(
1
2121
?
?
?
?
???
?????????
i
k
i
yk
yhL
yyLytfytf
ybtaytD
yytf
的解收敛于方程
求出的序列、式由迭代公式时则当
上满足李普希兹条件
在凸区域关于变量若函数定理
?
.
,,,,
)7(,2,1,0)],,(),([
2
1
;
1
)(
1
)1(
1
)0(
1
)(
11
)1(
1
?
???
??
?
?
????
i
k
iii
k
iiiii
k
i
y
yyy
kytfytfyy
程的解在一定条件下收敛于方
求出迭代序列
用迭代公式作为初值
??
?
);,( 21
.22~21,,2,1
.;;;
.,,2,1,01
.;;,
.
)(
),(
1
ythfykS
SSni
yt
yat
n
ab
h
niyihatn
nba
ay
btaytf
dt
dy
ii
??
?
??
?
?
????
?
?
?
?
?
?
???
做对
和输出
令
)(上解近似值个等距节点
初值区间等分个数区间端点
的近似解
校正法计算初值问题—用预估
?
?
?
?
?
三、预估 — 校正算法
目标
输入
输出
S1
S2
步骤
.
,22
).(
2
1
);,(;
21
12
停机
和输出 ytS
kky
kthfyk
ihat
??
??
??
S3
作业
教材 P198 习题 1
法
一、改进的欧拉公式
( 1 )
)(
),(
0
??
?
?
?
?
???
yay
btaytf
dt
dy
已知初值问题
??? 2
1
))(,()()( 21
t
t
dttytftyty
等价于
1111
1211
) ) ](,())(,([
2
)()(
,,,
????
??
????
????
iiiiiii
iiii
Rtytftytfhtyty
tthtttt 得且令
) ) ](,())(,([
2
)()( 111 ??? ??? iiiiii tytftytf
h
tyty
h 充分小时当
( 4)
)],(),([
2
111
0
??
?
?
?
???
?
??? iiiiii ytfytf
h
yy
y ?
并满足
的近似值与作为别用所谓改进欧拉法是指分
),,,2,1(
)()(,11
ni
tytyyy iiii
??
??
],[ )(
12
)]()([
2
)(
))](,())(,([
2
))(,(
1
3
1
111i
1
1
?
?
???
??????
??????
???
?
?
?
?
iiii
ii
t
t
iiii
t
t
tty
h
tyty
h
dtty
tytftytf
h
dttytfR
i
i
i
i
??
其中
.,2,1,0,)(,
)
2
1(12
,1
2
],,[)(
,},|),{(
),(),( 1
m ax
1
3
1
3
??????
??
?
?
?
?
?????????
?
?
??
?
iyM
f
h
Mh
f
h
baCty
ybtaytD
ytfytf
i
tt
y
i
y
y
iii
?
?
?
其中
局部截断误差满足
则改进欧拉法的时且当
解内连续
在凸区域及若函数定理
二、预估 — 校正法
)6( ),(
,)4(
)0(
1 iiii ythfyy ???
取的多次迭代解法考虑方程
.)4(),1,0(
}{)7()6(,1
2
1
),(),(
},|),{(
),( 2
1
)(
1
2121
?
?
?
?
???
?????????
i
k
i
yk
yhL
yyLytfytf
ybtaytD
yytf
的解收敛于方程
求出的序列、式由迭代公式时则当
上满足李普希兹条件
在凸区域关于变量若函数定理
?
.
,,,,
)7(,2,1,0)],,(),([
2
1
;
1
)(
1
)1(
1
)0(
1
)(
11
)1(
1
?
???
??
?
?
????
i
k
iii
k
iiiii
k
i
y
yyy
kytfytfyy
程的解在一定条件下收敛于方
求出迭代序列
用迭代公式作为初值
??
?
);,( 21
.22~21,,2,1
.;;;
.,,2,1,01
.;;,
.
)(
),(
1
ythfykS
SSni
yt
yat
n
ab
h
niyihatn
nba
ay
btaytf
dt
dy
ii
??
?
??
?
?
????
?
?
?
?
?
?
???
做对
和输出
令
)(上解近似值个等距节点
初值区间等分个数区间端点
的近似解
校正法计算初值问题—用预估
?
?
?
?
?
三、预估 — 校正算法
目标
输入
输出
S1
S2
步骤
.
,22
).(
2
1
);,(;
21
12
停机
和输出 ytS
kky
kthfyk
ihat
??
??
??
S3
作业
教材 P198 习题 1