【導讀】MATLAB會將運算結果直接存入一變量ans，代表MATLAB運算后的答案，并顯示其數(shù)值于屏幕上。（為簡便起見，在下述各例中，我們不再印出MATLAB的。不同，此提示符號常會消失不見，但這并不會影響到MATLAB的運算結果。由上例可知，MATLAB認識所有一般常用到的加。MATLAB同時也會自動進行內存的使用和回收，而不必像C語言，必須由使用者一一指定。這些功能使的MATLAB易學易用，使用者可專心致力于撰。寫程序，而不必被軟件枝節(jié)問題所干擾。若不想讓MATLAB每次都顯示運算結果，只需在表達式最后加上分號（；）即可，下表即為MATLAB常用的基本數(shù)學函數(shù)及三角函數(shù)：。y=[]%刪除第四個元素，文字均可視為程序的注解。MATLAB亦可取出向量的一個元素或一部份。例如已知inv是用來計算反矩陣，鍵入help. inv即可得知有關inv命令的用法。inverse，MATLAB即會列出所有和關鍵詞inverse相關的指令。后，即可用help進一步找出其用法。（lookfor事實上是對所有在搜尋路徑下的。M檔案進行關鍵詞對第一注解行的比對，詳見后敘。

　　

【正文】 注 解 title(39。Two Trigonometric Functions39。)。 % 圖 形 標題 legend(39。y = sin(x)39。,39。y = cos(x)39。)。 % 圖 形注解 grid on。 % 顯 示網(wǎng)格 線 我們 可用 subplot 來 同 時畫 出 數(shù)個 小 圖 形于同一 個 窗口之中： subplot(2,2,1)。 plot(x, sin(x))。 subplot(2,2,2)。 plot(x, cos(x))。 subplot(2,2,3)。 plot(x, sinh(x))。 第 23 頁 共 36 頁 subplot(2,2,4)。 plot(x, cosh(x))。 MATLAB 還 有其它各種二 維繪圖 函 數(shù) ，以適合不同的 應 用， 詳見下 表。 小整理：其它各種二 維繪圖 函 數(shù) bar 長條圖 errorbar 圖 形加上 誤 差范 圍 fplot 較 精確的函 數(shù)圖 形 polar 極坐 標圖 hist 累 計圖 rose 極坐 標 累 計圖 stairs 階 梯 圖 stem 針狀圖 fill 實 心 圖 feather 羽毛 圖 pass 羅盤圖 quiver 向量 場圖 以下 我們針對 每 個 函 數(shù)舉 例。 當數(shù) 據(jù) 點數(shù) 量不 多時 ， 長條圖 是很適合的表示方式： close all。 % 關閉 所有的 圖 形窗口 x=1:10。 y=rand(size(x))。 bar(x,y)。 第 24 頁 共 36 頁 如果已知 資 料 的誤 差量，就可用 errorbar 來 表示。下例以 單 位 標準差來 做 數(shù) 據(jù) 的誤 差量： x = linspace(0,2*pi,30)。 y = sin(x)。 e = std(y)*ones(size(x))。 errorbar(x,y,e) 對 于 變 化 劇 烈的函 數(shù) ，可用 fplot 來進 行 較 精確的 繪圖 ， 會對劇 烈 變 化 處進 行 較 密集的取 樣 ，如下例： fplot(39。sin(1/x)39。, [ ])。 % [ ]是 繪圖 范 圍 第 25 頁 共 36 頁 若要 產 生極坐 標圖 形，可用 polar： theta=linspace(0, 2*pi)。 r=cos(4*theta)。 polar(theta, r)。 對 于大量的 數(shù) 據(jù)， 我們 可用 hist 來顯 示 數(shù) 據(jù)的分 情 況 和 統(tǒng)計 特性。下面幾 個 命令可用 來驗證 randn 產 生的高斯 隨 機 數(shù) 分 ： x=randn(5000, 1)。 % 產 生 5000 個 ?=0， ?=1 的高斯 隨 機 數(shù) hist(x,20)。 % 20 代表 長條 的 個數(shù) 第 26 頁 共 36 頁 rose 和 hist 很接近，只不 過是將數(shù) 據(jù)大小 視為角 度， 數(shù) 據(jù) 個數(shù)視為 距離，并用極坐 標繪 制表示： x=randn(1000, 1)。 rose(x)。 stairs 可 畫 出 階 梯 圖 ： x=linspace(0,10,50)。 y=sin(x).*exp(x/3)。 stairs(x,y)。 第 27 頁 共 36 頁 stems 可 產 生 針狀圖 ，常被用 來繪 制 數(shù) 字 訊號 ： x=linspace(0,10,50)。 y=sin(x).*exp(x/3)。 stem(x,y)。 stairs 將資 料 點視為多邊 行 頂點 ，并 將 此 多邊 行涂上 顏 色： x=linspace(0,10,50)。 y=sin(x).*exp(x/3)。 fill(x,y,39。b39。)。 % 39。b39。為藍 色 第 28 頁 共 36 頁 feather 將每 一 個資 料 點視 復 數(shù) ，并以箭 號畫 出： theta=linspace(0, 2*pi, 20)。 z = cos(theta)+i*sin(theta)。 feather(z)。 pass 和 feather 很接近，只是每 個 箭 號 的起 點 都在 圓點 ： theta=linspace(0, 2*pi, 20)。 z = cos(theta)+i*sin(theta)。 pass(z)。 第 29 頁 共 36 頁 第 30 頁 共 36 頁 XYZ立體 繪圖 命令 在科 學目視 表示（ Scientific visualization）中，三度空 間 的立體 圖 是一 個 非常重要的技巧。本章 將 介 紹 MATLAB 基本 XYZ 三度空 間 的各 項繪圖 命令。 mesh和 plot 是三度空 間立 體 繪圖 的基本命令， mesh 可 畫 出立體網(wǎng) 狀圖 ， plot 則可畫 出三 維 曲面 圖 ， 兩 者 產 生的 圖 形都 會 依高度而有不同 顏 色。下列命令可 畫 出由函數(shù) 形成的立體網(wǎng) 狀圖 : x=linspace(2, 2, 25)。 % 在 x 軸 上取 25 點 y=linspace(2, 2, 25)。 % 在 y 軸 上取 25 點 [xx,yy]=meshgrid(x, y)。 % xx 和 yy 都是 21x21 的矩 陣 zz=xx.*exp(xx.^2yy.^2)。 % 計 算函 數(shù) 值， zz 也是 21x21 的矩 陣 mesh(xx, yy, zz)。 % 畫 出立體網(wǎng) 狀圖 surf 和 mesh 的用法 類 似： x=linspace(2, 2, 25)。 % 在 x 軸 上取 25 點 y=linspace(2, 2, 25)。 % 在 y 軸 上取 25 點 [xx,yy]=meshgrid(x, y)。 % xx 和 yy 都是 21x21 的矩 陣 zz=xx.*exp(xx.^2yy.^2)。 % 計 算函 數(shù) 值， zz 也是 21x21 的矩 陣 surf(xx, yy, zz)。 % 畫 出三 維 曲面 圖 第 31 頁 共 36 頁 為 了方便 測試 立體 繪圖 ， MATLAB 提供了一 個 peaks 函 數(shù) ，可 產 生一 個 凹凸有致的曲面，包含了三 個 局部極大 點 及三 個 局部極小 點 ，其方程式 為： 要 畫 出此函 數(shù) 的最快方法即是直 接鍵入 peaks： peaks z = 3*(1x).^2.*exp((x.^2) (y+1).^2) ... 10*(x/5 x.^3 y.^5).*exp(x.^2y.^2) ... 1/3*exp((x+1).^2 y.^2) 我們 亦可 對 peaks 函 數(shù)取點 ，再以各種不同方法 進 行 繪圖 。 meshz 可將 曲面加上 圍裙： 第 32 頁 共 36 頁 [x,y,z]=peaks。 meshz(x,y,z)。 axis([inf inf inf inf inf inf])。 waterfall 可在 x方向或 y 方向 產 生水流效果： [x,y,z]=peaks。 waterfall(x,y,z)。 axis([inf inf inf inf inf inf])。 下列命令 產 生在 y方向的水流效果： [x,y,z]=peaks。 waterfall(x39。,y39。,z39。)。 axis([inf inf inf inf inf inf])。 第 33 頁 共 36 頁 meshc 同 時畫 出網(wǎng) 狀圖 與等高 線： [x,y,z]=peaks。 meshc(x,y,z)。 axis([inf inf inf inf inf inf])。 surfc 同 時畫 出曲面 圖 與等高 線： [x,y,z]=peaks。 surfc(x,y,z)。 axis([inf inf inf inf inf inf])。 第 34 頁 共 36 頁 contour3 畫 出曲面在三度空 間 中的等高 線： contour3(peaks, 20)。 axis([inf inf inf inf inf inf])。 contour 畫 出曲面等高 線 在 XY平面的投影： contour(peaks, 20)。 第 35 頁 共 36 頁 plot3 可 畫 出三度空 間 中的 曲線： t=linspace(0,20*pi, 501)。 plot3(t.*sin(t), t.*cos(t), t)。 亦可同 時畫 出 兩條 三度空 間 中的 曲線： t=linspace(0, 10*pi, 501)。 plot3(t.*sin(t), t.*cos(t), t, t.*sin(t), t.*cos(t), t)。 第 36 頁 共 36 頁