MathLabs

Phương trình vi phân và Hệ động lực

Phương trình vi phân tuyến tính

Phương trình cấp cao có nguyên lý chồng chất, giải bằng phương trình đặc trưng.

Trực giácTrực giác: dao động và nguyên lý chồng chất

Vật nặng gắn lò xo, con lắc dao động biên độ nhỏ, dòng điện trong mạch RLC đều tuân theo cùng một dạng quy luật: hàm số chưa biết cùng các đạo hàm của nó kết hợp tuyến tính với nhau. Nếu biết hai chuyển động riêng lẻ đều thỏa mãn quy luật đó thì tổng (hay tổ hợp có hệ số) của chúng cũng thỏa mãn — đó là nguyên lý chồng chất. Đây là lý do ta chỉ cần hai nghiệm cơ bản là dựng được mọi chuyển động có thể xảy ra.

Parabol của đa thức đặc trưng r^2-3r+2, cắt trục r tại r=1 và r=2.
Dao động điều hòa tắt dần y′′+2ζy′+y=0y'' + 2\zeta y' + y = 0: nghiệm dao động bên trong đường bao hàm mũ giảm dần ±e−ζt\pm e^{-\zeta t}.

Phổ thôngPhát biểu chính xác: phương trình tuyến tính cấp hai hệ số hằng

Định nghĩa: Phương trình vi phân tuyến tính thuần nhất

Phương trình cấp hai dạng ay′′+by′+cy=0ay'' + by' + cy = 0, với a,b,ca,b,c là hằng số và a≠0a\neq 0, được gọi là tuyến tính thuần nhất vì vế phải bằng 0 và yy chỉ xuất hiện bậc nhất ở mọi số hạng.

ay′′+by′+cy=0ay'' + by' + cy = 0

Vì đạo hàm của erxe^{rx} chỉ là chính nó nhân với rr, ta thử nghiệm dạng y=erxy=e^{rx}. Thay vào phương trình rồi chia cho thừa số chung erx≠0e^{rx}\neq 0, phương trình vi phân trở thành một phương trình bậc hai đơn giản theo rr:

ar2+br+c=0ar^2+br+c=0
Ba trường hợp của phương trình đặc trưng ar2+br+c=0ar^2+br+c=0 và nghiệm tổng quát tương ứng
Biệt thứcNghiệmNghiệm tổng quát
Δ>0\Delta>0Hai nghiệm thực phân biệt r1≠r2r_1\neq r_2y=c1er1x+c2er2xy = c_1 e^{r_1 x} + c_2 e^{r_2 x}
Δ=0\Delta=0Một nghiệm thực kép rry=(c1+c2x)erxy = (c_1 + c_2 x) e^{r x}
Δ<0\Delta<0Nghiệm phức liên hợp r=α±iβr = \alpha \pm i\betay=eαx(c1cos⁡(βx)+c2sin⁡(βx))y = e^{\alpha x}\left(c_1 \cos(\beta x) + c_2 \sin(\beta x)\right)

Đại họcCấp đại học: định lý và chứng minh

Nếu y1y_1 và y2y_2 đều thỏa mãn ay′′+by′+cy=0ay'' + by' + cy = 0 thì với mọi hằng số c1,c2c_1,c_2, tổ hợp y=c1y1+c2y2y=c_1y_1+c_2y_2 cũng thỏa mãn phương trình. Hơn nữa, nếu định thức Wronski W(y1,y2)=y1y2′−y2y1′W(y_1,y_2) = y_1y_2' - y_2y_1' khác 0 tại một điểm nào đó thì mọi nghiệm của phương trình đều có dạng này.

Vì sao đúng?

Phép lấy đạo hàm là một phép toán tuyến tính: đạo hàm của tổng bằng tổng các đạo hàm, hằng số được đưa ra ngoài. Vì vế trái của phương trình chỉ gồm yy, y′y', y′′y'' nhân hằng số rồi cộng lại, nên khi thay một tổ hợp của hai nghiệm vào, ta chỉ cộng hai số 0 lại với nhau.

Chứng minh

