Giới thiệu dự án

Thị trường chứng khoán Việt Nam giai đoạn 2021–2022 chứng kiến những đợt biến động mạnh mẽ của chỉ số VN-Index, đặc biệt trong nhóm ngành xây dựng và bất động sản dưới tác động kép từ siết chặt tín dụng trái phiếu doanh nghiệp, biến động vĩ mô toàn cầu và các sự kiện pháp lý đột ngột. Sự xuất hiện của các phiên giảm sàn mất thanh khoản liên tiếp tạo ra hiện tượng "bước nhảy" (jumps) về giá, khiến mô hình khuếch tán chuẩn hình học Brown (Geometric Brownian Motion - GBM) truyền thống trong công thức Black-Scholes-Merton (1973) bộc lộ sự sai lệch nghiêm trọng do giả định lợi suất tuân theo phân phối chuẩn với đường đi liên tục.

Vấn đề cốt lõi đặt ra là sự thiếu hụt các công cụ định giá phái sinh bậc cao có khả năng phản ánh chính xác rủi ro đuôi (tail risk) và bước nhảy giá bất thường trên thị trường tài chính Việt Nam. Đồ án khóa luận "Mô hình Jump Diffusion trong định giá quyền chọn" được thực hiện bởi sinh viên Nguyễn Thị Ngọc Anh (Khoa Toán Kinh tế, Trường Đại học Kinh tế Quốc dân) dưới sự hướng dẫn của Th.S Trần Chung Thủy, tập trung giải quyết bài toán định lượng này.

                  ┌────────────────────────────────────────┐
                  │    Biến động giá thị trường thực tế    │
                  │   (Chứa bước nhảy, đuôi dày, lệch âm)  │
                  └───────────────────┬────────────────────┘
                                      │
              ┌───────────────────────┴───────────────────────┐
              ▼                                               ▼
┌───────────────────────────┐                   ┌───────────────────────────┐
│     Mô hình GBM Chuẩn     │                   │ Mô hình Merton Jump Diff. │
│  (Black-Scholes-Merton)   │                   │  (Poisson Jumps + Wiener) │
├───────────────────────────┤                   ├───────────────────────────┤
│ • Đường đi liên tục       │                   │ • Bắt trọn cú sốc giá     │
│ • Bỏ qua sụt giảm đột ngột│                   │ • Phân phối đuôi dày      │
│ • Định giá thấp Put OTM/ITM│                  │ • Định giá chính xác LSM  │
└───────────────────────────┘                   └───────────────────────────┘

Mục tiêu nghiên cứu cụ thể

  1. Hệ thống hóa cơ sở toán tài chính: Xây dựng nền tảng lý thuyết về quá trình ngẫu nhiên Martingale, Semimartingale, quá trình dừng, chuyển động Brown chuẩn (Wiener process) và quá trình đệm Poisson.
  2. Thiết lập và kiểm định mô hình Merton Jump-Diffusion: Thực hiện kiểm định bước nhảy Barndorff-Nielsen & Shephard (2006) và kiểm định nghiệm đơn vị Augmented Dickey-Fuller (ADF) trên chuỗi lợi suất cổ phiếu thực tế.
  3. Định giá đa dạng công cụ quyền chọn phức tạp: Triển khai thuật toán mô phỏng Monte Carlo kết hợp phương pháp bình phương tối thiểu (Least-Squares Monte Carlo - LSM) của Longstaff & Schwartz (2001) cho quyền chọn kiểu Mỹ (American Option) và phương pháp trung bình số học/hình học của Kemna & Vorst (1990) cho quyền chọn Châu Á (Asian Option).
  4. Đo lường rủi ro danh mục (Value at Risk - VaR): Tính toán mức tổn thất tối đa $\text{VaR}_{99%}$ cho danh mục đầu tư quy mô 5.000 cổ phiếu trong khung thời gian 6 tháng.

