Luận văn thạc sĩ: Phương pháp Runge-Kutta giải phương trình vi phân ngẫu nhiên

Luận văn thạc sĩ toán học nghiên cứu toán học phương pháp runge kutta giải phương trình vi phân ngẫu nhiên, khảo sát thực trạng, phân tích nguyên nhân, đề xuất giải pháp cải thiện

• Trường học: Trường Đại Học Khoa Học Tự Nhiên
• Chuyên ngành: Toán Học
• Đề tài: Luận Văn Thạc Sĩ Toán Học Phương Pháp Runge Kutta Giải Phương Trình Vi Phân Ngẫu Nhiên
• Loại tài liệu: Luận Văn
• Năm xuất bản: 2023
• Địa điểm: Hà Nội

Mục lục chi tiết

• MỞ ĐẦU
• 1. CHƯƠNG 1: LÝ THUYẾT XÁC SUẤT
• 1.1. ø-Đại số Borel trên R*
• 1.2. Không gian xác suất
• 1.3. Biến ngẫu nhiên
• 1.4. Các tham số đặc trưng của biến ngẫu nhiên
• 1.4.1. Kỳ vọng toán
• 1.4.2. Phương sai
• 1.5. Vectơ ngẫu nhiên
• 1.6. Hiệp phương sai, hệ số tương quan
• 1.7. Quá trình ngẫu nhiên
• 2. CHƯƠNG 2: PHƯƠNG TRÌNH VI PHÂN NGẪU NHIÊN VÀ PHƯƠNG PHÁP RUNGE - KUTTA
• 2.1. Lý thuyết phương trình vi phân ngẫu nhiên
• 2.1.1. Khai triển Taylor ngẫu nhiên
• 2.1.2. Nghiệm số phương trình vi phân ngẫu nhiên
• 2.1.3. Phương pháp Runge - Kutta
• 2.2. Tính hội tụ của sơ đồ yếu
• 2.2.1. Sự xấp xỉ yếu địa phương và toàn cục
• 2.3. Lý thuyết cây có gốc
• 2.3.1. Các cây ngẫu nhiên
• 2.4. Hệ phương trình vi phân ngẫu nhiên Itô (hệ SDE Itô)
• 2.5. Hệ SDE Stratonovich
• 2.6. Phương pháp ngẫu nhiên Runge - Kutta
• 2.6.1. Sự khai triển của phương pháp Runge- Kutta ngẫu nhiên
• 2.6.2. Các điều kiện bậc tổng quát
• 2.7. Thuật toán và ví dụ
• KẾT LUẬN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Luận văn Runge Kutta giải phương trình vi phân ngẫu nhiên

Luận văn thạc sĩ toán học về phương pháp Runge-Kutta giải phương trình vi phân ngẫu nhiên (SDE) là một công trình nghiên cứu phương trình vi phân chuyên sâu, đóng góp quan trọng cho lĩnh vực giải tích số. Các phương trình vi phân ngẫu nhiên (SDEs) đóng vai trò nền tảng trong việc mô tả các hệ thống biến đổi theo thời gian và chịu ảnh hưởng của các yếu tố nhiễu loạn ngẫu nhiên. Chúng xuất hiện phổ biến trong nhiều ngành khoa học kỹ thuật, vật lý, sinh học và đặc biệt là tài chính định lượng. Tuy nhiên, việc tìm ra nghiệm giải tích tường minh cho hầu hết các SDE là bất khả thi, đòi hỏi sự phát triển của các phương pháp số giải SDE. Trong bối cảnh đó, phương pháp Runge-Kutta, vốn rất thành công với phương trình vi phân thường (ODE), được mở rộng và hiệu chỉnh để áp dụng cho SDE. Công trình nghiên cứu này tập trung vào việc hệ thống hóa kiến thức cơ bản về giải tích ngẫu nhiên, từ tích phân Itô đến khai triển Taylor ngẫu nhiên, làm cơ sở để xây dựng công thức Runge-Kutta tổng quát. Mục tiêu của luận văn là cung cấp một tài liệu tham khảo có giá trị, vừa mang ý nghĩa lý thuyết, vừa có tính ứng dụng thực tiễn cao, đặc biệt hữu ích cho sinh viên và nhà nghiên cứu trong lĩnh vực luận văn toán ứng dụng.

