基于MATLAB实现RAS法投入产出表.zip

利用matlab实现RAS法快速求解计算期或规划期投入产出表的目标年直接消耗系数。本压缩包有《RAS法的matlab计算代码》和《基年IO表——数据格式》两个文件。可以正常运行。

function RAS（）
%假设条件已满足了R·A·S法的基本条件，可以具体进行了。
%“R”的含义是“行乘数”，而“S”的含义是“列乘数”。
%R·A·S法的基本思路就是，计算出“行乘数”和“列乘数”，然后用它们来不断调整报告期的直接消耗系数矩阵，直到满意的结果为止。
clcclose all
format long g
warning（‘off‘）
%已知报告期投入产出表的直接消耗系数矩阵为A0=[0.2 0.2 0.1;0 0.1 0.2;0.2 0.2 0.2];
%（1）已知：基年投入产出表 1‘D4:AA33‘
A=xlsread（‘基年IO表.xlsx‘）;

%（2）根据基年投入产出表，计算基年直接消耗系数A0=Z0/X0
Z0=A（[1:24][1:24]）; %中间投入
X0=A（26[1:24]）; %总投入
[mn]=size（Z0）;
for i=1:m
    for j=1:n
        A0（ij）=Z0（ij）/X0（j）;
    end
end
Z0（1:51）
Y1=A（29[1:24]）; %最终需求
Va1=A（28[1:24]）; %增加值
%已知计划期（预测期）各部门的总产出向量 
X1=A（27[1:24]）; %目标年总投入

%（3）已知目标年总投入行向量X1，增加值行向量Va1和最终需求合计行向量Y1，计算行和U1=X1-Y1、列和V1=X1-Va1
u1=X1-Y1;  %计算行和   中间投入合计
v1=X1-Va1; %计算列和   中间需求合计

%按基年中间投入结构预测的目标年中间流量矩阵
[mn]=size（A0）;
for i=1:m
    for j=1:n
        A0X1（ij）=A0（ij）*X1（j）;
    end
end

for k=1:30
%计算出计划期各部门物资消耗的合计数和中间产品的合计数  A0Xt=ceil（A0Xt） %floor round ceil fix
disp（[‘第‘num2str（k）‘次迭代：‘]）
%行合计数ui
for i=1:m
    ui（i）=sum（A0X1（i:））;
    ri（i）=u1（i）/ui（i）;
end
ui;
ri;
%列合计数vi
for i=1:m
    for j=1:n
        vi_s（ij）=ri（i）*A0X1（ij）;
    end
end
vi=sum（vi_s）;

for i=1:length（vi）
sj（i）=v1（i）/vi（i）; 
end
sj;

for i=1:m
    for j=1:n
        A0X1i（ij）=ri（i）*A0X1（ij）*sj（j）;
    end
end
%[A0X1（1:31）  A0X1i（1:31）]
A0X1=A0X1i;

%（7）计算误差：（V1-V1‘）^2/Xi^2+（U1-U1‘）^2/Xi^2
for i=1:m
    u1T（i）=sum（A0X1（i:））;
end
v1T=sum（A0X1）;
disp（‘行合计数ui前后比较‘）
A=[u1T‘u1‘];
disp（A）
disp（‘列合计数vi前后比较‘）
A=[v1T‘v1‘];
disp（A）

err=0;
for i=1:m
    err=err+（u1T（i）-u1（i））^2/X1（i）^2+（v1T（i）-v1（i））^2/X1（i）^2;
end
err;
err0=0.00145;
if err<=err0
    disp（[‘计算误差Δ=‘num2str（err）‘< ‘num2str（err0）]）
    disp（‘现年投入产出系数矩阵‘）
    A1=A0X1*inv（diag（X1））;
    disp（A1）
    break
end
end

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       64760  2019-05-30 10:14  基年IO表.xlsx
     文件        2139  2019-05-30 10:34  RAS.m