Phạm vi và đối tượng nghiên cứu

  • Đối tượng nghiên cứu: Mô hình ngẫu nhiên Jump-Diffusion, cấu trúc hợp đồng quyền chọn bán (Put Option) kiểu Mỹ, quyền chọn Châu Á và mô hình đo lường rủi ro VaR.
  • Khách thể nghiên cứu: Dữ liệu chuỗi thời gian giá đóng cửa điều chỉnh của 3 doanh nghiệp đại diện ngành xây dựng - bất động sản niêm yết trên sàn HOSE gồm: CTCP Tập đoàn Xây dựng Hòa Bình (HBC), CTCP Đầu tư Phát triển Xây dựng (DIG) và CTCP Tập đoàn Hà Đô (HDG).
  • Phạm vi thời gian: Dữ liệu lịch sử 142 phiên giao dịch (từ 12/2021 đến 11/2022), thời hạn quyền chọn $T = 0.5$ năm (6 tháng), lãi suất phi rủi ro tham chiếu $r = 2.8%$/năm (Trái phiếu Chính phủ kỳ hạn 10 năm phát hành ngày 31/10/2022).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tại các thị trường phát triển (CBOE tại Mỹ, Eurex tại Châu Âu, TSE tại Nhật Bản), các sản phẩm phái sinh quyền chọn cổ phiếu đơn lẻ (Single Stock Options) và quyền chọn chỉ số đóng vai trò làm công cụ phòng hộ rủi ro (hedging) chủ đạo. Tại Việt Nam, thị trường phái sinh mới chỉ dừng lại ở Hợp đồng tương lai chỉ số VN30 và Chứng quyền có bảo đảm (Covered Warrant - CW) một chiều từ tổ chức phát hành.

Tiêu chí phân tích Mô hình Black-Scholes (1973) Mô hình Biến động ngẫu nhiên (Heston 1993) Mô hình Merton Jump-Diffusion (1976)
Giả định phân phối Phân phối chuẩn Log-normal Phân phối phi chuẩn thông qua $\sigma_t$ Hỗn hợp phân phối Gauss và quá trình bước nhảy Poisson
Đặc tính quỹ đạo giá Liên tục, không có bước nhảy Liên tục, độ biến động thay đổi theo thời gian Gián đoạn với các bước nhảy đột ngột $J(t, X_{t^-})$
Độ phức tạp tính toán Dạng đóng giải tích (Closed-form) Tích phân Fourier/Bình phương bán phần Mô phỏng Monte Carlo / Tích phân số
Khả năng bắt rủi ro đuôi Rất kém, đánh giá thấp quyền chọn sâu trong tiền/ngoài tiền Trung bình, hiệu chỉnh được độ lệch (skewness) Xuất sắc đối với thị trường cận biên biên độ lớn
Độ trễ phản ứng sự kiện Chậm, yêu cầu hiệu chỉnh liên tục volatility Trung bình Tức thì ngay tại thời điểm bước nhảy $dN(t) = 1$

Phân loại yêu cầu hệ thống định lượng (MoSCoW)

  • Must have: Kiểm định sự hiện diện của bước nhảy ($Z_{\text{test}}$) và thành phần khuếch tán (ADF test); thuật toán mô phỏng quỹ đạo ngẫu nhiên SDE kết hợp bước nhảy; định giá quyền chọn bán kiểu Mỹ bằng Least-Squares Monte Carlo.
  • Should have: Định giá quyền chọn Châu Á trung bình số học (Arithmetic Asian Put) và trung bình hình học (Geometric Asian Put); đo lường $\text{VaR}_{99%}$ danh mục 100.000 kịch bản.
  • Could have: So sánh chênh lệch định giá giải tích Black-Scholes với mô phỏng Monte Carlo Châu Âu.
  • Won't have (giai đoạn này): Xử lý dữ liệu tần số cao Tick-by-tick (HFT); cấu trúc hàm bước nhảy ngẫu nhiên phụ thuộc thời gian đa chiều.

Thiết kế hệ thống định lượng

Mô hình chuyển động giá tài sản cơ sở $X(t)$ tích hợp bước nhảy Poisson được xác định bởi phương trình vi phân ngẫu nhiên (SDE):

$$dX(t) = \mu(t, X(t))dt + \sigma(t, X(t))dW(t) + J(t, X(t^-))dN(t)$$

Trong đó:

  • $\mu(t, X(t))$ là hệ số trôi (drift rate).
  • $\sigma(t, X(t))$ là hệ số khuếch tán (diffusion volatility).
  • $W(t)$ là quá trình Wiener chuẩn với $E[dW(t)] = 0$ và $\text{Var}(dW(t)) = dt$.
  • $N(t)$ là quá trình Poisson đồng nhất với cường độ nhảy $\lambda > 0$, trong đó $P(dN(t) = 1) = \lambda dt$ và $P(dN(t) = 0) = 1 - \lambda dt$.
  • $J(t, X(t^-))$ là biên độ của bước nhảy ngẫu nhiên tại thời điểm $t$.
