matlabl實(shí)驗(yàn)數(shù)據(jù)分析

...wd......wd......wd...實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)課程名稱《金融MATLAB》開(kāi)課系部金融學(xué)院班級(jí)12金融1班學(xué)號(hào)20125633姓名宋歌指導(dǎo)教師洪振木2014年5月12日實(shí)驗(yàn)名稱MATLAB數(shù)據(jù)分析學(xué)院金融學(xué)院學(xué)號(hào)20125633姓名宋歌實(shí)驗(yàn)準(zhǔn)備實(shí)驗(yàn)?zāi)康膶W(xué)會(huì)MATLAB與Excel的數(shù)據(jù)交換；學(xué)會(huì)使用MATLAB進(jìn)展現(xiàn)金流分析；學(xué)會(huì)使用MATLAB計(jì)算概率密度、分布函數(shù)值，學(xué)會(huì)生成各種隨機(jī)數(shù)；掌握MATLAB中的CFTOOL數(shù)據(jù)擬合工具；學(xué)會(huì)如何使用MATLAB來(lái)求方程和方程組的數(shù)值解。實(shí)驗(yàn)設(shè)計(jì)方案1.利用MATLAB讀取Chapter3_code文件夾中StockPriceData.xls的數(shù)據(jù)，分別找出所有最高價(jià)中的最大值和所有最低價(jià)中的最小值，將最大值和最小值分別寫(xiě)入StockPriceData.xls文件的sheet2工作表中；并在同一圖像中畫(huà)出最高價(jià)曲線和最低價(jià)曲線。2.投資人貸款65W買(mǎi)房，還款期20年，每月還4000元，則貸款利率為多少

假設(shè)改為每月還4500元，貸款利率不變，則還貸期限為多長(zhǎng)〔提示：年金現(xiàn)金流計(jì)算〕3.求參數(shù)為n=50，p=0.4的二項(xiàng)分布在各點(diǎn)處的密度函數(shù)值與分布函數(shù)值，然后畫(huà)出該二項(xiàng)分布的密度函數(shù)曲線；調(diào)用random函數(shù)生成10000×1的二項(xiàng)分布〔n=50，p=0.4〕的隨機(jī)數(shù)向量，然后作出頻率直方圖。4.使用cftool工具對(duì)下面的數(shù)據(jù)進(jìn)展擬合，并畫(huà)出擬合曲線?！仓辽偈褂脙煞N曲線進(jìn)展擬合〕。Xdata17901800181018201830184018501860187018801890Ydata3.95.37.29.612.917.123.131.438.650.262.9Xdata1900191019201930194019501960197019801990Ydata76.092.0105.7122.8131.7150.7179.0205.0226.5248.75.〔1〕求方程的零點(diǎn)；〔2〕用fsolve解下面的非線性方程組在〔1,1,1〕附近的根1.利用MATLAB讀取Chapter3_code文件夾中StockPriceData.xls的數(shù)據(jù)，分別找出所有最高價(jià)中的最大值和所有最低價(jià)中的最小值，將最大值和最小值分別寫(xiě)入StockPriceData.xls文件的sheet2工作表中；并在同一圖像中畫(huà)出最高價(jià)曲線和最低價(jià)曲線。解：1〕讀取數(shù)據(jù)[data,textdate]=xlsread('C:\Users\Administrator\Desktop\Chapter3_code\StockPriceData.xls')2〕找出所有最高價(jià)中的最大值和所有最低價(jià)中的最小值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.在圖像中畫(huà)出最高價(jià)曲線和最低價(jià)曲線。>>textdate=[0:509]>>Y1=xlsread('StockPriceData.xls',1,'C2:C511')>>Y2=xlsread('StockPriceData.xls',1,'D2:D511')>>plot(textdate,Y1,'r',textdate,Y2,'g')實(shí)驗(yàn)結(jié)論2.投資人貸款65W買(mǎi)房，還款期20年，每月還4000元，則貸款利率為多少假設(shè)改為每月還4500元，貸款利率不變，則還貸期限為多長(zhǎng)〔提示：年金現(xiàn)金流計(jì)算〕解：第一問(wèn)是年金現(xiàn)金流計(jì)算，可以用annurate函數(shù)，NumPeriods=20*12Payment=4000PresentValue=650000FutureValue=0;Due=0結(jié)果如下rate=annurate(20*12,4000,650000,0,0)rate=0.0035年利率a=rate*12a=0.0418第二問(wèn)是計(jì)算年金周期函數(shù)，可以用annuterm函數(shù)計(jì)算Rate=0.0035Payment=4000PresentValue=650000FutureValue=0Due=0結(jié)果如下numperiods=annuterm(0.0035,-4500,650000,0,0)numperiods=201.5864運(yùn)行結(jié)果numperiods=201.5864〔月〕=16.7989〔年〕3.求參數(shù)為n=50，p=0.4的二項(xiàng)分布在各點(diǎn)處的密度函數(shù)值與分布函數(shù)值，然后畫(huà)出該二項(xiàng)分布的密度函數(shù)曲線；調(diào)用random函數(shù)生成10000×1的二項(xiàng)分布〔n=50，p=0.4〕的隨機(jī)數(shù)向量，然后作出頻率直方圖。解：〔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('二項(xiàng)分布〔n=50,p=0.4〕隨機(jī)數(shù)');>>ylabel('f(x)');4.使用cftool工具對(duì)下面的數(shù)據(jù)進(jìn)展擬合，并畫(huà)出擬合曲線?！仓辽偈褂脙煞N曲線進(jìn)展擬合〕。Xdata17901800181018201830184018501860187018801890Ydata3.95.37.29.612.917.123.131.438.650.262.9Xdata1900191019201930194019501960197019801990Ydata76.092.0105.7122.8131.7150.7179.0205.0226.5248.7解：第一種曲線擬合x(chóng)=[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〕求方程的零點(diǎn)；〔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)計(jì)算結(jié)果x=1.1424-1.00000.0000fval=1.0e-011*0.0013-0.0001-