【正文】
。系統(tǒng)會暫停，等待用戶鍵入命令。會給調(diào)試帶來很大困難，解決方法： ? （ 1）把某些分號改為逗號，使中間結(jié)果能顯示在屏幕上，作為查錯依據(jù)。 ? 程序調(diào)試 ? 主程序不需要專門的調(diào)試命令，需要用調(diào)試命令的主要是函數(shù)程序。 nargout表示輸出變元數(shù)目的變量。 end %d2為 pi時的設(shè)置 y=(10).^[d1+(0:n2)*(d2d1)/(n1), d2]。 function y=logspace (d1, d2, n) % LOGSPACE 對數(shù)均分?jǐn)?shù)組 % LOGSPACE（ d1, d2）在 10d1與 10d2之間生成長度為 50的對數(shù)均分?jǐn)?shù)組 % 如果 d2為 pi，則這些點(diǎn)在 10d1和 pi之間 % LOGSPACE（ d1, d2， n）生成的數(shù)組長度為 n， n的默認(rèn)值為 50 if nargin ==2 n=50。假如它是一個已知向量或矩陣 Z，可以寫成 mean（ Z），該變量 Z通過變元替換傳遞給 mean函數(shù)后，在子程序內(nèi)，它就變成了局部變量 x。 If m==1 M=n。如果在執(zhí)行該文件前，工作空間中已有同名的變量，系統(tǒng)會把兩者看作各自無關(guān)的變量。對于向量， mean(x)返回該向量 x中各元素的平均值 % 對于矩陣， mean(x)是一個包含各列元素平均值的行向量 [m, n]=size (x)。（這條語句可以區(qū)分程序文件和函數(shù)文件）輸入變元、輸出變元可以有若干個，必須在第一條語句中列出。 （ 3）除非用 global聲明，程序中的變量均為局部變量，不保存在工作空間中。與主程序文件的主要區(qū)別有 3點(diǎn)： （ 1）由 function起頭，后面跟的函數(shù)名，函數(shù)名必須與文件名相同。menu 用來產(chǎn)生人際交互的備選擇菜單的命令。^c (controlc) 強(qiáng)行停止程序運(yùn)行的命令。數(shù)據(jù)輸入后，程序繼續(xù)運(yùn)行。 end, x N = 44 x =2 3 5 7 11 13 17 19 23 29 31 37 41 43 人機(jī)交互命令： x~=u)。), x=2:N。 %求素數(shù) (prime number)的程序 clear, close all N=input(39。 [例 2]求素數(shù)的程序。 i=i+1。 i=1。輸入最大數(shù)值范圍 N=39。 [例 1]列出一個求 fibonnaci數(shù)的程序，它是一個數(shù)列，從 [1， 1]開始，由數(shù)列的最后兩個元素之和生成新的元素，依次遞推。 ? （ 3） M文件的文件名、路徑名不能用中文，要按MATLAB的標(biāo)識符編制，因?yàn)?M文件也就是MATLAB的調(diào)用命令。不同行以分號隔開。選項(xiàng)可以自動對程序進(jìn)行縮進(jìn)排版。 ? 全局變量語句： global 變量名 1 變量名 2 …… ? 程序必須用半角英文字母和符號編制（包括標(biāo)點(diǎn)符號），只有 %號后面的注釋可以用中文。 （ 2）程序的主體 ? 全局變量：在子程序中和主程序中共用的變量。 MATLAB不執(zhí)行 %號后面的任何內(nèi)容。 ? 1．主程序文件 ? 主程序文件格式特征： ? （ 1）用 clear 、 close all等語句開始，清除原有的變量和圖形。 ? M文件分為兩種：一種是主程序，為用戶解決特定的問題編制的；一種是子程序，函數(shù)文件，必須由其他 M文件來調(diào)用，函數(shù)文件可以遞歸調(diào)用（自己調(diào)用自己）。 ()可以用任何編輯器來建立，可直接閱讀。方向照射時的 peak曲面。 人工設(shè)定的命令 caxis([cmin cmax]) 彩色條中的最小值、最大值 光源在方位角 30186。 1:65]’) hot(8) pcolor([1:9。 z=x.*exp(x.^2y.^2)。 end movie(m,30) eye單位矩陣 %因?yàn)楫a(chǎn)生的圖形是圓形，將坐標(biāo)設(shè)成相等比例 %為變量 M預(yù)留 16幅圖的存儲空間 %做 16次循環(huán) %畫圖 %依次存入 M中 16個畫面放在矩陣 M中，每幅 16858*8=134864個字節(jié)， 運(yùn)行 movie(m,30)，將 M中圖形播放 30遍 。 for j=1:16 plot(fft(eye(j+16)))。 用 movie命令把它們連起來重放。surfc(z1), view(20,0)39。),shading flat, %colormp(gray) subplot(2,2,4), surfc(z1),view(20,0)。 surfc(z1),pause title(39。),shading interp subplot(2,2,3), R=abs(X)+abs(Y)+eps。 mesh(z),pause title(39。),shading flat Subplot(2,2,2), R=sqrt(X.^2+Y.*Y)+eps。 meshc(z), pause title(39。 rotate3d contour coutour3 可用鼠標(biāo)拖動立體圖作空間連續(xù)轉(zhuǎn)動 把曲面的等高線投影在 XY平面 在三維立體圖中畫出等高線 subplot(2,2,1), R=sqrt(X.^2+Y.*Y)。 abs(X)+abs(Y)稱為一范數(shù) （ 3）其他三維繪圖命令 view(20, 0) view(37,30) view（方位角， 俯仰角）可以變換立體圖的視角 默認(rèn)值 shading flat shading interp shading faceted 可將曲面上的小格平滑掉，使曲面成為光滑表面。 figure(3), surf(z1) 在 R=0（原點(diǎn)）處出現(xiàn) 0/0運(yùn)算，得 NaN結(jié)果。 figure(2),mesh(z) R=abs(X)+abs(Y)+eps。 R=sqrt(X.*X+Y.*Y)+eps。 X、 Y方向各有 33個樣本點(diǎn)， size多維矩陣的各維長度。 z=sin(R)./R。 Y=y*ones(size(x))。 y=x39。r39。 y=sin(z)。) z=0::4*pi。) v = 0 1 0 1 5. 三維曲線和曲面 （ 1）空間曲線繪制 plot3 plot3(x,y,z,39。) axis(39。 控制圖形的縱橫比，使屏幕上 x與 y的比例尺相同 . v=axis axis(39。) 用戶自行設(shè)定坐標(biāo)比例，選擇圖形邊界范圍。 subplot(1,2,1),plot(x,y) subplot(1,2,2), plot(x,y) axis(39。 x=sin(z)。square39。equal39。 要在復(fù)平面內(nèi)繪多條曲線，必須用 hold命令，或?qū)⒍鄺l曲線的實(shí)部和虛部明確地寫出，作為函數(shù)的輸入變元，即Plot(real(z1),imag(z1),real(z2),imag(z2)) 子圖 1畫出了復(fù)數(shù)圖形； 子圖 2只畫出了 z的實(shí)部隨 t 的變化規(guī)律， 子圖 3是用極坐標(biāo)繪制的復(fù)數(shù)曲線； 子圖 4說明了半對數(shù)坐標(biāo)繪圖的結(jié)果。即相當(dāng)于 plot(real(z),imag(z))。semilogx(t,z)39。polar(angle(z),abs(z))39。復(fù)數(shù)繪圖 plot(t,z)39。復(fù)數(shù)繪圖 plot(z)39。 極坐標(biāo)繪圖，角度 theta為一個坐標(biāo)，單位是弧度，另一坐標(biāo)是矢徑 rho。) subplot(1,1,1) loglog semilogx semilogy polar(theta,rho) 取消子圖 繪出以 log10log10為坐標(biāo)刻度的對數(shù)圖 使用半對數(shù)刻度繪圖， x軸為 log10刻度， y軸為線性刻度。39。39。 title(39。r39。bar(t,y)39。) subplot (2,2,3), bar(t,y)。 title(39。stem(t,y)39。) 繪脈沖圖 繪階梯圖 繪條形圖 繪誤差條形圖 繪直方圖 顏色標(biāo)注符 r subplot (2,2,1), stem(t,y)。 figure 創(chuàng)建圖形窗口 shg 顯示圖形 圖形 窗的 控制 gcf 獲取當(dāng)前圖形窗的句柄 refresh 刷新圖形 clf 清除當(dāng)前圖形窗 close 關(guān)閉圖形窗 axes 在任意位置創(chuàng)建坐標(biāo)系 ishold 保持當(dāng)前圖形狀態(tài)為真 軸系 的 控制 gca 獲取當(dāng)前坐標(biāo)系的句柄 box 形成軸系方向 cla 清除當(dāng)前坐標(biāo)系 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 獲得當(dāng)前對象的句柄 copyobj 為圖形對象及其子項(xiàng)作硬拷貝 工具 closereq 請求關(guān)閉圖形窗 ishandle 是圖形句柄時為真 newplot 說明 Nextplot的 M文件 ginput 從鼠標(biāo)作圖形輸入 uiputfile 給出存儲文件的對話框 雜項(xiàng) 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。 關(guān)閉所有圖形窗口。 figure clf hold hold on hold off close close all subplot(n,m,p) 打開圖形窗口； 清除當(dāng)前圖形窗的內(nèi)容； 保持當(dāng)前圖形窗的內(nèi)容； 再次用 hold就解除保持狀態(tài)。只要每張圖的自變量和因變量長度相同就可以。 hold off 畫完一張圖用命令保持住，再畫下一條曲線。 y2=exp(*t2).*sin(5*t2+1)。)。 plot(t,y), hold on ,plot(t,y1,39。) Plotyy設(shè)有兩個縱坐標(biāo)，可以繪制兩個 y尺度不同的變量， x仍只用同一比例尺。),gtext(39。 gtext(39。t, t239。,t2,y3)。 為曲線圖加標(biāo)題 y3=5*y2。Y39。時間 39。)。 title(39。:r39。+g39。不便于用戶自行設(shè)定線型和顏色。 y的行（列）數(shù)就是曲線的根數(shù)。 plot(t,[y。 y=exp(*t).*sin(t)。+r39。:y39。*b39。) t是橫坐標(biāo)， y為縱坐標(biāo) y1相位超前 1弧度。 加標(biāo)題， 坐標(biāo)軸說明 加坐標(biāo)網(wǎng)格線 t=0::4*pi y=exp(*t).*sin(t) plot(t,y) y1=exp(*t).*sin(t+1) plot(t,y1,39。y39。x39。my first plot39?？?39。偶 39。奇 39。n=39?？?39。偶 39。奇 39。n=39。 case語句可以有 N1個，如果沒有任何一個 case值能與表達(dá)式值相符，則執(zhí)行 otherwise后面的語句組N。 break 是中止循環(huán)的命令，在多重循環(huán)中， break只能使程序跳出包含它的最內(nèi)部的那個循環(huán)。 end, end, h n =5 h = 1 1/2 1/3 1/4 1/5 1/2 1/3 1/4 1/5 1/6 1/3 1/4 1/5 1/6 1/7 1/4 1/5 1/6 1/7 1/8 1/5 1/6 1/7 1/8 1/9 增加可讀性 format rat, n=input(‘n=’) for i=1:n for j=1:n, h(i, j)=1/(i+j1)。n=39。 （ 3） for語句 for k= 初值：增量：終值 語句組 A， end 將語句組 A反復(fù)執(zhí)行 N次，每次執(zhí)行時程序中的 k值不同。 while 1+y1, y1=y y=y/2。 X不