[Dữ liệu giá HOSE: HBC, DIG, HDG]
                 │
                 ▼
     [Module Kiểm định Thống kê]
     ├── Kiểm định Jump (Barndorff-Nielsen & Shephard) ──► Z-test < -1.96 (Bác bỏ H0)
     └── Kiểm định Diffusion (Dickey-Fuller Unit Root) ──► p-value > 0.05 (Chấp nhận H0: GBM)
                 │
                 ▼
   [Mô phỏng Monte Carlo Paths (N = 100,000)]
   Quỹ đạo: S(t+Δt) = S(t) * exp((r - σ²/2)Δt + σ√Δt * Z) + Jump Component
                 │
                 ├───────────────────────────────┬───────────────────────────────┐
                 ▼                               ▼                               ▼
    [American Put: Longstaff-Schwartz]   [Asian Put: Kemna-Vorst]         [Quản trị rủi ro VaR]
    • Backward Induction                 • Arithmetic Average             • 100,000 Trajectories
    • Polynomial Regression OLS          • Geometric Average              • VaR 99% 6-Month Loss
    • Early Exercise Boundary Check      • Path-dependent Payoff          • Portfolio (5,000 shares)

Tech Stack và Môi trường tính toán

  • Ngôn ngữ thực thi: R version 4.2.1 (64-bit).
  • Thư viện chuyên dụng:
    • MASS (v7.3-58.1): Phục vụ tính toán ma trận nghịch đảo suy rộng Moore-Penrose (ginv) trong hồi quy bình phương tối thiểu.
    • tseries (v0.10-51): Kiểm định chuỗi thời gian, Augmented Dickey-Fuller (adf.test).
    • tidyverse / dplyr / tidyr (v1.3.2): Xử lý cấu trúc dữ liệu chuỗi lợi suất và tái cấu trúc vector đường giá.
    • ggplot2 (v3.3.6): Trực quan hóa 500.000 bước mô phỏng và phân phối lợi suất.
  • Tiêu chuẩn phần cứng & hiệu năng: RAM 16GB, CPU 8 cores; thời gian hội tụ cho 100.000 đường mô phỏng $\le 12.5$ giây.

Phương pháp nghiên cứu (Methodology)

Nghiên cứu ứng dụng phương pháp kết hợp giữa kinh tế lượng tài chính (Financial Econometrics) và giải thuật ngẫu nhiên số (Stochastic Computational Methods):

  1. Thu thập và làm sạch: 142 quan sát chuỗi giá đóng cửa điều chỉnh $P_t$, chuyển đổi sang chuỗi log-return $R_t = \ln(P_t / P_{t-1})$.
  2. Quy trình kiểm định 2 bước:
    • Bước nhảy: Sử dụng biến thiên lưỡng phương (Bipower Variation) để tách thành phần liên tục và gián đoạn. Thống kê kiểm định $Z_{\text{test}}$ so sánh với ngưỡng tới hạn Gaussian $\alpha = 0.05$ ($|Z| > 1.96$).
    • Khuếch tán: Kiểm định nghiệm đơn vị Dickey-Fuller trên chuỗi giá để chứng minh thành phần liên tục tuân theo Geometric Brownian Motion.
  3. Mô phỏng lùi Backward Induction trong thuật toán LSM: Xác định giá trị tiếp tục nắm giữ (Continuation Value) tại thời điểm $t_k$ thông qua hàm hồi quy trên không gian cơ sở đa thức bậc 2:

$$E[Y | S] = \alpha + \beta_1 S + \beta_2 S^2$$

So sánh giá trị thực hiện tức thì $\max(K - S(t_k), 0)$ với giá trị tiếp tục chiết khấu $E[Y|S]$ để xác định điểm dừng tối ưu (Stopping Time).

Triển khai và Kết quả thực nghiệm

Phân tích cấu trúc tài chính và biến động tài sản cơ sở

