第3章matlab與excel數(shù)據(jù)連結(jié)及繪圖

金融數(shù)值計算與應(yīng)用軟件

--MATLAB王輝中央財經(jīng)大學金融學院E-mail：Tel：第3章MATLAB與Excel數(shù)據(jù)連接及繪圖3.1MATLAB與Excel數(shù)據(jù)連接3.1.1MATLAB和Excel的接口安裝：.ExcelLink的安裝和和設(shè)置ExcelLink的設(shè)置

啟動Excel，選擇菜單“工具”項下的“加載宏”項，彈出如下對話框若該項不存在，則通過瀏覽目錄，在目錄

MATLAB/toolbox/exlink下找到

excllink.xla文件在Excel中多了一個ExcelLink工具條，如下圖示:3.1.2.ExcelLink管理函數(shù)

手動啟動matlab

要不想在每次啟動時都啟動matlab，則可以在電子表格中或函數(shù)輸入框中輸入以下命令=MLAutoStart(“no”)手動啟動matlab(1)在excel電子表格中輸入命令=MLOpen()(2)利用excel的宏工具，在宏對話框中輸入Matlabinit，然后單擊執(zhí)行按鈕。

(3)按左上按鈕"startmatlab"。1、ExcelLink連接管理函數(shù)

(1)Matlabinit

該函數(shù)只能在宏子例程中使用。初始化ExcelLink和啟動Matlab進程。只有在MLAutoStart函數(shù)中使用“no”參數(shù)，才需要手動使用Matlabinit來初始化ExcelLink和啟動Matlab進程，

如果使用參數(shù)“yes”，則Matlabinit是自動執(zhí)行的。

使用語法：Matlabinit

(2)MLAutoStart設(shè)置自動啟動Matlab和ExcelLink。

在工作表中的使用語法：

MLAutoStart("yes")

MLAutoStart("no")

使用“yes”參數(shù)，則當Excel啟動時，自動啟動Matlab和ExcelLink；如果使用參數(shù)“no”，則當Excel啟動時，不啟動Matlab和ExcelLink。如果在此之前它們已經(jīng)啟動，則無任何影響。

(3)MLClose終止Matlab進程并刪除Matlab工作空間的所有變量。并通知Excel，Matlab不再運行。

在工作表中的使用語法：

MLClose（）

在宏中的使用語法：

MLClose

(4)MLOpen—啟動Matlab進程如果Matlab進程已經(jīng)啟動，則MLOpen函數(shù)不進行任何操作。在使用MLClose關(guān)閉Matlab進程后使用MLOpen來重新啟動Matlab。

在工作表中的使用語法：

MLOpen（）

在宏中的使用語法：

MLOpen

(3)MLClose終止Matlab進程并刪除Matlab工作空間的所有變量。并通知Excel，Matlab不再運行。

在工作表中的使用語法：

MLClose（）

在宏中的使用語法：

MLClose

Excel時間與Matlab的轉(zhuǎn)換EXCEL的默認日期數(shù)開始與1900年1月1日，而MATLAB的日期數(shù)開始于0000年1月1日。如果用戶使用MATLAB中的日期數(shù)計算,常數(shù)693960起關(guān)鍵作用:EXCEL日期數(shù)加上它進入MATLAB,或從MATLAB日期數(shù)中減去它進入EXCEL

2、ExcelLink數(shù)據(jù)管理函數(shù)

(1)matlabfcn

根據(jù)給定的Excel數(shù)據(jù)執(zhí)行Matlab命令。

在工作表中使用時的語法：

matlabfcn(command，inputs)=matlabfcn(“sum”,B3:B10)；把從B3到B10的單元格中數(shù)據(jù)相加例：在Matlab中生成12維均勻隨機向量aa，利用matlabfcn在Excel中求其分量之和。(2)Matlabsub根據(jù)給定的Excel數(shù)據(jù)執(zhí)行Matlab命令，并將結(jié)果返回到指定的單元格中。

在工作表中的使用語法：

matlabsub（command，edat，inputs）

