量化投资策略与技术期末报告样本

《量化投资方略与技术》结课报告学科专业：研究方向：姓名：学号：7月9日

目录1目标股票的选取 31.1备选对象 31.2择股标准 42投资组合的构建 42.1目标股票的收益率 42.2最优投资组合 63预期收益和风险 74风险对冲策略 8附录A：原始数据 10表A1:财务指标 10表A2：投资组合股票收益率（节选） 11表A3：投资组合日收益率（节选） 11表A4：资产组合预期收益率（节选） 12表A5：沪深300股指期货日收盘价和收益率（节选） 12附录B：代码 13B1：择股 13B2：投资组合风险-收益的计算 13B3：最优投资组合求解 14

1目的股票选用1.1备选对象从上海主板、深圳主板、中小板和创业板分别随机选取10只股票，列示于表1.1中：表1.1投资组合备选对象所属板块股票名称股票代码上海主板武钢股份600005中华人民共和国国贸600007首创股份600008包钢股份600010华能国际600011日照港600017上港集团600018上海电力600021中华人民共和国石化600028三一重工600031深圳主板万科A000002.SZ世纪星源000005.SZ深振业A000006.SZ零七股份000007.SZ中华人民共和国宝安000009.SZ南玻A000012.SZ沙河股份000014.SZ深康佳A000016.SZ深华发A000020.SZ深赤湾A000022.SZ中小板新和成00.SZ鸿达兴业00.SZ伟星股份00.SZ华邦颖泰00.SZ德豪润达00.SZ精功科技00.SZ大族激光00.SZ天奇股份00.SZ传化股份00.SZ盾安环境00.SZ创业板特锐德300001南风股份300004莱美药业300006汉威电子300007立思辰300010鼎汉技术300011新宁物流300013亿纬锂能300014银江股份300020大禹节水3000211.2择股原则本文使用财务选股办法选取投资组合目的股票，所使用财务指标来源于上市公司财务报表，详细原则为：①每股收益（EPS）排名前16位；②净资产收益率（ROE）排名前16位；③销售净利率（NPM）排名前20位同步满足以上三个条件成为备选股票。使用R软件对原始数据进行解决（解决过程见附录B1），选取出满足条件10只股票，剔除停牌股票后，共有8只股票进入投资组合。（表1.2）表1.2投资组合股票股票代码股票名称所属市场备注600007中华人民共和国国贸上海主板600011华能国际600018上港集团600021上海电力000002.SZ万科A深圳主板000022.SZ深赤湾A停牌00.SZ大族激光中小板00.SZ传化股份00.SZ新和成00.SZ伟星股份停牌2投资组合构建2.1目的股票收益率8只股票日收益率和合计收益率如图2.1和2.2所示：图2.1日收益率图2.2合计收益率股票有关系数矩阵如下：万科A华能国际上海电力新和成大族激光传化股份中华人民共和国国贸上港集团万科A1.0000.4030.3360.2920.2850.3090.4840.346华能国际0.4031.0000.5100.2810.2870.2820.3540.353上海电力0.3360.5101.0000.3050.3390.2870.3820.397新和成0.2920.2810.3051.0000.4170.3330.3700.333大族激光0.2850.2870.3390.4171.0000.4060.3830.332传化股份0.3090.2820.2870.3330.4061.0000.2870.244中华人民共和国国贸0.4840.3540.3820.3700.3830.2871.0000.445上港集团0.3460.3530.3970.3330.3320.2440.4451.0002.2最优投资组合依照马科维茨资产组合理论拟定最优投资组合，该理论核心是在给定收益水平下，拟定资产组合中每只股票权重使风险最小；或者在给定风险承受能力状况下，选取权重使组合盼望收益最大。假设市场上仅有种资产，其收益率向量记为,投资者投资此种风险资产资产组合向量记为。两种资产收益率协方差记为，，其相应协方差矩阵记为。相应地，该资产组合收益率为，总风险为。最优投资组合含义是在给定预期收益水平下，风险最小投资方略，即：s.t.其中代表预期收益。数据解决过程见附录B2、B3，计算所得最优投资组合（表2.1）收益率方差（风险）为0.。表2.1最优投资组合（节选）股票名称万科A华能国际上海电力新和成大族激光传化股份中华人民共和国国贸上港集团组合权重0.11310.4434-0.02570.2366-0.13750.07890.14570.14571日收益率/1/50.0390.0270.0340.0540.0360.0520.0440.0363.917%/1/60.030.0130.0330.0120.0530.0370.0350.0231.504%/1/7-0.006-0.003-0.013-0.016-0.0050.02-0.0040.009-0.246%/1/80.006-0.008-0.019-0.049-0.0380.035-0.007-0.042-1.336%/1/9-0.0010.0240.0560.0190.0610.0680.0260.0091.555%/1/12-0.0120-0.0060.0030.012-0.021-0.0220.023-0.352%/1/13-0.035-0.022-0.025-0.043-0.037-0.056-0.058-0.02-3.384%注：①文中使用是无买空限制投资组合模型，因而也许浮现权重为负状况。②由于篇幅所限，上表仅列示了日收益率某些计算成果，完整时间跨度为1月日—7月2日。该投资组合日收益率和合计收益率分别如图2.1和图2.2所示：图2.1投资组合日收益率图2.2投资组合合计收益率3预期收益和风险报告采用CAPM模型计算投资组合预期收益率和风险。其中，为资产组合盼望收益率，为无风险资产收益率，代表资产组合系数，表达市场组合预期收益率。其中，是股票（组合）收益率与市场组合收益率协方差，为市场组合收益率方差。也即代表了股票（组合）收益率变动对市场组合收益率变动敏感度，因而用来衡量该股票（组合）系统性风险大小。本文使用沪深300指数代表市场组合，经计算可得，投资组合与市场组合协方差，市场组合收益率方差，投资组合，预期收益率见附录A表A4。图3.1组合收益率与预期收益率拟合图3.1显示了组合收益率和运用CAPM模型计算收益率预测值关系，可以看出两者存在较强有关性，两者有关系数为0.75954，拟合效果较好。下面探讨当大盘趋势不同步，股票组合收益率状况。（Ⅰ）大盘上升以3月16日到4月23日这段时间为例，沪深300指数收盘价从3705.671上涨到4740.892，上涨了31.049%。同一时间段，股票组合上涨幅度为26.9606%，略弱于大盘涨幅，但也有可观收益（Ⅱ）大盘下降以6月13日到7月2日这段时间为例，沪深300指数下跌幅度为-23.001%，同一时期股票组合收益率为-13.3968%，在一定限度上规避了系统性风险，有效减损。（Ⅲ）大盘振荡在5月7日到7月11日之间，沪深300指数处在震荡阶段。最高收盘价2191.855，最低2115.143，涨幅0.5%，收益率方差。同步期股票组合收益率为14.48%，几何平均收益率0.202%，可见在大盘振荡阶段，股票组合仍能维持稳定上涨态势。4风险对冲方略本文选取沪深300股指期货进行风险对冲，选取最小方差风险对冲模型计算套期保值比率，计算过程如下：，其中，为股票组合T日收盘价，为股票组合T+1日收盘价，为股指期货T日收盘价，股指期货为T+1日收盘价；为现货头寸，为期货头寸，代体现货收益率，为期货收益率，为套期保值比率。、和分别为现货组合、股指期货和风险对冲后组合收益率方差，为现货组合日收益率和期货日收益率协方差。为使最小化，有最优套期保值比率。经计算，沪深300股指期货日收益率方差,现货组合日收益率方差，两者协方差，有关系数，最优套期保值比率，即投资者在建立交易头寸时，沪深300股指期货合约价值与现货（股票组合）价值之比约为0.26，可以将风险控制到最小。（计算所需原始数据见附录A表A3、A4）

