MathLabs

Xác suất và Thống kê

Chuyển động Brown, giải tích ngẫu nhiên

Một quá trình ngẫu nhiên liên tục mô hình hóa chuyển động bất định, có phép giải tích riêng dùng trong tài chính và vật lý.

Trực giácTại sao một hạt phấn hoa nhảy múa trong nước lại cần một loại giải tích riêng?

Năm 1827, nhà thực vật học Robert Brown quan sát các hạt phấn hoa lơ lửng trong nước dưới kính hiển vi và thấy chúng rung động không ngừng, đổi hướng ở mọi khoảnh khắc mà không có quy luật nào ông nhận ra được. Lời giải thích, chỉ được hiểu nhiều thập kỷ sau, là mỗi hạt bị hàng tỷ phân tử nước vô hình va đập từ mọi phía mỗi giây, và cú hích ngẫu nhiên ròng ở mỗi khoảnh khắc đẩy hạt đi thêm một chút theo hướng hoàn toàn không đoán trước được. Chuyển động Brown là sự lý tưởng hóa toán học của sự rung động không ngừng, không định hướng này: một đường đi ngẫu nhiên liên tục — nó không bao giờ nhảy cóc — nhưng lại gồ ghề tới mức không có vận tốc xác định rõ ràng ở bất kỳ khoảnh khắc nào, không bao giờ. Chính sự thật đó, rằng đường đi liên tục nhưng không trơn ở đâu cả, buộc các nhà toán học phải xây dựng một loại giải tích hoàn toàn mới, vì giải tích thông thường giả định rằng bạn có thể phóng to một đường cong và thấy nó trông như một đường thẳng.

Đường cong mật độ xác suất Gauss hình chuông của vị trí Brown tại một thời điểm cố định.
Tại mọi thời điểm t>0t > 0, vị trí WtW_t của chuyển động Brown chuẩn có mật độ hình chuông Gauss N(0,t)\mathcal{N}(0, t), với xác suất ≈68%\approx 68\% nằm trong [−t,+t][-\sqrt{t}, +\sqrt{t}].

Đường phân phối phía trên mô tả vị trí có thể có của hạt tại một thời điểm cố định. Chuyển động Brown hai chiều gồm hai bản sao độc lập của quá trình này; thêm thời gian làm tọa độ thứ ba, ta được một quỹ đạo không gian–thời gian gồ ghề.

Quỹ đạo mẫu rời rạc 3D tương tác của hai tọa độ Brown độc lập theo thời gian, có nút điều chỉnh số bước.
Một quỹ đạo mẫu xác định gồm hai tọa độ Brown độc lập theo thời gian. Tọa độ x, y biểu diễn vị trí, z là thời gian; đổi số bước để so sánh các mẫu rời rạc có độ phân giải khác nhau.

Đại họcQuá trình Wiener: một định nghĩa chặt chẽ

Định nghĩa: Quá trình Wiener chuẩn (chuyển động Brown)

Một quá trình Wiener chuẩn (hay chuyển động Brown) là một quá trình ngẫu nhiên thời gian liên tục (Wt)t≥0(W_t)_{t\ge0} sao cho: (1) W0=0W_0=0 hầu chắc chắn; (2) nó có gia số độc lập — với mọi 0≤t1<t2<t3<⋯0\le t_1<t_2<t_3<\cdots, các gia số Wt2−Wt1, Wt3−Wt2,…W_{t_2}-W_{t_1},\ W_{t_3}-W_{t_2},\dots độc lập với nhau; (3) mỗi gia số phân phối chuẩn, Wt−Ws∼N(0,t−s)W_t-W_s\sim N(0,t-s) với s<ts<t; và (4) hầu như mọi quỹ đạo mẫu t↦Wtt\mapsto W_t liên tục. Chỉ riêng tính chất (3) đã buộc quá trình trông giống hệt nhau về phân phối ở mọi thang thời gian — tính tự đồng dạng này là lý do chuyển động Brown xuất hiện như giới hạn phổ quát của rất nhiều bước đi ngẫu nhiên rời rạc khác nhau.

