Tổng quan luận án

Luận án tiến sĩ "Models of Selected Problems in Mathematical Finance and Numerical Methods for Stochastic Differential Equations" (Các mô hình cho các bài toán chọn lọc trong tài chính toán học và các phương pháp số cho phương trình vi phân ngẫu nhiên) được thực hiện bởi nghiên cứu sinh Timothy L. Seaman dưới sự hướng dẫn của Phó Giáo sư Harbir Lamba tại Khoa Khoa học Toán học, thuộc chương trình Tiến sĩ Khoa học Tính toán và Tin học (Computational Sciences and Informatics) của Đại học George Mason (Fairfax, Virginia), hoàn thành vào học kỳ mùa thu năm 2006. Hội đồng chấm luận án gồm Tiến sĩ Harbir Lamba (Trưởng ban), Tiến sĩ Timothy Sauer, Tiến sĩ James Gentle, Phó Trưởng khoa Becker và Trưởng khoa Menas Kafatos.

Về tính cấp thiết và khoảng trống nghiên cứu, công trình giải quyết hai vấn đề độc lập trong khoa học tính toán và mô hình hóa toán học:

  1. Trong tài chính toán học: Lịch sử kinh tế ghi nhận nhiều chu kỳ bất thường, bong bóng và sụp đổ thị trường như hội chứng hoa tulip thế kỷ 17, bong bóng South Sea thế kỷ 18, cùng các cuộc khủng hoảng thị trường chứng khoán năm 1929, 1987 và 2000. Giả thuyết thị trường hiệu quả (Efficient Market Hypothesis - EMH) dựa trên giả định phân phối chuẩn của Louis Bachelier (1900) đã bộc lộ giới hạn nghiêm trọng khi dữ liệu thực nghiệm cho thấy sự xuất hiện nhất quán của các đặc tính thống kê phi Gauss: phân phối lợi suất có đuôi dày (fat tails), độ nhọn cao (leptokurtosis) và hiện tượng cụm biến động (volatility clustering). Các mô hình vi mô dựa trên tác tử trước đây thường quá phức tạp, chứa nhiều tham số khó ước tính hoặc phụ thuộc vào các giả định nhân tạo như rào chắn phản xạ để tránh rơi vào trạng thái hấp thụ. Luận án đặt ra yêu cầu xây dựng một lớp mô hình dựa trên tác tử đơn giản, có số lượng tham số cơ bản rất nhỏ và có thể quy về mô hình EMH ở các giới hạn tham số nhất định để làm rõ mối quan hệ nhân quả của các hiện tượng phi chuẩn.
  2. Trong giải số phương trình vi phân ngẫu nhiên (SDE): SDE chứa số hạng cưỡng bức ngẫu nhiên Brownian đóng vai trò nền tảng trong nhiều ngành khoa học. Tuy nhiên, các phương pháp số bậc cao cho SDE trở nên cực kỳ phức tạp và tốn kém tài nguyên tính toán do sự phức tạp của các khai triển Itô-Taylor và Stratonovich-Taylor. Việc áp dụng trực tiếp các kỹ thuật thích nghi bước thời gian biến thiên từ phương trình vi phân thường (ODE) sang SDE không đem lại hiệu quả tương đương. Thuật toán kiểm soát sai số kép của Harbir Lamba (2003) trên phương pháp Milstein bậc 1 gặp phải các vấn đề kỹ thuật: sự sụp đổ của ước lượng sai số khuếch tán đối với SDE có nhiễu cộng (additive noise), hạn chế của phương pháp hiện đối với các bài toán cứng (stiff SDEs), và cơ chế lựa chọn bước thời gian còn mang tính đặc thù.

Mục tiêu và câu hỏi nghiên cứu của luận án được xác định cụ thể:

  • Mục tiêu 1: Xây dựng và mô phỏng một lớp mô hình thị trường tài chính dựa trên tác tử có hành vi ngưỡng, kết hợp giữa tính duy lý và tính duy lý giới hạn, nhằm tái tạo các đặc trưng thực nghiệm (stylized facts) của chuỗi thời gian tài chính thực tế.
  • Mục tiêu 2: Phân tích tính ổn định trung bình bình phương (mean-square stability) của thuật toán thích nghi bước thời gian sử dụng chiến lược kiểm soát sai số kép trên các bài toán kiểm tra tuyến tính có nhiễu nhân; mở rộng thuật toán sang phương pháp Milstein cân bằng (Balanced Milstein Method) cho SDE cứng; và khảo sát các bộ kiểm soát sai số ngẫu nhiên bổ sung trên SDE có nhiễu cộng.
  • Giả thuyết và câu hỏi nghiên cứu:
    1. Liệu hai sai lệch phi Gauss quan trọng nhất trong lợi suất giá—gồm hiện tượng đuôi dày và hiện tượng cụm biến động—có bắt nguồn từ hai cơ chế nội sinh riêng biệt hay không?
    2. Thuật toán bước thời gian biến thiên kiểm soát sai số kép có bảo toàn tính ổn định trung bình bình phương đối với lớp bài toán kiểm tra có nhiễu nhân hay không?
    3. Việc tích hợp phương pháp Milstein cân bằng và các bộ kiểm soát sai số ngẫu nhiên bổ sung có giúp cải thiện độ chính xác và hiệu năng tính toán so với các phương pháp bước cố định tương ứng hay không?

