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）‘????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