電-Matlab基礎及應用3-091130

1、第3 章 Matlab數(shù)值計算 多項式 插值和擬合 數(shù)值微積分 線性方程組的數(shù)值解 稀疏矩陣 常微分方程的數(shù)值解3.1 多項式1、多項式的建立 Matlab 采用列向量表示多項式，將多項式的系數(shù)按降冪次序存放在列向量中。2、多項式的四則運算(1).多項式的加減運算： 相同次數(shù)的，直接相加減 不同次數(shù)的，在較低次冪的多項式系數(shù)前補0(2).多項式乘法運算函數(shù)conv(p1,p2)用于求多項式p1和p2的乘積。這里，p1、p2是兩個多項式系數(shù)向量。(3).多項式除法函數(shù)Q,r=deconv(p1,p2)用于對多項式p1和p2作除法運算。其中Q返回多項式p1除以p2的商式，r返回p1除以p2的余式。

2、這里，Q和r仍是多項式系數(shù)向量。deconv是conv的逆函數(shù)，即有p1=conv(p2,Q)+r。3、多項式求值函數(shù)：polyval(p,x)p：多項式的系數(shù)向量x：為一個向量或者矩陣，要求解的x點4、多項式求根函數(shù)：r=roots(p)p:多項式的系數(shù)向量r:方程p(x)=0的全部根(重根、復根)；r(1),r(2),r(n)分別代表多項式的n個根。由給定的根創(chuàng)建多項式：p=poly(r)r:方程的根p:多項式的系數(shù)向量5、多項式的求導函數(shù)是p=polyder(a)p=polyder(a)：求多項式a的導函數(shù)p=polyder(a,b)：求ab的導函數(shù)p,q=polyder(a,b)：求a

3、/b的導函數(shù)，導函數(shù)的分子存入p，分母存入q。上述函數(shù)中，參數(shù)a,b是多項式的向量表示，結果p,q也是多項式的向量表示。3.2 插值和擬合 如果測量值是準確的，沒有誤差，一般用插值；如果測量值與真實值有誤差，一般用曲線擬合。(1) 多項式的插值和擬合數(shù)值插值的原理：參考數(shù)值計算方法的有關內容1、一維數(shù)值插值：被插值函數(shù)是一個單變量函數(shù)。 Y1=interp1(X,Y,X1,method) 根據(jù)X,Y的值，計算函數(shù)在X1處的值。X,Y：分別描述采樣點X和采樣點函數(shù)值YX1： 欲插值點的向量或標量;method：插值方法：linear(默認),nearest,spline,cubic2、二維數(shù)值插

