【導讀】本課研究物流管理的關鍵問題：物流管理決策。著名管理大師西蒙曾說：“管。系統(tǒng)設計上，側重于用戶界面設計、數(shù)據(jù)文件設計、數(shù)據(jù)分。析及圖形化界面、線性規(guī)劃、函數(shù)優(yōu)化及圖形分析等知識的學習。推薦參考書：1、張瑞豐，精通MATLAB，中國水利水電出版社，2020年；在業(yè)務中使用定量方法，常稱之為管理科學或運籌學。決策就是作決定、選滿意方案；廣義的決策還包括選擇前設計、構造，提出多種可行方案，事后對選定的方案進行實施、監(jiān)測、評價和鑒定。定性分析是基礎，定量分析則是定性分析的深化與精確化。析與定量分析兩者各有所長，不可偏廢。[MATLAB:取自英文matrix+laboratory，矩陣實驗室]. 其可視化編程能力較強，如運用DDE和OLE，實現(xiàn)了MATLAB. 宇處理于一體的高水準環(huán)境。邁入20世紀90年代后，MATLAB已經(jīng)成為國際公認的。標準計算軟件，至今沒有一個別的計算軟件可以與MATLAB匹敵。工具箱和學科性工具箱]。

　　

【正文】 矩陣的乘積，若出錯，則自動轉去求兩矩陣的點乘。 程序如下： 17 A=[1,2,3。4,5,6]。 B=[7,8,9。10,11,12]。 try C=A*B。 catch C=A.*B。 end C lasterr %顯示出錯原因 （ 3）循環(huán)結構 for 語句 for 語句的格式為： for 循環(huán) 變量 =表達式 1:表達式 2:表達式 3 循環(huán)體語句 end 其中表達式 1的值為循環(huán)變量的初值，表達式 2的值為步長，表達式 3的值為循環(huán)變量的終值。步長為 1時，表達式 2可以省略。 for 語句更一般的格式為： for 循環(huán)變量 =矩陣表達式 循環(huán)體語句 end 執(zhí)行過程是依次將矩陣的各列元素賦給循環(huán)變量，然后執(zhí)行循環(huán)體語句，直至各列元素處理完畢。 寫出下列程序的執(zhí)行結果。 s=0。 a=[12,13,14。15,16,17。18,19,20。21,22,23]。 for k=a s=s+k。 end disp(s39。)。 while 語句 while 語句的一般格式為： while (條件 ) 循環(huán)體語句 end 其執(zhí)行過程為：若條件成立，則執(zhí)行循環(huán)體語句，執(zhí)行后再判斷條件是否成立，若不成立則跳出循環(huán)。 例 從鍵盤輸入若干個數(shù)，當輸入 0 時結束輸入，求這些數(shù)的平均值和它們之和。 18 sum=0。t=0。val=input(39。Enter a number (end in 0):39。)。 while (val~=0) sum=sum+val。 t=t+1。 val=input(39。Enter a number (end in 0):39。)。 end if t 0 sum； mean=sum/t， disp(sum),end 與循環(huán)結構相關的語句還有 break 語句和 continue 語句。它們一般與 if 語句配合使用。 break 語句用于終止循環(huán)的執(zhí)行。當在循環(huán)體內(nèi)執(zhí)行到該語句時，程序將跳出循環(huán)，繼續(xù)執(zhí)行循環(huán)語句的下一語句。 continue 語句控制跳過循環(huán)體中的某些語句。當在循環(huán)體內(nèi)執(zhí)行到該語句時，程序將跳過循環(huán)體中所有剩下 的語句，繼續(xù)下一次循環(huán)。 例 求 [100， 200]之間第一個能被 21 整除的整數(shù)。 for n=100:200， if rem(n,21)~=0， continue end Break end， n 程序調(diào)試 一般來說，應用程序的錯誤有兩類，一類是語法錯誤，另一類是運行時的錯誤。 語法錯誤包括詞法或文法的錯誤，例如函數(shù)名的拼寫錯、表達式書寫錯等。 程序運行時的錯誤是指程序的運行結果有錯誤，這類錯誤也稱為程序邏輯錯誤。 程序調(diào)試器： Debug、 Breakpoints [即 M文件編輯器 ] Debug 菜單項：用 于程序調(diào)試，需要與 Breakpoints 菜單項配合使用。 Breakpoints 菜單項：該菜單項共有 6 個菜單命令，前兩個是用于在程序中設置和清除斷點的，后 4個是設置停止條件的，用于臨時停止 M文件的執(zhí)行，并給用戶一個檢查局部變量的機會。 實驗三 MATLAB 程序設計 MATLAB 繪圖 MATLAB 數(shù)據(jù)圖形化通過圖形的線型、平面、色彩、光線、視角等屬性的控制，能把數(shù)據(jù)的內(nèi)在特征表現(xiàn)得淋漓盡致，包括二維、三維，甚至四維的圖形。 二維數(shù)據(jù)曲線圖 19 三維圖形 圖像處理與動畫制作 （ 1）二維數(shù)據(jù)曲線圖 ①繪制二維曲線 plot 函數(shù)的基本調(diào)用格式為： plot(x,y) 其中 x 和 y為長度相同的向量，分別用于存儲 x坐標和 y坐標數(shù)據(jù)；當 x 是實向量時， plot(x)以該向量元素的下標為橫坐標，元素值為縱坐標畫出一條折線圖；當 x,y 是同維矩陣時，則以 x,y 對應列元素為橫、縱坐標分別繪制曲線，曲線條數(shù)等于矩陣的列數(shù)。 例 在 0≤ x≤ 2p 區(qū)間內(nèi)，繪制曲線 y=(4πx) 程序如下： x=0:pi/100:2*pi。 y=2*exp(*x).*cos(4*pi*x)。 plot(x,y) plot(x1,y1,x2,y2,?,xn,yn) *當輸入?yún)?shù)都為向量時， x1和 y1， x2和 y2， ? ， xn 和 yn 分別組成一組向量對，每一組向量對的長度可以不同。每一向量對可以繪制出一條曲線，這樣可以在同一坐標內(nèi)繪制出多條曲線。 *當輸入?yún)?shù)有矩陣形式時，配對的 x,y 按對應列元素為橫、縱坐標分別繪制曲線，曲線條數(shù)等于矩陣的列數(shù)。 例 分析下列程序繪制的曲線。 x1=linspace(0,2*pi,100)。x2=linspace(0,3*pi,100)。 y1=sin(x1)。y2=1+sin(x2)。 x=[x1。x2]39。y=[y1。y2]39。 plot(x,y,x1,y1) plotyy 繪圖函數(shù)可繪制出具有不同縱坐標標度的兩個圖形： plotyy(x1,y1,x2,y2) 其中 x1,y1 對應一條曲線， x2,y2 對應另一條曲線。橫坐標的標度相同，縱坐標有兩個，左縱坐標用于 x1,y1 數(shù)據(jù)對，右縱坐標用于 x2,y2 數(shù)據(jù)對。 例 用不同標度在同一坐標內(nèi)繪制曲線 y1=(4π x)和 y2=(π x)。 程序如下： x=0:pi/100:2*pi y1=*exp(*x).*cos(4*pi*x) y2=2*exp(*x).*cos(pi*x) 20 plotyy(x,y1,x,y2)。 hold on/off 命令控制是保持原有圖形還是刷新原有圖形，不帶參數(shù)的 hold命令在兩種狀態(tài)之間進行切換。 設置曲線樣式： “b .” 表示藍色點劃線， “y:d” 表示黃色虛線并用菱形符標記數(shù)據(jù)點。當選項省略時，線型一律用實線，顏色將根據(jù)曲線的先后順序依次。要設置曲線樣式可以在 plot 函數(shù)中加繪圖選項，其調(diào)用格式為： plot(x1,y1,選項 1,x2,y2,選項 2,? ,xn,yn,選項 n) ②圖形標注與坐標控制 圖形標注： title(圖形名稱 )、 xlabel(x軸說明 )、 ylabel(y軸說明 )、 text(x,y,圖形說明 )， legend(圖例 1,圖例 2,?) 例 在 0≤ x≤ 2p 區(qū)間內(nèi)繪制曲線 y1=，并添加圖形標注。 程序如下： x=0:pi/100:2*pi。y1=2*exp(*x)。 plot(x,y1) title(39。x from 0 to 2{\pi}39。)。 %加圖形標題 xlabel(39。Variable X39。)。 %加 X軸 說明 ylabel(39。Variable Y39。)。 %加 Y軸說明 text(,‘ 曲線 y1=2e^{}’ )。 %在指定位置添加圖形說明 legend(‘y1’) %加圖例 axis 函數(shù)的調(diào)用格式為： axis([xmin xmax ymin ymax zmin zmax]) axis 函數(shù)功能豐富，常用的格式還有： axis equal：縱、橫坐標軸采用等長刻度。 axis square：產(chǎn)生正方形坐標系 (缺省為矩形 )。 axis auto：使用缺省設置。 axis off：取消坐標軸。 axis on：顯示坐標軸。 grid on/off 命令控制是畫還是不畫網(wǎng)格線，不帶參數(shù)的 grid 命令在兩種狀態(tài)之間進行切換。 box on/off 命令控制是加還是不加邊框線，不帶參數(shù)的 box 命令在兩種狀態(tài)之間進行切換。 在同一坐標中，可以繪制 3 個同心圓，并加坐標控制。 程序如下： t=0::2*pi。x=exp(i*t)。y=[x。2*x。3*x]39。plot(y) grid on。 %加網(wǎng)格線 box on。 %加坐標邊框 21 axis equal %坐標軸采用等刻度 ③圖形的可視化編輯、圖形窗口的分割 利用圖形窗口菜單欄或工具欄中的有關命令可以完成對窗口中各種圖形對象的編輯處理。在圖形窗口上有一個菜單欄和工具欄。菜單欄包含 File、Edit、 View、 Insert、 Tools、 Window 和 Help 共 7個菜單項，工具欄包含11 個命令按鈕。 subplot 函數(shù)的調(diào)用格式為： subplot(m,n,p) 該函數(shù)將當前圖形窗口分成 m179。 n 個繪圖區(qū)，即每行 n 個，共 m 行，區(qū)號按行優(yōu)先編號，且 選定第 p個區(qū)為當前活動區(qū)。在每一個繪圖區(qū)允許以不同的坐標系單獨繪制圖形。 ④二維統(tǒng)計分析圖 ? 條形圖：是用水平條形或垂直條形來表示向量或矩陣中的值。 bar（ Y）：如果 Y為向量，則分別顯示每個分量的高度，橫坐標為 1到 LENGTH（ Y）；如果 Y 為矩陣，則 bar 函數(shù)對每一行中分量創(chuàng)建的條形進行分組，再分別繪出，橫坐標為 1 到 SIZE（ Y， 1），即矩陣的行數(shù)。 bar（ X， y）：在 X 指定的位置上繪制出 Y 中每一個分量的條形，其中 X為嚴格單調(diào)遞增向量。如果 Y為矩陣，則 bar 將矩陣分解成幾個行向量，在指定的橫坐標處分 別繪制出條形。 bar（ 178。178。178。 ， WIDTH）：設置條形的相對寬度并控制在一組內(nèi)條形的間距，WIDTH 默認值為 。 ? 二維散點圖、餅圖 SCATTER（ X， Y， S）：在向量 X 和 Y 指定的位置繪制彩色圓圈。其中 X和 Y必須具有相同的大小，參數(shù) S設置標記尺寸的大小。 pie（ X）：應用 X 中的數(shù)據(jù)繪制餅圖。如果輸入?yún)?shù)為向量，則繪制各個元素在向量所有元素之中所占比例；當輸入?yún)?shù)為矩陣時，繪制每個元素在矩陣所有元素之中所占的比例。 Pie（ X， explode）：將一個扇區(qū)從餅形圖中分離出，參數(shù) explode 為對 應于X 的為零或非零的向量或矩陣。與 explode 的非零值對應的部分將從餅形圖中心分離出來， explode 必須與 X具有相同大小。 例 ： x=[1 3 6 8 20],explode=[0 0 0 0 1] Pie(x, explode) 22 （ 2）三維圖形 ? 三維曲線 plot3 函數(shù)與 plot 函數(shù)用法十分相似： plot3(x1,y1,z1,選項 1,? ,xn,yn,zn,選項 n) 其中每一組 x,y,z 組成一組曲線的坐標參數(shù)，選項的定義和 plot 函數(shù)相同。當 x,y,z 是同維向量時，則 x,y,z 對應元素構成一 條三維曲線。當 x,y,z是同維矩陣時，則以 x,y,z 對應列元素繪制三維曲線，曲線條數(shù)等于矩陣列數(shù)。 例 繪制三維曲線。 程序如下： t=0:pi/100:20*pi。 x=sin(t)。y=cos(t)。z=t.*sin(t).*cos(t)。 plot3(x,y,z)。title(39。Line in 3D Space39。)。 xlabel(39。X39。)。ylabel(39。Y39。)。zlabel(39。Z39。)。grid on。 條形圖、餅圖等也可以三維形式出現(xiàn)： 常用格式為： bar3(y)、 pie3(x) 例 已知 x=[2347,1827,2043,3025]，繪制餅圖 x=[2347,1827,2043,3025] pie3(x) ? 視點處理： view(az,el) 其中 az 為方位角， el 為仰角，它們均以度為單位。系統(tǒng)缺省的視點定義為方位角 176。 ,仰角 30176。 （ 3）圖像處理與動畫制作 ①圖像處理 imread 和 imwrite 函數(shù)分別用于將圖像文件讀入 MATLAB 工作空間，以及將圖像數(shù)據(jù)和色圖數(shù)據(jù)一起寫入一定格式的圖像文件。 MATLAB 支持多種圖像文件格式，如 .bmp、 .jpg、 .jpeg、 .tif 等。 image 和 imagesc 函數(shù)用于圖像顯示。為了保證圖像的顯示效果，一般還應使用 colormap 函數(shù)設置圖像色圖。 例 有一圖像文件 ，在圖形窗口顯示該圖像。 程序如下： [x,cmap]=imread(39。39。)。 %讀取圖像的數(shù)據(jù)陣和色圖陣 image(x)。colormap(cmap)。 axis image off %保持寬高比并取消坐標軸 23 ②動畫制作 MATLAB 提供 moviein 和 movie 函數(shù)進行動畫制作。 moviein(n)函數(shù)用來建立一個足夠大的 n 列矩 陣。該矩陣用來保存 n 幅畫面的數(shù)據(jù)，以備播放。之所以要事先建立一個大矩陣，是為了提高程序運行速度。 movie(m,n)函數(shù)播放由矩陣 m所定義的畫面 n 次，缺省時播放一次。 例 繪制了 peaks 函數(shù)曲面并且將它繞 z 軸旋轉。 [X,Y,Z]=peaks(30)。su