Lược Đồ Sai Phân Khác Thường Giải Phương Trình Vi Phân

Luận văn thạc sĩ phân tích lược đồ sai phân khác thường giải một số phương trình vi phân, đánh giá thực trạng, chỉ ra hạn chế, đề xuất giải pháp khả thi cho thực tiễn.

• Tác giả: Hoàng Mạnh Tuấn
• Người hướng dẫn: GS. Đặng Quang Á
• Trường học: Đại Học Quốc Gia Hà Nội
• Chuyên ngành: Toán Ứng Dụng
• Đề tài: Lược Đồ Sai Phân Khác Thường Giải Phương Trình Vi Phân
• Loại tài liệu: Luận Văn Thạc Sĩ Khoa Học
• Năm xuất bản: 2015
• Địa điểm: Hà Nội

Tóm tắt

I. Tổng Quan Lược Đồ Sai Phân Khác Thường Bí Quyết Tiếp Cận

Bài viết này khám phá lược đồ sai phân khác thường, một phương pháp mạnh mẽ để giải phương trình vi phân. Phương pháp này nổi lên như một giải pháp thay thế cho các phương pháp sai phân hữu hạn truyền thống khi các phương pháp này gặp khó khăn trong việc bảo toàn các tính chất quan trọng của nghiệm, dẫn đến tính không ổn định. Lược đồ sai phân khác thường được thiết kế để khắc phục những hạn chế này bằng cách sử dụng các công thức rời rạc hóa đặc biệt, thường là phi địa phương, để xấp xỉ các đạo hàm. Cách tiếp cận này cho phép bảo toàn tốt hơn các tính chất như tính ổn định, tính đơn điệu và tính bảo toàn năng lượng, đặc biệt quan trọng trong việc mô phỏng các mô hình hóa toán học vật lý và sinh học. Luận văn này trình bày một cách hệ thống các kết quả từ các tác giả nước ngoài trong 20 năm qua.

1.1. Mục Tiêu Phạm Vi Nghiên Cứu Lược Đồ Sai Phân

Luận văn tập trung vào việc xây dựng lược đồ sai phân khác thường cho bài toán giá trị ban đầu đối với phương trình vi phân thường. Mục tiêu là trình bày chi tiết cách các lược đồ này có thể được thiết kế để bảo toàn các tính chất định tính quan trọng của nghiệm. Điều này đặc biệt quan trọng khi giải các bài toán trên miền thời gian lớn, nơi các phương pháp truyền thống có thể dẫn đến kết quả không chính xác. Các nội dung chính bao gồm việc trình bày lại các kết quả tiêu biểu của các tác giả nước ngoài trong vòng 20 năm trở lại đây. Cấu trúc luận văn bao gồm ba chương.

1.2. Kiến Thức Nền Tảng Quan Trọng Về Phương Trình Vi Phân

Để hiểu rõ lược đồ sai phân khác thường, cần nắm vững kiến thức cơ bản về phương trình vi phân. Điều này bao gồm các khái niệm như sự tồn tại và duy nhất nghiệm, tính liên tục của nghiệm theo dữ liệu ban đầu và các định nghĩa liên quan đến điểm bất động (điểm cân bằng) và tính ổn định. Một điểm bất động là một giá trị mà tại đó đạo hàm bằng không, và tính ổn định xác định xem một nghiệm gần điểm bất động có xu hướng tiến về nó (ổn định) hay rời xa nó (không ổn định). Các định nghĩa này cung cấp nền tảng để phân tích tính ổn định của các lược đồ sai phân.

II. Thách Thức Vì Sao Cần Lược Đồ Sai Phân Khác Thường