Đặt L[y]=ay′′+by′+cyL[y] = ay''+by'+cy. Vì đạo hàm thỏa mãn (u+v)′=u′+v′(u+v)'=u'+v' và (ku)′=ku′(ku)'=ku' với hằng số kk, toán tử LL là tuyến tính: L[c1y1+c2y2]=a(c1y1+c2y2)′′+b(c1y1+c2y2)′+c(c1y1+c2y2)=c1(ay1′′+by1′+cy1)+c2(ay2′′+by2′+cy2)=c1L[y1]+c2L[y2]L[c_1y_1+c_2y_2] = a(c_1y_1+c_2y_2)'' + b(c_1y_1+c_2y_2)' + c(c_1y_1+c_2y_2) = c_1(ay_1''+by_1'+cy_1) + c_2(ay_2''+by_2'+cy_2) = c_1L[y_1]+c_2L[y_2].

Vì y1,y2y_1,y_2 là nghiệm của phương trình nên L[y1]=0L[y_1]=0 và L[y2]=0L[y_2]=0, do đó L[c1y1+c2y2]=c1⋅0+c2⋅0=0L[c_1y_1+c_2y_2]=c_1\cdot 0+c_2\cdot 0=0: tổ hợp cũng là một nghiệm. Điều này chứng minh nửa đầu.

Với nửa sau, lý thuyết phương trình vi phân tuyến tính khẳng định không gian nghiệm của phương trình cấp hai có đúng hai chiều (một điều kiện đầu y(x0)=y0y(x_0)=y_0, y′(x0)=y0′y'(x_0)=y_0' xác định duy nhất một nghiệm, cho hai tham số tự do). Định thức Wronski khác 0, W(y1,y2)=y1y2′−y2y1′W(y_1,y_2)=y_1y_2'-y_2y_1', tại một điểm x0x_0 nghĩa là hệ tuyến tính c1y1(x0)+c2y2(x0)=y0c_1y_1(x_0)+c_2y_2(x_0)=y_0, c1y1′(x0)+c2y2′(x0)=y0′c_1y_1'(x_0)+c_2y_2'(x_0)=y_0' có nghiệm (c1,c2)(c_1,c_2) duy nhất với mọi dữ kiện đầu. Vậy mọi nghiệm đều trùng với một c1y1+c2y2c_1y_1+c_2y_2 nào đó, tức là {y1,y2}\{y_1,y_2\} trải hết toàn bộ không gian nghiệm hai chiều.

Nghiệm tổng quát của ay′′+by′+cy=0ay'' + by' + cy = 0 hoàn toàn được xác định bởi nghiệm của phương trình đặc trưng ar2+br+c=0ar^2+br+c=0: nghiệm thực phân biệt cho y=c1er1x+c2er2xy = c_1 e^{r_1 x} + c_2 e^{r_2 x}, nghiệm thực kép cho y=(c1+c2x)erxy = (c_1 + c_2 x) e^{r x}, nghiệm phức liên hợp r=α±iβr = \alpha \pm i\beta cho y=eαx(c1cos⁡(βx)+c2sin⁡(βx))y = e^{\alpha x}\left(c_1 \cos(\beta x) + c_2 \sin(\beta x)\right).

Vì sao đúng?

Hàm mũ là khối xây dựng tự nhiên vì đạo hàm của erxe^{rx} chỉ là chính nó nhân rr; đưa nghiệm thử này vào phương trình tuyến tính hệ số hằng biến giải tích thành đại số — tìm nghiệm một đa thức.

Chứng minh

Thử nghiệm dạng y=erxy=e^{rx}. Khi đó y′=rerxy'=re^{rx} và y′′=r2erxy''=r^2e^{rx}; thay vào ay′′+by′+cy=0ay'' + by' + cy = 0 ta được (ar2+br+c)erx=0(ar^2+br+c)e^{rx}=0. Vì erxe^{rx} không bao giờ bằng 0 nên buộc phải có phương trình đặc trưng ar2+br+c=0ar^2+br+c=0. Theo công thức nghiệm bậc hai, r=−b±b2−4ac2ar = \dfrac{-b \pm \sqrt{b^2-4ac}}{2a}.

Trường hợp Δ>0\Delta>0: hai nghiệm r1≠r2r_1\neq r_2 cho hai nghiệm riêng er1xe^{r_1x} và er2xe^{r_2x}. Định thức Wronski của chúng là (r2−r1)e(r1+r2)x≠0(r_2-r_1)e^{(r_1+r_2)x}\neq 0, nên theo định lý chồng chất, nghiệm tổng quát là y=c1er1x+c2er2xy = c_1 e^{r_1 x} + c_2 e^{r_2 x}.

Trường hợp Δ=0\Delta=0: có một nghiệm kép r=−b/(2a)r=-b/(2a), chỉ cho một nghiệm mũ erxe^{rx}. Để tìm nghiệm độc lập thứ hai, thử y2=v(x)erxy_2=v(x)e^{rx} (hạ bậc) rồi thay vào phương trình; vì rr là nghiệm kép của ar2+br+car^2+br+c nên các số hạng chứa vv và v′v' triệt tiêu, chỉ còn av′′erx=0av''e^{rx}=0, tức v′′=0v''=0, nên v=c1+c2xv=c_1+c_2x. Điều này cho nghiệm thứ hai xerxxe^{rx}, và nghiệm tổng quát y=(c1+c2x)erxy = (c_1 + c_2 x) e^{r x}.

Trường hợp Δ<0\Delta<0: nghiệm là hai số phức liên hợp r=α±iβr = \alpha \pm i\beta với α=−b/(2a)\alpha=-b/(2a), β=−Δ/(2a)\beta=\sqrt{-\Delta}/(2a). Công thức Euler eiβx=cos⁡(βx)+isin⁡(βx)e^{i\beta x}=\cos(\beta x)+i\sin(\beta x) biến hai nghiệm phức e(α±iβ)xe^{(\alpha\pm i\beta)x} thành hai nghiệm thực độc lập eαxcos⁡(βx)e^{\alpha x}\cos(\beta x) và eαxsin⁡(βx)e^{\alpha x}\sin(\beta x) (bằng cách lấy phần thực và phần ảo, cũng là nghiệm vì LL có hệ số thực), cho nghiệm tổng quát y=eαx(c1cos⁡(βx)+c2sin⁡(βx))y = e^{\alpha x}\left(c_1 \cos(\beta x) + c_2 \sin(\beta x)\right).

Đại họcỨng dụng thực tiễn và Ví dụ minh họa

Ba trường hợp của phương trình đặc trưng không chỉ là đại số — chúng mô tả ba dáng điệu vật lý khác nhau của hệ lò xo-vật nặng có ma sát my′′+cy′+ky=0my'' + cy' + ky = 0 (khối lượng mm, hệ số cản cc, độ cứng kk): dao động tắt dần rung lắc rồi tắt (Δ<0\Delta<0), quá tắt dần trở về vị trí cân bằng mà không rung lắc chút nào (Δ>0\Delta>0), và tắt dần tới hạn, cách trở về nhanh nhất mà không rung lắc (Δ=0\Delta=0) — chính là mục tiêu thiết kế của giảm xóc ô tô hay bộ đóng cửa tự động.

Ví dụ: Tìm nghiệm tổng quát: nghiệm thực phân biệt

Giải phương trình y′′−5y′+6y=0y'' - 5y' + 6y = 0.

Lời giải

Phương trình đặc trưng là r2−5r+6=0r^2-5r+6=0.

Phân tích thành nhân tử, (r−2)(r−3)=0(r-2)(r-3)=0, nên các nghiệm là r1=2r_1=2 và r2=3r_2=3: hai nghiệm thực phân biệt.

Theo định lý ở trên (trường hợp Δ>0\Delta>0), nghiệm tổng quát là y=c1e2x+c2e3xy = c_1e^{2x}+c_2e^{3x}, với c1,c2c_1,c_2 là các hằng số tùy ý xác định bởi điều kiện đầu.

Ví dụ: Hệ số bất định cho phương trình có vế phải

Tìm một nghiệm riêng của y′′−y=e2xy'' - y = e^{2x}.

Lời giải

Phương trình thuần nhất y′′−y=0y''-y=0 có phương trình đặc trưng r2−1=0r^2-1=0, nghiệm r=±1r=\pm 1; số hạng cưỡng bức e2xe^{2x} có số mũ 22, không trùng với các nghiệm này, nên không có cộng hưởng.

Thử nghiệm dạng yp=Ae2xy_p = Ae^{2x}. Khi đó yp′′=4Ae2xy_p''=4Ae^{2x}, thay vào ta được 4Ae2x−Ae2x=3Ae2x4Ae^{2x}-Ae^{2x}=3Ae^{2x}, phải bằng e2xe^{2x}.

So sánh hệ số, 3A=13A=1, nên A=13A=\tfrac{1}{3}, cho nghiệm riêng yp=13e2xy_p=\tfrac{1}{3}e^{2x}. Nghiệm tổng quát đầy đủ là y=c1ex+c2e−x+13e2xy=c_1e^{x}+c_2e^{-x}+\tfrac{1}{3}e^{2x}.

Nghiệm tổng quát của y′′−7y′+12y=0y''-7y'+12y=0 là gì?

Với y′′+2y′+5y=0y''+2y'+5y=0, phương trình đặc trưng có loại nghiệm nào, và nghiệm tổng quát là gì?

Giảm xóc ô tô, mô hình hóa bởi my′′+cy′+ky=0my'' + cy' + ky = 0, được thiết kế để đưa thân xe về vị trí cân bằng nhanh nhất mà không hề nảy lên. Kỹ sư nên hướng tới trường hợp nào của phương trình đặc trưng?

Để tìm nghiệm riêng của y′′−y=e2xy''-y=e^{2x} bằng phương pháp hệ số bất định, nên dùng nghiệm thử dạng nào?

Tài liệu tham khảo

  1. William E. Boyce, Richard C. DiPrima (2017). Elementary Differential Equations and Boundary Value Problems
  2. Lawrence Perko (2001). Differential Equations and Dynamical Systems