Về đối tượng và phạm vi nghiên cứu: Đối tượng nghiên cứu là các chuỗi thời gian giá và lợi suất tài sản tài chính trong mô phỏng vi mô, cùng sai số cục bộ và tính ổn định số của nghiệm SDE vô hướng. Phạm vi không gian của mô hình tài chính là một thị trường đóng với $N$ tác tử có quy mô thời gian dài hạn (khảo sát tại $M = 100$ và $M = 89$); phạm vi thời gian mô phỏng tương đương chu kỳ 40 năm giao dịch. Phạm vi của phần SDE giới hạn trong các phương trình vi phân ngẫu nhiên vô hướng chịu sự kích thích của một quá trình ngẫu nhiên vô hướng đơn lẻ trên khoảng thời gian tích phân $[0, 1]$.

Tổng quan tài liệu và vị trí của luận án

Luận án điểm lại hệ thống tài liệu khoa học theo hai nhánh nội dung:

Nhánh 1: Các mô hình thị trường tài chính và giả thuyết thị trường hiệu quả:

  • Khởi đầu từ Louis Bachelier (1900, "Théorie de la Spéculation"), người đầu tiên mô hình hóa biến động giá cổ phiếu bằng chuyển động Brown và phân phối chuẩn trước Albert Einstein (1905). Bachelier cho rằng tính thanh khoản và việc triệt tiêu cơ hội kinh doanh chênh lệch giá (arbitrage) làm cho giá cổ phiếu bước đi ngẫu nhiên. Nhược điểm giá có thể âm được khắc phục sau đó bằng chuyển động Brown hình học.
  • M. F. M. Osborne (1959) tiếp cận thị trường theo góc độ cơ học thống kê, xem xét giá cổ phiếu là kết quả của một tập hợp quyết định ở trạng thái cân bằng thống kê, trong đó logarit của thay đổi giá tương đối có phân phối dừng tương tự hạt chuyển động Brown.
  • Benoit Mandelbrot (1960, 1963) áp dụng luật Pareto-Lévy vào phân phối thu nhập và biến động giá đầu cơ (dữ liệu giá bông), chứng minh rằng phân phối giá thực tế có đỉnh nhọn và đuôi dày hơn phân phối chuẩn, dẫn đến độ nhọn (kurtosis) vượt xa giá trị 3 của phân phối Gauss (hiện tượng leptokurtosis). Mandelbrot đề xuất thay thế phân phối chuẩn bằng phân phối Paretian ổn định.
  • Eugene Fama (1963, 1965) củng cố các phát hiện của Mandelbrot qua khảo sát sai phân bậc một logarit giá cổ phiếu của các công ty lớn, xác nhận số mũ đặc trưng xấp xỉ 2 nhưng nhỏ hơn 2. Fama (1970) hệ thống hóa ba dạng của giả thuyết thị trường hiệu quả: dạng yếu, bán mạnh và mạnh; sau đó Fama (1991) khẳng định tính hiệu quả của thị trường không thể kiểm định độc lập mà phải kiểm định đồng thời với một mô hình cân bằng định giá tài sản.
  • Fischer Black và Myron Scholes (1972) xây dựng công thức định giá quyền chọn dựa trên giả định lợi suất hình học tuân theo phân phối chuẩn. McCauley chỉ ra rằng giả định chuẩn này sai lệch so với dữ liệu thực nghiệm vốn có đuôi dày.
  • Adrian Pagan (1996) tổng kết các dữ kiện kinh tế lượng tài chính: bình phương lợi suất có tương quan bền vững với sự phân rã chậm của hệ số tự tương quan; tồn tại hiệu ứng đòn bẩy (leverage effect); tồn tại mô-men bậc hai nhưng thiếu bằng chứng về sự tồn tại của mô-men bậc bốn.
  • Rama Cont (2001) tổng hợp danh sách 10 đặc trưng thực nghiệm được cách điệu hóa (stylized empirical facts): không có tự tương quan tuyến tính trong lợi suất; đuôi nặng dạng luật lũy thừa hoặc Pareto; cụm biến động; tính bất đối xứng lãi/lỗ; tính chuẩn khi tổng hợp thời gian; tính gián đoạn của lợi suất; đuôi nặng có điều kiện; hàm tự tương quan của trị tuyệt đối lợi suất phân rã chậm; hiệu ứng đòn bẩy; và tương quan giữa khối lượng giao dịch hoặc độ biến động với lợi suất.
  • Sergei Maslov (1999) đánh giá 4 mô hình vi mô dựa trên tiêu chí đơn giản:
    1. Mô hình Bak, Paczuski, Shubik (BPS, 1997): Kết hợp nhà giao dịch duy lý (tối ưu hóa hàm thỏa dụng cổ tức) và nhà giao dịch nhiễu. Khi có 2% nhà giao dịch duy lý, mô hình tạo ra bong bóng giá và đuôi dày; khi nâng lên 20%, các biến động lớn biến mất. Nhược điểm là giới hạn giá nhân tạo và chỉ tái hiện được 2 đặc trưng thực nghiệm.
    2. Mô hình Cont và Bouchaud (CB, 2000): Liên kết đuôi nặng với hành vi bầy đàn trên đồ thị ngẫu nhiên dựa trên lý thuyết thấm Flory-Stockmayer. Phương trình biến động giá tỉ lệ nghịch với độ sâu thị trường $\lambda$. Stauffer và Penna (1996) mở rộng mô hình CB lên mạng tinh thể $d$ chiều ($d = 2$ đến 7). Nhược điểm là chỉ tái hiện được đuôi dày và sự chuyển tiếp sang phân phối chuẩn ở thang thời gian lớn.
    3. Mô hình Caldarelli, Marsili, Zhang (CMZ, 1997): Dựa trên bài toán El Farol Bar của Brian Arthur (1994), mô phỏng đầu cơ ngoại hối trong hệ kín, tái tạo đuôi lũy thừa nhưng số mũ nhỏ hơn thực tế và không có chuyển tiếp sang phân phối chuẩn.
    4. Mô hình Trò chơi thiểu số (Minority Game) của Damien Challet và Yi-Cheng Zhang (CZ, 1997): Người chơi dùng điểm ảo để lựa chọn chiến lược thuộc nhóm thiểu số; việc đặt ngưỡng giúp giảm tần suất chuyển đổi chiến lược. Các mở rộng tiếp theo gồm Challet, Marsili, Zhang (2000) đưa vào nhà sản xuất và nhà đầu cơ; Marsili, Challet, Zecchina (2000) dùng vật lý thống kê giải mô hình cho vô hạn tác tử; Marsili và Challet (2001) chỉ ra việc chấp nhận giá và học nhanh gây ra biến động dư thừa.
  • Các mô hình khác: Jefferies, Hart, Hui, Johnson (2001) với cơ chế bù trừ thị trường; Didier Sornette (2003) với mô hình điều khiển bởi rủi ro và giá; Maslov (2003) với mô hình thị trường lệnh giới hạn có số mũ Hurst $H \approx 0.5$; chuỗi mô hình của Thomas Lux gồm Lux (1995) dựa trên thuyết hiệp đồng (synergetics), Lux (1998) cho phép chuyển đổi giữa chiến lược cơ bản và kỹ thuật, Lux và Marchesi (LM, 1999) với 3 nhóm tác tử và hệ phương trình vi phân có 3 trạng thái dừng (nhược điểm là cần tới 10 tham số hiệu chuẩn), và Lux và Alfarano (2003) đơn giản hóa nhưng phải dùng rào chắn phản xạ nhân tạo.

