2026/4/18 12:00:38
网站建设
项目流程
wordpress网站生成app应用,手机百度 网站提交,信阳网站建设制作公司,淘宝上面如何做网站策略功能与风险说明
本策略通过Johansen协整检验筛选具有长期均衡关系的标的资产组合#xff0c;构建指数期权统计套利头寸。核心功能包括#xff1a;1) 多变量时间序列协整关系验证#xff1b;2) 动态跟踪误差控制#xff1b;3) 期权希腊字母对冲。主要风险包含模型误设风…策略功能与风险说明本策略通过Johansen协整检验筛选具有长期均衡关系的标的资产组合构建指数期权统计套利头寸。核心功能包括1) 多变量时间序列协整关系验证2) 动态跟踪误差控制3) 期权希腊字母对冲。主要风险包含模型误设风险、流动性风险及极端市场条件下的协整关系断裂风险。Johansen协整检验理论基础向量误差修正模型(VECM)框架在VAR§模型基础上引入协整关系约束建立VECM模型ΔXₜ ΠXₜ₋₁ Γ₁ΔXₜ₋₁ ... Γₚ₋₁ΔXₜ₋ₖ₊₁ εₜ其中协整矩阵Π分解为αβᵀβ矩阵列向量构成协整空间。特征值检验方法通过极大似然估计获取特征值λᵢ构造迹统计量Trace(H₁/H₀) -T∑(ln(1-λ̂ᵢ)) ~ χ²(r)最大特征值统计量Max(H₁/H₀) -Tln(1-λ̂ᵣ₊₁) ~ χ²(1)标的筛选与组合构建流程数据预处理阶段importnumpyasnpimportpandasaspdfromstatsmodels.tsa.vector_ar.vecmimportcoint_johansendefpreprocess_data(price_df,window252):计算对数收益率并处理缺失值returnsnp.log(price_df).diff().dropna()returnreturns.rolling(window).mean().dropna()协整检验实施defjohansen_test(returns,lag1,alpha0.05):执行Johansen协整检验resultcoint_johansen(returns,det_order0,k_ar_difflag)eigenvaluesresult.eigvals.round(4)trace_statresult.trace_stat.round(4)max_statresult.max_stat.round(4)# 临界值判断95%置信度cv_trace[15.87,19.60,23.42,27.04]cv_max[14.07,18.63,22.00,25.32]cointegration_rank0foriinrange(len(trace_stat)-1,-1,-1):iftrace_stat[i]cv_trace[i]:cointegration_rank1else:breakreturn{eigenvalues:eigenvalues,trace_stat:trace_stat,max_stat:max_stat,rank:cointegration_rank,critical_values:{trace:cv_trace,max:cv_max}}组合权重优化defcalculate_hedge_ratio(series_a,series_b):计算最小方差对冲比率cov_matrixnp.cov(series_a,series_b)deltacov_matrix[0,1]/cov_matrix[1,1]returndeltadefbuild_portfolio(candidates,test_result):根据协整结果构建投资组合iftest_result[rank]1:returnNoneselectedcandidates.loc[candidates.index.isin(test_result[eigenvectors][:test_result[rank]].index)]weights{}foriinrange(len(selected.columns)-1):weights[selected.columns[i]]calculate_hedge_ratio(selected.iloc[:,i],selected.iloc[:,i1])returnpd.Series(weights)期权头寸构建逻辑Greeks参数计算fromscipy.statsimportnormdefblack_scholes_greeks(S,K,T,r,sigma,option_typecall):BS模型希腊字母计算d1(np.log(S/K)(r0.5*sigma**2)*T)/(sigma*np.sqrt(T))d2d1-sigma*np.sqrt(T)ifoption_typecall:deltanorm.cdf(d1)gammanorm.pdf(d1)/(S*sigma*np.sqrt(T))theta-(S*norm.pdf(d1)*sigma)/(2*np.sqrt(T))-r*K*np.exp(-r*T)*norm.cdf(d2)else:deltanorm.cdf(d1)-1gammanorm.pdf(d1)/(S*sigma*np.sqrt(T))theta-(S*norm.pdf(d1)*sigma)/(2*np.sqrt(T))r*K*np.exp(-r*T)*norm.cdf(-d2)return{delta:delta,gamma:gamma,theta:theta}动态对冲策略defdynamic_hedging(portfolio,option_chain,rebalance_freq5):基于希腊字母的动态对冲holdings{}forassetinportfolio.index:# 获取最优行权价期权best_optionselect_best_option(asset,option_chain)# 计算对冲比例hedge_ratioportfolio[asset]*best_option[delta]# 确定头寸规模position_sizeint(abs(hedge_ratio)*100)holdings[asset]{option_contract:best_option[contract],position_size:position_size,hedge_ratio:hedge_ratio}returnholdings