例如：=matlabsub（”sum”，”A1”，B1:B10）；把工作表中從單元格B1到B10的數(shù)據(jù)相加，并將結(jié)果返回到單元格A1中。

注意：edat指定的位置不能包含matlabsub所在的位置。

(3)MLAppendMatrix將Excel工作表中的數(shù)據(jù)追加到Matlab中指定的矩陣中如果該矩陣不存在，則創(chuàng)建矩陣。

在工作表中使用的語法：

MLAppendMatrix(var_name，mdat)注意要追加的數(shù)據(jù)維數(shù)要和原矩陣中的維數(shù)相匹配，否則出錯

假設(shè)matlab中，矩陣a是個2行3列的矩陣，如下圖:

將Excel表中A1：A2中的數(shù)據(jù)追加到矩陣后，如下圖示成為矩陣的第四列，如下圖示

(4)MLDeleteMatrix

刪除Matlab空間中指定的矩陣

在工作表中使用的語法：

MLDeleteMatrix（var_name）；Var_name，是要刪除的矩陣名，如果矩陣名在引號內(nèi)“var_name”的形式，則直接指定變量名，如果var_name不用引號引起來，則通過計算獲得實際矩陣的名字。

例如，單元格B1中的內(nèi)容為a，MLDeleteMatrix（B1），則表示刪除Matlab工作空間的矩陣a，它等價于MLDeleteMatrix（”a”）4.xlswrite與xlsread

MATLAB提供的函數(shù)xlswrite，具有將MATLAB中的數(shù)據(jù)寫入EXCEL的功能。MATLAB提供的函數(shù)xlsread，具有將EXCEL中的數(shù)據(jù)讀入MATLAB的功能。xlswrite的調(diào)用方式

xlswrite('filename',M)xlswrite('filename',M,sheet)xlswrite('filename',M,'range')xlswrite('filename',M,sheet,'range')xlswrite的輸入?yún)?shù)filename%EXCEL文件名M%MATLAB工作空間中的大小為m×n的數(shù)組，該數(shù)組可以是字符型的，數(shù)值型的也可以是單元數(shù)組，其中m<65536且n<256sheet%EXCEL中的工作簿range%EXCEL工作簿中的數(shù)據(jù)區(qū)域例如：將一個向量寫入EXCEL文件v.xls，在MATLAB中執(zhí)行如下命令：xlswrite('v',[15-83.90])例3.1.1把元胞數(shù)組d={'Time','Temp';1298;1399;1497}寫入到EXCEL文件m.xls的一個規(guī)定的工作簿當中。答案見ex3_1_1.mxlsread的調(diào)用方式

num=xlsread(filename)num=xlsread(filename,-1)num=xlsread(filename,sheet)num=xlsread(filename,'range')num=xlsread(filename,sheet,'range')[num,txt]=xlsread(filename,sheet,'range')xlsread的輸入與輸出參數(shù)輸入?yún)?shù)：

filename%EXCEL文件名

sheet%EXCEL中的工作簿名

range%EXCEL工作簿中的數(shù)據(jù)區(qū)域輸出參數(shù)：

num%讀入Excel文件數(shù)據(jù)

txt%保存文本內(nèi)容

例3.1.2將m.xls復制到matlab當前工作目錄下，并在matlab命令窗口輸入以下命令：答案見ex3_1_2.m3.2MATLAB繪圖

3.2.1二維數(shù)據(jù)曲線圖1、繪制單根二維曲線plot函數(shù)最簡單的調(diào)用格式是只包含一個輸入?yún)?shù)：plot(x)在這種情況下，當x是實向量時，以該向量元素的下標為橫坐標，元素值為縱坐標畫出一條連續(xù)曲線，這實際上是繪制折線圖。plot函數(shù)的基本調(diào)用格式為：

plot(x,y)

其中x和y為長度相同的向量，分別用于存儲x坐標和y坐標數(shù)據(jù)。二維繪圖的例子例3.2.1在0≤x≤2

區(qū)間內(nèi)，繪制曲線

