disp (步骤1:加载和可视化的市场数据……”);

histData =负载(“HistoricalYTMData.mat”);%为每个系列的到期时间tsYTMMats = histData.tsYTMMats;%日期利率观察tsYTMObsDates = histData.tsYTMObsDates;%观察利率tsYTMRates = histData.tsYTMRates;%可视化收益率曲线如矿坑的= min (tsYTMRates (:));maxy = max (tsYTMRates (:));图;h =情节(tsYTMMats tsYTMRates,' k o ');轴([0,32岁,如矿坑的,maxy]);包含(“到期时间”);ylabel (“收益”);传奇(“历史收益率曲线”,“位置”,“本身”);网格在;为我= 1:50:长度(tsYTMObsDates)组(h,“ydata”tsYTMRates(我:));标题(datestr (tsYTMObsDates (i)));暂停(0.1);结束

disp (步骤2:计算市场不变量……”);

%不变量则认为是日常到期利率利率变化tsYTMRateDeltas = diff (tsYTMRates);%约1/3的30年期利率(列9)从最初的失踪%数据集。而不是扔掉所有这些观察,一个期望%最大化例程用于估计的均值和协方差%不变量。默认选项(NaN跳过初始估计,等等)。[tsInvMu, tsInvCov] = ecmnmle (tsYTMRateDeltas);%计算标准差和相关系数[tsInvStd, tsInvCorr] = cov2corr (tsInvCov);%的投资期限是1个月。(21工作日9/1/2010之间%和10/1/2010)。自从不变量是可和和手段%方差的正态分布是正常的,我们可以计划%不变量对投资期限如下hrznInvMu = 21 * tsInvMu ';hrznInvCov = 21 * tsInvCov;[hrznInvStd, hrznInvCor] = cov2corr (hrznInvCov);%显示结果disp (市场不变量将地平线以下统计数据的);

disp (”的意思是:“);

disp (hrznInvMu);

disp (的标准差:);

disp (hrznInvStd);

disp (的关系:);

disp (hrznInvCor);

disp (步骤3:模拟市场不变量在地平线……”);

%进行PCA分解利用不变量的协方差[pcaFactorCoeff, pcaFactorVar pcaFactorExp] = pcacov (hrznInvCov);%将pca decompositon的所有组件numFactors = 9;%创建PCA因素协方差矩阵pcaFactorCov = corr2cov (√pcaFactorVar)、眼睛(numFactors));%的数量模拟(随机画)numSim = 10000;%解决随机种子可重复的结果流= RandStream (“mrg32k3a”);RandStream.setGlobalStream(流);%从多元正态分布与零均值随机了%和对角协方差(numFactors pcaFactorSims = mvnrnd (0, 1), pcaFactorCov, numSim);%地平线变换不变量和计算统计数据hrznInvSims = pcaFactorSims * pcaFactorCoeff ' + repmat (hrznInvMu numSim 1);hrznInvSimsMu =意味着(hrznInvSims);hrznInvSimsCov = x (hrznInvSims);[hrznInvSimsStd, hrznInvSimsCor] = cov2corr (hrznInvSimsCov);%显示结果disp (模拟不变量有非常相似的统计原始不变量的);

disp (”的意思是:“);

disp (hrznInvSimsMu);

disp (的标准差:);

disp (hrznInvSimsStd);

disp (的关系:);

disp (hrznInvSimsCor);

disp (第四步:计算分布的安全返回地平线……”);