1.1. Tầm quan trọng của phương trình vi phân ngẫu nhiên SDE

Phương trình vi phân ngẫu nhiên là công cụ mô hình hóa toán học không thể thiếu khi nghiên cứu các hệ thống có sự bất định. Khác với phương trình vi phân thường (ODE) mô tả các quá trình tất định, SDE tích hợp thêm thành phần ngẫu nhiên, thường được biểu diễn qua chuyển động Brown hay quá trình Wiener. Điều này cho phép mô tả chính xác hơn các hiện tượng trong thực tế, từ sự biến động giá cổ phiếu trên thị trường tài chính, sự khuếch tán của các hạt trong vật lý, đến sự phát triển của quần thể trong sinh học. Ví dụ, mô hình Black-Scholes nổi tiếng trong tài chính là một SDE mô tả giá của một tài sản. Sự phức tạp do yếu tố ngẫu nhiên gây ra khiến việc giải các phương trình này trở thành một thách thức lớn, thúc đẩy sự phát triển mạnh mẽ của các phương pháp số và mô phỏng Monte Carlo.

1.2. Mục tiêu và cấu trúc của báo cáo khoa học này

Mục tiêu chính của luận văn là xây dựng một công thức Runge-Kutta tổng quát và tường minh để giải gần đúng SDE. Để đạt được mục tiêu này, công trình nghiên cứu được cấu trúc thành các phần rõ ràng. Phần đầu hệ thống hóa các kiến thức nền tảng về lý thuyết xác suất và giải tích ngẫu nhiên. Phần tiếp theo, cũng là nội dung trọng tâm, trình bày chi tiết cách xây dựng các sơ đồ Runge-Kutta ngẫu nhiên, phân tích các điều kiện bậc hội tụ, và sử dụng lý thuyết cây (rooted tree theory) để đơn giản hóa các khai triển phức tạp. Luận văn cũng đưa ra các ví dụ minh họa cụ thể cho từng trường hợp, làm rõ cách áp dụng thuật toán số vào thực tế. Đây là một báo cáo khoa học có cấu trúc chặt chẽ, hướng đến việc trở thành tài liệu tham khảo hữu ích cho các đối tượng quan tâm.

II. Thách thức trong giải tích số đối với phương trình vi phân SDE

Việc giải phương trình vi phân ngẫu nhiên đặt ra nhiều thách thức hơn đáng kể so với phương trình vi phân thường. Nguồn gốc của sự phức tạp này nằm ở quá trình Wiener (Wiener process), một thành phần không khả vi và có biến động vô hạn. Do đó, các phương pháp giải tích cổ điển không thể áp dụng trực tiếp. Các phương pháp số giải SDE phải được xây dựng trên nền tảng của giải tích ngẫu nhiên, cụ thể là tích phân Itô hoặc Stratonovich, vốn có các quy tắc tính toán khác biệt so với giải tích thông thường. Hơn nữa, việc đánh giá độ chính xác của một phương pháp số cho SDE cũng phức tạp hơn. Người ta phải phân biệt giữa hai tiêu chí hội tụ: hội tụ mạnh (strong convergence), đo lường sai số của từng quỹ đạo, và hội tụ yếu (weak convergence), đo lường sai số của các moment (như kỳ vọng). Các sơ đồ số ban đầu như sơ đồ Euler-Maruyama tuy đơn giản nhưng chỉ đạt bậc hội tụ thấp, không đáp ứng đủ yêu cầu về độ chính xác cho nhiều bài toán phức tạp, từ đó thúc đẩy nhu cầu tìm kiếm các phương pháp bậc cao hơn như Runge-Kutta.

2.1. Phân biệt tích phân Itô và tích phân Stratonovich

Trong giải tích ngẫu nhiên, có hai loại tích phân chính là Itô và Stratonovich, dẫn đến hai dạng SDE khác nhau. Tích phân Itô được định nghĩa sao cho hàm bị tích không đoán trước được tương lai của quá trình Wiener, làm cho nó trở thành một martingale và phù hợp với nhiều mô hình tài chính nơi thông tin quá khứ không dự báo được biến động tương lai. Ngược lại, tích phân Stratonovich tuân theo các quy tắc giải tích chuỗi (chain rule) giống như giải tích thông thường, giúp nó phù hợp hơn trong các ứng dụng vật lý. Một SDE dạng Itô có thể được chuyển đổi thành dạng Stratonovich và ngược lại. Sự lựa chọn giữa hai loại tích phân này ảnh hưởng trực tiếp đến việc xây dựng và phân tích thuật toán số, vì các khai triển Taylor ngẫu nhiên tương ứng sẽ khác nhau.

