常微分方程初值问题计算问题6

常微分计算问题4_百度知道
常微分计算问题4
常微分计算问题4如图,dp^-1除以-p^2要怎么算?
我有更好的答案
如图所示:这是伯努利方程,如果设u=1/p的话更容易化简出来。
采纳率:89%
来自团队:
为您推荐:
其他类似问题
换一换
回答问题,赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。您所在位置: &
&nbsp&&nbsp&nbsp&&nbsp
第6章节常微分方程数值解法.ppt 70页
本文档一共被下载:
次 ,您可全文免费在线阅读后下载本文档。
下载提示
1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值,立即自动返金币,充值渠道很便利
第6章节常微分方程数值解法
你可能关注的文档:
··········
··········
常微分方程数值解法
欧拉法和改进的欧拉法 §3
龙格-库塔法 §4
阿当姆斯方法
在高等数学里我们已经接触过常微分方程,对于一些典型的常微分方程,有求解析解的基本方法,但多数情况下,遇到的问题比较复杂,此时,只能利用近似方法求解,一般有两种近似方法 。
设f(x,y)在带形区域R:{a≤x≤b,-∞<y<+∞}上为x,y的连续函数,且对任意的y满足李普希茨(Libusize)条件
|f(x ,y1)-f(x ,y2)|≤L|y1-y2|
其中( x ,y1)、( x ,y2)∈R,L为正常数。
上的准确解y(x)的近似值
y0,y1,y2,…,yn §2
欧拉法和改进的欧拉法
一、欧拉方法 1.
区间[a,b]上给定n+1个点x0,x1,x2,……xn
再用yi近似地代替y(xi),则初值问题(6-1)(6-2) 就化为 2.几何意义
欧拉公式有很明显的几何意义。我们知道初值问题(6.1) 中的微分方程的解是xoy平面上的一簇积分曲线 这簇积分曲线上任意点(x,y)的斜率为f(x,y), 而初值问题(6.1) (6.2)的解是过点(x0,y0)的一条特定的积分曲线。
二、欧拉方法的误差分析 定义1(p146)
对于初值问题,当假设yi是准确的时,用某种方法求yi+1时所产生的截断误差称为该方法的局部截断误差 。 我们来看在第i+1步使用欧拉方法所得yi+1的局部截断误差y(xi+1)-yi+1
y(xi+1)=y(xi+h)= y(xi) + hy′(xi) + y″(ξ)/2*h2 yi+1=yi+hf(xi,yi)=y(xi) + hf(xi,y(xi)) = y(xi) + hy′(xi)
定义2 (p147)
设yi是用某种方法计算初值问题(6-1)(6-2)在xi点的近似解,而y(xi)是它的精确解,则称
为该方法的整体截断误差,也称为该方法的精度。 补:若某方法的局部截断误差为O(hp+1),则该方法的精度为p阶的。
三、改进的欧拉法
欧拉法虽然形式简单,计算方便,但比较粗糙,精度也低。特别当y=y(x)的曲线曲率较大时,欧拉法的效果更差。为了构造较高精度的数值解法,对初值问题
再做分析 为了得到更精确的方法我们可使用梯形公式
例2: 证明解常微分方程初值问题的梯形方法精度是二阶的 证明:
四、预报-校正方法
我们看到梯形法虽然提高了精度,但其算法复杂,每算一点,都需进行反复迭代,为了控制计算量,通常只迭代一两次就转入下一步计算,以简化算法。
具体地说,我们先用欧拉公式求一个初步的近似值
用预报-校正公式求解初值问题
练习1: 用Euler法解
练习2: 用改进Euler法(梯形公式)解
练习3 :用预报—校正方法求
龙格库塔法
一、基本思想
对初值问题 若用这种方法研究欧拉公式,可以发现欧拉公式仅取一个点(xi,yi)的斜率f(xi,yi)代替k*,梯形公式却是利用xi与xi+1两个点的斜率值 它也启发我们若设法在[xi,xi+1]上多找几个点的斜率值,然后将它们加权平均作为K*的近似值,则有可能构造出精度更高的计算公式,这就是Runge-Kutta方法的基本思想。 二、Runge-Kutta公式
我们以两点xi,xi+p=xi+αh为例,说明(6-14)式中系数的求法 例 4 用经典的龙格-库塔法计算
取步长h=0.2 解:由x0=0,y0=1,h=0.2,利用公式(6-20)可计算出 表6-3 补充-----单步法的收敛性和稳定性 计算过程中舍入误差总是存在的,以Euler法为例。
假设,由于舍入误差的影响,实际得到的是…..
单步法的收敛性 补充定义1:
设y(x)是初值问题(6.2)的精确解,对单步法
补充定义2:用单步法解模型方程
得到解满足稳定性方程
,就称此方法是绝对稳定的
在复平面上,所有满足
所围成的区域称为方法的绝对稳定区域. 例5:求Euler法的绝对稳定区间
阿达姆斯方法
我们已经知道,初值问题(6―1)等价于积分方程,即
对积分式分别采用矩形和梯形
正在加载中,请稍后...【图文】第6章 常微分方程数值解法_百度文库
您的浏览器Javascript被禁用,需开启后体验完整功能,
赠送免券下载特权
10W篇文档免费专享
部分付费文档8折起
每天抽奖多种福利
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
第6章 常微分方程数值解法
&&数值分析算讲稿
阅读已结束,下载本文到电脑
登录百度文库,专享文档复制特权,积分每天免费拿!
你可能喜欢常微分计算问题6_百度知道
常微分计算问题6
常微分计算问题6如图,求这一步的详细过程
我有更好的答案
u+根号(1+u^2)=y=根号(1+u^2)=y-u;两边平方,消去u^21=-2yu+y^2;u=(y^2-1)/(2y).
采纳率:90%
来自团队:
为您推荐:
其他类似问题
换一换
回答问题,赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。

我要回帖

更多关于 微分方程初值问题 的文章

 

随机推荐