Nhánh 2: Nghiệm số của phương pháp vi phân ngẫu nhiên:

  • Kloeden và Platen cung cấp nền tảng toán học cho các lược đồ số SDE.
  • Các kỹ thuật ODE kinh điển được điểm lại gồm: thiết bị Milne dùng hai phương pháp đa bước cùng bậc; ngoại suy Richardson dùng bước $2h$ và hai bước $h$; cặp Runge-Kutta nhúng bậc $p$ và $p+1$ (Hairer, Norsett, Wanner) với công thức bước tối ưu $h_{opt} = h \cdot (1/err)^{1/(p+1)}$ kết hợp hệ số an toàn $fac_{safe}$ (0.8 - 0.9) và hệ số giới hạn $fac_{max}$ (1.5 - 5).
  • Gaines và Lyons (1997) xây dựng thuật toán giải SDE bước biến thiên sử dụng cây nhị phân lưu trữ quỹ đạo Brown.
  • Harbir Lamba (2003) đề xuất thuật toán thích nghi bước thời gian dựa trên kiểm soát sai số kép áp dụng cho phương pháp Milstein bậc 1, tách riêng sai số trôi (drift) và sai số khuếch tán (diffusion), cho phép từ chối bước tính ngay sau khi sinh số gia Brown.

Khoảng trống luận án lựa chọn giải quyết: Trong mô hình hóa tài chính, luận án xây dựng mô hình ngưỡng tập trung vào các nhà đầu tư dài hạn với số lượng tham số tối thiểu, loại bỏ các giả định nhân tạo (như rào chắn phản xạ) và cho phép thu hồi mô hình EMH tại các giới hạn tham số cụ thể. Trong giải số SDE, luận án giải quyết vấn đề sụp đổ của kiểm soát sai số khuếch tán khi gặp nhiễu cộng, mở rộng phương pháp thích nghi sang phương pháp Milstein cân bằng ẩn để giải SDE cứng, và chứng minh tính ổn định trung bình bình phương trên các bài toán kiểm tra tuyến tính.

Cơ sở lý thuyết và phương pháp nghiên cứu

Lý thuyết, khái niệm và khung phân tích:

  • Khái niệm tính duy lý giới hạn (bounded rationality) của Herbert Simon (1957, 1982): Khi khối lượng thông tin vượt quá khả năng xử lý, tác tử dừng thu thập thông tin và đưa ra lựa chọn dựa trên lượng thông tin hiện có tại thời điểm đó.
  • Khái niệm ngưỡng tâm lý và cấu trúc thị trường: Tác tử sẵn sàng đánh đổi lợi nhuận kỳ vọng rủi ro cao để lấy lợi nhuận nhỏ bảo đảm; các nhà quản lý quỹ chuyên nghiệp chịu áp lực từ các tiêu chí đánh giá hiệu quả và cơ chế thưởng nội bộ.
  • Khung phân tích 10 đặc trưng thực nghiệm của Rama Cont (2001).
  • Phương pháp Milstein cấp 1 hiện (Explicit Milstein Method) và phương pháp Milstein cân bằng (Balanced Milstein Method - Implicit) cho nghiệm số SDE vô hướng.
  • Lý thuyết ổn định trung bình bình phương (mean-square stability) cho các phương trình vi phân ngẫu nhiên tuyến tính có nhiễu nhân: $$dX(t) = \lambda X(t) dt + \mu X(t) dW(t)$$
  • Chiến lược kiểm soát sai số kép (dual error control): Sai số trôi $E_d$ tương ứng với sai số xác định và sai số khuếch tán $E_s$ tương ứng với sai số tích phân ngẫu nhiên. Khi độ khuếch tán tiến về 0, thuật toán tiệm cận về một bộ giải ODE thích nghi thông thường.