Các phương pháp số truyền thống, như phương pháp Euler và phương pháp Runge-Kutta, có thể gặp khó khăn trong việc giải phương trình vi phân, đặc biệt là khi đối mặt với các bài toán biên hoặc các hệ động lực phức tạp. Một trong những thách thức chính là bảo toàn các tính chất định tính quan trọng của nghiệm, chẳng hạn như tính ổn định, tính đơn điệu và tính bảo toàn năng lượng. Khi các phương pháp sai phân hữu hạn thông thường không bảo toàn được các tính chất này, chúng có thể dẫn đến các giải pháp không chính xác, không ổn định hoặc thậm chí không có ý nghĩa vật lý. Lược đồ sai phân khác thường ra đời như một giải pháp để giải quyết những hạn chế này.

2.1. Hiện Tượng Không Ổn Định Số Trong Giải Phương Trình

Hiện tượng không ổn định số xảy ra khi nghiệm của phương trình sai phân (thu được từ quá trình rời rạc hóa) không phản ánh chính xác các tính chất của nghiệm của phương trình vi phân tương ứng. Điều này có thể dẫn đến các giải pháp dao động không thực tế, hoặc nghiệm số có xu hướng tăng lên vô hạn, ngay cả khi nghiệm thực tế của phương trình vi phân là ổn định và bị chặn. Các phân tích cũng cho thấy rằng, hiện tượng không ổn định số cũng xảy ra khi ta sử dụng các kỹ thuật tinh vi hơn để xây dựng các lược đồ sai phân bình thường, chẳng hạn sử dụng phương pháp Taylor hoặc phương pháp Runge-Kutta.

2.2. Hạn Chế Của Phương Pháp Sai Phân Hữu Hạn Thông Thường

Các phương pháp sai phân hữu hạn thông thường thường dựa trên việc xấp xỉ các đạo hàm bằng các công thức sai phân đơn giản, chẳng hạn như sai phân tiến, sai phân lùi hoặc sai phân trung tâm. Mặc dù các công thức này có thể chính xác đối với các bước lưới nhỏ, nhưng chúng có thể không bảo toàn được các tính chất quan trọng của nghiệm khi sử dụng các bước lưới lớn hơn. Điều này đặc biệt đúng đối với các phương trình vi phân có tính phi tuyến cao, hoặc các hệ động lực phức tạp, nơi các tính chất định tính của nghiệm có thể nhạy cảm với các sai số nhỏ.

III. Phương Pháp Xây Dựng Lược Đồ Sai Phân Khác Thường Hiệu Quả

Lược đồ sai phân khác thường được xây dựng dựa trên một tập hợp các quy tắc cụ thể, được đề xuất bởi R.E. Mickens, nhằm giải quyết các vấn đề về tính ổn định số. Các quy tắc này bao gồm việc sử dụng các công thức rời rạc hóa phức tạp hơn cho các đạo hàm và việc rời rạc hóa phi địa phương các số hạng phi tuyến. Rời rạc hóa phi địa phương nghĩa là gì? Đó là rời rạc hóa hàm số dựa trên giá trị của hàm tại một số điểm trên lưới rời rạc thay vì rời rạc hóa địa phương trong các lược đồ sai phân hữu hạn.

3.1. Quy Tắc Mickens Nền Tảng Của Lược Đồ Sai Phân

Quy tắc Mickens là nền tảng cho việc xây dựng lược đồ sai phân khác thường. Hai quy tắc quan trọng nhất là: (1) các đạo hàm nên được rời rạc hóa bằng các công thức phức tạp hơn các công thức sai phân thông thường. (2) Các số hạng phi tuyến nên được rời rạc hóa một cách không địa phương, tức là bằng cách sử dụng các giá trị của hàm tại các điểm khác nhau trên lưới.

3.2. Rời Rạc Hóa Phi Địa Phương Yếu Tố Then Chốt

Rời rạc hóa phi địa phương là yếu tố then chốt trong việc xây dựng lược đồ sai phân khác thường. Thay vì xấp xỉ các giá trị hàm tại một điểm lưới duy nhất, phương pháp này sử dụng thông tin từ các điểm lân cận để tăng cường độ chính xác và tính ổn định. Điều này giúp bảo toàn các tính chất quan trọng của nghiệm, chẳng hạn như tính đơn điệu và tính bảo toàn năng lượng, đặc biệt quan trọng trong việc mô phỏng các hệ vật lý.

