基于回归分析应用的长江水质预测与控制.docVIP

基于回归分析的长江水质预测与控制 2013年12月22日星期日 摘要 水资源的保护问题是当今全球范围内讨论的焦点问题之一。我们使用eviews和matlab做回归模型对长江后10年废水排放总量预测及模拟并且得出了长江水质变化发展的特点:l)长江近10年主要特点: Ⅱ类，Ⅲ类水具有缓冲性。2)长江后10年主要特点: Ⅰ类，劣V类水具有两极性。然后我们、线性回归模型在预测长江水质发展趋势，最后得出了对长江未来10年水质发展趋势估计(2005、2014)不可饮用水占总流量的百分比(%)，通过对回归模型进行控制，保证了一年中长江劣V类水消失的概率为95%，并且使得长江干流的Ⅳ类和V类水的比例在20%以内。 关键字：回归方程 百分比 eviews matlab 问题三 为了在宏观上正确地把握长江未来水质污染的发展趋势，本文依照过去10年的主要统计数据首先对长江的水质特性作了一些简单的探索。我们取1995~2004年水文年全流域水占河长的百分比进行研究。 表1 1994~2004年长江全流域个类水百分比 年份T Ⅰ类Y1 Ⅱ类Y2 Ⅲ类Y3 Ⅳ类Y4 Ⅴ类Y5 劣Ⅴ类Y6 百分比 百分比 百分比 百分比 百分比 百分比 1995 25.8 42.6 24.7 3.9 3 0 1996 15.3 20.2 49.8 9.7 1.9 3.1 1997 12.2 24.9 43.6 13.3 2.6 3.4 1998 11.5 24.1 52.8 8.3 1.7 1.6 1999 5.2 39.8 35.2 9.5 6.2 4.1 2000 5.6 32.8 35.6 16.6 4.4 5.3 2001 5.9 33.1 34.7 14 5.5 6.8 2002 4.4 44 28.3 10 3.2 10 2003 4.7 41.5 31.3 6.4 5.8 10.3 2004 1.2 26.9 39.9 14.8 5.9 11.3 图1 图2 过观察统计一图1和图2，可以发现污染类Ⅱ与Ⅲ的线图呈现出此消彼长的走向并且很像一只蝴蝶结，具有极强的共线性，如果把两者相合并就会显示出非常平稳的走势，我们称之为“缓冲区”。进一步地，将Ⅱ与Ⅲ类合并，Ⅳ与Ⅴ类合并后作图（见附件），则可以看出，I类水质与劣V类水质，合并后的Ⅱ、Ⅲ类与合并后的Ⅳ、V类水质同样显示出了比较明显的共线性。其中，缓冲区的作用在过去10年中至关重要。随着时间的推移，废水排放量的增多，在近几年中I类水锐减，到2004年几乎消失，而缓冲区中的水(即Ⅱ类和Ⅲ类水)却保持相对较平缓的变化趋势，与此同时，Ⅳ类水及劣Ⅴ类水百分比会急剧升高，Ⅴ类水相对比较平缓。即优质水消失，污染严重的水剧增，这就是我们说的“两极性”。 对表1数据进行拟合，通过在Matlab中编程得到如下方程和图3中的曲线。各类水占的河长百分比y与时间t关系分别为： 图3 通过计算得到未来十年各类水占河长的百分比为下表： 表2 2005~2014各类水占河长的百分比 水的类别 Ⅰ类（%） Ⅱ类（%） Ⅲ类（%） Ⅳ类（%） Ⅴ类（%） 劣Ⅴ类（%） 水文年 2005 4.07 32.59 38.71 12.35 4.64 7.66 2006 3.88 32.58 38.75 12.42 4.67 7.74 2007 3.71 32.57 38.79 12.47 4.69 7.8 2008 3.57 32.56 38.82 12.52 4.70 7.86 2009 3.45 32.55 38.84 12.56 4.72 7.91 2010 3.35 32.54 38.87 12.59 4.73 7.95 2011 3.26 32.53 38.89 12.62 7.74 7.99 2012 3.17 32.53 38.90 12.65 4.75 8.02 2013 3.10 32.52 38.92 12.67 4.76 8.05 2014 3.03 32.51 38.93 12.70 4.77 8.08 由此，我们得到如下结论: l)长江10年内主要特点: Ⅱ类，Ⅲ类水具有缓冲性。 2)长江后10年主要特点:I类，劣Ⅴ类水具有两极性。 检验：通过表2和图三，我们发现未来十年，Ⅰ类、Ⅱ类水占河长的百分比呈逐年下降趋势，而其它类的水占河长的百分比则呈逐年上升趋势如果不采取措施长江水质将趋于恶化。 图3的模型中的各类水占河长百分比是平滑曲线，只是代表大致趋向，而实际上其受天气、水温等影响是成波动态势，这也是本模型的缺点。 问题4 Zi(0i11)：过去十年污水年排放量; Zj(10j21)：未来十