多因子算法文档

表城市地表饮水水源地水质安全评价指数和标准限值单位：mg/L项目安全评价指数及标准限值12345DO≥7.5≥6≥5≥3≥2COD≤2≤4≤6≤10≤15氨氮（NH3-N）≤0.15≤0.5≤1.0≤1.5≤2.0表长江干流断面水质检测报告攀枝花龙洞宜昌南津关岳阳城凌安庆皖河南京林山DO9.478.178.248.628.8COD2.01.94.23.01.2NH30.420.330.360.320.27算法名称：多因子综合评价（综合加权法）功能描述：综合分析一定时期内水环境质量，有污染特别严重的因子。对外接口：城市地表饮用水水源地水质安全评价指数和标准限值二维数组A（根据物质种类可以扩充）；长江干流断面水质检测报告二维数组B（根据物质种类可以扩充）；算法流程：1隶属函数的确定：112piiiijiisctss（计算DO）121piiiijiicstss（计算除DO外的其它物质）pii指标的隶属度，i=1,2…,n;ic指标i的监测值；is为指标i的各级标准值，it第i个指标的分级值；ijt取1，2，…6等值，j=1,2…6如下例：水流截面污染物质各水质指标隶属度P计算公式一类二类三类四类五类六类DO13.9113.97.5ic7.527.56ic6365ic5453ic3532ic26icCOD12ic2242ic4364ic64106ic1051510ic615icNH310.15ic0.1520.50.15ic0.5310.5ic141.51ic1.5521.5ic26iC2计算综合评价指标水质的综合指标由所有评价项目的连续的评价分类平均值P与其最大值Pmax的加权平均得。综合指数计算公式为：11max22Ppp式中P-综合指数；pmax-隶属度Pi中的最大值；11niippn各个隶属度的算术平均值。3水质分级标准指数范围1~22~33~44~55~6》6水质级别123456对应情况清洁尚清洁轻污染中污染重污染严重污染程序结果：2.372类水文献结果：2.32类水程序matlab源码：A=[9.478.178.248.628.80;21.94.23.01.2;0.420.330.360.320.27];B=[7.56532;2461015;0.150.51.01.52.0];[M,N]=size(A);fori=1:NifA(1,i)=B(1,1)temp(1,i)=1+(13.9-A(1,i))/(13.9-7.5);elseifA(1,i)=B(1,2)temp(1,i)=2+(7.5-A(1,i))/(7.5-6);elseifA(1,i)=B(1,3)temp(1,i)=3+(6-A(1,i))/(6-5);elseifAA(1,i)=B(1,4)tempA(1,i)=4+(5-AA(1,i))/(5-3);elseifA(1,i)=B(1,5)tempA(1,i)=5+(3-AA(1,i))/(3-2);elsetempA(1,i)=12/A(1,i);endendendendendendforj=2:Mfori=1:NifA(j,i)=B(j,1)temp(j,i)=1+A(j,i)/B(j,1);elseifA(j,i)=B(j,2)temp(j,i)=2+(A(j,i)-B(j,1))/(B(j,2)-B(j,1));elseifA(j,i)=B(j,3)temp(j,i)=3+(A(j,i)-B(j,2))/(B(j,3)-B(j,2));elseifA(j,i)=B(j,4)temp(j,i)=4+(A(j,i)-B(j,3))/(B(j,4)-B(j,3));elseifA(j,i)=B(j,5)temp(j,i)=5+(A(j,i)-B(j,4))/(B(j,5)-B(j,4));elsetemp(j,i)=A(j,i)*6/B(j,5);endendendendendendendindex=mean(temp,2);average=0;pmax=index(1);[N,T]=size(index);fori=1:Naverage=average+index(i);endaverage=average/N;fori=1:N-1ifindex(i)index(i+1)pmax=index(i);index(i)=index(i+1);index(i+1)=pmax;endendP=(average+index(N))/2ifP1&P2disp('清洁');endifP2&P3disp('尚清洁');endifP3&P4disp('轻污染');endifP4&P5disp('中污染');endifP5&P6disp('重污染');endifP6disp('严重污染');end