ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS

Similar documents
Chapter Gauss-Seidel Method

12 Iterative Methods. Linear Systems: Gauss-Seidel Nonlinear Systems Case Study: Chemical Reactions

Lecture 3-4 Solutions of System of Linear Equations

MATRIX AND VECTOR NORMS

CS537. Numerical Analysis

CURVE FITTING LEAST SQUARES METHOD

ME 501A Seminar in Engineering Analysis Page 1

Advanced Algorithmic Problem Solving Le 3 Arithmetic. Fredrik Heintz Dept of Computer and Information Science Linköping University

SOLVING SYSTEMS OF EQUATIONS, DIRECT METHODS

CS321. Introduction to Numerical Methods

DATA FITTING. Intensive Computation 2013/2014. Annalisa Massini

INTRODUCTION ( ) 1. Errors

Level-2 BLAS. Matrix-Vector operations with O(n 2 ) operations (sequentially) BLAS-Notation: S --- single precision G E general matrix M V --- vector

Asynchronous Analogs of Iterative Methods and Their Efficiency for Parallel Computers

CS473-Algorithms I. Lecture 3. Solving Recurrences. Cevdet Aykanat - Bilkent University Computer Engineering Department

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America

i+1 by A and imposes Ax

MTH 146 Class 7 Notes

Area and the Definite Integral. Area under Curve. The Partition. y f (x) We want to find the area under f (x) on [ a, b ]

Chapter 2 Solving Linear Equation

Optimality of Strategies for Collapsing Expanded Random Variables In a Simple Random Sample Ed Stanek

Approximate Integration

POWERS OF COMPLEX PERSYMMETRIC ANTI-TRIDIAGONAL MATRICES WITH CONSTANT ANTI-DIAGONALS

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

In Calculus I you learned an approximation method using a Riemann sum. Recall that the Riemann sum is

The linear system. The problem: solve

Rendering Equation. Linear equation Spatial homogeneous Both ray tracing and radiosity can be considered special case of this general eq.

G S Power Flow Solution

PubH 7405: REGRESSION ANALYSIS REGRESSION IN MATRIX TERMS

The z-transform. LTI System description. Prof. Siripong Potisuk

Methods to Invert a Matrix

Computer Programming

Chapter 2 Intro to Math Techniques for Quantum Mechanics

On Several Inequalities Deduced Using a Power Series Approach

Soo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7: Figure 8: Figure 9: Figure 10: Figure 11:

Sequences and summations

Autar Kaw Benjamin Rigsby. Transforming Numerical Methods Education for STEM Undergraduates

An Alternative Method to Find the Solution of Zero One Integer Linear Fractional Programming Problem with the Help of -Matrix

