clear,clc;
global NUM_PATH NUM_T VOLA S0
NUM_PATH=1000;%股票路径数
NUM_T=25;%时间点数目
DELT_T=3/25;%相邻时间点间隔时间(年)
S0=8.30;%初始股价
SZ=13.63;%转股价
V0LA=0.324;%股价年波动率
RF=0.0225;%无风险利率
START_T=0;%转换权允许行使...clear,clc;
global NUM_PATH NUM_T VOLA S0
NUM_PATH=1000;%股票路径数
NUM_T=25;%时间点数目
DELT_T=3/25;%相邻时间点间隔时间(年)
S0=8.30;%初始股价
SZ=13.63;%转股价
V0LA=0.324;%股价年波动率
RF=0.0225;%无风险利率
START_T=0;%转换权允许行使起点
SELL_CON_DAY=20;%回售(即卖回)条款中规定的允许回售的连续天数
SELL_P=103;%回售价格
BUY_CON_DAY=20;%赎回条款中规定的赎回的连续天数
BUY_P=105;%赎回价格
s=get_paths(NUM_PATH,NUM_T,S0,VOLA,RF,DELT_T);%生成股价路径
sell_back_mat=get_sell_back_mat(s,SELL_CON_DAY,SELL_P,valid_d);
buy_back_mat =get_sell_back_mat(s,BUY_CON_DAY,BUY_P,valid_d);
save data;
%==========================================================================
function s_star=get_paths(NUM_PATH,NUM_T,S0,VOLA,RF,DELT_T)
%函数功能:生成股价路径
%输入参数:NUM_PATH:生成的路径数目,NUM_T:股票价格期数,S0:初始股价
%VOLA:股价波动率,RF:无风险利率(年),DETL_T:每期时间间隔(年)
%返回:s_star,股价路径矩阵(矩阵大小NUM_PATH*NUM_T)
eps=normrnd(0,l,NUM_PATH,NUM_T);%产生符合标准正态分布的乱数
s(l:NUM_PATH,1)=S0;%第1列全部置为初始股价
for i=l:NUM_PATH
for k=2:NUM_T
s(i,k)=s(i,k-l)*exp((RF-0.5*VOLA^2)*DETL_T+VOLA*eps(i,k)*DELT_T^0.5);
end
end
% EMS for stock price paths
s_tar(:,l)=s(:,l);
for tj=2:NUM_T
ztj=s_star(:,tj-l).*(s(:,tj)./s(:,tj-l));
z0tj=mean(ztj)*exp(-RF*(tl+(tj-2)*t2));
s_star(:,tj)=s_star(:,l).*(ztj/z0tj);
end
function sell_back_mat=get_sell_back_mat(s,continue_dates,sellback_sp,valid_d)
%函数功能:按照回售条款:回售生效日后，若某日前连续若干天的股价都低于约定股价，则该日有权卖回
%本函数根据日股价路径矩阵，确定哪些日期满足回售条件
%输入参数s:股价路径矩阵
%continue_dates:约定的连续天数
%sellback_sp:约定的股价
%valid_d:允许回售生效日
%注意:本函数假定valid_d总是大于continue_dates+l
%返回:返回值se!l_back_mat为一大小与股价路径矩阵相同的矩阵，
%其行表示路径，列表示交易日，各矩阵表示是否满足回售条件，O-不满足 1-满足
valid_d=max(valid_d,continue_dates+l);
rows=size(s,l);%计算行数(路径数)
cols=size(s,2);%计算列数(交易日数)
sell_back_mat(l:rows,l:valid_d)=0;%生效日前一律为0
for j=valid_d:cols
sell_back_mat(:,j)=sum(s(:,j-continue_dates:j-l)=sellback_sp,2)==continue_dates;
end展开 |
|
转播