Wt−Ws∼N(0,t−s)W_t-W_s\sim N(0,t-s)

Một đặc điểm định nghĩa, gần như nghịch lý của các quỹ đạo Brown là biến phân bậc hai của chúng: cộng các bình phương gia số của WW trên một phân hoạch mịn của [0,t][0,t] hội tụ (theo xác suất) về chính tt, ∑i(Wti+1−Wti)2→t\sum_{i}(W_{t_{i+1}}-W_{t_i})^2 \to t, bất kể phân hoạch mịn tới đâu — các dao động bình phương không bao giờ triệt tiêu, chúng tích lũy đúng thành tt. So sánh điều này với bất kỳ hàm khả vi thông thường nào ff, có biến phân bậc hai trên bất kỳ khoảng nào luôn là 00, vì gia số của nó co lại như Δt\Delta t chứ không phải Δt\sqrt{\Delta t}. Sự khác biệt duy nhất về tỉ lệ này — gia số Brown có bậc Δt\sqrt{\Delta t}, không phải Δt\Delta t — là nguyên nhân gốc rễ của mọi điều bất thường trong giải tích ngẫu nhiên.

∑i(Wti+1−Wti)2→t\sum_{i}(W_{t_{i+1}}-W_{t_i})^2 \to t
Giải tích thông thường so với giải tích ngẫu nhiên (Itô)
Khía cạnhHàm trơn f(t)f(t)Quỹ đạo Brown WtW_t
Tính khả viKhả vi, tồn tại tiếp tuyếnKhông khả vi ở đâu, không có tiếp tuyến ở đâu
Kích thước gia số trên Δt\Delta tBậc Δt\Delta tBậc Δt\sqrt{\Delta t}
Biến phân bậc hai trên [0,t][0,t]00tt (không bao giờ bằng 0)
Quy tắc dây chuyền cho F(f(t))F(f(t))dF=F′(f)f′(t) dtdF=F'(f)f'(t)\,dtCần thêm một số hạng đạo hàm bậc hai (bổ đề Itô)

Với một quá trình Wiener chuẩn WW và bất kỳ t>0t>0, ∑i(Wti+1−Wti)2→t\sum_{i}(W_{t_{i+1}}-W_{t_i})^2 \to t trong đó tổng lấy trên một phân hoạch của [0,t][0,t] có độ mịn tiến về 00, giới hạn đúng theo xác suất (thực ra hầu chắc chắn dọc theo các phân hoạch nhị phân). Hệ quả, hầu như mọi quỹ đạo mẫu s↦Wss\mapsto W_s không khả vi tại điểm nào s∈[0,t]s\in[0,t].

Vì sao đúng?

Đây là định lý khiến giải tích ngẫu nhiên trở nên cần thiết ngay từ đầu: nó loại trừ việc coi dWtdW_t như một vi phân thông thường theo cách dtdt được xử lý trong giải tích cổ điển, và nó biện minh cho quy tắc trực giác (dWt)2=dt(dW_t)^2=dt dẫn dắt bổ đề Itô dưới đây.

Chứng minh

Biến phân bậc hai. Phân hoạch [0,t][0,t] thành nn mảnh bằng nhau độ dài t/nt/n và đặt Qn=∑i=1n(Wit/n−W(i−1)t/n)2Q_n=\sum_{i=1}^n(W_{it/n}-W_{(i-1)t/n})^2. Mỗi số hạng (Wit/n−W(i−1)t/n)2(W_{it/n}-W_{(i-1)t/n})^2 có kỳ vọng t/nt/n (vì gia số là N(0,t/n)N(0,t/n)) và, dùng mômen bậc bốn của biến chuẩn, phương sai 2(t/n)22(t/n)^2. Cộng nn số hạng độc lập như vậy, E[Qn]=tE[Q_n]=t chính xác và Var(Qn)=2t2/n→0\mathrm{Var}(Q_n)=2t^2/n\to0 khi n→∞n\to\infty, nên Qn→tQ_n\to t theo xác suất (và trung bình bình phương) — điều này chứng minh đẳng thức biến phân bậc hai.

