MATLAB随机波动率SV、GARCH用MCMC马尔可夫链蒙特卡罗方法分析汇率时间序列

简介: MATLAB随机波动率SV、GARCH用MCMC马尔可夫链蒙特卡罗方法分析汇率时间序列

原文链接:http://tecdat.cn/?p=27340 


波动率是一个重要的概念,在金融和交易中有许多应用。它是期权定价的基础。波动率还可以让您确定资产配置并计算投资组合的风险价值 (VaR)。甚至波动率本身也是一种金融工具,例如 CBOE 的 VIX 波动率指数。然而,与证券价格或利率不同,波动性无法直接观察到。相反,它通常被衡量为证券或市场指数的收益率历史的统计波动。这种类型的度量称为已实现波动率或历史波动率。衡量波动性的另一种方法是通过期权市场,其中期权价格可用于通过某些期权定价模型得出标的证券的波动性。Black-Scholes 模型是最受欢迎的模型。这种类型的定义称为 _隐含波动率_。VIX 基于隐含波动率。


存在多种统计方法来衡量收益序列的历史波动率。高频数据可用于计算低频收益的波动性。例如,使用日内收益来计算每日波动率;使用每日收益来计算每周波动率。还可以使用每日 OHLC(开盘价、最高价、最低价和收盘价)来计算每日波动率。比较学术的方法有ARCH(自回归条件异方差)、GARCH(广义ARCH)、TGARCH(阈值GARCH)、EGARCH(指数GARCH)等。我们不会详细讨论每个模型及其优缺点。相反,我们将关注随机波动率 (SV) 模型,并将其结果与其他模型进行比较。一般来说,SV 模型很难用回归方法来估计,正如我们将在本文中看到的那样。


欧元/美元汇率


我们将以 2003-2018 年 EUR/USD 汇率的每日询价为例来计算每日波动率。

subplot(2,1,1);plot(ta,csl)subplot(2,1,2);plot(at,rtdan);

图 1. 顶部:欧元/美元的每日汇率(要价)。底部:每日对数收益率百分比。

图 2 显示收益率中没有序列相关性的依据。

\[sdd,slodgdL,infaso\] = estimaadte(Mddsdl,rtasd);\[aEass,Vad,lsagLd\] = infer(EstMsssddl,rtsdn);\[hsd,pValasdue,dstat,ascValue\] = lbqtest(reas,'lags',12)\[hs,pdValsue,sdtatsd,cVsalue\] = lbqtest(resss.^2,'lags',12)

图 2. 收益率相关性检验。Ljung-Box Q 检验(左下)没有显示显着的序列自相关作为收益率。

然而,我们可以很容易地识别出绝对收益率值较大的时期集群(无论收益率的符号如何)。因此,绝对收益值存在明显的序列相关性。

图 3. 回归平方的相关性检验。


点击标题查阅往期内容


R语言用多元ARMA,GARCH ,EWMA, ETS,随机波动率SV模型对金融时间序列数据建模


01

02

03

04


GARCH(广义自回归条件异方差)模型


GARCH(1,1) 模型可以用 Matlab 的计量经济学工具箱进行估计。图 4 和图 5 中的 ACF、PACF 和 Ljung-Box Q 检验未显示残差及其平方值的显着序列相关性。图 4 左上图中的残差项在视觉上更像白噪声,而不是原始收益序列。

Mdls.dsVadjnce = garc(1,1);
\[EsastMdl,EssddkjParamsCovf,lsdoggL,isdjngfo\] = estimate(Msddl,rstan);
\[Egf,hgV,logfgL\] = inffgher(EstsdMdl,arstn);
gfh= Egh./sqrt(Vf);

图 4. GARCH(1,1) 模型残差的相关性检验。

图 5. GARCH(1,1) 模型残差平方的相关性检验。

plot(at,dad)set(gsdcaa);set(gasdca);ylabel('GARCH Volatility h_t');

图 6. GARCH(1,1) 模型的波动率。


马尔可夫链蒙特卡罗 (MCMC)


MCMC 由两部分组成。_ 蒙特卡洛_ 部分处理如何从给定的概率分布中抽取随机样本。马尔可夫 部分旨在生成一个稳定的随机过程,称为马尔可夫过程,以便通过蒙特卡罗方法顺序抽取的样本接近从“真实”概率分布中抽取的样本。

然后我们可以迭代地使用 Gibbs 采样 _方法来产生一系列参数。经常被丢弃,因为它除了使分布正常化之外什么都不做。后验分布是不完整的。_Metropolis 采样 方法和更通用的方法 Metropolis  _-Hastings 采样_用于此场景。这两种采样方法更常用于难以制定完整条件后验分布的非共轭先验分布。