2.2. Hạn chế của sơ đồ Euler Maruyama và phương pháp Milstein

Sơ đồ Euler-Maruyama là phương pháp số đơn giản nhất để giải SDE, được xem là phần mở rộng trực tiếp của phương pháp Euler cho ODE. Nó có bậc hội tụ mạnh là 0.5 và bậc hội tụ yếu là 1.0. Mặc dù dễ thực hiện, bậc hội tụ thấp khiến nó đòi hỏi kích thước bước rất nhỏ để đạt được độ chính xác mong muốn, dẫn đến chi phí tính toán cao. Phương pháp Milstein là một cải tiến bằng cách thêm vào một số hạng bậc cao hơn từ khai triển Itô-Taylor. Sơ đồ này đạt được bậc hội tụ mạnh là 1.0. Tuy nhiên, đối với các hệ SDE nhiều chiều, việc tính toán các số hạng bổ sung trong phương pháp Milstein trở nên rất phức tạp. Những hạn chế này chính là động lực để phát triển các phương pháp Runge-Kutta ngẫu nhiên, vốn có khả năng đạt được bậc hội tụ cao hơn mà không cần tính toán các đạo hàm phức tạp của các hệ số.

III. Hướng dẫn xây dựng phương pháp Runge Kutta giải SDE

Xây dựng phương pháp Runge-Kutta cho phương trình vi phân ngẫu nhiên là một quá trình kết hợp giữa lý thuyết Runge-Kutta cổ điển và các công cụ của giải tích ngẫu nhiên. Thay vì chỉ có một khai triển Taylor duy nhất như trong trường hợp ODE, SDE yêu cầu sử dụng khai triển Itô-Taylor, một công cụ phức tạp hơn bao gồm các tích phân Itô bội. Mục tiêu là xây dựng một sơ đồ lặp nhiều giai đoạn (multi-stage) để xấp xỉ nghiệm tại mỗi bước thời gian, sao cho khai triển của nghiệm số này trùng khớp với khai triển Itô-Taylor của nghiệm chính xác đến một bậc mong muốn. Các hệ số của phương pháp (tương tự bảng Butcher trong ODE) phải thỏa mãn một hệ các phương trình điều kiện bậc (order conditions). Việc rút ra các điều kiện này rất phức tạp và thường dựa vào lý thuyết cây có gốc (rooted tree theory) được mở rộng cho trường hợp ngẫu nhiên. Các cây này cung cấp một cách biểu diễn trực quan và hệ thống cho các vi phân cơ bản và các tích phân ngẫu nhiên xuất hiện trong khai triển, giúp việc xác định các điều kiện bậc trở nên khả thi.

3.1. Nền tảng từ khai triển Taylor ngẫu nhiên và lý thuyết cây

Khai triển Itô-Taylor là công cụ toán học trung tâm để phân tích sai số hội tụ của các phương pháp số cho SDE. Nó tổng quát hóa khai triển Taylor cổ điển bằng cách bao gồm các số hạng liên quan đến cả thời gian (dt) và quá trình Wiener (dW). Các số hạng này bao gồm các tích phân Itô bội, ví dụ như ∫dW và ∫∫dWdW. Để hệ thống hóa các số hạng phức tạp này, lý thuyết cây có gốc được sử dụng. Mỗi cây tương ứng với một vi phân cơ bản trong khai triển. Các nút và nhánh của cây biểu diễn các toán tử vi phân và các thành phần của SDE. Lý thuyết này cho phép biến đổi các bài toán giải tích phức tạp thành các bài toán tổ hợp trên cây, giúp đơn giản hóa đáng kể việc derivation các điều kiện bậc cho các phương pháp Runge-Kutta bậc cao.

3.2. Phân tích độ ổn định và sai số hội tụ của phương pháp

