clc, clear allglobal s d1 beta1 r K d2 beta2 a p1 p2 c eta q m b s=10; d1=0.01; beta1=0.015; r=0.03; K=0.2; d2=0.01; beta2=0.015; a=1.1; p1=0.01; p2=0.01; c=0.8; eta=1; q=1.5; m=0.05; b=0.25; %% 闃堝??I1=(s.*beta1-d1.*a)./(beta1.*a);u1=-r.*beta1^2./(beta2.*K.*(beta2-p2)); v1=beta1.*r.*a./(beta2.*K.*(beta2-p2))-beta1.*(r-d2)./beta2-d1; w1=s;Delta1=v1^2-4.*u1.*w1; T2=(-v1-Delta1^0.5)./(2.*u1); I2=(s-d1.*T2)./(beta1.*T2); N2=(a-beta1.*T2)./(beta2-p2);R01=s.*beta1./(d1.*a);R11=1-(beta2-p2).*K.*(r-d2)./(a.*r);R21=1+beta1.*(r-d2)./(d1.*beta2);Rc1=1+beta1./d1.*(b./m./eta)^0.5;c1=roundn(m+b.*eta+2.*(b.*m.*eta)^0.5,-10);c1star=m+b.*eta+b./I1+m.*eta.*I1;u2=m.*eta+q.*eta.*K.*beta2./r; v2=m+b.*eta-c+q.*K.*beta2./r-q.*eta.*K.*(r-d2)./r; w2=b-q.*K.*(r-d2)./r;Delta2=v2^2-4.*u2.*w2;c2=roundn(m+b.*eta+q.*K.*beta2./r-q.*eta.*K.*(r-d2)./r+2.*(u2.*w2)^0.5,-10);c2star=m+b.*eta+m.*eta.*(r-d2)./beta2+b.*beta2./(r-d2); b1=q.*K.*(r-d2)./r;b1star=q.*K.*(r-d2)./r+u2.*(r-d2)^2./(beta2)^2;b2star=u2.*I2^2+q.*K.*(r-d2)./r;% q1=b.*r./K./(r-d2);% q1star=q1-u2.*(r-d2)./K./beta2./beta2;B=m+b.*eta-c;c3=m+(b-q.*N2).*eta+m.*eta.*I2+(b-q.*N2)./I2;%% 鏇茬嚎鏂圭▼f1p = @(c) (-(m+b.*eta-c)+((m+b.*eta-c)^2-4.*b.*m.*eta)^0.5)./(2.*m.*eta);f1n = @(c) (-(m+b.*eta-c)-((m+b.*eta-c)^2-4.*b.*m.*eta)^0.5)./(2.*m.*eta);line([0,c1star],[I1,I1],'linestyle','-','color','b','linewidth',2);hold online([c1star,2],[I1,I1],'linestyle','--','color','r','linewidth',2);hold onfplot(f1p,[0,c1star],'ShowPoles','off','linewidth',2,'linestyle',':','Color','[0.69,0.09,0.12] ');%涓嶇ǔ瀹氱殑姝e钩琛$偣Q3hold onfplot(f1n,[0, 2],'ShowPoles','off','linewidth',2,'linestyle','-','Color','[1,0.6,0.07]');%绋冲畾鐨勬骞宠 鐐筈3hold on% line([0,c2star],[(r-d2)./beta2,(r-d2)./beta2],'linestyle','--','color','k','linewidth',0.5);% hold online([c1,c1],[-1,(b./m./eta)^0.5],'linestyle','--','color','k','linewidth',0.2);hold online([c1star,c1star],[-1,I1],'linestyle','--','color','k','linewidth',0.2);hold on xlim ([0.3 0.9]);ylim([-1,10]) box on text(c1, -1.4,'c_1') text(c1star, -1.4,'c_1^*')

想得到下面这种图 !问问怎么修改 无能为力了!谢谢