y=2e-0.5xcos(4πx)答案見ex3_2_1.m例3.2.2假如你的帳戶中有一筆10年的$1000存款，年利率為r,那么10年后你的錢數(shù)與利率的關(guān)系是,畫出該函數(shù)的圖像。答案見ex3_2_2.m

1、當x,y是同維矩陣時，則以x,y對應(yīng)列元素為橫、縱坐標分別繪制曲線，曲線條數(shù)等于矩陣列數(shù)。

2、對只包含一個輸入?yún)?shù)的plot函數(shù)，當輸入?yún)?shù)是實矩陣時，按列繪制每列元素值相對其下標的曲線，曲線條數(shù)等于矩陣的列數(shù)。3、當輸入?yún)?shù)是復數(shù)矩陣時，則按列分別以元素實部和虛部為橫、縱坐標繪制多條曲線。2、繪制多根二維曲線plot函數(shù)的輸入?yún)?shù)是矩陣形式例3.2.3畫出S&P從2008年1月到2008年10.20開盤價、最高價、最低價和收盤價的時間圖。

(1)橫坐標是觀測個數(shù)，縱坐標是價格，將四種價格畫在一幅圖中（2）橫坐標是時間，縱坐標是價格，將四種價格畫在一幅圖中（3）橫坐標是時間，縱坐標是價格，畫出其收盤價的時間圖答案見ex3_2_3.m

調(diào)用格式為：plot(x1,y1,x2,y2,…,xn,yn)(1)當輸入?yún)?shù)都為向量時，x1和y1，x2和y2，…，xn和yn分別組成一組向量對，每一組向量對的長度可以不同。(2)當輸入?yún)?shù)有矩陣形式時，配對的x,y按對應(yīng)列元素為橫、縱坐標分別繪制曲線，曲線條數(shù)等于矩陣的列數(shù)。含多個輸入?yún)?shù)的plot函數(shù)例3.2.4分析下列程序繪制的曲線。x1=linspace(0,2*pi,100);x2=linspace(0,3*pi,100);x3=linspace(0,4*pi,100);y1=sin(x1);y2=1+sin(x2);y3=2+sin(x3);x=[x1;x2;x3]';y=[y1;y2;y3]';plot(x,y,x1,y1-1)

具有兩個縱坐標標度的圖形

調(diào)用格式為：

plotyy(x1,y1,x2,y2)其中x1,y1對應(yīng)一條曲線，x2,y2對應(yīng)另一條曲線。橫坐標的標度相同，縱坐標有兩個，左縱坐標用于x1,y1數(shù)據(jù)對，右縱坐標用于x2,y2數(shù)據(jù)對。

例3.2.5用不同標度在同一坐標內(nèi)繪制曲線y1=0.2e-0.5xcos(4πx)和y2=2e-0.5xcos(πx)。答案見ex3_2_5.m例3.2.6畫出S&P從2008年1月到2008年10.20收盤價及其對數(shù)收益率的時間圖。

答案見ex3_2_6.m圖形保持holdon/off命令控制是保持原有圖形還是刷新原有圖形，不帶參數(shù)的hold命令在兩種狀態(tài)之間進行切換。例3.2.7采用圖形保持，在同一坐標內(nèi)繪制曲線y1=0.2e-0.5xcos(4πx)和y2=2e-0.5xcos(πx)。

答案見ex3_2_7.m設(shè)置曲線樣式

plot(X1,Y1,LineSpec1,X2,Y2,LineSpec2…)

將按順序分別畫出由三參數(shù)定Xi,Yi,LineSpeci的線條。其中參數(shù)LineSpeci指明了線條的類型，標記符號，和畫線用的顏色。

例3.2.8采用不同顏色、不同線形以及不同標記在同一圖形中畫出y=cos(x)，y=sin(x-

/2)exp{x/100}和y=sin(x-

)的圖像答案見ex3_2_8.m例3.2.9分析如下程序plot(t,sin(2*t),'-mo','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor',[.491.63],'MarkerSize',12)在原來程序的基礎(chǔ)上plot(mydate,sp08(:,3),'y-*',mydate,sp08(:,3),'m-.')例3.2.10畫出S&P從2008年1月到2008年10.20最高價、和最低價的時間圖，不同價格以不同的線形、顏色和標記畫出。圖形標注

