为什么用matlab解这个微分方程组解不了?function matlab secondprogramoptions = odeset('RelTol',1e-4,'AbsTol',[1e-5 1e-5]);[Z,Y] = ode45(@factor2,[0 2],[1 0],options);plot(Z,real(Y(:,1)),'-r',Z,real(Y(:,2)),'-.g');h = legend('y1','y2',3);f

来源:学生作业帮助网 编辑:作业帮 时间:2024/10/03 05:04:43
为什么用matlab解这个微分方程组解不了?function matlab secondprogramoptions = odeset('RelTol',1e-4,'AbsTol',[1e-5 1e-5]);[Z,Y] = ode45(@factor2,[0 2],[1 0],options);plot(Z,real(Y(:,1)),'-r',Z,real(Y(:,2)),'-.g');h = legend('y1','y2',3);f
xSMo@+{AmMU ;Z>&x-E8*p@T‡ZUINֻjEj\J#͛fgh͠q׉=ǝ$;f;d-p>磝;M3(h d鲠(EXF4F>0ks#Y"]Ea[IjKni,ԉU-jy_l* thX![@lP`AElM7RW$FFHuܐkZ"5WruD$X%ڣ!..3-&6>xH mŰت[ Rp)A%z@H+\'ESpA"MlOG9!O

为什么用matlab解这个微分方程组解不了?function matlab secondprogramoptions = odeset('RelTol',1e-4,'AbsTol',[1e-5 1e-5]);[Z,Y] = ode45(@factor2,[0 2],[1 0],options);plot(Z,real(Y(:,1)),'-r',Z,real(Y(:,2)),'-.g');h = legend('y1','y2',3);f
为什么用matlab解这个微分方程组解不了?
function matlab secondprogram
options = odeset('RelTol',1e-4,'AbsTol',[1e-5 1e-5]);
[Z,Y] = ode45(@factor2,[0 2],[1 0],options);
plot(Z,real(Y(:,1)),'-r',Z,real(Y(:,2)),'-.g');
h = legend('y1','y2',3);
function dy = factor2(t,y)
dy = zeros(2,1); % a column vector
k=[1 1 1 1];
b=[0.7 1];
dy = [ -i*k(1)*y(1) -i*k(2)*y(2)*exp(-i*(b(1)-b(2))*z);
-i*k(3)*y(2) -i*k(4)*y(1)*exp(-i*(b(2)-b(1))*z);]
求大虾们帮一下,就是普通的ode问题,跑不出来。

为什么用matlab解这个微分方程组解不了?function matlab secondprogramoptions = odeset('RelTol',1e-4,'AbsTol',[1e-5 1e-5]);[Z,Y] = ode45(@factor2,[0 2],[1 0],options);plot(Z,real(Y(:,1)),'-r',Z,real(Y(:,2)),'-.g');h = legend('y1','y2',3);f
dy = [ -i*k(1)*y(1) -i*k(2)*y(2)*exp(-i*(b(1)-b(2))*z);
-i*k(3)*y(2) -i*k(4)*y(1)*exp(-i*(b(2)-b(1))*z);]
其中的z 是什么啊