Không khả vi ở đâu. Giả sử, để phản chứng, rằng WW khả vi tại một điểm ss nào đó với đạo hàm W′(s)=LW'(s)=L hữu hạn. Khi đó với hh nhỏ, Ws+h−Ws≈LhW_{s+h}-W_s\approx Lh, nên bình phương gia số trên một khoảng nhỏ độ dài hh sẽ có bậc h2h^2 — nhưng tính toán biến phân bậc hai ở trên cho thấy bình phương gia số điển hình trên một khoảng độ dài hh có bậc hh (lớn hơn nhiều h2h^2 với hh nhỏ), và cộng các số hạng bậc-hh này trên t/ht/h khoảng chính là điều tạo ra một tổng hội tụ về số hữu hạn, dương tt thay vì về 00.

Sự không tương thích này — khả vi sẽ buộc biến phân bậc hai bằng 00, nhưng nó được chứng minh là t>0t>0 — không thể hòa giải, nên không thể tồn tại điểm khả vi nào; một lập luận đầy đủ theo lý thuyết độ đo (của Paley, Wiener và Zygmund năm 1933) làm điều này chặt chẽ bằng cách chặn trực tiếp xác suất mà bất kỳ thương sai phân nào ở gần bất kỳ điểm nào vẫn bị chặn, cho thấy xác suất này đúng bằng không đồng thời cho mọi điểm trên quỹ đạo.

Định lý: Bổ đề Itô

Cho WtW_t là một quá trình Wiener chuẩn và F(t,x)F(t,x) khả vi liên tục hai lần theo xx và một lần theo tt. Khi đó quá trình F(t,Wt)F(t,W_t) thỏa dF=(∂tF+12∂x2F)dt+∂xF dWtdF=\Big(\partial_tF+\tfrac12\partial_x^2F\Big)dt+\partial_xF\,dW_t — một quy tắc dây chuyền ngẫu nhiên với một số hạng đạo hàm bậc hai thêm vào ("hiệu chỉnh Itô") so với quy tắc dây chuyền thông thường.

Vì sao đúng?

Đây là công cụ được dùng nhiều nhất trong giải tích ngẫu nhiên: nó cho bạn biết chính xác cách lấy vi phân một hàm của một đường đi ngẫu nhiên, điểm khởi đầu để suy ra động lực của bất kỳ đại lượng nào (giá quyền chọn, một đại lượng vật lý quan sát được) phụ thuộc vào chuyển động Brown.

Chứng minh

Khai triển Taylor. Với FF trơn, khai triển Taylor bậc hai thông thường theo cả hai biến trên một bước nhỏ Δt\Delta t với ΔW=Wt+Δt−Wt\Delta W=W_{t+\Delta t}-W_t có dạng ΔF≈∂tF Δt+∂xF ΔW+12∂x2F (ΔW)2+12∂t2F (Δt)2+∂t∂xF Δt ΔW+⋯\Delta F\approx \partial_tF\,\Delta t+\partial_xF\,\Delta W+\tfrac12\partial_x^2F\,(\Delta W)^2+\tfrac12\partial_t^2F\,(\Delta t)^2+\partial_t\partial_xF\,\Delta t\,\Delta W+\cdots — đây thuần túy là giải tích, đúng cho bất kỳ đường trơn nào.