4、值：Z1=interp2(X,Y,Z,X1,Y1,method)3、三維數(shù)值插值：V1=interp3(X,Y,Z,V,X1,Y1,Z1,method)注意：X1的取值范圍不能超出X的給定范圍，否則，會給出“NaN”錯誤。MATLAB中有一個專門的3次樣條插值函數(shù)Y1=spline(X,Y,X1)，其功能及使用方法與函數(shù)Y1=interp1(X,Y,X1,spline)完全相同。曲線擬合：目的是根據(jù)給定區(qū)間或者區(qū)域上的有限個采樣點的函數(shù)值，構造一個較簡單的函數(shù)去逼近一個復雜或者未知的函數(shù)。一般可以采用最小二乘法，實際上是求一個多項式的系數(shù)向量，通過多項式實現(xiàn)擬合。p,s=polyfit(x,y

5、,m)：根據(jù)采樣點x和采樣點函數(shù)值y，產生一個m次多項式系數(shù)向量p及其在采樣點上的誤差sy=polyval(p,x)：求擬合多項式函數(shù)在采樣點x處的函數(shù)值，其中p為多項式的系數(shù)向量。3.3 數(shù)值微積分數(shù)值微分在MATLAB中，沒有直接提供求數(shù)值導數(shù)的函數(shù)，只有計算向前差分的函數(shù)diff，其調用格式為：dx=diff(X,n):計算X的n階向前差分。例如，diff(X,2)=diff(diff(X)。dx=diff(A,n,dim):計算矩陣A的n階差分，dim=1時(缺省狀態(tài))，按列計算差分；dim=2，按行計算差分。dx=diff(X):計算向量X的向前差分，x(i)=x(i+1)-x(i)

6、，i=1,2,n-1。數(shù)值積分1. 數(shù)值積分基本原理求解定積分的數(shù)值方法多種多樣，如簡單的梯形法、辛普生(Simpson)法、牛頓柯特斯(Newton-Cotes)法等都是經常采用的方法。它們的基本思想是無限元法：將整個積分區(qū)間a,b分成n個子區(qū)間xi,xi+1，i=1,2,n，其中x1=a，xn+1=b。這樣求定積分問題就分解為求和問題。f： 被積函數(shù)a,b： 積分區(qū)間a,btol： 計算精度默認為0.001trace：非零時畫出積分圖形2. 求解數(shù)值積分的函數(shù)：quad(f,a,b,tol,trace)注意：調用quad函數(shù)時，要先建立一個描述被積函數(shù)f的函數(shù)文件或者語句函數(shù)用兩種不同的方

7、法求定積分先建立一個函數(shù)文件ex.m：function ex=ex(x)ex=exp(-x.2);I=quad(ex,0,1)I = 0.7468I=quad (ex,0,1)I = 0.7468然后在MATLAB命令窗口，輸入命令：也可不建立關于被積函數(shù)的函數(shù)文件，而使用語句函數(shù)(內聯(lián)函數(shù)inline)求解，命令如下：g=inline(exp(-x.2); %定義一個語句函數(shù)I=quadl(g,0,1) %注意函數(shù)名不加號I = 0.74683.二重積分數(shù)值求解 使用MATLAB提供的dblquad函數(shù)就可以直接求出上述二重定積分的數(shù)值解。該函數(shù)的調用格式為： I=dblquad(f,a,b

8、,c,d,tol,trace) 該函數(shù)求f(x,y)在a,bc,d區(qū)域上的二重定積分。參數(shù)tol，trace的用法與函數(shù)quad完全相同。3.5 稀疏矩陣稀疏矩陣的建立稀疏矩陣創(chuàng)建函數(shù)：sparseS=sparse(A):可以將一個全元素存儲矩陣A轉換為稀疏存儲矩陣S A=full(S):將一個稀疏存儲矩陣S轉換為全元素存儲矩陣AS=sparse(i,j,s,m,n):直接創(chuàng)建稀疏矩陣。 i和j： 矩陣非零元素的行和列的下標向量 s： 對應行列處的非零元素值向量 m和n：矩陣的行數(shù)和列數(shù)稀疏矩陣的查看和運算nnz(s):返回稀疏矩陣非零元素的個數(shù)nonzeros(s):返回稀疏矩陣中所有非零元

9、素的值nzmax(s):返回分配給稀疏矩陣中非零元素的總存儲空間i,j,k=find(s):返回非零元素的下標和數(shù)值 i：非零元素的行下標向量 j：非零元素的列下標向量 k：非零元素的數(shù)值第4章 結構數(shù)組與細胞數(shù)組 結構數(shù)組的創(chuàng)建 結構數(shù)組的操作 細胞數(shù)組的創(chuàng)建 細胞數(shù)組的操作4.1 結構數(shù)組結構數(shù)組的格式:.結構數(shù)組的創(chuàng)建賦值法 struct_name(record#).field_name=data 創(chuàng)建1*1結構數(shù)組時可省略記錄號。 1*2的結構數(shù)組，并且當結構數(shù)組包含兩個以上的元素時，查看結構數(shù)組不顯示各個元素的值，而是顯示數(shù)組的結構信息。利用函數(shù)struct創(chuàng)建 函數(shù)基本形式為： s

10、truct_name = struct(field1,field2, .) struct_name = struct(field1,val1,field2,val2, .) 注意：續(xù)行符“” 如果域沒有值，則一定要賦空值，不能空著。 多個元素域值相同，可以只賦一次值。結構數(shù)組的操作函數(shù)名說明struct創(chuàng)建結構數(shù)組isstruct判定是否為結構數(shù)組fieldnames獲取結構數(shù)組域名setfield設定域值getfield獲取域值isfield判斷域是否在結構數(shù)組中rmfield刪除結構數(shù)組中的域orderfields域排序4.2 細胞數(shù)組細胞數(shù)組的創(chuàng)建通過賦值 cell_namei,j=value利用函數(shù)cell創(chuàng)建細胞數(shù)組 cell_name=cell(n) cell_name=cell(m,n)或cell_name=cell(m n) cell_name=cell(m,n,p,)或cell_name=cell(m n p) cell_name=cell(size(A)細胞數(shù)組的操作函數(shù)名說明celldisp顯示細胞數(shù)組所有元素的內容iscell判定是否為細胞數(shù)組is