%购买和投资期限日期settleDate =“9/1/2010”;hrznDate =“10/1/2010”;%购买新发行国债的到期日期结算日期treasuryMaturities = {“9/1/2020”,“9/1/2030”,“9/1/2040”};%的观察收益率证券结算日期上的利益treasuryYTMAtSettle = tsYTMRates(结束,七章);%初始化数组以备后用treasuryYTMAtHorizonSim = 0 (numSim, 3);treasuryPricesAtSettle = 0 (1、3);treasuryPricesAtHorizonSim = 0 (numSim, 3);%使用实际/实际年化收益率计算基础基础= 8;%价格国债结算日期使用已知的到期收益率%注意:为简单起见,我们假定这些证券%包括利息。希望是,虽然价格可能不是%准确的总体结构/价值之间的关系%保存资产配置过程。为j = 1:3 treasuryPricesAtSettle (j) = bndprice (treasuryYTMAtSettle (j), 0, settleDate,…treasuryMaturities (j),“基础”、基础);结束%的美国国债价格,我们需要知道到期收益率% 9岁11个月,19年11个月,29年11个月%的模拟。我们使用三次样条插值近似这些%变换模拟不变量在地平线到期利率hrznYTMRatesSims = repmat (tsYTMRates(最终,:),numSim, 1) + hrznInvSims;hrznYTMMaturities = {“4/1/2011”,“10/1/2011”,“10/1/2012”,“10/1/2013”,…“10/1/2015”,“10/1/2017”,“10/1/2020”,“10/1/2030”,…“10/1/2040”};%将日期转换成数字串行日期x = datenum (hrznYTMMaturities);ξ= datenum (treasuryMaturities);%的数值精度,改变x值从0开始minDate = min (x);x = x - minDate;ξ= xi - minDate;%为每个模拟和成熟度近似收益率接近10,20年,30年%的节点。注意,花键配合的影响与线性适合有一个%显著影响产生的理想配置。这是由于%显著低估收益当使用一个线性的健康%点的已知节点为i = 1: numSim treasuryYTMAtHorizonSim(我,:)= interp1 (x, hrznYTMRatesSims(我,:),,样条的);结束% 1可视化模拟收益率曲线插值图;情节(x, hrznYTMRatesSims (1:)' k o '习,treasuryYTMAtHorizonSim (1:)“罗”);包含(的时间(天));ylabel (“收益”);传奇({“模拟收益率曲线”,“窜改收益率”},“位置”,“本身”);网格在;标题(“看样条与线性Interpolants变焦”);

%公债价格为每个模拟到期收益率在地平线%就是做了相同的假设在调用bndprice之上基础= 8 * 1 (numSim 1);为j = 1:3 treasuryPricesAtHorizonSim (:, j) = bndprice (treasuryYTMAtHorizonSim (:, j), 0,…hrznDate treasuryMaturities (j),“基础”、基础);结束%计算线性分布的回报treasuryReturns = (treasuryPricesAtHorizonSim - repmat (treasuryPricesAtSettle numSim 1))。/ repmat (treasuryPricesAtSettle numSim 1);%计算返回数据retsMean =意味着(treasuryReturns);retsCov = x (treasuryReturns);[retsStd, retsCor] = cov2corr (retsCov);%可视化结果30年期国债图;嘘(treasuryReturns, 100);标题(的收益分配10年,20年,30年国债的);网格在;传奇({“十年”,“20年”,“30年”});

disp (第五步:优化资产配置……”);

%计算10分的投影有效边界% std /回报空间。[portStd, portRet portWts] = portopt (retsMean retsCov 10);%可视化图;次要情节(2,1,1)情节(portStd portRet,' k o ');包含(“组合性病”);ylabel (“投资回报”);标题(“有效边界投影”);传奇(“最佳组合”,“位置”,“本身”);网格在;次要情节(2,1,2)栏(portWts,“堆叠”);包含(组合数的);ylabel (的组合权重);标题(的比例在每个财政投资);传奇({“十年”,“20年”,“30年”});

%夏普比率计算使用0无风险利率,因为我们是投资%的美国国债夏普= portRet. / portStd;基于近%的投资者选择投资组合意味着夏普比率sharpeTarget =意味着(夏普);investorChoice =找到(min (abs (sharpe-sharpeTarget)) = = abs (sharpe-sharpeTarget));investorPortfolioWts = portWts (investorChoice:);disp (的投资者比例分配在10年,20年,30年国债的);

disp (investorPortfolioWts);