Fitting a Polynomial to Heat Capacity as a Function of Temperature for Ag. Mathematical Background Document

Similar documents
Rank One Update And the Google Matrix by Al Bernstein Signal Science, LLC

Lecture 4: Piecewise Cubic Interpolation

Solubilities and Thermodynamic Properties of SO 2 in Ionic

Applied Statistics Qualifier Examination

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II

INTRODUCTION TO COMPLEX NUMBERS

Statistics 423 Midterm Examination Winter 2009

6 Roots of Equations: Open Methods

Jens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers

Quiz: Experimental Physics Lab-I

APPROXIMATE SOLUTION OF HIGHER ORDER LINEAR DIFFERENTIAL EQUATIONS BY MEANS OF A NEW RATIONAL CHEBYSHEV COLLOCATION METHOD

Exponents and Powers

Principle Component Analysis

CISE 301: Numerical Methods Lecture 5, Topic 4 Least Squares, Curve Fitting

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Lecture 3 Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab

Katholieke Universiteit Leuven Department of Computer Science

Chapter 5 Supplemental Text Material R S T. ij i j ij ijk

The Number of Rows which Equal Certain Row

Course Review Introduction to Computer Methods

OXFORD H i g h e r E d u c a t i o n Oxford University Press, All rights reserved.

GAUSS ELIMINATION. Consider the following system of algebraic linear equations

Least Squares Fitting of Data

Two Coefficients of the Dyson Product

System in Weibull Distribution

Numbers Related to Bernoulli-Goss Numbers

Least Squares Fitting of Data

Lecture 8: Camera Calibration

Chapter 2 Introduction to Algebra. Dr. Chih-Peng Li ( 李 )

Special Relativity and Riemannian Geometry. Department of Mathematical Sciences

The Schur-Cohn Algorithm

Least squares. Václav Hlaváč. Czech Technical University in Prague

Review of linear algebra. Nuno Vasconcelos UCSD

Applied Mathematics Letters

Discussion Introduction P212, Week 1 The Scientist s Sixth Sense. Knowing what the answer will look like before you start.

Definition of Tracking

1 Linear Least Squares

An Ising model on 2-D image

ECONOMETRIC THEORY. MODULE IV Lecture - 16 Predictions in Linear Regression Model

Model Fitting and Robust Regression Methods

Remember: Project Proposals are due April 11.

A Family of Multivariate Abel Series Distributions. of Order k

Chapter 12 Lyes KADEM [Thermodynamics II] 2007

The graphs of Rational Functions

Demand. Demand and Comparative Statics. Graphically. Marshallian Demand. ECON 370: Microeconomic Theory Summer 2004 Rice University Stanley Gilbert

MA 131 Lecture Notes Calculus Sections 1.5 and 1.6 (and other material)

arxiv:math/ v2 [math.ho] 16 Dec 2003

Department of Mechanical Engineering, University of Bath. Mathematics ME Problem sheet 11 Least Squares Fitting of data

r = cos θ + 1. dt ) dt. (1)

INTERPOLATION(1) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

Denote the function derivatives f(x) in given points. x a b. Using relationships (1.2), polynomials (1.1) are written in the form

MATHEMATICAL MODEL AND STATISTICAL ANALYSIS OF THE TENSILE STRENGTH (Rm) OF THE STEEL QUALITY J55 API 5CT BEFORE AND AFTER THE FORMING OF THE PIPES

6.6 The Marquardt Algorithm

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Creative Practicing. By Jimmy Wyble edited by David Oakes

CURVE FITTING LEAST SQUARES METHOD

The Regulated and Riemann Integrals

8. INVERSE Z-TRANSFORM

Math 8 Winter 2015 Applications of Integration

4.4 Areas, Integrals and Antiderivatives

Research Article On the Upper Bounds of Eigenvalues for a Class of Systems of Ordinary Differential Equations with Higher Order

Overview of Calculus I

Variable time amplitude amplification and quantum algorithms for linear algebra. Andris Ambainis University of Latvia

Many-Body Calculations of the Isotope Shift

8 Transportation Problem Alpha-Beta

NUMERICAL MODELLING OF A CILIUM USING AN INTEGRAL EQUATION

The Atwood Machine OBJECTIVE INTRODUCTION APPARATUS THEORY

COMPLEX NUMBERS INDEX

econstor Make Your Publication Visible

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

Lesson 1: Quadratic Equations

Lecture 36. Finite Element Methods