Sắp thứ tự các số hạng theo kích thước. Vì ΔW\Delta W có bậc Δt\sqrt{\Delta t} (không phải Δt\Delta t, như đã chỉ ra trong định lý trước), các số hạng có tỉ lệ: ∂tF Δt\partial_tF\,\Delta t bậc Δt\Delta t; ∂xF ΔW\partial_xF\,\Delta W bậc Δt\sqrt{\Delta t} (số hạng ngẫu nhiên chủ đạo, bậc dẫn đầu); (Δt)2(\Delta t)^2 và Δt ΔW\Delta t\,\Delta W bậc (Δt)2(\Delta t)^2 và (Δt)3/2(\Delta t)^{3/2} tương ứng, không đáng kể so với Δt\Delta t; nhưng (ΔW)2(\Delta W)^2 bậc Δt\Delta t — cùng bậc với ∂tF Δt\partial_tF\,\Delta t, hoàn toàn không đáng bỏ qua, khác với giải tích thông thường nơi (Δx)2(\Delta x)^2 luôn không đáng kể so với Δx\Delta x.

Thay (ΔW)2(\Delta W)^2 bằng kỳ vọng của nó. Theo định lý biến phân bậc hai, cộng trên nhiều bước nhỏ (ΔW)2(\Delta W)^2 hành xử như Δt\Delta t (kỳ vọng của nó, với các dao động quanh kỳ vọng đó triệt tiêu khi các bước co lại và được cộng lại), đó là sự biện minh không chính thức cho quy tắc thay thế trực giác (dWt)2=dt(dW_t)^2=dt trong giới hạn các bước vô cùng nhỏ.

Lấy giới hạn. Bỏ các số hạng bậc cao không đáng kể và thay (ΔW)2→dt(\Delta W)^2\to dt trong giới hạn Δt→0\Delta t\to0 biến khai triển Taylor thành dạng vi phân dF=(∂tF+12∂x2F)dt+∂xF dWtdF=\Big(\partial_tF+\tfrac12\partial_x^2F\Big)dt+\partial_xF\,dW_t chính xác như đã khẳng định — số hạng 12∂x2F dt\tfrac12\partial_x^2F\,dt chính xác là đóng góp mà giải tích thông thường bỏ qua nhưng giải tích ngẫu nhiên phải giữ lại.

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

Ngoài mô tả sự khuếch tán vật lý (phấn hoa trong nước, nhiệt lan qua vật rắn, các phân tử khí hòa trộn), chuyển động Brown và giải tích Itô trở thành nền tảng toán học của tài chính định lượng hiện đại. Giá cổ phiếu StS_t thường được mô hình hóa như chuyển động Brown hình học, dSt=μSt dt+σSt dWtdS_t=\mu S_t\,dt+\sigma S_t\,dW_t, trong đó μ\mu là tốc độ tăng trưởng trung bình và σ\sigma là độ biến động. Áp dụng bổ đề Itô cho giá trị V(t,St)V(t,S_t) của một sản phẩm phái sinh tài chính (như quyền chọn) dẫn tới phương trình đạo hàm riêng nổi tiếng Black–Scholes, ∂tV+12σ2S2∂S2V+rS∂SV−rV=0\partial_tV+\tfrac12\sigma^2S^2\partial_S^2V+rS\partial_SV-rV=0, mà nghiệm của nó cho giá công bằng của các quyền chọn giao dịch trên mọi sàn lớn. Cùng toán học đó — quá trình Wiener và bổ đề Itô — cũng làm nền tảng cho các mô hình điện thế màng neuron trong khoa học thần kinh, mô hình lãi suất trong kinh tế học, và sự phân tán rối trong động lực học chất lỏng.

Ví dụ: Kỳ vọng và phương sai của một độ dịch chuyển Brown

Vị trí của một hạt tuân theo một quá trình Wiener chuẩn WtW_t trong một chiều, xuất phát từ W0=0W_0=0. Tìm E[W5]E[W_5] và Var(W5−W2)\mathrm{Var}(W_5-W_2).

Lời giải

Với đại lượng đầu, dùng Wt−Ws∼N(0,t−s)W_t-W_s\sim N(0,t-s) với s=0s=0: W5−W0=W5W_5-W_0=W_5 là N(0,5)N(0,5), nên E[W5]=0E[W_5]=0 trực tiếp — chuyển động Brown không có xu hướng, nên vị trí kỳ vọng của nó không bao giờ rời xa điểm xuất phát.

