MATLAB-常微分方程求解
MATLAB中可以用来求解常微分方程(组)的函数有ode23、 ode23s、 ode23t、 ode23tb 、ode45、ode15s和odel13等,见下表。它们的具体调用方法类似,为了方便后面的描述,
在后面的介绍中将使用solver统一代替它们。
![](/assets/blank.gif)
函数的具体调用方法如下。
[T,Y] =solver( odefun, tspan,y0)
[T,Y] = solver( odefun,tspan, y0,options)
[T,Y]=soIver( odefun, tspan,yO,options,pl ,p2,...)
在区间tspan=[t0,tf].上,使用初始条件y0,求解常微分方程y'=f(t,y)。其中解向量Y中的每行结果对应于时间向量T中的每个时间点。利用传递给函数odefun的p1, p2, .参数进行求解。solver 函数中涉及的参数见表。
![](/assets/blank.gif)
odeset函数的用法与optimset函数类似,利用所设置的可选参数进行求解,solver 中options的参数见表。
![](/assets/blank.gif)
![](/assets/blank.gif)
MATLAB可以求解3种类型一阶常微分方程,即显式常微分方程、线性隐式常微分方程和完全隐式常微分方程。
显式常微分方程
显式常微分方程的形式为
![](/assets/blank.gif)
例:显式常微分方程求解实例。
已知微分方程为y"-u(1-y^2)y'+y=0 (y(0)=0,y'(0)=2;t∈[0,30]),分别取μ=3和μ=5求解该方程。
首先对微分方程进行变换得到形式为
![](/assets/blank.gif)
然后对方程组进行函数描述,并保存为myfun6. 9. m,其内容如下所示。
function output = myfun6.9(t,y ,mu)
output = zeros(2,1);
output(1)=y(2);
output(2)= mu*(1-y(1)^2) *y(2)-y(1);
对方程组进行求解,在命令行窗口中输入如下语句。
>>[t1,y1] =ode45( @ myfun6_ 9,[0 30],[0;2],[],3); %mu=3
[t2,y2]=ode45(@ myfun6 9,[0 30],[0;2],[],5); %mu=5
plot(t1,y1(:,1),'-',t2,y2(:,2),'--')
title( '显式常微分方程的解');
xlabel('t');
ylabel('y ');
legend( 'mu=3',' mu=5');
图形窗口中的输出结果如图所示。
![](/assets/blank.gif)
线性隐式常微分方程
线性隐式常微分方程的形式为
![](/assets/blank.gif)
例:线性隐式常微分方程求解实例。
已知微分方程(ty^2+1)y'=3y"+y+4(t∈[0,10];y(0)=2),求解该方程。
首先根据微分方程(ty^2+1)y'=3y^3 +y+4和通式4M(1,y)y'=f(t,y),得到:
![](/assets/blank.gif)
然后对f(1,y)进行函数描述,并保存为myfun6_ 10f. m,其内容如下所示。
function output= myfun6_ 10f(t,y)
output=3* y. ^3+y+4;
对M(t,y)进行函数描述,并保存为myfun6_ 10M. m,其内容如下所示。
function output = myfun6_10M(t,y)
output=t*y. ^2+1;
最后对方程进行求解,在命令行窗口中输人如下语句。.
>>options = odeset( ' RelTol ',le-6, 'OutputFen ',' odeplot ',' Mass ' ,@ myfun6_ 10M);
[t,y] = ode45( @ myfun6_ 10f,[0 10] ,2 ,options);
xlabel('t');
ylabel('y ');
title( '线性隐式常微分方程的解')
图形窗口中的输出结果如图所示。
![](/assets/blank.gif)
完全隐式常微分方程
完全隐式常微分方程的形式为:
![](/assets/blank.gif)
例:完全隐式常微分方程求解实例。
已知方程如下求解该方程。
![](/assets/blank.gif)
首先对方程进行函数描述,并保存为myfun6_ 11. m,其内容如下所示。
function output = myfun6_ 11(t,y ,dydt)
output=t*y. ^2*dydt. ^3-2* y.^3*dydt. ^2+3*t*(t^2+1) *dydt -t^2*y;
其次对方程进行求解,在命令行窗口中输人如下语句。
>>t0=1;
y0=sqrt(3/2);
yp0=0;
[y0,yp0] =decie( @ myfun6_ 11 ,t0,y0,1 ,yp0,0) ;
[t,y]=ode15i( @ myfun6 11,[1 20] ,y0,yp0);
plot(t,y);
xlabel('t');
ylabel('y');
title(完全隐式常微分方程的解');
图形窗口中的输出结果如图所示。
![](/assets/blank.gif)
MATLAB-常微分方程求解相关推荐
- 视频教程-三课时精通matlab常微分方程求解-Matlab
三课时精通matlab常微分方程求解 图像和算法等领域有多年研究和项目经验:指导发表科技核心期刊经验丰富:多次指导数学建模爱好者参赛. 宋星星 ¥12.00 立即订阅 扫码下载「CSDN程序员学院AP ...
- matlab常微分方程数值求解
本节将介绍常微分方程初值问题的数值求解,主要内容分为三个部分:常微分方程数值求解的概念.求解函数及刚性问题. 一.常微分方程数值求解的一般概念 首先,凡含有参数,未知函数和未知函数导数 (或微分) 的 ...
- 微分方程求解 matlab,4MATLAB常微分方程求解.ppt
4MATLAB常微分方程求解 MATLAB微分方程 1 求简单微分方程的解析解 2 求微分方程的数值解 3 建模实例 1 求简单微分方程的解析解 求微分方程(组)的解析解命令: dsolve('方程1 ...
- Matlab:Matlab编程语言应用之数学计算(求极限/渐近线求导数常微分方程求解求微分方程组的解求临界阻尼系数的解)的简介、案例实现之详细攻略
Matlab:Matlab编程语言应用之数学计算(求极限/渐近线&求导数&常微分方程求解&求微分方程组的解&求临界阻尼系数的解)的简介.案例实现之详细攻略 目录 三.极 ...
- Simulink基础【1】-弹簧-阻尼模型的常微分方程求解
Simulink基础[1]-弹簧-阻尼模型的常微分方程求解 0. Simulink模块是什么?能干什么? 1. 弹簧阻尼模型简介 1.1 受常力的弹簧阻尼模型 1.2 动力学方程 2. simulin ...
- 微分方程Gear方法MATLAB,常微分方程的数值解
常微分方程的数值解Tag内容描述: 1.常微分方程的数值解法,引入,微分方程的数值解法是动态系统仿真的基础. 思考:数值分析课程-计算机求解数学问题? 仿真软件的实现(具体执行步骤) ?,常微分方程的 ...
- 常微分方程数值解matlab欧拉,matlab 常微分方程数值解法 源程序代码
matlab 常微分方程数值解法 源程序代码 所属分类:其他 开发工具:matlab 文件大小:16KB 下载次数:41 上传日期:2019-02-13 11:03:29 上 传 者:XWLYF 说明 ...
- matlab中predictor怎么填,在MATLAB中求解非線性有限元
我嘗試在MATLAB中求解帶有節點熱源的四面體有限元的熱擴散問題,這個節點取決於解矢量.非線性方程系統如下:在MATLAB中求解非線性有限元 乙U」 + A U = Q(T) 與B是熱capactiy ...
- 【生产优化】基于matlab遗传算法求解帐篷工序问题【含Matlab源码 2145期】
⛄一.获取代码方式 获取代码方式1: 完整代码已上传我的资源: [生产优化]基于matlab遗传算法求解帐篷工序问题[含Matlab源码 2145期] 点击上面蓝色字体,直接付费下载,即可. 获取代码 ...
- 【GA MTSP】基于matlab遗传算法求解多旅行商问题(多且同始终点)【含Matlab源码 1339期】
一.获取代码方式 获取代码方式1: 通过订阅紫极神光博客付费专栏,凭支付凭证,私信博主,可获得此代码. 获取代码方式2: 完整代码已上传我的资源:[MTSP]基于matlab遗传算法求解多旅行商问题[ ...
最新文章
- 数据蒋堂 | 大数据技术的4个E
- XML配置文件的命名空间与Spring配置文件中的头
- conda切换环境_Python多环境管理,你需要知道的几种方法
- CF 961E Tufurama
- 第四十八期:你太菜了,竟然不知道Code Review?
- 【机器学习】Pima数据集的可视化
- 汉典速查: 一个简易的国学阅读工具
- 2014版C 语言程序设计基础,清华大学出版社-图书详情-《C/C++程序设计基础与实践教程》...
- oracle wire protocol,WMB 8.0.0.2 在linux下面怎么配置连接ORACLE数据库
- 前端学习路线(详细)
- 山型组合数c语言,(人教版)高中数学选修2-3课件:组合与组合数公式自主学习...
- Download.Accelerator.Plus.v10.0.4.3_Patch-REPT
- 使用ssh关联github
- 光纤快速连接器如何安装使用?与冷接子有何区别?
- Vivado仿真报错合集(更新中)
- 二二、方程Ax=b的行空间中的解
- Python网络爬虫入门案例
- 计算机音乐里面的歌怎么关掉,电脑qq音乐怎么退出登录
- java利用poi生成pptx格式的ppt
- Java对象的持久化,将对象入库保存以及取出转换成对象实例