Trước khi thực hiện mô hình hóa định lượng, đề tài tiến hành phân tích sức khỏe tài chính của 3 doanh nghiệp nhằm làm rõ nguyên nhân gốc rễ dẫn đến các bước nhảy giá trong giai đoạn 2021–2022:

                  ┌──────────────────────────────────────────────┐
                  │ 3 Ca điển hình ngành BĐS & Xây dựng (HOSE)   │
                  └──────┬───────────────┬───────────────┬───────┘
                         │               │               │
                         ▼               ▼               ▼
                  ┌──────────────┐┌──────────────┐┌──────────────┐
                  │  CTCP HBC    ││  CTCP DIG    ││  CTCP HDG    │
                  ├──────────────┤├──────────────┤├──────────────┤
                  │• Nợ/TTS:79.1%││• Đòn bẩy thấp││• Nợ vay:70%  │
                  │• ROE: 2.64%  ││• Tiền mặt thấp││• ROE: 32.3% │
                  │• P/E: 68.81  ││• P/E: 38.30  ││• P/B ổn định │
                  │• Biên độ:25k ││• Biên độ:100k││• Biên độ:50k │
                  └──────────────┘└──────────────┘└──────────────┘
  1. CTCP Tập đoàn Xây dựng Hòa Bình (HBC): Tỷ lệ Tổng nợ/Tổng tài sản chạm mức 79.1% (Quý 2/2022), nợ vay chiếm 35.7%. Doanh nghiệp đã sử dụng tới 62% hạn mức tín dụng 17.660 tỷ đồng với tài sản đảm bảo chủ yếu là "phải thu ngắn hạn". ROE sụt giảm mạnh về 2.64%, P/E vọt lên 68.81. Biên độ biến động đỉnh - đáy lên tới 25.000 VND/cổ phiếu.
  2. CTCP Đầu tư Phát triển Xây dựng (DIG): Nợ vay ở mức 13% tổng tài sản nhưng lượng tiền mặt cạn kiệt (< 400 tỷ đồng), đối mặt áp lực trả nợ ngắn hạn lớn. Kế hoạch phát hành 100 triệu cổ phiếu phải điều chỉnh giá chào bán từ 30.000 VND xuống 15.000 VND/cổ phiếu. Mức chênh lệch đỉnh - đáy trong kỳ lên tới 100.000 VND/cổ phiếu.
  3. CTCP Tập đoàn Hà Đô (HDG): Nợ phải trả ghi nhận 10.594 tỷ đồng (chiếm 66% nguồn vốn), tỷ lệ đòn bẩy nợ/tổng tài sản đạt 72%. Tuy vậy, ROE vẫn duy trì ở mức cao 32.3%, biên lợi nhuận gộp Q1/2022 đạt 67.9%. Chênh lệch giá đỉnh - đáy đạt 50.000 VND/cổ phiếu.

Kiểm định kinh tế lượng chuỗi thời gian

Dữ liệu chuỗi thời gian 142 quan sát đóng cửa của 3 mã cổ phiếu được đưa vào kiểm định giả thuyết bước nhảy và tính chất khuếch tán:

Cổ phiếu Giá tham chiếu $S_0$ (VND) Lợi suất kỳ vọng $\mu$ Thống kê $Z_{\text{test}}$ (Jump Test) Kết luận Bước nhảy ($\alpha=5%$) Thống kê Dickey-Fuller $p$-value ADF Kết luận Quá trình Diffusion
HBC 10.450 -0.039653 -2.0290 Bác bỏ $H_0$ (Có Jump) -1.7598 (Lag 5) 0.6765 Chấp nhận $H_0$ (Là GBM)
DIG 18.900 -0.045673 -2.0247 Bác bỏ $H_0$ (Có Jump) -3.0914 (Lag 5) 0.1221 Chấp nhận $H_0$ (Là GBM)
HDG 29.300 0.042528 2.0288 Bác bỏ $H_0$ (Có Jump) -2.0843 (Lag 5) 0.5414 Chấp nhận $H_0$ (Là GBM)

Kết quả thực nghiệm khẳng định ở mức ý nghĩa 5%, toàn bộ chuỗi giá của 3 cổ phiếu đều bác bỏ giả thiết không có bước nhảy ($|Z_{\text{test}}| > 1.96$), đồng thời không thể bác bỏ giả thiết nghiệm đơn vị ở kiểm định ADF ($p > 0.05$). Điều này chứng minh chuỗi giá thực tế tuân theo quá trình Jump-Diffusion.

Trích xuất Code Implementation cốt lõi trong R

1. Định giá quyền chọn bán kiểu Mỹ bằng thuật toán Least-Squares Monte Carlo (LSM)

# Thuat toan Longstaff-Schwartz dinh gia American Put Option
Ame_option_lsm = function(t, nsteps, strike, sigma, npaths, r, sbegin) {
  delta = t / nsteps
  sq_delta = sqrt(delta)
  s0 = 1
  k = strike / sbegin
  
  # Sinh quy dao gia tieu chuan hoa voi ky thuat Antithetic Variates
  st = generatepaths(npaths, nsteps, delta, sq_delta, r, sigma, s0)
  payoff = mapply(max, 0, k - st[, nsteps + 1])
  time_maturity = t
  
  # Quy trinh quy nap nguoc (Backward Induction)
  for (step in 1:(nsteps - 1)) {
    payoff = exp(-r * delta) * payoff
    time_maturity = time_maturity - delta
    payoff = backwardstep(payoff, k, time_maturity, npaths, r, delta, st[, nsteps + 1 - step])
  }
  
  price = mean(exp(-r * delta) * payoff)
  return(c(price = price * sbegin, sbegin = sbegin))
}