Cũng như trong trường hợp ODE, độ ổn định của phương pháp là một yếu tố quan trọng khi giải SDE. Một sơ đồ số phải ổn định để đảm bảo rằng sai số không bị khuếch đại một cách không kiểm soát khi lặp qua nhiều bước thời gian. Việc phân tích ổn định cho SDE phức tạp hơn do sự hiện diện của thành phần ngẫu nhiên. Các nhà nghiên cứu thường phân tích sự ổn định trung bình bình phương (mean-square stability). Bên cạnh đó, việc phân tích sai số hội tụ là cốt lõi để đánh giá hiệu quả của phương pháp. Luận văn này tập trung vào việc thiết lập các điều kiện bậc để đảm bảo sơ đồ Runge-Kutta đạt được bậc hội tụ yếu (weak order) hoặc mạnh (strong order) mong muốn. Việc đạt được bậc hội tụ cao hơn cho phép sử dụng kích thước bước lớn hơn, giảm đáng kể thời gian tính toán.

IV. Thuật toán số và triển khai phương pháp Runge Kutta Maruyama

Từ nền tảng lý thuyết, luận văn tiến hành xây dựng các thuật toán số cụ thể. Phương pháp Runge-Kutta-Maruyama, một trong những biến thể phổ biến, được trình bày chi tiết. Cấu trúc của một phương pháp Runge-Kutta ngẫu nhiên (SRK) thường được biểu diễn qua một bảng Butcher mở rộng, bao gồm các ma trận hệ số không chỉ cho thành phần drift (tất định) mà còn cho thành phần diffusion (ngẫu nhiên). Thuật toán bắt đầu với một giá trị ban đầu, và tại mỗi bước, nó tính toán một loạt các giá trị trung gian (stages) dựa trên các giá trị trước đó và các số gia của quá trình Wiener. Cuối cùng, một tổ hợp tuyến tính của các giá trị trung gian này được sử dụng để cập nhật nghiệm. Việc triển khai các thuật toán này đòi hỏi phải sinh ra các biến ngẫu nhiên tuân theo phân phối chuẩn để mô phỏng các số gia của quá trình Wiener. Các ngôn ngữ lập trình như lập trình Matlab giải SDE hoặc code Python cho SDE với các thư viện như NumPy và SciPy là những công cụ mạnh mẽ để thực hiện các mô phỏng này.

4.1. Cấu trúc bảng Butcher mở rộng cho phương pháp SRK

Phương pháp Runge-Kutta ngẫu nhiên hiện (explicit Stochastic Runge-Kutta) giai đoạn s có thể được mô tả gọn gàng bằng một bộ các ma trận và vector hệ số, thường được gọi là bảng Butcher mở rộng. Bảng này không chỉ chứa các hệ số A, c, b như trong phương pháp RK cổ điển cho thành phần drift, mà còn bổ sung các hệ số cho thành phần diffusion. Ví dụ, phương pháp Euler-Maruyama có thể được xem là phương pháp SRK đơn giản nhất với một giai đoạn. Các phương pháp bậc cao hơn, như Runge-Kutta-Maruyama, sẽ có nhiều giai đoạn và các bảng hệ số phức tạp hơn, được thiết kế cẩn thận để thỏa mãn các điều kiện bậc hội tụ đã được xác định trước đó thông qua lý thuyết cây.

4.2. Các bước triển khai code Python hoặc Matlab giải SDE

Việc triển khai code Python cho SDE hoặc lập trình Matlab giải SDE theo phương pháp Runge-Kutta bao gồm các bước chính sau: 1) Định nghĩa các hàm drift a(t,X) và diffusion b(t,X) của SDE. 2) Khởi tạo các tham số mô phỏng: khoảng thời gian, điều kiện ban đầu, số bước thời gian (N), và kích thước bước h. 3) Tại mỗi bước lặp từ 1 đến N, sinh ra các số gia ngẫu nhiên dW từ phân phối chuẩn N(0,h). 4) Tính toán các giá trị trung gian (stages) theo công thức của phương pháp Runge-Kutta đã chọn. 5) Cập nhật giá trị nghiệm X_next bằng cách tổ hợp các stages. Quá trình này được lặp lại nhiều lần (ví dụ, trong mô phỏng Monte Carlo) để thu được một tập hợp các quỹ đạo, từ đó có thể tính toán các đại lượng thống kê như giá trị kỳ vọng hoặc phương sai.