title(圖形名稱)xlabel(x軸說明)ylabel(y軸說明)text(x,y,圖形說明)legend(圖例1,圖例2,…)例3.2.11用不同顏色及線形畫出HSI指數(shù)從1990年1月2日到1992年12月31日的收盤價和對數(shù)收益率的1萬倍的時間圖。圖形名稱：HIS收盤價和對數(shù)收益率X軸標注：觀測Y軸標注：HIS收盤價和對數(shù)收益率的1萬倍圖形說明：時間是1990.1.2-1992.12.31加圖例答案見ex3_2_11.m坐標控制

axis([xminxmaxyminymaxzminzmax])axisequal：縱、橫坐標軸采用等長刻度。axissquare：產(chǎn)生正方形坐標系(缺省為矩形)。axisauto：使用缺省設(shè)置。axisoff：取消坐標軸。axison：顯示坐標軸。網(wǎng)格線和邊框給坐標加網(wǎng)格線用grid命令來控制

gridon/off命令控制是畫還是不畫網(wǎng)格線，不帶參數(shù)的grid命令在兩種狀態(tài)之間進行切換。給坐標加邊框用box命令來控制

boxon/off命令控制是加還是不加邊框線，不帶參數(shù)的box命令在兩種狀態(tài)之間進行切換。例3.2.12在同一坐標中，可以繪制3個同心圓，并加坐標控制。答案見ex3_2_12.m圖形的可視化編輯MATLAB在圖形窗口中提供了可視化的圖形編輯工具，利用圖形窗口菜單欄或工具欄中的有關(guān)命令可以完成對窗口中各種圖形對象的編輯處理。在圖形窗口上有一個菜單欄和工具欄。菜單欄包含F(xiàn)ile、Edit、View、Insert、Tools、Window和Help共7個菜單項，工具欄包含11個命令按鈕。利用可視化的圖形編輯工具處理前面的例子圖形窗口的分割

subplot(m,n,p)該函數(shù)將當前圖形窗口分成m×n個繪圖區(qū)，即每行n個，共m行，區(qū)號按行優(yōu)先編號，且選定第p個區(qū)為當前活動區(qū)。在每一個繪圖區(qū)允許以不同的坐標系單獨繪制圖形。期權(quán)的利潤函數(shù)Call多頭：L=-Cer(T-t)+max(ST-X,0)Put多頭：L=-Per(T-t)+max(X-ST,0)例3.2.13假設(shè)投資者1和投資者2在3月份分別購買了一份12月份到期的歐式看漲期權(quán)和看跌期權(quán)，執(zhí)行價格都是X=100元，看漲期權(quán)價格C=8元，看跌期權(quán)價格為6元，無風險利率為3%，試給出該兩位投資者及其交易對手的利潤函數(shù)并用四幅子圖畫出利潤函數(shù)圖。答案見ex3_2_13.m3.2.2三維圖形

三維曲線

plot3(x1,y1,z1,選項1,x2,y2,z2,選項2,…,xn,yn,zn,選項n)其中每一組x,y,z組成一組曲線的坐標參數(shù)，選項的定義和plot函數(shù)相同。當x,y,z是同維向量時，則x,y,z對應(yīng)元素構(gòu)成一條三維曲線。當x,y,z是同維矩陣時，則以x,y,z對應(yīng)列元素繪制三維曲線，曲線條數(shù)等于矩陣列數(shù)。例3.2.14繪制三維曲線。程序如下：t=0:pi/100:20*pi;x=sin(t);y=cos(t);z=t.*sin(t).*cos(t);plot3(x,y,z);title('Linein3-DSpace');xlabel('X');ylabel('Y');zlabel('Z');gridon;

三維曲面

產(chǎn)生三維數(shù)據(jù)

在MATLAB中，利用meshgrid函數(shù)產(chǎn)生平面區(qū)域內(nèi)的網(wǎng)格坐標矩陣。其格式為：

x=a:d1:b