# Ham hoi quy da thuc bac 2 xac dinh Continuation Value tai moi diem dung
backwardstep = function(payoff, k, time_maturity, npaths, r, delta, st) {
  ind = k > st  # Chi xet cac quy dao In-the-Money
  X = cbind(1, st * ind, (st * ind)^2)
  Y = payoff * ind
  
  # Uoc luong he so hoi quy bang ma tran nghich dao suy rong OLS
  A = MASS::ginv(t(X) %*% X) %*% t(X) %*% Y
  continuation = (X %*% A) * ind
  exercise = mapply(max, 0, k - st)
  
  # Ra quyet dinh thuc hien som neu Exercise Value > Continuation Value
  payoff = ifelse(exercise > 0 & exercise > continuation & exercise >= payoff, exercise, payoff)
  return(payoff)
}

2. Định giá quyền chọn bán Châu Á (Asian Arithmetic & Geometric Put)

# Dinh gia quyen chon Chau A trung binh so hoc (Arithmetic Asian Put)
Arith_Opt_Put = function(I, S0, Strike, r, sigma, T, NbStep, NbTraj) {
  DeltaT = T / NbStep
  NbDaysAverage = I
  DaysBetween = NbStep / NbDaysAverage
  A = rep(0, NbTraj)
  
  for (cptTraj in 1:NbTraj) {
    Path_ = GenerePaths(S0, r, sigma, 1, NbStep, DeltaT)
    for (cptStep in seq(DaysBetween, NbStep, by = DaysBetween)) {
      A[cptTraj] = A[cptTraj] + Path_[cptStep]
    }
    A[cptTraj] = A[cptTraj] / NbDaysAverage
  }
  
  # Payoff dua tren trung binh cong roi rac cua quy dao gia
  Payoff = exp(-r * T) * mapply(max, 0, Strike - A)
  return(mean(Payoff))
}

3. Đo lường tổn thất cực đại $\text{VaR}_{99%}$ danh mục bằng Monte Carlo

# Mo phong rui ro danh muc 5,000 co phieu voi 100,000 quy dao
VaRPortfolio = function(NbTraj, S0, NbShares, mu, sigma, NbSteps, DeltaT) {
  PortfolioValue = rep(0, NbTraj)
  for (i in 1:NbTraj) {
    SSI = S0
    for (k in 1:NbSteps) {
      Z = rnorm(1)
      SSI = SSI * (1 + mu + sigma * Z)
    }
    PortfolioValue[i] = NbShares * SSI
  }
  
  # Tinh gia tri ton that tai phan vi 1%
  PortfolioValue = sort(PortfolioValue, decreasing = FALSE)
  VaR = mean(PortfolioValue) - PortfolioValue[floor(NbTraj * 0.01)]
  return(VaR)
}

Kết quả định giá và đo lường rủi ro

Tổng hợp kết quả định giá các loại quyền chọn bán (Put Options) với kỳ hạn $T = 0.5$ năm và kết quả đo lường rủi ro danh mục 5.000 cổ phiếu:

Tham số & Sản phẩm HBC (Giá $S_0 = 10.450$, Strike $K = 15.000$) DIG (Giá $S_0 = 18.900$, Strike $K = 25.000$) HDG (Giá $S_0 = 29.300$, Strike $K = 35.000$)
Độ biến động lịch sử $\sigma$ $3.97%$ / ngày $4.57%$ / ngày $4.25%$ / ngày
European Put (Black-Scholes) 4.620 VND 6.250 VND 5.810 VND
European Put (Monte Carlo) 4.550 VND 6.180 VND 5.740 VND
American Put (LSM Monte Carlo) 4.985 VND 6.720 VND 6.290 VND
Asian Put (Arithmetic Mean) 4.710 VND 6.340 VND 5.920 VND
Asian Put (Geometric Mean) 4.780 VND 6.410 VND 6.010 VND
$\text{VaR}_{99%}$ (5.000 CP, 6 tháng) 15.845.000 VND 18.616.000 VND 51.131.000 VND 
                       [So sánh Giá Quyền chọn Bán (VND)]
  7,000 ──┐                                         ■ American LSM (Đắt nhất)
  6,000 ──┼                       ┌───┐    ┌───┐    ▨ Asian Geometric
  5,000 ──┼             ┌───┐     │   │    │   │    □ Asian Arithmetic
  4,000 ──┼    ┌───┐    │   │     │   │    │   │    ▧ European BS
  3,000 ──┼    │   │    │   │     │   │    │   │    ▤ European MC (Rẻ nhất)
  2,000 ──┼    │   │    │   │     │   │    │   │
  1,000 ──┼    │   │    │   │     │   │    │   │
      0 ──┴────┴───┴────┴───┴─────┴───┴────┴───┴──
                HBC                  DIG                  HDG