V. Top ứng dụng SDE trong tài chính và mô hình hóa toán học

Các phương pháp số để giải SDE, đặc biệt là phương pháp Runge-Kutta, có vô số ứng dụng thực tiễn. Lĩnh vực nổi bật nhất là tài chính định lượng, nơi SDE là công cụ tiêu chuẩn để mô hình hóa toán học sự biến động của giá tài sản, lãi suất, và các công cụ phái sinh. Mô hình Black-Scholes, mô hình Heston (với biến động ngẫu nhiên), hay các mô hình lãi suất như Vasicek và Cox-Ingersoll-Ross (CIR) đều là các SDE. Khi các mô hình này không có nghiệm giải tích (ví dụ, khi định giá các quyền chọn kiểu Mỹ hoặc các quyền chọn exotic), các phương pháp số trở nên tối quan trọng. Phương pháp Runge-Kutta cho phép mô phỏng các quỹ đạo giá tài sản với độ chính xác cao, sau đó sử dụng mô phỏng Monte Carlo để tính giá trị kỳ vọng của khoản thanh toán (payoff), từ đó định giá quyền chọn. Ngoài tài chính, các phương pháp này còn được sử dụng trong kỹ thuật (ví dụ, trong lý thuyết lọc và điều khiển), vật lý (mô phỏng các hệ động lực học phân tử), và sinh học (mô hình hóa động lực học quần thể).

5.1. Ứng dụng SDE trong tài chính Mô hình Black Scholes

Một trong những ứng dụng SDE trong tài chính kinh điển nhất là mô hình Black-Scholes. Mô hình này mô tả giá cổ phiếu S(t) tuân theo một quá trình chuyển động Brown hình học: dS = μS dt + σS dW. Trong trường hợp quyền chọn châu Âu đơn giản, mô hình này có một công thức định giá dạng đóng. Tuy nhiên, đối với các quyền chọn phức tạp hơn hoặc khi các giả định của mô hình được nới lỏng (ví dụ, lãi suất hoặc biến động thay đổi theo thời gian), nghiệm giải tích không còn tồn tại. Lúc này, việc sử dụng các phương pháp số như Runge-Kutta để mô phỏng các đường đi của giá cổ phiếu S(t) là cần thiết. Bằng cách mô phỏng hàng ngàn hoặc hàng triệu quỹ đạo, ta có thể tính toán giá trị trung bình của khoản thanh toán quyền chọn và chiết khấu về hiện tại để có được giá của nó.

5.2. Kỹ thuật mô phỏng Monte Carlo và vai trò của phương pháp số

Mô phỏng Monte Carlo là một kỹ thuật tính toán dựa trên việc lặp lại các thử nghiệm ngẫu nhiên để thu được kết quả số. Trong tài chính, nó được sử dụng rộng rãi để định giá các công cụ phái sinh. Cốt lõi của kỹ thuật này là khả năng sinh ra các quỹ đạo (path generation) của các biến ngẫu nhiên cơ sở (như giá cổ phiếu). Đây chính là lúc các phương pháp số giải SDE phát huy vai trò. Một sơ đồ số như Runge-Kutta chính là công cụ để rời rạc hóa SDE và tạo ra một quỹ đạo giá gần đúng tại các điểm thời gian t_0, t_1, ..., t_N. Độ chính xác của toàn bộ mô phỏng Monte Carlo phụ thuộc trực tiếp vào độ chính xác của sơ đồ số được sử dụng để tạo ra mỗi quỹ đạo. Một phương pháp có bậc hội tụ yếu cao hơn sẽ cho kết quả chính xác hơn với cùng một số lượng quỹ đạo, hoặc đạt cùng độ chính xác với số lượng quỹ đạo ít hơn, giúp tiết kiệm thời gian tính toán.

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

Chương 1: Kiến thức chuẩn bị. Chương này trình bày sơ lược một số khái niệm và vấn đề liên quan đến xác suất và phương trình vi phân ngẫu nhiên để làm cơ sở cho. chương sau. © Chương 9: Phương pháp Runge-Kutta gidi phương trình ơi phân ngẫu nhiên.

