Ứng dụng chuỗi số giải phương trình vi phân

Một số phương trình vi phân rất khó (nếu không muốn nói là không thể) tìm nghiệm ở dạng tổ hợp của các hàm số sơ cấp xác định. Điều này cũng xảy ra ngay cả khi các phương trình vi phân có dạng rất đơn giản. Ví dụ như phương trình sau:

y''- 2x.y' + y = 0 (1)

Đây là phương trình vi phân cấp hai, hệ số hàm nhưng ta không thể tìm được 1 nghiệm riêng dưới dạng hàm số sơ cấp. Tuy nhiên, việc giải các phương trình như dạng phương trình (1) là rất quan trọng vì nó nảy sinh từ các vấn đề, các bài toán của vật lý, cụ thể, nó liên quan đến phương trình Schrödinger trong cơ học lượng tử. Vì vậy, ta cần thiết phải xây dựng các phương pháp nhằm tìm nghiệm cho các phương trình dạng này. Một trong các phương pháp thông dụng là ứng dụng lý thuyết chuỗi để tìm nghiệm của phương trình dưới dạng chuỗi lũy thừa:

y = \sum\limits_{n=0}^{\infty}c_n.x^n = c_0 + c_1.x + c_2.x^2 + c_3.x^3 + ... + c_n.x^n (2)

Cơ sở Toán học của phương pháp này là ta thay thế biểu thức (2) vào phương trình vi phân và từ đó xác định giá trị của các hằng số c_0,c_1,c_2,... sao cho nó nghiệm đúng phương trình vi phân. Tuy nhiên, việc này chỉ có giá trị khi chuỗi lũy thừa ứng với các hệ số tìm được phải là chuỗi hội tụ.

Ta nhắc lại một số điều thường gặp đối với chuỗi lũy thừa:

Trong khoảng hội tụ của chuỗi, ta có thể lấy đạo hàm và tích phân từng số hạng của chuỗi, chuỗi mới nhận được (sau khi lấy đạo hàm hoặc tích phân) cũng có bán kính hội tụ như chuỗi ban đầu.

Trước khi sử dụng chuỗi lũy thừa để giải phương trình 1, chúng ta sẽ minh họa phương pháp này bằng một ví dụ đơn giản hơn phương trình (1). Xét phương trình:

y'' + y = 0 (3)

Theo phương pháp sơ cấp, ta đã biết nghiệm của phương trình (3) có dạng:

y = C_1.cosx + C_2.sinx

Bây giờ, sử dụng phương pháp chuỗi số. Ta giả sử nghiệm của phương trình (3) có dạng:

y = c_0 + c_1.x + c_2.x^2 + c_3.x^3 + ... + c_n.x^n = \sum\limits_{n=0}^{\infty}c_n.x^n (2)

Khi đó, ta có:

y'' = 2c_2 + 2.3c_3x + 3.4c_4x^2 + ... + n(n-1)c_nx^{n-2} \\ \qquad = \sum\limits_{n=2}^{\infty}n(n-1)c_nx^{n -2} (4)

Ta lưu ý: Chỉ số của chuỗi thứ (4) bắt đầu với n = 2, do đó để thuận tiện cho việc so sánh các hệ số của y và y” được dễ dàng, ta sẽ viết chuỗi (4) với chỉ số bắt đầu với n = 0, nghĩa là:

y'' = \sum\limits_{n=0}^{\infty}(n+2)(n+1)c_{n+2}x^{n} (5)

Đây là kỹ thuật quan trọng nhất trong việc sử dụng phương pháp chuỗi số. Sở dĩ, ta làm được việc này là do ta đặt m = n -2, hay n = m + 2. Khi đó ta có:

y'' = \sum\limits_{m=0}^{\infty}(m+2)(m+1)c_{m+2}x^{m} (6)

Khi đó, theo lý thuyết chuỗi số ta có hai chuỗi (5) và (6) là như nhau. Như vậy, ta thế biểu thức (2) và (5) vào phương trình (3) ta có:

{ \sum\limits_{n=0}^{\infty}(n+2)(n+1)c_{n+2}x^{n}} + { \sum\limits_{n=0}^{\infty}c_n.x^n} = 0

Hay ta có:

\sum\limits_{n=0}^{\infty}[(n+2)(n+1)c_{n+2}+c_n]x^n = 0 (7)

Theo phương pháp hệ số bất định, hai chuỗi số muốn bằng nhau thì từng hệ số tương ứng phải bằng nhau. Vì vậy, hệ số của x^n ở biểu thức (7) phải bằng 0. Hay:

(n+2)(n+1).c_{n+2} + c_n = 0

Từ đây ta có:

c_{n+2} = - { \dfrac{c_n}{(n+1)(n+2)}}

Như vậy, ta có được công thức truy hồi. Do đó:

Với n = 0: c_2 = - { \dfrac{c_0}{1.2}}

Với n = 1: c_3 = - { \dfrac{c_1}{2.3}}

Với n = 2: c_4 = - { \dfrac{c_2}{3.4}} = { \dfrac{c_0}{1.2.3.4}} = { \dfrac{c_0}{4!}}

Với n = 3: c_5 = - { \dfrac{c_3}{4.5}} = { \dfrac{c_1}{2.3.4.5}} = { \dfrac{c_1}{5!}}

Với n = 4: c_6 = - { \dfrac{c_4}{5.6}} = - { \dfrac{c_0}{4!.5.6}} = - { \dfrac{c_0}{6!}}

Với n = 5: c_7 = - { \dfrac{c_5}{6.7}} = - { \dfrac{c_1}{5!.6.7}} = - { \dfrac{c_1}{7!}}

Vì vậy, theo quy luật trên chúng ta có:

Với các hệ số chẵn: c_{2k} = (-1)^{k} { \dfrac{c_0}{(2k)!}}

Với các hệ số lẻ: c_{2k+1} = (-1)^{k} { \dfrac{c_1}{(2k+1)!}}

Do đó, thế vào chuỗi (2) ta có:

y = c_0 + c_1.x + c_2.x^2 + c_3.x^3 + ... + c_n.x^n \\ = c_0(1- { \dfrac{x^2}{2!}} + { \dfrac{x^4}{4!}} - { \dfrac{x^6}{6!}} + ... + (-1)^{n} { \dfrac{x^{2n}}{(2n)!}} + ... ) + \\ c_1(x- { \dfrac{x^3}{3!}} + { \dfrac{x^5}{5!}} - { \dfrac{x^7}{7!}} + ... + (-1)^{n} { \dfrac{x^{2n+1}}{(2n+1)!}} + ... )

Hay:

y = c_0{ \sum\limits_{n=0}^{\infty}(-1)^{n}{ \dfrac{x^{2n}}{(2n)!}}} + c_1{ \sum\limits_{n=0}^{\infty}(-1)^{n}{ \dfrac{x^{2n+1}}{(2n+1)!}}}

Như vậy, ta có kết quả của phương trình phụ thuộc vào 2 chuỗi số với c0 và c1 là 2 hằng số tùy ý. Và, rõ ràng, cả hai chuỗi số đều hội tụ với mọi x thuộc R. Nhận xét:

  • Nếu chúng ta chú ý các dạng của chuỗi Taylor – Maclaurin, thì hai chuỗi số vừa tìm được ở trên chính là khai triển Maclaurin của hai hàm số cosx và sinx. Do đó, nghiệm của phương trình sẽ là:

y = c_0.cosx + c_1.sinx

  • Tuy nhiên, thông thường, ít khi nào ta có thể biểu diễn các chuỗi lũy thừa trong công thức nghiệm dưới dạng các hàm số sơ cấp đã biết.

Trang: 1 2

Thảo luận

