Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2014, Article ID 634852, 7 pages http://dx.doi.org/10.1155/2014/634852
Research Article A Novel Parameter Estimation Method for Muskingum Model Using New Newton-Type Trust Region Algorithm
Zhou Sheng,1 Aijia Ouyang,2,3 Li-Bin Liu,4 and Gonglin Yuan1
1 College of Mathematics and Information Science, Guangxi University, Nanning, Guangxi 530004, China 2Department of Electrical and Information Engineering, Hunan Institute of Traffic Engineering, Hengyang, Hunan 421001, China 3School of Information Science and Engineering, Hunan City University, Yiyang, Hunan 413000, China 4Department of Mathematics and Computer Science, Chizhou College, Chizhou, Anhui 247000, China
Correspondence should be addressed to Li-Bin Liu; [email protected]
Received 28 August 2014; Accepted 4 December 2014; Published 21 December 2014
Academic Editor: Valder Steffen Jr.
Copyright Β© 2014 Zhou Sheng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Parameters estimation of Muskingum model is very significative in both exploitation and utilization of water resources and hydrological forecasting. The optimal results of parameters directly affect the accuracy of flood forecasting. This paper considers the parameters estimation problem of Muskingum model from the following two aspects. Firstly, based on the general trapezoid formulas, a class of new discretization methods including a parameter π to approximate Muskingum model is presented. The accuracy of these methods is second-order, when π =1/3ΜΈ .Particularly,ifwechooseπ=1/3, the accuracy of the presented method can be improved to third-order. Secondly, according to the Newton-type trust region algorithm, a new Newton-type trust region algorithm is given to obtain the parameters of Muskingum model. This method can avoid high dependence on the initial parameters. The average absolute errors (AAE) and the average relative errors (ARE) of the proposed algorithm of parameters estimation for Muskingum model are 8.208122 and 2.462438%, respectively, where π=1/3. It is shown from these results that the presented algorithm has higher forecasting accuracy and wider practicability than other methods.
1. Introduction where π(π‘) represents the channel storage at time π‘; πΌ(π‘) and π(π‘) represent the rates of inflow and outflow at time π‘, Flood routing in open channels is a very important tool in respectively. The linear Musikingum model is the design of flood protection measures to estimate how the proposedmeasureswillaffectthebehavioroffloodwaves π (π‘) =πΎ[π₯πΌ (π‘) + (1βπ₯) π (π‘)] , (2) in rivers so that enough protection and economic solutions can be found [1]. In general, flood routing procedures can where πΎ is a storage time constant for the river reach and π₯ be classified as either hydrologic or hydraulic. Hydrologic is a weighting factor commonly varying between 0.0 and 0.3 routing method is on the basis of the storage-continuity fortheriverchannel. equation, whereas hydraulic routing method is on the basis of It is worth mentioning that the parameters πΎ and π₯ in (2) both continuity and momentum equations. One of the hydro- are graphically estimated by a trial-and-error procedure [2]. If logic routing approaches, the Muskingum method, was first π₯ is obtained, the values of [π₯πΌ(π‘)+(1βπ₯)π(π‘)] are calculated developed by McCarthy for flood control studies in the by using observed data in both upstream and downstream Muskingum river basin in Ohio. and plotted against π. The particular value which generates The following continuity and storage equations are used the loop is accepted as the best estimation of π₯.Theslopeof to describe the Muskingum model: the straight line fits through the loop derives πΎ. Although the trial-and-error procedure has been used for many years, it is dπ (π‘) =πΌ(π‘) βπ(π‘) , (1) time consuming and subjective interpretation. Therefore, in dπ‘ order to avoid subjective interpretations of observed data in 2 Mathematical Problems in Engineering estimating πΎ and π₯, the following finite difference scheme [3] are that the global search ability is strong and the program is used to the resulting ordinary differential equation (1): is relatively simple to design. However, these methods have unavoidable disadvantages of slow convergence precision, π(π‘π)βπ(π‘πβ1) low precision solution in limited generations, and premature Ξπ‘ convergence. (3) πΌ(π‘)+π(π‘ ) π(π‘)+π(π‘ ) In a word, almost all of the above methods used to =[ π πβ1 ]β[ π πβ1 ], estimate parameters πΎ and π₯ are based on (3).Weknowthat 2 2 2 the truncation error of (3) is only π((Ξπ‘) ).Thus,itisvery necessary to design a higher accuracy method to discrete where differential equation (1). Based on the generalized trapezoidal
π(π‘π) = πΎ [π₯πΌπ (π‘ )+(1βπ₯) π(π‘π)] , formula [16], this paper will first develop a class of new differ- ence scheme which contains a parameter π to approximate the π(π‘πβ1) = πΎ [π₯πΌπβ1 (π‘ )+(1βπ₯) π(π‘πβ1)], π=2,...,π, differential equation (1).Then,similartotheaboveoptimiza- (4) tion problem (7), we can obtain a class of new unconstrained nonlinear optimization problems which contain parameter and π(π‘π) represents observed outflow discharges at time π‘π, π. It is noted that the accuracy of the presented difference πΌ(π‘π) represents the inflow discharge at time interval π‘π, Ξπ‘ is schemes to approximate the differential equation (1) can be thetimestep,andπ is the total time number. Substituting (4) improved to third order, when π=1/3. In other words, we into (3),wehave can get a higher accuracy parameter estimation model, if π= 1/3. In addition, based on the Newton-type trust region algo- π (π‘π+1) =π0πΌ (π‘π+1) +π1πΌ (π‘π) +π2π (π‘π) , (5) rithm [17], this paper will design a new Newton-type trust region algorithm (NN-TTRA) for estimating the Muskingum where model parameters. Briefly speaking, this algorithm can avoid βπΎπ₯ + (Ξπ‘/2) high dependence on the initial parameters and find the global π = , 0 πΎ (1βπ₯) + (Ξπ‘/2) optimization solution quickly. πΎπ₯ + (Ξπ‘/2) π = , 1 πΎ (1βπ₯) + (Ξπ‘/2) (6) 2. A New Parameter Estimation of Muskingum Model πΎ (1βπ₯) β (Ξπ‘/2) π = . 2 πΎ (1βπ₯) + (Ξπ‘/2) 2.1. Model Description. At first, (1) can be rewritten as follows:
By minimizing the sum of the square of the deviations dπ (π‘) between observed and calculated outflows, we obtain the =π(π‘) +π(π‘) , (8) dπ‘ following objective function:
πβ1 π(π‘) = πΌ(π‘) β π(π‘) βπ(π‘) 2 where . min π=β [π(π‘π+1)βπ0πΌ(π‘π+1)βπ1πΌ(π‘π)βπ2π(π‘π)] . (7) Then, applying the generalized trapezoidal formula [16] π=1 GTF (π)ofChawlaetal.to(8),wecanobtain In the past two decades, in order to obtain the parameters π(π‘ ) πΎ and π₯, some researchers adopted many optimization Μ d π+1 π(π‘π)=π(π‘π+1)βΞπ‘ , (9) methods to solve the above optimization problem (7).Until dπ‘ now, these methods can be classified into traditional opti- Ξπ‘ π(π‘) mization methods and intelligent algorithms. The existing d π π(π‘π+1)=π(π‘π)+ [(1βπ) traditional methods include nonlinear programming method 2 dπ‘ (NPM) [4], the least-square method (L-SM) [5, 6], method (10) π(π‘Μ ) π(π‘ ) of trial-and-error (TAE) [2], the minimum area method +πd π + d π+1 ]. (MAM) [1, 7], the Broydene-Fletchere-Goldfarbe-Shanno dπ‘ dπ‘ (BFGS) technique [8], test-method and least residual square method [9], and Nelder-Mead simplex method [10]. These Combining (8), (9),and(10),weget traditional methods have their own advantages, but there also exist complex calculations and poor generality disadvantages, Μ and some of these methods are related to the selection π(π‘π)=(1βΞπ‘) π(π‘π+1)βΞπ‘π(π‘π+1), (11) of initial point, which is easy to fall into local optimum. Ξπ‘ The existing intelligent algorithms include genetic algorithm π(π‘π+1)=π(π‘π)+ [π π+1(π‘ )+π(π‘π)+(1βπ) π(π‘π) [11], harmony search [12], Gray-encoded accelerating genetic 2 algorithm [9], particle swarm optimization [13], immune +ππ(π‘Μ )+π(π‘ )] . clonal selection algorithm [14], and differential evolution π π+1 algorithm [15]. The advantages of those intelligent algorithms (12) Mathematical Problems in Engineering 3
Substituting (11) into (12),wefinallyobtain 3. Algorithm Construction Ξπ‘ 3.1. Newton-Type Trust Region Algorithm. As far as we know, (1 β π) ππ+1 (π‘ ) 2 the basic trust region algorithm used to solve the following Ξπ‘ Ξπ‘ unconstrained optimization problem =(1β π) π (π‘ )+ [πΌ (π‘ )βπ(π‘)] 2 π 2 π π (13) minπ (π₯) Ξπ‘ π₯βπ π (20) + (1βΞπ‘π) [πΌ (π‘ )βπ(π‘ )] . 2 π+1 π+1 π=0 was first presented clearly in [19]. Recently, many researchers Obviously, when , (13) becomes (3). have proposed some improved trust region algorithm to solve At last, substituting (4) into (13),wehave (20); see for example Esmaeili and Kimiaei [20]andAmini π(π‘π+1)=π0πΌ(π‘π+1)+π1πΌ(π‘π)+π2π(π‘π), (14) and Ahookhosh [21]. Here, we assume that π₯π is the πth iterative point, ππ = where π(π₯π), ππ =βπ(π₯π),andπ΅π is the πth iteration of Hesse matrix βπΎπ₯ (1β(Ξπ‘/2) π) + (Ξπ‘/2)(1βΞπ‘π) β2π(π₯ ) π π = , π ; the trust region subproblem, which is at the th 0 πΎ (1βπ₯)(1β(Ξπ‘/2) π) + (Ξπ‘/2)(1βΞπ‘π) iterative step of problem (20), can be formulated as follows: πΎπ₯ (1β(Ξπ‘/2) π) + (Ξπ‘/2) π 1 π π1 = , (15) π (π) =π π+ π π΅ π, πΎ (1βπ₯)(1β(Ξπ‘/2) π) + (Ξπ‘/2)(1βΞπ‘π) min π π 2 π (21) πΎ (1βπ₯)(1β(Ξπ‘/2) π) β (Ξπ‘/2) π = . s.t. βπβ β€Ξπ, 2 πΎ (1βπ₯)(1β(Ξπ‘/2) π) + (Ξπ‘/2)(1βΞπ‘π) π By minimizing the residual sum of squares between where Ξ π is a trust region radius and πβπ is an iterative step observed and calculated outflows, the objective function can and ββ βdenotes the Euclidian norm of vectors or its induced be given as follows: matrix norm. Denote πβ1 2 min π=β [π(π‘π+1)βπ0πΌ(π‘π+1)βπ1πΌ(π‘π)βπ2π(π‘π)] . (16) π=1 Ξππ =π(π₯π)βπ(π₯π +ππ), (22)
Ξπ =π (0) βπ (π ), 2.2. Truncation Errors Analysis. Here, we analyze the local π π π π (23) truncation error of (13) at time direction [16]. Firstly, from Ξπ (9),wehave π = π ; π Ξπ (24) Μ 2 π dπ(π‘π) dπ(π‘π+1) d π(π‘π+1) = βΞπ‘ . (17) π‘ π‘ π‘2 d d d the general trust region algorithm [17]forsolvingtheuncon- Then, by using Taylor expansion, (10) can be written as strained optimization problem (20) is given as follows. follows: Algorithm 1. (Ξπ‘)3 π3π(π‘) (Ξπ‘)4 π4π(π‘) β ( π )β ( π )+π((Ξπ‘)5) 4 4 π 12 ππ‘ 24 ππ‘ Step 0.Givenastartingpointπ₯0 βπ , Ξ 0 >0is the initial trust region radium, π>0, π΅0 = Hesse(π₯0).Setπ:=0. Μ (Ξπ‘) π dπ(π‘π) dπ(π‘π) 5 = [ β ]+π((Ξπ‘) ). π βπ ββ€π 2 dπ‘ dπ‘ Step 1.Calculate π,if π , and stop iteration; otherwise, (18) go to Step 2. Μ Substituting for dπ(π‘π)/dπ‘ from (17) into (18), the truncation Step 2. Utilize the smoothing Newton method to obtain the π errors of generalized trapezoidal formula can be obtained as solution π of the subproblem (21). follows: Step 3.Calculatetheππ of formula (24). π 1 π3π(π‘) πΈ=( β ) (Ξπ‘)3 ( π ) 4 12 ππ‘3 Step 4. Regulate trust region radius. If ππ <0.25,letΞ π+1 := 0.5Ξ π > 0.75 βπ β=Ξ Ξ := 2Ξ (19) π;if π and π π,let π+1 π;otherwise, 4 π 1 π π(π‘) let Ξ π+1 := Ξ π. +( β ) (Ξπ‘)4 ( π )+π((Ξπ‘)5). 6 24 ππ‘4 Step 5.Ifππ > 0.25,letπ₯π+1 := π₯π +ππ,updateπ΅π to be From (19), it is clear that the order of our presented π΅π+1 = Hesse (π₯π+1),letπ:=π+1,andgotoStep1;otherwise, 2 π₯ =π₯ π:=π+1 Ξ = method in (13) is π((Ξπ‘) ) if π =1/3ΜΈ .Inparticular,for let π+1 π, , and go to Step 2, where 0 β6 π=1/3, our presented method in (13) is third-order accuracy. (1/10)βπ(π₯0)β, π=10 . 4 Mathematical Problems in Engineering
Obviously, the above trust region subproblem (21) is Start an unconstrained optimization problem whose objective function is a quadratic. For the sake of convenience, we first denote Initialize parameters π§=(π, π,) π , (25) π Termination criterion is met? π (π§) := ( 2 ), 2 2 β 2 2 2 πββπβ2 +Ξπ β (π + βπβ2 βΞπ) +4π Yes
(π΅π +ππΈ)πβππ Use Algorithm 2 to solve (26) subproblem (19) π½ (π§) =πΎβπ (π§)β min {1, βπ (π§)β} , (27) respectively. Then, the following smoothing Newton algo- Calculate formula (22) rithm [17] is given to solve the subproblem (21).
Algorithm 2. No π Μ Formula (22) criterion Step 0.Givenπ0 βπ , π§0 =(π0,π0,π0), π§0 =(π0,0,0).Select is met? parameters πΏ, π, πΎ β (0,1), π0πΎ<1,andπΎβπ(π§0)β < 1;set β:=0. Yes No Step 1.Calculateβπ(π§β)β,ifβπ(π§β)β = 0,andterminate algorithm; otherwise, calculate π½β =π½(π§β). Use (23) to act as the adjustment formula σΈ Step 2.Obtainthesolutionforequationsπ(π§β)+π (π§β)Ξπ§β = π½βπ§Μ,andthesolutionisΞπ§β =(Ξπβ,Ξπβ,Ξπβ). End Step 3.Supposethatπβ is the minimum nonnegative integer π π βπ(π§ +πΏ β Ξπ§ )β β€ [1 β π(1 β π½π )πΏ β ]βπ(π§ )β meeting β β 0 β .Let Figure 1: Flow chart of NN-TTRA. πβ πΌβ := πΏ and π§β+1 =π§β +πΌβΞπ§β.
Step 4.Letβ:=β+1,andgotoStep1. It is proved that the method is of super-linear convergence [22]. The flow chart of NN-TTRA is given in Figure 1. 3.2. A New Newton-Type Trust Region Algorithm. Optimizing abilityofNewton-typetrustregionalgorithm(N-TTRA) highly depends on initial parameters. And the algorithm 4. Numerical Experiments and is short of global optimizing capability although its local Results Analysis optimizing is fast. In order to enhance the global optimization capability and get rid of dependence on initial parameters, In this paper, we provide actual observed data of flood we construct a new Newton-type trust region algorithm by runoff process between Chenggouwan and Linqing segment combining a new BFGS updating formula with N-TTRA. in Nanyunhe River of Haihe River Basin. (Length of the reach is 83.8 km, where there is no tributary, but a levee control on Thebasicideaistoupdateπ΅π+1 with new BFGS formulas at both sides. There may occur lifting irrigation during the water Step 5 in Algorithm 1. This idea guarantees positive definite- delivery, and flood water may discharge into the reach when ness of π΅π+1 to improve the global optimization capability; rainfall is high. But these situations have little effect on flood, meanwhile, it gets rid of dependence on the initial parameter where the routing time interval Ξπ‘ = 12 h.) The detailed selection. New adjustment formula of BFGS is as follows: datacanbeseenin[23]. Here, we will give the numerical β βπ π experiments from the following three aspects. π¦π π¦π π΅ππ ππ π π΅π π΅ =π΅ + β , (28) Firstly, in order to verify the advantage of our new π+1 π π ππ¦β π ππ΅ π π π π π π parameter estimation models (16) by using the generalized trapezoid formula to approximate (1), we use the NN-TTRA π¦β =(π¦ππ /|π¦ππ |)π¦ π =π₯ βπ₯ π¦ =π βπ where π π π π π π, π π+1 π,and π π+1 π. tosolvetheaboveunconstrainedoptimizationproblem(16) β Here, π¦π guarantees positive definiteness of π΅π+1,sothenew by choosing different π. Furthermore, we get the correspond- Newton-type trust region algorithm (NN-TTRA) uses (28) as ing parameters πΎ and π₯ for each π. From these parameters, Μ the new adjustment formula of π΅π+1 at Step 5 in Algorithm 1. we use (14) to obtain the calculated outflow π(π‘π),where Mathematical Problems in Engineering 5
Table 1: Numerical results with different π for flood routing in 1960.
π AAE ARE (%) π AAE ARE (%) π=0 8.208160 2.462451 π=1/3 8.208122 2.462438 π = 1/10 8.208159 2.462451 π=1/2 8.208163 2.462452 π=1/7 8.208154 2.462449 π=2/3 8.208162 2.462452 π=1/5 8.208169 2.462454 π=1 8.208162 2.462451
Table 2: Numerical results with different π for flood routing in 1961.
π AAE ARE (%) π AAE ARE (%) π=0 4.011335 0.998421 π=1/3 4.011335 0.998421 π = 1/10 4.011335 0.998421 π=1/2 4.011336 0.998421 π=1/7 4.011336 0.998421 π=2/3 4.011336 0.998421 π=1/5 4.011336 0.998421 π=1 4.011335 0.998421
Table 3: Numerical results with different π for flood routing in 1964.
π AAE ARE (%) π AAE ARE (%) π=0 10.949204 2.662713 π=1/3 10.949203 2.662712 π = 1/10 10.949205 2.662713 π=1/2 10.949204 2.662713 π=1/7 10.949208 2.662713 π=2/3 10.949203 2.662712 π=1/5 10.949205 2.662713 π=1 10.949204 2.662712
Μ π(π‘1)=π(π‘1). At last, the average absolute errors (AAE) and Table 4: Comparison of performance of the two algorithms. the average relative errors (ARE) can be given as follows: Algorithm Initial value Iterations πΎπ₯ π 1 σ΅¨ σ΅¨ (1, 1) 5 0 β1.938700 = β σ΅¨π(π‘Μ )βπ(π‘)σ΅¨ , AAE σ΅¨ π π σ΅¨ N-TTRA (2, 2) β β β π π=1 (29) (11, 1) 3 11.190740 0.998456 π σ΅¨ Μ σ΅¨ 1 σ΅¨π(π‘π)βπ(π‘π)σ΅¨ (1, 1) 27 11.190740 0.998456 ARE = β σ΅¨ σ΅¨ . π σ΅¨ π(π‘) σ΅¨ NN-TTRA (2, 2) 25 11.190740 0.998456 π=1 σ΅¨ π σ΅¨ (11, 1) 4 11.190740 0.998456 From (29),Tables1, 2,and3 list the AAE and ARE for different π for flood routing in 1960, 1961, and 1964, respectively. It is shown from Tables 1β3 thattheAAEandAREarethesmallest Meanwhile Table 5 also gives the numerical results obtained when π=1/3. In other words, the numerical results given by using method of trial-and-error (TAE) [2], the least-square in Tables 1β3 confirm the theoretical analysis presented in method (L-SM) [6], and direct optimal method (DOM) Section 2.2. [18], respectively. Here, the observed data is the outflow of Secondly, to illustrate the efficiency of the NN-TTRA flood runoff between Chenggouwan and Linqing segment presented in this paper, for different parameter initial values, in Nanyunhe River of Haihe Basin in 1960. It can be seen we employ the N-TTRA and the NN-TTRA to solve the above from Table 5 that the AAE and ARE of NN-TTRA are all less unconstrained optimization problem (16), respectively, where than those of other estimation methods. Thus, we conclude the parameter π is set to 1/3 and a maximum number of that it is very effective by using NN-TTRA to estimate the iterations of the N-TTRA and the NN-TTRA are set to 150. parameters of Muskingum model. In addition, the calculated The numerical results are shown in Table 4.FromTable 4, valuesandobservedvaluesoftheflowsin1960,1961,and1964 we can indicate that the NN-TTRA not only gets rid of are plotted in Figures 2, 3,and4, respectively. From Figures dependence on the initial values of parameters, but also 2β4, it can be seen that the calculated flow data via the NN- guarantees that the number of iterations is finite. Meanwhile, TTRA highly coincides with the observed flow data. In brief, through this algorithm the global optimum value is obtained. the accuracy of the method is satisfactory. However, the N-TTRA highly depends on initial values of parameters. When the selected initial value is far from the optimal ones, the algorithm hardly finds the global optimum 5. Conclusions and even cannot complete iterations in finite times. Finally, Table 5 lists the numerical results calculated by By combining a new BFGS adjustment formula with N- using NN-TTRA to solve problem (16),whereπ=1/3. TTRA, this paper implements parameter estimation of 6 Mathematical Problems in Engineering
Table 5: Comparison of results via several parameter estimation 600 methods. 550 Algorithm πΎπ₯AAE ARE (%) 500 TAE [2] 12.400000 0.100000 10.230000 3.150000 β 450
L-SM [6] 11.790000 0.325000 9.800000 2.960000 s) / β 3 DOM [18] 12.440000 0.260000 9.700000 2.920000 400 NN-TTRA 11.190740 0.998456 8.208122 2.462438 350
500 (m Outο¬ow 300
250 450 200 400 150 s)
/ 051015 20 25 30 3 350 Time (h)
300 Calculated
Outο¬ow (m Outο¬ow Observed 250 Figure 3: Comparison between calculated values and observed 200 values in 1961.
150 0 5 10 15 20 25 30 650 Time (h) 600
Calculated 550 Observed 500 s)
Figure 2: Comparison between calculated values and observed / 3 values in 1960. 450 400
Muskingum model by selecting different initial values for (m Outο¬ow 350 parameters and comparing it with N-TTRA. The results show 300 that the NN-TTRA can get rid of the dependence on initial value selection for parameters when searching solutions for 250 Muskingum model parameters, which avoids the influence 200 of initial value selection for parameters on the optimization 0 5 10 15 20 25 30 35 results and averts local optimum. What is more, this algo- Time (h) rithm is of application value in flood disaster management and should be generalized. In addition, this algorithm can be Calculated extended to other similar parameter estimation problems to Observed help obtain excellent results. Figure 4: Comparison between calculated values and observed values in 1964. Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper. Province (Grant no. 2014GK3043), and the Guangxi Natural Science Foundation (2012GXNSFAA053002). The authors wouldliketoacknowledgeGuangxiCollegesandUniversities Acknowledgments Key Laboratory of Mathematics and Its Applications. ThisworkissupportedbytheNationalNaturalScienceFoun- dation of China (nos. 11301044, 11261006, and 11161003), the References Key Projects of Excellent Young Talents Fund in universities of Anhui Province (2013SQRL095ZD), the Project Supported [1] E. M. Wilson, Engineering Hydrology, Macmillan Education, by Scientific Research Fund of Hunan Provincial Education 1990. Department (Grant no. 13C333), the Project Supported by [2] S. E. Serrano, βThe Theis solution in heterogeneous aquifers,β the Science and Technology Research Foundation of Hunan Ground Water,vol.35,no.3,pp.463β467,1997. Mathematical Problems in Engineering 7
[3] L. Kang, C. Wang, and T. B. Jiang, βA new genetic simulated [22] Y. Yuan and W. Sun, Theory and Methods of Optimization,Sci- annealing algorithm for flood routing model,β Journal of Hydro- ence Press, Beijing, China, 1999. dynamics B,vol.16,no.2,pp.233β239,2004. [23]A.Ouyang,Z.Tang,K.Li,A.Sallam,andE.Sha,βEstimating [4]J.JinandJ.Ding,Genetic Algorithm and Its Applications for parameters of Muskingum model using an adaptive hybrid PSO Water Science, Sichuan University Press, 2000. algorithm,β International Journal of Pattern Recognition and [5] M. A. Gill, βFlood routing by the Muskingum method,β Journal Artificial Intelligence,vol.28,no.1,ArticleID1459003,29pages, of Hydrology,vol.36,no.3-4,pp.353β363,1978. 2014. [6] A. A. Aldama, βLeast-squares parameter estimation for Musk- ingum flood routing,β Journal of Hydraulic Engineering,vol.116, no. 4, pp. 580β586, 1990. [7]V.T.Chow,Handbook of Applied Hydrology, Macmillan Educa- tion, 1990. [8] Z. W. Geem, βParameter estimation for the nonlinear Musk- ingum model using the BFGS technique,β Journal of Irrigation and Drainage Engineering,vol.132,no.5,pp.474β478,2006. [9] J. Chen and X. Yang, βOptimal parameter estimation for Muskingum model based on Gray-encoded accelerating genetic algorithm,β Communications in Nonlinear Science and Numeri- cal Simulation,vol.12,no.5,pp.849β858,2007. [10] R. Barati, βParameter estimation of nonlinear Muskingum models using nelder-mead simplex algorithm,β Journal of Hydrologic Engineering, vol. 16, no. 11, pp. 946β954, 2011. [11] S. Mohan, βParameter estimation of nonlinear Muskingum models using genetic algorithm,β Journal of Hydraulic Engineer- ing,vol.123,no.2,pp.137β142,1997. [12] J. H. Kim, Z. W. Geem, and E. S. Kim, βParameter estimation of the nonlinear Muskingum model using Harmony Search,β Journal of the American Water Resources Association,vol.37,no. 5, pp. 1131β1138, 2001. [13] H.-J. Chu and L.-C. Chang, βApplying particle swarm opti- mization to parameter estimation of the nonlinear muskingum model,β Journal of Hydrologic Engineering,vol.14,no.9,pp. 1024β1027, 2009. [14] J. Luo and J. Xie, βParameter estimation for nonlinear Musk- ingum model based on immune clonal selection algorithm,β Journal of Hydrologic Engineering,vol.15,no.10,pp.844β851, 2010. [15] D.-M. Xu, L. Qiu, and S.-Y. Chen, βEstimation of nonlinear Muskingum model parameter using differential evolution,β Journal of Hydrologic Engineering,vol.17,no.2,pp.348β353, 2012. [16] M. M. Chawla, M. A. Al-Zanaidi, and D. J. Evans, βGeneralized trapezoidal formulas for parabolic equations,β International JournalofComputerMathematics,vol.70,no.3,pp.429β443, 1999. [17] C. Ma, Optimization Method and the Matlab Programming,Sci- ence Press, Beijing, China, 2010. [18] D. Stephenson, βDirect optimization of Muskingum routing coefficient,β Journal of Hydrologic Engineering,vol.36,pp.353β 363, 1979. [19] M. J. D. Powell, βA new algorithm for unconstrained optimiza- tion,β Report T.P. 393, Atomic Energy Research Establishment, Oxfordshire, UK, 1970. [20] H. Esmaeili and M. Kimiaei, βA new adaptive trust-region method for system of nonlinear equations,β Applied Mathemat- ical Modelling,vol.38,no.11-12,pp.3003β3015,2014. [21] K. Amini and M. Ahookhosh, βA hybrid of adjustable trust- region and nonmonotone algorithms for unconstrained opti- mization,β Applied Mathematical Modelling,vol.38,no.9-10,pp. 2601β2612, 2014. Advances in Advances in Journal of Journal of Operations Research Decision Sciences Applied Mathematics Algebra Probability and Statistics Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
The Scientific International Journal of World Journal Differential Equations Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
Submit your manuscripts at http://www.hindawi.com
International Journal of Advances in Combinatorics Mathematical Physics Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
Journal of Journal of Mathematical Problems Abstract and Discrete Dynamics in Complex Analysis Mathematics in Engineering Applied Analysis Nature and Society Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014
International Journal of Journal of Mathematics and Mathematical Discrete Mathematics Sciences
Journal of International Journal of Journal of Function Spaces Stochastic Analysis Optimization Hindawi Publishing Corporation Hindawi Publishing Corporation Volume 2014 Hindawi Publishing Corporation Hindawi Publishing Corporation Hindawi Publishing Corporation http://www.hindawi.com Volume 2014 http://www.hindawi.com http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014 http://www.hindawi.com Volume 2014