14.2 Line Integrals. determines a partition P of the curve by points Pi ( xi, y

Stats & Summary

Numerical differentiation (Finite difference)

International Journal of Scientific and Research Publications, Volume 3, Issue 5, May ISSN

Bond Additive Modeling 5. Mathematical Properties of the Variable Sum Exdeg Index

An Indian Journal FULL PAPER ABSTRACT KEYWORDS. Trade Science Inc. Matrix eigenvalues of jacobi iterative method for solving

Numerical Analysis Topic 4: Least Squares Curve Fitting

Section IV.6: The Master Method and Applications

10.5 Power Series. In this section, we are going to start talking about power series. A power series is a series of the form

The Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 11

UNIVERSITI KEBANGSAAN MALAYSIA PEPERIKSAAN AKHIR SEMESTER I SESI AKADEMIK 2007/2008 IJAZAH SARJANAMUDA DENGAN KEPUJIAN NOVEMBER 2007 MASA : 3 JAM

Chapter 7 Infinite Series

4 Linear Homogeneous Recurrence Relation 4-1 Fibonacci Rabbits. 组合数学 Combinatorics

Taylor Polynomials. The Tangent Line. (a, f (a)) and has the same slope as the curve y = f (x) at that point. It is the best

10.2 Series. , we get. which is called an infinite series ( or just a series) and is denoted, for short, by the symbol. i i n

Systems of second order ordinary differential equations

UNIT 2 SOLUTION OF ALGEBRAIC AND TRANSCENDENTAL EQUATIONS

Numerical Differentiation and Integration

APPLICATION OF THE CHEBYSHEV POLYNOMIALS TO APPROXIMATION AND CONSTRUCTION OF MAP PROJECTIONS

Math 1313 Final Exam Review

Introduction to Numerical Integration Part II

Chapter 3 Supplemental Text Material

Preliminary Examinations: Upper V Mathematics Paper 1

means the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever.

Homework 2 solutions

Mathematically, integration is just finding the area under a curve from one point to another. It is b

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

CS321. Numerical Analysis

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

Introduction to mathematical Statistics

CS321. Numerical Analysis

Log1 Contest Round 2 Theta Complex Numbers. 4 points each. 5 points each

Answer: First, I ll show how to find the terms analytically then I ll show how to use the TI to find them.

Chapter Unary Matrix Operations

6 Roots of Equations: Open Methods

Regression. By Jugal Kalita Based on Chapter 17 of Chapra and Canale, Numerical Methods for Engineers

SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS

St John s College. UPPER V Mathematics: Paper 1 Learning Outcome 1 and 2. Examiner: GE Marks: 150 Moderator: BT / SLS INSTRUCTIONS AND INFORMATION

Chapter Linear Regression

A Damped Guass-Newton Method for the Generalized Linear Complementarity Problem

Limit of a function:

A Technique for Constructing Odd-order Magic Squares Using Basic Latin Squares

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

SUM PROPERTIES FOR THE K-LUCAS NUMBERS WITH ARITHMETIC INDEXES

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

Pre-Calculus - Chapter 3 Sections Notes

Chapter 2 Intro to Math Techniques for Quantum Mechanics

Chapter 1. Introduction. Fundamental Concepts. Introduction. Historical background. Historical background. Fundamental Concepts

n 2 + 3n + 1 4n = n2 + 3n + 1 n n 2 = n + 1

I1 = I2 I1 = I2 + I3 I1 + I2 = I3 + I4 I 3

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

6. Chemical Potential and the Grand Partition Function

Mu Sequences/Series Solutions National Convention 2014

GAUSS ELIMINATION. Consider the following system of algebraic linear equations

Minimal DFA. minimal DFA for L starting from any other

Module 2: Introduction to Numerical Analysis

D. L. Bricker, 2002 Dept of Mechanical & Industrial Engineering The University of Iowa. CPL/XD 12/10/2003 page 1

More Regression Lecture Notes CE 311K - McKinney Introduction to Computer Methods Department of Civil Engineering The University of Texas at Austin

Solutions Manual for Polymer Science and Technology Third Edition

Mathematical models for computer systems behaviour

Patterns of Continued Fractions with a Positive Integer as a Gap

Chapter Simpson s 1/3 Rule of Integration. ( x)

Transcription:

Numercl Alyss for Egeers Germ Jord Uversty ITERATIVE METHODS FOR SOLVING SYSTEMS OF LINEAR ALGEBRAIC EQUATIONS Numercl soluto of lrge systems of ler lgerc equtos usg drect methods such s Mtr Iverse, Guss Elmto, Guss-Jord Elmto, d LU-Decomposto my clude lrge roud-off errors, especlly for ll-codtoed systems. Aother pproch for solvg such systems s usg tertve methods. The pproch s to strt wth tl soluto vector, d tertvely coverge to the true soluto. The tertos re the stopped whe pre-specfed tolerce s reched. Itertve methods re preferred over drect methods some cses such s whe the coeffcet mtr s sprse hs too my zeros). The tertve methods hve my dvtges such s My e more rpd. My e ecoomcl storge d clcultos. Self correctg. Reduce roud-off errors. JACOBI METHOD Cosder the followg system of equtos: stes.google.com/ste/zydmsoud/umercl

Numercl Alyss for Egeers Germ Jord Uversty The ove system c e rewrtte s ) Strtg wth tl guess vector, ew ppromto c e clculted y susttutg the tl guess the ove system s ) ) ) After ech terto, the reltve error of the lst ppromto c e determed s k k k % where k s the curret terto, k s the prevous terto, d. The reltve error s the compred wth commo tolerce for ll vrles, or specfc tolerce for ech vrle. Emple Solve the followg system of ler lgerc equtos usg Jco method d strtg wth the tl guess vector,, T. Use tolerce of.% for ll vrles. stes.google.com/ste/zydmsoud/umercl

Numercl Alyss for Egeers Germ Jord Uversty The ove system c e rewrtte s 6 7 5 5 5 6 7 5 Susttutg the tl guess vector the ove equtos, we c ot the frst ppromto of the system soluto 5 6 6.8 7 5 7.75 5 5. The error ths frst ppromto c e clculted s.8 % % %.8.75 % % %.75 %.. % % The secod ppromto of the system soluto c e oted y susttutg the frst ppromto ck the system of equtos s stes.google.com/ste/zydmsoud/umercl

Numercl Alyss for Egeers Germ Jord Uversty 5 6.75). 6.9 7 5.8).) 7.8 5.8.75) 5.856 The error ths frst ppromto c e clculted s % % %.9.8 %.5%.9.8.75 % 9.55%.8.856. % 76.5%.856 The secod ppromto s the used to determe the thrd d so o. The error s clculted every terto d compred wth the tolerce. Iterto results re show the tle elow. It #.8.75..9.8.856.5 9.55 76.5.86.5.8.5.5..5..8. 5.9. 5.995.99..5.79.59 6.997.9986.9955..7.659 7...9988.999.. stes.google.com/ste/zydmsoud/umercl