Phương pháp thực tế:

  • Phương pháp mô phỏng số dựa trên tác tử: Mô phỏng hành vi của các nhà đầu tư dài hạn có vị thế mua (long) hoặc bán (short). Tác tử liên tục điều chỉnh trạng thái nội tại; khi trạng thái này vượt ngưỡng chịu đựng, tác tử thực hiện giao dịch chuyển đổi vị thế. Tỷ lệ giữa vị thế mua và bán xác định tâm lý thị trường tổng thể. Nhiễu Brown ngoại sinh đại diện cho dòng tin tức mới được các nhà giao dịch kinh doanh chênh lệch giá và nhà đầu tư cơ bản xử lý tức thời ở thang thời gian ngắn.
  • Phương pháp phân tích chuỗi thời gian thống kê: Tính toán số mũ Hurst ($H$), hàm tự tương quan của lợi suất và bình phương lợi suất, phân phối tích lũy của log-lợi suất, ước lượng số mũ luật lũy thừa, phân tích tính bất đối xứng giữa thị trường giá lên (bull market) và giá xuống (bear market).
  • Phương pháp thử nghiệm số cho SDE: Thực hiện tích phân số với các mức dung sai $\tau$ khác nhau. Đo lường số bước tính trung bình, số bước được chấp nhận, số bước bị từ chối, sai số tuyệt đối cực đại trung bình và sai số trung bình tại thời điểm $t = 1$. Thử nghiệm trên Test Problem 1 (ổn định trung bình bình phương với $\lambda = -2$ và $\mu^2 = 4$) và Test Problem 2 (chia thành bốn trường hợp: không cứng - hiện, không cứng - cân bằng, cứng - hiện, cứng - cân bằng).

Nguồn tư liệu và dữ liệu: Luận án sử dụng dữ liệu mô phỏng tổng hợp (synthetic simulated data) từ mô hình tác tử trong thời gian tương đương 40 năm giao dịch; dữ liệu số từ các bài toán kiểm tra chuẩn của SDE với quỹ đạo chuyển động Brown được sinh ngẫu nhiên trên máy tính.

Nội dung chính theo từng chương

Chương 1: Giới thiệu (Chapter 1: Introduction)

Chương 1 thiết lập nền tảng nghiên cứu qua hai phân nhánh: mô hình hóa thị trường tài chính và giải số SDE. Tác giả phác thảo lịch sử tư tưởng kinh tế từ Adam Smith, Thomas Malthus, David Ricardo, Karl Marx đến Léon Walras với hệ thống cân bằng tổng thể Walrasian (Walrasian General Equilibrium). Luận án phân tích việc đưa toán học vào kinh tế tài chính từ luận án năm 1900 của Louis Bachelier, đặt nền móng cho chuyển động Brown và giả thuyết thị trường hiệu quả (EMH). Tác giả chỉ ra mâu thuẫn giữa giả định phân phối chuẩn của EMH với các phát hiện thực nghiệm về đuôi dày của Benoit Mandelbrot (1963) và Eugene Fama (1965). Về triết lý mô hình hóa, tác giả lựa chọn tập trung vào các tác tử có tầm nhìn dài hạn (tháng đến năm), mô phỏng trạng thái nội tại và hành vi ngưỡng thay vì gán xác suất chuyển đổi ngẫu nhiên kiểu Markov. Về phần SDE, tác giả trình bày sự phức tạp của các khai triển Itô-Taylor và Stratonovich-Taylor, chi phí tính toán cao của các phương pháp bậc cao, và sự cần thiết của thuật toán thích nghi bước thời gian dựa trên kiểm soát sai số kép trên phương pháp Milstein bậc 1.

Chương 2: Tổng quan tài liệu (Chapter 2: Review of the Literature)

Chương 2 cung cấp bức tranh chi tiết về lịch sử nghiên cứu:

  • Mục 2.1 đánh giá các mô phỏng thị trường dựa trên tác tử: Phân tích chi tiết mô hình Bak, Paczuski, Shubik (1997); mô hình Cont-Bouchaud (2000) và mở rộng của Stauffer-Penna (1996); mô hình Caldarelli, Marsili, Zhang (1997); mô hình Trò chơi thiểu số (Minority Game) của Challet-Zhang (1997) cùng các mở rộng của Marsili, Challet, Zecchina (2000, 2001); mô hình thị trường lệnh giới hạn của Maslov (2003); và đặc biệt là hệ thống mô hình của Thomas Lux (1995, 1998, 1999, 2003). Tác giả nhấn mạnh mô hình Lux-Marchesi (1999) dù tái hiện tốt số mũ Hurst ($H = 0.48$ cho lợi suất thô, $H = 0.85$ cho trị tuyệt đối lợi suất) nhưng đòi hỏi tới 10 tham số hiệu chuẩn, trong khi mô hình Lux-Alfarano (2003) phải dựa vào rào chắn phản xạ nhân tạo.
  • Mục 2.2 tổng quan các kỹ thuật số giải SDE: Điểm lại các kỹ thuật kiểm soát sai số cục bộ trong ODE (thiết bị Milne, ngoại suy Richardson, cặp Runge-Kutta nhúng bậc $p$ và $p+1$ của Hairer et al., công thức tính bước tối ưu $h_{opt}$). Chuyển sang SDE, tác giả phân tích thuật toán biểu diễn quỹ đạo Brown bằng cây nhị phân của Gaines và Lyons (1997) và thuật toán kiểm soát sai số kép của Harbir Lamba (2003).

