Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl
|
|
- Lora Underwood
- 5 years ago
- Views:
Transcription
1 RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept , HA1N, SE Gothenburg Sweden. Emal: A. Forgo Volvo Car Corporaton Engne Desgn and Development Dept , HA1S SE Gothenburg Sweden Emal: Keywords: Splne nterpolaton method, Observers, Estmaton, Combuston qualty montorng, Spark gnton engnes. Abstract In ths paper new computatonally effcent algorthms for splne nterpolaton method are explaned. Theoretcal comparatve analyss of the splne nterpolaton method wth combned hgh-gan observer and splne nterpolaton method s presented. New splne nterpolaton algorthms are mplemented for estmaton of the engne angular acceleraton from the crankshaft angle measurements. 1 Introducton Numercal calculaton of the dervatves of a sgnal s an old problem n numercal analyss and dgtal sgnal processng. The backward dfference method gves one of the smplest numercal dfferentators. Despte the fact that t s qute common n engneerng applcatons the behavor of the dervatvesvery often accompaned wth peakng phenomena. Splne nterpolaton method proposed n [2] s based on on-lne least-squares polynomal fttng over the movng n tme wndow of a sze w. The advantage of ths method over the backward dfference method s ts good transent behavor. The dea for the splne nterpolaton method s to ft a polynomal of a certan order as a functon of tme n least squares sense and take the dervatves analytcally. Propertes of ths method are descrbed n [2],[1]. However, several practcal problems reman. Relatvely large wndow sze w requres more on-lne computatons and makes practcal mplementaton of the method dffcult. Ths necesstates the development of computatonally effcent recursve algorthms. Moreover, n [2],[1] only constant ntersamplng tme s consdered, whch n many practcal applcatons s not constant. For example, crankshaft angle measurements n automotve engnes are based on the measurements of the crankwheel tooth number and therefore the elapsed tme between two teeth passng a fxed pont vares. Ths necesstates the development of recursve computatonally effcent algorthms for varable dscretzaton step. The contrbutons of ths paper are the followng : new computatonally effcent recursve splne nterpolaton algorthms for varable dscretzaton step ; theoretcal comparatve analyss of the splne nterpolaton method and combned splne method wth hgh gan observer; real-tme mplementaton of the splne nterpolaton method for the crankshaft acceleraton estmaton. 2 Estmaton of the dervatves of sgnal The frst step s to choose the nterpolatng polynomal as ^ff + c 1 t + ::: + c n t n (1) where ^ff s an estmate of the measured sgnal ff, t s contnuous tme, c, 0; :::; n are coeffcents to be found. The estmates of the dervatves are obtaned by dfferentatng (1) analytcally.
2 Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to place the orgn at the pont k (w 1), where k s the step number, then there s a wndow ofdatafff 0 ; :::; ff g measured at f0; :::; t g respectvely. Varable dscretzaton step s presented as t t t 1, 1;:::;w 1. The sum to be mnmzed at every step s S (ff ( + c 1 t + ::: + c n t n ))2 (2) where t, 0;:::;w 1 s a dscrete tme whch corresponds to the sgnal measurements ff, t 0 0. Mnmum of S s acheved when equatng to zero partal dervatves of S wth respect to c, 0 Equatons (3) can be wrtten as follows w + c 1 ::: :::::: t + ::: + c n t + c 1 t 2 + ::: + c n t n t n+1 ff (4) ff t (5) t n + c 1 t n+1 + ::: + c n t 2n ff t n (6) where equaton (4) 0, equaton 0 and fnally equaton n 1 The system (4) - (6) should be resolved at every step n order to fnd the coeffcents c, 0; :::; n. It s clear that for a suffcently large wndow sze w, the calculaton of sums n the system (4) - (6) requres a lot of computatonal power and our next step s to present recursve computatonally effcent algorthms to compute these sums. Let us consder one step of the movng wndow of a sze w. Suppose that at step k 1 there s the followng data fff 0 ; :::; ff g measured at f0;:::;t g, and at step k there s fff 1 ; :::; ff w g measured at f0; :::; t w g. The new value ff w measured at the tme t w enters the wndow ( buffer) and the value ff 0 leaves the wndow. Our problem statement s to fnd computatonally effcent recursve algorthms to compute the sums P S mk va the sums on the prevous step S mk 1 tm, where m 1; :::2n. Defne the sum of order m at step k as follows S mk w 1 P w 1 tm t m t m 2 +( t 2 + t 3 ) m + ::: + ( t 2 + t 3 + ::: + t w ) m (7) where t 1 0. The sum (7) should be computed usng the same sum on a prevous step k 1 whch can wrtten as S mk 1 t m t m 1 +( t 1 + t 2 ) m + ::: + ( t 1 + t 2 + ::: + t ) m (8) and the sums of lower order at step k whch are defned as S (m j)k w 1 t m j t m j 2 +( t 2 + t 3 ) m j + ::: + ( t 2 + t 3 + ::: + t w ) m j (9) where 1» j» (m 1), m>1. Startng wth (8) and usng the followng dentty (x + y) m m j0 where m 0; 1; 2; :::, one gets Cj m xj y m j ; Cj m m! (10) j!(m j)!
3 S mk 1 (w 1) tm 1 + tm 2 +( t 2 + t 3 ) m + ::: +( t 2 + t 3 + ::: + t ) m + Cj m ( tj 1 tm j 2 + ::: + t j 1 ( t 2 + ::: + t ) m j ) (11) Notce that and t m 2 +( t 2 + t 3 ) m + ::: +( t 2 + t 3 + ::: + t ) m S mk ( t 2 + t 3 + ::: + t w ) m (12) C m j ( tj 1 tm j 2 + ::: + t j 1 ( t 2 + ::: + t ) m j ) C m j tj (S (m j)k ( t 2 + ::: + t w ) m j ) (13) Substtutng (12) and (13) n (11) one gets where S ffmk S ffmk 1 (ff 1 + ff 2 + ::: + ff ) t m 1 C m j tj 1 (S ff(m j)k ff w ( t 2 + t 3 + ::: + t w ) m j ) +ff w ( t 2 + t 3 + ::: + t w ) m (15) S ffmk 1 ff t m ff 1 t m 1 + ff 2( t 1 + t 2 ) m + ::: +ff ( t 1 + t 2 + ::: + t ) m (16) s the sum on the step (k 1), and S ff(m j)k w 1 ff t m j +ff 3 ( t 2 + t 3 ) m j + ff 2 t m j 2 ::: + ff w ( t 2 + t 3 + ::: + t w ) m j (17) where 1» j» (m 1), m>1. The order of the nterpolatng polynomal should be so selectedtobeaslow as possble n order to reduce the computatonal burden and to flter out measurement nose. In the next Secton the detaled soluton of the nterpolaton problem for the second order polynomal s presented. Ths example s used n Secton 4 for the crankshaft acceleraton estmaton. S mk S mk 1 (w 1) tm 1 C m j tj (S (m j)k ( t 2 + ::: + t w ) m j ) +( t 2 + t 3 + ::: + t ) m (14) where Cj m m!,1» j» (m 1), m>1. j!(m j)! Our next step s to calculate the sums on the rght hand sde of the equatons (4) - (6). The sums S ffmk P w 1 ff t m should be P calculated va the sums on the prevous step S ffmk 1 ff t m, where m 2; :::n. Usng smlar arguments one can show that 3 Second Order Example 3.1 Splne Interpolaton Method For the second order polynomal (1), where n 2, equatons (4) - (6) can be wrtten as follows w + c 1 t + c 2 t + c 1 t 2 + c 2 t 2 ff (18) t 3 ff t (19)
4 t 2 + c 1 t 3 + c 2 t 4 ff t 2 (20) problem, repeatable ntalzaton of the algorthms s requred. where the sums are computed usng recursve formulas (14) and (17). Presentng equatons (18) - (20) n matrx form yelds where A Ac b (21) a 11 a 12 a 13 1 a 21 a 22 a 23 A, c T ( ;c 1 ;c 2 ), b T a 31 a 32 a 33 (b 1 ;b 2 ;b 3 ), where 3.2 Combnaton of hgh gan observer and splne nterpolaton method and ther comparatve analyss In [2], [1] t s mentoned that the splne nterpolaton method can be combned wth hgh-gan observers. The combned scheme shows mproved transent behavor. In ths Secton the comparatve analyss of the splne nterpolaton method and combned method s presented. The sgnal s estmated va the second order polynomal descrbed above a 11 w; a 12 S 1k ; a 13 S 2k a 21 S 1k ; a 22 S 2k ; a 23 S 3k a 31 S 2k ; a 32 S 3k ; a 33 S 4k b 1 S ffk ; b 2 S ff1k ; b 3 S ff2k Notce that, a 12 a 21 and a 22 a 31 a 13, a 23 a 32. In order to fnd splne coeffcents c, 0; 1; 2 the matrx equaton (21) should be solved wth respect to c at every step ( c A 1 b). To ths end matrx A s nverted analytcally. Remark 1. Ths algorthm has only one parameter to be optmzed, ths beng the sze of the movng wndow w. If the dervatve of the measured sgnal changes slowly t s advsable to have a relatvely large wndow sze to flter out measurement nose. If the dervatve changes quckly, the wndow sze should be suffcently small to capture correspondng fast changes n the dervatve. The dsadvantage of a small wndow sze s the nose n the estmated sgnal. Ideally, the wndow sze should be adjustable so that t s small enough durng transents to capture fast changes n the dervatve of the sgnal, and large enough under steady-state condtons to flter out measurement and space-dscretzaton nose. Remark 2. For constant dscretzaton step the sums S m do not change wth tme, matrx A s constant and computatonal burden s mnmal. Remark 3. The recursve computatons presented above accumulate an approxmaton error, whch ncreases wth tme. To avod ths error accumulaton ^ff(t) + c 1 t + c 2 t 2 (22) where ^ff(t) s the estmate of the sgnal ff. A combnaton of the splne nterpolaton algorthm wth a smple hgh gan observer s presented as follows y(t) 1 f ff + y 1(t) (23) _y 1 (t) 1 f y(t)+^ ff; y 1 (t 0 )^_ff(t 0 ) 1 f ff(t 0)(24) where y(t) s the estmate of the _ff, f > 0 s the algorthm parameter, ^_ff(t 0 ) s the estmate of the dervatve from splne nterpolaton method evaluated at tme t 0, ^ ff s splne estmate of the second dervatve. The condton y 1 (t 0 )^_ff(t 0 ) 1 f ff(t 0) corresponds to the followng ntal value assgnment y(t 0 )^_ff(t 0 ). Ths means that hgh gan observer (23), (24) s ntalzed to the splne estmate. Moreover, the second order dervatve ^ ff s used drectly as an nput to the hgh gan observer. Notce that, ^_ff(t 0 )c 1 +2c 2 t 0 and ^ ff 2c 2. Remark 4. Equatons (23) and (24) represent the followng estmator of the dervatve y(p) p ff(p) (25) fp +1 f ^ ff 0, where p jw s a Laplace varable. Estmates from the splne nterpolaton method are used as a feedforward part n the estmator (25).
5 Assumng a splne approxmaton error, then ths error s expressed n terms of the error n splne coeffcents and the true sgnal can be expressed as follows e 2 (t 1 )( c 1 +2 c 2 t 0 2 c 2 f)e t 1 t 0 f +2 c 2 f (34) ff(t) Λ + c 1Λ t + c 2Λ t 2 (26) where ff(t) s the sgnal, c Λ, 0; 1; 2 are true coeffcents. Assumng a constant error n the coeffcents w.r.t. tme one gets c c Λ + c (27) where c, 0; 1; 2 are constant errors n splne coeffcents. Our task s to compare the followng estmaton errors where t 1 s a fxed tme. Frst we evaluate e 1 (t 1 ) e 1 (t 1 ) ^_ff(t 1 ) _ff(t 1 ) (28) e 2 (t 1 ) y(t 1 ) _ff(t 1 ) (29) e 1 (t 1 ) c 1 +2 c 2 t 1 (30) The next step s to evaluate e 2 (t 1 ). Dfferentatng (23) and takng nto account that ^ ff ff +2 c 2 one gets _y ff 1 f (y _ff)+2 c 2 (31) The soluton of (31) s the followng y(t) _ff(t) (y(t 0 ) _ff(t 0 ) 2 c 2 f)e t t 0 f +2 c 2 f (32) Intalzaton of hgh gan observer gves the followng: y(t 0 )^_ff(t 0 ). Takng nto account that ^_ff(t 0 ) _ff(t 0 ) c 1 +2 c 2 t 0 one obtans y(t) _ff(t) ( c 1 +2 c 2 t 0 2 c 2 f)e t t 0 f Evaluatng e 2 (t 1 ) one gets +2 c 2 f (33) Comparng (34) and (30), t can be seen that je 2 (t 1 )j can always be made smaller than je 1 (t 1 )j by reducng the desgn parameter f. However, n the presence of nose the reducton of f leads to a deteroraton of the sgnal qualty. The advantages of ths combned scheme can be shown f the parameter f can be reduced to a suffcently small value wthout any sgnal qualty deteroraton. As the smulaton results dd not show any sgnfcant superorty of the combned method over the splne nterpolaton method and the latter s chosen for mplementaton. 4 Implementaton Results In ths Secton the mplementaton results are presented for a second order splne fttng method, as descrbed n Secton 3.1. A Volvo S80 passenger car equpped wth a crankshaft angle sensor s used for experments. The method descrbed n Secton 3.1 s used to estmate the crankshaft acceleraton,whch gves vtal nformaton about the qualty of the combuston n the engne, from the crankshaft angle measurements. The car s equpped wth a specal Engne Control Unt, whch s called Volvo Rapd Prototypng System. Ths system conssts of the PM (Power Module) and the AM (Applcaton Module) connected together va a 1 Mbt/s CAN channel. The PM provdes nformaton based on a number of engne sensor sgnals. The crankshaft wheel has 58 teeth spaced every 6 degrees and a gap correspondng to two mssng teeth. Engne speed and cylnder poston can be calculated from the resultng crankshaft sensor sgnal. The PM calculates the engne speed based on the measured tme between the tooth events va the backward dfference method. Ths engne speed s transmtted along wth other measurements every 4th mllsecond to the AM. The 4 mllsecond sampled crankshaft angle calculated from the tooth number sgnal s the nput to the crankshaft acceleraton estmaton, whch runs n the
6 AM. The accuracy of the estmaton s verfed by comparng the engne speed, calculated by the splne nterpolaton method, wth the engne speed calculated n the PM. Implementaton results are presented n Fg.1 whch shows engne speeds as a functon of tme. The accuracy of the estmaton s verfed by comparng fast engne speed measured n the PM module ( dotted lne ) and engne speed calculated va splne nterpolaton method ( sold lne ). Crankshaft Acceleraton [rev/sec 2 ] Tme [sec] Engne Speed [rev/sec] Fgure 2 Expermental results. Two engne cycles of 5 cylnder engne. The msfre s generated on the second cycle n the cylnder 4. Tooth number sgnal s plotted wth dashdot lne. Crankshaft acceleraton estmated by the splne nterpolaton method s plotted as a functon of tme wth a sold lne Tme [sec] Fgure 1 Expermental results. The comparson of two engne speeds as functons of tme. Engne speed computed n PM s plotted wth dotted lne. Engne speed calculated wth splne nterpolaton method s plotted wth sold lne. 5 Concluson Computatonally effcent algorthms for splne nterpolaton method allowng a reducton n computatonal burden as well as the mplementaton of the splne nterpolaton method are presented. A further challenge s to mprove the performance of these algorthms by adjustng the sze of the movng wndow. References Fgure 2 shows an applcaton of the proposed recursve splne nterpolaton method to the msfre dagnostcs. In ths Fgure two engne cycles of a 5 cylnder engne are shown, and the msfre s generated on the second cycle of cylnder N 4. The tooth number sgnal s plotted wth a dashdot lne. Crankshaft acceleraton estmated by the splne nterpolaton method s plotted wth a sold lne. In the event of a msfre, the crankshaft acceleraton changes dramatcally, whch permts the cylnder ndvdual msfre detecton. [1] Dabroom A., Khall H. (1999). Dscrete-Tme Implementaton of Hgh-Gan Observers for Numercal Dfferentaton, Internatonal Journal of Control 72 (17), p.p [2] Dop S., Grzzle J., Moraal P., Stefanopoulou A. (1994). Interpolaton and Numercal Dfferentaton for Observer Desgn, Proc. Amercan Control Conference, Baltmore, Maryland, p.p
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationCurve Fitting with the Least Square Method
WIKI Document Number 5 Interpolaton wth Least Squares Curve Fttng wth the Least Square Method Mattheu Bultelle Department of Bo-Engneerng Imperal College, London Context We wsh to model the postve feedback
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationTime-Varying Systems and Computations Lecture 6
Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy
More informationParameter Estimation for Dynamic System using Unscented Kalman filter
Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,
More informationPolynomial Regression Models
LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance
More informationLab 2e Thermal System Response and Effective Heat Transfer Coefficient
58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationLeast squares cubic splines without B-splines S.K. Lucas
Least squares cubc splnes wthout B-splnes S.K. Lucas School of Mathematcs and Statstcs, Unversty of South Australa, Mawson Lakes SA 595 e-mal: stephen.lucas@unsa.edu.au Submtted to the Gazette of the Australan
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationAnnexes. EC.1. Cycle-base move illustration. EC.2. Problem Instances
ec Annexes Ths Annex frst llustrates a cycle-based move n the dynamc-block generaton tabu search. It then dsplays the characterstcs of the nstance sets, followed by detaled results of the parametercalbraton
More informationAn efficient algorithm for multivariate Maclaurin Newton transformation
Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,
More informationLecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationCubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two
World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationLecture 2: Numerical Methods for Differentiations and Integrations
Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationStatistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA
14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton
More informationDepartment of Electrical & Electronic Engineeing Imperial College London. E4.20 Digital IC Design. Median Filter Project Specification
Desgn Project Specfcaton Medan Flter Department of Electrcal & Electronc Engneeng Imperal College London E4.20 Dgtal IC Desgn Medan Flter Project Specfcaton A medan flter s used to remove nose from a sampled
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationTracking with Kalman Filter
Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
More informationThis column is a continuation of our previous column
Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationNON-LINEAR CONVOLUTION: A NEW APPROACH FOR THE AURALIZATION OF DISTORTING SYSTEMS
NON-LINEAR CONVOLUTION: A NEW APPROAC FOR TE AURALIZATION OF DISTORTING SYSTEMS Angelo Farna, Alberto Belln and Enrco Armellon Industral Engneerng Dept., Unversty of Parma, Va delle Scenze 8/A Parma, 00
More informationMultivariate Ratio Estimator of the Population Total under Stratified Random Sampling
Open Journal of Statstcs, 0,, 300-304 ttp://dx.do.org/0.436/ojs.0.3036 Publsed Onlne July 0 (ttp://www.scrp.org/journal/ojs) Multvarate Rato Estmator of te Populaton Total under Stratfed Random Samplng
More informationTOPICS MULTIPLIERLESS FILTER DESIGN ELEMENTARY SCHOOL ALGORITHM MULTIPLICATION
1 2 MULTIPLIERLESS FILTER DESIGN Realzaton of flters wthout full-fledged multplers Some sldes based on support materal by W. Wolf for hs book Modern VLSI Desgn, 3 rd edton. Partly based on followng papers:
More informationComparison of Wiener Filter solution by SVD with decompositions QR and QLP
Proceedngs of the 6th WSEAS Int Conf on Artfcal Intellgence, Knowledge Engneerng and Data Bases, Corfu Island, Greece, February 6-9, 007 7 Comparson of Wener Flter soluton by SVD wth decompostons QR and
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
More informationNice plotting of proteins II
Nce plottng of protens II Fnal remark regardng effcency: It s possble to wrte the Newton representaton n a way that can be computed effcently, usng smlar bracketng that we made for the frst representaton
More informationDO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.
EE 539 Homeworks Sprng 08 Updated: Tuesday, Aprl 7, 08 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. For full credt, show all work. Some problems requre hand calculatons.
More informationOn an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1
On an Extenson of Stochastc Approxmaton EM Algorthm for Incomplete Data Problems Vahd Tadayon Abstract: The Stochastc Approxmaton EM (SAEM algorthm, a varant stochastc approxmaton of EM, s a versatle tool
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationDr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur
Analyss of Varance and Desgn of Experment-I MODULE VII LECTURE - 3 ANALYSIS OF COVARIANCE Dr Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Any scentfc experment s performed
More informationThe Quadratic Trigonometric Bézier Curve with Single Shape Parameter
J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma
More informationHomework Assignment 3 Due in class, Thursday October 15
Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationModeling curves. Graphs: y = ax+b, y = sin(x) Implicit ax + by + c = 0, x 2 +y 2 =r 2 Parametric:
Modelng curves Types of Curves Graphs: y = ax+b, y = sn(x) Implct ax + by + c = 0, x 2 +y 2 =r 2 Parametrc: x = ax + bxt x = cos t y = ay + byt y = snt Parametrc are the most common mplct are also used,
More informationPop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing
Advanced Scence and Technology Letters, pp.164-168 http://dx.do.org/10.14257/astl.2013 Pop-Clc Nose Detecton Usng Inter-Frame Correlaton for Improved Portable Audtory Sensng Dong Yun Lee, Kwang Myung Jeon,
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationMAE140 - Linear Circuits - Winter 16 Final, March 16, 2016
ME140 - Lnear rcuts - Wnter 16 Fnal, March 16, 2016 Instructons () The exam s open book. You may use your class notes and textbook. You may use a hand calculator wth no communcaton capabltes. () You have
More informationCOEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN
Int. J. Chem. Sc.: (4), 04, 645654 ISSN 097768X www.sadgurupublcatons.com COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN R. GOVINDARASU a, R. PARTHIBAN a and P. K. BHABA b* a Department
More informationCONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION
CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING N. Phanthuna 1,2, F. Cheevasuvt 2 and S. Chtwong 2 1 Department of Electrcal Engneerng, Faculty of Engneerng Rajamangala
More informationReview of Taylor Series. Read Section 1.2
Revew of Taylor Seres Read Secton 1.2 1 Power Seres A power seres about c s an nfnte seres of the form k = 0 k a ( x c) = a + a ( x c) + a ( x c) + a ( x c) k 2 3 0 1 2 3 + In many cases, c = 0, and the
More informationAdiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram
Adabatc Sorpton of Ammona-Water System and Depctng n p-t-x Dagram J. POSPISIL, Z. SKALA Faculty of Mechancal Engneerng Brno Unversty of Technology Techncka 2, Brno 61669 CZECH REPUBLIC Abstract: - Absorpton
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationCOMPUTATIONALLY EFFICIENT WAVELET AFFINE INVARIANT FUNCTIONS FOR SHAPE RECOGNITION. Erdem Bala, Dept. of Electrical and Computer Engineering,
COMPUTATIONALLY EFFICIENT WAVELET AFFINE INVARIANT FUNCTIONS FOR SHAPE RECOGNITION Erdem Bala, Dept. of Electrcal and Computer Engneerng, Unversty of Delaware, 40 Evans Hall, Newar, DE, 976 A. Ens Cetn,
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationStudy on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI
2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,
More informationDEGREE REDUCTION OF BÉZIER CURVES USING CONSTRAINED CHEBYSHEV POLYNOMIALS OF THE SECOND KIND
ANZIAM J. 45(003), 195 05 DEGREE REDUCTION OF BÉZIER CURVES USING CONSTRAINED CHEBYSHEV POLYNOMIALS OF THE SECOND KIND YOUNG JOON AHN 1 (Receved 3 August, 001; revsed 7 June, 00) Abstract In ths paper
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationOn the correction of the h-index for career length
1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat
More informationA PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS
HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,
More informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More informationComposite Hypotheses testing
Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationFor now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.
Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson
More informationVQ widely used in coding speech, image, and video
at Scalar quantzers are specal cases of vector quantzers (VQ): they are constraned to look at one sample at a tme (memoryless) VQ does not have such constrant better RD perfomance expected Source codng
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More information( ) [ ( k) ( k) ( x) ( ) ( ) ( ) [ ] ξ [ ] [ ] [ ] ( )( ) i ( ) ( )( ) 2! ( ) = ( ) 3 Interpolation. Polynomial Approximation.
3 Interpolaton {( y } Gven:,,,,,, [ ] Fnd: y for some Mn, Ma Polynomal Appromaton Theorem (Weerstrass Appromaton Theorem --- estence ε [ ab] f( P( , then there ests a polynomal
More informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationOn balancing multiple video streams with distributed QoS control in mobile communications
On balancng multple vdeo streams wth dstrbuted QoS control n moble communcatons Arjen van der Schaaf, José Angel Lso Arellano, and R. (Inald) L. Lagendjk TU Delft, Mekelweg 4, 68 CD Delft, The Netherlands
More informationAssortment Optimization under MNL
Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.
More informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationBeyond Zudilin s Conjectured q-analog of Schmidt s problem
Beyond Zudln s Conectured q-analog of Schmdt s problem Thotsaporn Ae Thanatpanonda thotsaporn@gmalcom Mathematcs Subect Classfcaton: 11B65 33B99 Abstract Usng the methodology of (rgorous expermental mathematcs
More information