1、1安徽财经大学金融证券实验室实 验 报 告实验课程名称 金融 MATLAB 开 课 系 部 金融学院 班 级 12 金融 1 班 学 号 20125633 姓 名 宋歌 指 导 教 师 洪振木 2014 年 5 月 12 日2实验名称 MATLAB 数据分析学院 金融学院 学号 20125633 姓名 宋歌实验准备实 验 目 的 学会 MATLAB 与 Excel 的数据交换； 学会使用 MATLAB 进行现金流分析； 学会使用 MATLAB 计算概率密度、分布函数值，学会生成各种随机数； 掌握 MATLAB 中的 CFTOOL 数据拟合工具； 学会如何使用 MATLAB 来求方程和方程组的数
2、值解。实 验 设 计 方 案1. 利用 MATLAB 读取 Chapter3_code 文件夹中 StockPriceData.xls 的数据，分别找出所有最高价中的最大值和所有最低价中的最小值，将最大值和最小值分别写入StockPriceData.xls 文件的 sheet2 工作表中；并在同一图像中画出最高价曲线和最低价曲线。2. 投资人贷款 65W 买房，还款期 20 年，每月还 4000 元，则贷款利率为多少？ 若改为每月还 4500 元，贷款利率不变，则还贷期限为多长？（提示：年金现金流计算）3. 求参数为 n=50，p=0.4 的二项分布在各点处的密度函数值与分布函数值，然后画出该
3、二项分布的密度函数曲线；调用 random 函数生成 100001 的二项分布（n=50，p=0.4 ）的随机数向量，然后作出频率直方图。4. 使用 cftool 工具对下面的数据进行拟合，并画出拟合曲线。 （至少使用两种曲线进行拟合） 。Xdata 1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 1890Ydata 3.9 5.3 7.2 9.6 12.9 17.1 23.1 31.4 38.6 50.2 62.9Xdata 1900 1910 1920 1930 1940 1950 1960 1970 1980 1990Ydata 76.
4、0 92.0 105.7 122.8 131.7 150.7 179.0 205.0 226.5 248.75. （1）求方程 的零点；（2）用 fsolve 解下面的非线性方程组在（ 1,1,1）附近的根004sin2xyzex)3)(cos)(2tettft31. 利用 MATLAB 读取 Chapter3_code 文件夹中 StockPriceData.xls 的数据，分别找出所有最高价中的最大值和所有最低价中的最小值，将最大值和最小值分别写入StockPriceData.xls 文件的 sheet2 工作表中；并在同一图像中画出最高价曲线和最低价曲线。解：1）读取数据data,tex
5、tdate=xlsread(C:UsersAdministratorDesktopChapter3_codeStockPriceData.xls)2）找出所有最高价中的最大值和所有最低价中的最小值a=max(xlsread(StockPriceData.xls, 1,C2:C511)a = 3.2926e+003 status1,message1=xlswrite(StockPriceData.xls, a, sheet2,A1)status1 = 1message1 = message: identifier: b=min(xlsread(StockPriceData.xls,1, D2:D
6、511)b = 1.0005e+003 status2,message2=xlswrite(StockPriceData.xls, b, sheet2,A2)status2 = 1message2 = message: identifier: 3.在图像中画出最高价曲线和最低价曲线。textdate=0:509Y1=xlsread(StockPriceData.xls, 1, C2:C511) Y2=xlsread(StockPriceData.xls, 1,D2:D511) plot(textdate,Y1,r,textdate,Y2,g)4实验结论 52. 投资人贷款 65W 买房，还款期
7、 20 年，每月还 4000 元，则贷款利率为多少？ 若改为每月还 4500 元，贷款利率不变，则还贷期限为多长？（提示：年金现金流计算）解：第一问是年金现金流计算，可以用 annurate 函数，NumPeriods=20*12 Payment=4000PresentValue=650000FutureValue=0; Due=0结果如下rate=annurate(20*12,4000,650000,0,0)rate = 0.0035年利率 a=rate*12a = 0.0418第二问是计算年金周期函数，可以用 annuterm 函数计算Rate=0.0035 Payment=4000 Pr
8、esentValue=650000FutureValue=0Due=0结果如下6numperiods=annuterm(0.0035,-4500,650000,0,0)numperiods = 201.5864运行结果 numperiods = 201.5864（月）=16.7989（年）3. 求参数为 n=50，p=0.4 的二项分布在各点处的密度函数值与分布函数值，然后画出该二项分布的密度函数曲线；调用 random 函数生成 100001 的二项分布（n=50，p=0.4）的随机数向量，然后作出频率直方图。解：（1）x=0:10;y=binopdf(x,50,0.4)y = 0.0000
9、 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0002 0.0005 0.0014（2） x=random(bino,50,0.4,10000,1); fp,xp=ecdf(x); ecdfhist(fp,xp,50);xlabel(二项分布（ n=50,p=0.4）随机数); ylabel(f(x);74. 使用 cftool 工具对下面的数据进行拟合，并画出拟合曲线。 （至少使用两种曲线进行拟合） 。Xdata 1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 1890Ydata 3
10、.9 5.3 7.2 9.6 12.9 17.1 23.1 31.4 38.6 50.2 62.9Xdata 1900 1910 1920 1930 1940 1950 1960 1970 1980 1990Ydata 76.0 92.0 105.7 122.8 131.7 150.7 179.0 205.0 226.5 248.7解：第一种曲线拟合x=1790,1800,1810,1820,1830,1840,1850,1860,1870,1880,1890,1900,1910,1920,1930,1940,1950,1960,1970,1980,1990;y=3.9,5.3,7.2,9.6
11、,12.9,17.1,23.1,31.4,38.6,50.2,62.9,76,92,105.7,122.8,131.7,150.7,179,205,226.5,248.7Cftool89第二种曲线拟合10115. （1）求方程 的零点；（2）用 fsolve 解下面的非线性方程组在（ 1,1,1）附近的根004sin2xyzex解：（1） fzero(t)(cos(t)2*exp(-2*t)-3*abs(t),0)ans = 0.2097解：（2）syms x edit myfun （在 mufun 中输入题目中的方程组） x0 = 1,1,1; options=optimset(Display, iter); x,fval,exitflag,output = fsolve(myfun,x0,options)计算结果x = 1.1424 -1.0000 0.0000fval = 1.0e-011 *0.0013-0.0001-0.1232exitflag = 1output = iterations: 9funcCount: 37algorithm: trust-region doglegfirstorderopt: 1.4079e-012)3)(cos)(2tettft12message: 1x76 char实验结果评价13成绩： 指导教师签字：
