volatility-modeling

Compare original and translation side by side

🇺🇸

Original

English

🇨🇳

Translation

Chinese

Volatility Modeling

波动率建模

Purpose

用途

Model, forecast, and interpret volatility using time-series models and options-implied measures. This skill covers EWMA and GARCH(1,1) for volatility forecasting, implied volatility extraction, volatility smile/skew/surface construction, the volatility term structure, the realized-vs-implied volatility gap (volatility risk premium), and the VIX index. These tools are foundational for options pricing, risk management, and trading strategy development.

使用时间序列模型和期权隐含指标对波动率进行建模、预测与解读。该技能涵盖用于波动率预测的EWMA和GARCH(1,1)、隐含波动率提取、波动率微笑/偏斜/曲面构建、波动率期限结构、已实现与隐含波动率差值（波动率风险溢价）以及VIX指数。这些工具是期权定价、风险管理和交易策略开发的基础。

Layer

层级

1b — Forward-Looking Risk

1b — 前瞻性风险

Direction

方向

Prospective

前瞻性

When to Use

适用场景

Forecasting future volatility for risk management or position sizing
Building EWMA or GARCH models to capture volatility clustering and mean reversion
Extracting implied volatility from option prices using Black-Scholes or other models
Constructing or interpreting volatility smiles, skews, and surfaces
Analyzing the volatility term structure across different maturities
Comparing realized (historical) volatility to implied volatility to assess the volatility risk premium
Understanding VIX and its relationship to market sentiment and expected risk

为风险管理或头寸规模预测未来波动率
构建EWMA或GARCH模型以捕捉波动率聚集和均值回归特性
使用Black-Scholes或其他模型从期权价格中提取隐含波动率
构建或解读波动率微笑、偏斜与曲面
分析不同到期期限的波动率期限结构
对比已实现（历史）波动率与隐含波动率以评估波动率风险溢价
理解VIX及其与市场情绪和预期风险的关系

Core Concepts

核心概念

EWMA (Exponentially Weighted Moving Average)

EWMA（指数加权移动平均）

A simple volatility model that gives more weight to recent observations. RiskMetrics popularized this approach with a standard decay factor.

sigma^2_t = lambda * sigma^2_{t-1} + (1 - lambda) * r^2_{t-1}

where:

lambda = decay factor (RiskMetrics standard: 0.94 for daily, 0.97 for monthly)
r_{t-1} = return in period t-1 (typically demeaned, but for daily returns the mean is often assumed to be zero)
sigma^2_{t-1} = previous period's variance estimate

Properties:

Assigns exponentially decaying weights to past squared returns.
Effective window is approximately 1/(1 - lambda) observations. For lambda = 0.94, effective window is approximately 17 days.
No mean reversion: the model is equivalent to IGARCH (integrated GARCH) where alpha + beta = 1. Volatility shocks persist indefinitely.
Simple to implement and requires only one parameter.

一种简单的波动率模型，赋予近期观测值更高权重。RiskMetrics推广了这种方法，并设定了标准衰减因子。

sigma^2_t = lambda * sigma^2_{t-1} + (1 - lambda) * r^2_{t-1}

其中：

lambda = 衰减因子（RiskMetrics标准：日度数据为0.94，月度数据为0.97）
r_{t-1} = t-1期的收益率（通常去均值，但日度收益率的均值常假设为0）
sigma^2_{t-1} = 上一期的方差估计值

特性：

对过去的平方收益率赋予指数衰减权重。
有效窗口约为1/(1 - lambda)个观测值。当lambda=0.94时，有效窗口约为17天。
无均值回归：该模型等价于IGARCH（集成GARCH），其中alpha + beta = 1。波动率冲击会无限持续。
实现简单，仅需一个参数。

GARCH(1,1)

The Generalized Autoregressive Conditional Heteroskedasticity model adds a constant term that induces mean reversion in volatility.

sigma^2_t = omega + alpha * r^2_{t-1} + beta * sigma^2_{t-1}

where:

omega > 0: constant term (determines long-run variance level)
alpha >= 0: reaction coefficient (sensitivity to recent shocks)
beta >= 0: persistence coefficient (memory of past variance)

