matlab实验报告模板.docx
- 文档编号:23228825
- 上传时间:2023-05-15
- 格式:DOCX
- 页数:22
- 大小:45.32KB
matlab实验报告模板.docx
《matlab实验报告模板.docx》由会员分享,可在线阅读,更多相关《matlab实验报告模板.docx(22页珍藏版)》请在冰豆网上搜索。
matlab实验报告模板
MATLAB高数实验
主题:
投资组合问题
指导教师:
阮小娥老师
制作时间:
—
学院:
机械学院
小组成员:
机自07赵磊80
机自07周策81
机自07邹业兵82
一、引言:
现在社会,随着社会的发展和生产力的提升,人们的生活水平日益提高。
人们不仅仅满足于吃饱穿暖,人们对财富的追求也愈来愈强烈,此时股票作为社会化大生产的产物,股票至今已有将近似400年的历史,它的发展与人类经济的发展可以说是和谐统一的,而它集惊险刺激和可能存在的巨额经济效益的特点更是牢牢的牵住了现代人的心,毕竟,在已经不愁衣食的前提下,谁不想像股神巴菲特一样,动动脑筋便收拢全世界的财富呢?
所以,研究股票自然而然就成了现代人生活中极为重要的一部分。
二、实验问题:
企业(或政府)或个人进行投资的目的是为了取得利润,在现代市场经济中,由于投资环境瞬息万变,任何一项投资的收益总是不能确定的,也就是说具有一定的风险,如何度量收益和风险呢?
下面给出一个实际问题进行分析:
表9-1给出3只股票(A、B、C股)49个周末的收盘价,
表9-1股票49个周末的收盘价
周数
A股
B股
C股
周数
A股
B股
C股
1
25
2
26
3
27
4
28
5
29
6
11
30
7
31
8
32
9
11
1
33
10
34
11
35
12
36
13
37
14
14
38
15
39
16
40
17
12
41
12
14
18
42
19
15
43
12
20
44
12
21
45
12
22
14
46
23
47
12
24
12
48
12
49
问题:
(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=,Eb=,Ec=
Da=,Db=,Dc=
Dab=、Dac=、Dbc=;
由概率论知识可知
covAC=Dac*Da^*Dc^
covAB=Dab*Da^*Db^
covBC=Dbc*Db^*Dc^
投资的总期望收益率为:
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)=*x1^2+*x2^2+*x3^2
+*x1*x2++*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=[,,,,,,,,11,,,,,,,,12,,,,,,,12,,,,,,,,,,,,,,,,,12,,12,12,12,,12,12,];
plot(x,y)
2、B股:
x=1:
1:
49;
y=[,,,,,,,,1,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,]
plot(x,y)
3、C股:
x=1:
1:
49;
y=[,,,,,11,,,,,,,,,,,,,15,,,14,,,,,,,,,,,,,,,14,,,,14,,,,,,,,]
plot(x,y)
第二、三、四步:
PA=[111212121212121212];
PB=[.9.93.961];
PC=[1115141414];
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)^*((Db)^;
covAC=Dab*((Da)^*((Dc)^;
covBC=Dab*((Db)^*((Dc)^;
第五步:
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=
Eb=
Ec=
Dab=
Dac=
Dbc=
Da=
Db=
Dc=
covAB=
covAC=
covBC=
第五步:
x=
fv=
ef=
1
out=
iterations:
1
constrviolation:
algorithm:
'medium-scale:
active-set'
firstorderopt:
cgiterations:
[]
message:
'Optimizationterminated.'
lag=
lower:
[3x1double]
upper:
[3x1double]
eqlin:
ineqlin:
0
F=
S=
由此可知
(1)A、B、C的平均收益率分别为Ea=,Eb=,Ec=
收益率分布的方差分别为Da=,Db=,Dc=
(2)、投资A、B、C股票的投资比例分别为:
%、%、%,
此时风险最小为%,收益率为%
六、问题的进一步拓展与实验
类似地,还可以引入四种以上数量的股票及股票不同时间的投资问题。
一般来说,投资人希望投资风险小而收益大,是一个多目标的优化问题。
我们可以将风险和收益的变化关系图像画出来。
给投资者最感官的认识。
最佳投资组合
示例5如表8-7给出4只股票在同一时期内的平均收益率r%,购买股票是交易率为p%,风险损失率q%,投资越分散,总的风险越小。
假定总风险用投资中最大一种股票的风险来度量,且同期银行存款利率是r%=5%,既无交易费有无风险,由于投资者承担风险的程度不一样,若给定风险一个界限a,季最大风险不超过a的情况下,为投资者建议一种投资策略,使其尽可能获得最大收益。
表8-7股票信息
股票
r
P
q
1
28
1
2
21
2
3
23
4
25
实验过程
为了建立最佳投资组合模型。
我们分以下三个步骤来完成次实验。
1:
条件假设
假设投资四种股票的资金分别为x,总投资金额为M.并设四种购票之间是相互独立的,且在投资的同一时间内都为定值。
不受意外因素影响。
2:
建立模型
投资四种股票的风险度为,i=1,2,3,4;购买四种股票时所付交易费分别为,i=1,2,3,4,则购买四种股票的收益分别为(),i=1,2,3,4,为使投资者获的更大利益,在总风险度不超过a的情况下,可建立如下模型:
max
;
.
=M
/M≤a,i=1,2,3,4
≥0,i=0,1,2,3
3:
模型求解
模型简化为
max
,
.
=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=
.
+
+
+
+
=1
≤a
≤a
≤a
≤a
≥0,i=0,1,2,3,4
如何给定风险度a没有一定的准则,不同的投资者承受风险度的能力不同,本次试验从a=0开始,以步长△a=编写下面的程序进行循环搜索。
程序设计
a=0;
while>1
c=[,,,,];
aeq=[1,,,,];
beq=[1];
A=[0,,0,0,0;0,0,,0,0;0,0,0,,0;0,0,0,0,];
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,])
holdon
a=a+;
end
xlabel(‘a’),ylabel(‘Q’)
程序运行结果
>>a=0;
while>1
c=[,,,,];
aeq=[1,,,,];
beq=[1];
A=[0,,0,0,0;0,0,,0,0;0,0,0,,0;0,0,0,0,];
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,])
holdon
a=a+;
end
xlabel('a'),ylabel('Q')
b=
0
0
0
0
Optimizationterminated.
a
0
x=
Q=
b=
*
Optimizationterminated.
a=
x=
Q=
b=
Optimizationterminated.
a=
x=
Q=
b=
Optimizationterminated.
a=
x=
Q=
b=
Optimizationterminated.
a=
x=
Q=
b=
Optimizationterminated.
a=
x=
Q=
b=
Optimizationterminated.
a=
x=
Q=
b=
Optimizationterminated.
a=
x=
Q=
b=
Optimizationterminated.
a=
x=
Q=
……………
………………………………………………………………….
……………………………………………………………………
…………………………………………………………………………
b=
Optimizationterminated.
a=
x=
Q=
b=
Optimizationterminated.
a=
x=
Q=
b=
Optimizationterminated.
a=
x=
Q=
风险与投资的关系图像
1)、风险大,收益也大;
2)、当投资越分散时,投资者承担的风险越小,这与题意一致,及冒险的投资者会出现集中投资的情况,保守的投资者则尽量分散投资。
3)上图曲线上的点表示当风险水平给定时能获得的最大收益率。
4)…………………..
七、实验的总结与体会:
赵磊:
通过这次实验,我觉得收获有两点:
一、掌握了Matlab的基本操作知识以及其在生活中的很简单应用。
二、感觉到三个人合作的重要性。
在做实验的过程中,我们各有分工,但同时我们又相互联系确保实验的同步进行,最后终于完成了,很有成就感!
周策:
数学是一门伟大的学科!
我想,了解数学的人是不会反对的。
在这次实验中,我们选择了比较具有实际意义的这个关于股票的问题,总的来说,这个题目与实际的契合度比较高,所以在做的过程中我们也遇到了不少困难。
首先,为了提高模型的准确度我们必须使用一定数量的数据。
再者就是数据拟合时的编程较为困难。
但是,在我们的努力下,在我们花费了大量的时间和努力之后终于做出来了。
在做成之后我们都很开心,我们坚信这是团队的力量,是合作的结果。
我想,只要我们不惧困难,团结一心,我们就能发挥百分之百的力量把事情做好。
邹业兵:
总结:
在整个程序中,根据不同的问题,我们将程序分成了三个模块(实验的第一步;第二、三、四步;第五步),分别交给本组的三个成员来编写。
我负责最为基础也是最为复杂的第二、三、四步。
编写过程中,一开始根据所学知识,在心中定下了程序的框架,但在实际编程时发现数据太多,一个一个编程计算很麻烦。
后来经过思考、查阅有关matlab使用的书籍,并和其他成员集体讨论,发现可以将数据用矩阵的形式来表示并用相关的语句顺利地解决了这个问题,大大地简化了程序的过程并减少了编程的工作量。
体会:
和小组成员一起完成的实验,可以说也算是一个项目,也就是说我们完成了一个项目,这对我们的自信心有一个很大的提高。
从一开始对编程的一知半解到基本掌握一些编程技巧,让我真正体会到了什么叫进步,什么叫提高。
现代的社会可以说是团队的社会,一个人想要干出一番大事业几乎是不可能的是,通过这次大作业我真正体会到了团队的作用,团队的凝聚力。
此外我还意识到在学习的过程中遇到问题要积极思考,不怕困难,要学会通过各种手段(例如:
和他人讨论,查阅书籍、询问老师学长)来解决问题。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 实验 报告 模板
![提示](https://static.bdocx.com/images/bang_tan.gif)