- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 199 个
- 通用积分
- 20.0051
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 5490 点
- 帖子
- 103
- 精华
- 0
- 在线时间
- 303 小时
- 注册时间
- 2014-5-24
- 最后登录
- 2023-11-3
博士生
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
 - 199 个
- 通用积分
- 20.0051
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 5490 点
- 帖子
- 103
- 精华
- 0
- 在线时间
- 303 小时
- 注册时间
- 2014-5-24
- 最后登录
- 2023-11-3
 | 开心 2023-3-5 08:29:50 |
|---|
签到天数: 225 天 连续签到: 1 天 [LV.7]常住居民III
|
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)
|
|