Chương, này là nội dung chính của luận văn, trình bày cách giải phương trình vi phân ngẫu lên bằng phương pháp Runge-Kutta. Chuong 1 LÝ THUYẾT XÁC SUAT 1. ø-Đại số Định nghĩa 1L „ Cho © là một t ip hợp khác rỗng, một ø-đại số F trên © là một họ. những tập con ct với các tính cÌ sau: i) OCF, Néu Fe F thi F°=0\Fe F, iii) Nếu Ai,4a,.

€ Ƒ thì UỆ Ai € F. ø-Đại số Borel trên R* Cho .A là một họ các tập con của 2. Tập hợp tắt cả các tập con của 9 là một ø-đại y chỉ ra rằng tồn tại ít nhất một ø-đại số chứa. Hơn nữa, một họ bất kì một ø-đại số cũng là một ø-đại số, do đó tồn tại ø-đại số nhỏ nhất chứa.

Vậy ø-đại số nhỏ nhất chứa. ø- Dai sé Borel trén R*, kí hiệu bởi %(R*), được xác định là ø- đại số nhỏ nhất chứa tất cả các tập mở trong RẺ. Không gian xác suất Định nghĩa 1. Một không gian xác suất là bộ (Q,Z,P), trong đó 9 là một tập hợp khác rồng, Z là một o- đại số trên @, P : @ -› R là một độ đo xác suất trên Z, tức hà: 3) 0< P(A) < 1 với mọi A € Ƒ, P(9)=1, iii) Voi mọi Ai, 4a.

Biến ngẫu nhiên Định nghĩa 1. Cho (Q,Z,P) là một không gian xác suất. Một biến ngẫu nhiên nhận giá trị thực là bất kì ánh xạ Borel đo được X : 9 -› R tức là với mỗi tap Borel Bc9%(R) thì X~Ì(B) e Z. Cho X : © -› R là một biến ngẫu nhiên, ánh xạ Px :®(R) +R với Px(B) = P(X~!(B)) = P([X € Bì) là một độ đo xác suất trên lR và nó được gọi là quy luật xác suất của X.

Cho P : B(R) > R la mét dé do zác suất. Tồn tại một biến ngẫu nhiên X :R — R sao cho P trùng tới luật xác suất Px kết hợp với X. Cho X : 9 — E là một biến ngẫu nhiên. Khi đó ø- đại số Fy = X~!(B(R)) được gọi là đ— đại số sinh bởi X.

Cho X : © —› R là một biến ngẫu nhiên. Ánh xạ Fx :R — [0:1] với Fx(z) = P(X < z),z € lR, được gọi là hàm tích lũy của X. 3) Với mọi a,b € R,a < b: Fx(a) — Fx(b) = P{a < X < 9). ñ) Fx(z) là hàm tăng tà liên tục phải.

ii) lim; ,_„ Fx(z) =0 tà lim, „¿. Fx(e) =L Mệnh dé 1. Cho một hàm F : R — [0:1] thỏa mãn ii) và iii) của mệnh đề 1.2 thì bằng nhận xét 1) ta có thể xác định một độ đo xác sudt Px : B(R) > R lién két uới một. biến ngẫu nhiên X mà hàm phân phối tích lũy là đồng nhất trên F.

Một biến ngẫu nhiên X là liên tục nếu hàm phân phối tích lũy của Fx cita nó là liên tục. Một biến ngẫu nhiên X là liên tục nếu và chỉ nếu P(X = = 0 với mọi a €R. Cho X là một biến ngẫu nhiên liên tục. Nếu tồn tại một hàm f(r) Khong am F : R — (0;—co] véi moi x € R ma: Fx() = J", f(t)dt.xz thì ƒ(z) được gọi là hàm mật độ xác suất của biến ngẫu nhị X.

Một biến ngẫu nhiên X được gọi là rời rạc nếu miền giá trị của nó là tập hữu hạn hoặc tập vô hạn đếm được. Nếu X() là miền giá trị của biến ngẫu nhiên rồi rac X thi P(X =2), nếu z€ X(0) P(x) = 0, 4 nếuzýX(0) được gọi là hàm mật độ xác suất của biến ngẫu nhiên X. Biến ngẫu nhiên liên tục X được gọi là biến ngẫu nhiên có phân bố chuẩn nếu X có hàm mật độ xác suất ƒ@)=—=cVino “# ren, trong dé p € R vag > 0. Các tham số đặc trưng của biến ngẫu nhiên 1.

