數(shù)學(xué)實(shí)驗(yàn)2-MATLAB操作基礎(chǔ)

1、重要通知第三次課為教材第1-3章的課堂作業(yè)，占總成績(jī)的20%。望大家提前預(yù)習(xí)，將書上例題全部學(xué)習(xí)一遍，方可做好課堂作業(yè)。周末實(shí)驗(yàn)室開放，有數(shù)學(xué)建模國(guó)賽隊(duì)員輔導(dǎo)，望大家前往學(xué)習(xí)。實(shí)驗(yàn)二目的正確區(qū)分命令文件和函數(shù)文件的區(qū)別，掌握函數(shù)文件的組成，函數(shù)調(diào)用等掌握程序結(jié)構(gòu)組成，循環(huán)結(jié)構(gòu)，選擇結(jié)構(gòu)掌握Matlab 符號(hào)運(yùn)算 所謂M文件就是由MATLAB語(yǔ)言編寫的可在MATLAB語(yǔ)言環(huán)境下運(yùn)行程序源代碼文件。 M文件可以根據(jù)調(diào)用方式的不同分為兩類： 命令文件(Script File)：自動(dòng)重復(fù)執(zhí)行的一組MATLAB命令和函數(shù)組合，不需輸出輸入?yún)?shù)。 函數(shù)文件(Function File)：M文件的第一個(gè)可執(zhí)

2、行以function開始，便是函數(shù)文件，每一個(gè)函數(shù)文件定義一個(gè)函數(shù)。 2.1.1 M文件1. M文件概述2.1 MATLAB語(yǔ)言的流程結(jié)構(gòu)function output，output2=function-name（input1，input2） 2. M 文件的建立與打開 M文件是一個(gè)文本文件，擴(kuò)展名為“m”，可以用任何文本編輯器來(lái)建立和編輯，通常使用 Matlab 自帶的 M 文件編輯器。（1）新建一個(gè) M 文件方法一 菜單操作 ( File New M-File )方法二 命令操作 ( edit M 文件名 )方法三 快捷欄按鈕（2） 打開已有的 M 文件方法一 菜單操作 ( File Op

3、en )方法二 命令操作 ( edit M文件名 )方法三 快捷欄按鈕方法四 雙擊 M 文件 3. M 文件的運(yùn)行（1）從命令窗口 輸入文件名即可運(yùn)行。（2）從編輯器窗口方法一 在編輯器菜單Debug-save and run方法二 快捷欄按鈕方法三 快捷鍵F5注1 命令文件產(chǎn)生的變量是全局變量，函數(shù)文件內(nèi)部的變量（非輸入?yún)?shù)）是局部變量。注2 運(yùn)行若出現(xiàn)死循環(huán)或其它需要中斷程序運(yùn)行：1. 命令窗口使用組合鍵ctrl+c。2. 組合鍵ctrl+alt+del調(diào)出任務(wù)管理器，結(jié)束matlab進(jìn)程。例2.1.1 分別建立命令文件和函數(shù)文件，將華氏溫度f(wàn)轉(zhuǎn)換為攝氏溫度c。命令文件建立： （以文件名f

4、2c.m存盤） clear; %清除工作空間中的變量f=input(temperature：);c=5*(f-32)/9 f2cInput Fahrenheit temperature：73c = 22.7778命令窗口輸入：函數(shù)文件建立：以文件名f2c_fun.m存盤function c=f2c_fun(f)c=5*(f-32)/9; c=f2c_fun(20)c = -6.6667命令窗口輸入：注1：函數(shù)調(diào)用的一般格式是： 輸出實(shí)參表=函數(shù)名(輸入實(shí)參表)實(shí)參傳遞給形參例2.1.2 編寫函數(shù)文件求半徑為r的圓的面積和周長(zhǎng)。 函數(shù)文件如下： function s,p=fcircle(r) %

