早教吧作业答案频道 -->其他-->
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解微分方程组,这...的网友还看了以下:
提先谢谢了,越快越好1.求下列函数的值:(1)已知f(x)=|x-2|分之x+1,求f(0),f( 2020-04-27 …
1.求f(x)=x²-2x-3在下列区间上的值域①R②[-3,0]③[2,3]④[0,3]2.已知 2020-05-02 …
1.函数y=根号下x=4/x=2定义域是?2.若函数f(x+1)定义域为(-1/2,2),求f(x 2020-05-13 …
已知f(x+x/1)=x^2+(1/x^2)+3,求f(x)已知f(x/x+1)=x^2+1/x^ 2020-06-07 …
(1)已知函数f(x)满足:2f(x)-f(1/x)=x+1,求f(x)(2)已知f(x-1)=x 2020-06-12 …
什么叫配凑法?1.用配凑法做一下这道题:f(x+1)=x²+2x-3求f(x)2.(1)已知f(x 2020-06-20 …
limx->无穷,x*(e^(1/x)-1)我知道答案是用洛比达法则求导分子e^(1/x)-1求导 2020-06-27 …
设全集U={(x,y)│x∈R,y∈R},M={(x,y)│y-3/x-2=1},P={(x,y) 2020-07-13 …
1.若X∧2+X=-1,求2X(X+1)+3的值.2.已知2X∧2-3X+2=0,求X∧2+1/X 2020-07-19 …
高一函数数学,逆求法求值域的题个人已经算过,不过想求下标准答案.另外,那个格式,那个话,该怎么写.题 2020-12-31 …