getlb.m

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

function Lb=getLb(A);
[m,n]=size(A);
Lb=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
    %Lb(i,1)=0.5*0.63169*A(i,9)*A(i,3); 
    Lb(i+20,1)=0.5*0.56856*A(i,8)*A(i,5);
end