Chương 3: Phương pháp luận (Chapter 3: Methodology)

Chương 3 trình bày chi tiết cấu trúc toán học của hai mô hình:

  • Mục 3.1 thiết lập mô hình ngưỡng về hành vi nhà đầu tư (A Threshold Model of Investor Behavior):
    • Tiểu mục 3.1.1 trình bày cơ chế xác định giá tài sản: Giá tài sản phản ánh tác động đồng thời của thành phần thông tin ngoại sinh (chuyển động Brown) thông qua các nhà đầu tư cơ bản/kinh doanh chênh lệch giá ở thang thời gian ngắn và thành phần mất cân bằng cung cầu do các tác tử dài hạn tạo ra.
    • Tiểu mục 3.1.2 mô hình hóa vị thế của nhà đầu tư: Trạng thái nội tại của mỗi tác tử chịu tác động của các yếu tố duy lý và tâm lý. Khi trạng thái vượt quá một ngưỡng giới hạn, tác tử chuyển đổi vị thế.
    • Bảng 3.1 so sánh số mũ Hurst tính toán được từ các phương pháp khác nhau so với số mũ Hurst của chuyển động Brown phân số đầu vào nhằm xác minh độ chính xác của các thuật toán đo lường chuỗi thời gian.
  • Mục 3.2 phát triển các phương pháp số cho SDE:
    • Tiểu mục 3.2.1 mô tả phương pháp Milstein hiện cấp 1.
    • Tiểu mục 3.2.2 phát triển phương pháp Milstein cân bằng (Balanced Milstein Method) ẩn để xử lý các bài toán cứng.
    • Tiểu mục 3.2.3 thiết lập cơ chế thích nghi sử dụng hai bộ kiểm soát sai số $E_d$ (thành phần trôi) và $E_s$ (thành phần khuếch tán).
    • Tiểu mục 3.2.4 đề xuất các bộ kiểm soát sai số ngẫu nhiên bổ sung ($E_{d2}$, $E_{s2}$) để xử lý các bài toán SDE có nhiễu cộng khi ước lượng sai số khuếch tán thông thường bị mất tác dụng.

Chương 4: Kết quả và thảo luận (Chapter 4: Results and Discussion)

Chương 4 trình bày các phát hiện thực nghiệm và phân tích số liệu:

  • Mục 4.1 trình bày kết quả mô phỏng thị trường tài chính:
    • 4.1.1 Hiệu chuẩn mô hình.
    • 4.1.2 Mô phỏng khi loại bỏ áp lực bầy đàn và đặt $a = 0$: Kết quả tại Hình 4.1 ($M = 100$) và Hình 4.2 ($M = 89$) cho thấy chuỗi giá không xuất hiện các biến động phi chuẩn cực đoan.
    • 4.1.3 Mô phỏng có bầy đàn nhưng $a = 0$ trong chu kỳ 40 năm (Hình 4.3).
    • 4.1.4 Mô phỏng mô hình đầy đủ (có bầy đàn và $a \neq 0$) trong chu kỳ 40 năm (Hình 4.4): Tái hiện đầy đủ các đặc trưng thực nghiệm. Phân phối tích lũy của log-lợi suất giá (dương, âm và trị tuyệt đối) thể hiện rõ đuôi dày tuân theo luật lũy thừa (Hình 4.5). Hàm tự tương quan của độ biến động phân rã chậm theo dạng luật lũy thừa (Hình 4.6). Phân phối lợi suất logarit chuẩn hóa tách riêng phần tăng và giảm (Hình 4.7). Luận điểm then chốt: kết quả chỉ ra hiện tượng đuôi dày và hiện tượng cụm biến động bắt nguồn từ các cơ chế nội sinh khác nhau trong mô hình.
    • 4.1.5 Đưa tính bất đối xứng giá vào mô hình: Phân phối tích lũy trị tuyệt đối chuẩn hóa của thay đổi giá tương đối trong thị trường giá lên (Hình 4.8) và thị trường giá xuống (Hình 4.9) phản ánh chính xác tính bất đối xứng lãi/lỗ thực tế.
  • Mục 4.2 trình bày các kết quả thực nghiệm số cho SDE:
    • 4.2.1 Phân tích tính ổn định trung bình bình phương: Bảng 4.1 trình bày các giá trị tính toán của $E(X_n^2)$ và bước thời gian $\Delta t$ tương ứng khi cố định $\lambda = -2$; Bảng 4.2 ghi nhận kết quả khi cố định $\mu^2 = 4$ ($\mu = 2$), xác nhận thuật toán bảo toàn tính ổn định trung bình bình phương trên lớp bài toán kiểm tra.
    • 4.2.2 Cải tiến phương pháp tiếp cận kiểm soát sai số kép:
      • Bảng 4.3 và Bảng 4.4 so sánh số bước thực hiện trung bình, nguyên nhân chấm dứt tăng kích thước bước, sai số cực đại trung bình trên mỗi mẻ và sai số trung bình giữa hệ thống 2 bộ kiểm soát ($E_d, E_s$) và 4 bộ kiểm soát ($E_d, E_s, E_{d2}, E_{s2}$) dưới các mức dung sai $\tau$.
      • Bảng 4.5 cung cấp sai số tuyệt đối cực đại trung bình và sai số tại $t = 1$ của nghiệm giải Test Problem 2 sử dụng phương pháp Milstein bước cố định $h_{fix}$.
      • Bảng 4.6 đến Bảng 4.9 đánh giá hiệu năng giải Test Problem 2 của các biến thể thích nghi: Không cứng, Hiện (Bảng 4.6); Không cứng, Cân bằng (Bảng 4.7); Cứng, Hiện (Bảng 4.8); Cứng, Cân bằng (Bảng 4.9). Các bảng ghi nhận chi tiết sai số tuyệt đối cực đại trung bình, sai số tại $t = 1$, số bước thực hiện trung bình, số bước được chấp nhận và số bước bị từ chối.
      • Bảng 4.10 thống kê số bước mà từng bộ kiểm soát sai số thực hiện chấm dứt việc tăng kích thước bước thời gian.