Stationarity condition: alpha + beta < 1. This ensures the process is covariance-stationary and mean-reverting.

Long-run (unconditional) variance:

V_L = omega / (1 - alpha - beta)

Long-run annualized volatility: sigma_L = sqrt(V_L * 252).

Persistence: The quantity alpha + beta measures how quickly volatility reverts to its long-run level. Higher persistence means slower mean reversion.

Half-life of volatility shocks: The number of periods for a volatility shock to decay by half:

h = -ln(2) / ln(alpha + beta)

Since alpha + beta < 1, ln(alpha + beta) < 0, and h is positive.

Multi-step forecasts: The h-step-ahead GARCH(1,1) forecast:

E[sigma^2_{t+h}] = V_L + (alpha + beta)^h * (sigma^2_t - V_L)

The forecast converges to V_L as h approaches infinity.

广义自回归条件异方差模型添加了一个常数项，使波动率具备均值回归特性。

sigma^2_t = omega + alpha * r^2_{t-1} + beta * sigma^2_{t-1}

其中：

omega > 0：常数项（决定长期方差水平）
alpha >= 0：反应系数（对近期冲击的敏感度）
beta >= 0：持续系数（对过去方差的记忆性）

平稳性条件： alpha + beta < 1。这确保过程是协方差平稳且均值回归的。

长期（无条件）方差：

V_L = omega / (1 - alpha - beta)

长期年化波动率：sigma_L = sqrt(V_L * 252)。

持续性： alpha + beta的值衡量波动率回归至长期水平的速度。持续性越高，均值回归越慢。

波动率冲击的半衰期： 波动率冲击衰减一半所需的周期数：

h = -ln(2) / ln(alpha + beta)

由于alpha + beta < 1，ln(alpha + beta) < 0，因此h为正数。

多步预测： GARCH(1,1)的h步前瞻预测：

E[sigma^2_{t+h}] = V_L + (alpha + beta)^h * (sigma^2_t - V_L)

当h趋近于无穷大时，预测值收敛至V_L。

Implied Volatility

隐含波动率

The volatility value that, when plugged into an option pricing model (typically Black-Scholes), produces a theoretical price equal to the observed market price.

For a European call under Black-Scholes:

C = S * N(d1) - K * exp(-rT) * N(d2)

d1 = [ln(S/K) + (r + sigma^2/2) * T] / (sigma * sqrt(T))
d2 = d1 - sigma * sqrt(T)

Implied volatility is the sigma that solves C_model(sigma) = C_market. There is no closed-form solution; it must be found numerically (e.g., Newton-Raphson, bisection).

将其代入期权定价模型（通常为Black-Scholes）时，能使理论价格与观测到的市场价格相等的波动率值。

对于Black-Scholes框架下的欧式看涨期权：

C = S * N(d1) - K * exp(-rT) * N(d2)

d1 = [ln(S/K) + (r + sigma^2/2) * T] / (sigma * sqrt(T))
d2 = d1 - sigma * sqrt(T)

隐含波动率是满足C_model(sigma) = C_market的sigma值。没有闭式解，必须通过数值方法求解（如牛顿-拉夫逊法、二分法）。

Volatility Smile and Skew

波动率微笑与偏斜

In practice, implied volatility varies by strike price, contradicting the constant-volatility assumption of Black-Scholes.

Volatility smile: IV is higher for both deep in-the-money and deep out-of-the-money options, forming a U-shape. Common in FX markets.
Volatility skew (smirk): IV increases for lower strikes (OTM puts have higher IV than OTM calls). This is the dominant pattern in equity markets and reflects demand for downside protection and the reality of fat left tails.
Skew is often quantified as the difference in IV between a 25-delta put and a 25-delta call, or between 90% moneyness and 110% moneyness strikes.

实际中，隐含波动率随行权价格变化，这与Black-Scholes模型的常数波动率假设相矛盾。

波动率微笑： 深度实值和深度虚值期权的隐含波动率更高，形成U型。常见于外汇市场。
波动率偏斜（ smirk）： 行权价格越低，隐含波动率越高（虚值看跌期权的隐含波动率高于虚值看涨期权）。这是股票市场的主导模式，反映了对下行保护的需求以及肥左尾的现实。
偏斜通常量化为25-delta看跌期权与25-delta看涨期权的隐含波动率差值，或行权价为标的价格90%与110%时的隐含波动率差值。

