matlab求解一阶微分方程组

在matlab里运用ode计算一阶微分方程组
可以根据需求选择合适的ode

举个例

程序为：function dy = odefun(t,y)
dy = zeros(3,1); % 一个列向量
dy(1)=[-21y(1)+19y(2)-20y(3)];
dy(2)=[19y(1)-21y(2)+20y(3)];
dy(3)=[40y(1)-40y(2)-40*y(3)];
end
在命令窗口输入
tspan = [0,0.01];
y0 = [1,0,-1];
[T,Y] = ode113( @odefun, tspan, y0 )

其中tspan为步长，y0表示初值，T:求解区域内离散数据，Y:求解区域内离散数据的对应数值解

运行结果为
在举一个例子

程序为：
function dy = odefun2(t,y)
dy = zeros(3,1); % a column vector
dy(1)=[y(2)];
dy(2)=[-1y(1)];
dy(3)=[-1y(3)];
end
在命令窗口输入
tspan = [0:1];;y0 = [-1,0,1];
[T,Y] = ode23s( @odefun2, tspan, y0 )
出现的结果为

matlab求解一阶微分方程组相关推荐

Matlab求解常微分方程组
求解这个常微分方程组. 初始条件为其中ε取0.01,a是有上限的参数,求解方程的目的其实是找出a的临界值. syms y(t) for i = [0:0.5:1.5,1.7 ...
matlab解二阶微分方程组ode,MATLAB解含参数方程、矩阵方程、二阶微分方程组
1.如下公式: 其中x.y.m.n为参数,a.b为未知数,利用MATLAB求解方程: syms x y m n a b; [a,b]=solve('x=m*cos(a)+n*cos(a+b)','y= ...
matlab常系数线性矩阵微分方程组,基于Matlab常系数线性微分方程组的求解
·基础数学· 基于 Matlab 常系数线性微分方程组的求解* 严水仙 (赣南师范大学数学与计算机科学学院,江西赣州 341000) 摘要: 在常微分方程课程教学中,常系数线性微分方程组可以通过 ...
Matlab：Matlab编程语言应用之数学计算(求极限/渐近线求导数常微分方程求解求微分方程组的解求临界阻尼系数的解)的简介、案例实现之详细攻略
Matlab:Matlab编程语言应用之数学计算(求极限/渐近线&求导数&常微分方程求解&求微分方程组的解&求临界阻尼系数的解)的简介.案例实现之详细攻略目录三.极 ...
8 精通matlab入门求微分方程组的通解特解数值解看这篇文章就够了
%% 学习目标:使用MATLAB求解常微分方程通解特解数值解 %% 求微分方程的通解 clear all; f1=dsolve('Dy-y=sin(x)') %默认情况下以t为自变量求解求出的 ...
二阶偏微分方程组龙格库塔法_1、经典四阶龙格库塔法解一阶微分方程组
1.经典四阶龙格库塔法解一阶微分方程组陕西科技大学数值计算课程设计任务书理学院信息与计算科学/应用数学专业信息08/数学08 班级学生: 题目:典型数值算法的C++语言程序设计课程 ...
四阶龙格库塔法的基本思想_经典四阶龙格库塔法解一阶微分方程组讲义.doc
1.经典四阶龙格库塔法解一阶微分方程组 1.1运用四阶龙格库塔法解一阶微分方程组算法分析 , 经过循环计算由推得 -- 每个龙格-库塔方法都是由一个合适的泰勒方法推导而来,使得其最终全局误差为,一种 ...
第二十三讲解一阶微分方程组
例题 {u1′=−u1+2u2u1′=u1−2u2\left\{\begin{matrix}{u_{1}}'=-u_{1}+2u_{2}\\ {u_{1}}'=u_{1}-2u_{2}\end{mat ...
高等工程数学--求解非齐次线性微分方程组
今天终于弄懂了怎么求解非齐次线性微分方程组,大家都要学会哦,非常重要的,不懂的留言哦! 嘻嘻第一张图片的标题写错了,实在抱歉,是非齐次些许抽象,来看例题吧刚开始卡在了积分那里,后来灵光一现,会了 ...

matlab求解一阶微分方程组

matlab求解一阶微分方程组相关推荐

最新文章

热门文章