免积分超清晰可编辑微分方程求解

上传人:痛*** 文档编号:216329062 上传时间:2023-06-05 格式:PPT 页数:27 大小:609KB
收藏 版权申诉 举报 下载
免积分超清晰可编辑微分方程求解_第1页
第1页 / 共27页
免积分超清晰可编辑微分方程求解_第2页
第2页 / 共27页
免积分超清晰可编辑微分方程求解_第3页
第3页 / 共27页
资源描述:

《免积分超清晰可编辑微分方程求解》由会员分享,可在线阅读,更多相关《免积分超清晰可编辑微分方程求解(27页珍藏版)》请在装配图网上搜索。

1、免积分超清晰可编辑微分方程求解 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望微分方程的解析解方法例题1:求解微分方程组初值问题clearsyms t xx1,x2,x3=dsolve(Dx1=-x1+x3+1,Dx2=x1+2*x2-1,.Dx3=-4*x1+3*x3+2,x1(0)=1,x2(0)=0,x3(0)=1)t1=0:.1:5;xt1=subs(x1,t,t1);xt2=subs(x2,t,t1);xt3=subs(x3,t,t1);plot(t1,

2、xt1,t1,xt2,t1,xt3)微分方程的数值积分微分方程的数值积分MATLAB中用来进行常微分方程数值积分的函数有好多种,例如ode23,ode45,等,ode是常微分方程(ordinary differential equation)的缩写。它们都用来解形如的一阶微分方程组在给定初始值y0时的解。对入门者而言,会一种ode函数就行。微分方程求解函数clear,close allf1=(t,x)-1,0,1;1,2,0;-4,0,3*x+1;-1;2;t_final=5;x0=1;0;1;t,x=ode23(f1,0,t_final,x0);plot(t,x)【例】Simulink 简介

3、1990 年前后出现最早的 Simulink,当时名为SimuLAB,1992 年改为 SimulinkSimulink 的名字有两重含义仿真(simu)与模型连接(link)odegroup 命令可以打开自定义模块集 Simulink 相关模块常用的模块:微分方程的Simulink建模与求解建立起微分方程的 Simulink 模型可以用 sim()函数对其模型直接求解得出微分方程的数值解【例7-30】Simulink模型图微分方程转换单个高阶常微分方程处理方法微分方程数值积分【例微分方程数值积分【例5-3-7】用数值积分法求解微分方程 设初始时间t0=0;终止时间tf=3;初始条件y(0)=

4、1,y(0)=0.解:先将方程化为两个一阶微分方程的方程组,其左端为两维变量的一阶导数。微分方程化为标准形式写成矩阵形式为 其中 为取代变量y的变量向量,为x的导数,在程序中用xdot表示。x的初始条件为 这就是待积分的微分方程组的标准形式。用MATLAB语句表述为:xdot=0,1;-t,0*x+0;1*(1-t2/pi2);【例【例5-3-7】数值解的程序】数值解的程序将微分方程的右端写成一个exn547f.m函数程序,内容如下:function xdot=exn547f(t,x)u=1-(t.2)/(pi2);xdot=0,1;-t,0*x+0;1*u;%向量导数方程主程序exn547如

5、下,它调用MATLAB中的现成的数值积分函数ode23进行积分。clf,t0=0;tf=3*pi;x0=1;0;%给出初始值t,x=ode23(exn547f,t0,tf,x0)%此处显示结果y=x(:,1);%y为x的第一列plot(t,y),grid%绘曲线xlabel(t),ylabel(y(t)数值解程序数值解程序exn547的运行结果的运行结果程序运行的结果见图5-37。这个数值积分函数是按精度要求自动选择步长的。它的默认精度为1.e-3,因此图中的积分结果是可靠的。若要改变精度要求,可在调用命令中增加备选变元,具体做法可键入help ode23查找。Simulink模型图【例】Simulink模型图【例】Simulink模型图

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

copyright@ 2023-2025  zhuangpeitu.com 装配图网版权所有   联系电话:18123376007

备案号:ICP2024067431-1 川公网安备51140202000466号


本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!