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
来源:学生作业帮助网 编辑:作业帮 时间:2024/08/13 16:49:40
![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](/uploads/image/z/1157097-57-7.jpg?t=matlab%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%2C%E8%BF%99%E4%B8%AA%E7%A8%8B%E5%BA%8F%E6%8A%A5%E9%94%99%2C%E6%B1%82%E6%8C%87%E5%AF%BC%E5%A6%82%E4%BD%95%E6%94%B9%7Efunction+x_dot%3Dcha721lotka_volterra%28t%2Cx%29x_dot%3Dzeros%282%2C1%29%3Bx%3Dzeros%282%2C1%29%3Bx_dot%3D%5B4%2Ax%281%29-2%2Ax%281%29%2Ax%282%29%3Bx%281%29%2Ax%282%29-3%2Ax%282%29%5D%3Bt_final%3D100%3Bx0%3D%5B2%2C3%5D%3B%5Bt%2Cx%5D%3Dode45%28%27cha721lotka_volterra%27%2C%5B0%2Ct_final%5D%2Cx0%29%3Bp)
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)
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下输入: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)');
%具体的可以看一下论坛:
%