第頁實驗01講評、參考答案講評未交實驗報告的同學名單數學：6人（11、12級）信科：12-04,12-22,13-47批改情況：問題1：不仔細，式子中出錯。問題2：提交的過程不完整。問題3：使用語句尾分號(;)不當，提交的過程中不該顯示的結果顯示。問題4：截屏窗口沒有調整大小。

附參考答案：《MATLAB軟件》課內實驗《MATLAB軟件》課內實驗王平實驗01MATLAB運算基礎（第2章MATLAB數據及其運算）一、實驗目的1.熟悉啟動和退出MATLAB的方法。2.熟悉MATLAB命令窗口的組成。3.掌握建立矩陣的方法。4.掌握MATLAB各種表達式的書寫規則以及常用函數的使用。二、實驗內容1.數學表達式計算先求下列表達式的值，然后顯示MATLAB工作空間的使用情況并保存全部變量。1.1計算三角函數（注意：度要轉換成弧度，e2如何給出）示例：點擊CommandWindow窗口右上角的，將命令窗口提出來成懸浮窗口，適當調整窗口大小。命令窗口中的執行過程：1.2計算自然對數，其中（提示：clc命令擦除命令窗口，clear則清除工作空間中的所有變量，使用時注意區別，慎用clear命令。應用點乘方）命令窗口中的執行過程：1.3求數學表達式的一組值提示：利用冒號表達式生成a向量，求各點的函數值時用點乘運算。命令窗口中的執行過程：1.4求分段函數的一組值，其中t=0:0.5:2.5提示：用邏輯表達式求分段函數值。命令窗口中的執行過程：1.5對工作空間的操作接著顯示MATLAB當前工作空間的使用情況并保存全部變量提示：用到命令who,whos,save,clear,load，請參考教材相關內容。命令窗口中的執行過程：>>who%顯示當前工作空間中的變量名Yourvariablesare:atxz1z2z3z4>>whos%顯示當前工作空間中的變量名及信息NameSizeBytesClassAttributesa1x61488doublet1x648doublex2x264doublecomplexz11x18doublez22x264doublecomplexz31x61976doublecomplexz41x648double>>savevarfile%生成文件var，并將當前工作空間的變量存入其中>>clearall%清除當前工作空間中的全部變量>>who>>whos>>loadvarfile%將文件var裝入>>whoYourvariablesare:atxz1z2z3z4>>whosNameSizeBytesClassAttributesa1x61488doublet1x648doublex2x264doublecomplexz11x18doublez22x264doublecomplexz31x61976doublecomplexz41x648double2.矩陣運算及數組運算已知：2.1矩陣加、減和數乘運算A+6*B和A-B+I（其中I為單位矩陣。注意：如何得到單位矩陣）命令窗口中的執行過程：2.2矩陣乘積和點乘積A*B和A.*B命令窗口中的執行過程：2.3矩陣的乘方和點乘方A^3和A.^3命令窗口中的執行過程：2.4矩陣的右除和左除A/B及B\A命令窗口中的執行過程：2.5拼接矩陣成大的矩陣[A,B]和[A([1,3],:);B^2]命令窗口中的執行過程：3.矩陣乘積、矩陣的子矩陣設有矩陣A和B3.1求矩陣A和B的乘積求它們的乘積并賦給C。（提示：可簡化A的輸入，用冒號表達式、reshape、矩陣轉置）命令窗口中的執行過程：3.2求矩陣C的子矩陣將矩陣C的右下角3×2子矩陣賦給D。命令窗口中的執行過程：3.3查看MATLAB工作空間的使用情況命令窗口中的執行過程：4.完成下列數據操作4.1求[100,999]之間能被21整除的數的個數提示：先利用冒號表達式，再利用find和length函數。命令窗口中的執行過程：4.2建立一個字符串向量，刪除其中的大寫字母提示：利用find函數和空矩陣。命令窗口中的執行過程：三、實驗提示1.1計算三角函數提示（注意：度要轉換成弧度，e2如何給出）提示1：示例：點擊CommandWindow窗口右上角的，將命令窗口提出來成懸浮窗口，適當調整窗口大小。命令窗口中的執行過程：提示2：e2的表示請查表“2.3常用數學函數及其含義”。用exp(2)。1.2計算自然對數提示，其中提示1：clc命令擦除命令窗口，clear則清除工作空間中的所有變量，使用時注意區別，慎用clear命令。提示2：ln和開方的表示請查“表2.3常用數學函數及其含義”。用log,sqrt。提示3：x2是數組運算。用x.^2或x.*x。提示4：i為復數的虛數單位（j也是）。若將i作變量則其虛數單位無效，恢復用cleari。慎用i,j做變量。1.4求分段函數的一組值提示，其中t=0:0.5:2.5提示：用邏輯表達式求分段函數值。提示1：參考例2.4。提示2：對于第1個分段條件表達式0≤t<1，有以下結果：對于第2、3個分段條件表達式，有類似的結果。于是，可以得到：分段函數值＝第1個表達式*第1個分段條件表達式+第2個表達式*第2個分段條件表達式+第3個表達式*第3個分段條件表達式+…注意：分段區間是不重疊的。對于在定義區間中任給的變量值，它只滿足一個分段條件表達式（值為1），其它不滿足（值為0），滿足的對應項的值是其表達式的值，不滿足的對應項值為0。代入本題得：z4=t^2*(t>=0&&t<1)+(t^2-1)*(t>=1&&t<2)+(t^2+2*t+1)*(t>=2&&t<3)繼續做下去可以完成本題。若變量取100個值，求對應的函數值，還按這種方法來求做就有問題啦。MATLAB提供了一種簡單的方法：t不是一個個取值，而是把所有值放到t中成一個行向量，表達式自動逐個將t中的值代入表達式計算，得到對應的一組函數值。實現這種功能需要將上面的表達式中的乘方號(^)和乘號(*)前加點(.)，即改成(.^)和(.*)，相當于循環，它不同于線性代數中相應的矩陣乘方(^)和乘法(*)運算。用冒號表達式更簡單：以下是分解式：至此，能理解表達式的計算結果了嗎？1.5對工作空間的操作提示接著顯示MATLAB當前工作空間的使用情況并保存全部變量提示1：用到命令who,whos,save,clear,load，請參考教材相關內容。步驟：查看工作空間；保存工作空間；清除工作空間；再查看工作空間；裝入工作空間；再查看工作空間。2.1矩陣加、減和數乘運算提示A+6*B和A-B+I（其中I為單位矩陣。注意：如何得到單位矩陣）提示1：3階單位矩陣用eye(3)或eye(3,3)。4.完成下列數據操作4.1求[100,999]之間能被21整除的數的個數提示提示：先利用冒號表達式，再利用find和length函數。提示1：步驟：用冒號表達式把[100,999]之間的整數放在一個向量里（如A）；用mod求得一個向量（如B），使A中能被21整除的數在B中對應位置為1，B的其它位置為0；用find求得B中不為0的元素的序號，放入一個向量中（如K）；顯示A中對應K的元素；%即能被21整除的數用length顯示K的元素個數。%即能被21整除的數的個數4.2建立一個字符串向量，刪除其中的大寫字母提示提示：利用find函數和空矩陣。提示1：利用find函數和空矩陣。注意使用關系運算、邏輯運算。注意！本題中大寫字母所在位置是隨機的，要用關系運算和邏輯運算，以及find找到它們所在位置的序號。提示2：四、教程：第2章MATLAB數據及其運算2.1MATLAB數據的特點p17矩陣是MATLAB最基本、最重要的數據對象。MATLAB的大部分運算或命令都是在矩陣（復數域上）運算的意義下執行的。標量對應僅含一個元素的矩陣。數值數據雙精度數 占64位，轉換函數double單精度數 占32位，轉換函數single帶符號整數 轉換函數int8,int16,int32無符號整數 轉換函數uint8,uint16,uint32字符數據 轉換函數char結構體類型 (Structure)單元類型 (Cell)邏輯型 非0為true，0為false2.2變量及其操作2.2.1變量及賦值1．變量命名以字母開頭，后接字母、數字或下劃線的字符序列，最多63個字符。區分字母的大小寫。MATLAB提供的標準函數名以及命令名必須用小寫字母。2．賦值語句(1)變量=表達式(2)表達式（值賦給預定義變量ans）其中，表達式是用運算符將有關運算量連接起來的式子，其結果是一個矩陣。例2.1計算表達式的值p18將的計算結果賦給變量x，然后顯示出結果。>>x=(5+cos(47*pi/180))/(1+sqrt(7)-2*i)x=1.1980+0.6572i其中，pi和i都是MATLAB預先定義的變量，分別代表圓周率π和虛數單位。3.預定義變量p18表2.1常用的預定義變量及其含義預定義變量含義ansepspii,jinf,InfNaN,nannarginnargoutrealmaxrealminlasterrlastwar默認賦值變量機器零閾值2.2204e-016π近似值虛數單位無窮大，如1/0的結果非數，如0/0,inf/inf的結果函數輸入參數個數函數輸出參數個數最大正實數最小正實數存放最新的錯誤信息存放最新的警告信息預定義變量有特定的含義，應盡量避免重新賦值。i=2 %對預定義變量重新賦值2*icleari %恢復預定義變量2*ii=2ans=4ans=0+2.0000i2.2.2變量的管理p191．內存變量的顯示及刪除顯示工作空間中駐留的變量名清單及信息命令：who %顯示變量名清單whos %顯示變量名清單及信息刪除工作空間中的變量命令：clear變量名工作空間窗口專門用于內存變量的管理。在該窗口中顯示所有內存變量的屬性。當選中某些變量后，再單擊Delete按鈕，將刪除這些變量。雙擊變量后，將進入變量編輯器。可以觀察變量中的具體元素，也可修改具體元素。輸入較大矩陣時，可采用變量編輯器。2．內存變量文件用MAT文件可把當前工作空間中的一些有用變量長久地保留下來，擴展名是.mat。MAT文件的生成和裝入，常用格式為：save文件名[變量名表][-append][-ascii]load文件名[變量名表][-ascii]文件名可帶路徑，不需帶擴展名.mat，默認對.mat文件進行操作。變量名表變量名以空格分隔。省略時，保存或裝入全部變量。-ascii選項使文件以ASCII格式處理，省略時以二進制格式處理。-append選項使變量追加到MAT文件中。2.2.3數據的輸出格式p20用十進制數表示一個常數，可采用日常記數法和科學記數法。一般情況下，內部每一個數據元素都用雙精度數表示和存儲。設置或改變數據輸出格式format命令的格式為：format格式符格式符決定數據的輸出格式表2.2控制數據輸出格式格式符及含義格式符含義short（默認）小數點后4位，不超過7位有效數字。大于1000的實數用5位有效數字的科學記數法long15位有效數字shorte5位有效數字科學記數法longe15位有效數字科學記數法shortg從short和shorte選擇最佳方式longg從long和longe選擇最佳方式rat近似有理數表示hex十六進制表示+正數、負數、零分別用+、-、空格表示bank銀行格式，元、角、分表示compact輸出變量之間無空行loose輸出變量之間有空行可用Helpformat查詢例輸出格式format>>4/3 %當前輸出格式ans=1.3333>>formatlong %15位有效數字形式輸出>>4/3ans=1.333333333333333>>formatrat %近似有理數表示>>4/3ans=4/3>>formatcompact %輸出變量之間沒有空行>>4/3ans=4/3>>formatloose %輸出變量之間有空行>>4/3ans=4/3>>formatshort %short為默認輸出格式2.3矩陣的表示p21表矩陣操作函數及其含義函數名含義eyeoneslinspacesub2indind2subsizelengthreshapeend[]單位矩陣全1矩陣生成行向量下標轉換成序號序號轉換成下標給出矩陣的行數和列數給出矩陣行數和列數中較大者矩陣重排預定義變量，某一維末尾下標空矩陣2.3.1矩陣的建立1.直接輸入法從鍵盤直接輸入矩陣的元素。方法如下：將矩陣的元素用方括號括起來，輸入元素；同一行的元素間用空格或逗號分隔；不同行的元素間用分號（或回車）分隔。例>>A=[1,2,3;456;7,89]A=1234567892.利用M文件建立矩陣對于比較大且復雜的矩陣，可專門建立一個M文件。例2.2利用M文件建立矩陣(1)啟動文本編輯器，輸入：MYMAT=[101,102,103,104,105; 201,202,203,204,205; 301,302,303,304,305](2)存盤(文件名為mymatrix.m)。(3)在命令窗口中輸入mymatrix，即運行該M文件，就建立一個名為MYMAT的矩陣。3.建立大矩陣（矩陣拼接）大矩陣可由方括號中的小矩陣或向量建立。>>A=[1,2;3,4]A=1234>>eye(2)ans=1001>>ones(2)ans=1111>>C=[A,eye(2),ones(2),A]C= 1 2 1 0 1 1 1 2 3 4 0 1 1 1 3 42.3.2冒號表達式用冒號表達式產生行向量，一般格式：e1:e2:e3其中e1為初始值，e2為步長，e3為終止值上限。用linspace函數產生行向量。調用格式：linspace(a,b,n)其中a和b是第一個和最后一個元素，n是元素總數。2.3.3矩陣的拆分1．矩陣元素通過下標引用矩陣的元素，例>>A=[1,2,3;4,5,6]A=123456>>A(2,3)ans=6>>A(2,2)=22;>>AA=1234226>>A(4,5)=10;>>AA=1230042260000000000010采用矩陣元素的序號來引用矩陣元素。元素的序號是相應元素在內存中的排列順序。在MATLAB中，矩陣按列存儲，先第一列，再第二列，依次類推。>>A=[1,2,3;4,5,6]A=123456>>A(3)ans=2>>A(5)ans=3>>A(8)???AttemptedtoaccessA(8);indexoutofboundsbecausenumel(A)=6.numel(A)返回A中元素個數。序號(Index)及下標(Subscript)一一對應。以m×n矩陣A為例，矩陣元素A(i,j)的序號為(j-1)*m+i。互轉換關系用sub2ind和ind2sub函數求得。>>A=ones(3,4)A=111111111111>>size(A) %給出矩陣的行數和列數ans=34>>sub2ind(size(A),1,2)%下標(1,2)轉換為序號ans=4>>[i,j]=ind2sub(size(A),3)%序號轉換為下標i=3j=1>>length(A) %給出矩陣行數和列數中較大者ans=4reshape(A,m,n) 在矩陣總元素保持不變的前提下，將矩陣重新排成m×n矩陣。>>x=1:1:12;>>y=reshape(x,3,4)y=147102581136912>>z=reshape(y,2,7)???Errorusing==>reshapeToRESHAPEthenumberofelementsmustnotchange.2.矩陣拆分(1)利用冒號表達式獲得子矩陣已知A為矩陣A(:,j)第j列全部元素；A(i,:)第i行的全部元素；A(i,j)第i行、第j列的元素；A(i:i+m,:)第i～i+m行的全部元素；A(:,k:k+m)第k～k+m列的全部元素；A(i:i+m,k:k+m)第i～i+m行內，并在第k～k+m列中的所有元素；A(:)將每一列元素堆疊起來，成為一個列向量。end表示某一維的末尾元素下標。A(end,:)最后一行A(:,end)最后一列(2)利用空矩陣刪除矩陣的元素給變量X賦空矩陣的語句：X=[]注意，X=[]及clearX不同：clear是將X從工作空間中刪除，而空矩陣則存在于工作空間中，只是維數為0。2.4MATLAB數據的運算p262.4.1算術運算1．基本算術運算MATLAB的基本算術運算有：＋ 加－ 減* 乘/ 右除\ 左除^ 乘方注意，運算是在矩陣意義下進行的。(1)矩陣加減運算假定矩陣A和B，則A+B和A-B實現矩陣的加減運算。運算規則是：若A和B的維數相同，則可以執行矩陣的加減運算，A和B的相應元素相加減。若A及B的維數不相同，將給出錯誤信息。注意，A、B均不是標量。(2)矩陣乘法若A為m×n矩陣，B為n×p矩陣，則C=A*B為m×p矩陣。(3)矩陣除法有兩種矩陣除法運算：\ 左除/ 右除若A矩陣是非奇異方陣，則A\B和B/A運算可以實現。A\B等效于inv(A)*BB/A等效于B*inv(A)對于含有標量的運算，兩種除法運算的結果相同。如3/4和4\3有相同的值，都等于0.75。又如設a=[10.5,25]，則a/5=5\a=[2.10005.0000]對于矩陣運算，一般A\B≠B/A。(4)矩陣的乘方一個矩陣的乘方運算可表示成A^x要求A為方陣，x為標量。2．點運算有一種特殊的運算，因為其運算符是在有關算術運算符前面加點，所以叫點運算。點運算符有.* 點乘./ 點右除.\ 點左除.^ 點乘方兩矩陣進行點運算是指它們的對應元素進行相關運算，要求兩矩陣的維數相同。3.MATLAB常用數學函數p29函數的自變量規定為矩陣變量。運算法則是將函數逐項作用于矩陣的元素上。結果是一個及自變量同維數的矩陣。表2.3常用數學函數及其含義p29函數名含義三角函數sqrtloglog10log2exppow2absanglerealimagconjremmodfixfloorceilroundsigngcdlcm略平方根函數自然對數函數常用對數函數以2為底的對數函數自然指數函數2的冪絕對值函數復數的幅角復數的實部復數的虛部復數共軛運算求余數或模運算模運算向零方向取整不大于自變量的最大整數不小于自變量的最小整數四舍五入到最鄰近的整數符號函數最大公因子最小公倍數函數使用說明：(1)三角函數以弧度為單位計算。(2)abs函數可以求實數的絕對值、復數的模、字符串的ASCII碼值。>>abs(-10)ans=10>>abs(3+4i)ans=5>>abs('abc')ans=979899(3)用于取整的函數有fix 向0方向取整floor 向-∞方向取整ceil 向+∞方向取整round 四舍五入到最鄰近的整數(4)rem(x,y)和mod(x,y)要求x,y必須為相同大小的實矩陣或為標量。rem 求余數或模運算mod 模運算當y≠0時rem(x,y)=x-y.*fix(x./y)mod(x,y)=x-y.*floor(x./y)當y=0時rem(x,0)=NaN（非數）mod(x,0)=x2.4.2關系運算p306種關系運算符：< 小于<= 小于或等于> 大于>= 大于或等于== 等于～= 不等于運算法則：兩個標量的比較若關系成立，結果為1，否則為0。兩個同維數矩陣的比較對相同位置的元素比較。標量及矩陣的比較標量及矩陣的每一個元素比較。例2.3建立5階方陣，判斷A的元素是否能被3整除。A=[24,35,13,22,63;23,39,47,80,80;... 90,41,80,29,10;45,57,85,62,21;... 37,19,31,88,76]P=rem(A,3)==0 %判斷元素是否被3整除A=24351322632339478080904180291045578562213719318876P=10001010001000011001000002.4.3邏輯運算p31邏輯運算符：& 及| 或～ 非運算法則：(1)在邏輯運算中，確認非零元素為真，用1表示，零元素為假，用0表示。(2)a和b是兩個標量a&ba,b全為非零時，結果為1，否則為0。a|ba,b中只要有一個非零，結果為1，否則為0。～a當a是零時，結果為1；當a非零時，結果為0。(3)兩個運算量是同維矩陣相同位置上的元素按標量規則逐個進行。(4)一個是標量，另一個是矩陣標量及矩陣中的每個元素之間按標量規則逐個進行。(5)邏輯非是單目運算符，也服從矩陣運算規則。(6)在算術、關系、邏輯運算中，算術運算優先級最高，邏輯運算優先級最低。例2.4在[0,3π]區間，求y=sinx的值。要求：①消去負半波，即(π,2π)內的函數值置0。②和內取均值為。%方法1x=0:pi/100:3*pi;y=sin(x);y1=(x<pi|x>2*pi).*y;%①消去負半波q= (x>pi/3&x<2*pi/3)|... (x>7*pi/3&x<8*pi/3);qn=~q;y2=q*sin(pi/3)+qn.*y1; %處理②plot(x,[y;y1;y2]);%方法2x=0:pi/100:3*pi;y=sin(x);y1=(y>=0).*y; %①消去負半波p=sin(pi/3);y2=(y>=p)*p+(y<p).*y1; %處理②plot(x,[y;y1;y2]);表2.4關系邏輯運算函數及其含義p32函數名含義allanyexistfindisemptyisglobalisinfisnanisfiniteissparseisstrxor向量的所有元素非0，結果為1向量中任何一個非0，結果為1變量在工作空間，結果為1找出矩陣中非0元素的位置變量是空矩陣，結果為1變量是全局變量，結果為1元素是±∞，對應位置取1元素不是一個數，對應位置取1元素值有限，對應位置取1變量是稀疏矩陣，結果為1變量是字符串，結果為1兩矩陣對應元素同為0或非0，對應位置取0例2.5建立矩陣A，然后找出在[10,20]區間的元素的位置（元素序號）。>>A=[4,15,-45,10,6;56,0,17,-45,0]A=415-4510656017-450>>find(A>=10&A<=20)ans=3672.5字符串p33字符串是用單撇號括起來的字符序列。將字符串當作一個行向量，每個元素對應一個字符，其標識方法和數值向量相同。可以建立多行字符串矩陣。字符串是以ASCII碼形式存儲的。abs和double函數都可以用來獲取字符串矩陣所對應的ASCII碼數值矩陣。相反，char函數可以把ASCII碼矩陣轉換為字符串矩陣。例2.6建立一個字符串向量，然后對該向量做如下處理：(1)取第1～5個字符組成的子字符串。(2)將字符串倒過來重新排列。(3)將字符串中的小寫字母變成相應的大寫字母，其余字符不變。(4)統計字符串中小寫字母的個數。clc;ch='ABc123d4e56Fg9';subch=ch(1:5)%取子字符串revch=ch(end:-1:1)%將字符串倒排k=find(ch>='a'&ch<='z')%找小寫字母的位置ch(k)='A'+(ch(k)-'a') %小寫字母變成大寫chchar(ch) %ASCII碼矩陣轉換為字符矩陣length(k) %統計小寫字母的個數char([0:29;30:59;60:89;90:119;120:149])subch=ABc12revch=9gF65e4d321cBAk=37913ch=ABC123D4E56FG9ch=ABC123D4E56FG9ans=ABC123D4E56FG9ans=4ans=




