电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

5常微分方程组与高阶常微分方程的数值解法

5常微分方程组与高阶常微分方程的数值解法_第1页
1/10
5常微分方程组与高阶常微分方程的数值解法_第2页
2/10
5常微分方程组与高阶常微分方程的数值解法_第3页
3/10
§5 常微分方程组与 高阶常微分方程的数值解法 为了简单起见,在此仅研究含有两个一阶常微分方程的初值问题 1112101 02212202 0( ,,),(),( ,,),()yf x yyyxyyfx yyyxy    (5 .1 ) 若记 12yy y 12( , )( , )( , )f xxfx yfyy 1 002 0yy y 10020()()()yxxyx y 则(5.1)可以改写为向量的形式 00( , )()xx yfyyy (5.2) 相应地,将前述求解常微分方程初值问题的公式写成向量的形式,就得到常微分方程组的数值解法。 例如,求解(5.1)的四阶经典 R-K 方法为 1 ,111 31 11 21 42 ,122 32 12 22 4226kkkkyyKKKKhyyKKKK 其中 1121 12122 1(,,)(,,)kkkkkkf xyyKfxyyK   111 122 11 22 2211 122 1(,,)222(,,)222kkkkkkhhhf xyKyKKKhhhfxyKyK   111 222 21 32 3211 222 2(,,)222(,,)222kkkkkkhhhf xyKyKKKhhhfxyKyK   111 322 31 4211 322 32 4(,,)(,,)kkkkkkf xh yhKyhKKfxh yhKyhKK   注:教材中1122(),()kkkkyy xyy x不正确,它们仅当0k 时成立(即初始条件)。 对于高阶微分方程,可以化成微分方程组后再求解。例如二阶常微分方程初值问题 0000( , ,),(),()yf x y yy xyy xy  (5 .3 ) 令 12,yy yy, 则(5 .3 )可以化为 121001021220020,()( ,,),()()yy y xyyyf x y yy xy xy   §6 边值问题的数值解法 本节用最简单的二阶常微分方程问题为例,说明边值问题的数值解法。 ( , ,),[ , ],( ),( ).yf x y yxa by ay b 用中心差商代替导数,上述问题可以转化为方程组,然后求解。这种求解方法称为差分法。过程如下: 记 00,(1,2,, ),kbahxxkh kn xan, 将一阶和二阶中心差商公式 211()()()()2kkky xy xy xO hh 2112()2()()()()kkkky xy xy xy xO hh 中的余项2(...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

5常微分方程组与高阶常微分方程的数值解法

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部