早教吧作业答案频道 -->其他-->
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解微分方程组,这...的网友还看了以下:
如图,一面积为S的单匝矩形线圈处于一个周期性变化的匀强磁场中,磁感应强度的变化规律为:B=B0si 2020-05-17 …
确定常数a,b,c使lim(x→0)(ax-sinx)/∫(上标x,下标b)ln(1+t^2)dt 2020-07-18 …
1.已知集合A={(x,y)|x^2-y^2-y=4},B={(x,y)|x^2-xy-2y^2= 2020-07-30 …
如图所示质量为m的小球在水平面内作匀速圆周运动,细线长L,偏角θ,线的拉力为F,小球作圆周运动的角 2020-07-31 …
已知集合A={2,3,4},B={x│x=t^2,t∈A},用列举法表示集合B= 2020-08-01 …
(2006•眉山)有MgSO4、Fe2(SO4)3、K2SO4三种的混合物,测得其中硫元素的质量分数 2020-11-12 …
在线求指导:听录音,选择正确的答听录音,选择正确的答语.1.Isthisyourhouse?[]AY 2020-12-22 …
一、听录音,选出你所听到的内容。(听一遍)(5分)()1.A./b/B./p/C./t/()2.A. 2020-12-22 …
重要推论1,若有f(x)的2个对称轴x=a,x=b.则T=2|a+b|2,若有f(X)的2个对称中心 2020-12-23 …
如图16-6-5所示,一个质量m=1kg的小球被长为=0.5m的细线系住,在光滑水平面上做匀速圆周运 2020-12-25 …