《污染负荷指数法程序》由会员分享,可在线阅读,更多相关《污染负荷指数法程序(2页珍藏版)》请在金锄头文库上搜索。
1、0污染负荷指数法程序污染负荷指数法程序function yy,count,count1=I_plz(A,domain) %污染负荷指数法 C=A(:,6:13);C_O=domain(:,2);y=D_count(A);for i=1:5sum1=1;if i=1a=1;b=y(i); elsea=y(i-1)+1;b=y(i);endfor j=a:bsum2=1;for k=1:8F(k)=C(j,k)/C_O(k);sum2=F(k)*sum2;FF(j,k)=F(k); end Ipl(j-a+1)=sum2(1/8);sum1=sum1*Ipl(j-a+1);endIplz(i)=s
2、um1(1/(b-a+1);endcount=zeros(5,8);for j=1:3191for k=1:8if FF(j,k)=2.5if jy(1)count(2,k)=count(2,k)+1;elseif jy(2)count(3,k)=count(3,k)+1;elseif jy(3)count(4,k)=count(4,k)+1;elseif jy(4)count(5,k)=count(5,k)+1;endendendend for i=1:5for j=1:8if i=1count1(i,j)=count(i,j)/y(i);elsecount1(i,j)=count(i,j)/(y(i)-y(i-1);endendendyy=Iplz;