摘要
随着应用领域的日益广泛,现有的灰色预测理论与方法往往难以解决实际应用中不断出现的新问题。本文按照“提出问题、分析问题、解决问题”的逻辑思路,主要研究含可变参数的缓冲算子和GM(1,1)幂模型的理论与应用问题,以期进一步完善灰色预测理论,扩大灰色预测理论与方法的应用范围。本文的主要工作体现在以下几个方面:
1.提出了缓冲算子光滑性的概念,给出了缓冲算子能否提高序列光滑性的判别条件,在此基础上,分别证明了弱化缓冲算子和强化缓冲算子的光滑性特点,并通过经典算例加以验证。从提高建模序列光滑性和消除冲击扰动因素两个方面,揭示了缓冲算子提高灰色模型对冲击扰动系统预测精度的原因。
2.针对传统缓冲算子不能实现作用强度的微调,从而导致缓冲算子的作用效果过强或过弱的问题,将可变参数引入到缓冲算子的构造中,分别提出了两类含可变参数的缓冲算子:幂缓冲算子和变权缓冲算子。构造了若干实用的幂弱化缓冲算子和幂强化缓冲算子,研究了各类幂缓冲算子之间的关系,并分析了幂缓冲算子的参数与其作用强度的定性关系;研究了变权缓冲算子对原始序列作用强度的定量测算与有效控制的问题,从根本上解决了冲击扰动系统的预测过程中常常出现的定量预测结果与定性分析结论不符的问题。在此基础上,提出了“单调性不变”公理,进一步完善了缓冲算子的公理体系。
3.提出了GM(1,1)幂模型幂指数的白化方法,给出了其白化微分方法的求解过程,并分析了GM(1,1)幂模型解的性质;基于矩阵求逆条件数的定义,分三种情形研究了GM(1,1)幂模型参数估计过程可能出现的矩阵病态性问题,总结了影响GM(1,1)幂模型病态性的主要影响因素。基于非线性优化方法,依次对GM(1,1)幂模型的幂指数、背景值插值系数,以及初始条件进行优化,并通过实例验证了本文的优化效果,实现了灰色预测方法对小样本振荡序列的高精度预测。此外,基于GM(1,1)幂模型的基本定义,本文提出了GM(1,1)幂模型的五类派生模型,进一步完善了灰色系统幂模型体系。
4.基于GM(1,1)幂模型的误差分析,提出了无偏GM(1,1)幂模型,并从理论上和实例应用两个方面证明了该模型的无偏性;在此基础上,考虑参数β1和β2已知和未知两种情形下,分别研究了无偏GM(1,1)幂模型初始条件的优化问题,并通过实例验证模型的有效性。
5.以本文提出的变权缓冲算子和优化的GM(1,1)幂模型为基本工具,研究了我国31个省、市、自治区的工业废水排放达标率的预测和预警问题,并提出了相应的政策建议。
With the wider application of Grey forecast theory, there are increasingly many new issues that need to be tackled with this theory. Based on ideas of“questions raising, question analysis, and question solution”, this paper focuses on the theories and application of parameter-contained buffer operators and GM(1,1) exponential model. This research aims at improving Grey Systems Theory as well as expanding its application. The study can be concluded as the following parts.
1. The smoothness of buffer operator is defined by initiating criterion of improving the smoothness of series as well as analyzing the smoothness of weakening buffer operator and reinforcement buffer operator. The reason of buffer operator to improve the forecast accuracy of shock disturbed system has been revealed in aspects of improving smoothness of series and eliminating shock disturbed factors.
2. It is known that function strength cannot be slightly adjusted in traditional buffer operator so that its function effect is hard to co-ordinate. To solve this problem, the variable parameter is introduced in the construction of buffer operators. Such operators can be divided into two categories, exponential buffer operators and weight variable buffer operators. Several practical exponential weakening and reinforcement buffer operators are constructed. In addition, the relation of various kinds of buffer operators has been studied as well as the quantificational relation between parameters and their function strength in exponential buffer operators. Besides, the function strength is measured quantitatively and effectively controlled in weight variable buffer operator to original series. This can perfectly solve the discrepancy of results in quantitative forecast and qualitative analysis in forecast process of shock disturbed system. Based on this, the axiom of“monotony invariance”is put forward to further improve the axiom system of buffer operators.
3. The whitenization methods of power exponent in GM(1,1) Power Model is put forward with the solution in whitenized differential method and properties of solution. Based on definition of condition matrix in inverse problem, the issues of possible ill-condition of matrix are studied in the parameter estimation. The main factors affecting the ill-condition of GM(1,1) Power Model are concluded. Based on non-linear optimization methods, the power exponents, interpolation coefficient of background value, and initial condition are optimized. The optimal results are verified in case study to achieve high accuracy in fluctuant series with small sample. In addition, based on the basic definition of GM(1,1) Power Model, five derived models are put forward to perfect the system of Grey Exponential Models.
4. Based on the error analysis of GM(1,1) Power Model, unbiased GM(1,1) Power Model is put forward. Furthermore, the unbiasedness of model is proved in theory and in case study. Based on this, the optimization of initial condition in unbiased GM(1,1) Power Model is analyzed separately taking account of given parametersβ1 andβ2. The effectiveness of model is studied in case study.
5. The forecast and warning of water discharge rate in 31 provinces and administrative regions are studied by method of weight variable buffer operator and optimized GM(1,1) Power Model, followed by policy analysis.
引文
[1] Deng J L. The Control problem of grey systems [J]. System Control Letter, 1982, 1(5): 288~294.
[2]徐伟宣,李建平.我国管理科学与工程学科的新进展[J].中国科学院院刊,2008,23(2) : 162~167.
[3] Wang J, Yan R, Hollister K, Zhu D. Ahistoric reviewof management science research in China [J]. Omega, 2008, (36) :919~932.
[4]邓聚龙.灰理论基础[M].武汉:华中科技大学出版社,2002.
[5]邓聚龙.灰预测与灰决策[M].武汉:华中科技大学出版社,2002.
[6]邓聚龙.灰色数理资源科学导论[M].武汉:华中科技大学出版社,2007.
[7] Liu S F, Lin Y. Grey Information Theory and Practical Applications [M]. Springer-Verlag, London, 2006.
[8]刘思峰,党耀国,方志耕.灰色系统理论及其应用[M].北京:科学出版社.2004.
[9]党耀国,刘思峰,王正新,林益.灰色预测与决策模型研究[M].科学出版社,2009.
[10]肖新平,宋中民,李峰.灰技术基础及其应用[M].北京:科学出版社,2005.
[11]张岐山.灰朦胧集的差异性信息理论[M].北京:石油工业出版社,2002.
[12]王清印.灰色数学基础[M].武汉:华中理工大学出版社,1996.
[13] Wen K L. Grey systems: modeling and prediction [M]. Arizona: Yang’s, 2004.
[14]邓聚龙.灰色系统理论教程[M].武汉:华中理工大学出版社,1990.
[15]陈涛捷.灰色预测模型的一种拓广[J].系统工程,1990,8(4):50~52.
[16]李群.灰色预测模型的进一步拓广[J].系统工程理论与实践,1993,13(1):64~66.
[17]黄福勇.灰色建模的变换方法[J].系统工程理论与实践,1996,14(16):35~38.
[18]何斌,蒙清.灰色预测模型的拓广方法研究[J].系统工程理论与实践,2002,22(9):138~141.
[19]戴文战,李俊峰.基于函数x ?α的GM(1,1)模型及我国农村人均住房面积建模中的应用[J].系统工程理论与实践,2004,24(11):63~67.
[20]李翠凤,戴文战.基于函数cot x变换的灰色建模方法[J].系统工程,2005,23(3):110~114.
[21] Lin Y H, Lee P C. Novel high-precision grey forecasting model [J]. Automation in Construction, 2007, (16):771~77.
[22] Li G D, Yamaguchi D and Nagai M. New methods and accuracy improvement of GM according to Laplace Transform. The Journal of Grey System. 2005, 8(1):13~24.
[23]钱吴永.数据变换技术与GM(1,1)模型研究[D].南京:南京航空航天大学,2009.
[24] Wei Y, Zhang Y. A Criterion of Comparing the Function Transformations to Raise the SmoothDegree of Grey Modeling Data [J]. The Journal of Grey System.2007, 19 (1):91~98.
[25] Wei Y, Zhang Y. An Essential Characteristic of the Discrete Function Transformation to Increase the Smooth Degree of Data [J]. The Journal of Grey System, 2007,19(3): 293~300.
[26] Wei Y, Zheng F. Studying grey modeling properties in special series transformation cot x , x ?α(α> 0) and cos x [J]. The Journal of Grey System.2008, 20(1): 79~86.
[27] Wang S H, Wei Y. One kind of effective discrete transformation on converting the high growth series into the growth one [J]. Journal of System Science and Information, 2008, (2):181~188.
[28]魏勇,胡大红.光滑性条件的缺陷及其弥补办法[J].系统工程理论与实践,2009,29(8):165~170.
[29]刘思峰.冲击扰动系统预测陷阱与缓冲算子[J].华中理工大学学报,1997, 25(1):25~27.
[30]党耀国,刘思峰,刘斌等.关于弱化缓冲算子的研究[J].中国管理科学, 2004, 12 (2) : 108~111.
[31]吴正朋,刘思峰,米传民等.基于反向累积法的弱化缓冲算子序列研究[J].中国管理科学,2009,17 (3) : 136~141.
[32]吴正朋,刘思峰,崔立志.基于不动点的新弱化缓冲算子的研究[J].控制与决策,2009,24(12):1805~1809.
[33]刘以安,陈松灿,张明俊等.缓冲算子及数据融合技术在目标跟踪中的应用[J].应用科学学报,2006,24(2):154~158.
[34]尹春华,顾培亮.基于灰色序列生成中缓冲算子的能源预测[J].系统工程学报, 2003,18(2):189~192.
[35]党耀国,刘斌,关叶青.关于强化缓冲算子的研究[J].控制与决策,2005,20(12):1332~1336.
[36]党耀国,刘思峰,米传民.强化缓冲算子性质的研究[J].控制与决策,2007,22(7):730~734.
[37]吴正朋,刘思峰,米传民等.基于单调函数的若干实用强化缓冲算子的构造[J].系统工程,2009,27(5):124~128.
[38]谢乃明,刘思峰.强化缓冲算子的性质与若干实用强化算子的构造[J].统计与决策,2006,(4): 9~10.
[39]崔杰,党耀国.一类新的SBO及其在GM(1,1)中的应用研究[J].管理工程学报,2009,(4):153~156.
[40]崔杰,党耀国,刘思峰等.一类新的强化缓冲算子及其数值仿真[J].中国工程科学,2010,(2):108~112.
[41]魏勇,孔新海.几类强弱缓冲算子的构造方法及其内在联系[J].控制与决策,2010,25(2):196~202.
[42]关叶青,刘思峰.基于不动点的强化缓冲算子序列及其应用[J].控制与决策,2007,22(10):1189~1192.
[43]王正新,党耀国,刘思峰.变权缓冲算子及缓冲算子公理的补充[J].系统工程,2009,27(1):113~117.
[44]王正新,党耀国,刘思峰.变权缓冲算子及其作用强度的研究[J].控制与决策,2009,24(8):1218~1222.
[45]张庆,刘思峰,王正新.几何变权缓冲算子及其作用强度研究[J].系统工程,2009,27(11):105~108.
[46] Carlin J B, Dempster A P, Jonas A B. On models and methods for Bayesian time series analysis [J].Journal of Econometrics, 1985(30): 67 ~90.
[47] Ben-Haim Y. Info-gap forecasting and the advantage of sub-optimal models [J]. European Journal of Operational Research, 2009, 197: 203 ~213.
[48]刘思峰,邓聚龙. GM(1,1)模型的适用范围[J].系统工程理论与实践,2000,20(5):121 ~124.
[49] Li X C. On parameters in Grey Model GM(1,1) [J]. The Journal of Grey System, 1998, 10(2):155 ~162.
[50]李留藏,许文天.灰色系统GM(1,1)模型的讨论[J].数学的实践与认识,1993,(1):15~22.
[51]李炳乾.原始序列特征与灰色模型的选择[J].灰色系统理论与实践,1992,2(1):63~66.
[52]胡大红,魏勇.灰模型对单调递减序列的适应性与参数估计[J].系统工程与电子技术,2008,30(11):63~66.
[53] Kuo C Y, Ching J L, Yen T H. Generalized admissible region of class ratio for GM(1,1) [J]. The Journal of Grey System, 2000, 12(2):153~156.
[54] Yeh J F, Lu H C. On some of the basic feature of GM(1,1) model [J]. The Journal of Grey System, 1996, 8(1):19~36.
[55]吉培荣,黄巍松,胡翔勇.无偏灰色预测模型[J].系统工程与电子技术,2000,(6):6~8.
[56]王文平,邓聚龙.灰色系统中GM(1,1)模型的混沌特性研究[J].系统工程,1997,15(2):13~15.
[57]王正新,党耀国,刘思峰.无偏GM(1,1)模型的混沌特性分析[J].系统工程理论与实践, 2007,27(11):153~158.
[58]郑照宁,武玉英,包涵龄. GM模型的病态性问题[J].中国管理科学,2001,(5):38~44
[59]党耀国,王正新,刘思峰.灰色模型的病态问题研究[J].系统工程理论与实践,2008,28(1):156~160.
[60]王正新. GM(1,1)模型的特性与优化研究[D].南京:南京航空航天大学,2007.
[61]张彤,王子才,吴建伟. GM(1,1)模型参数的改进计算方法[J].系统工程与电子技术,1998, 20(8):58~60.
[62]王义闹,刘开第,李应川.优化灰导数白化值的GM(1,1)建模法[J].系统工程理论与实践,2001,21(05):124~128 .
[63] Wang Y N, Chen Z J , Gao Z Q , et al. A generalization of the GM(1 ,1) direct modeling method with a step by step optimizing grey derivative’s whiten values and it s applications [J]. Kybernetes, 2004, 33 (2): 382~389.
[64]穆勇.优化灰导数白化值的无偏灰色GM(1,1)模型[J].数学的实践与认识,2003, (3):13~16.
[65]李玻,魏勇.优化灰导数后的新GM(1,1)模型[J].系统工程理论与实践,2009,29(2):100~105.
[66]张凌霜,王丰效.逐步优化灰导数的非等间距GM(1,1)模型[J].数学的实践与认识,2010, (11):63~67.
[67]谭冠军. GM(1,1)模型的背景值构造方法和应用(Ⅰ)[J].系统工程理论与实践,2000,20(5):70~97.
[68]李俊蜂,戴文战.基于插值和Newton-Cotes公式的GM(1,1)模型的背景值构造新方法及应用[J].系统工程理论与实践,2002,22(10):122~126.
[69]唐万梅,向长合.基于二次插值的GM(1,1)模型预测方法的改进[J].中国管理科学,2006,14(6):109~112.
[70]向跃霖.废气排放量灰色建模初探[J].环境科学研究,1995,8(6):45~48.
[71]谢开贵,李春燕.基与遗传算法的GM(1,1,λ)模型[J].系统工程学报, 2000,15(2): 168~172.
[72]王义闹,刘光珍,刘开第. GM(1,1)模型的一种逐步优化直接建模法[J].系统工程理论与实践,2000,20(9):99~104.
[73]罗党,刘思峰,党耀国.灰色模型GM(1,1)优化[J].中国工程科学,2003,5(8):50~53.
[74] Lin Y H, Lee P C, Chang T P. Adaptive and high-precision grey forecasting model [J]. Expert Systems with Applications, 2009, (36): 9658~9662.
[75]王正新,党耀国,刘思峰.基于离散指数函数优化的GM(1,1)模型[J].系统工程理论与实践,2008,28(2):61~67.
[76]穆勇.灰色预测模型参数估计的优化方法[J].青岛大学学报,2003,16(3):95~98.
[77] Xiao X P. A New Modified GM(1,1) Model: Grey optimization Model [J]. Journal of Systems Engineering and Electronics, 2001, 12(2): 1~5.
[78]何文章,宋国乡,吴爱弟.估计GM(1,1)模型中参数的一族算法[J].系统工程理论与实践,2005,25(1):68~75.
[79] Wang C H, Hsu L C. Using genetic algorithms grey theory to forecast high technology industrial output [J]. Applied Mathematics and Computation, 2008, (195) : 256~263.
[80]张岐山.提高灰色GM(1,1)模型精度的微粒群方法[J].中国管理科学, 2007,15(5):126~129.
[81] Hsu C C, Chen C Y. A modified grey forecasting model for long-term prediction [J]. Journal of The Chinese Institute of Engineers, 2003, 26(3): 301~308.
[82] Shih N Y, Shih N J, Yu W C. The Parameter Estimation of Time-varying GM(1,1) [J]. Journal of Grey System, 2006, 9(1):51~56.
[83]王义闹,吴利丰.基于平均相对误差绝对值最小的GM(1 ,1)建模[J].华中科技大学学报(自然科学版),2009,37(10):29~31.
[84] Dang Y G, Liu S F. The GM Models That x(n) be Taken as Initial Value [J]. Kybernetes, 2004, 33(2):247~255.
[85]刘斌,刘思峰,翟振杰等. GM(1,1)模型时间响应函数的最优化[J].中国管理科学, 2003,(4):54~57.
[86] Ren H B, Dang Y G, Wang Z X. GM(1,1) Model of Time Sequence Coefficient and Application [C]. Proceedings of 2007 IEEE International Conference on Grey Systems and Intelligent Services, 2007, 18~20.
[87] Wang Y H, Dang Y G, Li Y Q, Liu S F. An approach to increase prediction precision of GM(1,1) model based on optimization of the initial condition [J]. Expert Systems with Applications, 2010, (37): 5640~5644.
[88]宋中民,同小军,肖新平.中心逼近式灰色GM(1,1)模型[J].系统工程理论与实践,2001,21(5):110 ~113.
[89] C.K. Chen, T.L. Tien, The indirect measurement of tensile strength by the deterministic grey dynamic model DGDM(1, 1, 1) [J]. International Journal of Systems Science, 1997, 28 (7) :683~690.
[90] C.K. Chen, T.L. Tien, A new forecasting method of discrete dynamic system[J]. Applied Mathematics and Computation, 1997, 86 (1): 61~84.
[91] Hsu C C, Chen C Y. Applications of improved grey prediction model for power demand forecasting [J]. Energy Conversion and Management, 2003, 44(14):2241~2249.
[92]谢乃明,刘思峰.离散GM(1,1)模型与灰色预测模型建模机理[J].系统工程理论与实践,2005,25(1):93~98.
[93] Xie N M, Liu S F. Discrete Grey Forecasting Model and its Optimization. Applied Mathematical Modelling, 2009, (33):1173~1186.
[94]姚天祥,刘思峰,党耀国.初始值优化的离散灰色预测模型[J].系统工程与电子技术,2009,31(10):2394~2398.
[95] Wang Z X, Dang YG, Liu B. Recursive Solution and Approximating Optimization to Grey Models with High Precision. Journal of Grey System, 2009, 21(2): 185~194.
[96]王正新,党耀国,刘思峰.两阶段灰色模型及其应用[J].系统工程理论与实践,2008,28(11):109~114.
[97]钱吴永,党耀国.基于振荡序列的GM(1,1)模型[J].系统工程理论与实践2009,29(3): 149~154.
[98]罗佑新,周继荣.非等间距GM(1,1)模型及其在疲劳试验数据处理和疲劳试验在线监测中的应用[J].机械强度,1996,18(3):60~63.
[99] Shi B Z. Modeling of the non-equigap GM(1,1) modeling [J]. The Journal of Grey System.1993,5(2):105~114.
[100] Deng J L. A novel GM(1,1) model for non-equigap series [J]. The Journal of Grey System.1997,9(2):111~116.
[101]王钟羡,吴春笃,史学荣.非等间距序列的灰色模型[J].数学的实践与认识,2003,33(10):16~20.
[102]史雪荣,王作雷,张正娣.变参数非等间距GM(1,1)模型及应用[J].数学的实践与认识,2006,36(6):216~220.
[103]蒋卫东,李夕兵,赵国彦.非等时序列预测的时数分离研究[J].系统工程理论与实践,2003,23(1):68~72
[104]戴文战,李俊峰.非等间距GM(1,1)模型建模研究[J].系统工程理论与实践,2005,9(9):89~93.
[105]王叶梅,党耀国,王正新.非等间距GM(1,1)模型背景值的优化[J].中国管理科学,2008,16(4):159~162.
[106] Lin C B, Su S F, Hsu Y T. High-precision forecast using grey models [J]. International Journal of Systems Science, 2001, 32(5):609~619.
[107] Chang L C. Applying the grey prediction model to the global integrated circuit industry. Technol Forecast Soc Change [J]. 2003, 70(6):563~574.
[108] Wang Y F. Predicting stock price using grey prediction system [J]. Expert Systems with Applications, 2002(22): 33~39.
[109] Tsaur R C. Forecasting analysis by using fuzzy grey regression model for solving limited time series data [J]. Soft Comput, 2008, (12):1105~1113.
[110] Tseng F M, Yu H C, Tzeng G H. Applied hybrid Grey Model to forecast seasonal time series [J]. Technol Forecast Soc Change, 2001, 67(2-3):291~302.
[111] Kung L M, Yu S W. Prediction of index futures returns and the analysis of financial spillovers - A comparison between GARCH and the grey theorem [J]. European Journal of Operational Research, 2009, 186(3): 1184 ~1200.
[112]耿立艳.非线性金融波动率模型及其实证研究[D].天津大学,2009.
[113]齐长鑫,汪树玉.灰色系统模型在坝基位移预测中的应用[J].水利学报,1996,(9):49~52.
[114] Li G D, Yamaguchi D, Nagai M. The development of stock exchange simulation predictionmodeling by a hybrid grey dynamic model [J]. Int J Adv Manuf Technol, 2008, (36):195~204.
[115] Chi C T, Shih C L. Efficient path planning of manipulator based on grey prediction. The Journal of Grey System, 2000, 12(1): 19~27.
[116] Hao Y H, Yeh T C J, Gao Z Q, Wang Y R, Zhao Y. A Gray System Model for Studying the Response to Climatic Change: the Liulin Karst Springs, China [J]. Journal of Hydrology, 2006, 328(3-4): 668~676.
[117] Zhang Y G, Xu Y, Wang Z P. GM(1,1) grey prediction of Lorenz chaotic system [J]. Chaos, Solitons and Fractals, 2009, (42): 1003~1009.
[118] Li D C, Yeh C W, Chang C J. An improved grey-based approach for early manufacturing data forecasting [J]. Computers & Industrial Engineering, 2009, (57):1161~1167.
[119] Lu L J, Lewis C, Lin S J. The forecast of motor vehicle, energy demand and CO2 emission from Taiwan’s road transportation sector [J]. Energy Policy, 2009, (37): 2952~2961.
[120]王正新,党耀国,刘思峰,练郑伟. GM(1,1)幂模型求解方法及其解的性质[J].系统工程与电子技术,2009,31(10):2380~2383.
[121]宋彦辉,聂德新.基础沉降预测的灰色Verhulst模型[J].岩土力学,2003,24(1):13~126.
[122]施玉民,陈雯,赵长利.运用灰色Verhulst模型对我国车险业保费收入的预测[J].统计与决策,2008,(12):11~13.
[123]何文章,吴爱弟.估计Verhulst模型中参数的线性规划方法及应用[J].系统工程理论与实践,2006,26(8):141~144.
[124] Dang Y G, Liu S F. The GM Models That x(n) be Taken as Initial Value [J]. Kybernetes 2004,33(2):247~255.
[125] Wang Z X, Dang Y G, Wang Y M. A new grey Verhulst model and its application [C]. Proceedings of the 2007 IEEE International Conference on Grey Systems and Intelligent Services 2007,571~574.
[126] Kayacan E, Ulutas B, Kaynak O. Grey system theory-based models in time series prediction [J]. Expert Systems with Applications, 2010, (37):1784~1789.
[127]王正新,党耀国,刘思峰.无偏灰色Verhulst模型及其应用[J].系统工程理论与实践,2009,29(10):138~144.
[128]许秀莉,罗健. GM(1,1)模型的改进方法及其应用[J].系统工程与电子技术,2002,(4):105~108.
[129]安宁宁,韩兆洲.几种统计平均数大小关系证明的新方法及其推广[J].统计与决策,2006,(7):141~142.
[130]玄光男,程润伟.遗传算法与工程优化[M] .北京:清华大学出版社,2004.
[131]戴华.矩阵论[M] .北京:科学出版社,2001.
[132] [前南斯拉夫]D.S.密特利诺维奇著.张小萍,王龙译.解析不等式[M].北京:科学出版社,1987.
[133]於方,张强,过孝民.中国主要工业废水排放行业的污染特征与行业治理重点[J].环境保护,2003,(14):38~43.
[134]谢红彬,刘兆德,陈雯.工业废水排放的影响因素量化分析[J].长江流域资源与环境, 2004,13(4):394~398.
[135] Hipel K W and McLeod A I. Time Series Modelling of Water Resources and Environmental Systems [M]. Elsevier, Armsterdam, The Netherlands, 1994.