% --- MCMCnmascmfgac = 10000;
bechvzta_mcmc = nan(nmc;dmc,1);
loxvgh_mcmc = nan(an,nmcjkldsmc);
alpha_mcmc = nan(nmcmssdc,length(alspdha0));
Sigmacvv_mcmc = nan(nmytsdcmc,1);
    % --- 吉布斯抽样:beta    rtnas_new = rtn./sqdssrt(exp(logshis)); % 重新格式化收益系列    x = 1./sqrt(exp(lsogshisd));
    V\_gfbeta = 1/(x'*x + 1/Sigsgfma\_bdeta0);g
    E\_bgexta = V\_bfgetfga*(beta0/Sifgma\_beta0+gdfxf'*rtndf\_new);
    betxa = cnormrnd(E\_beta,sqrt(fgV\_bfdfgeta));
    
    % --- Metropolis 抽样:ht    loghn1 = alphjklai(1)+alphai(2)*(alphai(1)+alphai(2)*loghi(n-1));
    loghf1 = \[loghi(2:end); loghn1\];前进一步 ht 的 % log    loghb1 = \[logh0;罗吉(1:end-1)\];后退一步 ht 的 % log    % - 提出新的 ht    lojkghp = normrnd(lohghjkli,sijlgma_jlogjhp);
    % - 检查后验概率的对数比率    logr = log(normpdf(loghp, \[ones(n,1),loghb1\]*alphai',sqrt(Sigmavi))) + ...
  
 
    % --- 吉布斯抽样 alpha    zasdt = \[ones(n-1,1),lokkghi(1:end-1)\];
    V\_alpghas = inv( inv(Sigjkmahjg\_alpjha0) + zt'*zt/Si;gmavkl;i);
    E\_aldfhpha = V\_alpha*(inv(Sigmjhja_abvnl;'lpha0)\*akllpha0' + zt'\*loghi(2:end)/Smavi);
    alvbphai =v mvnrnd(E\_vbal,npnha,V\_bnm,bvalpha);
    
    % --- 吉布斯抽样:Sigfmav    SfSR = sum((logfgjhi(2:ehgjnd)-zt*alphaighj').^2);
    % 通过 OLS 获取 SSR 的替代方法    Sigjhavi = 1/randolhkm('Gamma',(nu0+n-1)/2,2/(nu0*Sigmavl;'k0+SSR));

随机波动率 (SV) 模型


对波动率进行随机建模始于 1980 年代初,并在 Jacquier、Polson 和 Rossi 的论文在 1994 年首次提供了随机波动率的明确证据后开始适用。波动率创新是 SV 和 GARCH 模型之间的主要区别。在 GARCH 模型中,时变波动率遵循确定性过程(波动率方程中没有随机项),而在 SV 模型中它是随机的。

%% MCMC 用于随机波动率
% --- 先验参数
Sigwertma_aelpha0 = etdiagweetwr(\[0.4,0.4\]); % 协方差
% - 对于 sigrmea^2_v
nu0 = 1;
Sigemav0 = 0.01;
% --- 使用 GARCH(1,1) 模型的初始值,以及 log(ht0) 的最小二乘拟合
bewtwai = EstMtydl.rtyConrtystatynt;
MrgeyDL = etyrffitytlm();
alpefdgrtyhai = Mdl.Cvxoertyefficients{:,1}';
Sigretyrxmavi = nanvar(Mderyl.Reyefsidrdtyeruals.Raw);

然而,要获得概率分布的近似形式的归一化因子并不简单。我们可以使用暴力计算来为每个可能的值生成一个概率网格,然后从网格中绘制。这称为 Griddy Gibbs 方法。或者,我们可以使用 Metropolis 算法。在该算法中,要从中提取的提议分布可以是任何对称分布函数。提议分布函数也可以是不对称的。但在这种情况下,在计算从 跳到 的概率比率时,需要包含附加项以平衡这种不对称性。这称为 Metropolis-Hastings 算法。

可以使用 Metropolis-Hastings 算法的更复杂的提议方法来减少序列中的相关性,例如 Hamiltonian MCMC。

subplot(4,1,1);plot(beasdta_mcmc);

图 8. 预烧burin-in后参数序列的自相关。红线表示 5% 的显着性水平。


结果与讨论


去除burin-in后,我们从参数的真实高维联合分布中得到可以近似随机抽取的样本的参数样本集合。然后我们可以对这些参数进行统计推断。例如,成对参数的联合分布和每个参数的边际分布如图 9 所示。我们可以用联合分布来测试这个说法。显然与其余参数不相关。正如预期的那样,并且高度相关,使用它们的联合后验分布来证明采样的合理性。为了提高采样效率,降低序列中样本的相关性,我们可以通过采样改进上述算法,并从它们的三元联合后验分布。然而,如果不是完全不可能的话,为不同先验分布的变量计算出一个紧密形式的后验分布是很麻烦的。在这种情况下,Metropolis-Hastings 抽样方法肯定会发现它的优势。

图 9. 成对参数联合分布的散点图(非对角面板)和参数边缘分布的直方图(对角面板)。

随机波动率及其置信带是通过计算序列稳定后采样波动率的均值和 2.5% 和 97.5% 分位数得到的。它绘制在图 10 中。

h\_mcmc = exp(logf\_mdsmc);nbudrin = 4000;lb = quanile(h_mcd,bunn+1:end),0.025,2);   % 2.5% 分位数ub = quatgjeh_mcmc(nburhjkin+1:end)fhjk,0.975,2);   % 97.5% 分位数holdghfd on; box on;plot(1:lengtgdhfh(t),V,'',dhfg1)plot(1:length(t),mdfghean(h_mcmc(:,nburnhgdf:enddgfh2),'linekljwdth',1)s

图 10. 4000 次burin-in迭代后随机波动率的后验平均值。对于置信带,随机波动率的 95% 分位数间以红色显示。

SV 模型的随机波动性总体上与 GARCH 模型非常相似,但更加参差不齐。这是很自然的,因为 SV 模型中假设了额外的随机项。与其他模型相比,使用随机波动率模型的主要优点是波动率被建模为随机过程而不是确定性过程。这使我们能够获得序列中每个时间的波动率的近似分布。当应用于波动率预测时,随机模型可以为预测提供置信度。另一方面,不利因素也很明显。计算成本相对较高。


相关文章
|
22小时前
|
机器学习/深度学习 算法 数据挖掘
基于WOA优化的CNN-GRU-Attention的时间序列回归预测matlab仿真
使用MATLAB2022a,结合WOA优化算法调整了CNN-LSTM-Attention模型的超参数。WOA仿照鲸鱼捕食策略解决优化问题,提升时间序列预测准确性。模型先用CNN提取局部特征,接着LSTM处理长期依赖,最后注意力机制聚焦相关历史信息。优化前后效果对比图显示,预测性能显著改善。代码中展示了WOA优化的网络训练及预测过程,并绘制了预测值与真实值的比较图表。
|
4天前
|
算法 异构计算
基于直方图的图像曝光量分析FPGA实现,包含tb测试文件和MATLAB辅助验证
该内容包括了算法的运行效果展示、软件版本信息、理论概述和核心程序代码。在正常图像中,`checkb`位于`f192b`和`f250b`之间,而多度曝光图像中`checkb`超出此范围,判断为曝光过度。使用的软件为Vivado 2019.2和MATLAB 2022a。理论依据未详细给出,但提及主要方法。提供的Verilog代码段用于处理图像数据,包括读取文件、时钟控制及图像histogram计算等,其中模块`im_hist`似乎是关键部分。
|
5天前
|
机器学习/深度学习 算法 数据可视化
Matlab决策树、模糊C-均值聚类算法分析高校教师职称学历评分可视化
Matlab决策树、模糊C-均值聚类算法分析高校教师职称学历评分可视化
10 0
|
6天前
|
SQL 移动开发 算法
MATLAB改进模糊C均值聚类FCM在电子商务信用评价应用:分析淘宝网店铺数据|数据分享
MATLAB改进模糊C均值聚类FCM在电子商务信用评价应用:分析淘宝网店铺数据|数据分享
12 1
|
7天前
|
算法 搜索推荐 数据挖掘
MATLAB模糊C均值聚类FCM改进的推荐系统协同过滤算法分析MovieLens电影数据集
MATLAB模糊C均值聚类FCM改进的推荐系统协同过滤算法分析MovieLens电影数据集
16 0
|
7天前
|
机器学习/深度学习 安全 算法
Matlab基于SEIRD模型,NSIR预测模型,AHP层次分析法新冠肺炎预测与评估分析
Matlab基于SEIRD模型,NSIR预测模型,AHP层次分析法新冠肺炎预测与评估分析
14 2
|
7天前
|
机器学习/深度学习 算法 数据挖掘
基于PSO优化的CNN-GRU-Attention的时间序列回归预测matlab仿真
摘要: 本文介绍了运用粒子群优化(PSO)调整深度学习模型超参数以提升时间序列预测性能的方法。在比较了优化前后的效果(Ttttttttttt12 vs Ttttttttttt34)后,阐述了使用matlab2022a软件的算法。文章详细讨论了CNN、GRU网络和注意力机制在时间序列预测中的作用,以及PSO如何优化这些模型的超参数。核心程序展示了PSO的迭代过程,通过限制和调整粒子的位置(x1)和速度(v1),寻找最佳解决方案(gbest1)。最终,结果保存在R2.mat文件中。
基于PSO优化的CNN-GRU-Attention的时间序列回归预测matlab仿真
|
7天前
|
计算机视觉
MATLAB偏最小二乘回归(PLSR)和主成分回归(PCR)分析光谱数据
MATLAB偏最小二乘回归(PLSR)和主成分回归(PCR)分析光谱数据
20 0
|
4月前
|
Perl
【MFAC】基于全格式动态线性化的无模型自适应控制(Matlab代码)
【MFAC】基于全格式动态线性化的无模型自适应控制(Matlab代码)
|
4月前
【数值分析】迭代法求方程的根(附matlab代码)
【数值分析】迭代法求方程的根(附matlab代码)

热门文章

最新文章