ch2、matlab語(yǔ)言基礎(chǔ)-在線瀏覽
2024-11-04 14:58本頁(yè)面
【正文】 .39。 b = + 》 b=[1+2i 27i].39。3 4]； b=[ 3 5。 29 51] 》 a/b=[ 。2 3] 》 a^3=[37 54。15 36] 》 a./b=[ 。 ] 》 a.^3= [1 8。 ? 注意只有當(dāng)兩個(gè)矩陣中前一個(gè)矩陣的列數(shù)和后一個(gè)矩陣的行數(shù)相同時(shí)，才可以進(jìn)行乘法運(yùn)算。只有方陣才可以求冪。 逆矩陣與行列式計(jì)算 求逆： inv(A)； 求行列式： det(A) 要求矩陣必須為方陣 了解矩陣超越函數(shù) ?在 MATLAB中 exp、 sqrt等命令也可以作用到矩陣上，但這種運(yùn)算是定義在矩陣的單個(gè)元素上的，即分別對(duì)矩陣的每一個(gè)元素進(jìn)行計(jì)算。矩陣的超越函數(shù)要求運(yùn)算矩陣為方陣。 4 5 6。 》 b=inv(a) b = 》 det(a) ans = 3 ? A(m,n)：提取第 m行，第 n列元素 ? A(:,n)：提取第 n列元素 ? A(m,:)：提取第 m行元素 ? A(m1:m2,n1:n2)：提取第 m1行到第 m2行和第 n1列到 第 n2列的所有元素（提取子塊）。 ? 矩陣擴(kuò)展：如果在原矩陣中一個(gè)不存在的地址位置上設(shè)定一個(gè)數(shù)（賦值），則該矩陣會(huì)自動(dòng)擴(kuò)展行列數(shù)，并在該位置上添加這個(gè)數(shù)，而且在其他沒(méi)有指定的位置補(bǔ)零。 三、矩陣的操作 矩陣下標(biāo) ?MATLAB通過(guò)確認(rèn)矩陣下標(biāo)，可以對(duì)矩陣進(jìn)行插入子塊，提取子塊和重排子塊的操作。 ? length(A)=max(size(A))：返回行數(shù)或列數(shù)的最大值。3 4 5]。若已知多項(xiàng)式等于 0的根，函數(shù) poly可以求出相應(yīng)多項(xiàng)式。 b=[1 2] c=conv(a,b)=1 4 7 6 conv指令可以嵌套使用，如 conv(conv(a,b),c) ?相除 deconv [q,r]=deconv(c,b) q=1 2 3 ％商多項(xiàng)式 r=0 0 0 ％余多項(xiàng)式 ?求多項(xiàng)式的微分多項(xiàng)式 polyder polyder(a)=2 2 ?求多項(xiàng)式函數(shù)值 polyval(p,n)：將值 n代入多項(xiàng)式求解。這在分析實(shí)驗(yàn)數(shù)據(jù)，將實(shí)驗(yàn)數(shù)據(jù)做解析描述時(shí)非常有用。 ? 例 （ 4） *多項(xiàng)式插值 ? 多項(xiàng)式插值是指根據(jù)給定的有限個(gè)樣本點(diǎn)，產(chǎn)生另外的估計(jì)點(diǎn)以達(dá)到數(shù)據(jù)更為平滑的效果。 ? 所用指令有一維的 interp二維的 interp三維的interp3。 Help polyfun可以得到更詳細(xì)的內(nèi)容。 [u,s,v]=svd(a) u = s = 0 0 v = （ 2）特征值分解 [V,D]=eig(A) 例： a = 9 8 6 8 [v,d]=eig(a) v = d = 0 0 求矩陣 A的特征向量 V及特征值D，滿足 A*V=V*D。如果 D=eig(A)則只返回特征值。 norm(q) ans = 1 q*r ans = （ 4）三角分解 [L,U]=lu(A) ? 將 A做對(duì)角線分解，使得 A=L*U,其中 L為下三角矩陣， U為上三角矩陣。 [L1,U1,P]=lu(A) 例： a=[1 2 3。7 8 9] 比較： [l1,u1,p]=lu(a) [l,u]=lu(a) l1 = 0 0 0 u1 = 0 0 0 p = 0 0 1 1 0 0 0 1 0 l = 0 0 0 u = 0 0 0 可以驗(yàn)證： u1=u， inv(p)*l1=l a=l*u p*a=l1*u1 2*、數(shù)據(jù)分析 （ 1） 繪制函數(shù)圖形： fplot() （ 2） 求極值： fmin,fmins （ 3） 求零點(diǎn)：尋找一維函數(shù)的過(guò)零點(diǎn) fzero() （ 4） 頻譜分析 （ fft） ： y=FFT(x)； unwrap()； abs；angle畫(huà)出幅頻和相頻曲線 （ 5） 了解數(shù)據(jù)分析函數(shù)： max,min,mean,sum,prod等 （ 6） 了解積分運(yùn)算： trap2,quad,quad8 3*、 常微分方程數(shù)值解 [t,x]=ode23(‘xfun’,t0,tf,x0,tol) [t,x]=ode45(‘xfun’,t0,tf,x0,tol) 第三節(jié)、繪圖簡(jiǎn)介 ? MATLAB提供了豐富的繪圖功能 help graph2d可得到所有畫(huà)二維圖形的命令 help graph3d可得到所有畫(huà)三維圖形的命令 下面介紹常用的二維圖形命令 基本的繪圖命令 plot（ x1,y1,option1,x2,y2,option2,… ） x1,y1給出的數(shù)據(jù)分別為 x,y軸坐標(biāo)值， option1為選項(xiàng)參數(shù)，以逐點(diǎn)連折線的方式繪制 1個(gè)二維圖形；同時(shí)類似地繪制第二個(gè)二維圖形， ……等。比如： plot(x,y)； plot(x,y,option) 選項(xiàng)參數(shù) option定義了圖形曲線的顏色、線型及標(biāo)示符號(hào)，它由一對(duì)單引號(hào)括起來(lái)。 grid on：在所畫(huà)出的圖形坐標(biāo)中加入柵格 grid off：除去圖形坐標(biāo)中的柵格 hold on：把當(dāng)前圖形保持在屏幕上不變，同時(shí) 允許在這個(gè)坐標(biāo)內(nèi)繪制另外一個(gè)圖形。 ? gtext(‘字符串’ ) 利用鼠標(biāo)在圖形的某一位置標(biāo)示字符串。 ? xlabel(‘字符串’ )， ylabel(‘字符串’ ) 設(shè)置 x， y坐標(biāo)軸的名稱。 legend(‘字符串 1’,‘字符串 2’,…,‘ 字符串 n’) ? 在屏幕上開(kāi)啟一個(gè)小視窗，然后依據(jù)繪圖命令的先后次序，用對(duì)應(yīng)的字符串區(qū)分圖形上的線。 semilogx：繪制以 x軸為對(duì)數(shù)坐標(biāo)（以 10為底）， y軸為線性坐標(biāo)的半對(duì)數(shù)坐標(biāo)圖形。 例 subplot（ mnk）： 分割圖形顯示窗口 m:上下分割個(gè)數(shù)， n:左右分割個(gè)數(shù)， k:子圖編號(hào) 例 了解應(yīng)用型繪圖指令：可用于數(shù)值統(tǒng)計(jì)分析或離散
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1