Volatility Term Structure

波动率期限结构

Implied volatility varies across option expiration dates.

Normal (upward-sloping): Longer-dated options have higher IV. Reflects uncertainty increasing over time.
Inverted (downward-sloping): Near-term IV exceeds long-term IV. Common during market stress when short-term uncertainty spikes (e.g., around earnings, elections, crises).
Humped: IV peaks at an intermediate maturity. May occur around a specific anticipated event.

隐含波动率随期权到期日变化。

正常型（向上倾斜）： 长期期权的隐含波动率更高。反映不确定性随时间增加。
反转型（向下倾斜）： 短期隐含波动率高于长期隐含波动率。常见于市场压力时期，此时短期不确定性飙升（如财报季、选举、危机期间）。
驼峰型： 隐含波动率在中期到期日达到峰值。可能出现在特定预期事件前后。

Volatility Surface

波动率曲面

The two-dimensional surface of implied volatility across both strike (or delta/moneyness) and maturity. The volatility surface is the most complete representation of the options market's view of future uncertainty.

Practitioners interpolate the surface to price options at arbitrary strike/maturity combinations. Surface dynamics (how the surface shifts, tilts, and bends) are critical for options portfolio risk management.

隐含波动率在行权价（或delta/价内程度）和到期日两个维度上构成的二维曲面。波动率曲面是期权市场对未来不确定性看法的最完整体现。

从业者会对曲面进行插值，以对任意行权价/到期日组合的期权定价。曲面动态（曲面如何移动、倾斜和弯曲）对期权投资组合风险管理至关重要。

Realized vs Implied Volatility: The Volatility Risk Premium

已实现波动率 vs 隐含波动率：波动率风险溢价

Implied volatility systematically exceeds subsequent realized volatility on average. This gap is the volatility risk premium (VRP).

VRP = IV - RV_subsequent

The VRP exists because investors are willing to pay a premium for options (insurance), and option sellers demand compensation for bearing tail risk. The VRP is typically positive and has been a persistent source of return for volatility sellers.

Key considerations:

The VRP varies over time and is larger during periods of market stress.
Selling volatility (harvesting VRP) earns a steady premium but is exposed to large, infrequent losses.
The VRP can turn negative during extreme events.

平均而言，隐含波动率系统性地高于后续的已实现波动率。这个差值就是波动率风险溢价（VRP）。

VRP = IV - RV_subsequent

波动率风险溢价的存在是因为投资者愿意为期权（保险）支付溢价，而期权卖方要求承担尾部风险的补偿。波动率风险溢价通常为正，一直是波动率卖方的稳定收益来源。

关键考量：

波动率风险溢价随时间变化，在市场压力时期更大。
卖出波动率（获取波动率风险溢价）能获得稳定溢价，但面临大额、偶发的损失风险。
在极端事件期间，波动率风险溢价可能转为负数。

VIX Index

VIX指数

The CBOE Volatility Index measures the market's expectation of 30-day forward volatility, derived from S&P 500 option prices.

VIX is quoted in annualized percentage points (e.g., VIX = 20 means approximately 20% expected annualized vol).
VIX is computed from a wide strip of OTM put and call options, not from the Black-Scholes model.
VIX levels: 12-15 = low/complacent, 15-20 = normal, 20-30 = elevated, 30+ = high stress, 40+ = crisis.
VIX has strong mean-reverting properties and tends to spike during market selloffs ("fear gauge").

CBOE波动率指数衡量市场对未来30天波动率的预期，由标普500期权价格推导而来。

VIX以年化百分点报价（如VIX=20意味着预期年化波动率约为20%）。
VIX由一系列虚值看跌和看涨期权计算得出，而非基于Black-Scholes模型。
VIX水平：12-15 = 低/市场自满，15-20 = 正常，20-30 = 升高，30+ = 高压力，40+ = 危机。
VIX具有很强的均值回归特性，在市场抛售期间往往飙升（“恐慌指标”）。

