【文章內(nèi)容簡(jiǎn)介】 ]來(lái)輸入數(shù)組，也可以用以下提供的函數(shù)來(lái)生成矩陣。 ones( ) 創(chuàng)建一個(gè)所有元素都為 1的矩陣，其中可以制定維數(shù) zeros() 創(chuàng)建一個(gè)所有元素都為 0的矩陣 eye() 創(chuàng)建對(duì)角元素為 1，其他元 素為 0的矩陣 diag() 根據(jù)向量創(chuàng)建對(duì)角矩陣，即以向量的元素為對(duì)角元素 magic() 創(chuàng)建魔方矩陣 rand() 創(chuàng)建隨機(jī)矩陣，服從均勻分布 randn() 創(chuàng)建隨機(jī)矩陣，服從正態(tài)分布 課件之家精心整理資料 歡迎你的欣賞 課件之家精心整 理資料 歡迎你的欣賞 randperm() 創(chuàng)建隨機(jī)行向量 horcat C=[A,B]，水平聚合矩陣，還可以用 cat(1,A,B) vercat C=[A。B]，垂直聚合矩陣 , 還可以用 cat(2,A,B) repmat(M,v,h) M在垂直方向上聚合 v次 水平方向上 聚合 h次 blkdiag（ A， B） 以 A，和 B 為塊創(chuàng)建塊對(duì)角矩陣 length 返回矩陣最長(zhǎng)維的的長(zhǎng)度 ndims 返回維數(shù) numel 返回矩陣元素個(gè)數(shù) size 返回每一維的長(zhǎng)度， [rows,cols]=size(A) reshape 重塑矩陣 rot90 旋轉(zhuǎn)矩陣 90 度，逆時(shí)針?lè)较? fliplr 沿垂 軸翻轉(zhuǎn)矩陣 flipud 沿水平軸翻轉(zhuǎn)矩陣 transpose 沿主對(duì)角線翻轉(zhuǎn)矩陣 ctranspose 轉(zhuǎn)置矩陣，也可用 A’或 A.’，這僅當(dāng)矩陣為復(fù)數(shù)矩陣時(shí)才有區(qū)別 inv 矩陣的逆 det 矩陣的行列式值 trace 矩陣對(duì)角元素的和 norm 矩陣或矢量的范數(shù)， norm（ a， 1）， norm（ a，Inf）?? . normest 估計(jì)矩陣的最大范數(shù)矢量 chol 矩陣的 cholesky 分解 cholinc 不完全 cholesky 分解 lu LU 分解 luinc 不完全 LU 分解 qr 正交分解 kron（ A， B） A 為 m n， B 為 p q，則生成 mp nq 的矩陣， A的每一個(gè)元素都會(huì)乘上 B，并占據(jù) p q大小的空間 rank 求出矩陣的秩 pinv 求偽逆矩陣 多項(xiàng)式運(yùn)算 多項(xiàng)式 Matlab 里面的多項(xiàng)式是以向量來(lái)表示的， 其具體操作函數(shù)如下： conv 多項(xiàng)式的乘法 deconv 多項(xiàng)式的除法，【 a， b】＝ deconv（ s），返回商和余數(shù) poly 求多項(xiàng)式的系數(shù)（由已知根求多項(xiàng)式的系數(shù)） polyeig 求多項(xiàng)式的特征值 課件之家精心整理資料 歡迎你的欣賞 課件之家精心整 理資料 歡迎你的欣賞 Polyfit（ x， y， n） 多項(xiàng)式的曲線擬合， x， y為被擬合的向量， n為擬合多 項(xiàng)式階數(shù) polyder 求多項(xiàng)式的一階導(dǎo)數(shù)， polyder（ a， b）返回 ab的導(dǎo)數(shù) [a,b]＝ polyder（ a， b）返回 a/b 的導(dǎo)數(shù)。 polyint 多項(xiàng)式的積分 polyval 求多項(xiàng)式的值 polyvalm 以矩陣為變量求多項(xiàng)式的值 residue 部分分式展開(kāi)式 roots 求多項(xiàng)式的根（返回所有根組成的武漢理工大學(xué)《 Matlab 課程設(shè)計(jì)》說(shuō)明書(shū) 課件之家精心整 理資料 歡迎你的欣賞 3. matlab 基本運(yùn)算 基礎(chǔ)微積分計(jì)算 微積分的計(jì)算可以說(shuō)是數(shù)學(xué)的基礎(chǔ)，關(guān)于基礎(chǔ)的微機(jī)分計(jì)算主要有求函數(shù)的極限，微分運(yùn) 算，積分運(yùn)算，級(jí)數(shù)的計(jì)算，代數(shù)的求解，以及常微分方程的求解，下面分別介紹 。 極限的基本運(yùn)算 MATLAB 中極限函數(shù) limit 格式如下： limit(F,x,a) ； limit(F,a)； limit(F)； limit(F,x,a,39。right39。)；limit(F,x,a,39。left39。) 其中 F表示的是函數(shù)式， x 表示求極限的變量， a表示的是變量取的值， left 或right 表示是取左極限還是右極限。 例 ： 求2201)ln(lim yx ex yyx ????； 分析：這是一個(gè)典型的求極限的題目，題目中涉及了兩個(gè)變 量，故在程序開(kāi)始應(yīng)該先定義兩個(gè)變量，然后直接利用函數(shù) limit 編程計(jì)算。 示例程序如下： syms x y limit(limit(log(x+exp(y))/sqrt(x*x+y*y),x,1),y,0) 運(yùn)行結(jié)果如下： 武漢理工大學(xué)《 Matlab 課程設(shè)計(jì)》說(shuō)明書(shū) 課件之家精心整 理資料 歡迎你的欣賞 微分的計(jì)算 MATLAB 中微分函數(shù) diff 格式如下： Y = diff(X)； Y = diff(X,n)； Y = diff(X,n,dim) 其中 X表示待微分的變量， n表示 n 次微分，第三式表示沿著定維 dim 的 n階微分。 例 ： 求 ?y 23 sin(4x+2)的微分 。 分析：微分運(yùn)算能直觀的顯示函數(shù)值的變化快慢，在平時(shí)計(jì)算中常常反映變量的變化對(duì)函數(shù)值的影響大小這是一個(gè)簡(jiǎn)單的一階微分式子 ,可以直接寫(xiě)定義變量然后求結(jié)果。 示例程序如下： syms x y dy y=sqrt(3)/2*sin(4*x+2) dy=diff(y) 運(yùn)行結(jié)果如下： 積分的計(jì)算 MATLAB 中積分常用函數(shù)為 int，格式如下： int(f,x,a,b) 其中 f 表示待積分的函數(shù)， x 表 示積分變量，而 a， b 則分別表示積分起始終止點(diǎn)。 例： 求下列變上限積分 ? ?2 1xxdxx 分析：積分運(yùn)算大量運(yùn)用于求面積體積等，此處選用的是一個(gè)變上限積分，武漢理工大學(xué)《 Matlab 課程設(shè)計(jì)》說(shuō)明書(shū) 課件之家精心整 理資料 歡迎你的欣賞 屬于有些典型的積分例子，但任然是根據(jù)函數(shù)格式就可以寫(xiě)出程序。 示例程序如下： syms k x y=(1+k)^(1/2) f=int(y,k,x,x^2) 運(yùn)行結(jié)果如下： 級(jí)數(shù)的計(jì)算 MATLAB 中級(jí)數(shù)常用函數(shù)為 symsum，格式如下： r = symsum(s)； r = symsum(s,v)； r = symsum(s,a,b)； r = symsum(s,v,a,b) 函數(shù)表達(dá)的意義是表達(dá)式 s 關(guān)于變量 v從 a到 b求和。 例： 計(jì)算級(jí)數(shù) .11 2???n n 的值 , 分析：對(duì)于級(jí)數(shù)我們最先接觸的就是常數(shù)項(xiàng)無(wú)窮級(jí)數(shù)，因?yàn)榧夹g(shù)結(jié)果一般較復(fù)雜，所以用 simple 求解最簡(jiǎn)形式 , 其中還用到 Inf 為無(wú)窮大。 示例程序如下： syms k simple(symsum(1/k^2,1,Inf)) %simple 求解最簡(jiǎn)形式 ,Inf 為無(wú)窮大 運(yùn)行結(jié)果如下： 求解代數(shù)微分方程 代數(shù)方程分為很多種，有簡(jiǎn)單有復(fù)雜 ，方法也有很多。一般多項(xiàng)式方程的根武漢理工大學(xué)《 Matlab 課程設(shè)計(jì)》說(shuō)明書(shū) 課件之家精心整 理資料 歡迎你的欣賞 可為實(shí)數(shù)，也可為復(fù)數(shù)， 可用 MATLAB 符號(hào)工具箱中的 solve( )函數(shù) ， MATLAB中求解代數(shù)方程常用函數(shù) solve 格式如下： [x,? ]=solve(eqn1,eqn2,? ,eqnn，’ x,?’ ) 其中 eqn 表示的是式子， x 等表示的是變量。 例 ：求解方程組 3x+2=20 分析：為了計(jì)算的簡(jiǎn)單，同時(shí)為了避免需要驗(yàn)證的麻煩，這里選用了基礎(chǔ)的一元一次方程，直接用語(yǔ)句描述即可。 示例程序如下： Syms x [x]=solve(‘3*x+2=20’ ,’x’) 運(yùn)行結(jié)果如下 ： 求解常微分方程 未知的函數(shù)以及它的某些階的導(dǎo)數(shù)連同自變量都由一已知方程聯(lián)系在一起的方程稱(chēng)為微分方程。如果未知函數(shù)是一元函數(shù)，稱(chēng)為常微分方程 。 MATLAB 中主要用 dsolve 求符號(hào)解析解 。 dsolve 的基本格式是： s=dsolve(‘方程 1’ , ‘方程 2’ ,? ,’初始條件 1’，’初始條件 2’ ? ,’自變量’ ) 例 ： 求解微分方程 ,1)0(,139。 ????? ytyy 求解析解 。 分析：求解此微分方程的解析解直接采用 dsolve 即可。 Dsolve 示例程序如下 : 武漢理工大學(xué)《 Matlab 課程設(shè)計(jì)》說(shuō)明書(shū) 課件之家精心整 理資料 歡迎你的欣賞 syms y t %定義變量 s=dsolve(39。dy=y+t+139。,39。y(0)=139。,39。t39。) %程序主體部分 simplify(s)%化簡(jiǎn)函數(shù) 運(yùn)行結(jié)果如下： 矩陣的基本運(yùn)算 矩陣計(jì)算是線性代數(shù)中的核心內(nèi)容，其對(duì)于整個(gè)數(shù)學(xué)系統(tǒng)的計(jì)算方面的意義是十分巨大的，集中它的基本運(yùn)算包括最大值、最小值、均值、方差、轉(zhuǎn)置、逆、行列式、特征值的計(jì)算、矩陣的相乘、右除、左除、冪運(yùn)算等等，下面將具體介紹。 矩陣的運(yùn)算都是要以矩陣為基礎(chǔ)的，本報(bào)告中決定選用一組矩陣來(lái)完成幾乎全部可以完成的計(jì)算，那么首先就得生成矩陣了。 矩陣的定義和分配可以 有多種方法。最簡(jiǎn) 單的方法是有方括號(hào) []包圍的逐行給定元素。若定義一個(gè)標(biāo)量，則方括號(hào)就不需要了。 相同行中的元素是由一行或多個(gè)空格 ‘’ 或一個(gè)逗號(hào) ‘ ， ’ 分隔，列由分號(hào) ‘。 ’ 或回車(chē)鍵分隔。沒(méi)有結(jié)尾分號(hào)的每個(gè)命令在屏幕上顯示出其結(jié)果。若結(jié)尾帶分號(hào)，就執(zhí)行計(jì)算，但計(jì)算結(jié)果并不顯示。 如定義 3 3矩陣如下： ???????????987654321A ?????????? ???13519661291B 則在命令窗口輸入： A=[1 2 3。4 5 6。 7 8 9] 武漢理工大學(xué)《 Matlab 課程設(shè)計(jì)》說(shuō)明書(shū) 課件之家精心整 理資料 歡迎你的欣賞 屏幕顯示結(jié)果為： 同上依次輸入： B=[1 9 12。6 6 9。1 5 13] 就可以得到 B矩陣（如下）。 矩陣的最大最小值 MATLAB 中 max 函數(shù)可以表示求每一列的最大值，那么經(jīng)過(guò)分析可以知道，先求出每一列的最大值然后求出這些最大值里面的最大值，下面以 A 矩陣為例。 示例程序如下： y=max(A), x=max(y) 運(yùn)行結(jié)果如下： 對(duì)比 A中數(shù)值發(fā)現(xiàn)結(jié)果是正確的。 MATLAB 中求最小值的函數(shù)為 min，求解思路與求最大值思路類(lèi)似，仍然以矩武漢理工大學(xué)《 Matlab 課程設(shè)計(jì)》說(shuō)明書(shū) 課件之家精心整 理資料 歡迎你的欣賞 陣 A 為例。 示例程序如下： y=min(A), x=min(y) 運(yùn)行結(jié)果如下： 對(duì)比 A中 數(shù)值發(fā)現(xiàn)結(jié)果是正確的。 矩陣的均值方差 MATLAB 中求解矩陣均值的函數(shù)是 mean，它的具體用法如下： mean(A,1)表示對(duì)列取平均， mean(A,2）表示對(duì)行取平均， mean(A）則默認(rèn)為 mean(A,1)。 下面以矩陣 B 分別舉例，程序示例如下： a=mean(B),b=mean(B,2) 運(yùn)行結(jié)果如下： 武漢理工大學(xué)《 Matlab 課程設(shè)計(jì)》說(shuō)明書(shū) 課件之家精心整 理資料 歡迎你的欣賞 觀察可知， a， b分別顯示出了矩陣行列的均值。如果想求矩陣的均值可以進(jìn)行 2次操作。 示例程序如下： c=mean（ a） 運(yùn)行結(jié)果如下： 可以觀察到 c 的值就是矩陣 b所有值的均值。 MATLAB 中求解矩陣方差的函數(shù)是 var，它的常用格式是 V = var(X)， 如果 X是一個(gè)矩陣， var(X)返回一個(gè)包含矩陣 X 每一列方差的行向量。 下面還是以矩陣 B 來(lái)示例，程序如下： var(B) d=var(var(B)) 運(yùn)行結(jié)果如下： 矩陣的轉(zhuǎn)置 矩陣的一個(gè)重要的運(yùn)算是轉(zhuǎn)置，如果 A是一個(gè)實(shí)數(shù)矩陣，那么它被轉(zhuǎn)置時(shí)，第 1 行變成第 1 列，第 2 行變成第 2 列，依此類(lèi)推，一個(gè) m n 矩陣變?yōu)橐粋€(gè) n m 矩陣。如果矩陣是方陣，那么這個(gè)矩陣在主對(duì)角線反映出來(lái)。 MAT