5、r 圓半徑%s 圓面%p 圓周長(zhǎng)s=pi*r*r;p=2*pi*r;2.2 程序控制結(jié)構(gòu) 程序控制結(jié)構(gòu)有三種：順序結(jié)構(gòu)、選擇結(jié)構(gòu)和循環(huán)結(jié)構(gòu)。任何復(fù)雜的程序都由這三種基本結(jié)構(gòu)組成。2.2.1 順序結(jié)構(gòu) 按排列順序依次執(zhí)行各條語(yǔ)句，直到程序的最后。 這是最簡(jiǎn)單的一種程序結(jié)構(gòu)，一般涉及數(shù)據(jù)的輸入輸出、數(shù)據(jù)的計(jì)算或處理等。1. 數(shù)據(jù)的輸入 從鍵盤輸入數(shù)據(jù)： A=input(提示信息，s)；加上s選項(xiàng)，則允許用戶輸入一個(gè)字符串。 xm=input(Whats your name?,s);Whats your name?張三2. 數(shù)據(jù)的輸出disp(輸出項(xiàng))其中輸出項(xiàng)既可以為字符串，也可以為矩陣。例2.2

6、.1 輸入x,y的值，并將它們的值互換后輸出。 x=input(Input x please:);y=input(Input y please:);z=x;x=y;y=z;disp(x);disp(y);Input x please:1Input y please:2 2 1例2.2.2 求一元二次方程ax2 +bx+c=0的根。 a=input(a=?);b=input(b=?); c=input(c=?); d=b*b-4*a*c; x=(-b+sqrt(d)/(2*a),(-b-sqrt(d)/(2*a);disp(x1=,num2str(x(1),x2=, num2str(x(2);a

7、=?1b=?3c=?2x1=-1,x2=-23. 程序的暫停pause(延遲秒數(shù))如果省略延遲時(shí)間，直接使用 pause，則將暫停程序，直到用戶按任一鍵后程序繼續(xù)執(zhí)行。2.2.2 選擇結(jié)構(gòu) 選擇結(jié)構(gòu)是根據(jù)給定的條件成立或不成立，分別執(zhí)行不同的語(yǔ)句。 MATLAB 用于實(shí)現(xiàn)選擇結(jié)構(gòu)的語(yǔ)句有 if 語(yǔ)句和 switch 語(yǔ)句。1. if語(yǔ)句 在MATLAB中，if語(yǔ)句有單分支，雙分支，多分枝3種格式。if 條件 語(yǔ)句組end當(dāng)條件1成立時(shí)，執(zhí)行語(yǔ)句組1，否則，當(dāng)條件2成立，執(zhí)行語(yǔ)句組2，否則，,否則，當(dāng)條件m成立，執(zhí)行語(yǔ)句組m，否則執(zhí)行語(yǔ)句組n，之后繼續(xù)執(zhí)行 if 語(yǔ)句的后繼語(yǔ)句，若條件不成立，則

8、直接執(zhí)行if語(yǔ)句的后繼語(yǔ)句。if 條件語(yǔ)句組1else語(yǔ)句組2endif 條件1 語(yǔ)句組1elseif 條件2 語(yǔ)句組2elseif 條件m 語(yǔ)句組melse 語(yǔ)句組nend例2.2.3 計(jì)算分段函數(shù)的值。 x=input(請(qǐng)輸入x的值:); if x c=input(請(qǐng)輸入一個(gè)字符,s); if c=A & c=a& c=0& c=9 disp(abs(c)-abs(0); else disp(c); end2. switch語(yǔ)句switch 表達(dá)式 %計(jì)算表達(dá)式的值 case 表達(dá)式1 %等于表達(dá)式1的值時(shí)，執(zhí)行語(yǔ)句組1 語(yǔ)句組1 case 表達(dá)式2 %等于表達(dá)式1的值時(shí)，執(zhí)行語(yǔ)句組2 語(yǔ)

9、句組2 case 表達(dá)式m %等于表達(dá)式m的值時(shí)，執(zhí)行語(yǔ)句組m 語(yǔ)句組m otherwise %不為上述值時(shí)，執(zhí)行語(yǔ)句組n 語(yǔ)句組n end %當(dāng)任意一個(gè)分支的語(yǔ)句執(zhí)行完后，直接執(zhí)行switch語(yǔ)句的下一句。例2.2.5 某商場(chǎng)對(duì)顧客所購(gòu)買的商品實(shí)行打折銷售，標(biāo)準(zhǔn)如下(商品價(jià)格用price來(lái)表示)： price200 沒(méi)有折扣 200price500 3%折扣 500price1000 5%折扣 1000price2500 8%折扣 2500price price=input(請(qǐng)輸入商品價(jià)格);switch fix(price/100) case 0,1 %價(jià)格小于200 rate=0; ca

10、se 2,3,4 %價(jià)格大于等于200但小于500 rate=3/100; case num2cell(5:9) %價(jià)格大于等于500但小于1000 rate=5/100; case num2cell(10:24) %價(jià)格大于等于1000但小于2500 rate=8/100; case num2cell(25:49) %價(jià)格大于等于2500但小于5000 rate=10/100; otherwise %價(jià)格大于等于5000 rate=14/100;endprice=price*(1-rate) %輸出商品實(shí)際銷售價(jià)格2.2.3 循環(huán)結(jié)構(gòu) 循環(huán)結(jié)構(gòu)是按照給定的條件，重復(fù)執(zhí)行指定的語(yǔ)句。 MATL

11、AB 用于實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的語(yǔ)句有 for 語(yǔ)句和 while 語(yǔ)句。1. for語(yǔ)句 for 循環(huán)變量=表達(dá)式1:表達(dá)式2:表達(dá)式3循環(huán)體語(yǔ)句end 其中表達(dá)式 1 的值為循環(huán)變量的初值，表達(dá)式2的值為步長(zhǎng)，表達(dá)式3的值為循環(huán)變量的終值。步長(zhǎng)為1時(shí)，表達(dá)式2可以省略。for 循環(huán)變量=矩陣表達(dá)式循環(huán)體語(yǔ)句end 執(zhí)行過(guò)程是依次將矩陣的各列元素賦給循環(huán)變量，然后執(zhí)行循環(huán)體語(yǔ)句，直至各列元素處理完畢。例2.2.6 已知 ，當(dāng)n=100時(shí)，求y的值。 y=0; n=100; for i=1:n y=y+1/(2*i-1); end y例2.2.7 寫出下列程序的執(zhí)行結(jié)果。 s=0; a=12,13,1

12、4;15,16,17;18,19,20;21,22,23; for k=a s=s+k; end disp(s);2. while語(yǔ)句 while語(yǔ)句的一般格式為： while (條件) 循環(huán)體語(yǔ)句 end 其執(zhí)行過(guò)程為：若條件成立，則執(zhí)行循環(huán)體語(yǔ)句，執(zhí)行后再判斷條件是否成立，如果不成立則跳出循環(huán)。例2.2.8 從鍵盤輸入若干個(gè)數(shù)，當(dāng)輸入0時(shí)結(jié)束輸入，求這些數(shù)的平均值和它們之和。 sum=0;cnt=0;val=input(Enter a number (end in 0):);while (val=0) sum=sum+val; cnt=cnt+1; val=input(Enter a nu

13、mber (end in 0):);endif (cnt 0) sum mean=sum/cntend3. break語(yǔ)句和continue語(yǔ)句 與循環(huán)結(jié)構(gòu)相關(guān)的語(yǔ)句還有 break 語(yǔ)句和 continue語(yǔ)句。 它們一般與if語(yǔ)句配合使用。 break 語(yǔ)句用于終止循環(huán)的執(zhí)行。當(dāng)在循環(huán)體內(nèi)執(zhí)行到該語(yǔ)句時(shí)，程序?qū)⑻鲅h(huán)，繼續(xù)執(zhí)行循環(huán)語(yǔ)句的下一語(yǔ)句。 continue 語(yǔ)句控制跳過(guò)循環(huán)體中的某些語(yǔ)句。當(dāng)在循環(huán)體內(nèi)執(zhí)行到該語(yǔ)句時(shí)，程序?qū)⑻^(guò)循環(huán)體中所有剩下的語(yǔ)句，繼續(xù)下一次循環(huán)。例2.2.9 求100，200之間第一個(gè)能被21整除的整數(shù)。 for n=100:200 if rem(n,21)=

14、0 continue end break end n4. 循環(huán)的嵌套 如果一個(gè)循環(huán)結(jié)構(gòu)的循環(huán)體又包括一個(gè)循環(huán)結(jié)構(gòu)，就稱為循環(huán)的嵌套，或稱為多重循環(huán)結(jié)構(gòu)。例2.2.10 若一個(gè)數(shù)等于它的各個(gè)真因子之和，則稱該數(shù)為完數(shù)，如6=1+2+3，所以6是完數(shù)。求 1,500之間的全部完數(shù)。 for m=1:500 s=0; for k=1:m/2 if rem(m,k)=0 s=s+k; end end if m=s disp(m); end end MathWorks公司在1993年收購(gòu)了主要針對(duì)符號(hào)計(jì)算的MAPLIE V 軟件的使用，以MAPLIE內(nèi)核為符號(hào)計(jì)算的引擎，開發(fā)了符號(hào)數(shù)學(xué)工具箱（ Symb

15、olic Math Toolbox ）。 MAPLIE占據(jù)符號(hào)軟件的主導(dǎo)地位，符號(hào)運(yùn)算工具包通過(guò)調(diào)用Maple軟件實(shí)現(xiàn)符號(hào)計(jì)算的。 在符號(hào)計(jì)算中，參與運(yùn)算的是符號(hào)變量而不是數(shù)值，使用字符串進(jìn)行分析。數(shù)值運(yùn)算中必須先對(duì)變量賦值，然后才能參與運(yùn)算。符號(hào)運(yùn)算無(wú)須事先對(duì)獨(dú)立變量賦值，運(yùn)算結(jié)果以標(biāo)準(zhǔn)的符號(hào)形式表達(dá)，可以獲得任意精度的解。2.3 MATLAB符號(hào)計(jì)算工具箱2.3.1 建立符號(hào)對(duì)象1. 建立符號(hào)變量和符號(hào)常量 創(chuàng)建一個(gè)符號(hào)變量x，可以是字符、字符串、表達(dá)式或字符表達(dá)式。 x=sym(符號(hào)字符串) syms var1 var2 varn函數(shù)定義符號(hào)變量var1,var2,varn注1 sym定

16、義符號(hào)變量一次只能定義一個(gè)。注2 var1varn之間用空格隔開2.3.2 建立符號(hào)表達(dá)式及符號(hào)矩陣?yán)?.3.1 用2種方法建立符號(hào)表達(dá)式。方法一 利用單引號(hào)來(lái)生成符號(hào)表達(dá)式方法二 用sym函數(shù)建立符號(hào)表達(dá)式方法三 使用已經(jīng)定義的符號(hào)變量組成符號(hào)表達(dá)式 U=sym(3*x2+5*y+2*x*y+6) %定義符號(hào)表達(dá)式Usyms x y; %建立符號(hào)變量x、yV=3*x2+5*y+2*x*y+6 %定義符號(hào)表達(dá)式V2*U-V+6 %求符號(hào)表達(dá)式的值例2.3.2 構(gòu)造符號(hào)矩陣。 S=sym(l,x,y,z;n,u,v,w;a,b,c,d;g,h,j,p) S = l, x, y, z n, u,

17、v, w a, b, c, d g, h, j, p H=sym(cos(t), -sin(t); sin(t), cos(t) H= cos(t), -sin(t) sin(t), cos(t)2.3.3 符號(hào)表達(dá)式運(yùn)算例2.3.3 求兩個(gè)符號(hào)表達(dá)式的和與差。 syms x fx gx % 定義符號(hào)變量 fx = 2*x2+3*x-5 fx = 2*x2+3*x-5 gx= x2-x+7； fx+gx ans= 3*x2+2*x+2 fx-gx ans= x2+4*x-12 syms a b c d e f g h; A = a b; c d; B = e f; g h; C1 = A.*B

18、 C1 = a*e, b*f c*g, d*h C2 = A*B/A C2 = (d*a*e+d*b*g-c*a*f-c*b*h)/(d*a-c*b), -(b*a*e+b2*g-a2*f-b*h*a)/(d*a-c*b) (d2*g+d*c*e-c2*f-c*d*h)/(d*a-c*b), -(d*b*g+b*c*e-c*a*f-d*h*a)/(d*a-c*b) 例 2.3.4 符號(hào)矩陣與數(shù)組的乘除示例。2.4 獨(dú)立變量與表達(dá)式化簡(jiǎn) 2.4.1 符號(hào)表達(dá)式中變量的確定MATLAB中的符號(hào)可以表示符號(hào)變量和符號(hào)常數(shù)。findsym可以幫助用戶查找一個(gè)符號(hào)表達(dá)式中的符號(hào)變量。該函數(shù)的調(diào)用格式為：findsym(S,n) 函數(shù)返回符號(hào)表達(dá)式S中的n個(gè)符號(hào)變量，若沒(méi)有指定n，則返回S中的全部符號(hào)變量。在求函數(shù)的極限、導(dǎo)數(shù)和積分時(shí)，如果用戶沒(méi)有明確指定自變量，MATLAB將按缺省原則確定主變量并