Xem mẫu

  1. CHƢƠNG VI : PHƢƠNG TRÌNH SAI PHÂN I. Các khái niệm cơ bản 1. Hàm số đối số nguyên Hàm có tập xác định thuộc Z gọi là hàm số có đối số nguyên. Ký hiệu y = f(n). f(n) = n2 + n – 1 Ví dụ: f(n) = n3 + 1 f(n) = sina (a là hằng số) 2. Định nghĩa sai phân: Sai phân của hàm số Un là chênh lệch giá trị của hàm số tại hai giá trị kế tiếp nhau. Ký hiệu: ΔUn = Un +1 - Un Sai phân cấp m của hàm số Un là sai phân của sai phân cấp m-1 của hàm số đó : ΔmUn = Δ(Δm-1Un )= Δm-1Un +1 - Δm-1Un Chẳng hạn sai phân cấp 2 được tính : Δ2Un = Δ(ΔUn )= ΔUn +1 – ΔUn= (Un +2 - Un+1 )- (Un +1 – Un ) = Un +2 -2 Un +1 + Un Tương tự ta có thể biểu diễn ΔmUn qua Un , Un+1,..., Un+m I. Phƣơng trình sai phân Định nghĩa : là PT với hàm số phải tìm là 1 hàm đối số rời rạc f (n) = Un có mặt dưới dạng sai phân các cấp. PT sai phân cấp m có dạng tổng quát : G(n, Un, ΔUn, Δ2Un,..., ΔmUn) = 0 Hay có thể viết dưới dạng : F(n, Un, Un+1,..., Un+m) = 0 Nghiệm của PT sai phân là hàm số đối số rời rạc Un =f(n) mà khi thay Un = f(n), Un+1 =f(n+1),..., Un+m =f(n+m) ta được một đồng nhất thức trên tập hợp các số nguyên n 0. Nghiệm tổng quát của một PT sai phân cấp n có dạng : Un =f(n, C1, C2,...,Cn) trong đó C1, C2,...,Cn là các hằng số bất kì, khi gán cho mỗi kí tự C1, C2,...,Cn một số xác định ta được một nghiệm riêng của PT. PT sai phân Ôtônôm là PT có dạng Un+m = f(Un, Un+1,..., Un+m-1) 1
  2. II. Phƣơng trình sai phân tuyến tính 1. Phương trình sai phân tuyến tính cấp 1 Định nghĩa: Là phương trình có dạng: anUn+1 + bnUn = fn (1) Trong đó an, bn, fn là các hàm đối số nguyên. Un và Un+1 là hai giá trị kề nhau của hàm Un đối số nguyên cần tìm. Nếu an và bn là các hằng số thì ta có phương trình sai phân hệ số hằng. Phương trình anUn+1 + bnUn = 0 (2) gọi là phương trình thuần nhất tương ứng của (1). Ví dụ: Một khách hàng có số tiền là A đồng, đem gửi tiết kiệm, lãi xuất mỗi tháng là 1%. Lập mô hình về tình hình tiền vốn của khách hàng. 1 Ta có un+1 = un + 100 un = 1,01.un un+1 – 1,01.un = 0, u0 = A 2. Phương trình sai phân cấp cao a. Phương trình sai phân cấp 2 Dạng : an.un+2 + bn.un+1 + cn.un = fn Nếu an, bn và cn là các hằng số thì ta có phương trình sai phân hệ số hằng. Nếu fn = 0 thì ta có phương trình thuần nhất liên kết an.un+2 + bn.un+1 + cn.un = 0 Nếu U*n là một nghiệm của PT sai phân tuyến tính không thuần nhất và U1n, U2n là 2 nghiệm độc lập tuyến tính của PT thuần nhất liên kết thì nghiệm tổng quát của PT là : U = U*n+ C1U1n + C2 U2n Ví dụ: Ngày 01/ 01/ 1202, Giáo hoàng La Mã cho Fibonacci một bài toán như sau: “Hôm nay, người ta tặng tôi một cặp thỏ. Biết thỏ hai tháng tuổi bắt đầu đẻ và sau đó mỗi tháng đẻ một lứa, mỗi lứa là một cặp thỏ. Hết năm, tôi có bao nhiêu cặp thỏ ?” Giải: Gọi Fn là số cặp thỏ có được ở tháng thứ n. Tháng trước có Fn-1 cặp, trong đó chỉ có số thỏ tháng trước nữa là đẻ Fn = Fn-1 + Fn-2 với F1 = 1, F2 = 1. b. Phương trình sai phân cấp k Là phương trình có dạng: ak.Un+k + ak-1.Un+k-1 + … + a0.Un = fn 2
  3. III. Phƣơng trình sai phân tuyến tính cấp 1 hệ số hằng 1. Phương trình sai phân tuyến tính thuần nhất Nghiệm tổng quát : Un = C(- p) n Dạng Un+1 + pUn = 0 Un+1 = - pUn Ví dụ: Năm 1990 dân số Hà Nội là 1,6 triệu người, tốc độ tăng dân số là 1% một năm. Hỏi dân số Hà Nội năm 2050 là bao nhiêu? Giải: Gọi un là dân số Hà Nội năm thứ n + 1990 1 Ta có un+1 = un + 100 un = 1,01.un un = u0.(1,01)n. Có u0 = 1,6 triệu u60 = 1,6.(1,01)60 2.91 triệu. 2. Phương trình sai phân tuyến tính không thuần nhất Dạng Un+1 + pUn = q (1) với q 0. PT thuần nhất liên kết Un+1 + pUn = 0 (2). Định lý : Nếu U*n là một nghiệm của PT sai phân tuyến tính không thuần nhất (1) và U1n là một nghiệm của PT thuần nhất liên kết (2) thì U1n+ U*n là nghiệm của PT (1). Nghiệm tổng quát của (1) dạng Un= U*n + C(- p) n Ta tìm nghiệm riêng của (1) : q +) Nếu p -1 nghiệm riêng là U*n = 1p U*n +) Nếu p = -1 nghiệm riêng là = qn. IV. Phƣơng trình sai phân tuyến tính cấp 2 hệ số hằng 1. Phương trình sai phân tuyến tính thuần nhất : Xét phương trình: Un+2 + pUn+1 + qUn = 0 (3) Bổ đề 1: Nếu xn, yn là nghiệm của (3) thì A.xn + B.yn (A, B : const) cũng là nghiệm của (3). Chứng minh: Ta có: (A.xn+2 + B.yn+2) + p.(A.xn+1 + B.yn+1) + q.(A.xn + B.yn) = A(xn+2 + p.xn+1 + q.xn ) + B(yn+2 + p.yn+1 + q.yn ) = 0 3
  4. Định nghĩa: x0 x1 Nếu 0 thì xn và yn độc lập tuyến tính y0 y1 Bổ đề 2: Nếu xn, yn là nghiệm riêng độc lập tuyến tính của (3) thì Un = A.xn + B.yn là nghiệm tổng quát của (3). Chứng minh: Gọi Un là một nghiệm bất kỳ của (3). Ta chứng minh rằng tồn tại Au và Bu sao cho Un = Au.xn + Bu.yn (Au, Bu là các hằng số phụ thuộc un). Ax0 + By0 = U0 Hệ phương trình Ax1 + By1 = U1 Có nghiệm duy nhất Au và Bu. U2 = p.U1 + q.U0 = Aux2 + Buy2. Chứng minh bằng quy nạp, ta có Un = Au.xn + Bu.yn mọi nghiệm của (3) đều biểu diễn qua xn và yn đ.p.c.m Ta tìm nghiệm riêng dưới dạng xn = λn (λ 0). Thay vào (3), ta có: λn+2 + p.λn+1 + q.λn = 0 λ2 + pλ + q = 0 (4). Phương trình (4) gọi là phương trình đặc trưng của (3). Trường hợp 1: Nếu (4) có hai nghiệm thực phân biệt λ1 và λ2 (3) có hai nghiệm riêng độc lập tuyến tính xn = λ1n và yn = λ2n . Nghiệm tổng quát Un = C1 λ1n + C2 λ2n Trường hợp 2: Nếu (4) có nghiệm kép là λ0, (3) có hai nghiệm riêng độc lập tuyến tính xn= λ0n và yn = n.λ0n . Nghiệm tổng quát Un = (C1+ nC2) λ0n p .i Trường hợp 3: Nếu (4) có hai nghiệm phức λ1,2 = =A Bi 2 B p ) và với r = A2 + B2 và α = arctgA . (A = ,B= 2 2 λ1,2 = r(cosα i.sinα) PT (3) có hai nghiệm riêng độc lập tuyến tính là xn = rn.cosnα và yn = rn.sinnα Nghiệm tổng quát Un = rn [C1 cosnα +C2 sinnα]. 4
  5. Ví dụ 1: Tìm nghiệm un+2 = 5un+1 + 6un biết u0 = 1, u1 = 0 Bài làm: Phương trình đặc trưng: λ2-5λ + 6 = 0 λ1 =1 và λ2 = 2 Vậy nghiệm tổng quát un = A + B.2n. u0 = A + B = 1 Hệ phương trình u 1 = A + 2B = 0 A = 2 và B = -1. n Vậy nghiệm riêng thoả mãn là un = 2 – 2 5 Ví dụ 2: Tìm nghiệm un+2 = 2 un+1 - un biết u0 = 0, u1 = 1 5 1 Bài làm: Phương trình đặc trưng: λ2- 2 λ+1 = 0 λ1 = 2 và λ2 = 2 1 Vậy nghiệm tổng quát un = A 2n + B.2n. u0 = A + B = 0 Hệ phương trình A 2 2 u1 = 2 + 2B = 1 A = -3 v à B = 3 . 2 Vậy nghiệm riêng cần tìm là un = 3 (2-n – 2n) Ví dụ 3: Tìm nghiệm un+2 = 10un+1 - 25un Bài làm: Phương trình đặc trưng: λ2- 10λ + 25 = 0 λ1 = λ2 = 5 Vậy nghiệm tổng quát un = (A + Bn)5n Ví dụ 4: Tìm nghiệm un+2 - 2un+1 + un = 0 biết u0 = 1, u1 = 2 Bài làm: Phương trình đặc trưng: λ2- 2λ+1 = 0 λ1 = λ2 = 1 Vậy nghiệm tổng quát un = A + Bn u0 = A = 1 Hệ phương trình u1 = A + B = 2 A = B = 1. Vậy nghiệm riêng cần tìm là un = 1 + n Ví dụ 5: Tìm nghiệm un+2 - un+1 + un = 0 Bài làm: Phương trình đặc trưng: λ2- λ+1 = 0 3 2 1 i3 1 3 (2)2 + ( 2 )2 = 1, tgα = 1 = 3 λ1,2 = ,r= 2 2 5
  6. α=3 λ1,2 = cos 3 i.sin 3 n. n. Vậy nghiệm tổng quát un = Acos 3 + Bsin 3 Ví dụ 6: Tìm nghiệm un+2 - 2un+1 + 4un = 0, u0 = u1 = 1 Bài làm: Phương trình đặc trưng: λ2- 2λ+4 = 0 12 +( 3 )2 = 2, tgα = 3 λ1,2 = 1 α=3 λ1,2 = 2(cos3 i. 3 , r = i.sin3 ) n. n. Vậy nghiệm tổng quát un = 2n(Acos 3 + Bsin 3 ) u0 = A = 1 Hệ phương trình u1 = 2(cos3 + Bsin3 ) = 1 A = 1 và B = 0. n. Vậy nghiệm riêng cần tìm là un = 2n.cos 3 2. Phương trình sai phân tuyến tính không thuần nhất Dạng Un+2 + pUn+1 + qUn = r (5) (r 0) Ta tìm nghiệm riêng U*n của (5) : ? r +) Nếu p+q -1 thì nghiệm riêng là : U*n = 1pq +) Nếu p+q = -1 rn Khi p -2 thì nghiệm riêng là : U*n = p2 rn 2 * Khi p = -2 thì nghiệm riêng là : U n = 2 Từ nghiệm của PT thuần nhất liên kết ta suy ra nghiệm tổng quát của (5). Trường hợp Un+2 + pUn+1 + qUn = f(n) ta xét ở dạng tổng quát cho PT sai phân tuyến tính hệ số hằng cấp k. V. Phƣơng trình sai phân tuyến tính cấp k hệ số hằng. 1. Phương trình sai phân tuyến tính thuần nhất cấp k hệ số hằng: Là phương trình có dạng: ak.Un+k + ak-1.Un+k-1 + … + a0.Un = 0 (6) Trong đó a0, a1, …, ak là các số thực. 6
  7. Ta tìm nghiệm riêng dưới dạng Un = λn, thay vào (6) ta có phương trình đặc trưng: ak.λk + ak-1.λk-1 + … + a0.λ = 0 (7) Trường hợp 1: Nếu (7) có k nghiệm thực phân biệt λ1, λ2, … λk ta có k nghiệm riêng độc lập tuyến tính x1n = λ1n, … xkn = λkn . Nghiệm tổng quát : Un = C1. λ1n + C2. λ2n + … + Ck. λkn Trường hợp 2: Nếu (7) có nghiệm bội, chẳng hạn λ1 có bội s và k-s nghiệm thực phân biệt: λ1 = λ2 = … = λs , ta thay thế s nghiệm riêng x1n, x2n, …, xsn tương ứng bằng: x1n = λ1n, x2n = nλ1n, … , xsn = ns-1.λ1n. Nghiệm tổng quát : Un = (C1+n C2 + … + ns-1Cs) λ1n + Cs+1 λ1n+...+ Ck. λkn Trường hợp 3: Nếu phương trình (7) có nghiệm phức, chẳng hạn λ1 = r(cosα +i.sinα) thì sẽ có nghiệm phức liên hợp λ2 = r(cosα - i.sinα) và k-2 nghiệm thực phân biệt, khi đó tương ứng ta thay thế x1n = rn.cosnα và x2n = rn.sinnα trong nghiệm tổng quát. Nghiệm tổng quát : Un = rn [C1. cosnα + C2. sinnα]+ C3. λ3n … + Ck. λkn Ví dụ 1: Tìm nghiệm un+3 – 10un+2 + 31un+1 - 30un = 0. Bài làm: Phương trình đặc trưng: λ3 -10λ2 + 31λ -30 = 0 λ1 =2, λ2 = 3 và λ3 = 5 Vậy nghiệm tổng quát un = A1.2n + A2.3n + A3.5n Ví dụ 2: Tìm nghiệm un+3 – 7un+2 + 16un+1 - 12un biết u0 = 0, u1 = 1, u2 = -1 Bài làm: Phương trình đặc trưng: λ3 - 7λ2 + 16λ -12 = 0 λ1 = λ2 = 2 và λ3 = 3 Vậy nghiệm tổng quát un = (A + n.B)2n + C.3n u0 = A + C = 0 Có hệ phương trình u1 = 2A + 2B + 3C = 1 u2 = 4(A + 2B) + 9C = -1 A = 5, B = 3 và C = -5. Vậy nghiệm riêng thoả mãn là un = (5 + 3n).2n – 5.3n Ví dụ 3: Tìm nghiệm un+3 – un = 0 Bài làm: Phương trình đặc trưng: λ3 -1= 0 1 i3 λ1 = 1, λ2,3 = 2 = cos3 i.sin3 n. n. Vậy nghiệm tổng quát un = A + Bcos 3 + Csin 3 7
  8. 2. Phương trình sai phân tuyến tính không thuần nhất cấp k hệ số hằng Là phương trình dạng: ak.Un+k + ak-1.Un+k-1 + … + a0.Un = fn (8) Trong đó a0, a1, …, ak là các số thực, fn 0n. Phương trình thuần nhất tương ứng ak.Un+k + ak-1.Un+k-1 + … + a0.Un = 0 (6). Bổ đề: Nghiệm tổng quát của phương trình (8) bằng nghiệm tổng quát của phương trình (6) cộng với nghiệm riêng bất kỳ của (8). Chứng minh: Giả sử vn là nghiệm tổng quát của (6) và xn là nghiệm riêng của (8). Đặt un = vn + xn. Ta có: ak.Un+k + ak-1.Un+k-1 + … + a0.Un = ak(vn+k + xn+k) + ak-1(vn+k-1 + xn+k-1) … + a0(vn + xn) = (ak.vn+k + ak-1.vn+k-1 + … + a0.vn)+(ak.xn+k + ak-1.xn+k-1+…+ a0.xn) = 0 + fn = fn un = vn + xn. Ngược lại hiệu 2 nghiệm riêng bất kỳ của (8) cũng là nghiệm riêng của (6). Vậy nghiệm tổng quát của (8) bằng nghiệm tổng quát của phương trình (6) cộng với nghiệm riêng bất kỳ của (8). Cách tìm nghiệm riêng xn fn = Pm(n) = bmnm + bm-1nm-1 + … + b1n + b0 Trường hợp 1: Nếu λ = 1 là nghiệm cấp s của phương trình đặc trưng ( s có thể nhận giá trị 0) thì nghiệm riêng có dạng xn= ns(cmnm + cm-1nm-1+…+ c1n + c0) và tìm ci bằng phương pháp hệ số bất định. Nếu λ = 1 không là nghiệm của phương trình đặc trưng thì nghiệm riêng có dạng xn= Cmnm + Cm-1nm-1+…+ C1n + C0 và tìm Ci bằng phương pháp hệ số bất định. fn = Pm(n).βn Trường hợp 2: Nếu λ = β là nghiệm cấp s của phương trình đặc trưng (s có thể nhận giá trị 0) thì nghiệm riêng có dạng xn= Qm(n).ns.βn, thay vào phương trình tìm Qm(n) bằng phương pháp hệ số bất định. Nếu λ = β không là nghiệm của phương trình đặc trưng thì nghiệm riêng có dạng xn= Qm(n).βn, thay vào phương trình tìm Qm(n) bằng phương pháp hệ số bất định. fn = Rl(n) + Pm(n).βn Trường hợp 3: Ta tìm nghiệm riêng dạng xn = x1n + x2n. 8
  9. Trong đó x1n là nghiệm riêng ứng với f1(n) = Rl(n) (đưa về trường hợp 1) và x2n là nghiệm riêng ứng với f2(n) = Pm(n).βn (đưa về trường hợp 2). 5 Ví dụ 1: Tìm một nghiệm riêng của phương trình un+2 – 2 un+1 + un = n2 + n + 1 5 1 Bài làm: Phương trình đặc trưng λ2 –2 λ+1 = 0 λ1= 2 và λ2 = 2 λ = 1 không là nghiệm ta tìm nghiệm riêng dạng xn= an2 + bn+ c Thay vào phương trình, ta có: 5 a(n+2)2+b(n+2)+c - 2 [a(n+1)2 +b(n+1) +c] + an2+bn+c = n2+ n+1. xn = -2n2 + 2n - 10 Đồng nhất hệ số a = -2, b =2 và c = -10 Ví dụ 2: Tìm một nghiệm riêng của phương trình un+2 – un = 6n2 + 12n + 8 Bài làm: Phương trình đặc trưng λ2 –1 = 0 λ1= 1 và λ2 = -1 λ = 1 là nghiệm đơn ta tìm nghiệm riêng dạng xn= n(an2+bn+c) x n = n3 Thay vào phương trình a = 1, b = c = 0 5 Ví dụ 3: Tìm một nghiệm riêng của phương trình un+2 – 2 un+1 + un = 3n 5 1 Bài làm: Phương trình đặc trưng λ2 –2 λ+1 = 0 λ1= 2 và λ2 = 2 ta tìm nghiệm riêng dạng xn= A.3n λ = 3 không là nghiệm 5 2 2 Thay vào phương trình, ta có: A.3n+2 - 2 A.3n+1 + A.3n = 3n A = 5 xn = 5 .3n un+2 – un+1 - 2un = -3. 2n Ví dụ 4: Tìm một nghiệm riêng của phương trình Bài làm: Phương trình đặc trưng λ2 – λ - 2 = 0 λ1= 2 và λ2 = -1 λ = 2 là nghiệm đơn ta tìm nghiệm riêng dạng xn= A.n.2n 1 -n Thay vào PT, ta có: A(n+2)2n+2 – A(n+1)2n+1 – 2A.n.2n = -3.2n A = - 2 xn = 2 .2n Ví dụ 5: Tìm một nghiệm riêng của phương trình 5 un+2 – 2 un+1 + un = n2 + n + 1 + 3n 2 Bài làm: Áp dụng ví dụ 1 và ví dụ 3 nghiệm riêng xn = -2n2 + 2n – 10 + 5 .3n 6. Ứng dụng của phƣơng trình sai phân 9
nguon tai.lieu . vn