附录A：原始数据表A1：财务指标股票名称每股收益（EPS）净资产收益率（ROE）销售净利率（NPM）武钢股份0.763.481.09中华人民共和国国贸1.0810.6410.33首创股份0.169.896.02包钢股份0.191.061.51华能国际2.3116.1115.28日照港0.365.95.68上港集团0.4512.9611.4上海电力1.26156.67中华人民共和国石化1.278.147.27三一重工0.162.981.07万科A3.7819.0817.67世纪星源-0.016.67-8.75深振业A-0.8112.356.4零七股份0.05-8.6-7.6中华人民共和国宝安0.329.36-3.06南玻A0.6810.665.26沙河股份-1.69.548.65深康佳A-0.531.29-11.59深华发A0.422.770.83深赤湾A1.2710.3610.15新和成0.9911.911.16鸿达兴业-0.1113.0811.75伟星股份0.9813.4113.07华邦颖泰0.679.167.47德豪润达0.240.29-4.68精功科技0.121.5-10.65大族激光0.9118.3715.55天奇股份0.66.932.96传化股份0.6911.7910.64盾安环境0.173.672.15特锐德0.2811.7910.3南风股份0.735.623.5莱美药业-0.240.890.8汉威电子0.356.393.38立思辰0.319.357.6鼎汉技术0.2813.459.2新宁物流0.081.91.7亿纬锂能0.19.051.99银江股份-0.9413.048.36大禹节水0.054.363.32表A2：投资组合股票收益率（节选）日期万科A华能国际上海电力新和成大族激光传化股份中华人民共和国国贸上港集团/1/50.0390.0270.0340.0540.0360.0520.0440.036/1/60.0300.0130.0330.0120.0530.0370.0350.023/1/7-0.006-0.003-0.013-0.016-0.0050.020-0.0040.009/1/80.006-0.008-0.019-0.049-0.0380.035-0.007-0.042/1/9-0.0010.0240.0560.0190.0610.0680.0260.009/1/12-0.0120.000-0.0060.0030.012-0.021-0.0220.023/1/13-0.035-0.022-0.025-0.043-0.037-0.056-0.058-0.020/1/140.0490.0110.0230.0370.0630.0480.0240.026/1/15-0.006-0.0190.009-0.0020.0730.0150.000-0.006/1/160.0100.0100.0000.044-0.026-0.0360.024-0.017/1/190.014-0.0140.019-0.004-0.046-0.0170.0070.099/1/20-0.014-0.0040.0000.0220.0330.0210.011-0.011/1/210.026-0.020-0.0120.075-0.0190.0260.066-0.024/1/220.0040.0030.0060.0300.0220.004-0.0160.011/1/23-0.013-0.025-0.006-0.012-0.023-0.016-0.020-0.008注：①收益率计算方式：

