11.1 Polynomial Least-Squares Curve Fit
|
|
- Brooke Holt
- 6 years ago
- Views:
Transcription
1 11.1 Polynomial Least-Squares Curve Fit A. Purpose This subroutine determines a univariate polynomial that fits a given disrete set of data in the sense of minimizing the weighted sum of squares of residuals. The fitted polynomial an be onstrained to math some data points exatly by appropriate setting of the a priori standard deviations of the data errors. Auxiliary subroutines, desribed in Chapter 11.2, may be used to evaluate, differentiate, or integrate a polynomial produed by this fitting proedure. B. Usage B.1 Program Prototype, Single Preision INTEGER M, NMAX, NDEG REAL X( M),Y( M),SD( M or =1),P( NMAX+3), SIGFAC, W( (NMAX+3) (NMAX+3)) LOGICAL SEEKN, COMTRN, CHBBAS Assign values to M, X(), Y(), SD(), NMAX, SEEKN, COMTRN and CHBBAS. If COMTRN =.FALSE. values must also be assigned to P(1) and P(2). CALL SPFIT (M, X, Y, SD, NMAX, SEEKN, COMTRN, CHBBAS, P, NDEG,SIGFAC, W) Following the CALL to SPFIT one may wish to use SCPVAL or SMPVAL to evaluate the fitted polynomial at speifi points, SCPDRV or SMPDRV to differentiate the polynomial, or SCPINT or SMPINT to integrate the polynomial. See Chapter 11.2 for desriptions of these subprograms. B.2 Argument Definitions M [in] Number of data points. X() [in] Array of values of the independent variable. Values need not be ordered and may be repeated. Y() [in] Array of values of the dependent variable indexed to orrespond to the X() values. SD() [in] If SD(1) > 0., eah SD(i) must be positive and will be assumed to be the user s a priori estimate of the standard deviation of the unertainty (observational errors, et.) in the orresponding data value Y(i). If SD(1) < 0., SD(1) will be used as the a priori standard deviation of eah data value Y(i). In this ase the array SD() an be dimensioned SD(1) sine loations following SD(1) will not be used Calif. Inst. of Tehnology, 2015 Math à la Carte, In. NMAX [in] Speifies the maximum degree polynomial to be onsidered. Require NMAX 0. If NMAX > M 1 the subroutine will funtion as though NMAX = M 1. SEEKN [in] If SEEKN =.TRUE. the subroutine will determine the optimum value of NDEG NMAX in the sense desribed below in Setion D. If SEEKN =.FALSE. the subroutine will set NDEG = NMAX unless this is a singular or nearly singular problem, in whih ase NDEG will be redued as neessary. COMTRN [in] If COMTRN =.TRUE. the subroutine will determine the transformation parameters P(1) and P(2) as desribed in Setion D. If COMTRN =.FALSE. the subroutine will use the values of P(1) and P(2) that have been set by the user. CHBBAS [in] If CHBBAS =.TRUE. the subroutine will use the Chebyshev basis. If CHBBAS =.FALSE. the subroutine will use the monomial basis. The Chebyshev basis is reommended for better numerial stability, when NMAX is greater than about six. P() [inout] Parameter vetor defining a polynomial aording to Eqs. (1) and (2) or (3) and (4). Values of P(1) and P(2) may either be set by the user, or by the subroutine, depending upon the value of COMTRN. Values of P(j+3), j = 0,..., NMAX will be omputed by the subroutine. If NDEG < NMAX, the subroutine will set P(i+3) = 0, for i = NDEG + 1,..., NMAX. NDEG [out] Set by the subroutine to indiate the degree of the fitted polynomial. On return NDEG will satisfy 0 NDEG NMAX if the omputation was suessful and will be set to 1 if an error ourred. SIGFAC [out] Number omputed by the subroutine. See disussion of SIGFAC in Setions C and D following. W() [srath] Array of working storage. B.3 Modifiations for Double Preision For double preision usage hange the REAL type statements to DOUBLE PRECISION and hange the subroutine name from SPFIT to DPFIT. Change the names of the auxiliary subroutines mentioned above to DCP- VAL, DMPVAL, DCPDRV, DMPDRV, DCPINT, and DMPINT, respetively. July 11, 2015 Polynomial Least-Squares Curve Fit
2 C. Examples and Remarks C.1 Example Given a set of 12 data pairs (x i, y i ), ompute the weighted least-squares polynomial fit to these data letting the subroutine determine the preferred polynomial degree not to exeed 8. After omputing the leastsquares polynomial p(x), ompute and tabulate the quantities x i, y i, p(x i ), and r i = y i p(x i ). The program DRSPFIT arries out this omputation using subroutines SPFIT and SCPVAL. The output is shown in ODSPFIT. Note that the value of NDEG seleted by the subroutine is 7. C.2 Transformation of the Independent Variable For the best numerial auray, both in determining the best fitting polynomial and in the later evaluation of the polynomial, it is generally advisable to use transformation parameters P(1) and P(2) that ause the transformed variable u of Eqs. (1) or (3) to range from 1 to +1. This ondition an be obtained by setting COMTRN =.TRUE. When fitting with polynomials of degree higher than about six, the Chebyshev basis generally gives better numerial auray than the monomial basis. To obtain the potential numerial advantages of the Chebyshev basis it is essential to ause the transformed variable, u, to range from 1 to +1. To fore an identity transformation (i.e., effetively no transformation) of the independent variable, set COMTRN =.FALSE., P(1) = 0., and P(2) = 1. C.3 Interpretation of SIGFAC If the user has set all SD(i) = 1.0, or equivalently, set SD(1) = 1.0, SIGFAC an be interpreted as an a posteriori estimate of the standard deviation of the error in eah Y(i) value. More generally, whatever values the user has assigned to the SD(i) s, the a posteriori estimate of the standard deviation of the error in Y(i) is SIGFAC SD(i). C.4 Equality Constraints If the user wishes the fitted polynomial to agree to mahine auray with one or more of the data points, this an be aomplished by setting the SD() value for suh points muh smaller than for the other points. Let ɛ denote the mahine preision, i.e., ɛ = R1MACH(4) in single preision or ɛ = D1MACH(4) in double preision, see Chapter For some value of in the range , we suggest setting SD(i) = ɛ for the points at whih an exat fit is desired and SD(i) = 1.0 for all other points. For example if ɛ = 10 8 set SD(i) in the range 10 4 through 10 6 for the exat fit points. Using < 0.5 may not produe suffiiently small residuals at the desired points of exat fit, whereas setting > 0.75 may trip the near-singularity test in these subroutines leading to an unwanted alteration of the problem. Note that it is not mathematially reasonable to attempt to fore exat fits at more than NMAX+1 points. D. Funtional Desription The user provides data (x i, y i, s i ), for i = 1,..., M, and an integer NMAX. If SEEKN =.FALSE. and the problem is not rank-defiient or extremely ill-onditioned, the subroutine simply determines the polynomial p n (x) of degree n = NMAX that minimizes the following weighted sum of squares of residuals: M [ ] 2 ρ 2 yi p n (x i ) n = i=1 s i The subroutine also omputes [ ρ 2 n σ n = max(1, M n 1) ] 1/2 and returns this value as SIGFAC. The subroutine makes use of subroutines SROTMG and SROTM (or DROTMG and DROTM) from the BLAS1 (See Chapter 6.3.) to implement a sequential Modified Givens orthogonal transformation method of reduing the matrix of the least-squares problem to triangular form [1], and then solves the triangular system. This method has exellent numerial stability. After triangularizing the matrix of the problem for degree NMAX, the quantities ρ 2 n and σ 2 n may be omputed inexpensively for all degrees n from zero through NMAX. Thus, if the user sets SEEKN =.TRUE. the subroutine omputes these values of ρ 2 n and omputes CMIN = min{σ 2 n : 0 n NMAX}. The quantity ρ n is a stritly noninreasing funtion of n whereas σ n typially dereases initially with inreasing n but then osillates and slowly inreases as n beomes so large that the polynomial p n (x) is fitting noise rather than a true polynomial signal in the data. The subroutine sets NDEG to be the smallest value of n for whih σn CMIN, and it sets SIGFAC = σ NDEG. This method of seleting NDEG is generally satisfatory when the ratio NMAX/M is reasonably small, say less than 1/3. If this ratio is larger there is a tendeny for the NDEG seleted by this method to be larger than one might pik if one viewed graphs of solution polynomials of different degrees Polynomial Least-Squares Curve Fit July 11, 2015
3 Whatever the setting of SEEKN, if the problem for degree NMAX is rank-defiient, or very ill-onditioned, the subroutine will restrit onsideration to lower degrees for whih the problem is not extremely ill-onditioned. In partiular this redution of degree will ertainly our if NMAX+1 > M. Parameterization of the polynomial, p n (x) In this subroutine an n th degree polynomial p n (x) is represented in one of two speial parametri forms (monomial basis or Chebyshev basis), both of whih inlude the speifiation of a linear transformation of the independent variable in the parameterization. The parameters defining an n th degree polynomial p n (x) are stored in the array P(i), i = 1,..., n + 3. Using the monomial basis these parameters define the polynomial p n (x) as follows: u = x P (1) (1) P (2) n p n (x) = P (i + 3)u i (2) i=0 whereas, if the Chebyshev basis is requested, p n (x) is defined as: u = x P (1) (3) P (2) n p n (x) = P (i + 3)T i (u) (4) i=0 where T i (u) denotes the Chebyshev polynomial of degree i. The Chebyshev polynomials may be defined as follows: T 0 (u) = 1, T 1 (u) = u, T 2 (u) = 2u 2 1 T i (u) = 2u T i 1 (u) T i 2 (u) i = 3, 4,... The parameters P(1) and P(2) that our in the transformation formula (1) or (3) may be set by the user (COMTRN =.FALSE.) or else omputed by this subroutine (COMTRN =.TRUE.). In the latter ase P(1) and P(2) will be omputed as P(1) = (x max + x min )/2 P(2) = (x max x min )/2 where x max and x min are respetively the maximum and minimum values of x in the data array X(i), i = 1,..., M. This auses the transformed variable u of Eq. (1) or (3) to range from 1 to +1 as x ranges from x min to x max. Referenes 1. Charles L. Lawson and Rihard J. Hanson, Solving Least-Squares Problems, Prentie-Hall, Englewood Cliffs, N. J. (1974) 340 pages. E. Error Proedures and Restritions The use of the Chebyshev basis is effetive in improving the onditioning of the problem only if a transformation of the variable is used that maps the interval of interest for both fitting and evaluation to an interval approximately oinident with [ 1., 1.]. The automati proedure for seleting the degree of the fitted polynomial, used when SEEKN =.TRUE., tends to selet the degree somewhat higher in some ases than one would hoose by looking at plots of fits of different degrees. This subroutine will issue an error message using the error proessor in Chapter 19.2 at level 0, set NDEG = 1, and return if any of the following erroneous onditions exists: 1. SD(1) = 0, or M 0, or NMAX < COMTRN =.FALSE. and P(2) = SD(1) > 0. and some SD(i) 0. for 1 < i M. F. Supporting Information The soure language is Fortran 77. Entry Required Files DPFIT AMACH, DERM1, DERV1, DPFIT, DROTM, DROTMG, ERFIN, ERMSG, IERM1, IERV1 SPFIT AMACH, ERFIN, ERMSG, IERM1, IERV1, SERM1, SERV1, SPFIT, SROTM, SROTMG Initially designed and programmed by C. L. Lawson, JPL, Modified by Lawson and S. Y. Chiu, 1984, to use Modified Givens rather than Householder transformations to redue the storage requirement for W() and the omplexity of the program. Also, adapted the ode to Fortran 77. July 11, 2015 Polynomial Least-Squares Curve Fit
4 DRSPFIT program DRSPFIT >> DRSPFIT Krogh S p e i a l ode f o r C onversion. >> DRSPFIT Krogh Added e x t e r n a l statement. >> DRSPFIT Krogh Delare a l l v a r i a b l e s. >> DRSPFIT Krogh Remove 0 i n format ( ag ai n?) >> DRSPFIT Krogh Changes t o use M77CON >> DRSPFIT WVS Remove 0 i n format >> DRSPFIT CLL E d i t e d f o r Fortran 90 >> DRSPFIT Lawson I n i t i a l Code. S r e p l a e s? : DR?PFIT,?PFIT,?CPVAL Demonstration d r i v e r f o r SPFIT. ++ Code f o r.c. i s i n a t i v e %% l o n g i n t j ; ++ End external SCPVAL real X( 1 2 ),Y( 1 2 ),P( 1 1 ),SIGFAC,W( ),SCPVAL real SIG ( 1 ), R, YFIT integer I, M, NDEG, NP3 data X / 2. E0, 4. E0, 6. E0, 8. E0, 1 0. E0, 1 2. E0, 1 4. E0, 1 6. E0, 1 8. E0, 2 0. E0, 2 2. E0, 2 4. E0 / data Y / 2. 2 E0, 4. 0 E0, 5. 0 E0, 4. 6 E0, 2. 8 E0, 2. 7 E0, 3. 8 E0, 5. 1 E0, 6. 1 E0, 6. 3 E0, 5. 0 E0, 2. 0 E0 / data SIG ( 1 ) / 1.E0 / data M / 12 / 1001 format ( 1X/ I X Y YFIT R=Y YFIT /1X) 1002 format (1X, I2, F6. 0, 2 F9. 3, F10. 3 ) a l l SPFIT(M,X,Y, SIG, 8,.TRUE.,.TRUE.,.TRUE.,P,NDEG, SIGFAC,W) NP3 = NDEG+3 ++ Code f o r.c. i s a t i v e print (/ NDEG =, I2, 1 0X, SIGFAC =, F8. 4 / / P( 1 ),P( 2 ) =,9X, 2 F15. 8 / / P ( 3 ),..., P(NDEG+3) =,3 F15. 8 / (21X, 3 F15. 8 ) ), NDEG, SIGFAC, ( P( I ), I =1,NP3) ++ Code f o r.c. i s i n a t i v e %% p r i n t f ( %% \n NDEG =%2l d SIGFAC =%8.4 f, ndeg, s i g f a ) ; %% p r i n t f ( %% \n\n P( 1 ),P(2) = %15.8 f %15.5 f \n\n P ( 3 ),..., P(NDEG+3) =, %% p [ 0 ], p [ 1 ] ) ; %% f o r ( i = 2 ; i < ( np3 ) ; i +=3){ %% f o r ( j = i ; j < ( i < np3 2? i+3 : np3 ) ; j++) %% p r i n t f ( %15.8 f, p [ j ] ) ; %% i f ( i < np3 2) p r i n t f ( \n ) ; } %% p r i n t f ( \n ) ; ++ End write (, ) do 20 I = 1,M YFIT = SCPVAL(P,NDEG,X( I ) ) R = Y( I ) YFIT write (, ) I,X( I ),Y( I ),YFIT,R 20 ontinue stop end Polynomial Least-Squares Curve Fit July 11, 2015
5 ODSPFIT NDEG = 7 SIGFAC = P( 1 ),P( 2 ) = P ( 3 ),..., P(NDEG+3) = I X Y YFIT R=Y YFIT July 11, 2015 Polynomial Least-Squares Curve Fit
7.1 Roots of a Polynomial
7.1 Roots of a Polynomial A. Purpose Given the oeffiients a i of a polynomial of degree n = NDEG > 0, a 1 z n + a 2 z n 1 +... + a n z + a n+1 with a 1 0, this subroutine omputes the NDEG roots of the
More information3.2 Gaussian (Normal) Random Numbers and Vectors
3.2 Gaussian (Normal) Random Numbers and Vetors A. Purpose Generate pseudorandom numbers or vetors from the Gaussian (normal) distribution. B. Usage B.1 Generating Gaussian (normal) pseudorandom numbers
More information4.3 Singular Value Decomposition and Analysis
4.3 Singular Value Deomposition and Analysis A. Purpose Any M N matrix, A, has a Singular Value Deomposition (SVD) of the form A = USV t where U is an M M orthogonal matrix, V is an N N orthogonal matrix,
More information2.9 Incomplete Elliptic Integrals
.9 Inomplete Ellipti Integrals A. Purpose An integral of the form R t, P t) 1/) dt 1) in whih Pt) is a polynomial of the third or fourth degree that has no multiple roots, and R is a rational funtion of
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More information6.3 Basic Linear Algebra Subprograms (BLAS1)
6.3 Basi Linear Algebra Subprograms (BLAS1) A. Purpose This is a set of subroutine and funtion subprograms for basi mathematial operations on a single vetor or a pair of vetors. The operations provided
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationAn iterative least-square method suitable for solving large sparse matrices
An iteratie least-square method suitable for soling large sparse matries By I. M. Khabaza The purpose of this paper is to report on the results of numerial experiments with an iteratie least-square method
More informationMethods of evaluating tests
Methods of evaluating tests Let X,, 1 Xn be i.i.d. Bernoulli( p ). Then 5 j= 1 j ( 5, ) T = X Binomial p. We test 1 H : p vs. 1 1 H : p>. We saw that a LRT is 1 if t k* φ ( x ) =. otherwise (t is the observed
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationComputer Science 786S - Statistical Methods in Natural Language Processing and Data Analysis Page 1
Computer Siene 786S - Statistial Methods in Natural Language Proessing and Data Analysis Page 1 Hypothesis Testing A statistial hypothesis is a statement about the nature of the distribution of a random
More informationSingular Event Detection
Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More information8.2 Solve System of Nonlinear Equations
8.2 Solve System of Nonlinear Equations A. Purpose This subroutine solves a system of n nonlinear equations in n unknowns. The problem can be expressed as f(x) = 0 where x denotes an n-vector of unknowns
More informationDanielle Maddix AA238 Final Project December 9, 2016
Struture and Parameter Learning in Bayesian Networks with Appliations to Prediting Breast Caner Tumor Malignany in a Lower Dimension Feature Spae Danielle Maddix AA238 Final Projet Deember 9, 2016 Abstrat
More informationUTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.
UTC Engineering 329 Proportional Controller Design for Speed System By John Beverly Green Team John Beverly Keith Skiles John Barker 24 Mar 2006 Introdution This experiment is intended test the variable
More informationThe Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationAcoustic Waves in a Duct
Aousti Waves in a Dut 1 One-Dimensional Waves The one-dimensional wave approximation is valid when the wavelength λ is muh larger than the diameter of the dut D, λ D. The aousti pressure disturbane p is
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationTaste for variety and optimum product diversity in an open economy
Taste for variety and optimum produt diversity in an open eonomy Javier Coto-Martínez City University Paul Levine University of Surrey Otober 0, 005 María D.C. Garía-Alonso University of Kent Abstrat We
More informationController Design Based on Transient Response Criteria. Chapter 12 1
Controller Design Based on Transient Response Criteria Chapter 12 1 Desirable Controller Features 0. Stable 1. Quik responding 2. Adequate disturbane rejetion 3. Insensitive to model, measurement errors
More informationElectromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.
arxiv:physis/99536v1 [physis.lass-ph] 15 May 1999 Eletromagneti radiation of the travelling spin wave propagating in an antiferromagneti plate. Exat solution. A.A.Zhmudsky November 19, 16 Abstrat The exat
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationCase I: 2 users In case of 2 users, the probability of error for user 1 was earlier derived to be 2 A1
MUTLIUSER DETECTION (Letures 9 and 0) 6:33:546 Wireless Communiations Tehnologies Instrutor: Dr. Narayan Mandayam Summary By Shweta Shrivastava (shwetash@winlab.rutgers.edu) bstrat This artile ontinues
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationA Queueing Model for Call Blending in Call Centers
A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl
More informationUNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS. V. N. Matveev and O. V. Matvejev
UNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS V. N. Matveev and O. V. Matvejev Joint-Stok Company Sinerta Savanoriu pr., 159, Vilnius, LT-315, Lithuania E-mail: matwad@mail.ru Abstrat
More informationChapter 8 Hypothesis Testing
Leture 5 for BST 63: Statistial Theory II Kui Zhang, Spring Chapter 8 Hypothesis Testing Setion 8 Introdution Definition 8 A hypothesis is a statement about a population parameter Definition 8 The two
More informationNonreversibility of Multiple Unicast Networks
Nonreversibility of Multiple Uniast Networks Randall Dougherty and Kenneth Zeger September 27, 2005 Abstrat We prove that for any finite direted ayli network, there exists a orresponding multiple uniast
More informationA Functional Representation of Fuzzy Preferences
Theoretial Eonomis Letters, 017, 7, 13- http://wwwsirporg/journal/tel ISSN Online: 16-086 ISSN Print: 16-078 A Funtional Representation of Fuzzy Preferenes Susheng Wang Department of Eonomis, Hong Kong
More informationLightpath routing for maximum reliability in optical mesh networks
Vol. 7, No. 5 / May 2008 / JOURNAL OF OPTICAL NETWORKING 449 Lightpath routing for maximum reliability in optial mesh networks Shengli Yuan, 1, * Saket Varma, 2 and Jason P. Jue 2 1 Department of Computer
More informationThe Impact of Information on the Performance of an M/M/1 Queueing System
The Impat of Information on the Performane of an M/M/1 Queueing System by Mojgan Nasiri A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master
More informationSimple FIR Digital Filters. Simple FIR Digital Filters. Simple Digital Filters. Simple FIR Digital Filters. Simple FIR Digital Filters
Simple Digital Filters Later in the ourse we shall review various methods of designing frequeny-seletive filters satisfying presribed speifiations We now desribe several low-order FIR and IIR digital filters
More informationOptimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach
Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)
More informationVerka Prolović Chair of Civil Engineering Geotechnics, Faculty of Civil Engineering and Architecture, Niš, R. Serbia
3 r d International Conferene on New Developments in Soil Mehanis and Geotehnial Engineering, 8-30 June 01, Near East University, Niosia, North Cyprus Values of of partial fators for for EC EC 7 7 slope
More information2. The Energy Principle in Open Channel Flows
. The Energy Priniple in Open Channel Flows. Basi Energy Equation In the one-dimensional analysis of steady open-hannel flow, the energy equation in the form of Bernoulli equation is used. Aording to this
More informationmax min z i i=1 x j k s.t. j=1 x j j:i T j
AM 221: Advaned Optimization Spring 2016 Prof. Yaron Singer Leture 22 April 18th 1 Overview In this leture, we will study the pipage rounding tehnique whih is a deterministi rounding proedure that an be
More information7 Max-Flow Problems. Business Computing and Operations Research 608
7 Max-Flow Problems Business Computing and Operations Researh 68 7. Max-Flow Problems In what follows, we onsider a somewhat modified problem onstellation Instead of osts of transmission, vetor now indiates
More informationModel-based mixture discriminant analysis an experimental study
Model-based mixture disriminant analysis an experimental study Zohar Halbe and Mayer Aladjem Department of Eletrial and Computer Engineering, Ben-Gurion University of the Negev P.O.Box 653, Beer-Sheva,
More informationIMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION
09-1289 Citation: Brilon, W. (2009): Impedane Effets of Left Turners from the Major Street at A TWSC Intersetion. Transportation Researh Reord Nr. 2130, pp. 2-8 IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationCONDITIONAL CONFIDENCE INTERVAL FOR THE SCALE PARAMETER OF A WEIBULL DISTRIBUTION. Smail Mahdi
Serdia Math. J. 30 (2004), 55 70 CONDITIONAL CONFIDENCE INTERVAL FOR THE SCALE PARAMETER OF A WEIBULL DISTRIBUTION Smail Mahdi Communiated by N. M. Yanev Abstrat. A two-sided onditional onfidene interval
More informationarxiv:hep-ph/ v2 30 May 1998
Ref. SISSA 31/98/EP hep ph/9805262 8 May, 1998 Diffrative-Like (or Parametri-Resonane-Like?) Enhanement of the Earth (Day-Night) Effet arxiv:hep-ph/9805262v2 30 May 1998 for Solar Neutrinos Crossing the
More informationPart G-4: Sample Exams
Part G-4: Sample Exams 1 Cairo University M.S.: Eletronis Cooling Faulty of Engineering Final Exam (Sample 1) Mehanial Power Engineering Dept. Time allowed 2 Hours Solve as muh as you an. 1. A heat sink
More informationA NETWORK SIMPLEX ALGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM
NETWORK SIMPLEX LGORITHM FOR THE MINIMUM COST-BENEFIT NETWORK FLOW PROBLEM Cen Çalışan, Utah Valley University, 800 W. University Parway, Orem, UT 84058, 801-863-6487, en.alisan@uvu.edu BSTRCT The minimum
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationMATHEMATICAL AND NUMERICAL BASIS OF BINARY ALLOY SOLIDIFICATION MODELS WITH SUBSTITUTE THERMAL CAPACITY. PART II
Journal of Applied Mathematis and Computational Mehanis 2014, 13(2), 141-147 MATHEMATICA AND NUMERICA BAI OF BINARY AOY OIDIFICATION MODE WITH UBTITUTE THERMA CAPACITY. PART II Ewa Węgrzyn-krzypzak 1,
More informationSimplified Buckling Analysis of Skeletal Structures
Simplified Bukling Analysis of Skeletal Strutures B.A. Izzuddin 1 ABSRAC A simplified approah is proposed for bukling analysis of skeletal strutures, whih employs a rotational spring analogy for the formulation
More informationCMSC 451: Lecture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017
CMSC 451: Leture 9 Greedy Approximation: Set Cover Thursday, Sep 28, 2017 Reading: Chapt 11 of KT and Set 54 of DPV Set Cover: An important lass of optimization problems involves overing a ertain domain,
More informationOrthogonal Complement Based Divide-and-Conquer Algorithm (O-DCA) for Constrained Multibody Systems
Orthogonal Complement Based Divide-and-Conquer Algorithm (O-DCA) for Constrained Multibody Systems Rudranarayan M. Mukherjee, Kurt S. Anderson Computational Dynamis Laboratory Department of Mehanial Aerospae
More informationProcess engineers are often faced with the task of
Fluids and Solids Handling Eliminate Iteration from Flow Problems John D. Barry Middough, In. This artile introdues a novel approah to solving flow and pipe-sizing problems based on two new dimensionless
More informationSensitivity Analysis in Markov Networks
Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores
More informationDIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS
CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install
More informationRelativity in Classical Physics
Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of
More informationInternational Journal of Advanced Engineering Research and Studies E-ISSN
Researh Paper FINIE ELEMEN ANALYSIS OF A CRACKED CANILEVER BEAM Mihir Kumar Sutar Address for Correspondene Researh Sholar, Department of Mehanial & Industrial Engineering Indian Institute of ehnology
More informationDirectional Coupler. 4-port Network
Diretional Coupler 4-port Network 3 4 A diretional oupler is a 4-port network exhibiting: All ports mathed on the referene load (i.e. S =S =S 33 =S 44 =0) Two pair of ports unoupled (i.e. the orresponding
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationCombined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2
Sensor and Simulation Notes Note 53 3 May 8 Combined Eletri and Magneti Dipoles for Mesoband Radiation, Part Carl E. Baum University of New Mexio Department of Eletrial and Computer Engineering Albuquerque
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationRobust Flight Control Design for a Turn Coordination System with Parameter Uncertainties
Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo
More informationSufficient Conditions for a Flexible Manufacturing System to be Deadlocked
Paper 0, INT 0 Suffiient Conditions for a Flexile Manufaturing System to e Deadloked Paul E Deering, PhD Department of Engineering Tehnology and Management Ohio University deering@ohioedu Astrat In reent
More informationSpeed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion
Speed-feedbak Diret-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion Y. Yamamoto, T. Nakamura 2, Y. Takada, T. Koseki, Y. Aoyama 3, and Y. Iwaji 3
More informationThe universal model of error of active power measuring channel
7 th Symposium EKO TC 4 3 rd Symposium EKO TC 9 and 5 th WADC Workshop nstrumentation for the CT Era Sept. 8-2 Kosie Slovakia The universal model of error of ative power measuring hannel Boris Stogny Evgeny
More informationWhat are the locations of excess energy in open channels?
Leture 26 Energy Dissipation Strutures I. Introdution Exess energy should usually be dissipated in suh a way as to avoid erosion in unlined open hannels In this ontext, exess energy means exess water veloity
More informationEE 321 Project Spring 2018
EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,
More informationModeling of discrete/continuous optimization problems: characterization and formulation of disjunctions and their relaxations
Computers and Chemial Engineering (00) 4/448 www.elsevier.om/loate/omphemeng Modeling of disrete/ontinuous optimization problems: haraterization and formulation of disjuntions and their relaxations Aldo
More informationREFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS. 1. Introduction
Version of 5/2/2003 To appear in Advanes in Applied Mathematis REFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS MATTHIAS BECK AND SHELEMYAHU ZACKS Abstrat We study the Frobenius problem:
More informationMicroeconomic Theory I Assignment #7 - Answer key
Miroeonomi Theory I Assignment #7 - Answer key. [Menu priing in monopoly] Consider the example on seond-degree prie disrimination (see slides 9-93). To failitate your alulations, assume H = 5, L =, and
More informationINFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE BEHAVIOR OF HYDRAULIC CYLINDER SUBJECTED TO JERKY MOTION
Proeedings of ICFDP 8: 8 th International Congress of Fluid Dynamis & Propulsion Deember 14-17, 006, Sharm El-Shiekh, Sinai, Egypt ICFDP8-EG-154 INFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE
More informationMath 151 Introduction to Eigenvectors
Math 151 Introdution to Eigenvetors The motivating example we used to desrie matrixes was landsape hange and vegetation suession. We hose the simple example of Bare Soil (B), eing replaed y Grasses (G)
More informationTHE METHOD OF SECTIONING WITH APPLICATION TO SIMULATION, by Danie 1 Brent ~~uffman'i
THE METHOD OF SECTIONING '\ WITH APPLICATION TO SIMULATION, I by Danie 1 Brent ~~uffman'i Thesis submitted to the Graduate Faulty of the Virginia Polytehni Institute and State University in partial fulfillment
More informationarxiv: v2 [math.pr] 9 Dec 2016
Omnithermal Perfet Simulation for Multi-server Queues Stephen B. Connor 3th Deember 206 arxiv:60.0602v2 [math.pr] 9 De 206 Abstrat A number of perfet simulation algorithms for multi-server First Come First
More informationEXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION
Journal of Mathematial Sienes: Advanes and Appliations Volume 3, 05, Pages -3 EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION JIAN YANG, XIAOJUAN LU and SHENGQIANG TANG
More informationSURFACE WAVES OF NON-RAYLEIGH TYPE
SURFACE WAVES OF NON-RAYLEIGH TYPE by SERGEY V. KUZNETSOV Institute for Problems in Mehanis Prosp. Vernadskogo, 0, Mosow, 75 Russia e-mail: sv@kuznetsov.msk.ru Abstrat. Existene of surfae waves of non-rayleigh
More informationSINCE Zadeh s compositional rule of fuzzy inference
IEEE TRANSACTIONS ON FUZZY SYSTEMS, VOL. 14, NO. 6, DECEMBER 2006 709 Error Estimation of Perturbations Under CRI Guosheng Cheng Yuxi Fu Abstrat The analysis of stability robustness of fuzzy reasoning
More informationOn the Bit Error Probability of Noisy Channel Networks With Intermediate Node Encoding I. INTRODUCTION
5188 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 54, NO. 11, NOVEMBER 2008 [8] A. P. Dempster, N. M. Laird, and D. B. Rubin, Maximum likelihood estimation from inomplete data via the EM algorithm, J.
More informationSurging in Coil Springs
Purdue University Purdue e-pubs nternational Compressor Engineering Conferene Shool of Mehanial Engineering 1996 Surging in Coil Springs R. A. Simmons Purdue University W. Soedel Purdue University Follow
More informationOn Component Order Edge Reliability and the Existence of Uniformly Most Reliable Unicycles
Daniel Gross, Lakshmi Iswara, L. William Kazmierzak, Kristi Luttrell, John T. Saoman, Charles Suffel On Component Order Edge Reliability and the Existene of Uniformly Most Reliable Uniyles DANIEL GROSS
More informationThe shape of a hanging chain. a project in the calculus of variations
The shape of a hanging hain a projet in the alulus of variations April 15, 218 2 Contents 1 Introdution 5 2 Analysis 7 2.1 Model............................... 7 2.2 Extremal graphs.........................
More informationAnalysis of discretization in the direct simulation Monte Carlo
PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 16 Aug 2004
Computational omplexity and fundamental limitations to fermioni quantum Monte Carlo simulations arxiv:ond-mat/0408370v1 [ond-mat.stat-meh] 16 Aug 2004 Matthias Troyer, 1 Uwe-Jens Wiese 2 1 Theoretishe
More informationError Bounds for Context Reduction and Feature Omission
Error Bounds for Context Redution and Feature Omission Eugen Bek, Ralf Shlüter, Hermann Ney,2 Human Language Tehnology and Pattern Reognition, Computer Siene Department RWTH Aahen University, Ahornstr.
More informationFour-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field
Four-dimensional equation of motion for visous ompressible substane with regard to the aeleration field, pressure field and dissipation field Sergey G. Fedosin PO box 6488, Sviazeva str. -79, Perm, Russia
More informationCollinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b
International Frontier Siene Letters Submitted: 6-- ISSN: 9-8, Vol., pp -6 Aepted: -- doi:.8/www.sipress.om/ifsl.. Online: --8 SiPress Ltd., Switzerland Collinear Equilibrium Points in the Relativisti
More informationA Characterization of Wavelet Convergence in Sobolev Spaces
A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat
More informationApplying CIECAM02 for Mobile Display Viewing Conditions
Applying CIECAM2 for Mobile Display Viewing Conditions YungKyung Park*, ChangJun Li*, M.. Luo*, Youngshin Kwak**, Du-Sik Park **, and Changyeong Kim**; * University of Leeds, Colour Imaging Lab, UK*, **
More informationA NORMALIZED EQUATION OF AXIALLY LOADED PILES IN ELASTO-PLASTIC SOIL
Journal of Geongineering, Vol. Yi-Chuan 4, No. 1, Chou pp. 1-7, and April Yun-Mei 009 Hsiung: A Normalized quation of Axially Loaded Piles in lasto-plasti Soil 1 A NORMALIZD QUATION OF AXIALLY LOADD PILS
More informationArray Design for Superresolution Direction-Finding Algorithms
Array Design for Superresolution Diretion-Finding Algorithms Naushad Hussein Dowlut BEng, ACGI, AMIEE Athanassios Manikas PhD, DIC, AMIEE, MIEEE Department of Eletrial Eletroni Engineering Imperial College
More informationNormative and descriptive approaches to multiattribute decision making
De. 009, Volume 8, No. (Serial No.78) China-USA Business Review, ISSN 57-54, USA Normative and desriptive approahes to multiattribute deision making Milan Terek (Department of Statistis, University of
More informationTutorial 4 (week 4) Solutions
THE UNIVERSITY OF SYDNEY PURE MATHEMATICS Summer Shool Math26 28 Tutorial week s You are given the following data points: x 2 y 2 Construt a Lagrange basis {p p p 2 p 3 } of P 3 using the x values from
More informationModeling Probabilistic Measurement Correlations for Problem Determination in Large-Scale Distributed Systems
009 9th IEEE International Conferene on Distributed Computing Systems Modeling Probabilisti Measurement Correlations for Problem Determination in Large-Sale Distributed Systems Jing Gao Guofei Jiang Haifeng
More informationTail-robust Scheduling via Limited Processor Sharing
Performane Evaluation Performane Evaluation 00 200) 22 Tail-robust Sheduling via Limited Proessor Sharing Jayakrishnan Nair a, Adam Wierman b, Bert Zwart a Department of Eletrial Engineering, California
More informationSubject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (
16.50 Leture 0 Subjet: Introdution to Component Mathing and Off-Design Operation At this point it is well to reflet on whih of the many parameters we have introdued (like M, τ, τ t, ϑ t, f, et.) are free
More informationONLINE APPENDICES for Cost-Effective Quality Assurance in Crowd Labeling
ONLINE APPENDICES for Cost-Effetive Quality Assurane in Crowd Labeling Jing Wang Shool of Business and Management Hong Kong University of Siene and Tehnology Clear Water Bay Kowloon Hong Kong jwang@usthk
More informationAharonov-Bohm effect. Dan Solomon.
Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that
More informationA Combination of BIM and BEM for efficiently analyzing optical. elements
A Combination of BIM and BEM for effiiently analyzing optial elements SUN Fang, LIU Juan, ZHANG Guo-ting, HU Chuan-fei, SU Xiaoxing Institute of Lightwave Tehnology, hool of eletronis and information engineering,
More information