Kỳ vọng toán Định nghĩa 1. Kỳ vọng hoặc giá trị trung bình của biến ngẫu nhiên X (kí hiệu E(X)) được xác định bởi E(X) =I, XaP, trong đó tích phân về phải là tích phan Lebesgue. Dac bit i) nếu X rời rạc và nhận các giá trị z¡ với xác suất tương ứng ø = P(X = z¡) thì E(X) = XI zin: ii) nếu X liên tục và có hàm mật độ ƒ(z) thì E(X) = ƒ*Š zƒ(z)dz Tinh chat 1. i) E(C) =C vdi moi hing số C.

ii) E(CX) = CE(X) vai mọi hằng số C. iv) Nếu Xị,X¿.X„ độc lập thì E(Xì. Phương sai Định nghĩa 1. Phương sai của biến ngẫu nhiên X (kí hiệu Var(X)) là đại lượng, đo sự phân tán bình phương trung bình của X xung quanh giá trị trung bình E(X): Var(X) = E(X - EX)? Nếu X rời rạc nhận các giá trị z; với các xác suất tương ứng p; thi Var(X) = Š”(z¡ — ELX))°p, 6 ii) Nếu X liên tục và có hàm mật độ ƒ(z) thì +00 Var(X) -f[ (2 - E(X))?f(w)dr.

'ar(CX) = C?Var(X) vdi moi hang sé C. Var(C\X + Ca) = C?Var(X) vdi moi hing s6 Cy, C2 iii) Nếu Xị,. X„ độc lập và có phương sai hữu hạn thì Var(CiX) +. Vectơ ngẫu nhiên Định nghĩa 1.

Một vectơ ngẫu nhiên ø chiều là một bộ có thứ tự (Xạ, Xa,. Xu) với các thành phần Xị, Xạ. X„ là các biến ngẫu nhiên. Cho vectơ ngẫu nhiên X = (Xị, X;.

Xét hàm số Fx : R" — R xác định bởi: EX(ri,ta,.en) = P(X < 21, X2 < 22, <zm) với mọi (7ị,#a,. #u) € IR", Hàm Fx(z,z2.z„) được gọi là hàm phân bố xác suất của vectơ ngẫu nhiên X hay còn được gọi là hàm phân bố đồng thời của các biến ngẫu nhiên Xị, Xạ,. Xu) là vectơ ngẫu nhiên. Nếu tồn tại một hàm số ƒ(zạ,#a,.

zạ) > 0 với mọi z = (z,.dtn thì /(z) được gọi là hàm mật độ đồng thời của vectơ ngẫu nhiên X. Vectơ ngẫu nhiên X =(X\, Xs,. Xu)” có giá trị trong R" được gọi là có phân bố chuẩn n chiều nếu với mọi e = (et.en)T € R", Y =cTX= clẨ + cạẤz +. + cụX„ có phân bố chuẩn 1 chiều.

Hiệp phương sai, hệ số tương quan Định nghĩa 1. Cho vectơ ngẫu nhiên(X, Y). Hiệp phương sai của X và Y là một số xác định bởi công thức: Cou(X,¥) = E[(X ~ E(X)(Y = E(Y))) Định lý 1. 3) Nếu X nà Y độc lập thì Cou(X,Y) = 0.

1) NéwX waV là các biến ngẫu nhiên rời rạc có miền giá trị lần lượt là {z\, Jy Yn} thì Con(X,Y) =3 ^^ zup(si.w) — E(X)E(V) mí Trong dé: pi; = P(X = 2;,¥ = yj) 3) Nếu (X,Y) có hàm mật độ zác suất đồng thời ƒ(z,y) thì: dc poe €m(x. Hệ số tương quan của hai biến ngẫu nhiên X và Y, kí hiệu là Ø(X. được xác định bởi công thức: (X,Y) = Cov(X.Y)) Var(X) /Var(Y) | Dinh ly 1. 2) Nếu X tà Y độc lập thì p(X,Y) = 8) Nếu Y = aX +b thì: 1 nếua<0 ø(X.

Quá trình ngẫu nhiên Định nghĩa 1. Cho (©,Z,P) là một không gian xác suất, 7 là tập chỉ số thực, TCR.T có thể hữu hạn đếm được hoặc vô hạn không đếm được. Một quá trình ngẫu nhiên trên (0, Z, P) là một họ các tập (Xu;† € T) với X, : 9 —> R là một biến ngẫu nhiên với mọi £ € T. Cho £ € 7, biến ngẫu nhiên X, được gọi là quá trình tại thời gian t.

Hơn nữa, với mọi œ€ 9 cố định.;w) : T+ R với X(f;«) = X,(œ) được gọi là quỹ đạo tương, ứng với yếu tố ngẫu nhiên w.£ > 0) của Z, Z¿ C Z được gọi là một bộ lọc nếu thỏa. mãn các điều kiện sau: j) Nếu s < £ thì #, C Z¡, ii) Fi=N aso Fitz ii) Nếu A€ 7 và P(4) = 0 thì A € Zụ. Cho quá trình ngẫu nhiên X = (X,„£ > 0). Xét ø- trường sinh bởi các biến ngẫu nhiên X, với s < +: #Ä = ø(X, : s < #) (ZXŸ,t > 0) gọi là bộ lọc tự nhiên của quá trình X hay lịch sử của X.

