Chapter Gauss-Seidel Method

Transcription

1 Chpter Guss-Sedel Method After redg ths hpter, you should be ble to:. solve set of equtos usg the Guss-Sedel method,. reogze the dvtges d ptflls of the Guss-Sedel method, d. determe uder wht odtos the Guss-Sedel method lwys overges. Why do we eed other method to solve set of smulteous ler equtos? I ert ses, suh s whe system of equtos s lrge, tertve methods of solvg equtos re more dvtgeous. Elmto methods, suh s Guss elmto, re proe to lrge roud-off errors for lrge set of equtos. Itertve methods, suh s the Guss-Sedel method, gve the user otrol of the roud-off error. Also, f the physs of the problem re well kow, tl guesses eeded tertve methods be mde more udously ledg to fster overgee. Wht s the lgorthm for the Guss-Sedel method? Gve geerl set of equtos d ukows, we hve If the dgol elemets re o-zero, eh equto s rewrtte for the orrespodg ukow, tht s, the frst equto s rewrtte wth o the left hd sde, the seod equto s rewrtte wth o the left hd sde d so o s follows

2 Chpter 04.08,,,,,, These equtos be rewrtte summto form s...,, Hee for y row,.,,,, Now to fd s, oe ssumes tl guess for the s d the uses the rewrtte equtos to lulte the ew estmtes. Remember, oe lwys uses the most reet estmtes to lulte the et estmtes,. At the ed of eh terto, oe lultes the bsolute reltve ppromte error for eh s 00 ew old ew where ew s the reetly obted vlue of, d old s the prevous vlue of.

3 Guss-Sedel Method Whe the bsolute reltve ppromte error for eh s less th the pre-spefed tolere, the tertos re stopped. Emple The upwrd veloty of roket s gve t three dfferet tmes the followg tble Tble Veloty vs. tme dt. Tme, t (s) Veloty, v (m/s) The veloty dt s ppromted by polyoml s v t t t, 5 t Fd the vlues of,, d usg the Guss-Sedel method. Assume tl guess of the soluto s 5 d odut two tertos. Soluto The polyoml s gog through three dt pots, v, t, v, d t v bove tble t 5, v 06.8 t 8, v 77. t, v 79. Requrg tht v t t t v t t v t t v t t vt vt vt Substtutg the dt, v, t, v, d t v or, t where from the, psses through the three dt pots gves t gves

4 Chpter The oeffets,, d for the bove epresso re gve by Rewrtg the equtos gves Iterto # Gve the tl guess of the soluto vetor s 5 we get () (5) The bsolute reltve ppromte error for eh the s % % % At the ed of the frst terto, the estmte of the soluto vetor s

5 Guss-Sedel Method d the mmum bsolute reltve ppromte error s 5.47%. Iterto # The estmte of the soluto vetor t the ed of Iterto # s Now we get ( 55.6) ( 55.6) = The bsolute reltve ppromte error for eh the s % % % At the ed of the seod terto the estmte of the soluto vetor s d the mmum bsolute reltve ppromte error s %. Codutg more tertos gves the followg vlues for the soluto vetor d the orrespodg bsolute reltve ppromte errors.

6 Chpter Iterto % % % As see the bove tble, the soluto estmtes re ot overgg to the true soluto of The bove system of equtos does ot seem to overge. Why? Well, ptfll of most tertve methods s tht they my or my ot overge. However, the soluto to ert lsses of systems of smulteous equtos does lwys overge usg the Guss-Sedel method. Ths lss of system of equtos s where the oeffet mtr [A] [ A][ X ] [ C] s dgolly domt, tht s for ll for t lest oe If system of equtos hs oeffet mtr tht s ot dgolly domt, t my or my ot overge. Fortutely, my physl systems tht result smulteous ler equtos hve dgolly domt oeffet mtr, whh the ssures overgee for tertve methods suh s the Guss-Sedel method of solvg smulteous ler equtos. Emple Fd the soluto to the followg system of equtos usg the Guss-Sedel method Use 0 s the tl guess d odut two tertos. Soluto The oeffet mtr

7 Guss-Sedel Method A 5 7 s dgolly domt s d the equlty s strtly greter th for t lest oe row. Hee, the soluto should overge usg the Guss-Sedel method. Rewrtg the equtos, we get Assumg tl guess of 0 Iterto # The bsolute reltve ppromte error t the ed of the frst terto s % % %

8 Chpter The mmum bsolute reltve ppromte error s 00.00% Iterto # At the ed of seod terto, the bsolute reltve ppromte error s % % % The mmum bsolute reltve ppromte error s 40.6%. Ths s greter th the vlue of 00.00% we obted the frst terto. Is the soluto dvergg? No, s you odut more tertos, the soluto overges s follows. Iterto % Ths s lose to the et soluto vetor of 4 Emple Gve the system of equtos 7 76 % %

9 Guss-Sedel Method fd the soluto usg the Guss-Sedel method. Use 0 s the tl guess. Soluto Rewrtg the equtos, we get Assumg tl guess of 0 the et s tertve vlues re gve the tble below. Iterto % % % You see tht ths soluto s ot overgg d the oeffet mtr s ot dgolly domt. The oeffet mtr 7 A 5 5 s ot dgolly domt s 7 0 Hee, the Guss-Sedel method my or my ot overge. However, t s the sme set of equtos s the prevous emple d tht overged. The oly dfferee s tht we ehged frst d the thrd equto wth eh other d tht mde the oeffet mtr ot dgolly domt.

10 Chpter Therefore, t s possble tht system of equtos be mde dgolly domt f oe ehges the equtos wth eh other. However, t s ot possble for ll ses. For emple, the followg set of equtos ot be rewrtte to mke the oeffet mtr dgolly domt. Key Terms: Guss-Sedel method Covergee of Guss-Sedel method Dgolly domt mtr