,..., k N. , k 2. ,..., k i. The derivative with respect to temperature T is calculated by using the chain rule: & ( (5) dj j dt = "J j. k i.

Expectation and Variance

Proof that if Voting is Perfect in One Dimension, then the First. Eigenvector Extracted from the Double-Centered Transformed

Math Lecture 23

COMPLEX NUMBER & QUADRATIC EQUATION

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Lecture 8: Camera Calibra0on

Chapters 4 & 5 Integrals & Applications

Analysis of Geometric, Zernike and United Moment Invariants Techniques Based on Intra-class Evaluation

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

Introduction to Numerical Integration Part II

Predict Global Earth Temperature using Linier Regression

LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER

Recitation 3: More Applications of the Derivative

Homework Solution - Set 5 Due: Friday 10/03/08

4. More general extremum principles and thermodynamic potentials

Chapter 0. What is the Lebesgue integral about?

Before we can begin Ch. 3 on Radicals, we need to be familiar with perfect squares, cubes, etc. Try and do as many as you can without a calculator!!!

4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve

MTH 122 Fall 2008 Essex County College Division of Mathematics Handout Version 10 1 October 14, 2008

CHOVER-TYPE LAWS OF THE ITERATED LOGARITHM FOR WEIGHTED SUMS OF ρ -MIXING SEQUENCES

THE COMBINED SHEPARD ABEL GONCHAROV UNIVARIATE OPERATOR

Improper Integrals, and Differential Equations

MAE140 - Linear Circuits - Fall 10 Midterm, October 28

V E L O C I T Y a n d V E L O C I T Y P R E S S U R E I n A I R S Y S T E M S

Transcription:

Fttng Polynol to Het Cpcty s Functon of Teperture for Ag. thetcl Bckground Docuent by Theres Jul Zelnsk Deprtent of Chestry, edcl Technology, nd Physcs onouth Unversty West ong Brnch, J 7764-898 tzelns@onouth.edu Theres Jul Zelnsk, 996, 997, 998. You re welcoe to use ths docuent n your own clsses but coercl use s not llowed wthout the persson of the uthor. Fttng thetcl functon to experentl dt s useful technque for representng dt wth n nlytcl functon. To do ths we need to dentfy the type of functon nd procedure for optzng the fttng preters n the functon. Polynol expressons y be sutble choce n ny cses. Polynols lso provde n esy begnnng for exnng the bsc procedure for fttng functon to set of dt. The obectve s to choose polynol expresson such tht the su of the squres of the devtons between the functon nd the dt re nzed. Here we wll use the flr lest squres procedure tht ost students lern bout n the physcl chestry lbortory. The thetcs for ths s reltvely strght forwrd nd cn be found n ost physcl chestry texts or lbortory nuls. Consder polynol f ( x) x x x x () where,,, up to re the fttng preters, nd s the order of the polynol. If you hve set of dt consstng of x nd y vlues, then you cn clculte n f(x) for ech x vlue,.e. f(x ). The f(x ) cn then be copred to the y vlues. If the f(x ) nd y gree for every x then you hve functon f(x) tht fts the dt. We use sttstcs to deterne the degree of goodnes of ths ft. A frst pproch s to use the ethod of lest squres. Ths s the ethod tht wll be developed n these notes. Exercse : For qudrtc functon the coeffcents re.5,., nd -4.5x -4. Copute f(x ) for the followng prs of x, y dt: 3, 4.8; 5,.7; 7, 6.3; 9, 9.. The ethod of lest squres requres tht the su of the squres of the devtons,.e. squre of the resduls, be nzed. Ths s defned s: [ ] R y f( x ) () where f(x ) s the vlue of f(x) for dt pont nd there re dt ponts n the experentl set. The experentl dt re tbulted n prs of x nd y vlues. f(x ) s the coputed vlue of the property y nd t s copred to the experentl vlue y to obtn the resdul. R s the resdul functon. It s the su of the squres of the devtons nd s soetes clled SSD n thcd docuents. It doesn t tke Creted: Februry 996 fttng.doc Author: Theres Jul Zelnsk odfed: July, 998 Pge

long to get used to the for of ths functon so tht you wll be ble to recognze t no tter wht ne s hs. The dscusson tht follows does not consder the experentl weght of the dt ponts y. Exercse : Copute R for the dt n Exercse.. nzton of the resduls s ccoplshed by stndrd pplcton of clculus. We wnt to fnd the best set of fttng preters, nely those whch nze R. Therefore we tke the dervtve wth respect to ech fttng preter nd set t equl to zero s shown n equton 3.,,, etc. (3) The result s set of sultneous equtons, one for ech coeffcent. Soe of these equtons re: ( y x x ) ( ) x x ( ) x x (4) where there s one equton for ech fttng preter,.e. equtons becuse the powers of x rnge fro to n the polynol. The suton s tken over the,, 3, dt ponts. In thcd ths y pper s,,, 3. - dependng on how the rry nuberng s set n the docuent you re usng. Exercse 3: Show tht the equtons n (4) coe fro the equtons n (3). Rerrngng (4) gves: ( x x ) 3 ( x x x ) 3 ( ) y y x x x x y x 4 (5) Systes of lner equtons lke ths re best represented by trx notton. Creted: Februry 996 fttng.doc Author: Theres Jul Zelnsk odfed: July, 998 Pge

Creted: Februry 996 fttng.doc Author: Theres Jul Zelnsk odfed: July, 998 Pge 3 A b (6) otce the bold font on the coponents of equton (6). Ths ndctes tht these re trces. trx A s squre syetrc rel trx. b nd re colun vectors. We see ther explct for here n equtons 7 nd 8. Proble: Check your skll wth trx npultons by showng tht the trx equton (6) corresponds to representton of the set of sultneous equtons shown n equton (5). Hnt: use the trx nd vectors n equton (7) nd (8). A 3 (7) nd b y (8) The eleents n trx A re esly wrtten by usng the followng expresson:

A,k ( x ) ( x ) k (9) where s the row nd k s the colun. nd k correspond to the subscrpt of the coeffcent used n the polynol expresson. Do not confuse ths wth the rnge nteger,, tht dentfes the dt ponts. nd k rnge fro to,.e. over the totl nuber of fttng preters. The defult ntl ndex for vectors nd rrys s zero n thcd. If our dt were represented n ths defult then would lso run fro to - nd the trx n equton (7) nd the vectors n equton (8) would need to be rewrtten. It s portnt to be flexble but ccurte wth our use of the rnge vrble. In lke fshon we wrte the vector b n equton () ( ) () b y x where b s row eleent n the colun vector b. There re eleents n the vector b nd eleents n the vector. The trx A s by squre. Through trx lgebr we cn solve for the vector A. Be creful bout how you ultply trces. A - A A b A b () The vector contns ll of the fttng preters n the polynol expnson. We obtn the vrnce of the ft nd the vrnce of the fttng preters by: R σ σ σ A () where,,... nd R s the su of the squres of the devtons t the end of the fttng process. Queston : Wht crter would you use to deterne whch polynol you would choose to represent set of dt? The docuent Slver.cd s thcd pleentton of polynol curve fttng. Ths docuent, the thcd docuent, nd sple dt re t the thcd WWW ste or dskette provded by your nstructor. The thcd docuent shows vrous opertons but you should frst experent wth vryng the order of the polynol nd observng the ft of the functon to the dt. Record your observtons n your note book. After you hve deonstrted copetence wth the requred obectves of ths lesson you Creted: Februry 996 fttng.doc Author: Theres Jul Zelnsk odfed: July, 998 Pge 4

ght lke to try other sets of dt or functons to lern ore bout curve fttng. Other thcd docuents on the WWW ste deonstrte lterntve curve fttng pproches. Queston : How does the goodness of ft prove wth ncresng order of the fttng polynol? Queston 3: How does one choose the best polynol to represent the dt? Queston 4: Use your textbook or nother text to copre the functons used to represent the het cpcty of terls s functon of teperture. Queston 5: Wht s the crter would you use for choosng to use one functon or nother to use when fttng polynol to het cpcty dt? References: Johnson, K. Jeffrey, uercl ethods n Chestry ; ew York: rcel Dekker, Inc., 98. Bevngton, Phlp R., Dt Reducton nd Error Anlyss for the Physcl Scences ; ew York: cgrw Hll Book Copny, 969. Hng, R. W., uercl ethods for Scentsts nd Engneers, nd Ed.;ew York: Dover Publctons, 986. Acknowledgent: Prtl support for ths work ws provded to TJZ by the tonl Scence Foundton's Dvson of Undergrdute Educton through grnt DUE #9354473 nd by the ew Trdtons proect t the Unversty of Wsconsn - dson through the tonl Scence Foundton's Dvson of Undergrdute Educton through grnt DUE #945598. Creted: Februry 996 fttng.doc Author: Theres Jul Zelnsk odfed: July, 998 Pge 5