二阶非线性常微分方程的打靶法_第1页
二阶非线性常微分方程的打靶法_第2页
二阶非线性常微分方程的打靶法_第3页
二阶非线性常微分方程的打靶法_第4页
二阶非线性常微分方程的打靶法_第5页
已阅读5页，还剩9页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、二阶非线性常微分方程的二阶非线性常微分方程的打靶法打靶法计算思路计算思路主要分为以下五步：给定容许误差，迭代初始值1，对k=1,2，.做：(1)用四阶Runge-Kutta 方法求解初值问，得出u1之 后 取 其 在 k的 值 ， 从 而 得 到F(k)=u1(b,k)-(2)若|F(k)|,则u1 即为所求，跳出循环。否则：否则：(3)用四阶Runge-Kutta 方法求解初值问题由此得到F(k)=v1(b,k)(4)用牛顿迭代计算k+1，即(5)置k+1k，转(1)直到误差在范围内。综上，实现了打靶法对非线性方程的拟合。接下来是每一步的代码：程序开头各变量的设置function ys=nd

2、bf(f,g,a,b,alfa,beta,n,eps,s0)%f 为二阶导数,y=f(x,y,y)，g 为f 对y 求偏导后的%(a,b)为自变量迭代区间%alfa,beta 为给定的边值条件%eps 题目规定的精度%n 为迭代次数%选取适当的s0 的初值循环第一步循环第一步x0=alfa,s0;%选取合适的迭代初值y0=RK4(f,a,b,h,x0);%龙格库塔算出u1(k)constant=y0(n,1)-beta;这里的y0(n,1)即是u1(k，b)，constant即为F(k)循环第二步循环第二步检验误差以跳出if abs(constant)eps%检验精度是否合题意break;循环

3、第三步循环第三步得到F函数的导数u0=0,1;u1=NRK4(g,a,b,h,u0,y0);通过偏导数求得的F。循环第四步循环第四步用牛顿迭代算出新的初值s0=s0-constant/u1(n,1);之后把所有的步奏放进一个while循环，跳出的条件即是constant小于给定的误差数值实验数值实验微分方程得到图像可以看出精度还是蛮高的与与matlabmatlab自带函数的比较自带函数的比较分别用Bvp4c 函数和打靶法解微分方程：得到结论：图像上的拟合度差不多，但是具体做数值误差之后Bvp4c函数的精确度比较高。Bvp4c大概为1.0e-10，打靶法只有1.0e-4再改变迭代方法进行比较再改变迭代方法进行比较用二分法迭代k此时不用再求F的导数，所以也不用输入偏导数g同样解上面的微分方程：图像的拟合度也比较高。数值误差方面和牛顿迭代在同一数量级，所以相差不大。算法总结算法总结牛顿打靶法精度

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论