matlab实验报告模板.docx
《matlab实验报告模板.docx》由会员分享，可在线阅读，更多相关《matlab实验报告模板.docx（25页珍藏版）》请在冰豆网上搜索。
matlab实验报告模板
MATLAB高数实验
主题：
投资组合问题
指导教师:
阮小娥老师
制作时间：
2010.12—2011.01
学院：
机械学院
小组成员：
机自07赵磊2010011180
机自07周策2010011181
机自07邹业兵2010011182
一、引言:
现在社会，随着社会的发展和生产力的提升，人们的生活水平日益提高。
人们不仅仅满足于吃饱穿暖，人们对财富的追求也愈来愈强烈，此时股票作为社会化大生产的产物，股票至今已有将近似400年的历史，它的发展与人类经济的发展可以说是和谐统一的，而它集惊险刺激和可能存在的巨额经济效益的特点更是牢牢的牵住了现代人的心，毕竟，在已经不愁衣食的前提下，谁不想像股神巴菲特一样，动动脑筋便收拢全世界的财富呢？
所以，研究股票自然而然就成了现代人生活中极为重要的一部分。
二、实验问题：
企业（或政府）或个人进行投资的目的是为了取得利润，在现代市场经济中，由于投资环境瞬息万变，任何一项投资的收益总是不能确定的，也就是说具有一定的风险，如何度量收益和风险呢？
下面给出一个实际问题进行分析：
表9-1给出3只股票（A、B、C股）49个周末的收盘价，
表9-1股票49个周末的收盘价
周数
A股
B股
C股
周数
A股
B股
C股
1
10.3
0.77
10.6
25
11.8
1.2
14.3
2
10.6
0.81
11.2
26
11.8
1.2
14.7
3
10.8
0.89
11.9
27
12.2
1.24
14.9
4
10.7
0.93
11.7
28
12.4
1.32
15.6
5
10.8
0.90
11.6
29
12.4
1.32
15.5
6
10.8
0.93
11
30
12.7
1.29
15.8
7
11.2
1.01
12.1
31
12.5
1.32
15.4
8
11.3
0.96
13.2
32
12.3
1.32
14.9
9
11
1
12.6
33
12.6
1.33
15.3
10
10.9
1.04
12.6
34
12.4
1.37
15.6
11
10.9
1.04
12.6
35
12.6
1.48
14.2
12
10.6
1.09
12.6
36
12.1
1.49
13.8
13
11.5
1.1
12.7
37
11.9
1.51
14
14
11.5
1.26
13.6
38
11.9
1.32
13.8
15
11.8
1.17
14.2
39
12.4
1.32
14.2
16
12.2
1.06
14.7
40
12.1
1.3
13.8
17
12
1.13
14.4
41
12
1.28
14
18
12.1
1.09
14.3
42
12.4
1.17
14.3
19
12.3
1.03
15
43
12
1.14
13.6
20
12.1
1.08
14.9
44
12
1.08
14.1
21
11.4
1.08
13.4
45
12
1.09
14.4
22
11.9
1.12
14
46
12.1
1.07
14.6
23
11.9
1.12
14.4
47
12
1.2
14.4
24
12
1.15
14.7
48
12
1.24
14.4
49
12.1
1.26
14.5
问题：
（1）分别计算这三只股票的周收益率的时间序列，平均收益率和他们收益率分布的方差；
（2）允许卖空的条件下给出这3只股票风险最小的投资策略。
三、问题分析
（1）、已知和股票有关的一些概念（实验所涉及的理论知识）：
记
为证券A在第t周末的价格，定义
为证券A当周的证券收益率，由于投资初期不可能知道将来的收益率，因此收益率是一个随机变量，
为证券A在某一时间段内收益率分布，且定义该收益率分布的均值
为证券A的期望收益率，并用该收益率分布的方差
度量该证券A期望收益的偏差称为风险，同样可以计算证券B的期望收益率
及方差
两证券收益率分布间的协方差为
协方差给出了证券A和证券B收益率分布相互影响的程度，协方差的大小依赖于收益所的
单位，相关系数
反映了证券A和证券B投资比率分
为多少时风险比率最小，受益最大，这就是证券A和证券B的投资组合问题。
（2）、数学建模：
用决策变量X1，X2，X3分别表示投资股票A，B，C的投资比率，记股票A，B，C的投资比率，记股票A，B，C的收益率分别为Ea、Eb、Ec，根据题意,Ea、Eb、Ec是随机变量，投资的总收益率为
S=X1*Ea+X2*Eb+X3*Ec
也是随机变量，用E和D分别表示随机变量的数学期望和方差，则有
（以下结果是经过第一、二、三、四步操作所得）
Ea=0.0036，Eb=0.0121,Ec=0.0073
Da=5.4449e-004,Db=0.0036,Dc=0.0016
Dab=3.8817e-005、Dac=5.6361e-004、Dbc=1.5600e-004；
由概率论知识可知
covAC=Dac*Da^0.5*Dc^0.5
covAB=Dab*Da^0.5*Db^0.5
covBC=Dbc*Db^0.5*Dc^0.5
投资的总期望收益率为：
y=E=x1*Ea+x2*Eb+x3*Ec
用总收益的方差衡量投资风险是，投资总收益的方差为：
Z=D（x1*Ea+x2*Eb+x3*Ec）
=Da*x1^2+Db*x2^2+Dc*x3^2
+2covAB*x1*x2+2covAC*x1*x3+2covBC*x2*x3
由于x1、x2、x3分别表示投资股票A、B、C的投资比例，所以
x1+x2+x3=1
数学模型：
Minz（x1,x2,x3）=5.4449e-004*x1^2+0.0036*x2^2+0.0016*x3^2
+11.8000e-008*x1*x2+7.2e-0088x1x3+18.39e-008*x2*x3
x1+x2+x3=1
（3）、求解方法：
第一步：
画出A、B、C三只股票49个周末收盘价的变化情况的图像（便于分析问题）；
第二步：
求出A、B、C三只股票49个周末收盘价的期望收益率Ea、Eb、Ec；
第三步：
求出A、B、C三只股票49个周末收盘价的方差Da、Db、Dc；
第四步：
求出A、B、C三只股票收益率分布直间的协方差：
Dab、Dac、Dbc；
第五步：
求出投资风险最小和此时的收益率。
四、程序设计：
第一步：
1、A股：
x=1:
1:
49;
y=[10.3,10.6,10.8,10.7,10.8,10.8,11.2,11.3,11,10.9,10.6,10.9,11.5,11.5,11.8,12.2,12,12.1,12.3,12.1,11.4,11.9,11.9,12,11.8,11.8,12.2,12.4,12.4,12.7,12.5,12.3,12.6,12.4,12.6,12.1,11.9,11.9,12.4,12.1,12,12.4,12,12,12,12.1,12,12,12.1];
plot（x,y）
2、B股：
x=1:
1:
49;
y=[0.77,0.81,0.89,0.93,0.9,0.93,1.01,0.96,1,1.04,1.04,1.09,1.1,1.26,1.17,1.06,1.13,1.09,1.03,1.08,1.08,1.12,1.12,1.15,1.2,1.2,1.24,1.32,1.32,1.29,1.32,1.32,1.33,1.37,1.48,1.49,1.51,1.32,1.32,1.3,1.28,1.17,1.14,1.08,1.09,1.07,1.2,1.24,1.26]
plot（x,y）
3、C股：
x=1:
1:
49;
y=[10.6,11.2,11.9,11.7,11.6,11,12.1,13.2,12.6,12.6,12.6,12.6,12.7,13.6,14.2,14.7,14.4,14.3,15,14.9,13.4,14,14.4,14.7,14.3,14.7,14.9,15.6,15.5,15.8,15.4,14.9,15.3,15.6,14.2,13.8,14,13.8,14.2,13.8,14,14.3,13.6,14.1,14.4,14.6,14.4,14.4,14.5]
plot（x,y）
第二、三、四步：
PA=[10.310.610.810.710.810.811.211.31110.910.610.911.511.511.812.21212.112.312.111.411.911.91211.811.812.212.412.412.712.512.312.612.412.612.111.911.912.412.11212.412121212.1121212.1];
PB=[0.770.810.890.93.9.931.1.9611.041.041.091.11.261.171.061.131.091.131.081.081.121.121.151.21.21.241.321.321.291.321.321.331.371.481.491.511.321.321.31.281.171.141.081.091.071.21.241.26];
PC=[10.611.211.911.711.61112.113.212.612.612.612.612.713.614.214.714.414.31514.913.41414.414.714.314.714.915.615.515.815.414.915.315.614.213.81413.814.213.81414.313.614.114.414.614.414.414.5];
fori=2:
1:
49
ra（i-1）=（PA（i）-PA（i-1））/（PA（i-1））;
rb（i-1）=（PB（i）-PB（i-1））/（PB（i-1））;
rc（i-1）=（PC（i）-PC（i-1））/（PC（i-1））;
end
rA=[ra];
rB=[rb];
rC=[rc];
T=48;
Ea=sum（rA（1:
T））/T;
Eb=sum（rB（1:
T））/T;
Ec=sum（rC（1:
T））/T;
Dab=（1/（T-1））*sum（（rA（1:
T）-Ea）.*（rB（1:
T）-Eb））;
Dac=（1/（T-1））*sum（（rA（1:
T）-Ea）.*（rC（1:
T）-Ec））;
Dbc=（1/（T-1））*sum（（rC（1:
T）-Ec）.*（rB（1:
T）-Eb））;
Da=sum（（rA（1:
T）-Ea）.^2）/T;
Db=sum（（rB（1:
T）-Eb）.^2）/T;
Dc=sum（（rC（1:
T）-Ec）.^2）/T;
covAB=Dab*（（Da）^0.5）*（（Db）^0.5）;
covAC=Dab*（（Da）^0.5）*（（Dc）^0.5）;
covBC=Dab*（（Db）^0.5）*（（Dc）^0.5）;
第五步：
h=[2*Da^22*covAB2*covAC;2*covAB2*Db^22*covBC;2*covAC2*covBC2*Dc^2];
c=[0,0,0];
a1=-[Ea,Eb,Ec];
b1=0;
A2=[1,1,1];
b2=1;
[x,fv,ef,out,lag]=quadprog（h,c,a1,b1,A2,b2）
F=sqrt（fv）
S=-a1*x
五、程序运行结果（实验结果）为：
第一步：
A股49个周末的收盘价的变化情况
.B股49个周末的收盘价的变化情况
C股49个周末的收盘价的变化情况
第二、三、四步：
Ea=
0.0036
Eb=
0.0121
Ec=
0.0073
Dab=
3.8817e-005
Dac=
5.6361e-004
Dbc=
1.5600e-004
Da=
5.4449e-004
Db=
0.0036
Dc=
0.0016
covAB=
5.4000e-008
covAC=
3.5989e-008
covBC=
9.1950e-008
第五步：
x=
0.8892
0.0167
0.0940
fv=
2.6792e-007
ef=
1
out=
iterations:
1
constrviolation:
-0.2826
algorithm:
'medium-scale:
active-set'
firstorderopt:
6.3527e-022
cgiterations:
[]
message:
'Optimizationterminated.'
lag=
lower:
[3x1double]
upper:
[3x1double]
eqlin:
-5.3584e-007
ineqlin:
0
F=
5.1761e-004
S=
0.0041
由此可知
（1）A、B、C的平均收益率分别为Ea=0.0036，Eb=0.0121,Ec=0.0073
收益率分布的方差分别为Da=5.4449e-004,Db=0.0036,Dc=0.0016
（2）、投资A、B、C股票的投资比例分别为：
88.92%、1.67%、9.40%，
此时风险最小为0.05176%，收益率为0.41%
六、问题的进一步拓展与实验
类似地，还可以引入四种以上数量的股票及股票不同时间的投资问题。
一般来说，投资人希望投资风险小而收益大，是一个多目标的优化问题。
我们可以将风险和收益的变化关系图像画出来。
给投资者最感官的认识。
最佳投资组合
示例5如表8-7给出4只股票在同一时期内的平均收益率r%，购买股票是交易率为p%，风险损失率q%，投资越分散，总的风险越小。
假定总风险用投资中最大一种股票的风险来度量，且同期银行存款利率是r%=5%，既无交易费有无风险，由于投资者承担风险的程度不一样，若给定风险一个界限a，季最大风险不超过a的情况下，为投资者建议一种投资策略，使其尽可能获得最大收益。
表8-7股票信息
股票
r
P
q
1
28
2.5
1
2
21
1.5
2
3
23
5.5
4.5
4
25
2.6
6.5
实验过程
为了建立最佳投资组合模型。
我们分以下三个步骤来完成次实验。
1：
条件假设
假设投资四种股票的资金分别为x,总投资金额为M.并设四种购票之间是相互独立的，且在投资的同一时间内都为定值。
不受意外因素影响。
2：
建立模型
投资四种股票的风险度为，i=1，2，3，4；购买四种股票时所付交易费分别为，i=1，2，3，4，则购买四种股票的收益分别为（），i=1，2，3，4，为使投资者获的更大利益，在总风险度不超过a的情况下，可建立如下模型：
max
；
s.t.
=M
/M≤a,i=1,2,3,4
≥0,i=0,1,2,3
3：
模型求解
模型简化为
max
，
s.t.
=1
≤a，i=0,1,2,3，4
≥0,i=0,1,2,3，4
其中
=
/M,i=0,1,2,3，4
将表8-7给定的数据带入模型的
minF=-0.05
-0.27
-0.19
-0.185
-0.185
s.t.
+1.01
+1.02
+1.045
+1.065
=1
0.025
≤a
0.015
≤a
0.055
≤a
0.026
≤a
≥0,i=0,1,2,3，4
如何给定风险度a没有一定的准则，不同的投资者承受风险度的能力不同，本次试验从a=0开始，以步长△a=0.001编写下面的程序进行循环搜索。
程序设计
a=0;
while（1.1-a）>1
c=[-0.05,-0.27,-0.19,-0.185,-0.185];
aeq=[1,1.01,1.02,1.045,1.065];
beq=[1];
A=[0,0.025,0,0,0;0,0,0.015,0,0;0,0,0,0.055,0;0,0,0,0,0.026];
b=[a;a;a;a]
vlb=[0,0,0,0,0];vub=[];
[x,val]=linprog（c,A,b,aeq,beq,vlb,vub）;
a%风险度
x=x’
Q=-val%风险度对应的收益率
plot（a,Q,’,’）
axis（[0,0.1,0,0.5]）
holdon
a=a+0.001;
end
xlabel（‘a’）,ylabel（‘Q’）
程序运行结果
>>a=0;
while（1.1-a）>1
c=[-0.05,-0.27,-0.19,-0.185,-0.185];
aeq=[1,1.01,1.02,1.045,1.065];
beq=[1];
A=[0,0.025,0,0,0;0,0,0.015,0,0;0,0,0,0.055,0;0,0,0,0,0.026];
b=[a;a;a;a]
vlb=[0,0,0,0,0];vub=[];
[x,val]=linprog（c,A,b,aeq,beq,vlb,vub）;
a
x=x'
Q=-val
plot（a,Q,'.'）
axis（[0,0.1,0,0.5]）
holdon
a=a+0.001;
end
xlabel（'a'）,ylabel（'Q'）
b=
0
0
0
0
Optimizationterminated.
a
0
x=
1.00000.00000.00000.00000.0000
Q=
0.0500
b=
1.0e-003*
1.0000
1.0000
1.0000
1.0000
Optimizationterminated.
a=
1.0000e-003
x=
0.83160.04000.06670.01820.0385
Q=
0.0755
b=
0.0020
0.0020
0.0020
0.0020
Optimizationterminated.
a=
0.0020
x=
0.66330.08000.13330.03640.0769
Q=
0.1011
b=
0.0030
0.0030
0.0030
0.0030
Optimizationterminated.
a=
0.0030
x=
0.49490.12000.20000.05450.1154
Q=
0.1266
b=
0.0040
0.0040
0.0040
0.0040
Optimizationterminated.
a=
0.0040
x=
0.32660.16000.26670.07270.1538
Q=
0.1521
b=
0.0050
0.0050
0.0050
0.0050
Optimizationterminated.
a=
0.0050
x=
0.15820.20000.33330.09090.1923
Q=
0.1776
b=
0.0060
0.0060
0.0060
0.0060
Optimizationterminated.
a=
0.0060
x=
0.00000.24000.40000.10910.2212
Q=
0.2019
b=
0.0070
0.0070
0.0070
0.0070
Optimizationterminated.
a=
0.0070
x=
0.00000.28000.46670.12730.1016
Q=
0.2066
b=
0.0080
0.0080
0.0080
0.0080
Optimizationterminated.
a=
0.0080
x=
0.00000.32000.53330.12710.0000
Q=
0.2112
……………
………………………………………………………………….
……………………………………………………………………
…………………………………………………………………………
b=
0.0970
0.0970
0.0970
0.0970
Optimizationterminated.
a=
0.0970
x=
0.00000.99010.00000.00000.0000
Q=
0.2673
b=
0.0980
0.0980
0.0980
0.0980
Optimizationterminated.
a=
0.0980
x=
0.00000.99010.00000.00000.0000
Q=
0.2673
b=
0.0990
0.0990
0.0990
0.0990
Optimizationterminated.
a=
0.0990
x=
0.00000.99010.00000.00000.0000
Q=
0.2673
风险与投资的关系图像
1）、风险大，收益也大；
2）、当投资越分散时，投资者承担的风险越小，这与题意一致，及冒险的投资者会出现集中投资的情况，保守的投资者则尽量分散投资。
3）上图曲线上的点表示当风险水平给定时能获得的最大收益率。
4）…………………..
七、实验的总结与体会：
赵磊：
通过这次实验，我觉得收获有两点：
一、掌握了Matlab的基本操作知识以及其在生活中的很简单应用。
二、感觉到三个人合作的重要性。
在做实验的过程中，我们各有分工，但同时我们又相互联系确保实验的同步进行，最后终于完成了，很有成就感！
周策：
数学是一门伟大的学科！
我想，了解数学的人是不会反
升级会员 