getwater.m

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

function W=getWater(A);
[m,n]=size(A);
W=zeros(4,20);

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
    W(1,i)=A(i,7)*A(i,1)*1.3011; 
    W(2,i)=A(i,9)*A(i,3)*0.63169; 
    W(3,i)=A(i,8)*A(i,5)*A(i,1)*0.56856;
    W(4,i)=W(1,i)+W(2,i)+W(3,i);
end