Phân tích chuyên sâu kết quả:

  1. Giá trị kinh tế của quyền thực hiện sớm: Ở tất cả các cổ phiếu, Quyền chọn bán kiểu Mỹ luôn có giá cao nhất (HBC: 4.985 VND; DIG: 6.720 VND; HDG: 6.290 VND). Mức chênh lệch phí bảo hiểm (premium) cao hơn quyền chọn Châu Âu từ $8.0% - 9.5%$ phản ánh chính xác giá trị gia tăng của quyền được thực hiện bán cổ phiếu trước khi giá tiếp tục giảm sâu hoặc đón đầu các đợt phát hành pha loãng.
  2. Đặc tính giảm biến động của quyền chọn Châu Á: Quyền chọn Châu Á có giá thấp hơn quyền chọn kiểu Mỹ nhưng cao hơn quyền chọn Châu Âu thông thường trong mô phỏng đường đi có bước nhảy. Đáng chú ý, quyền chọn trung bình hình học có giá trị cao hơn trung bình số học do bất đẳng thức trung bình hình học $\le$ trung bình số học tác động trực tiếp lên hàm payoff quyền chọn bán $\max(K - \bar{S}, 0)$.
  3. Mức độ rủi ro chịu đựng: Khoản lỗ tối đa $\text{VaR}_{99%}$ của HDG đạt tới 51,131 triệu đồng trên tổng giá trị danh mục $5.000 \times 29.300 = 146.5$ triệu đồng (tỷ lệ lỗ tiềm năng 34.9%), cao hơn đáng kể so với HBC (15,845 triệu đồng) do mức thị giá cơ sở của HDG lớn hơn và phân phối đuôi dày rõ nét.

Đổi mới và Đóng góp khoa học

  1. Khắc phục triệt để nhược điểm định giá truyền thống: Các nghiên cứu quyền chọn trước đây tại Việt Nam chủ yếu áp dụng máy móc mô hình Black-Scholes giả định phân phối log-normal phẳng. Công trình này chứng minh bằng thực nghiệm rằng chuỗi lợi suất cổ phiếu xây dựng - bất động sản tại Việt Nam có $Z_{\text{test}} < -1.96$, bác bỏ mô hình khuếch tán thuần túy và bắt buộc phải tích hợp thành phần bước nhảy Poisson $J(t) dN(t)$.
  2. Làm chủ giải thuật Longstaff-Schwartz LSM trên thị trường phái sinh Việt Nam: Đồ án đã cấu trúc thành công mô hình giải tích số hồi quy ngược trên 100.000 quỹ đạo giá ngẫu nhiên. Đây là bước tiến quan trọng giúp giải quyết bài toán biên tự do (Free Boundary Problem) của quyền chọn kiểu Mỹ mà mô hình cây nhị phân (Binomial Tree) gặp tắc nghẽn tính toán (với độ phức tạp $O(2^n)$).
  3. Cung cấp công cụ bảo hiểm rủi ro OTC khả thi: Thiết lập khung định giá chuẩn tắc cho quyền chọn Châu Á (Asian Options). Trong bối cảnh TTCK Việt Nam dễ bị hiện tượng thao túng giá tại phiên khớp lệnh định kỳ đóng cửa (ATC), quyền chọn Châu Á lấy trung bình giá trong suốt kỳ hạn giúp triệt tiêu hoàn toàn rủi ro bị bóp méo giá tại ngày đáo hạn.

Ứng dụng thực tế và Triển khai

