Matlab多項式運算.ppt

1、Matlab多項式運算與方程求根,Matlab多項式運算,在 Matlab 中，n 次多項式是用一個長度為 n+1的向量來表示，缺少的冪次項系數(shù)為 0。例如：,在 Matlab中表示為相應(yīng)的向量：,例：,注：系數(shù)中的零不能省！,多項式四則運算,多項式加減運算：Matlab沒有提供專門進行多項式加減運算的函數(shù)，事實上，多項式的加減就是其所對應(yīng)的系數(shù)向量的加減運算。,例：,對于次數(shù)相同的多項式，可以直接對其系數(shù)向量進行加減運算； 如果兩個多項式次數(shù)不同，則應(yīng)該把低次多項式中系數(shù)不足的高次項用0補足，然后進行加減運算。,多項式四則運算（續(xù)）,多項式乘法運算： k=conv(p,q),例：計算多項式

2、和 的乘積, p=2,-1,0,3; q=2,1; k=conv(p,q);,多項式除法運算：k,r=deconv(p,q),其中k返回的是多項式p除以 q的商，r是余式。,k,r=deconv(p,q),p=conv(q,k)+r,多項式的導(dǎo)數(shù)：polyder,k=polyder(p):多項式p的導(dǎo)數(shù)； k=polyder(p,q): p*q 的導(dǎo)數(shù)； k,d=polyder(p,q):p/q 的導(dǎo)數(shù)，k是分子，d是分母。, k1=polyder(2,-1,0,3); k2=polyder(2,-1,0,3,2,1); k2,d=polyder(2,-1,0,3,2,1);,例：已知 ， ，

3、 求,多項式求值, p=2,-1,0,3; x=2;polyval(p,x) x=-1, 2;-2,1;polyval(p,x),例：已知 ，分別取 x=2和一個2x2矩陣， 求 p(x)在 x處的值,代數(shù)多項式求值：,y=polyval(p,x):計算多項式p在x點的值,注：若 x 是向量或矩陣，則采用數(shù)組運算（點運算）！,多項式求值（續(xù)）, p=2,-1,0,3; x=-1, 2;-2,1;polyval(p,x) polyvalm(p,x),例：已知 ，則,矩陣多項式求值：,Y=polyvalm(p,X)：以方陣X為自變量， 計算多項式的值，采用矩陣運算。,polyvalm(p,A)=2

4、*A*A*A - A*A + 3*eye(size(A); polyval(P,A)=2*A.*A.*A - A.*A + 3*ones(size(A),多項式求根, p=2,-1,0,3; x=roots(p),例：已知 ，求p(x)的零點。,x=roots(p)：若p是n次多項式，則輸出x為包含p=0的n個根的n維向量。,若已知多項式的全部零點，則可用poly函數(shù)給出該多項式。,p=ploy(x),注：以上多項式運算中，使用的都是多項式 的 系數(shù)向量，不涉及符號計算！,Matlab非線性方程的數(shù)值求解,fzero(f,x0)：求方程f=0在x0附近的根。,(1)方程可能有多個根，但fzer

5、o之給出離x0最近的一個根；,(2)若x0是一個標(biāo)量，則fzero先找出一個包含x0的區(qū)間，使得f在這個區(qū)間兩個端點上的值異號，然后再在這個區(qū)間內(nèi)尋找方程f=0的根；如果找不到這樣的區(qū)間，則返回 NaN。,幾點說明：,(4)由于fzero是根據(jù)函數(shù)是否穿越橫軸來決定零點，因此它無法確定函數(shù)曲線僅觸及橫軸但不穿越的零點，如|sin(x)|的所有零點。,(3)若x0是一個2維向量，則表示在x0(1),x0(2)區(qū)間內(nèi)求方程的根，此時必須滿足f在這兩個端點上的值異號。,(5)函數(shù)中的f是一個函數(shù)句柄，可通過一下方式給出：,字符串形式：fzero(x3-3*x+1,2); 通過調(diào)用的函數(shù)句柄：fzer

6、o(sin,4);,(6) f不能用符號表達(dá)式！,例：,fzero(sin(x),10) fzero(sin,10),fzero(x3-3*x+1,1) fzero(x3-3*x+1,1,2) fzero(x3-3*x+1,-2,0),Matlab符號方程求解,s=solve(f,v)：求方程關(guān)于指定自變量的解； s=solve(f)：求方程關(guān)于默認(rèn)自變量的解。,其中 f可以是用字符串表示的方程，或符號表達(dá)式；若f 中不含等號，則表示解方程 f=0。,例：解方程 x3-3*x+1=0,syms x; f=x3-3*x+1; s=solve(f,x),s=solve(x3-3*x+1,x),s=

7、solve(x3-3*x+1=0,x),solve也可以用來解方程組,solve(f1,f2,.,fN,v1,v2,.,vN) 求解由 f1,f2,.,fN 確定的方程組關(guān)于 v1, v2,.,vN 的解。,例：解方程組,x,y,z=solve(x+2*y-z=27,x+z=3, . x2+3*y2=28,x,y,z),輸出變量的順序要書寫正確！,solve在得不到解析解時，會給出數(shù)值解。,線性方程組求解,linsolve(A,b)：解線性方程組,例：解方程組, A=1 2 1; 1 0 1; 1 3 0; b=2;3;8; X=linsolve(A,b),b是列向量！,求解方程函數(shù)小結(jié),roots(p)：多項式的所有零點，p是多項式系數(shù)向量。,fzero(f,x0)：求f=0在x0附近的根，f是函數(shù)句柄，可以由字符串給出或使用，但不能是符