getub.m

这个程序是为了实现我们中国的南水北调的应用

function Ub=getUb(A);
[m,n]=size(A);
Ub=zeros(40,1);

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
    %Ub(i,1)=1.5*0.63169*A(i,9)*A(i,3); 
    Ub(i+20,1)=1.5*0.56856*A(i,8)*A(i,5);
end