用matlab解微分方程组(很急)!解这个微分方程组dx(1)=-k1*(x(1)+n*x(2))+g1*x(1)*x(3);dx(2)=-k2*(x(2)-n*x(1))+g2*x(2)*x(4);dx(3)=-(1+i1+abs(x(1))^2)*x(3)+i1-1;dx(4)=-(1+i2+abs(x(2))^2)*x(4)+i2-1;其中k1=1000;k2=1000;i1=4;i2=4;g1=10500
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/21 09:56:22
![用matlab解微分方程组(很急)!解这个微分方程组dx(1)=-k1*(x(1)+n*x(2))+g1*x(1)*x(3);dx(2)=-k2*(x(2)-n*x(1))+g2*x(2)*x(4);dx(3)=-(1+i1+abs(x(1))^2)*x(3)+i1-1;dx(4)=-(1+i2+abs(x(2))^2)*x(4)+i2-1;其中k1=1000;k2=1000;i1=4;i2=4;g1=10500](/uploads/image/z/399357-45-7.jpg?t=%E7%94%A8matlab%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%EF%BC%88%E5%BE%88%E6%80%A5%EF%BC%89%21%E8%A7%A3%E8%BF%99%E4%B8%AA%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84dx%281%29%3D-k1%2A%28x%281%29%2Bn%2Ax%282%29%29%2Bg1%2Ax%281%29%2Ax%283%29%3Bdx%282%29%3D-k2%2A%28x%282%29-n%2Ax%281%29%29%2Bg2%2Ax%282%29%2Ax%284%29%3Bdx%283%29%3D-%281%2Bi1%2Babs%28x%281%29%29%5E2%29%2Ax%283%29%2Bi1-1%3Bdx%284%29%3D-%281%2Bi2%2Babs%28x%282%29%29%5E2%29%2Ax%284%29%2Bi2-1%3B%E5%85%B6%E4%B8%ADk1%3D1000%3Bk2%3D1000%3Bi1%3D4%3Bi2%3D4%3Bg1%3D10500)
用matlab解微分方程组(很急)!解这个微分方程组dx(1)=-k1*(x(1)+n*x(2))+g1*x(1)*x(3);dx(2)=-k2*(x(2)-n*x(1))+g2*x(2)*x(4);dx(3)=-(1+i1+abs(x(1))^2)*x(3)+i1-1;dx(4)=-(1+i2+abs(x(2))^2)*x(4)+i2-1;其中k1=1000;k2=1000;i1=4;i2=4;g1=10500
用matlab解微分方程组(很急)!
解这个微分方程组
dx(1)=-k1*(x(1)+n*x(2))+g1*x(1)*x(3);
dx(2)=-k2*(x(2)-n*x(1))+g2*x(2)*x(4);
dx(3)=-(1+i1+abs(x(1))^2)*x(3)+i1-1;
dx(4)=-(1+i2+abs(x(2))^2)*x(4)+i2-1;
其中k1=1000;k2=1000;i1=4;i2=4;g1=10500;g2=4350;n=0.2;
再画出x(1)和x(2)的关系图
用matlab解微分方程组(很急)!解这个微分方程组dx(1)=-k1*(x(1)+n*x(2))+g1*x(1)*x(3);dx(2)=-k2*(x(2)-n*x(1))+g2*x(2)*x(4);dx(3)=-(1+i1+abs(x(1))^2)*x(3)+i1-1;dx(4)=-(1+i2+abs(x(2))^2)*x(4)+i2-1;其中k1=1000;k2=1000;i1=4;i2=4;g1=10500
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dxdt=zhidao_rk4_t55262177(t,x)
k1=1000;k2=1000;i1=4;i2=4;g1=10500;g2=4350;n=0.2;
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;g...
全部展开
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dxdt=zhidao_rk4_t55262177(t,x)
k1=1000;k2=1000;i1=4;i2=4;g1=10500;g2=4350;n=0.2;
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)');
用dsolve命令
收起