Dưới đây là bảng tổng hợp các cấu trúc dữ liệu thực nghiệm số chính được tác giả xây dựng trong Chương 4:

Ký hiệu bảng trong luận án Nội dung khảo sát thực nghiệm số Các biến số và chỉ số đo lường chính
Bảng 4.1 Ổn định trung bình bình phương với $\lambda = -2$ cố định $E(X_n^2)$, bước thời gian $\Delta t$
Bảng 4.2 Ổn định trung bình bình phương với $\mu^2 = 4$ ($\mu = 2$) cố định $E(X_n^2)$, bước thời gian $\Delta t$
Bảng 4.3 Thuật toán thích nghi dùng 2 bộ kiểm soát ($E_d, E_s$) Dung sai $\tau$, số bước trung bình, lý do dừng tăng bước, sai số cực đại, sai số trung bình
Bảng 4.4 Thuật toán thích nghi dùng 4 bộ kiểm soát ($E_d, E_s, E_{d2}, E_{s2}$) Dung sai $\tau$, số bước trung bình, lý do dừng tăng bước, sai số cực đại, sai số trung bình
Bảng 4.5 Phương pháp Milstein bước cố định (Test Problem 2) Kích thước bước $h_{fix}$, sai số tuyệt đối cực đại trung bình, sai số tại $t = 1$
Bảng 4.6 Thuật toán thích nghi: Không cứng, Hiện (Test Problem 2) Dung sai $\tau$, sai số cực đại, sai số tại $t = 1$, số bước thực hiện, số bước chấp nhận/từ chối
Bảng 4.7 Thuật toán thích nghi: Không cứng, Cân bằng (Test Problem 2) Dung sai $\tau$, sai số cực đại, sai số tại $t = 1$, số bước thực hiện, số bước chấp nhận/từ chối
Bảng 4.8 Thuật toán thích nghi: Cứng, Hiện (Test Problem 2) Dung sai $\tau$, sai số cực đại, sai số tại $t = 1$, số bước thực hiện, số bước chấp nhận/từ chối
Bảng 4.9 Thuật toán thích nghi: Cứng, Cân bằng (Test Problem 2) Dung sai $\tau$, sai số cực đại, sai số tại $t = 1$, số bước thực hiện, số bước chấp nhận/từ chối
Bảng 4.10 Tần suất chấm dứt tăng kích thước bước của từng bộ kiểm soát Số bước mà mỗi tiêu chí kiểm soát ($E_d, E_s, E_{d2}, E_{s2}$) ra quyết định dừng tăng bước

Chương 5: Kết luận (Chapter 5: Conclusions)

Chương 5 tóm tắt các phát hiện của công trình theo hai nhánh:

  • Mục 5.1 tổng kết mô hình thị trường tài chính: Lớp mô hình tác tử dựa trên ngưỡng giải thích được sự xuất hiện của các đặc tính phi Gauss mà không cần các giả định gượng ép; việc phân tách nguyên nhân giữa đuôi dày và cụm biến động cung cấp cơ sở để thiết lập các mô hình định lượng chính xác hơn trong tương lai.
  • Mục 5.2 tổng kết phương pháp số cho SDE: Khẳng định tính khả thi của các thuật toán thích nghi bước thời gian dựa trên kiểm soát sai số kép và kiểm soát bổ sung, chứng minh sự vượt trội về độ ổn định và hiệu năng tính toán so với các phương pháp bước cố định.

Kết quả và những đóng góp mới

