【正文】 ? [例 3]函數(shù)文件 ? ? function y=mean(x) ? %文件的第一條語句定義了函數(shù)名、輸入變元及輸出變元。 ? echo on(off) 在執(zhí)行 M文件每行程序前先顯示其內(nèi)容。只能被自身和 1除凈的數(shù)。 ? %計算 fibonnaci數(shù)的 M文件 ? clear,close all ? N=input(39。 ? 要注意流程控制語句的縮進及與 end的對應(yīng)關(guān)系。 MATLAB軟件的大部分功能是來自其建立的函數(shù)集。1:9]’) caxis([cmin cmax]) surfl(peaks,[30,0]) 彩色條分為 64個等級，顏色排序的默認值為 HSV，用戶可以自行定義為 hot， cool， gray， copper， pink， jet，prism等， m=colormap將返回一個 64 3階的矩陣 65個分割點，產(chǎn)生 64根彩條 hot(8)將偽彩色板設(shè)置為 hot色，改為 8份。 for j=1:16 plot(fft(eye(j+16)))。surfc(z1),shading flat39。 meshc(z),shading flat39。 z1=sin(R)./R。 z=sin(R)./R。 y=sin(z)。 v=axis axis(39。 x=sin(z)。即相當于 plot(real(z),imag(z))。復(fù)數(shù)繪圖 plot(z)39。39。) subplot (2,2,3), bar(t,y)。 figure 創(chuàng)建圖形窗口 shg 顯示圖形 圖形 窗的 控制 gcf 獲取當前圖形窗的句柄 refresh 刷新圖形 clf 清除當前圖形窗 close 關(guān)閉圖形窗 axes 在任意位置創(chuàng)建坐標系 ishold 保持當前圖形狀態(tài)為真 軸系 的 控制 gca 獲取當前坐標系的句柄 box 形成軸系方向 cla 清除當前坐標系 line 創(chuàng)建直線 surface 創(chuàng)建曲面 圖形 對象 patch 創(chuàng)建圖形填充塊 light 創(chuàng)建照明 image 創(chuàng)建圖象 通用圖形函數(shù)（ graphics） (h) set 設(shè)置對象特性 gcbo 獲得回叫對象的曲柄 圖形 get 獲得對象特性 gcbf 獲得回叫圖形的曲柄 句柄 reset 復(fù)位對象特性 drawnow 直接等待圖形事件 操作 delet 刪除對象 findobj 尋找具有特定值的對象 gco 獲得當前對象的句柄 copyobj 為圖形對象及其子項作硬拷貝 工具 closereq 請求關(guān)閉圖形窗 ishandle 是圖形句柄時為真 newplot 說明 Nextplot的 M文件 ginput 從鼠標作圖形輸入 uiputfile 給出存儲文件的對話框 雜項 graymon 設(shè)定圖形窗灰度監(jiān)視器 uigetfile 給出詢問文件名的對話框 rbbox 涂抹塊 whitebg 設(shè)定圖形窗背景色 rotate 圍繞指定方向旋轉(zhuǎn)對象 zoom 二維圖形的放大和縮小 terminal 設(shè)定圖形終端類型 warndlg 警告對話框 （ 2）其他二維繪圖命令 stem stairs bar errorbar hist fill(t,y,39。y339。 plotyy(t,y39。)。)。不便于用戶自行設(shè)定線型和顏色。+r39。 加標題， 坐標軸說明 加坐標網(wǎng)格線 t=0::4*pi y=exp(*t).*sin(t) plot(t,y) y1=exp(*t).*sin(t+1) plot(t,y1,39。, end ? n = 8 ? a =偶 ? 判斷輸入數(shù) n的奇、偶、空的程序 ? mod（ x， m） x整除 m取正余數(shù)， ? rem（ a， b） a整除 b，求余數(shù) ? n= 負數(shù) ? n = 5 ? a =奇 ? n=8 ? a =偶 ? 基本繪圖方法 ? 1．直角坐標中的兩維曲線 ? （ 1） plot(y)輸入一個數(shù)組的情況 y=5*(rand(1,10) ) plot(y) title (39。, end ? n = 5 ? a =奇 ? n=input( 39。 ? n=input( 39。 N=1+（終值 初值） /增量 ? 用 for語句求三角函數(shù)表 ? for x=0: : pi/4 disp([x, sin(x), cos(x), tan(x)]), end ? 運行結(jié)果 ? x sin(x) cos(x) tan(x) ? 0 0 1 0 ? 1/10 839/8404 1195/1201 1499/14940 ? 1/5 209/1052 295/301 374/1845 ? 3/10 409/1384 1647/1724 275/889 ? 2/5 368/945 2882/3129 1777/4203 ? 1/2 501/1045 1699/1936 820/1501 ? 3/5 1153/2042 430/521 979/1431 ? 7/10 947/1470 992/1297 486/577 列出構(gòu)成 Hilbert矩陣的程序 format rat顯示形式是分數(shù)近似 n=input(39。 X不斷增大，直到無法表示它的值，只能用 inf表示為止。odd39。, end ? n = 7 a =odd。 與 或 非 異或 G整除 3，求余數(shù) p=(rem(G, 3)==0) u=p|~p all(p) all(u) any(p) 可以按行、按列判斷一群元素的邏輯值。顯示 x sin(x) cos(x) tan(x)39。 元素群的運算的兩個矩陣必須是同階的。 4 9 2 2] B1=fliplr(A) B2=flipud(A) B3=reshape(A,2,6) 提取矩陣中某些特殊結(jié)構(gòu)的元素， 組成新的矩陣，改變矩陣結(jié)構(gòu)。除此以外，其他 MATLAB函數(shù)都是對矩陣中的元素分別進行，英文直譯為數(shù)組運算（ Array Operations），譯為“ 元素群運算 ” S=[1 2。 3 6 0] ? 運算： a*b d\a ? a*b ? ??? Error using == * ? Inner matrix dimensions must agree. ? d\a ? ??? Error using == \ ? Matrix dimensions must agree. ? a39。 3 5 4。 （ 2）矩陣的除法及線性方程組的解 a =1 2 3 4 5 6 7 8 9 AV=I V=A1 V=inv(a) inv(a)*a V = +016 * n n階方陣 A和同階的方陣 V相乘，得出 n階單位矩陣 I。 x*y’ ans = 2 ans = 2 0 2 y‘*x X右乘 y’。 [n,m]=size(fb2) x=[1 0 1]。f4,f2] fb2=[fb1。 6,8]*i f=sqrt(1+2i) f*f 復(fù)數(shù)的虛數(shù)部分用 i或 j表示，如曾用過 i, j 作變量，用 clear i,j 復(fù)數(shù)矩陣有兩種賦值方法： ①將其元素逐個賦予復(fù)數(shù)； ②將其實部和虛部矩陣分別賦值。語句結(jié)尾用回車或逗號，會顯示結(jié)果，如果不想顯示結(jié)果，用分號。標識符可以是英文字母、數(shù)字和下劃線等符號。這種文件簡單、短小、高效，并且便于調(diào)試。 ? MATLAB語言有如下優(yōu)點： ? ? MATLAB的基本數(shù)據(jù)單元是既不需要指定維數(shù)、也不需要說明數(shù)據(jù)類型的矩陣，而且數(shù)學(xué)表達式和運算規(guī)則與通常的習(xí)慣相同。因此，在 MATLAB環(huán)境下，數(shù)組的操作與數(shù)的操作一樣簡單。 ? ? MATLAB具有二維和三維繪圖功能，使用方法十分簡便。標識符第 1個字符必須是英文字母， MATLAB對大、小寫敏感。 元素用（）中的數(shù)字（下標）來注明，一維用一個下標，二維用兩個下標，逗號分開。 w=z’ （共軛轉(zhuǎn)置） u=conj(z) （共軛） v=conj(z)’ （轉(zhuǎn)置） Z’復(fù)數(shù)矩陣共軛轉(zhuǎn)置：行列互換，各元素的虛部反號。f5] 全 1矩陣 全 0矩陣 魔方矩陣：元素由 1到 nn的自然數(shù)組成，每行、每列及兩對角線上的元素之和均等于 (n3+n)/2。 y=x1 y = 2 1 0 語句 size檢查矩陣階數(shù)，兩矩陣相加，階數(shù)必須相同。 1 0 1 0 0 0 eye(3)*a 左、右乘結(jié)果不同，只有單位矩陣例外。 I為 eye(n)。 7 8 9] ? x=[x1,x2,x3] ? b=[2。*b ans = ? 6 16 20 ? 9 23 25 ? 12 30 30 ? a*b39。 3 4] D=[1 4 7。 fliplr矩陣左右翻轉(zhuǎn) flipud矩陣上下翻轉(zhuǎn) reshape階數(shù)重組（元素總數(shù)不變） B4=rot90(A) B5=diag(A) B6=tril(A) B7=triu(A) B8=A(: )39。（標量會自動擴展為同階矩陣參與運算） x=[1, 2, 3] y=[4, 5, 6] z=x.*y z=x.\y z=x.^y z=x.^2 z=2.^[x y] z =4 10 18 x*y不能成立 z = 元素群沒有左除右除之分 z =1 32 729 x^y 能成立嗎？ z =1 4 9 x^2能成立嗎？ z =2 4 8 16 32 64 2^[x y] 能成立嗎？ d=[1 4 7。) ? disp([x, sin(x) cos(x) tan(x)]) ? 顯示 x sin(x) cos(x) tan(x) ? 0 0 0 ? ? ? ? ? ? ? 邏輯判斷及流程控制 1. 關(guān)系運算 a= 2+2==4 a=(2+2==4) a=(34) a=(43) a=(3=4) a=(4=3) a=(43) a=(34) 等于， a = 1 a = 1 小于， a = 1 a = 0 小于等于， a = 1 a = 0 大于， a =1 a =0 a=(4=3) a=(3=4) a=(3~=4) A=magic(6) rem(A,3) p=(rem(A, 3)==0) lp=find (p)39。 兩個對元素群運算的函數(shù)： 列中有一個元素為 0，即為