悬赏 2 个论坛币 未解决
抄写了对外贸易、地理优势与中国的地区差异中的代码,但是运行报错,请问大家,如何解决?
[size=13.3333px]出错 figureinHuang (第 10 行) abs (w0-w0old)./w0>0&&abs(w1-w1old)./w1>0&&abs(w2-w2old)./w2>0)
[size=13.3333px]/////
clear;clc
mu=0.55; sigma=7; L0=1; theta=1/3; A=0.45; gama=1/4;
T=1.67;T01=2.1; T02=T01;
L1=linspace(0.2, 0.8, 101);
L2=1-L1;
len=length(L1);
w0=1*ones(1,len); w1=1*ones(1,len); w2=1*ones(1,len);
w0old=0; w1old=0; w2o1d=0;
while( ...
abs (w0-w0old)./w0>0.00001&&abs(w1-w1old)./w1>0.0001&&abs(w2-w2old)./w2>0.00001)
w0old=w0; w1old=w1; w2old=w2;
LA1=(A./w1).^(1/theta);
LM1=L1-LA1;
LA2=(A./w2).^(1/theta);
LM2=L2-LA2;
Y0=L0.*w0;
Y1=L1.*w1;
Y2=L2.*w2;
G0=(L0.*w0.^(1-sigma)+LM1.*(w1.*T01).*(1-sigma)+LM2.* ...
(w2.*T02).^(1-sigma).^(1/(1-sigma)));
G1=(L0.*(w0.*T01).^(1-sigma)+LM1.*w1.^(1-sigma)+LM2.* ...
(w2.*T).^(1-sigma).^(1/(1-sigma)));
G2=(L0.*(w0.*T02).^(1-sigma)+LM1.*(w1.*T)^(1-sigma)+LM2.* ...
w2.^(1-sigma).^(1/(1-sigma)));
fail=1-(1-mu).*(1+gamma.*(G1./w1));
fail2=1-(1-mu).*(1+gamma.*(G2./w2));
w0=(Y0.*G0.^(sigma-1)+fail.*Y1.*(G1./T01).^(sigma-1)+fail2.*Y2.* ...
(G2./T02).^(sigma-1)).^(1/sigma);
w1=(Y0.*(G0./T01).^(sigma-1)+fail.*Y1.*G1.^(sigma-1)+fail2.*Y2.* ...
(G2./T).^(sigma-1)).^(1/sigma);
w2=(Y0.*(G0./T02).^(sigma-1)+fail.*Y1.*(G1./T).^(sigma-1)+fail2.*Y2.* ...
G2.^(sigma-1)).^(1/sigma);
end
omegal=((w1.*(1-fail)).^(1-mu)).*((fail.*w1./G1+gamma).^mu);
omega2=((w2.*(1-fai2)).^(1-mu)).*((fai2.*w2./G2+gamma).^mu);
ratio=omegal./omega2;
plot(L1, ratio, L1, 1)