Key Formulas

关键公式

Formula	Expression	Use Case
EWMA Variance	sigma^2_t = lambda * sigma^2_{t-1} + (1-lambda) * r^2_{t-1}	Simple volatility forecast
GARCH(1,1) Variance	sigma^2_t = omega + alpha * r^2_{t-1} + beta * sigma^2_{t-1}	Mean-reverting vol forecast
GARCH Long-Run Variance	V_L = omega / (1 - alpha - beta)	Unconditional variance level
GARCH Half-Life	h = -ln(2) / ln(alpha + beta)	Speed of mean reversion
GARCH h-Step Forecast	V_L + (alpha+beta)^h * (sigma^2_t - V_L)	Multi-period vol forecast
Black-Scholes Call	S * N(d1) - K * exp(-rT) * N(d2)	Option pricing (IV extraction)
Volatility Risk Premium	IV - RV_subsequent	Premium earned by vol sellers
EWMA Effective Window	approximately 1 / (1 - lambda)	Implicit lookback period

公式	表达式	使用场景
EWMA方差	sigma^2_t = lambda * sigma^2_{t-1} + (1-lambda) * r^2_{t-1}	简单波动率预测
GARCH(1,1)方差	sigma^2_t = omega + alpha * r^2_{t-1} + beta * sigma^2_{t-1}	均值回归波动率预测
GARCH长期方差	V_L = omega / (1 - alpha - beta)	无条件方差水平
GARCH半衰期	h = -ln(2) / ln(alpha + beta)	均值回归速度
GARCH h步预测	V_L + (alpha+beta)^h * (sigma^2_t - V_L)	多周期波动率预测
Black-Scholes看涨期权	S * N(d1) - K * exp(-rT) * N(d2)	期权定价（隐含波动率提取）
波动率风险溢价	IV - RV_subsequent	波动率卖方获取的溢价
EWMA有效窗口	approximately 1 / (1 - lambda)	隐含回溯周期

Worked Examples

实例演示

Example 1: EWMA Variance Update

实例1：EWMA方差更新

Given: Yesterday's variance estimate sigma^2_{t-1} = 0.0004 (daily vol = 2%), yesterday's return r_{t-1} = -3% (i.e., r = -0.03), and lambda = 0.94.

Calculate: Today's EWMA variance estimate and daily volatility.

Solution:

sigma^2_t = 0.94 * 0.0004 + (1 - 0.94) * (-0.03)^2
          = 0.94 * 0.0004 + 0.06 * 0.0009
          = 0.000376 + 0.000054
          = 0.000430

Daily volatility:

sigma_t = sqrt(0.000430) = 0.02074 = 2.074%

The large negative return (-3%) caused the volatility estimate to increase from 2.0% to 2.074%. The EWMA responded to the shock, but the high lambda (0.94) dampened the reaction.

已知： 昨日方差估计值sigma^2_{t-1}=0.0004（日度波动率=2%），昨日收益率r_{t-1}=-3%（即r=-0.03），lambda=0.94。

计算： 今日EWMA方差估计值和日度波动率。

解答：

sigma^2_t = 0.94 * 0.0004 + (1 - 0.94) * (-0.03)^2
          = 0.94 * 0.0004 + 0.06 * 0.0009
          = 0.000376 + 0.000054
          = 0.000430

日度波动率：

sigma_t = sqrt(0.000430) = 0.02074 = 2.074%

大幅负收益率（-3%）使波动率估计值从2.0%上升至2.074%。EWMA对冲击做出了反应，但高lambda值（0.94）削弱了反应幅度。

Example 2: GARCH(1,1) Long-Run Volatility and Half-Life

实例2：GARCH(1,1)长期波动率与半衰期

Given: GARCH(1,1) parameters estimated from daily S&P 500 returns: omega = 0.000002, alpha = 0.08, beta = 0.91.

Calculate: Long-run daily variance, long-run annualized volatility, and half-life of volatility shocks.

Solution:

Stationarity check: alpha + beta = 0.08 + 0.91 = 0.99 < 1 (stationary, but highly persistent).

Long-run variance:

V_L = 0.000002 / (1 - 0.99) = 0.000002 / 0.01 = 0.0002

Long-run daily volatility:

sigma_L = sqrt(0.0002) = 0.01414 = 1.414%

Annualized:

sigma_annual = 0.01414 * sqrt(252) = 22.45%

Half-life:

h = -ln(2) / ln(0.99) = -0.6931 / (-0.01005) = 68.97 ~ 69 trading days

Interpretation: After a volatility shock, it takes approximately 69 trading days (about 3 months) for the excess volatility to decay by half. This high persistence (alpha + beta = 0.99) is typical for equity index returns.

已知： 由日度标普500收益率估计的GARCH(1,1)参数：omega=0.000002，alpha=0.08，beta=0.91。

计算： 长期日度方差、长期年化波动率以及波动率冲击的半衰期。

解答：

平稳性检验： alpha + beta = 0.08 + 0.91 = 0.99 < 1（平稳，但持续性极高）。

长期方差：

V_L = 0.000002 / (1 - 0.99) = 0.000002 / 0.01 = 0.0002

长期日度波动率：

sigma_L = sqrt(0.0002) = 0.01414 = 1.414%

年化：

sigma_annual = 0.01414 * sqrt(252) = 22.45%

半衰期：

h = -ln(2) / ln(0.99) = -0.6931 / (-0.01005) = 68.97 ~ 69个交易日

解读：波动率冲击后，大约需要69个交易日（约3个月）才能使超额波动率衰减一半。这种高持续性（alpha + beta=0.99）在股票指数收益率中很典型。

Example 3: Implied Volatility Interpretation

实例3：隐含波动率解读

Given: A stock trades at $100. A 3-month ATM call (K = $100) trades at $6.50. The risk-free rate is 5%. Using Black-Scholes, the implied volatility is determined (via numerical solver) to be 30%.

Calculate: What does this tell us, and how does it compare to realized vol of 22%?

Solution:

The implied volatility of 30% represents the market's consensus forecast of annualized volatility over the next 3 months, as embedded in option prices.

Comparing to realized (historical) volatility of 22%:

VRP = IV - RV = 30% - 22% = 8%

The positive 8-percentage-point gap is the volatility risk premium. Possible interpretations:

The market expects volatility to rise above recent realized levels.
Option sellers are demanding a premium for bearing tail risk.
There may be an upcoming event (earnings, regulatory decision) that could cause a volatility spike.

A systematic vol-selling strategy would sell this option, expecting to profit from the VRP if realized vol remains near 22%. However, the seller bears the risk that realized vol could exceed 30%.

已知： 某股票价格为100美元。3个月平值看涨期权（行权价K=100美元）成交价为6.50美元。无风险利率为5%。使用Black-Scholes模型，通过数值求解器确定隐含波动率为30%。

分析： 这说明了什么？与22%的已实现波动率相比如何？

解答：

30%的隐含波动率代表了期权价格中所包含的市场对未来3个月年化波动率的一致预期。

与22%的已实现（历史）波动率对比：

VRP = IV - RV = 30% - 22% = 8%

8个百分点的正差值即为波动率风险溢价。可能的解读：

市场预期波动率将高于近期已实现水平。
期权卖方要求承担尾部风险的溢价。
可能有即将到来的事件（财报、监管决策）导致波动率飙升。

系统性波动率卖出策略会卖出该期权，如果已实现波动率维持在22%左右，有望从波动率风险溢价中获利。但卖方需承担已实现波动率超过30%的风险。

Common Pitfalls

常见误区