IV. Ứng Dụng Giải Phương Trình Vi Phân Với Lược Đồ Khác Thường

Lược đồ sai phân khác thường có thể được áp dụng để giải một loạt các phương trình vi phân, bao gồm cả phương trình đạo hàm riêng và phương trình đạo hàm thường. Các lược đồ này đặc biệt hữu ích khi giải các bài toán biên và bài toán giá trị ban đầu, nơi các phương pháp truyền thống có thể gặp khó khăn.Ưu thế của các lược đồ khác thường so với lược đồ bình thường là bảo toàn tính chất nghiệm của bài toán với mọi cỡ bước h > 0.

4.1. Ví Dụ Về Các Phương Trình Có Thể Giải Bằng Phương Pháp Này

Lược đồ sai phân khác thường có thể được sử dụng để giải các phương trình phân rã tuyến tính, phương trình Logistic, và các hệ động lực học khác. Trong mỗi trường hợp, các lược đồ được xây dựng để bảo toàn các tính chất quan trọng của nghiệm, chẳng hạn như tính ổn định, tính đơn điệu, và tính bảo toàn năng lượng.

4.2. Thử Nghiệm Số Đánh Giá Tính Hiệu Quả Của Lược Đồ

Để đánh giá tính hiệu quả của lược đồ sai phân khác thường, các thử nghiệm số được thực hiện để so sánh hiệu suất của chúng với các phương pháp sai phân hữu hạn truyền thống. Các thử nghiệm này cho thấy rằng, lược đồ sai phân khác thường có thể cung cấp các giải pháp chính xác và ổn định hơn, đặc biệt khi giải các bài toán trên miền thời gian lớn.

V. Kết Luận Ưu Nhược Điểm và Hướng Phát Triển Lược Đồ

Lược đồ sai phân khác thường cung cấp một phương pháp mạnh mẽ để giải phương trình vi phân, đặc biệt là khi các phương pháp truyền thống gặp khó khăn trong việc bảo toàn các tính chất quan trọng của nghiệm. Mặc dù có một số hạn chế, như độ phức tạp cao hơn và khả năng giảm độ chính xác, nhưng các ưu điểm của lược đồ sai phân khác thường khiến chúng trở thành một công cụ có giá trị cho các nhà toán học, kỹ sư và nhà khoa học. Nhược điểm của các lược đồ khác thường là khó có thể đưa ra các lược đồ có cấp chính xác cao như các lược đồ bình thường và thời gian thực hiện tính toán có thể lâu hơn.

5.1. Tóm Tắt Ưu Nhược Điểm Của Lược Đồ Sai Phân

Ưu điểm chính của lược đồ sai phân khác thường là khả năng bảo toàn các tính chất quan trọng của nghiệm, chẳng hạn như tính ổn định, tính đơn điệu, và tính bảo toàn năng lượng. Nhược điểm chính là độ phức tạp cao hơn và khả năng giảm độ chính xác so với các phương pháp sai phân hữu hạn truyền thống.

5.2. Hướng Nghiên Cứu và Phát Triển Lược Đồ Sai Phân Trong Tương Lai

Hướng nghiên cứu trong tương lai bao gồm việc phát triển các lược đồ sai phân khác thường hiệu quả hơn về mặt tính toán, cũng như việc mở rộng các phương pháp này để giải các loại phương trình vi phân phức tạp hơn. Đặc biệt, việc tích hợp các kỹ thuật phân tích Fourier và các phương pháp dựa trên học máy có thể mở ra những khả năng mới cho việc thiết kế và ứng dụng lược đồ sai phân khác thường.

Trích đoạn nội dung tài liệu

