全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
358 0
2022-03-21
悬赏 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)




二维码

扫码加我 拉你入群

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

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

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群