核心边缘模型(CP Model)Matlab程序

clear;clc

mu=0.4;sigma=5;T=1.75;fai=T^(1-sigma);

L1=linspace(0,1,21);len=length(L1);

w1=1*ones(1,len);w2=1*ones(1,len);

w1old=0;w2old=0;

val1=(w1-w1old)./w1>0.00001;val2=abs(w2-w2old)./w2>0.00001;ii=1; while val1(ii)&&val2(ii)||ii>length(val1)

w1old=w1;w2old=w2;

E1=(1-mu).*(0.5+mu.*w1.*L1./(1-mu));

delta1=L1.*w1.^(1-sigma)+fai.*(1-L1).*w2.^(1-sigma);

delta2=fai.*L1.*w1.^(1-sigma)+(1-L1).*w2.^(1-sigma);

P1=delta1.^(mu/(1-sigma));P2=delta2.^(mu/(1-sigma));

B1=E1./delta1+fai.*(1-E1)./delta2;

B2=fai.*E1./delta1+(1-E1)./delta2;

w1=B1.^(1/sigma);w2=B2.^(1/sigma);

val1=(w1-w1old)./w1>0.00001;val2=abs(w2-w2old)./w2>0.00001;

ii=ii+1;

end

omega1=w1./P1;omega2=w2./P2;

d=omega1-omega2

plot(L1,d)