23 bình luận về “Ứng dụng chuỗi số giải phương trình vi phân”

  1. các bạn ơi cho mình hỏi cách giải bài vi phân sau: y” – 4y’= e^x ((-4x+4)cosx-(2x+6)sinx)

    Thích

    Posted by mai thị chi | 15/05/2015, 11:03
    Reply to this comment

  2. cho em hỏi một bài này: Q=3L^(1/3)K^(2/3), sử dụng vi phân toàn phần tính gần đúng mức thay đổi của Q khi K tăng từ 1000 lên 1000,5 và L giảm từ 125 xuống 124

    Thích

    Posted by Dương Quế | 13/12/2014, 16:09
    Reply to this comment

  3. có ai giúp mình giải phương trình vi phân này với:
    huhu, mình làm k ra.
    ( x * y’ – y) arctan (y/x) = x
    thank nhiều!!!

    Thích

    Posted by tienchang1 | 30/09/2011, 11:31
    Reply to this comment

  4. Ai giúp em giải bài này với ạ: y’ = ycos(x-y)

    Thích

    Posted by quyen | 02/01/2010, 21:25
    Reply to this comment

  5. mình có bài này, không biết có thể sử dụng kiến thức mà bạn nói ở trên để giải không:

    for each integer n \ge 2 , we consider the iteration equation:
    u_n = u_{n -1} + 2u_{n-2} + (-1)^n
    if we put u_o = u_1 = 1 , we obtain a sequence (u_n)_n and for each complex z the power series:
    u(z) = \sum_{n = 0}^{\infty}u_n.z^n
    if we suppose that R,bán kính hội tụ của chuỗi lũy thừa với |z| < R, ta có: u(z) = \dfrac{1 + z +z^2}{(1-3z)(1+z)^2} là đúng hay sai

    Thích

    Posted by ntd | 14/08/2009, 22:21
    Reply to this comment

    • Bài này không phải là dạng ứng dụng chuỗi số để giải phương trình vi phân mà đây là dạng xây dựng chuỗi số phức khi biết dạng công thức truy hồi của nó. Bài này thuộc phần khai triển Taylor-Maclaurin cho chuỗi số phức. Để giải bài này, cần sử dụng lý thuyết thặng dư cho chuỗi phức.

      Thích

      Posted by 2Bo02B | 16/08/2009, 10:17
      Reply to this comment

  6. Rất vui khi biết bạn đang làm về lĩnh vực phương trình vi phân. Hướng của mình đang làm là lý thuyết bài toán biên nhiều điểm cho hệ phương trình vi phân hàm.nhưng cũng còn nhiều trở ngại. Thật bất ngờ khi biết giáo sư Ngô Ngọc Diệp đang nghiên cứu về lĩnh vực này, vì trước đây, GS Diệp chuyên nghiên cứu về lý thuyết nhóm Lie để khảo sát các đa tạp. Nếu bạn có những nguồn tài liệu nào về vấn đề này (tiếng Anh càng tốt), mong bạn share giúp mình.

    Cảm ơn bạn đã góp ý Về việc viết tên Tiếng Việt. Tuy nhiên, hiện tại, đa số các giáo trình đều viết tên các công thức, định lý theo nguyên bản tên của nhà Toán học. Điều này giúp các bạn sinh viên dễ dàng tiếp cận với các nguồn giáo trình tiếng Anh hoặc các tài liệu liên quan. Mặc dù vậy, mình sẽ thêm phiên âm tiếng Việt ngay sau các tên nhà Toán học để mọi người dễ theo dõi.

    Thân mến,

    P.S: Ngoài các cách trên, ta còn có cách ứng dụng dãy Fourier (Phu-ri-ê) và phép biến đổi Laplace để giải ptvp nữa

    Thích

    Posted by 2Bo02B | 13/03/2009, 16:42
    Reply to this comment

  7. Xin chào các bạn!
    Mình là Huấn đang học ở Liên bang Nga và cũng đang nghiên cứu về phương trình vi phân. Của mình là ứng dụng lý thuyết nhóm Li vào dể giải các phương trình vi phân. Lĩnh vực này ở Việt Nam có giáo sư Ngô Ngọc Diệp là làm gần nhất. Bạn viết rất cụ thể về phương pháp dùng dãy số để giải phương trình. Tuy nhiên khi nêu tên các nhà khoa học nên viết bằng tiếng Việt, tiếng anh sẽ có ít người hiểu lắm. Mục đích của bạn là viết lên để cho mọi người cùng đọc đúng không. Vì thế nên viết tất cả bằng tiếng Việt. Mình bật mí là cùng với chuỗi có thêm ứng dụng dãy Phu-ri-ê để giải phương trình đấy. Phương pháp này ứng dụng với các phương trình có nghiệm tuần hoàn ví dụ như bài toán về phương trình dao động của dây đàn ghi ta!

    Thích

    Posted by Hoàng Ngự Huấn! | 13/03/2009, 04:47
    Reply to this comment
« Bình luận cũ hơn

Bình luận về bài viết này

Translators & RSS

English French RussiaMaths 4 Physics (M4Ps)

Bạn hãy nhập địa chỉ email của mình để đăng ký theo dõi tin tức từ blog này và nhận những bài viết mới nhất qua địa chỉ email.

Join 2 787 other subscribers

Đôi lời

Bạn có thể theo dõi các lời bình liên quan đến lời bình của mình qua email bằng cách chọn dòng thông báo Báo cho bạn khi có người bình luận tiếp theo đề tài này bằng điện thư mỗi khi viết 1 lời bình.

Rất mong các bạn viết lời nhắn bằng tiếng việt có dấu nhé.

Để viết tiếng việt có dấu bạn dùng font chữ Unicode và bảng mã là Unicode UTF-8.

Để biết cách gõ công thức Toán học trong các lời nhắn ở trang web này, mời bạn đọc bài hướng dẫn tại đây hoặc bạn có thể xem bài hướng dẫn dùng MathType tại đây và bài tạo công thức trực tuyến tại đây

Get Well

Lời nhắn mới nhất

Dương Khánh Uyên trong Trang 2
Trần Thái An trong Trang 2
Chúc Chúc trong Xác suất có điều kiện
Hoang Anh trong Khai triển Taylor – Macl…
Trần Trung Đức trong Mẹo phân tích nhanh 1 phân…
Nhung Duong trong Trang 2
khoi trong Khai triển Taylor – Macl…
Minh pham trong Chuỗi Fourier Sine và Cos…
Minh Phạm trong Chuỗi Fourier
Anh Tuấn trong Cực trị (không điều kiện) của…