Các đóng góp mới về mặt lý luận và thực tiễn được luận án xác định bao gồm:

  1. Đóng góp trong mô hình hóa tài chính:
    • Phát triển một lớp mô hình thị trường tài chính dựa trên tác tử với cơ chế ngưỡng đơn giản, phản ánh thỏa đáng tính duy lý giới hạn và các yếu tố tâm lý hành vi của nhà đầu tư dài hạn.
    • Mô hình có số lượng tham số cơ bản rất nhỏ và dễ ước tính. Tại các giới hạn tham số cụ thể, mô hình quy tụ về mô hình hoạt động theo giả thuyết thị trường hiệu quả (EMH), cho phép phân tích rõ ràng mối quan hệ nhân quả giữa các vi phạm giả định EMH và sự hình thành các thuộc tính thống kê phi Gauss.
    • Tái hiện thành công các đặc trưng thực nghiệm cốt lõi của thị trường thực tế trong mô phỏng chu kỳ 40 năm: sự vắng mặt của tự tương quan trong lợi suất, phân phối đuôi dày tuân theo luật lũy thừa, hiện tượng cụm biến động với tốc độ suy giảm tự tương quan chậm, và tính bất đối xứng lãi/lỗ giữa thị trường giá lên và giá xuống.
    • Phát hiện ra rằng hai sai lệch phi Gauss lớn nhất—đuôi dày và cụm biến động—có thể xuất phát từ hai nguyên nhân cơ chế nội sinh khác nhau, tạo tiền đề lý thuyết cho các nghiên cứu mô hình hóa tiếp theo.
  2. Đóng góp trong giải số phương trình vi phân ngẫu nhiên:
    • Chứng minh về mặt lý thuyết rằng thuật toán thích nghi bước thời gian dựa trên chiến lược kiểm soát sai số kép bảo toàn tính ổn định trung bình bình phương trên một lớp bài toán kiểm tra tuyến tính có nhiễu nhân.
    • Triển khai thành công cơ chế từ chối bước tính thử nghiệm ngay sau khi sinh số gia Brown trước khi phải đánh giá các hệ số trôi và khuếch tán tiếp theo, giúp tiết kiệm chi phí tính toán.
    • Mở rộng thuật toán thích nghi sang phương pháp Milstein cân bằng (Balanced Milstein Method) dạng ẩn, giải quyết hiệu quả sự mất ổn định số khi tích phân các phương trình SDE cứng.
    • Đề xuất và thử nghiệm hệ thống 4 bộ kiểm soát sai số ngẫu nhiên ($E_d, E_s, E_{d2}, E_{s2}$), chứng minh tính khả thi của việc khắc phục sự cố sụp đổ ước lượng sai số khuếch tán trên các phương trình SDE có nhiễu cộng.

Kiến nghị và giải pháp kỹ thuật đề xuất: Luận án khuyến nghị sử dụng phương pháp Milstein cân bằng kết hợp chiến lược kiểm soát sai số thích nghi khi giải các bài toán SDE có tính chất cứng trong thực tế; đồng thời đề xuất tiếp tục nghiên cứu các cơ chế lựa chọn bước thời gian trực tiếp để thay thế cho quy trình tăng bước thử nghiệm ngẫu nhiên.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu:

  • Trong phần tài chính toán học: Mô hình chưa xây dựng được nghiệm giải tích dạng đóng (closed-form solutions), chủ yếu vận hành như một mô hình mô phỏng tính toán vi mô. Các nhà giao dịch ngắn hạn (nhà giao dịch nhiễu, nhà kinh doanh chênh lệch giá) không được mô hình hóa thành các thực thể độc lập mà chỉ được tích hợp gián tiếp thông qua tác động lũy kế lên thị trường.
  • Trong phần giải số SDE: Toàn bộ phân tích lý thuyết và thuật toán thích nghi mới chỉ áp dụng cho các SDE vô hướng được kích thích bởi một quá trình ngẫu nhiên vô hướng đơn lẻ (single stochastic scalar process); chưa mở rộng cho hệ phương trình nhiều chiều hoặc SDE có nhiều nguồn nhiễu độc lập. Quy trình chọn bước thời gian thích nghi trong thử nghiệm đòi hỏi chia nhỏ chuyển động Brown thành các bước cực nhỏ, dẫn đến số lần đánh giá chuyển động Brown còn rất lớn, chưa thể ứng dụng trực tiếp như một bộ phần mềm thương mại tối ưu hoàn chỉnh.

Hướng nghiên cứu tiếp theo:

  • Mở rộng các thuật toán thích nghi bước thời gian kiểm soát sai số kép và phương pháp Milstein cân bằng cho các hệ SDE nhiều chiều và SDE có nhiều quá trình ngẫu nhiên Brown độc lập.
  • Thiết kế các chiến lược lựa chọn bước thời gian hiệu quả và có tính hệ thống hơn nhằm giảm thiểu số lượng mẫu chuyển động Brown cần sinh ra.
  • Tiếp tục mở rộng mô hình thị trường dựa trên tác tử để phân lập sâu hơn các yếu tố cấu trúc gây ra hiện tượng đuôi dày và cụm biến động.

Giá trị tham khảo

Luận án có giá trị tham khảo chuyên môn đối với các đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Toán ứng dụng, Khoa học tính toán và Tài chính định lượng (Quantitative Finance): Chương 2 cung cấp một tài liệu tổng quan toàn diện và có hệ thống về lịch sử phát triển của các mô hình thị trường vi mô (từ mô hình BPS, Cont-Bouchaud, Minority Game đến các mô hình của Thomas Lux) và các phương pháp giải số thích nghi cho SDE.
  • Các nhà nghiên cứu và phát triển thuật toán số: Mục 3.2, Mục 4.2 và các Bảng từ 4.3 đến 4.10 cung cấp chi tiết thiết kế thuật toán kiểm soát sai số kép, phương pháp Milstein cân bằng và các dữ liệu kiểm tra nghiệm số chi tiết trên các bài toán cứng và không cứng, rất hữu ích cho việc xây dựng các thư viện tính toán khoa học cho SDE.
  • Các nhà phân tích định lượng tại các tổ chức tài chính: Cơ chế mô hình ngưỡng hành vi tại Mục 3.1 và phân tích tái hiện các đặc trưng thực nghiệm tại Mục 4.1 cung cấp góc nhìn tính toán thực tế về sự hình thành các rủi ro hệ thống và biến động phi Gauss trên thị trường tài chính.

