Phương pháp tính
Chương 4: Giải Hệ phương trình tuyến tính1. Đặt vấn đề: Trong chương này ta xét việc giải hệ phương trình đại số n phương trình, n ẩn số.=+++=+++=++++++12211122222121111212111 .nnnnnnnnnnnnnaxaxaxaaxaxaxaaxaxaxatrong đó aij ( i,j = n,1) gọi là hệ số của hệ phương trình ain+1( i = n,1) gọi là vế phải của hệ phương trình xi ( i = n,1) là các ẩn số phải tìmA = nnnnnnnaaaaaaaaaaaa .32122322211131211gọi là ma trận hệ số của hệ phương trình b = +++11211 nnnnaaa gọi là vectơ vế phải của hệ phương trình x = nxxx 21 gọi là vectơ ẩn số của hệ phương trình Hệ phương trình trên được viết gọn dưới dạng: Ax = bNếu ma trận A không suy biến, nghĩa là det(A) ≠ 0 thì hệ phương trình có nghiệm duy nhất.Thật vậy, vì det(A) ≠ o nên tồn tại A-1, ta suy ra A-1.A.x = A-1.b ⇒ x = A-1bPhương pháp Cramer, giải đúng hệ phương trình bằng công thức:xi = ),1( nii=∆∆ trong đó ∆ = det(A)∆i là định thức cấp n thu được bằng cách thay cột thứ i của ∆ bằng cột vế phải b của hệ phương trình.Khi n lớn, phương pháp Cramer với khối lượng tính toán tăng rất nhanh. Vì vậy, người ta phải xây dựng những phương pháp giải hệ phương trình sao cho khi n lớn, khối lượng tính toán không quá lớn, có thể thực hiện được.2. Phương pháp trực tiếp: ( phương pháp Gauss, phương pháp khử )Phương pháp Gauss là 1 phương pháp được dùng phổ biến để giải hệ phương trình, đặc biệt người ta dùng phương pháp này khi ma trận hệ số A “đầy” (nghĩa là những phần tử = 0 ít và cấp của ma trận không quá lớn)2.1. Nội dung phương pháp: Để đơn giản việc trình bày, xét hệ thống 4 phương trình, 4 ẩn số:=+++=+++=+++=+++)0(454)0(443)0(432)0(421)0(41)0(354)0(343)0(332)0(321)0(31)0(254)0(243)0(232)0(221)0(21)0(154)0(143)0(132)0(121)0(11axaxaxaxaaxaxaxaxaaxaxaxaxaaxaxaxaxa (I)Nội dung của phương pháp Gauss là khử dần các ẩn số để đưa hệ (I) về hệ tam giác tương đương ( ma trận hệ số của hệ là ma trận tam giác trên)==+=++=+++)4(454)3(354)3(343)2(254)2(243)2(232)1(154)1(143)1(132)1(121axaxaxaxaxaxaxaxaxax (II)Sau đó giải (II) từ dưới lên(I) → (II) : Quá trình thuậnGiải (II) gọi là quá trình ngượca) Quá trình thuận:* Khử x1: Giả sử )0(11a ≠ 0 ()0(11a gọi là trụ thứ nhất)Chia phương trình đầu của (I) cho )0(11a, ta được:)1(154)1(143)1(132)1(121axaxaxax =+++ với )0(11)0(1)1(1aaajj= (j= 2,3,4,5))1(1)0(1)0()1(jiijijaaaa −= (i=2,3,4 ; j= 2,3,4,5)*Khử x2: chú ý các công thức:)1(22)1(2)2(2aaajj=( j= 3,4,5))2(2)1(2)1()2(jiijijaaaa −=(i= 3,4; j = 3,4,5)*Khử x3: chú ý các công thức:)2(33)2(3)3(3aaajj=( j= 4,5))3(3)2(43)2(4)3(4 jjjaaaa −=( j = 3,4,5)*Khử x4: với )3(44)3(45)4(45aaa =b) Quá trình ngược:Giải hệ phương trình (II)x4 = )4(45ax3 = 4)3(34)3(35xaa −x2 = 4)2(243)2(23)2(25xaxaa −−x1 = 4)1(143)1(132)1(12)1(15xaxaxaa −−− 2.2 Sơ đồ tính và kiểm tra quá trình tính:Sơ đồ tính :x1x2x3x4Số hạng tự doΣQuá trìnhQuá trình thuậnx4x3x2x14x3x2x1xQuá trình ngượcĐể kiểm tra người ta dùng cột “tổng kiểm tra” (Σ)∑==51)0()0(6jijiaa (i= 1,2,3,4)nó như một vế phải mới của hệ phương trình :)0(651)0(ijjijaxa =∑= ( i= 1,2,3,4)rõ ràng : 1+=jjxxĐối với mỗi phần tử của cột Σ ta đều thực hiện những phép tính giống như đối vớI những phần tử của những cột ở bên trái cột Σ và nằm trong cùng một hàng vớI phần tử của cột Σ, nên nếu không có sai số tính toán thì những phần tử của cột Σ phải bằng tổng những phần tử tương ứng của những cột ở bên trái cột Σ. hiện tượng này được dùng để kiểm tra quá trình thuận. Quá trình ngược được kiểm tra bằng hệ thức 1+=jjxx2.3 Phương pháp Gauss có tìm trụ lớn nhất:Phương pháp Gauss sẽ không thực hiện được nếu một trong các trụ bằng 0, hoặc tuy khác 0 nhưng về trị tuyệt đối rất nhỏ so với các phần tử còn lại trong cùng hàng thì khi chia các phần tử ấy cho phần tử trụ, sai số làm tròn sẽ lớn, do đó giảm nhiều độ chính xác của nghiệm tìm được.Để khắc phục những hạn chế vừa nêu, người ta thường dùng phương pháp Gauss có tìm trụ lớn nhất. Chẳng hạn, khi khử x1, người ta chọn số lớn nhất về trị tuyệt đối trong các số hệ số của x1 trong các phương trình làm trụ thứ nhất và hoán vị hàng thứ nhất với hàng có chứa trụ vừa tìm được. Cứ thế tương tự khi khử x2, x3,…2.4 Chuẩn của ma trận và chuẩn của vectơ:Ba chuẩn thông dụng của ma trận:- Chuẩn cột :∑=iijjamaxA1- Chuẩn Euclide:∑=jiijaA,22- Chuẩn hàng:∑=∞jijiaA maxVectơ là ma trận chỉ có một cột, ta có 3 chuẩn:-∑==+++=niinxxxxx1211 .-∑==niixx122-iixx max=∞3. Phương pháp lặp đơn:Những phương pháp lặp trong đó có phương pháp lặp đơn là những phương pháp giải gần đúng hệ phương trình đại số tuyến tính. Chúng chỉ cho nghiệm gần đúng của hệ phương trình, dù các phép tính trung gian tính đúng hoàn toàn, nhưng ta có thể tính nghiệm gần đúng ấy với độ chính xác bất kỳ. Người ta dùng những phương pháp lặp khi ma trận hệ số A “thưa” (số phần tử ≠ 0 ít, cấp của ma trận lớn)3.1 Nội dung phương pháp:Xét hệ phương trình Ax = b, ta hệ phương trình về dạng tương:x = β + αxtrong đó : ==nnnnnnnααααααααααββββ .2122221112112Sau đó tự cho 1 vectơ gọi là vectơ xấp xỉ đầu, ký hiệu x(0) ( thường chọn x(0) = β), sau đó tính dần x(k) theo công thức:x(k+1) = β + αx(k) (k= 0, 1, 2, …)Vectơ x(k) gọi là vectơ lặp thứ kNếu dãy vectơ lặp x(0), x(1),…, x(k) có giới hạn thì giới hạn đó chính là nghiệm đúng x* của hệ phương trình.3.2 Sự hội tụ của phương pháp: Quá trình lặp đơn hội tụ nếu:1<pα3.3 Đánh giá sai số của nghiệm gần đúng:Ta có 2 công thức:pkkpppkxxxx)1()(*)(1−−−≤−ααvà ppkppkxxxx)0()1(*)(1)(−−≤−ααThí dụ: Dùng phương pháp lặp đơn tìm nghiệm gần đúng của hệ phương trình:=+−=−+=−+20408,004,0915,0309,0808,024,04321321321xxxxxxxxx . thức:)2(33)2(3)3(3aaajj=( j= 4, 5))3(3)2 (43 )2 (4) 3 (4 jjjaaaa −=( j = 3 ,4, 5)*Khử x4: với )3 (44 )3 (45 )4( 45aaa =b) Quá trình ngược:Giải hệ phương trình (II)x4 = )4( 45ax3 = 4) 3( 34) 3(35xaa. dung phương pháp: Để đơn giản việc trình bày, xét hệ thống 4 phương trình, 4 ẩn số:=+++=+++=+++=+++)0 (45 4)0 (44 3)0 (43 2)0 (42 1)0 (41 )0(3 54) 0( 343 )0(332)0(321)0(31)0(2 54) 0( 243 )0(232)0(221)0(21)0(1 54) 0( 143 )0(132)0(121)0(11axaxaxaxaaxaxaxaxaaxaxaxaxaaxaxaxaxa