GARCH stationarity: alpha + beta must be strictly less than 1 for the GARCH(1,1) process to be covariance-stationary. If alpha + beta >= 1, the long-run variance is undefined and the model is IGARCH (or explosive). Always check this condition after estimation.
EWMA has no mean reversion: EWMA is equivalent to IGARCH (alpha + beta = 1), so volatility shocks never decay. This makes EWMA unsuitable for long-horizon volatility forecasts where mean reversion is expected.
Implied volatility is model-dependent: IV extracted using Black-Scholes assumes log-normal returns, constant volatility, continuous trading, and no jumps. The "smile" and "skew" exist precisely because these assumptions are violated. IV is a quoting convention, not a true forecast.
Conflating historical vol with forward-looking vol: Historical (realized) volatility measures what has happened. Implied volatility reflects market expectations about the future. They address different questions and can diverge substantially.
Using daily GARCH for long-horizon forecasts without term structure adjustment: Daily GARCH forecasts converge to the unconditional variance at long horizons. Use the multi-step forecast formula E[sigma^2_{t+h}] = V_L + (alpha+beta)^h * (sigma^2_t - V_L) and aggregate appropriately.
Overfitting GARCH models: Higher-order GARCH(p,q) models or extensions (EGARCH, TGARCH, GJR-GARCH) can overfit in-sample. GARCH(1,1) is remarkably hard to beat out-of-sample for most financial return series. Start with GARCH(1,1) and justify added complexity.
Lambda selection for EWMA: The choice of lambda significantly affects responsiveness. lambda = 0.94 responds quickly to shocks (effective window approximately 17 days); lambda = 0.97 is smoother (effective window approximately 33 days). The choice should match the application's horizon.

GARCH平稳性： GARCH(1,1)过程要具备协方差平稳性，alpha + beta必须严格小于1。如果alpha + beta >=1，长期方差无定义，模型为IGARCH（或爆炸性）。估计后务必检查该条件。
EWMA无均值回归： EWMA等价于IGARCH（alpha + beta=1），因此波动率冲击永远不会衰减。这使得EWMA不适用于预期存在均值回归的长期波动率预测。
隐含波动率依赖模型： 使用Black-Scholes提取的隐含波动率假设收益率服从对数正态分布、波动率恒定、连续交易且无跳空。“微笑”和“偏斜”的存在正是因为这些假设不成立。隐含波动率是一种报价惯例，而非真实预测。
混淆历史波动率与前瞻性波动率： 历史（已实现）波动率衡量过去发生的情况。隐含波动率反映市场对未来的预期。它们解决的是不同问题，可能大幅偏离。
未调整期限结构就使用日度GARCH进行长期预测： 日度GARCH预测在长期会收敛至无条件方差。使用多步预测公式E[sigma^2_{t+h}] = V_L + (alpha+beta)^h * (sigma^2_t - V_L)并进行适当聚合。
过度拟合GARCH模型： 高阶GARCH(p,q)模型或扩展模型（EGARCH、TGARCH、GJR-GARCH）可能在样本内过度拟合。对于大多数金融收益率序列，GARCH(1,1)在样本外表现出色，很难被超越。从GARCH(1,1)开始，再证明增加复杂度的合理性。
EWMA的lambda选择： lambda的选择显著影响模型的响应速度。lambda=0.94对冲击响应迅速（有效窗口约17天）；lambda=0.97更平滑（有效窗口约33天）。选择应与应用场景的时间范围匹配。

Cross-References

交叉引用

historical-risk (wealth-management plugin, Layer 1a): Close-to-close, Parkinson, and Yang-Zhang volatility estimators provide the realized volatility benchmarks against which GARCH forecasts and implied volatility are compared.
forward-risk (wealth-management plugin, Layer 1b): Volatility forecasts from EWMA and GARCH are direct inputs to parametric and Monte Carlo VaR calculations.
performance-metrics (wealth-management plugin, Layer 1a): Volatility estimates affect the denominators of Sharpe, Sortino, and other risk-adjusted ratios. Using forward-looking (GARCH) volatility can produce conditional performance ratios.

historical-risk（财富管理插件，层级1a）：收盘价-收盘价、Parkinson和Yang-Zhang波动率估计器提供了已实现波动率基准，用于对比GARCH预测值和隐含波动率。
forward-risk（财富管理插件，层级1b）：EWMA和GARCH的波动率预测是参数法和蒙特卡洛VaR计算的直接输入。
performance-metrics（财富管理插件，层级1a）：波动率估计值影响Sharpe、Sortino等风险调整比率的分母。使用前瞻性（GARCH）波动率可生成条件绩效比率。

Reference Implementation

参考实现

See

scripts/volatility_modeling.py

for computational helpers.

计算辅助工具请参见

scripts/volatility_modeling.py

。