当天收益率②实际计算中使用收益率时间跨度为1月5日—7月2日，限于文章篇幅，表A2中仅列示每只股票首月日收益率。表A3：投资组合日收益率（节选）日期收益率日期收益率日期收益率/1/53.917%/2/23.415%/3/22.002%/1/61.504%/2/32.388%/3/3-0.328%/1/7-0.246%/2/40.974%/3/44.220%/1/8-1.336%/2/50.099%/3/52.727%/1/91.555%/2/63.246%/3/6-1.390%/1/12-0.352%/2/91.650%/3/9-4.050%/1/13-3.384%/2/100.976%/3/101.552%/1/142.090%/2/112.139%/3/11-1.568%/1/15-1.969%/2/12-0.973%/3/12-0.802%/1/161.788%/2/133.524%/3/13-0.967%/1/191.435%/2/161.307%/3/160.058%/1/20-0.125%/2/17-1.465%/3/172.809%/1/212.295%/2/18-3.818%/3/18-0.699%/1/220.524%/2/193.883%/3/192.035%/1/23-1.722%/2/203.637%/3/20-0.827%/2/231.071%/3/230.957%/2/24-4.958%/3/242.773%/2/25-1.507%/3/25-4.382%/2/26-4.853%/3/262.585%/2/271.260%/3/271.288%/3/300.970%/3/31-1.342%注：实际计算收益率时间跨度为1月5日—7月2日，上表仅列示组合前三个月日收益率。表A4：资产组合预期收益率（节选）日期市场组合收益率（沪深300）资产组合预期收益率/1/53.5891%3.0479%/1/63.1780%1.0361%/1/7-0.5981%-0.2517%/1/8-2.2365%-1.5427%/1/91.6091%1.4876%/1/120.1209%0.3150%/1/13-2.3169%-1.6060%/1/144.2133%3.5398%/1/15-0.0189%0.2047%/1/161.8078%1.6442%/1/191.1185%1.1010%/1/200.6321%0.7177%/1/21-0.1718%0.0842%/1/221.1297%1.1099%/1/23-0.5806%-0.2378%注：实际计算收益率时间跨度为1月5日—7月2日，上表仅列示首月日收益率。表A5：沪深300股指期货日收盘价和收益率（节选）日期收盘价（元）收益率-04-183512—-04-193262-7.118%-04-2032920.920%-04-2133481.701%-04-223316-0.956%-04-233304-0.362%-04-2633120.242%-04-273262.8-1.486%-04-283255-0.239%-04-293215-1.229%-04-303243.20.877%注：实际计算中使用收益率时间跨度为4月19日—7月2日，限于文章篇幅，上表仅列示首月日收益率。