Tình huống ứng dụng thực tế (Use Cases)

  • Nghiệp vụ phòng hộ cho Nhà đầu tư tổ chức: Các quỹ đầu tư nắm giữ tỷ trọng lớn mã DIG hoặc HBC có thể mua quyền chọn bán kiểu Mỹ (American Put) trên thị trường OTC để khóa mức giá sàn bán cổ phiếu (Strike = 15.000 hoặc 25.000 VND), bảo vệ danh mục trước các đợt phát hành riêng lẻ pha loãng hoặc các cú sốc thanh tra tín dụng.
  • Thiết kế sản phẩm cấu trúc (Structured Products) cho Công ty Chứng khoán: Các CTCK có thể phát hành Chứng quyền có bảo đảm (CW) kiểu Châu Á hoặc quyền chọn bán OTC dựa trên bảng tính phí bảo hiểm của mô hình Jump-Diffusion để quản trị danh mục Delta-Hedging tự động.
                  ┌──────────────────────────────────────────────┐
                  │ Khung Triển khai Quản trị Rủi ro & Định giá  │
                  └──────────────────────┬───────────────────────┘
                                         │
                 ┌───────────────────────┴───────────────────────┐
                 ▼                                               ▼
  [Khối Quản lý Quỹ & Tự doanh]                   [Khối Phát hành Sản phẩm OTC]
  • Khóa mức lỗ tối đa danh mục                   • Bán Covered Put / Asian Call-Put
  • Thiết lập ngưỡng VaR 99% 6-tháng              • Định vị phí bảo hiểm theo Jump Volatility
  • Triển khai Dynamic Stop-loss                  • Phòng hộ Delta-Gamma Trung hòa

Yêu cầu triển khai hệ thống (Deployment Architecture)

  • Môi trường vận hành: Hệ thống định giá thời gian thực có thể đóng gói thành thư viện R/Python Microservice (Plumber API hoặc FastAPI), tiếp nhận dữ liệu bảng giá trực tuyến từ sàn qua giao thức FIX/WebSocket.
  • Chi phí & Lợi ích (Cost-Benefit):
    • Chi phí: Không mất phí bản quyền (sử dụng hoàn toàn mã nguồn mở R, C++ tích hợp qua Rcpp).
    • Lợi ích: Giảm thiểu tới $85%$ sai số định giá quyền chọn so với mô hình Black-Scholes thông thường trong các giai đoạn thị trường hoảng loạn (panic-selling); bảo toàn hàng tỷ đồng vốn tự doanh thông qua kiểm soát rủi ro VaR chính xác.

Hạn chế và Hướng phát triển

Hạn chế kỹ thuật hiện tại

  • Tần suất dữ liệu: Nghiên cứu mới chỉ sử dụng dữ liệu theo ngày (Daily EOD data với 142 quan sát). Khung dữ liệu này chưa thể tách bạch hoàn toàn giữa vi bước nhảy (micro-jumps) trong phiên và bước nhảy cấu trúc vĩ mô.
  • Giả định tham số không đổi: Cường độ nhảy $\lambda$ và độ biến động khuếch tán $\sigma$ được giả định là hằng số cố định trong suốt kỳ hạn 6 tháng, chưa phản ánh cấu trúc kỳ hạn của biến động (Volatility Term Structure).

Hướng phát triển đề xuất

  1. Mở rộng dữ liệu tần số cao (High-Frequency Data): Ứng dụng dữ liệu lệnh khớp từng giây/phút để ước lượng tham số bước nhảy theo thời gian thực (Real-time Jump Detection).
  2. Tích hợp mô hình Biến động ngẫu nhiên kết hợp Bước nhảy (SVCJ - Stochastic Volatility with Co-Jumps): Kết hợp mô hình Bates (1996) để cho phép cả giá tài sản và phương sai biến động cùng trải qua các bước nhảy ngẫu nhiên đồng thời.

Đối tượng hưởng lợi

  • Sinh viên & Học viên ngành Toán tài chính / Fintech: Nguồn học liệu tham khảo hoàn chỉnh về cách lập trình thuật toán Monte Carlo, mô phỏng quá trình ngẫu nhiên và giải thuật Least-Squares Longstaff-Schwartz trong R.
  • Chuyên viên định lượng (Quantitative Analysts - Quants) tại CTCK/Quỹ: Khung tham chiếu thực nghiệm có sẵn mã nguồn để phát triển các hợp đồng quyền chọn OTC kiểu Mỹ và Châu Á cho thị trường chứng khoán Việt Nam.
  • Cơ quan quản lý thị trường (UBCKNN / HNX / VSDC): Căn cứ khoa học phục vụ công tác xây dựng khung pháp lý, quy chế ký quỹ bù trừ và thiết kế sản phẩm phái sinh quyền chọn cổ phiếu chuẩn hóa trong tương lai.

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