!"#$%&'()*+,-./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{|}~??????????????????????函數調用格式：eval(t)其中t為字符串。作用：把字符串的內容作為對應的MATLAB語句來執行。表2.5字符串處理函數及其含義函數名含義setstrcharmat2strnum2strint2strstr2numstrcatstrcmpabs將ASCII碼值轉換成字符同上，將代替setstr將矩陣轉換成字符串將數值轉換成字符串將整數轉換成字符串將字符串轉換成數值用于字符串的連接用于字符串的比較將字符轉換成ASCII碼(1)若字符串中含有單撇號，該撇號用兩個單撇號表示。(2)某些情況下用字符串向量表示，即用[]括起來。disp('I''mateacher.')f=70;c=(f-32)/1.8;disp(['Roomtemperatureis',...num2str(c),'degreesC.'])I'mateacher.Roomtemperatureis21.1111degreesC.2.6結構數據和單元數據p342.6.1結構數據1．結構矩陣的建立及引用結構矩陣其元素可以是不同的數據類型，將一組不同屬性的數據納入到一個統一的變量名下進行管理。建立結構矩陣可采用給結構成員賦值的辦法。格式為：結構矩陣名.成員名=表達式其中表達式應理解為矩陣表達式。例建立含有3個元素的結構矩陣>>a(1).x1=10;a(1).x2='liu'; a(1).x3=[11,21;34,78];>>a(2).x1=12; a(2).x2='wang';a(2).x3=[34,191;27,578];>>a(3).x1=14; a(3).x2='cai'; a(3).x3=[13,890;67,231];>>a(2).x3 %引用矩陣元素a(2)的成員x3ans=34191