getf.m

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

function f=getF(A);
[m,n]=size(A);
f=zeros(1,40);

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
    %f(1,i)=(1+A(i,4))^10/A(i,9); 
    f(1,i+20)=(1+A(i,6))^10*(1+A(i,2))^10/A(i,8);
end