長江水質的評價和預測摘要在長江水質日趨下滑的情況下，經過對長江水質的評價和預測，對長江水質污染問題提出合理的建議和意見。針對不同的問題，建立不同的模型，提出不同的解決方案。對于問題1，為不同的水質設置不同的權重值，經過分析對水質做出定量的綜合評價；對各地區水質的污染情況，在對每一個地區在近兩年多的28個月中水質情況統計的基礎上，計算出該地區污染的種類及各種類污染出現的頻率，通過計算分析每個地區的污染情況。對于問題2，用Matlab程序解出時間矩陣T，用污染物的含量減去降解后的剩余量得到自排量。比較每一個觀測點的自排量，可以得出長江干流近一年多主要污染物高錳酸鹽指數和氨氮的污染源主要在哪些地區。對于問題3，運用灰色理論的GM（1，1）模型，解得關聯度e =0.9993，灰色關聯度精度等級為1級（好），最終得出如不治理，長江水質將會越來越差。對于問題，取19952004年文年水的IV和V類水質分別所占百分比和長江的年總流量作為預測原始值，預測出20052014年的IV和V類水質分別所占百分比和長江的年總流量。設長江污水年排放量與長江的年總流量的比作為因變量Y，利用spss 軟件的分析功能，可以估算出未來10年的年處理污水量。問題，在問題1對長江近兩年多的水質情況定量的綜合評價，并分析各地區水質的污染狀況；問題2研究、分析長江干流近一年多主要污染物高錳酸鹽指數和氨氮的污染源主要在哪些地區；問題3對長江未來水質污染的發展趨勢做出預測分析的基礎上，對解決長江水質污染問題提出一些切實可行的建議和意見。關鍵詞：水質問題 權重 灰色模型預測法 變量剔除 線性回歸一、問題重述長江是我國第一、世界第三大河流，長江水質的污染程度日趨嚴重，已引起了相關政府部門和專家們的高度重視。2004年10月，由全國政協與中國發展研究院聯合組成“保護長江萬里行”考察團，從長江上游宜賓到下游上海，對沿線21個重點城市做了實地考察，揭示了一幅長江污染的真實畫面，其污染程度讓人觸目驚心。為此，專家們提出“若不及時拯救，長江生態10年內將瀕臨崩潰”，并發出了“拿什么拯救癌變長江”的呼喚。通常認為一個觀測站的水質污染主要來自于本地區的排污和上游的污水。一般說來，江河自身對污染物都有一定的自然凈化能力，即污染物在水環境中通過物理降解、化學降解和生物降解等使水中污染物的濃度降低。反映江河自然凈化能力的指標稱為降解系數。事實上，長江干流的自然凈化能力可以認為是近似均勻的，根據檢測可知，主要污染物高錳酸鹽指數和氨氮的降解系數通常介于0.10.5之間，比如可以考慮取0.2(單位：1/天)。現需討論下列問題：（1）對長江近兩年多的水質情況做出定量的綜合評價，并分析各地區水質的污染狀況。（2）研究、分析長江干流近一年多主要污染物高錳酸鹽指數和氨氮的污染源主要在哪些地區?（3）假如不采取更有效的治理措施，依照過去10年的主要統計數據，對長江未來水質污染的發展趨勢做出預測分析，比如研究未來10年的情況。（4）根據預測分析，如果未來10年內每年都要求長江干流的類和類水的比例控制在20%以內，且沒有劣類水,那么每年需要處理多少污水？（5）對解決長江水質污染問題有什么切實可行的建議和意見。二、問題分析 (一)問題一的分析問題一可以分解為兩個方面進行討論。對于近兩年多水質情況的定量評價方面，引進評價指數和權。將三類水質在不同流域的評價指數設為ij，各類水質的百分比設為i；把I、II、III、IV、V和劣V水質的權分別記為6、5、4、3、2、1。計算出i和j，通過分析比較，定量做出水質的綜合評價。對于各地區水質的污染狀況，在上述對每一個地區在近兩年多的28個月中水質情況統計的基礎上，再找出該地區污染的種類及各種類污染出現的頻率，建立表格，通過計算分析得出結論。(二)問題2的分析設7的觀測點分別為， Si為i到i+1觀測點的距離，模型里假設各段之間水流為勻速，以兩觀測點的平均速度作為Si的速度，則 ，。用Matlab程序解出時間矩陣T，高錳酸鹽量（或氨氮含量）Q=c*q，設降解后的剩余量為，則觀測點的自排量M=Q-P。比較每一個觀測點的自排量，可以得出長江干流近一年多主要污染物高錳酸鹽指數和氨氮的污染源主要在哪些地區。(三)問題3的分析運用灰色預測的原理，建立GM (1, 1) 模型，編寫程序，得出未來年的水質情況。(四)問題4的分析取19952004年文年水的IV和V類水質分別所占百分比和長江的年總流量作為預測原始值，預測出20052014年的IV和V類水質分別所占百分比和長江的年總流量。設長江污水年排放量與長江的年總流量的比作為因變量Y，利用spss 軟件的線性回歸分析功能，對19952004 年的每年長江干流水文年各類水質所占河長的百分比做為歷史數據進行回歸分析，得到各個因子的常量系數，即可得到Y的表達式，由此得出，且未來的年處理污水量=期望長江年總流量，經過上述方法可以估算出未來10年的年處理污水量。(五)問題5的分析在問題1對長江近兩年多的水質情況定量的綜合評價，并分析各地區水質的污染狀況；問題2研究、分析長江干流近一年多主要污染物高錳酸鹽指數和氨氮的污染源主要在哪些地區；問題3對長江未來水質污染的發展趨勢做出預測分析的基礎上，對解決長江水質污染問題提出一些切實可行的建議和意見。三、模型假設1. 每一個觀測站的水質污染物主要來自本地區的排污河上游的污水，河水在通過觀測點后立即混入下段區域的污水；2. 降解作用在整個流域中是平穩進行的，主要污染物高錳酸鹽指數和氨氮的降解系數取0.2（單位：1/天）；3.不考慮洪水干旱等特殊氣候對水質的影響；42005年2014年的評價河長與2004年的評價河長相等；5.任意兩觀測點之間的江水流速是所觀測的兩點流速的平均值。四、定義與符號說明：三類水質在不同流域的評價指數（）；：第i類水質評價指數的平均值；：第j類流域評價指數的平均值；：各類水質的百分比（）；：第i個觀測點（）；：第i個觀測點到i+1個觀測點的距離（）；：第i個觀測點的水流速度（）；：從i觀測點到i+1觀測點的時間（）；：第i段流域的長度（）；：第i個觀測點自排量（）；：降解系數=0.2/天；：濃度；：水流量；：年廢水排放總量。五、模型的建立與求解（一）問題一模型的建立和求解：對于問題1，根據附件4中03年、04年的報告表，運用權重的表達式： 計算出百分比如下表1： 年份時段全流域干流支流平均值2003年枯水期3.7223.7883.7653.758豐水期3.8354.0713.7873.898水文年4.0204.2223.9784.0752004年枯水期3.6733.5483.6993.640豐水期3.7033.6783.6833.688水文年3.6883.6973.6843.690表1 近兩年來長江流域水質報告圖表2 2003年長江流域水質折線圖表3 2004年長江流域水質折線圖對于問題2，根據附件3近兩年的28個月中的水質情況進行統計，再找出該地區污染的種類及該種類污染出現的頻率。統計結果如下表4：觀測站pHDOCODMnNH3-N四川攀枝花龍洞001/281/28重慶朱沱0000湖北宜昌南津關0000湖南岳陽城陵磯0000江西九江河西水廠0000安徽安慶皖河口0000江蘇南京林山0000四川樂山岷江大橋06/284/289/28四川宜賓涼姜溝001/283/28四川瀘州沱江二橋05/281/285/28湖北丹江口胡家嶺0000湖南長沙新港01/28010/28湖南岳陽岳陽樓003/280湖北武漢宗關01/2800江西南昌滁槎1/288/282/2821/28江西九江蛤蟆石004/280江蘇揚州三江營01/2800表4 各地污染頻率統計表從上表可以看出江西南昌滁槎出現了PH值超標情況；四川樂山岷江大橋，四川瀘州沱江二橋，江蘇揚州三江營，湖北武漢宗關，江西南昌滁槎，湖南長沙新港這6個地區出現了含氧量超標情況；四川攀枝花龍洞，四川樂山岷江大橋，四川宜賓涼姜溝，四川瀘州沱江二橋，湖南岳陽岳陽樓，江西南昌滁槎，江西九江蛤蟆石這7個地區出現了高錳酸鹽超標情況；四川攀枝花龍洞，四川樂山岷江大橋，江西南昌滁槎，四川瀘州沱江二橋，湖南長沙新港，四川宜賓涼姜溝這6給地區出現了氨氮超標情況。經過上述分析統計，結合圖表可以看出2003年到2004年的水質情況下降，而且在水文年的水質較好，枯水期的水質最差。四川攀枝花龍洞CODMn，NH3-NC未達標率分別為：1/28，1/28，即該地區在抽樣的28個月內有1個月出現高錳酸鹽超標情況，1個月出現氨氮超標情況。四川樂山岷江大橋DO，CODMn，NH3-N的未達標率分別為：6/28，4/28，9/28，即該地區在抽樣的28個月內有6個月出現含氧量超標情況，4個月出現高錳酸鹽超標情況，9個月出現氨氮超標情況。四川宜賓涼姜溝CODMn，NH3-NC未達標率分別為：1/28，3/28，即該地區即該地區在抽樣的28個月內有1個月出現高錳酸鹽超標情況，3個月出現氨氮超標情況。四川瀘州沱江二橋DO，CODMn，NH3-N的未達標率分別為：5/28，1/28，5/28，即該地區在抽樣的28個月內有5個月出現含氧量超標情況，1個月出現高錳酸鹽超標情況，5個月出現氨氮超標情況。湖南長沙新港DO，NH3-N的未達標率分別為：1/28，10/28，即該地區在抽樣的28個月內有1個月出現含氧量超標況，10個月出現氨氮超標情況。湖南岳陽岳陽樓CODMn，未達標率為：3/28，即該地區即該地區在抽樣的28個月內有3個月出現高錳酸鹽超標情況。湖北武漢宗關湖南岳陽岳陽樓DO未達標率為：1/28，即該地區即該地區在抽樣的28個月內有1個月出現含氧量超標情況。江西南昌滁槎PH，DO，CODMn，NH3-NC未達標率分別為：1/28，8/28，2/28，21/28即該地區在抽樣的28個月內有1個月出現PH超標情況，8個月出現含氧量超標情況，2個月出現高錳酸鹽超標情況，21個月出現氨氮超標情況。江西九江蛤蟆石CODMn，未達標率為：4/28，即該地區即該地區在抽樣的28個月內有4個月出現高錳酸鹽超標情況。湖北武漢宗關湖南岳陽岳陽樓DO未達標率為：1/28，即該地區即該地區在抽樣的28個月內有1個月出現含氧量超標情況。（二）問題二模型的建立和求解：模型里假設各段之間水流為勻速，以兩觀測點的平均速度作為Si的速度，則，記對應的流域段長度為，則。分別運用Matlab程序解出：；速度矩陣江水流過所需時間矩陣單位時間內高錳酸鹽量（單位：g）：，運用Matlab程序解得：單位時間內氨氮量（單位：g）：，運用Matlab程序解得：上游高錳酸鹽經過降解到下游時剩余高錳酸鹽量，上游氨氮經過降解到下游相鄰觀測點時剩余氨氮量，運用Matlab程序解得：Q1 = 19.0, 3.09, 13.0, 4.68, .443e4, 45.0 28.6, .627, 10.3, 6.62, .571e4, 5.88 27.3, 6.50, 213., 103., .977e4, 196. 44.1, 21.3, 456., 203., .985e4, 290. 53.6, 110., 385., 168., .946e4, 392. .108e4, .243e4, .305e4, .200e4, .456e5, .934e4 3.73, 12.6, 633., 230., .740e4, 538. 1.92, .234, 1.10, .542, .108e4, 4.20 .694, .487e-1, .968, .331, 669., .917 .458e-1, .268e-2, .135e-1, .251e-1, 334., .316 .415e-2, .624e-5, .369e-2, .154e-1, 407., .366 .223e-2, .173e-5, .445e-2, .478e-1, 568., 2.49 .156e-1, .922e-4, .962e-3, .384e-2, 467., .948Q2 = 1.24, .185, 1.62, .369, 412., 3.53 .466, .770e-1, 1.03, .589, 405., 4.63 .436, .390, 19.1, 10.3, 573., 20.8 .734, 1.16, 22.8, 17.4, 606., 37.5 9.24, 9.34, 23.6, 13.9, .113e4, 33.0 15.9, 88.2, 260., 154., .599e4, 817. .373, 1.81, 39.3, 23.7, 740., 43.5 .410e-1, .316e-1, .139, .751e-1, 108., .229 .347e-1, .191e-1, .645e-1, .293e-1, 34.9, .755e-1 .267e-2, .111e-2, .156e-2, .207e-2, 29.2, .243e-1 .691e-3, .191e-5, .332e-3, .206e-2, 17.1, .244e-1 .528e-3, .493e-6, .402e-3, .490e-2, 34.7, .188 .142e-2, .291e-4, .521e-4, .556e-3, 32.7, .619e-1任意地區自身排放高錳酸鹽和氮氨排放量分別為即為某地區高錳酸鹽和氮氨總量減去上游相鄰城市經過降解后剩下即可，分別如下：，利用matlab解得：如下所示：時間A1A2A3A4A5A6A720040484874828146197844678446778170446552004051599628791593997379092373509062094200406100254257377134788877529741590660042004071118454076726191.01E+056999744380575102004082169221146742901.03E+05640325400059278200409383082.08E+051.79E+052.07E+051.36E+05731201.77E+05200410260825916706577741746890577007242220041142001633820140311993211936320352762004121521.68514.31776037449303602330129799200501854.45226868341769261603164636000200502550.86485.4902025536319303658334730200503685.390061087828864257406450238698200504706.273001296020996302004599333149平均值8986.02337077.950708.2370086.4657375.7745568.8557428.85如下所示：時間A1A2A3A4A5A6A7200404553.52896.85759.86654.47586.66078592.47200405260.43536.55939.965596555.480953100.4200406160.42555.66901.67890.94414.748841965.2200407186.42951.33630.88653.24302.664121567.520040837401792.84550.78523.47690.142141722200409565.27600.115427158801804643985663200410260.83725.64391.27988.74688.345301492.5200411902205.925444319.93211.91932979.7720041276.083340.511843316.915841939.11489.920050149.842170.81005.43449.822892430.8316820050291.81981.6811.83431.413392448.94681200503161.982559.6984.22956.815734887.3429.8120050464.22299.57023040.821142997.34419.9平均值481.58463047.4314140.9546358.8625030.3544249.7232405.496通過比較平均值的大小，得出高錳酸鹽超標的污染A4污染最嚴重，A7次之，A5為第三，A1最輕。氨氮超標的污染A4最嚴重，A5次之，A3為第三，A1最輕。（三）問題三模型的建立和求解：GM (1, 1) 模型：灰色系統理論的實質是將無規律的原始數據進行累加生成，得到規律性較強的生成數列后再重新建模。由生成模型得到的數據再通過累加生成的逆運算累減生成得到還原模型，再還原模型作為預測模型。灰色模型是預測工作的基礎模型。對于采集的原始數據序列:記為原始序列，為由經過一次累加生成的序列，其中，表示的均值生成序列， 表5通過運用灰色原理，GM模型解得關聯度e =0.9993，灰色關聯度精度等級為1級（好），畫出圖像。最終得出未來10年的值分別為：4.3314 ，4.2580，4.1857，4.1147，4.0449，3.9763，3.9089，3.8425，3.7774，3.7133。從整體分析，如果不加治理，長江水質在未來10年會持續下降。（四）問題四模型的建立和求解：取19952004年文年水的IV和V類水質分別所占百分比和長江的年總流量作為預測原始值，預測出20052014年的IV和V類水質分別所占百分比和長江的年總流量。長江污水年排放量與長江的年總流量的比作為因變量Y，利用spss 軟件的線性回歸分析功能對19952004 年的每年長江干流水文年各類水質所占河長的百分比做為歷史數據進行回歸分析，得到各個因子的常量系數，分別為系統自動剔除第三變量，從而得到：其中由此得出，且未來的年處理污水量=期望長江年總流量年份期望長江總流量IV類河長所占百分比% ()V類河長所占百分比% ()需要處理的污水量（億噸）2005102657.651.570.004344.13952006101908.351.950.005354.0072007101169.112.400.006565.7542008100429.932.970.008181.34022009996910.833.670.0199.692010989611.824.530.012118.7522011982412.895.600.015147.362012975314.066.920.019185.3072013968215.348.540.023222.6862014961116.7310.550.029278.719表6 每年需要處理的污水量（五）問題五對解決長江水質污染問題的一些切實可行的建議和意見：（1）對保護長江立法，讓人們的意識受到法律的約束。一方面可以制定促進循環經濟發展的政策和法律法規；另一方面，還可以加大對違法排污行為的處罰力度。（2）對于生活污水也需要進行一定的控制。喚醒民眾、領導干部的環保意識。由于長江兩岸農田過量噴灑的農藥被雨水沖刷到河里，出現了氨氮含量過高、河水富營養化、藍藻瘋長等情況，嚴重破壞了水資源。迫切需要提高人民大眾的環保意識，在保證收成的情況下，盡量少使用農藥毋將農藥水隨便排放。同時需要政府撥款在長江沿岸多建水凈化廠，將需要排放到長江里的水進行凈化處理后再排放。（3）在長江上游地區植樹造林，據統計長江上游地區森林覆蓋率曾達到60%至80%，現如今卻只有5%至7%，因此每年流入長江泥沙近24億噸，長江源頭出現沙漠化及斷流，污染還殃及近海。2003年出現赤潮近40次，赤潮導致“氧虧”，進而影響漁業資源及品質，水生生態系統瀕于崩潰。因此應該鼓勵在長江上游地區植樹造林。（4）對于長江沿線的鋼鐵廠、化工廠、造紙廠、造船廠、拆船廠這樣的重污染企業禁止直接將污水排入長江。嚴格要求工廠將污水經過處理設備處理后，再對各種污染物含量進行檢測，檢測所要排放的水屬于可飲用水，才可排放。若發現工廠將檢驗未達標水排入長江則處以罰款或者是責令他們拆遷工廠。六、模型評價與推廣模型中運用了灰色模型預測法，在預測實踐中以中短期預測的效果為佳，弱化了預測對象的隨機性變化，突出了趨勢性。忽略了次要因素對水質的影響，簡潔、易懂，具有一定的合理性和實用性。模型還考慮到了變量剔除問題，正確的選擇了變量，減少異常點對模型結果的影響。但是模型對于水質的凈化只考慮了自然降解，而對于實際情況比如是否存在污水處理廠等認為凈化設施的情況沒有考慮。對于降解系數選取，認為全流域相等，與實際情況存在一定的誤差。對于污水在監測站之后馬上匯入江水的假設，與實際沿江不規律分布的污水排放不符，也會產生誤差。實際情況下，污染物排放量與諸多因素有關，僅僅估計其與時間的關系會造成偏差。用給出的前十年的水質數據來預測后十年的水質變化情況，數據不夠充足，容易造成誤差。七、參考文獻 1 鄧聚龍1 灰色預測與決策M 1 武漢: 華中理工大學出版社，1986.1 2 姜啟源，謝金星，葉俊，數學模型，北京：高等教育出版社，2004.4八、附件第二題程序1clear;clc;a1=2.9 2.8 3.0 3.2 2.95 4.95 2.7 2.3 2.3 1.8 1.5 1.4 1.65;a2=1.5 1.35 1.65 1.85 2.75 3.25 1.9 1.3 1.2 1.0 0.7 0.65 0.8;a3=0.9 0.85 1.25 1.45 1.4 1.8 1.55 0.75 0.75 0.55 0.5 0.5 0.45;a4=0.95 1.0 1.4 1.5 1.45 2.0 1.6 0.85 0.8 0.65 0.65 0.7 0.6;a5=1.05 1.1 1.5 1.55 1.6 2.75 1.65 0.9 0.8 0.7 0.7 0.8 0.75;a6=1.15 1.15 1.55 1.65 1.7 3.4 1.8 0.95 0.85 0.75 0.75 0.85 0.8;b1=950000;b2=778000;b3=395000;b4=500000;b5=164000;b6=464000;c1=b1./(86400.*a1);c2=b2./(86400.*a2);c3=b3./(86400.*a3);c4=b4./(86400.*a4);c5=b5./(86400.*a5);c6=b6./(86400.*a6);c1=vpa(c1,3)c2=vpa(c2,3)c3=vpa(c3,3)c4=vpa(c4,3)c5=vpa(c5,3)c6=vpa(c6,3)運行結果：第二題程序2clear;clc;a1=8487 48300 46200 84480 81490 82600;. 15996 28820 59400 73800 92380 10800;. 10025 42600 77140 79100 75400 51360;. 11184 54120 72640 101220 70200 54230;. 21692 21200 74400 103600 64200 63460;. 38308 209440 181900 209820 138320 118720;. 2608 25920 70670 78050 47120 65100;. 4200 16340 20140 31200 32120 37400;. 1521.6 8515 17760 37450 30360 23970;. 854.4 5226 8683 41769 26160 31980;. 550.8 6485.4 9020 25536 31930 36990;. 685.3 9006 10878 28864 25740 65070;. 706.2 7300 12960 20996 30200 46460;a2=553.5 2898 5760 6656 7587 6490;. 260.4 3537 5940 6560 6556 8500;. 160.4 2556 6902 7910 4425 5457;. 186.4 2952 3632 8676 4320 7018;. 3740 1802 4560 8547 7704 5344;. 565.2 7616 15515 16140 18200 10388;. 260.8 3726 4393 8028 4712 5270;. 90 2205.9 2544 4320 3212 2040;. 76.08 3340.5 1184 3317 1584 1974;. 49.84 2170.8 1005.4 3449.8 2289 2460;. 91.8 1981.65 811.8 3431.4 1339 2466;. 161.98 2559.6 984.2 2956.8 1573 4922;. 64.2 2299.5 702 3040.8 2114 3030;t=3.79 6.00 5.08 6.09 1.81 4.67;. 3.93 6.67 5.38 5.79 1.73 4.67;. 3.67 5.46 3.66 4.13 1.27 3.46;. 3.44 4.87 3.15 3.86 1.22 3.25;. 3.73 3.27 3.27 3.99 1.19 3.16;. 2.22 2.77 2.54 2.89 0.69 1.58;. 4.07 4.74 2.93 3.62 1.15 2.98;. 4.78 6.93 6.10 6.81 2.11 5.65;. 4.78 7.50 6.10 7.23 2.37 6.32;. 6.11 9.00 8.31 8.90 2.71 7.16;. 7.33 12.9 9.14 8.90 2.71 7.16;. 7.85 13.9 9.14 8.27 2.37 6.32;. 6.66 11.3 10.2 9.64 2.59 6.71;x=0.2.t;y1=a1.*x;y2=a2.*x;Q1=vpa(y1,3)Q2=vpa(y2,3)運行結果：Q1 = 19.0, 3.09, 13.0, 4.68, .443e4, 45.0 28.6, .627, 10.3, 6.62, .571e4, 5.88 27.3, 6.50, 213., 103., .977e4, 196. 44.1, 21.3, 456., 203., .985e4, 290. 53.6, 110., 385., 168., .946e4, 392. .108e4, .243e4, .305e4, .200e4, .456e5, .934e4 3.73, 12.6, 633., 230., .740e4, 538. 1.92, .234, 1.10, .542, .108e4, 4.20 .694, .487e-1, .968, .331, 669., .917 .458e-1, .268e-2, .135e-1, .251e-1, 334., .316 .415e-2, .624e-5, .369e-2, .154e-1, 407., .366 .223e-2, .173e-5, .445e-2, .478e-1, 568., 2.49 .156e-1, .922e-4, .962e-3, .384e-2, 467., .948 Q2 = 1.24, .185, 1.62, .369, 412., 3.53 .466, .770e-1, 1.03, .589, 405., 4.63 .436, .390, 19.1, 10.3, 573., 20.8 .734, 1.16, 22.8, 17.4, 606., 37.5 9.24, 9.34, 23.6, 13.9, .113e4, 33.0 15.9, 88.2, 260., 154., .599e4, 817. .373, 1.81, 39.3, 23.7, 740., 43.5 .410e-1, .316e-1, .139, .751e-1, 108., .229 .347e-1, .191e-1, .645e-1, .293e-1, 34.9, .755e-1 .267e-2, .111e-2, .156e-2, .207e-2, 29.2, .243e-1 .691e-3, .191e-5, .332e-3, .206e-2, 17.1, .244e-1 .528e-3, .493e-6, .402e-3, .490e-2, 34.7, .188 .142e-2, .291e-4, .521e-4, .556e-3, 32.7, .619e-1第二題程序3clear;clc;a1=48300 46200 84480 81490 82600 44700;. 28820 59400 73800 92380 10800 62100;. 42600 77140 79100 75400 51360 66200;. 54120 72640 101220 70200 54230 57800;. 21200 74400 103600 64200 63460 59670;. 209440 181900 209820 138320 118720 186300;. 25920 70670 78050 47120 65100 72960;. 16340 20140 31200 32120 37400 35280;. 8515 17760 37450 30360 23970 29800;. 5226 8683 41769 26160 31980 36000;. 6485.4 9020 25536 31930 36990 34730;. 9006 10878 28864 25740 65070 38700;. 7300 12960 20996 30200 46460 33150;a2=2898 5760 6656 7587 6490 596;. 3537 5940 6560 6556 8500 3105;. 2556 6902 7910 4425 5457 1986;. 2952 3632 8676 4320 7018 1605;. 1802 4560 8547 7704 5344 1755;. 7616 15515 16140 18200 10388 6480;. 3726 4393 8028 4712 5270 1536;. 2205.9 2544 4320 3212 2040 980;. 3340.5 1184 3317 1584 1974 1490;. 2170.8 1005.4 3449.8 2289 2460 3168;. 1981.65 811.8 3431.4 1339 2466 4681;. 2559.6 984.2 2956.8 1573 4922 430;. 2299.5 702 3040.8 2114 3030 4420;x1=19.0, 3.09, 13.0, 4.68, .443e4, 45.0;. 28.6, .627, 10.3, 6.62, .571e4, 5.88;. 27.3, 6.50, 213., 103., .977e4, 196.;. 44.1, 21.3, 456., 203., .985e4, 290.;. 53.6, 110., 385., 168., .946e4, 392.;. .108e4, .243e4, .305e4, .200e4, .456e5, .934e4;. 3.73, 12.6, 633., 230., .740e4, 538.;. 1.92, .234, 1.10, .542, .108e4, 4.20;. .694, .487e-1, .968, .331, 669., .917;. .458e-1, .268e-2, .135e-1, .251e-1, 334., .316;. .415e-2, .624e-5, .369e-2, .154e-1, 407., .366;. .223e-2, .173e-5, .445e-2, .478e-1, 568., 2.49;. .156e-1, .922e-4, .962e-3, .384e-2, 467., .948;x2 = 1.24, .185, 1.62, .369, 412., 3.53;. .466, .770e-1, 1.03, .589, 405., 4.63;. .436, .390, 19.1, 10.3, 573., 20.8;. .734, 1.16, 22.8, 17.4, 606., 37.5;. 9.24, 9.34, 23.6, 13.9, .113e4, 33.0;. 15.9, 88.2, 260., 154., .599e4, 817.;. .373, 1.81, 39.3, 23.7, 740., 43.5;. .410e-1, .316e-1, .139, .751e-1, 108., .229;. .347e-1, .191e-1, .645e-1, .293e-1, 34.9, .755e-1;. .267e-2, .111e-2, .156e-2, .207e-2, 29.2, .243e-1;. .691e-3, .191e-5, .332e-3, .206e-2, 17.1, .244e-1;. .528e-3, .493e-6, .402e-3, .490e-2, 34.7, .188;. .142e-2, .291e-4, .521e-4, .556e-3, 3