wuranzhishu.m

来自「Matlab中的一个小程序」· M 代码 · 共 45 行

function wuranzhishu(a,T)
%  DO:rong jie yang   CODMn:gao meng suan yan  NH3-N:an yan
%  a shu zu shu ju     a(1):PH  a(2):DO a(3):CODMn a(4):NH3-N a(5):wu ran ji bie
%clear;clc;
C=[2 4 6 8 10 15 inf];    %CODMn:gao meng suan yan
N=[0.15 0.5 1.0 1.5 2.0 inf];            %NH3-N:an yan
D=[7.5 6 5 3 2 0];         %DO:rong jie yang
%ph
if a(1)<=7.0
   sph=(7.0-a(1))/(7.0-6);
else
   if a(1)>7.0
      sph=(a(1)-7.0)/2.0;
   end
end
%DO
dof=468/(31.6+T);   %T: wen du============mei you fu zhi
if a(2)<D(a(5))
   sdo=10-9*(a(2))/D(a(5));
end
if a(2)>=D(a(5))
      sdo=(abs(dof-a(2)))/(dof-D(a(5)));
      %if (dof-a(2))>=0
      %   sdo=(dof-a(2))/(dof-D(a(5)));
      %else
      %   if (dof-a(2))<0
      %      sdo=(a(2)-dof)/(dof-D(a(5)));
      %   end
      %end
   %end
end
%gao meng suan yan
if a(3)==C(6)
   sc=0;
else
   sc=a(3)/C(a(5));
end
%an yan
if a(4)==N(6)
   sn=0;
else
   sn=a(4)/N(a(5));
end
sij=(sph+sdo+sc+sn)/4