長江水質(zhì)的量化評價與預(yù)測分析建模(參考版)
2025-07-02 05:34本頁面
【正文】 )m 計(jì)算結(jié)果：未來十年每年處理污水量為（單位：億噸）：m = Columns 1 through 7 Columns 8 through 10 8．《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB3838—2002）中4個主要項(xiàng)目標(biāo)準(zhǔn)限值（ 單位：mg/L） 序號 分 類 標(biāo)準(zhǔn)值 項(xiàng) 目Ⅰ類Ⅱ類Ⅲ類Ⅳ類Ⅴ類劣Ⅴ類1溶解氧(DO) ≥（或飽和率90%）6532 02高錳酸鹽指數(shù)(CODMn)≤2461015∞3氨氮（NH3N） ≤∞4PH值（無量綱）69 。enddisp(39。%需處理的劣5類水for i=1:10, x(i)=P(i)/(L2*x1(i)*)。%長江干流未來十年劣5類水含量（單位：%）l2=x220*ones(1,10)。%長江干流未來十年廢水總量（單位：萬噸）x2=[ ]。x1=[ ]。)g(11:20) 計(jì)算結(jié)果:對長江全流域未來10年可飲用水河長比例預(yù)測結(jié)果：ans = Columns 1 through 7 Columns 8 through 10 function wenti4%.預(yù)測未來10年污水處理量L1=。disp(39。for k=1:19, x3(k+1)=x2(k+1)x2(k)。for k=0:19,x2(k+1)=(x0(1)b/a)*exp(a*k)+b/a。a=deta_a/deta。deta=(n1)*FC^2。F=sum(z1.^2)。D=x1(n)x1(1)。endfor i=1:n1, z1(i)=(x1(i+1)+x1(i))/2。x1(1)=x0(1)。%x0(i)表示第i年可飲用水的比例x0=sqrt(g0)。)L_nh3計(jì)算結(jié)果：各江段分別排放CODMn總量（單位：萬噸）L_codmn = 各江段分別排放NH3N總量（單位：萬噸）L_nh3 = function wenti3%利用灰度預(yù)測對長江全流域未來10年可飲用水河長比例預(yù)測n=9。)L_codmndisp(39。endL_codmn=L_codmn/1000/1000/10000%單位換成萬噸L_nh3=L_nh3/1000/1000/10000%單位換成萬噸disp(39。for i=2:7, L_codmn(i)=k*sum(codmn(:,i).*Q(:,i)codmn(:,i1).*Q(:,i1).*exp(a.*t(:,i1)))。%Q(i,j)表示第i月j站點(diǎn)測定流量(單位：立方米/秒)L_codmn(1)=sum(codmn(:,1).*Q(:,1))*k。%L_codmn(i)表示第i個站點(diǎn)排放codmn總量L_nh3=zeros(1,7)。%時間換為天k=3600*24*30。s=s*1000%兩相臨觀測點(diǎn)間距離,單位換算為米；for i=1:13, t(i,:)=s./(jun_v(i,:))。 endendjun_v=jun_v39。)C計(jì)算結(jié)果:綜合模糊評判矩陣:C = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 function wenti2%分別求出長江干流各江段污染物總量load datav=[ ]。 endenddisp(39。 end B(i,j)=max(tj)。A=w(i,:)。mn(i,:)。%w(i,j)表示第i個地區(qū)第j項(xiàng)指標(biāo)的權(quán)重for i=1:3, for j=1:17, w(j,i)=a(j,i)/sum(a(j,:))。a(:,3)=c(:,3)./s(3)。a(:,1)=s(1)./c(:,1)。 endendc=c./28。c=zeros(17,3)。s(2)=mean([2 4 6 10 15])。endends=zeros(1,3)。mn(i,j)=mn(i,j)./sum(mn(i,:))。 elseif x32, nh3(j,5)=nh3(j,5)+1。x32, nh3(j,4)=nh3(j,4)+(2x3)/。 nh3(j,4)=nh3(j,4)+(x31)/。 elseif x31amp。x31, nh3(j,2)=nh3(j,2)+(1x3)/。 nh3(j,2)=nh3(j,2)+()/。 elseif x3amp。 elseif x215, mn(j,5)=mn(j,5)+1。x215, mn(j,4)=mn(j,4)+(15x2)/5。 mn(j,4)=mn(j,4)+(x26)/4。 elseif x26amp。x26, mn(j,2)=mn(j,2)+(6x2)/2。 mn(j,2)=mn(j,2)+(x22)/2。 elseif x22amp。 elseif x12, do(j,5)=do(j,5)+1。x1=3, do(j,4)=do(j,4)+x12。 do(j,4)=do(j,4)(x15)/2。 elseif x13amp。x1=6, do(j,3)=do(j,3)(x16)。 do(j,2)=()/+do(j,2)。 elseif x1=6amp。 x3=d((i1)*17+j,4)。for j=1:17,for i=1:28, x1=d((i1)*17+j,2)。
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1