資源簡介
利用matlab實現(xiàn)RAS法快速求解計算期或規(guī)劃期投入產(chǎn)出表的目標(biāo)年直接消耗系數(shù)。本壓縮包有《RAS法的matlab計算代碼》和《基年IO表——數(shù)據(jù)格式》兩個文件??梢哉_\行。

代碼片段和文件信息
function?RAS()
%假設(shè)條件已滿足了R·A·S法的基本條件,可以具體進(jìn)行了。
%“R”的含義是“行乘數(shù)”,而“S”的含義是“列乘數(shù)”。
%R·A·S法的基本思路就是,計算出“行乘數(shù)”和“列乘數(shù)”,然后用它們來不斷調(diào)整報告期的直接消耗系數(shù)矩陣,直到滿意的結(jié)果為止。
clcclose?all
format?long?g
warning(‘off‘)
%已知報告期投入產(chǎn)出表的直接消耗系數(shù)矩陣為A0=[0.2?0.2?0.1;0?0.1?0.2;0.2?0.2?0.2];
%(1)已知:基年投入產(chǎn)出表?1‘D4:AA33‘
A=xlsread(‘基年IO表.xlsx‘);
%(2)根據(jù)基年投入產(chǎn)出表,計算基年直接消耗系數(shù)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]);?%增加值
%已知計劃期(預(yù)測期)各部門的總產(chǎn)出向量?
X1=A(27[1:24]);?%目標(biāo)年總投入
%(3)已知目標(biāo)年總投入行向量X1,增加值行向量Va1和最終需求合計行向量Y1,計算行和U1=X1-Y1、列和V1=X1-Va1
u1=X1-Y1;??%計算行和???中間投入合計
v1=X1-Va1;?%計算列和???中間需求合計
%按基年中間投入結(jié)構(gòu)預(yù)測的目標(biāo)年中間流量矩陣
[mn]=size(A0);
for?i=1:m
????for?j=1:n
????????A0X1(ij)=A0(ij)*X1(j);
????end
end
for?k=1:30
%計算出計劃期各部門物資消耗的合計數(shù)和中間產(chǎn)品的合計數(shù)??A0Xt=ceil(A0Xt)?%floor?round?ceil?fix
disp([‘第‘num2str(k)‘次迭代:‘])
%行合計數(shù)ui
for?i=1:m
????ui(i)=sum(A0X1(i:));
????ri(i)=u1(i)/ui(i);
end
ui;
ri;
%列合計數(shù)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(‘行合計數(shù)ui前后比較‘)
A=[u1T‘u1‘];
disp(A)
disp(‘列合計數(shù)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(‘現(xiàn)年投入產(chǎn)出系數(shù)矩陣‘)
????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
評論
共有 條評論