早教吧作业答案频道 -->其他-->
matlab解微分方程组,这个程序报错,求指导如何改~function x_dot=cha721lotka_volterra(t,x)x_dot=zeros(2,1);x=zeros(2,1);x_dot=[4*x(1)-2*x(1)*x(2);x(1)*x(2)-3*x(2)];t_final=100;x0=[2,3];[t,x]=ode45('cha721lotka_volterra',[0,t_final],x0);p
题目详情
matlab解微分方程组,这个程序报错,求指导如何改~
function x_dot=cha721lotka_volterra(t,x)
x_dot=zeros(2,1);
x=zeros(2,1);
x_dot=[4*x(1)-2*x(1)*x(2);
x(1)*x(2)-3*x(2)];
t_final=100;
x0=[2,3];
[t,x]=ode45('cha721lotka_volterra',[0,t_final],x0);
plot(t,x)
function x_dot=cha721lotka_volterra(t,x)
x_dot=zeros(2,1);
x=zeros(2,1);
x_dot=[4*x(1)-2*x(1)*x(2);
x(1)*x(2)-3*x(2)];
t_final=100;
x0=[2,3];
[t,x]=ode45('cha721lotka_volterra',[0,t_final],x0);
plot(t,x)
▼优质解答
答案和解析
昨天不是回了一个吗?
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dxdt=zhidao_rk4_t55262177(t,x)
k1=1000;k2=1000;i1=4;i2=4;g1=10500;g2=4700;n=0.2; %把这个g2改成4700就可以了.
dxdt=[
-k1*(x(1)+n*x(2))+g1*x(1)*x(3);
-k2*(x(2)-n*x(1))+g2*x(2)*x(4);
-(1+i1+abs(x(1))^2)*x(3)+i1-1;
-(1+i2+abs(x(2))^2)*x(4)+i2-1;
];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab下面输入:
t_start=0;
t_end=10;
xx0=[0;1;0;0]; %初值
[t,x]=ode45('zhidao_rk4_t55262177',[0,t_end],xx0);
plot(x(:,1),x(:,2));
xlabel('x(1)');
ylabel('x(2)');
%具体的可以看一下论坛:
%
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dxdt=zhidao_rk4_t55262177(t,x)
k1=1000;k2=1000;i1=4;i2=4;g1=10500;g2=4700;n=0.2; %把这个g2改成4700就可以了.
dxdt=[
-k1*(x(1)+n*x(2))+g1*x(1)*x(3);
-k2*(x(2)-n*x(1))+g2*x(2)*x(4);
-(1+i1+abs(x(1))^2)*x(3)+i1-1;
-(1+i2+abs(x(2))^2)*x(4)+i2-1;
];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab下面输入:
t_start=0;
t_end=10;
xx0=[0;1;0;0]; %初值
[t,x]=ode45('zhidao_rk4_t55262177',[0,t_end],xx0);
plot(x(:,1),x(:,2));
xlabel('x(1)');
ylabel('x(2)');
%具体的可以看一下论坛:
%
看了 matlab解微分方程组,这...的网友还看了以下:
求向量组的秩 和 一个极大线性无关部分组求向量组,a1=(1,-2,2,3)^T,a2=(-2,4 2020-05-13 …
问一道线性代数的题目设α1,α2,α3是四元非齐次线性方程组Ax=b的三个解向量,且秩r(A)=3 2020-07-26 …
线性代数的题向量组a1=(1,0,2,1)T,a2=(1,2,0,1)T,a3=(2,1,3.0) 2020-07-26 …
设有向量组A:a1=(a,2,10)T,a2=(-2,1,5)T,a3=(-1,1,4)T以及向量 2020-08-02 …
线性代数求解!设α1,α2,α3是四元非齐次线性方程组AX=b的3个解向量,且R(A)=3,α1= 2020-08-03 …
非满秩矩阵的通解问题已知η1,η2,η3是四元方程组AX=b的三个解,其中A的秩R(A)=3,η1 2020-08-03 …
向量组问题,条件给足了向量组A=(a1,a2,a3).a1=(1,0,1)^T,a2=(0,1,1) 2020-10-31 …
一道线性代数题,与向量表示有关,见问题补充,确定常数a,使向量组a1=(1,1,a)t,a2=(1, 2020-10-31 …
设向量组Αα1=(1,2,1,3)T,α2=(4,-1,-5,-6)T,2)T向量组B:β1=(-1 2020-11-30 …
lingo求救急MODEL:SETS:ID/1..4/;NO(ID):a,b,n;endsetsma 2020-12-19 …