matlabl实验数据分析

...wd......wd......wd...实验报告实验课程名称《金融MATLAB》开课系部金融学院班级12金融1班学号20125633姓名宋歌指导教师洪振木2014年5月12日实验名称MATLAB数据分析学院金融学院学号20125633姓名宋歌实验准备实验目的学会MATLAB与Excel的数据交换；学会使用MATLAB进展现金流分析；学会使用MATLAB计算概率密度、分布函数值，学会生成各种随机数；掌握MATLAB中的CFTOOL数据拟合工具；学会如何使用MATLAB来求方程和方程组的数值解。实验设计方案1.利用MATLAB读取Chapter3_code文件夹中StockPriceData.xls的数据，分别找出所有最高价中的最大值和所有最低价中的最小值，将最大值和最小值分别写入StockPriceData.xls文件的sheet2工作表中；并在同一图像中画出最高价曲线和最低价曲线。2.投资人贷款65W买房，还款期20年，每月还4000元，则贷款利率为多少

假设改为每月还4500元，贷款利率不变，则还贷期限为多长〔提示：年金现金流计算〕3.求参数为n=50，p=0.4的二项分布在各点处的密度函数值与分布函数值，然后画出该二项分布的密度函数曲线；调用random函数生成10000×1的二项分布〔n=50，p=0.4〕的随机数向量，然后作出频率直方图。4.使用cftool工具对下面的数据进展拟合，并画出拟合曲线。〔至少使用两种曲线进展拟合〕。Xdata17901800181018201830184018501860187018801890Ydata3.95.37.29.612.917.123.131.438.650.262.9Xdata1900191019201930194019501960197019801990Ydata76.092.0105.7122.8131.7150.7179.0205.0226.5248.75.〔1〕求方程的零点；〔2〕用fsolve解下面的非线性方程组在〔1,1,1〕附近的根1.利用MATLAB读取Chapter3_code文件夹中StockPriceData.xls的数据，分别找出所有最高价中的最大值和所有最低价中的最小值，将最大值和最小值分别写入StockPriceData.xls文件的sheet2工作表中；并在同一图像中画出最高价曲线和最低价曲线。解：1〕读取数据[data,textdate]=xlsread('C:\Users\Administrator\Desktop\Chapter3_code\StockPriceData.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:D511'))b=1.0005e+003>>[status2,message2]=xlswrite('StockPriceData.xls',b,'sheet2','A2')status2=1message2=message:''identifier:''3.在图像中画出最高价曲线和最低价曲线。>>textdate=[0:509]>>Y1=xlsread('StockPriceData.xls',1,'C2:C511')>>Y2=xlsread('StockPriceData.xls',1,'D2:D511')>>plot(textdate,Y1,'r',textdate,Y2,'g')实验结论2.投资人贷款65W买房，还款期20年，每月还4000元，则贷款利率为多少假设改为每月还4500元，贷款利率不变，则还贷期限为多长〔提示：年金现金流计算〕解：第一问是年金现金流计算，可以用annurate函数，NumPeriods=20*12Payment=4000PresentValue=650000FutureValue=0;Due=0结果如下rate=annurate(20*12,4000,650000,0,0)rate=0.0035年利率a=rate*12a=0.0418第二问是计算年金周期函数，可以用annuterm函数计算Rate=0.0035Payment=4000PresentValue=650000FutureValue=0Due=0结果如下numperiods=annuterm(0.0035,-4500,650000,0,0)numperiods=201.5864运行结果numperiods=201.5864〔月〕=16.7989〔年〕3.求参数为n=50，p=0.4的二项分布在各点处的密度函数值与分布函数值，然后画出该二项分布的密度函数曲线；调用random函数生成10000×1的二项分布〔n=50，p=0.4〕的随机数向量，然后作出频率直方图。解：〔1〕x=0:10;y=binopdf(x,50,0.4)y=0.00000.00000.00000.00000.00000.00000.00000.00000.00020.00050.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)');4.使用cftool工具对下面的数据进展拟合，并画出拟合曲线。〔至少使用两种曲线进展拟合〕。Xdata17901800181018201830184018501860187018801890Ydata3.95.37.29.612.917.123.131.438.650.262.9Xdata1900191019201930194019501960197019801990Ydata76.092.0105.7122.8131.7150.7179.0205.0226.5248.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,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.7]Cftool第二种曲线拟合5.〔1〕求方程的零点；〔2〕用fsolve解下面的非线性方程组在〔1,1,1〕附近的根解：〔1〕>>fzero(@(t)(cos(t))^2*exp(-2*t)-3*abs(t),0)ans=0.2097解：〔2〕symsx>>editmyfun〔在mufun中输入题目中的方程组〕>>x0=[1,1,1];>>options=optimset('Display','iter');>>[x,fval,exitflag,output]=fsolve(@myfun,x0,options)计算结果x=1.1424-1.00000.0000fval=1.0e-011*0.0013-0.0001-