資源描述:
《matlab應(yīng)用matlab數(shù)據(jù)及運(yùn)算》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。
1、第2章MATLAB數(shù)據(jù)及運(yùn)算2.1變量及其操作一、變量命名規(guī)則1.變量名、函數(shù)名對(duì)字母大小寫是敏感的myfile與MyFile表示不同的變量sin是MATLAB定義的正弦函數(shù)名,但SIN、Sin都不是2.變量名的第一個(gè)字符必須是英文字母3.變量名最多可包含63個(gè)字符(英文、數(shù)字和下劃線)4.變量名中不能包含空格、標(biāo)點(diǎn)my_exemplel2是合法的變量名,12exemple、_exemplel2^myexemplel2>my.exemplel2是非法變量名二、MATLAB默認(rèn)的預(yù)定義變量每當(dāng)MATLAB啟動(dòng)時(shí),不經(jīng)定義和賦值就會(huì)產(chǎn)生一些變
2、量,稱為MATLAB默認(rèn)的預(yù)定義變量這些變量都口J以重新賦值。但最好不要對(duì)這些變量名重新賦值預(yù)定義變量含義偵定義變量含義ans計(jì)算結(jié)果的缺省變量名nargin函數(shù)輸入?yún)?shù)個(gè)數(shù)eps機(jī)器零閾值(機(jī)器的浮點(diǎn)運(yùn)算誤差限)narout函數(shù)輸出參數(shù)個(gè)數(shù)Inf或inf無(wú)窮大,如1/0realmax最大正浮點(diǎn)數(shù)1或J虛數(shù)單位realmin最小正浮點(diǎn)數(shù)Pi圓周率兀lasterr存放最新的錯(cuò)誤信息NaN或nan不是一個(gè)數(shù)(NotaNumber),如0/0,oo/oolastwarn存放最新的警告信息例1用鍵盤在MATLAB指令窗中輸入以下內(nèi)容epseps—
3、—機(jī)器的浮點(diǎn)運(yùn)算誤差限。PC機(jī)上eps的默認(rèn)值為2.2204x1O'16,若某個(gè)量的絕對(duì)值小于eps,則可以認(rèn)為這個(gè)量為0。例2用鍵盤在MATLAB指令窗中輸入以下內(nèi)容1/0,I.elOOO,log(0)Inf——無(wú)窮大量+oo的MATLAB表示,也可以寫成inf。同樣地,?oo可以表示為?Inf。在MATLAB程序執(zhí)行時(shí),即使遇到了以0為除數(shù)的運(yùn)算,也不會(huì)終止程序的運(yùn)行,而只給出一個(gè)“除0”警告,并將結(jié)果賦成Inf,這樣的定義方式符合IEEE的標(biāo)準(zhǔn)。從數(shù)值運(yùn)算編程角度看,這樣的實(shí)現(xiàn)形式明顯優(yōu)于C語(yǔ)言。例3用鍵盤在MATLAB指令窗屮輸入
4、以下內(nèi)容0/0,inf/inf,inf*O注意在MATLAB中,即使遇到以0為除數(shù)的運(yùn)算,程序也不會(huì)終止運(yùn)行。這時(shí)只給出一個(gè)警告,并將結(jié)果賦給inf或NaNNaN不定式(notanumber),通常由0/0運(yùn)算、Inf/Inf及其他可能的運(yùn)算得出。NaN是一個(gè)很奇特的量,如NaN與Inf的乘積仍為NaNo例4用鍵盤在MATLAB指令窗中輸入以下內(nèi)容n=realmax,n=1.7977e+309例5用鍵盤在MATLAB指令窗中輸入以下內(nèi)容m=realmin,m=2.2251e?309i和j——若i或j量不被改寫,則它們表示純虛數(shù)量j。但在M
5、ATLAB程序編寫過(guò)程屮經(jīng)常事先改寫這兩個(gè)變量的值,比如在循環(huán)過(guò)程屮常用這兩個(gè)變量來(lái)表示循環(huán)變量。如果想恢復(fù)該變量,則可以用語(yǔ)句i=sqrt(-l)設(shè)置,即對(duì)-1求平方根。三、MATLAB運(yùn)算符和表達(dá)式運(yùn)算符操作數(shù)學(xué)表達(dá)式MATLAB運(yùn)算符MATLAB表達(dá)式加a+b+a?b減a-b-a-b乘axb*a*b除a^b/或a/b或ba幕abAaAb注意所有運(yùn)算都定義在復(fù)數(shù)域上例1用鍵盤在MATLAB指令窗中輸入以下內(nèi)容u=sqrt(-4)注意對(duì)于方根,運(yùn)算只給岀處于笫1象限的解用左除()和右除(/)計(jì)算表達(dá)式5+cos47°1+V7-2/
6、的值用鍵盤在MATLAB指令窗中輸入以下內(nèi)容(5+cos(47*pi/180))/(l+sqrt(7)?2*i)、(l+sqrt(7)-2*i)(5+cos(18047*pi))注意對(duì)牛標(biāo)量,左除()或右除(/)所得到的結(jié)果是一樣的對(duì)于矩陣,左除()和右除(/)所得到的結(jié)果是不一樣的AX=B-A_,AX=A_1B-X=A_1B=AB,XA=BfXAA_1=BA_,->X=BA-」B/A表達(dá)式表達(dá)式由變量名、運(yùn)算符、函數(shù)名組成表達(dá)式按優(yōu)先級(jí)和自左向右的順序運(yùn)算優(yōu)先級(jí)從高向低的排序是:指數(shù)運(yùn)算、乘除運(yùn)算、加減運(yùn)算括號(hào)可以改變運(yùn)算的順
7、序賦值符“二”和運(yùn)算符的兩側(cè)允許有空格復(fù)數(shù)和復(fù)數(shù)矩陣MATLAB把復(fù)數(shù)看成一個(gè)整體來(lái)處理(不像其它程序,把實(shí)部和虛部分開處理)虛數(shù)單位用預(yù)定義變量i或j表示復(fù)數(shù)z=a+bi=“i°直角坐標(biāo)表示與極坐標(biāo)表示之間的轉(zhuǎn)換real(z):給出復(fù)數(shù)z的實(shí)部a=rcos0imag(z):給出復(fù)數(shù)z的虛部b=rsin0abs⑵:給岀復(fù)數(shù)z的模r=sqrt(a2+b2)angle(z):以弧度為單位給出復(fù)數(shù)z的相角0=arctg(b/a)兀.例1寫岀復(fù)數(shù)Zj=3+4z,z?=1+2z,z3=2^6表達(dá)式,計(jì)算z-hll.用鍵盤在MATLAB指令窗屮輸入以
8、下內(nèi)容zl=3+4i%書寫格式1z2=1+2*i%書寫格式2z3=2*exp(i*pi/6)%書寫格式2z=zl*z2/z3注意%是注釋號(hào),在它后面及冋車之前的輸入部分是非執(zhí)行的注釋在格式1屮