Numercl Alyss for Egeers Germ Jord Uversty stes.google.com/ste/zydmsoud/umercl 8....5.. 9...... GAUSS-SEIDEL METHOD The Guss-Sedel method s very smlr to the Jco method ecept tht ech ew ppromto s used mmedtely the followg clculto. Cosder the followg system of equtos: The ove system c e rewrtte s ) ) Strtg wth tl guess vector, ew ppromto c e clculted y susttutg the tl guess the ove system s

Numercl Alyss for Egeers Germ Jord Uversty stes.google.com/ste/zydmsoud/umercl 5 ) ) After ech terto, the reltve error of the lst ppromto c e determed s % k k k where k s the curret terto, k s the prevous terto, d. The reltve error s the compred wth commo tolerce for ll vrles, or specfc tolerce for ech vrle. Emple Solve the followg system of ler lgerc equtos usg Guss-Sedel method d strtg wth the tl guess vector T,,. Use tolerce of.% for ll vrles. 5 5 7 6 The ove system c e rewrtte s

Numercl Alyss for Egeers Germ Jord Uversty 5 6 7 5 Susttutg the tl guess vector the ove equtos, we c ot the frst ppromto of the system soluto 5 6 6.8 7 5.8) 7.9 5.8.9) 5.6 The error ths frst ppromto c e clculted s.8 % % %.8.9 % % %.9 %.6.6 % % The secod ppromto of the system soluto c e oted y susttutg the frst ppromto ck the system of equtos s 5 6.9).6 6.7 7 5.7).6) 7. 5.7.) 5.5 stes.google.com/ste/zydmsoud/umercl 6

Numercl Alyss for Egeers Germ Jord Uversty The error ths frst ppromto c e clculted s % % %.7.8 %.%.7..9 %.65%..5.6 %.6%.5 The secod ppromto s the used to determe the thrd d so o. The error s clculted every terto d compred wth the tolerce. Iterto results re show the tle elow. It #.8.9.6.7..5..65.6.999.9959.998.55.65.68....5.595.8 5....998.. We c oserve tht the Guss-Sedel method s much fster th the Jco method sce ewer d more ccurte ppromtos re lwys used mmedtely ll clcultos. stes.google.com/ste/zydmsoud/umercl 7

Numercl Alyss for Egeers Germ Jord Uversty stes.google.com/ste/zydmsoud/umercl 8 CONVERGENCE CRITERION FOR GAUSS-SEIDEL METHOD Recll tht for the Fed-Pot Iterto method ) g, the covergece crtero ws tht ) g The lgorthm of the Guss-Sedel method s As wth the Fed-Pot Iterto method, we c ssume tht g ) It c e prove tht the covergece of the Guss-Sedel method requres tht g g g g g, g,, g The

Numercl Alyss for Egeers Germ Jord Uversty or whch mes tht for gurteed covergece, the mgtude of the dgol coeffcet must e lrger th the sum of mgtudes of the off-dgol coeffcets. Emple The order of equtos of the pervous emple s ltered y terchgg equtos ) d ) of the system of equtos s show elow. Solve the system of ler lgerc equtos usg Guss-Sedel method d strtg wth the tl guess vector,, T. Use tolerce of.% for ll vrles. The ove system c e rewrtte s 6 7 5 5 stes.google.com/ste/zydmsoud/umercl 9

Numercl Alyss for Egeers Germ Jord Uversty 5 7 6 5 Susttutg the tl guess vector the ove equtos, we c ot the frst ppromto of the system soluto 5 7 6 5.5 6.5). 5.5).) 5 5.5 The error ths frst ppromto c e clculted s.5 % % %.5. % % %. % 5.5 5.5 % % The secod ppromto of the system soluto c e oted y susttutg the frst ppromto ck the system of equtos s 5 7 6 5 7.) 5.5) 5.5 5 6 5.5) 5.5) 68.8 5 5.5) 68.8) 5 78. stes.google.com/ste/zydmsoud/umercl 5

Numercl Alyss for Egeers Germ Jord Uversty The error ths frst ppromto c e clculted s % % % 5.5.5) % 95.% 5.5 68.8.) % 9.% 68.8 78. 5.5) % 9.95% 78. Crryg out more tertos, the soluto cotues to dverge wy from the system soluto. Ths s due to the fct tht the mgtude dgol coeffcet equtos ) d ) s smller th the sum of mgtudes of the off-dgol coeffcets of ech equto. stes.google.com/ste/zydmsoud/umercl 5