常微分方程数值解法_第1页
常微分方程数值解法_第2页
常微分方程数值解法_第3页
常微分方程数值解法_第4页
全文预览已结束

下载本文档

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

文档简介

1、常微分方程数值解法【作用】微分方程建模是数学建模的重要方法,因为许多实际问题的数学描述将导致求解微分方程的定解问题。把形形色色的实际问题化成微分方程的定解问题,大体上可以按以下几步:1. 根据实际要求确定要研究的量(自变量、未知函数、必要的参数等)并确定坐标系。2. 找出这些量所满足的基本规律(物理的、几何的、化学的或生物学的等等)。3. 运用这些规律列出方程和定解条件。 基本模型1.发射卫星为什么用三级火箭2.人口模型3.战争模型4.放射性废料的处理通常需要求出方程的解来说明实际现象,并加以检验。如果能得到解析形式的解固然是便于分析和应用的,但是我们知道,只有线性常系数微分方程,并且自由项是

2、某些特殊类型的函数时,才可以得到这样的解,而绝大多数变系数方程、非线性方程都是所谓“解不出来”的于是对于用微分方程解决实际问题来说,数值解法就是一个十分重要的手段。1. 改进Euler 法:2. 龙格库塔(RungeKutta)方法:【源程序】1. 改进Euler 法:function x,y=eulerpro(fun,x0,x1,y0,n);%fun为函数,(x0,x1)为x区间,y0为初始值,n为子区间个数if nargin<5,n=50;endh=(x1-x0)/n;x(1)=x0;y(1)=y0;for i=1:nx(i+1)=x(i)+h;y1=y(i)+h*feval(fun

3、,x(i),y(i);y2=y(i)+h*feval(fun,x(i+1),y1);y(i+1)=(y1+y2)/2;end调用 command窗口f=inline('-2*y+2*x2+2*x')x,y=eulerpro(f,0,0.5,1,10)求解函数y'=2y+2+2x ,(0 x 0.5), y(0) = 12. 龙格库塔(RungeKutta)方法:t,y=solver('F',tspan,y0)这里solver为ode45,ode23,ode113,输入参数 F 是用M文件定义的微分方程y'= f (x, y)右端的函数。tspan

4、=t0,tfinal是求解区间,y0是初值。注:ode45和ode23变步长的,采用Runge-Kutta算法。ode45表示采用四阶-五阶Runge-Kutta算法,它用4阶方法提供候选解,5阶方法控制误差,是一种自适应步长(变步长)的常微分方程数值解法,其整体截断误差为(x)5。解决的是Nonstiff(非刚性)常微分方程。ode45是解决数值解问题的首选方法,若长时间没结果,应该就是刚性的,可换用ode23试试例如: odefun=(t,y) (y+3*t)/t2; %定义函数tspan=1 4; %求解区间y0=-2; %初值t,y=ode45(odefun,tspan,y0

5、);plot(t,y) %作图title('t2y''=y+3t,y(1)=-2,1<t<4')legend('t2y''=y+3t')xlabel('t')ylabel('y')% 精确解dsolve('t2*Dy=y+3*t','y(1)=-2') ans = (3*Ei(1,-1/t)+(-3*exp(-1)*Ei(1,-1)-2)/exp(-1)*exp(-1/t)求解高阶常微分方程需要求解的高阶常微分方程: 求解的关键是将高阶转为一阶,

6、odefun的书写.F(y,y',y''.y(n-1),t)=0用变量替换,y1=y,y2=y'.注意odefun方程定义为行向量程序:function Testode45tspan=3.9 4.0; %求解区间y0=2 8; %初值t,x=ode45(odefun,tspan,y0);plot(t,x(:,1),'-o',t,x(:,2),'-*')legend('y1','y2')title('y'' ''=-t*y + et*y'' +3sin2t')xlabel('t')ylabel('y')endfunction y=odefun(t,x)y=zeros(2,1); % 列向量y(1)=x(2);y(2)=-t*x(1)+exp(t)*x(2)+3*sin(2*t);end【原理】改进Euler 法:式称为由Euler 公式和梯形公式得到的预测校正系统,也叫改进Eul

温馨提示

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

评论

0/150

提交评论