Câu hỏi thường gặp

  1. Giả thuyết thị trường hiệu quả (EMH) gặp phải những hạn chế thực nghiệm cốt lõi nào khiến tác giả phải xây dựng mô hình thay thế? Trả lời: Theo tác giả, EMH giả định giá tài sản biến động ngẫu nhiên theo phân phối chuẩn Gauss, thị trường không có cơ hội kinh doanh chênh lệch giá kéo dài và các tác tử có kỳ vọng duy lý. Tuy nhiên, dữ liệu thực nghiệm từ thời Mandelbrot (1963) và Fama (1965) đã chỉ ra rằng lợi suất thực tế có quá nhiều biến động nhỏ và biến động lớn (đuôi dày), dẫn đến hiện tượng leptokurtosis (độ nhọn vượt mức 3 của phân phối chuẩn) và hiện tượng cụm biến động mà EMH không thể giải thích được.

  2. Mô hình thị trường dựa trên tác tử của Timothy L. Seaman có điểm gì khác biệt so với mô hình Lux-Marchesi (1999) và các mô hình Trò chơi thiểu số (Minority Game)? Trả lời: Khác với mô hình Lux-Marchesi vốn đòi hỏi tới 10 tham số hiệu chuẩn phức tạp và phân chia tác tử thành ba nhóm với các xác suất chuyển đổi ngẫu nhiên kiểu Poisson, mô hình của Seaman chỉ sử dụng một số lượng tham số cơ bản rất nhỏ và dễ ước tính. Khác với Trò chơi thiểu số (vốn ban đầu không có biến số đo lường giá cụ thể), mô hình của Seaman xác định trực tiếp giá tài sản, tập trung mô phỏng trạng thái nội tại của các tác tử dài hạn có hành vi ngưỡng và có thể quy tụ trực tiếp về mô hình EMH tại các giới hạn tham số nhất định.

  3. Chiến lược kiểm soát sai số kép (dual error control) của Harbir Lamba (2003) vận hành như thế nào trong giải số SDE? Trả lời: Chiến lược này áp dụng hai ước lượng sai số cục bộ độc lập cho phương pháp Milstein bậc 1: ước lượng sai số $E_d$ dành cho thành phần trôi (deterministic drift) và ước lượng sai số $E_s$ dành cho thành phần khuếch tán (stochastic diffusion). Khi hệ số khuếch tán tiến về 0, thuật toán tự động vận hành tương tự như một bộ giải thích nghi ODE thông thường. Điểm đặc thù là thuật toán cho phép từ chối bước tính ngay sau khi tính số gia Brown nếu sai số ngẫu nhiên vượt ngưỡng, trước khi phải thực hiện các phép đánh giá hàm phức tạp tiếp theo.

  4. Tại sao phương pháp Milstein cân bằng (Balanced Milstein Method) lại được đưa vào nghiên cứu bên cạnh phương pháp Milstein hiện? Trả lời: Phương pháp Milstein hiện (Explicit Milstein Method) bị giới hạn nghiêm trọng về kích thước bước thời gian khi giải các phương trình SDE cứng (stiff SDEs) do yêu cầu ổn định số. Luận án đã tích hợp phương pháp Milstein cân bằng—một biến thể bán ẩn—vào khung kiểm soát sai số thích nghi nhằm duy trì tính ổn định và cho phép thuật toán thực hiện các bước thời gian lớn hơn trên các bài toán SDE cứng, như được kiểm chứng qua các kết quả số tại Bảng 4.8 và Bảng 4.9.

  5. Luận án đã đề xuất giải pháp gì khi bộ kiểm soát sai số khuếch tán bị sụp đổ trên các phương trình SDE có nhiễu cộng? Trả lời: Khi giải các phương trình SDE có nhiễu cộng, đạo hàm của hệ số khuếch tán triệt tiêu khiến ước lượng sai số khuếch tán $E_s$ thông thường không còn tác dụng. Tác giả đã đề xuất và thử nghiệm bổ sung hai bộ kiểm soát sai số ngẫu nhiên mới, nâng tổng số tiêu chí kiểm soát lên 4 bộ ($E_d, E_s, E_{d2}, E_{s2}$), giúp thuật toán duy trì khả năng kiểm soát sai số và ổn định kích thước bước khi tích phân SDE có nhiễu cộng.

Kết luận

Luận án của Timothy L. Seaman đã giải quyết thành công hai bài toán mô hình hóa và tính toán khoa học trong tài chính toán học và giải số phương trình vi phân ngẫu nhiên. Công trình đã xây dựng một mô hình thị trường dựa trên tác tử với cơ chế ngưỡng đơn giản, tái hiện đầy đủ các đặc trưng thực nghiệm phi Gauss của thị trường tài chính và chỉ ra sự khác biệt về mặt cơ chế giữa hiện tượng đuôi dày và hiện tượng cụm biến động. Đồng thời, luận án đã chứng minh tính ổn định trung bình bình phương của thuật toán thích nghi bước thời gian kiểm soát sai số kép, mở rộng thành công sang phương pháp Milstein cân bằng cho các phương trình SDE cứng và đề xuất các bộ kiểm soát sai số ngẫu nhiên bổ sung hiệu quả.