请问如何实现matlab解一阶微分方程组(用其自身命令如ode45等怎么写能运行出结果)
来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/05/24 13:41:02
请问如何实现matlab解一阶微分方程组(用其自身命令如ode45等怎么写能运行出结果)
我不太会用matlab
方程如下:
dx1/dt=a*(x2-x1)+x4;
dx2/dt=d*x1-x1*x3+c*x2;
dx3/dt=x1*x2-b*x3;
dx4/dt=x2*x3+r*x4;
其中a,b,c,d,r均已知.
我不太会用matlab
方程如下:
dx1/dt=a*(x2-x1)+x4;
dx2/dt=d*x1-x1*x3+c*x2;
dx3/dt=x1*x2-b*x3;
dx4/dt=x2*x3+r*x4;
其中a,b,c,d,r均已知.
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dxdt=zhidao_rk4_5(t,x)
global a b c d r
x1=x(1);
x2=x(2);
x3=x(3);
x4=x(4);
dxdt=[
a*(x2-x1)+x4;
d*x1-x1*x3+c*x2;
x1*x2-b*x3;
x2*x3+r*x4;
];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab下面输入:
global a b c d r
a=1;
b=2;
c=3;
d=4;
r=5;
t_end=10; %时间计算的终点,即算到这个时间为止
x0=[1;1;1;1]; %x的初始值
[t,x]=ode45('zhidao_rk4_5',[0,t_end],x0);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dxdt=zhidao_rk4_5(t,x)
global a b c d r
x1=x(1);
x2=x(2);
x3=x(3);
x4=x(4);
dxdt=[
a*(x2-x1)+x4;
d*x1-x1*x3+c*x2;
x1*x2-b*x3;
x2*x3+r*x4;
];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab下面输入:
global a b c d r
a=1;
b=2;
c=3;
d=4;
r=5;
t_end=10; %时间计算的终点,即算到这个时间为止
x0=[1;1;1;1]; %x的初始值
[t,x]=ode45('zhidao_rk4_5',[0,t_end],x0);
请问如何实现matlab解一阶微分方程组(用其自身命令如ode45等怎么写能运行出结果)
我在编matlab程序的过程中,需要用ode45解一个一阶微分方程组,
matlab解微分方程组ode45
matlab ode45解微分方程组
请问如何实现使用matlab对下面的一阶微分方程组进行拉普拉斯变换?
matlab ode45求解微分方程组
matlab ode45 求解常微分方程组
求解二阶微分方程组.ps:用matlab能解这个吗?怎么写?
用matlab解微分方程组,怎么写程序呢?
如何用matlab解出如下微分方程组并作图?
Matlab编程解一阶非线性常微分方程组,谁知道程序怎么修改?
求matlab解一阶微分方程组,并画相图.