›
›

Bài giảng Toán cao cấp - Chương 4: Dạng toàn phương - Nguyễn Phương

pdf 27 trang huongle 19830
Download
Bạn đang xem 20 trang mẫu của tài liệu "Bài giảng Toán cao cấp - Chương 4: Dạng toàn phương - Nguyễn Phương", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pdfbai_giang_toan_cao_cap_chuong_4_dang_toan_phuong_nguyen_phuo.pdf

Nội dung text: Bài giảng Toán cao cấp - Chương 4: Dạng toàn phương - Nguyễn Phương

  1. Chương 4: DẠNG TOÀN PHƯƠNG Th.S NGUYỄN PHƯƠNG Khoa Giáo dục cơ bản Trường Đại học Ngân hàng TPHCM Blog: Email: nguyenphuong0122@gmail.com Yahoo: nguyenphuong1504 Ngày 28 tháng 10 năm 2013 1
  2. 1 Giá trị riêng - vectơ riêng Các định nghĩa Cách tìm vectơ riêng, giá trị riêng 2 Chéo hóa ma trận. Chéo hóa trực giao Định nghĩa chéo hóa Các bước chéo hóa ma trận vuông Chéo hóa trực giao ma trận đối xứng 3 Dạng toàn phương. Đưa DTP về dạng chính tắc Dạng toàn phương Đưa dạng toàn phương về dạng chính tắc Phương pháp Lagrange Phương pháp Jacobi 4 Dạng toàn phương xác định dấu. Định lý Sylvester Dạng toàn phương xác định dấu Định lý Sylvester 2
  3. Giá trị riêng - vectơ riêng Các định nghĩa Định nghĩa - Cho ma trận A Mn(R). Số thực λ được gọi là trị riêng của A nếu tồn tại n ∈ vector 0 < x R nếu A[x] = λ[x]. ∈ - Vector x , 0 thỏa A[x] = λ[x] được gọi là vector riêng của ma trận A ứng với trị riêng λ. ! 4 2 Ví dụ 1:Với A = x = (2; 1), ta được 1− 1 , ! ! ! ! 4 2 2 6 2 A[x] = = = 3 = 3[x] 1− 1 1 3 1 Vậy x = (2; 1) là vector riêng của A ứng với trị riêng λ = 3. Tính chất Nếu A có vectơ riêng x ứng với trị riêng λ thì kx, k , 0 cũng là vectơ riêng ứng với trị riêng λ. Nếu A có trị riêng λ thì λm là trị riêng của Am. m m Nếu A có trị riêng λ và A , 0 thì λ− là trị riêng của A− . | | 3
  4. Giá trị riêng - vectơ riêng Các định nghĩa Định nghĩa Cho ma trận vuông A = (aij) cấp n, ma trận đơn vị cấp n: I. - Ma trận đặc trưng của A là  a a a   11 λ 12 1n   −   a21 a22 λ . . . a2n  A λI =  −   . . .  −  . . . .    an1 an2 ann λ − - Đa thức đặc trưng của A là định thức của ma trận đặc trưng (là một đa thức λ), det(A λI). - Phương trình− đặc trưng của ma trận A là det(A λI) = 0. − ! 1 2 Ví dụ 2: Cho A = , ta có đa thức đặc trưng 3 4 1 λ 2 2 det(A λI) = − = λ 5λ 2 − 3 4 λ − − − 4
  5. Giá trị riêng - vectơ riêng Cách tìm vectơ riêng, giá trị riêng Cách tìm vectơ riêng, giá trị riêng - Bước 1: Giải phương trình đặc trưng:det(A λI) = 0. Các nghiệm tìm được là các giá trị riêng cần− tìm. - Bước 2: Giả sử λ0 là giá trị riêng. Giải hệ phương trình tuyến tính thuần nhất (A λ0I)X = 0. Nghiệm không− tầm thường của phương trình này là vectơ riêng cần tìm. ! 4 2 Ví dụ 3: Cho A = . Tìm giá trị riêng và vector riêng của A 1− 1 . Giải. Phương trình đặc trưng là 4 λ 2 2 det(A λI) = 0 − − = 0 λ 5λ + 6 = 0 − ⇔ 1 1 λ ⇔ − − Suy ra λ1 = 2 và λ2 = 3 là hai trị riêng của A. + Ứng với λ1 = 2: + Ứng với λ2 = 3: 5
  6. Giá trị riêng - vectơ riêng Cách tìm vectơ riêng, giá trị riêng  0 0 1    Ví dụ 4: Cho A =  0 1 0 . Tìm giá trị riêng và vector riêng của A   .  1 0 0  Giải. Phương trình đặc trưng: λ 0 1 − 2 0 1 λ 0 = 0 (1 λ)(λ 1) = 0 − ⇔ − − 1 0 λ − λ1 = 1, λ2 = 1 là hai trị riêng của A. ⇒ − + Với λ1 = 1, ta có: −  1 0 1   1 0 1  (     x1+ x3 = 0 A I =  0 2 0   0 1 0  λ1     x = 0 −  1 0 1  −→  0 0 0  ⇒ 2 x = α(1; 0; 1)(α , 0) là vetor riêng của A. ⇒ − 6
  7. Giá trị riêng - vectơ riêng Cách tìm vectơ riêng, giá trị riêng + Với λ2 = 1, ta có:  1 0 1   1 0 1      A I =  −0 0 0   −0 0 0  x x = 0 λ2     1 3 −  1 0 1  −→  0 0 0  ⇒ − − x = (α; β; α) là vetor riêng của A. ⇒ Không gian riêng Giả sử λ là giá trị riêng của ma trận A. Gọi tập hợp các vector riêng ứng với λ và vector không là E(λ).E(λ) được gọi là không gian riêng ứng với λ. Ví dụ 4: E( 1) = (1; 0; 1) ; dimE( 1) = 1 − − − E(1) = (1; 0; 1), (0; 1; 0) ; dimE(1) = 2 7
  8. Chéo hóa ma trận. Chéo hóa trực giao Định nghĩa chéo hóa Định nghĩa - Hai ma trận vuông cùng cấp A và B được gọi là đồng dạng nếu tồn tại ma 1 trận khả nghịch P thỏa B = P− AP. - Ma trận vuông cấp n được gọi là chéo hóa được nếu A đồng dạng với ma 1 trận đường chéo D, tức là, tồn tại ma trận P khả nghịch sao cho P− AP = D. Khi đó, ta nói ma trận P làm chéo hóa ma trận A. ! ! 1 0 1 0 Ví dụ: Hai ma trận A = và B = đồng dạng nhau vì có 6 1 −0 1 ! − 0 1 ma trận P = khả nghịch thỏa B = P 1AP 1 3 − .  0 0 0    Ví dụ: Ma trận A =  0 1 0  chéo hóa được, vì có ma trận    1 0 1   1 0 0   0 0 0      P =  0 1 0  thỏa P 1AP =  0 1 0    −    1 0 1   0 0 1  − 8
  9. Chéo hóa ma trận. Chéo hóa trực giao Các bước chéo hóa ma trận vuông Các bước chéo hóa ma trận vuông - Bước 1: Giải phương trình đặc trưng:det(A λI) = 0 để tìm trị riêng thực của A. − + Trường hợp A có trị riêng phức thì ta kết luận A không chéo hóa được. + Trường hợp A có n trị riêng thực thì A chéo hóa được, ta làm tiếp bước 2. + Trường hợp A có k trị riêng thực λi(i = 1, , k) với λi là nghiệm bội ni của phương trình đặc trưng. i) dim E(λi) = ni, i, ta kết luận A chéo hóa được, ta làm tiếp bước 2. ∀ ii) tồn tại dim E(λi) < ni, ta kết luận A không chéo hóa được. - Bước 2: Tìm cơ sở của các không gian riêng E(λi). - Bước 3: Lập ma trận P có các cột là các vector cơ sở của E(λi). 1 Khi đó, P− AP = D với D là ma trận chéo có các phần tử trên đường chéo chính lần lượt là λi (mỗi λi xuất hiện liên tiếp ni lần). 9
  10. Chéo hóa ma trận. Chéo hóa trực giao Các bước chéo hóa ma trận vuông  1 1 0    Ví dụ: Chéo hóa ma trận sau trên R:A =  2 −1 0     1− 2 3  1 λ 1 0 − − Giải. det(A λI) = 2 1 λ 0 = (3 λ)(1 + λ2) − − − − 1 2 3 λ − Phương trình det(A λI) = 0 không đủ 3 nghiệm thực nên A không chéo hóa − được trên R. ! 3 2 Ví dụ: Chéo hóa ma trận sau trên R:A = −2 1 − Giải. Phương trình đặc trưng det(A λI) = 0 (λ + 3)(λ 1) + 4 = 0 (λ + 1)2 = 0 λ = 1 (bội 2) − ⇔ − ⇔ ⇔ − ! ! ! ( 2 2 x1 0 2x1 + 2x2 = 0 (A λI)X = 0 − = − x1 = x2 − ⇔ 2 2 x2 0 ⇔ 2x1 + 2x2 = 0 ⇔ − − dim E( 1) = 1 < 2 A không chéo hóa được. ⇒ − ⇒ 10
  11. Chéo hóa ma trận. Chéo hóa trực giao Các bước chéo hóa ma trận vuông ! 1 0 Ví dụ: Chéo hóa ma trận A = . 6 1 − Giải. Phương trình đặc trưng " 1 λ 0 λ = 1 det(A λI) = − = 0 − A chéo hóa được. − 6 1 λ ⇔ λ = 1 ⇒ − − ! ! 2 0 1 0 λ = 1 :(A λI) = x = (0; 1) − − 6 0 → 0 0 ⇒ ! ! 0 0 3 1 λ = 1 :(A λI) = − x = (1; 3) − 6 2 → 0 0 ⇒ − ! ! 0 1 1 0 suy ra P = và D = 1 3 −0 1 ! 1 0 Vậy P 1AP = − −0 1 11
  12. Chéo hóa ma trận. Chéo hóa trực giao Các bước chéo hóa ma trận vuông  4 2 1    Ví dụ: Chéo hóa ma trận A =  6 4− 3     −6 −6 5  − − Giải. Phương trình đặc trưng 4 λ 2 1 " − − 2 λ = 1 det(A λI) = 6 4 λ 3 = 0 (λ 1)(λ 2) = 0 − − − − ⇔ − − ⇔ λ = 2 6 6 5 λ − − − + Với λ = 1 (nghiệm đơn), ta có:      3 2 1   3 0 1  (  −    3x1 + x3 = 0 A λI =  6 5 3   0 1 1  −  − −  →  −  ⇒ x2 x3 = 0 6 6 4 0 0 0 − − − x = α(1; 3; 3) dim E(1) = 1 ⇒ − − ⇒ 12
  13. Chéo hóa ma trận. Chéo hóa trực giao Các bước chéo hóa ma trận vuông + Với λ = 1 (nghiệm bội 2), ta có:  2 2 1   2 2 1      A I =  6 6− 3   0 0− 0  2x + 2x x = 0 λ     1 2 3 −  −6 −6 3  →  0 0 0  ⇒ − − − x = α(1; 0; 2) + β(0; 1; 2) dim E(2) = 2 ⇒ ⇒  1 0 0   1 1 0      Vậy A chéo hóa được và P 1AP =  0 2 0  với P =  3 0 1 . −      0 0 2   −3 2 2  − 13
  14. Chéo hóa ma trận. Chéo hóa trực giao Chéo hóa trực giao ma trận đối xứng Định nghĩa Ma trận thực A vuống cấp n được gọi là trực giao nếu ATA = I, tức là 1 T A− = A .  1 1       − √2 √2  Ví dụ:   là ma trận trực giao.  1 1    √2 √2 Định lý Nếu A là ma trận đối xứng thì các vector riêng thuộc các không gian riêng khác nhau là trực giao. Định nghĩa 1 Nếu tồn tại ma trận trực giao P sao cho P− AP là ma trận chéo thì A là ma trận chéo hóa trực giao được và P là ma trận làm chéo hóa trực giao A. 14
  15. Chéo hóa ma trận. Chéo hóa trực giao Chéo hóa trực giao ma trận đối xứng Các bước chéo hóa ma trận đối xứng - Bước 1:Tìm các trị riêng. - Bước 2: Tìm các vector cơ sở của các không gian riêng. - Bước 3: Chuẩn hóa các vector cơ sở này. - Bước 4: Lập ma trận chéo hóa trực giao của A.  9 2 2    Ví dụ: Chéo hóa ma trận đối xứng A =  2 0 2 .    2 2 0  15
  16. Dạng toàn phương. Đưa DTP về dạng chính tắc Dạng toàn phương Định nghĩa - Dạng toàn phương của n biến x1, x2, , xn (hay là dạng toàn phương trong Kn) là hàm Q từ Kn đến K cho bởi biểu thức Xn Xn Q(x) = aijxixj với aij = aji i=1 j=1  a a a   11 12 1n       a21 a22 a2n  - Nếu ta đặt A = aij =   Mn(K) và n n  . . .  ×  . . . .  ∈   an1 an2 ann X = (x1, x2, , xn)thì dạng toàn phương được viết ở dạng ma trận Q(X) = [x]TA[x] A được gọi là ma trận của dạng toàn phương, A là một ma trận đối xứng. 16
  17. Dạng toàn phương. Đưa DTP về dạng chính tắc Dạng toàn phương ! 1 1 Ví dụ: Tìm dạng toàn phương q(x), biết ma trận của q là A = . 1− 2 − Giải. Ta có ! ! T 1 1 x1 2 2 q(x) = [x] A[x] = (x1 x2) − = x1 + 2x2 2x1x2 1 2 x2 − − Ví dụ: Tìm ma trận của dạng toàn phương q : R3 R sau → 2 2 2 q(x) = 2x + 3x x 4x1x2 + 6x2x3 1 2 − 3 − Giải.  2 2 0    A =  2− 3 3     −0 3 1  − 17
  18. Dạng toàn phương. Đưa DTP về dạng chính tắc Dạng toàn phương Định nghĩa Trong Rn, dạng toàn phương n 2 2 2 X 2 Q(X) = a11x1 + a22x2 + + anxn = aiixi ··· i=1 được gọi là dạng toàn phương chính tắc hay gọi tắt là dạng chính tắc.   Ma trận của dạng chính tắc là A = diag a11 a22 ann . ! 1 0 Ví dụ: Trong R2, cho dạng chính tắc có ma trận A = . 0 2 − Khi đó, biểu thức của q là ! ! T 1 0 x1 2 2 q(x) = [x] A[x] = (x1 x2) = x1 2x2 0 2 x2 − − 18
  19. Dạng toàn phương. Đưa DTP về dạng chính tắc Đưa dạng toàn phương về dạng chính tắc Phương pháp Lagrange: - Phương pháp Lagrange nhóm trực tiếp lần lượt theo từng biến xi dưới dạng tổng bình phương. Sau đó, ta dựa vào các tổng bình phương để đổi biến. Trường hợp 1: (Q(x) có hệ số aii , 0) Bước 1: Giả sử a11 , 0, ta tách tất cả các số hạng chứa x1 trong Q(x) và • thêm hoặc bớt để có dạng 2 Q(x) = a11(x1 + α12x2 + + α1nxn) + Q1(x2; ; xn) ··· Đổi biến y1 = x1 + α12x2 + + α1nxn , yi = xi (i = 2, , n) ··· , tức là x1 = y1 α12x2 α1nxn , xi = yi (i = 2, , n) − − · · · − 2 Khi đó, Q(x) trở thành Q(x) = a11(y1) + Q1(y2; ; yn) Bước 2: Tiếp tục làm như bước 1 cho Q1(y2; ; yn). Sau vài bước thì dạng toàn• phương Q(x) sẽ có dạng chính tắc. Trường hợp 2: (Q(x) có tất cả các hệ số aii = 0) Giả sử a12 , 0, ta thực hiện đổi biến: x1 = y1 + y2, x2 = y1 y2, xi = yi (i = 2, , n) − 2 2 2 Khi đó, dạng toàn phương Q(x) = 2a12y 2a12y + có hệ số của y , 0. 1 − 2 ··· 1 Ta thực hiện theo trường hợp 1. 19
  20. Dạng toàn phương. Đưa DTP về dạng chính tắc Đưa dạng toàn phương về dạng chính tắc Ví dụ: Đưa dạng toàn phương về dạng chính tắc: 2 2 Q(x) = x + 4x + 2x1x2 + 4x1x3 − 2 3 Bước 1. Biến đổi 2 2 2 2 2 2 2 Q(x) = (x 2x1x2 + x ) + x + 4x + 4x1x3 = (x1 x2) + x + 4x + 4x1x3 − 1 − 2 1 3 − − 1 3    y1 = x1  x1 = y1   Đổi biến  y2 = x1 x2  x2 = y1 y2  − ⇔  −  y3 = x3  x3 = y3 2 2 2 2 2 Khi đó, Q(x) = y2 + y1 + 4y3 + 4y1y3 với Q1(y) = y1 + 4y3 + 4y1y3. − 2 2 2 Bước 2. Biến đổi Q1(y) = y + 4y + 4y1y3 = (y1 + 2y3)  1 3   z1 = y1 +2y3  y1 = z1 2z3   − Đổi biến  z2 = y2  y2 = z2  ⇔   z3 = y3  y3 = z3   x1 = z1 2z3 2 2  − Vậy Q(x) = z1 z2 với  x2 = z1 z2 2z3 −  − −  x3 = z3 20
  21. Dạng toàn phương. Đưa DTP về dạng chính tắc Đưa dạng toàn phương về dạng chính tắc Phương pháp Jacobi: - Phương pháp Jacobi áp dụng cho dạng toàn phương có ma trận A thỏa điều kiện a11 a12 D1 = a11 , 0, D2 = , 0, , Dn = det(A) , 0. a21 a22 Khi đó, biến đổi theo công thức   x1 = y1 +b21y2 +b31y3 + +bn1yn  ···  x2 = y2 +b32y3 + +bn2yn  ···  .  .   xn = yn i+j Di 1,j với bij = ( 1) − (j < i);Di 1 j là định thức được tạo thành từ i 1 dòng Di 1 , − − − − đầu và i 1 cột đầu (sau khi đã bỏ cột thứ j)của ma trận A. Khi đó, ma trận −  1 b b   21 n1     0 1 bn2  đổi biến là P =    . . .   . . . .    0 0 1 2 D2 2 D3 2 Dn 2 và dạng chính tắc là : Q = D1y1 + y2 + y3 + + yn. D1 D2 ··· Dn 1 21 −
  22. Dạng toàn phương. Đưa DTP về dạng chính tắc Đưa dạng toàn phương về dạng chính tắc Ví dụ: Đưa dạng toàn phương về dạng chính tắc 2 2 2 Q(x) = 2x1 + 3x1x2 + 4x1x3 + x2 + x3 Giải:  2 3 2   2  Ma trận của Q: A =  3 1 0   2   2 0 1  3 2 2 1 17 Ta có: D1 = 2; D2 = 3 = 4 ; D3 = det(A) = 4 2 1 − − Ta tính: 3 2 2 3 D 3 D 1 0 2+1 1,1 − 2 3+1 2,1 b21 = ( 1) = = ; b31 = ( 1) = 1 = 8 − D1 2 −4 − D2 − 4 2 2 3 0 3+2 D2,2 2 b32 = ( 1) = 1 = 12 − D2 − − 4  3  x1 = y1 4 y2 + 8y3  − 2 1 2 2 Vậy đổi biến  x2 = y2 12y3 ta được Q = 2y1 8 y2 + 17y3  − −  x3 = y3 22
  23. Dạng toàn phương xác định dấu. Định lý Sylvester Dạng toàn phương xác định dấu Định lý (Luật quán tính) Số hệ số dương và hệ số âm trong dạng chính tắc là những đại lượng bất biến, không phụ thuộc vào phép biến đổi tuyến tính không suy biến đưa dạng toàn phương về dạng chính tắc. Định nghĩa Dạng toàn phương Q(x) được gọi là: - Xác định dương nếu Q(x) > 0, x Rn và x , 0. ∀ ∈n - Xác định âm nếu Q(x) < 0, x R và x , 0. - Không xác định dấu nếu nó∀ nhận∈ cả giá trị âm và giá trị dương. 23
  24. Dạng toàn phương xác định dấu. Định lý Sylvester Dạng toàn phương xác định dấu Định lý Dạng toàn phương Q(x) của Rn xác định dương khi và chỉ khi tất cả n hệ số trong dạng chính tắc của nó đều dương. Hệ quả Dạng toàn phương Q(x) của Rn xác định âm khi và chỉ khi tất cả n hệ số trong dạng chính tắc của nó đều âm. Định lý Dạng toàn phương Q(x) của Rn xác định dương khi và chỉ khi ma trận của nó có tất cả các giá trị riêng đều dương. Hệ quả -Q(x) xác định âm khi và chỉ khi ma trận của nó có tất cả các giá trị đều âm. -Q(x) không xác định dấu khi và chỉ khi ma trận của nó giá trị riêng dương và giá trị riêng âm. 24
  25. Dạng toàn phương xác định dấu. Định lý Sylvester Định lý Sylvester Định lý (Định lý Sylvester) - Dạng toàn phương Q(x) xác định dương khi và chỉ khi tất cả các định thức con chính của ma trận của nó đều dương, tức là D1 > 0, D2 > 0, , Dn > 0. - Dạng toàn phương Q(x) xác định âm khi và chỉ khi các định thức con chính k cấp chẵn dương, cấp lẻ âm, tức là ( 1) Dk > 0, k. − ∀ Ví dụ: Xét tính xác định dấu của dạng toàn phương trong R3: 2 2 2 Q(x) = 2x 4x 3x + 4x1x2. − 1 − 2 − 3  2 2 0    Giải. Ma trận của Q là A =  −2 4 0     0− 0 3  − Ta có các định thức con chính: 2 2 D1 = 2 0, D3 = A = 12 < 0 − 2 4 | | − − Vậy Q xác định âm. 25
  26. Dạng toàn phương xác định dấu. Định lý Sylvester Định lý Sylvester Ví dụ: Xét tính xác định dấu của dạng toàn phương trong R3: 2 2 2 Q(x) = 7x + 2x x + 6x1x3. 1 2 − 3  7 0 3    Giải. Ma trận của Q là A =  0 2 0     3 0 1  − Ta có các định thức con chính: 7 0 D1 = 7 > 0, D2 = = 14 > 0, D3 = A = 32 < 0 0 2 | | − Vậy Q không xác định dấu. 26
  27. Dạng toàn phương xác định dấu. Định lý Sylvester Định lý Sylvester Ví dụ: Xét tính xác định dấu của dạng toàn phương trong R3: 2 2 2 Q(x) = x + 4x + mx 2x1x2 + 8x1x3 + 4x2x3. 1 2 3 −  1 1 4    Giải. Ma trận của Q là A =  1− 4 2     −4 2 m  Dạng toàn phương xác định dương khi và chỉ khi tất cả các định thức con chính xác định dương. 1 1 D1 = 1 > 0, D2 = − = 3 > 0, D3 = A = 3(m 28) > 0 m > 28 1 4 | | − ⇔ − 27