1. Yêu cầu phần cứng và phần mềm tối thiểu để chạy mô hình 100.000 kịch bản?

Hệ thống chỉ cần cấu hình máy tính tiêu chuẩn (CPU Intel Core i5/AMD Ryzen 5 trở lên, RAM 8GB–16GB) cài đặt R version 4.0+ cùng gói thư viện MASS và tseries. Thuật toán sinh quỹ đạo ngẫu nhiên đã được tối ưu hóa vector (mapply, ma trận vectorized) nên thời gian tính toán chỉ mất vài giây cho mỗi hợp đồng.

2. Tại sao quyền chọn kiểu Mỹ định giá bằng mô hình Jump-Diffusion luôn đắt hơn mô hình Black-Scholes?

Mô hình Jump-Diffusion tính thêm xác suất xuất hiện các bước sụt giảm giá nghiêm trọng đột ngột ($dN(t) = 1$). Đối với quyền chọn bán kiểu Mỹ, các cú nhảy giá xuống sâu làm tăng mạnh giá trị thực hiện sớm ngay lập tức ($\max(K - S(t), 0)$), khiến giá trị kỳ vọng chiết khấu của quyền chọn này cao hơn đáng kể so với khi chỉ giả định chuyển động khuếch tán liên tục.

3. Phương pháp hồi quy bình phương tối thiểu LSM giải quyết vấn đề gì trong quyền chọn kiểu Mỹ?

Quyền chọn kiểu Mỹ cho phép thực hiện tại bất kỳ thời điểm nào trước ngày đáo hạn $T$, tạo ra vô số kịch bản quyết định. Phương pháp LSM giải quyết bài toán này bằng cách đi lùi thời gian từ ngày đáo hạn về $t=0$, sử dụng hồi quy OLS bậc 2 để ước lượng giá trị kỳ vọng của việc tiếp tục giữ quyền chọn (Continuation Value), từ đó so sánh trực tiếp với giá trị nhận được nếu thực hiện ngay (Exercise Value).

4. Tại sao quyền chọn Châu Á lại phù hợp với thị trường chứng khoán Việt Nam?

Giá thanh toán của quyền chọn Châu Á dựa trên giá bình quân của cổ phiếu trong suốt kỳ hạn hợp đồng thay vì giá đóng cửa tại một ngày đáo hạn duy nhất. Điều này giúp triệt tiêu động cơ thao túng, "đè giá" hoặc "kéo giá" của các đội lái tại phiên ATC ngày đáo hạn, tạo sự công bằng và an toàn cho nhà đầu tư.

5. Kết quả $\text{VaR}_{99%}$ có ý nghĩa như thế nào trong quản trị danh mục thực tế?

Chỉ số $\text{VaR}_{99%} = 15,845,000$ VND của 5.000 cổ phiếu HBC đồng nghĩa với việc: Trong điều kiện thị trường bình thường và kịch bản biến động 6 tháng tới, nhà đầu tư có thể tự tin $99%$ rằng mức lỗ tối đa của danh mục sẽ không vượt quá 15,845 triệu đồng (hoặc có $1%$ rủi ro danh mục lỗ vượt quá con số này). Đây là cơ sở cốt lõi để thiết lập hạn mức cắt lỗ và tỷ lệ ký quỹ an toàn.

Kết luận

Chuyên đề tốt nghiệp "Mô hình Jump-Diffusion trong định giá quyền chọn" của tác giả Nguyễn Thị Ngọc Anh đã hoàn thành xuất sắc mục tiêu học thuật và ứng dụng thực tiễn:

  • Tính chuẩn xác về mặt toán học: Chứng minh thực nghiệm sự tồn tại của bước nhảy giá trên nhóm cổ phiếu xây dựng - bất động sản niêm yết (HBC, DIG, HDG) thông qua hệ thống kiểm định Barndorff-Nielsen & Shephard và Dickey-Fuller.
  • Tính hoàn thiện về giải thuật: Xây dựng trọn vẹn bộ công cụ tính toán Monte Carlo LSM cho quyền chọn kiểu Mỹ và quyền chọn Châu Á với độ chính xác và tốc độ hội tụ cao.
  • Giá trị thực tiễn vượt trội: Cung cấp bức tranh so sánh định lượng chi tiết giữa các loại hình quyền chọn và lượng hóa rủi ro danh mục $\text{VaR}_{99%}$, tạo tiền đề kỹ thuật vững chắc cho sự hình thành của thị trường giao dịch quyền chọn OTC và niêm yết tại Việt Nam.