Với đại lượng thứ hai, áp dụng cùng tính chất định nghĩa với s=2s=2, t=5t=5: W5−W2∼N(0,5−2)W_5-W_2\sim N(0,5-2), tức N(0,3)N(0,3).

Phương sai của một biến ngẫu nhiên N(0,σ2)N(0,\sigma^2) là σ2\sigma^2 theo định nghĩa, nên đọc trực tiếp tham số, Var(W5−W2)=3\mathrm{Var}(W_5-W_2)=3.

Chú ý phương sai chỉ phụ thuộc vào thời gian trôi qua t−s=3t-s=3, không phụ thuộc vào bản thân thời điểm xuất phát s=2s=2 — điều này phản ánh tính thuần nhất theo thời gian được xây dựng sẵn trong định nghĩa của quá trình Wiener.

Ví dụ: Áp dụng bổ đề Itô cho F(t,x)=x2F(t,x)=x^2

Dùng bổ đề Itô để tìm d(Wt2)d(W_t^2), vi phân ngẫu nhiên của quá trình Wiener bình phương, và dùng kết quả để xác nhận E[Wt2]=tE[W_t^2]=t.

Lời giải

Đặt F(t,x)=x2F(t,x)=x^2, nên ∂tF=0\partial_tF=0, ∂xF=2x\partial_xF=2x, ∂x2F=2\partial_x^2F=2.

Thay vào bổ đề Itô dF=(∂tF+12∂x2F)dt+∂xF dWtdF=\Big(\partial_tF+\tfrac12\partial_x^2F\Big)dt+\partial_xF\,dW_t với các đạo hàm này (tính tại x=Wtx=W_t): d(Wt2)=(0+12⋅2)dt+2Wt dWtd(W_t^2)=\big(0+\tfrac12\cdot2\big)dt+2W_t\,dW_t, đơn giản hóa thành d(Wt2)=dt+2Wt dWtd(W_t^2)=dt+2W_t\,dW_t.

Tích phân cả hai vế từ 00 tới tt (dùng W02=0W_0^2=0): Wt2=t+2∫0tWs dWsW_t^2=t+2\int_0^tW_s\,dW_s.

Lấy kỳ vọng, và dùng sự thật then chốt rằng một tích phân Itô ∫0tWs dWs\int_0^tW_s\,dW_s luôn có kỳ vọng bằng không (nó được xây từ các gia số độc lập với quá khứ, nên không có xu hướng hệ thống nào để tích lũy): E[Wt2]=t+2E[∫0tWs dWs]=t+0=tE[W_t^2]=t+2E\big[\int_0^tW_s\,dW_s\big]=t+0=t, xác nhận trực tiếp điều ta đã biết từ Wt∼N(0,t)W_t\sim N(0,t) (có phương sai tt), nhưng lần này được suy ra thuần túy từ bộ máy giải tích ngẫu nhiên thay vì từ định nghĩa.

Với một quá trình Wiener chuẩn, gia số Wt−WsW_t-W_s (với s<ts<t) phân phối theo:

Biến phân bậc hai của chuyển động Brown trên [0,t][0,t] là:

So với quy tắc dây chuyền thông thường, bổ đề Itô dF=(∂tF+12∂x2F)dt+∂xF dWtdF=\Big(\partial_tF+\tfrac12\partial_x^2F\Big)dt+\partial_xF\,dW_t có thêm một số hạng vì:

Phương trình đạo hàm riêng Black–Scholes cho định giá quyền chọn được suy ra bằng cách áp dụng:

Tài liệu tham khảo

  1. Ioannis Karatzas, Steven E. Shreve (1991). Brownian Motion and Stochastic Calculus
  2. Kiyosi Itô (1944). Stochastic Integral
  3. Fischer Black, Myron Scholes (1973). The Pricing of Options and Corporate Liabilities