Chương 1 Lược đồ sai phân khác thường 1.1 Một số kiến thức chuẩn bị Trong phần trình bày của luận văn, ta chủ yếu nghiên cứu việc giải gần đúng bài toán giá trị ban đầu đối với phương trình vi phân cấp một, hay còn gọi là bài toán Cauchy dy  Dy =  = f (t, y), t0 ≤ t ≤ T, dt (1.1)  y(t0 ) = y0 , y, f ∈ Rn ,  trong đó hàm y(t) : [t0 , T ] → Rn là hàm số cần xác định, giá trị ban đầu y0 ∈ Rn và hàm vế phải f : [t0 , T ] × Rn → Rn cho trước. Ta giả thiết rằng thời gian ban đầu t0 là hữu hạn, nhưng thời gian T có thể tiến đến vô cùng đối với hệ động lực học. Để đơn giản, ta giả sử rằng t0 = 0. Trong trường hợp f = f (y) thì phương trình được gọi là dừng (au- tonomous).

Không mất tính tổng quát ta có thể giả thiết một phương trình là dừng. Vì nếu phương trình không ở dạng dừng thì ta đưa thêm biến phụ yn+1 = t và đặt ŷ = (y1 , y2 ,. Khi đó phương trình được viết lại dưới dạng
T ŷ 0 = fˆ(ŷ), fˆ(ŷ) = f (y), 1 .2) Các kết quả liên quan đến bài toán giá trị ban đầu (1.1) như sự tồn tại và duy nhất nghiệm, sự phụ thuộc liên tục của nghiệm vào dữ liệu ban đầu. được trình bày trong hầu hết các giáo trình về phương trình vi phân (xem [3, 9, 10]) nên chúng ta không trình bày lại ở đây.