附录B：代码B1：择股library(xts)library(xtsExtra)library(quantmod)library(RODBC)fin<-read.csv("C:/Users/lirj/Desktop/期末考试/财务.csv"，header=TRUE)sort(table(fin$Nnindcd))fin.lrj<-fin[fin$Nnindcd=="CSRC"，]stk.lrj.list<-fin.lrj[，"Stkcd"]stk.lrj.num<-length(stk.lrj.list)stk.lrj.numrank.lrj<-apply(fin.lrj，MARGIN=2，FUN=function(var)stk.lrj.num-rank(var)+1)rank.lrj<-data.frame(rank.lrj)k<-16win.select<-ifelse(rank.lrj$EPS<=k&rank.lrj$ROE<=k&rank.lrj$NPM<=k+4，1，0)win.selectfin.win.lrj<-fin.lrj[win.select==1，]fin.lose.lrj<-fin.lrj[win.select==0，]win.list<-fin.win.lrj$Stkcdlose.list<-fin.lose.lrj$Stkcdwin.nme<-fin.win.lrj$Stknmelose.nme<-fin.lose.lrj$Stknmewin.nmelose.nmeB2：投资组合风险-收益计算library(RODBC)library(xts)library(xtsExtra)library(quantmod)library(TTR)library(fPortfolio)ret<-read.csv("C:/Users/lirj/Desktop/期末考试/股票收益.csv"，header=TRUE)head(ret)ret<-xts(ret[，-1]，order.by=as.Date(ret[，1]))na.num<-apply(is.na(ret)，MARGIN=2，sum)na.numret<-na.omit(ret)cumret<-cumprod(1+ret)plot.xts(ret，main="DailyReturnof8Stocks(Y-Y)")plot.xts(cumret，screens=1，main="CumulativeReturnof8Stocks(Y-Y)")head(data.frame(ret))pairs(data.frame(ret)，pch=20，col='darkblue',main="Correlationsamong8StocksDailyReturns")cov.mat<-cov(ret)cov.matcorr.mat<-cor(ret)corr.matw.g<-c(0.1131,0.443,-0.0257,0.2366,-0.1375,0.0789,0.1456,0.1457)ret.g<-ret%*%w.gret.g<-xts(ret.g，order.by=index(ret))cumret.g<-cumprod(1+ret.g)mean(ret.g)var(ret.g)w.g%*%cov(ret)%*%w.gplot.xts(cumret.g)B3：最优投资组合求解min.var<-function(ret.exp，mu){n<-dim(ret)[2]Q<-cov(ret)r<-ret.expL1<-cbind(Q，rep(1，n)，r)L2<-rbind(c(rep(1，n)，0，0)，c(r，0，0))L<-rbind(L1，L2)b<-c(rep(0，n)，1，mu)solve.res<-solve(L，b)wt<-solve.res[1:n]ret.mu<-r%*%wtret.var<-wt%*%Q%*%wtreturn(c(ret.mu，ret.var，wt))}ret.exp<-apply(ret，MARGIN=2，FUN=mean)step<-seq(-0.0002，0.0012，by=0.000002)frontier.res<-t(sapply(step，FUN=funct