getbili.m

来自「这个程序是为了实现我们中国的南水北调的应用」· M 代码 · 共 13 行

function L=getbili(A);
[m,n]=size(A);
L=zeros(1,3);

for i=1:m
    A(i,1)=A(i,1)*1.0e+004; A(i,3)=A(i,3)*1.0e+004; A(i,7)=A(i,7)*366/1000;
end

for i=1:m
    L(1)=A(i,7)*A(i,1)+L(1); L(2)=A(i,9)*A(i,3)+L(2); L(3)=A(i,8)*A(i,5)*A(i,1)+L(3);
end

L(1)=110*1.0e+008*0.4/L(1); L(2)=(110*1.0e+008*0.38)/L(2); L(3)=110*1.0e+008*0.22/L(3);