Từ giờ cho tới hết phần (LUAN.phan 8 TIEU LUAN MOI download : [e-mail] trình bày của luận văn, chúng ta luôn giả thiết rằng nghiệm của bài toán Cauchy (1.1) là duy nhất toàn cục trên [0, T ]. Điểm ȳ ∈ Rn được gọi là điểm bất động (fixed point) hay điểm cân bằng (equilibrium point) hoặc nghiệm hằng (constant solution) của phương trình (1. Giả sử ȳ là một điểm bất động (1. Điểm (vị trí) ổn định (stable) nếu với mọi  > 0, tồn tại δ = δ() > 0 sao cho nếu y0 ∈ B(ȳ, δ) thì nghiệm y(t, y0 ) ∈ B(ȳ, ) với mọi t ≥ 0.

Điểm ổn định tiệm cận (asymptotically stable) nếu y là ổn định và ||y(t, y0 ) − y|| → 0 khi t → ∞ với mọi ||y0 − y|| đủ nhỏ. Điểm không ổn định nếu điều kiện 1 của định nghĩa không được thỏa mãn. Từ định nghĩa trên ta thấy rằng: Nếu ȳ là ổn định thì các lời giải với giá trị ban đầu đủ gần ȳ cũng sẽ không nằm ngoài hình cầu tâm ȳ với bán kính . Còn nếu ȳ là ổn định tiệm cận thì nó sẽ hút các lời giải với giá trị ban đầu gần nó.

Tức là các lời giải với giá trị ban đầu gần ȳ sẽ hội tụ về ȳ. Hiển nhiên, nếu ȳ là ổn định tiệm cận thì nó là ổn định, trong trường hợp ȳ là ổn định nhưng không ổn định tiệm cận thì ta nói ȳ là ổn định yếu (marginally stable). Giả sử rằng f : Rn → Rn là hàm thuộc lớp hàm C 1 , ta ký hiệu Jacobian của hàm f tại điểm bất động ȳ là J ≡ Jf (ȳ). Nếu ma trận J có tất cả các giá trị riêng λk đều có phần thực Re(λk ) 6= 0 thì ta gọi ȳ là điểm bất động hyperbolic.

Ngược lại, ta gọi ȳ là điểm bất động không hyperbolic (non - hyperbolic). Giả sử f là hàm số khả vi liên tục với ȳ là điểm bất động hyperbolic và y(t) là một nghiệm của (1. Ta ký hiệu (t) = y(t) − y là thay đổi (nhiễu) của nghiệm y(t) xung quanh điểm bất (LUAN.phan 9 TIEU LUAN MOI download : [e-mail] động hyperbolic ȳ. Sử dụng khai triển Taylor, bỏ qua thành phần bậc hai trở đi ta có 0 = y 0 (t) = f (y) = f (ȳ + ) ≈ f (ȳ) + Jf (ȳ) = J.

Phương trình 0 = J, (1.3) được gọi là phương trình tuyến tính hóa của phương trình nhiễu xung quanh điểm ổn định hyperbolic ȳ. Giả sử f là hàm số khả vi liên tục và ȳ là một điểm bất động hyperbolic. Khi đó ȳ là ổn định tiệm cận khi và chỉ khi đối với các trạng thái ban đầu 0 = y0 − y mà ||0 || đủ nhỏ thì nghiệm (t) = 0 etJ của (1. t→∞ Điều này tương đương với Re(λ) < 0, ∀λ ∈ σ(J), trong đó σ(J) là tập hợp các giá trị riêng của ma trận J.

Ngược lại, điểm bất động y là không ổn định khi và chỉ khi tồn tại λ ∈ σ(J) sao cho Re(λ) > 0 hoặc limt→∞ (t) = ∞. Định lý không áp dụng được trong trường hợp ȳ không phải điểm bất động hyperbolic. Trong trường hợp Re(λ) < 0, ∀λ ∈ σ(J), ta còn gọi điểm ȳ là điểm ổn định tuyến tính (linearly stable). Điều này tương đương sự kiện phương trình tuyến tính 0 = J là ổn định.

Xét hệ động lực xác định bởi phương trình Logistic y 0 = λy(1 − y), y(0) = y0 .4) Phương trình này có nghiệm chính xác y0 y(t) = .5) (1 − y0 )e−λt + y0 Phương trình (1.4) có hai điểm bất động là y = 0, ŷ = 1. Sử dụng Định lý 1.1, ta nhận được.phan 10 TIEU LUAN MOI download : skknchat@gmail. ȳ = 0 là ổn định tiệm cận nếu λ < 0 và không ổn định nếu λ > 0. ŷ = 1 là không ổn định nếu λ < 0 và ổn định tiệm cận nếu λ > 0.

λ = 0 thì mọi hằng số đều là ổn định nhưng không ổn định hyperbolic. Trong trường hợp λ < 0, tất cả các nghiệm của (1.4) với y0 < 1 đều hội tụ đơn điệu về ȳ = 0. Tất cả các nghiệm với y0 > 1 dẫn đến các quỹ đạo tăng trưởng không bị chặn trong thời gian hữu hạn. Tức là, các nghiệm bùng nổ (blow - up solution).

Thời gian bùng nổ là 1 1 y0 y0 − 1 T (y0 ; λ) = − ln = ln( ) λ > 0.6) λ y0 − 1 y0 Nếu λ > 0, tất cả các nghiệm với y0 > 0 đều hội tụ đơn điệu về ŷ = 1. Tất cả các nghiệm với y0 < 0 dẫn đến các quỹ đạo tiến về −∞ trong khoảng thời gian hữu hạn T , trong đó 1 1 y0 − 1 y0 − 1 T (y0 ; λ) = ln = ln( ) λ > 0.4 biểu diễn các quỹ đạo của phương trình Logistic. Tiếp theo, ta trình bày lại một số khái niệm và kết quả liên quan tới các phương pháp số giải phương trình vi phân. Chúng ta xét bài toán giá trị ban đầu (1.1)  0 y = f (t, y), 0 ≤ t ≤ T, y(0) = y0 , y, f ∈ Rn ,  với giả thiết hàm f đủ trơn có các đạo hàm riêng bị chặn sao cho lời giải của bài toán là tồn tại duy nhất và lời giải có đạo hàm bị chặn tới cấp cần thiết.

Ta lấy một phân hoạch không nhất thiết đều của đoạn [0, T ] (rời rạc hóa trục thời gian) π = {0 = t0 < t1 < t2 <. < tN = T }, hn = tn − tn−1 được gọi là bước lưới hay cỡ bước (stepsize). Ta cần tìm các xấp xỉ yn ≈ y(tn ) tại các nút lưới t0 ,. Ta gọi các giá trị {yn }n=N n=0 là nghiệm số xấp xỉ, hay ngắn gọn là nghiệm số của phương pháp trên lưới π .phan 11 TIEU LUAN MOI download : [e-mail] 12 TIEU LUAN MOI download : [e-mail] Ta xét lược đồ số một bước ở dạng Dh (yk ) = Fh (f ; yk ), (1.8) dy trong đó Dh ≈ là phần rời rạc hóa đạo hàm trên lưới.

Còn Fh là phần dt rời rạc hóa hàm vế phải. Tương tự như đối với phương trình vi phân, ta có kết quả phát biểu cho phương trình sai phân (1. Giả sử ȳ là một điểm bất động của lược đồ (1. Ký hiệu Jh là Jacobian của hàm F tại ȳ.

Khi đó, nhiễu k = yk − ȳ xung quanh điểm bất động y được xấp xỉ bằng nghiệm phương trình tuyến tính hóa (sử dụng khai triển Taylor) Dh k = Jh k .9) Điểm bất động y được gọi là ổn định tuyến tính khi và chỉ khi nghiệm {k } (với ||0 || đủ nhỏ) của (1. Điều này tương đương với điều kiện mọi giá trị riêng λi của Jh đều có |λi | < 1. Trong rất nhiều các giáo trình về phương trình sai phân, Định lý 1.2 được lấy làm định nghĩa điểm bất động ổn định tuyến tính cho phương trình sai phân. Các phương pháp sai phân bình thường chủ yếu được xây dựng dựa trên việc rời rạc hóa đạo hàm bằng các công thức sai phân hữu hạn.

Chẳng hạn, công thức sai phân tiến, sai phân lùi, sai phân trung tâm yn+1 − yn yn − yn−1 yn+1 − yn−1 y 0 (tn ) ≈ , y 0 (tn ) ≈ , y 0 (tn ) ≈. h h 2h Tương ứng ta nhận được: Công thức Euler hiển yn+1 = yn + hf (tn , yn ).10) Công thức Euler ẩn yn+1 = yn + hf (tn+1 , yn+1 ).11) Công thức hình thang ẩn hn   yn+1 = yn + f (tn , yn ) + f (tn+1 , yn+1 ) .phan 13 TIEU LUAN MOI download : [e-mail] Các công thức Euler hiển, Euler ẩn, hình thang ẩn là trường hợp riêng của công thức θ có dạng hn   yn+1 = yn + θf (tn , yn ) + (1 − θ)f (tn+1 , yn+1 ). 2 Đây là trường hợp riêng của phương pháp Runge - Kutta tổng quát. Phương pháp Runge - Kutta và phương pháp Taylor là các phương pháp tiêu biểu cho các phương pháp sai phân bình thường giải phương trình vi phân.

Chúng ta đã biết rằng các phương pháp trên đều là tương thích, hội tụ (xem [3, 9, 10]). Điều đó có nghĩa là khi bước lưới h dần về 0 thì sai số toàn cục |y(tn ) − yn | cũng hội tụ về 0. Tuy nhiên, nếu chúng ta chỉ quan tâm đến sự tương thích, tính hội tụ của phương pháp thôi thì chưa đủ. Trong nhiều trường hợp, các phương pháp cần bảo toàn được các tính chất của bài toán.

Chẳng hạn, tính dương, tính bị chặn, tính đơn điệu, tính tuần hoàn, tính chất ổn định của các điểm bất động, các tính chất bất biến trong các phương trình mô tả các hiện tượng Vật lý, Sinh học. Để minh họa cho điều này, chúng ta xét ví dụ đơn giản sau đây.