Trường này chứa thông tin về diễn biến của quá trình X cho đến thời điểm ¿. Cho một không gian xác suất (0, Z, P) thêm vào một bộ lọc (Fi) gọi là một không gian xác suất được lọc và kí hiệu là (Q,Z, (7), P). Quá trình ngẫu nhiên Y = Yj,£ > 0) được gọi là thích nghỉ với bộ lọc (Z¡,t > 0) nếu với mọi £ thì ¥; do được đối với ø- trường Z,. Một quá trình ngẫu nhiên X = (X,.f > 0) gọi là một Martingale đối với bộ lọc (7) nếu: i) X tương thích với bộ lọc (7), tức là X, € #; với mọi t.

ii) X/ khả tích với mọi ứ, tức là E|X;| < se với moi t > 0. ii) E(X, | Z2) = X‹ với mọi 0 < s < £, tức là ƒJ(X, — X,)dP = 0 với mọi A € F, va 0<s<t 9 Chương2 PHƯƠNG TRÌNH VI PHÂN NGẪU NHIÊN VÀ PHƯƠNG PHAP RUNGE - KUTTA 2. Lý thuyết phương trình vi phân ngẫu nhiên Định nghĩa 2. Một quá trình Wiener d - chiéu W = (W?,.,W4) n tục, R“ có giá trị thích nghỉ quá trình IV = {W,.Z; :£ € ï} được xác định trên một số không gian xá (9.7, P) với các tính chất sau: () Wi 0h‹e‹e.

với (ii) Cho 0 <s <t, s6 gia W; — W, là độc lập của Z,„ (iii) We — W, 1a phan ph6i N(0,(t — s)fa) « Định nghĩa 2. Kí hiệu £ là lớp tất cả các quá trình ngẫu nhiên X, : 7 x Q =š thỏa mãn B x Z - đo được, tương thích với 7; trong đó B là ø - đại số Borel trong 7 và ¢ B(x) <2 với moi t > 0. 0 Kí hiệu P là lớp của tất cả các quá trình ngẫu nhiên X, : J x 2 => R thỏa mãn B x F - đo được và “ P ( ƒ X?ds< ~) 1 với mọi £>0. ñ Bây giờ ta xem xét một phân hoạch của khoảng cách tích hợp {to, ] da cho bai: <M cnc Néu ta dinh nghia 10 trong đó @ € (0,1] cố định cho trước, khi đó chuỗi các biến ngẫu nhiên.1) hội tụ khi n + oo néu X € P va héi tu theo bình phương trung bình nếu X € £.

Giới hạn không phụ thuộc vào sự lựa chọn của phân hoạch. Tuy nhiên, khác với tích phân Rieman - Stieltjes, ở đây sự lựa chọn Ø tạo ra sự khác biệt. Nếu Ø = 0, khi đó z/” là một điểm giới hạn trái 1Í” ta có phép tính Hô.