麻烦帮我看下这程序为啥运行不了，在线等，谢谢各位大神咯

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

%函数
function dxdt=differential2(t,x,flag,X,Y,Z,T,M,N,O,P)
dxdt=[x(1)*(1-x(1))*[X+Y-Z-T];x(2)*(1-x(2))*[M+N-O-P]];

%y-x
E(1)=2000;r(1)=0.03;e(1)=600;
d(1)=0.08;d(2)=0.02;t(1)=0.8;k(1)=0.06;c(1)=0.4;
m=0.05;n=0.02;
E(2)=800;r(2)=0.01;e(2)=350;t(2)=0.7;k(2)=0.04;c2=0.2;
p=25;
X=x(2)*(E(1)+r(1)*e(1)*p+d(1)*e(1)+d(1)*e(2)-0.5*c(1)*e(1)*e(1)-(1-t(1))*k(1)*r(1)*e(1)*p)
Y=(1-x(2))*(E(1)+r(1)*e(1)*p-0.5*c(1)*e(1)*e(1)-m*E(1))
Z=x(2)*(E(1)+d(1)*e(2)-n*E(1))
T=(1-x(2))*(E(1)-m*E(1)-n*E(1))
M=x(1)*(E(2)+r(2)*e(2)*p+d(2)*e(2)+d(2)*e(1)-0.5*c(2)*e(2)*e(2)-(1-t(2))*k(2)*r(2)*e(2)*p)
N=(1-x(1))*(E(2)+r(2)*e(2)*p-0.5*c(2)*e(2)*e(2)-m*E(2))
O=x(1)*(E(2)+d(2)*e(1)-n*E(2))


for i=0:0.1:1
    for j=0:0.1:1

P=(1-x(1))*(E(2)-m*E(2)-n*E(2))
        [T,Y]=ode45('differential2',[0 5],[i j]);
        figure
        grid on
        plot(Y(:,1),Y(:,2));
        hold on
    end
end

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：在线等 Differential Different function DIFFER function 在线 程序

你这X.Z.M.N.O算出来并没有用到啊