[計(jì)算機(jī)軟件及應(yīng)用]第四講 數(shù)值計(jì)算基礎(chǔ)-文庫(kù)吧
2024-10-04 04:17 本頁(yè)面
【正文】 [2 5 。7 9]及標(biāo)量 5的值 。 ? 在命令窗口輸入多項(xiàng)式系數(shù) ， 并求值： ? p = [3 2 1]。 ? polyvalm(p,[2 5。7 9]) ? polyvalm(p,5) 123 2 ?? xx單擊此處編輯母版標(biāo)題樣式 ? 2. 求多項(xiàng)式的根 ? MATLAB進(jìn)行多項(xiàng)式的求根運(yùn)算時(shí)，可以有兩種方法， 一種是直接調(diào)用求根函數(shù) roots。另一種是先把多項(xiàng)式轉(zhuǎn)化為伴隨矩陣， 然后再求其特征值。 在 matlab中約定，多項(xiàng)式系數(shù)用行矢量表示，一組根用列矢量表示。 單擊此處編輯母版標(biāo)題樣式 ? 【 例 36】 分別用兩種方法求多項(xiàng)式 ? 的根 。 ? 在命令窗口輸入多項(xiàng)式系數(shù) ， 并求根： ? a=[1 5 3 6 4 10]。 ? r=roots(a) ? s=pan(a) ? r=eig(s) 104635 2345 ????? xxxxx單擊此處編輯母版標(biāo)題樣式 ? 3. 多項(xiàng)式的乘除運(yùn)算 ? 多項(xiàng)式的乘法由函數(shù) conv實(shí)現(xiàn)， conv也是矢量的卷積函數(shù)：多項(xiàng)式的除法由函數(shù) deconv實(shí)現(xiàn)， deconv也是矢量的卷積函數(shù)的逆函數(shù)。 單擊此處編輯母版標(biāo)題樣式 ? 【 例 37】 計(jì)算兩多項(xiàng)式 和 的乘法。 ? 在命令窗口輸入多項(xiàng)式系數(shù) ， 求兩多項(xiàng)式的乘法： ? a=[1 5 3 4 2]。 ? b=[1 2 5 3]。 ? c=conv(a,b) 2435 234 ???? xxxx352 23 ??? xxx單擊此處編輯母版標(biāo)題樣式 ? 【 例 38】 計(jì)算例 37中求得的乘積被? 除所得結(jié)果 。 ? 利用例 37中求得的結(jié)果 ， 求兩多項(xiàng)式的除法： ? c=[1 3 12 30 36 33 22 6]。 ? b=[1 2 5 3]。 ? d=deconv(c,b) 352 23 ??? xxx單擊此處編輯母版標(biāo)題樣式 ? 4. 多項(xiàng)式的微積分 ? 多項(xiàng)式的微分由函數(shù) polyder實(shí)現(xiàn) ， 多項(xiàng)式的積分由函數(shù) polyint實(shí)現(xiàn) 。 ? 【 例 39】 計(jì)算多項(xiàng)式 的微分。 ? 在命令窗口輸入多項(xiàng)式系數(shù) ， 求多項(xiàng)式的微分 ： ? p=[3 5 2 6 10]。 ? polyder(p) ? poly2sym(ans) 106253 234 ???? xxxx單擊此處編輯母版標(biāo)題樣式 ? 【 例 310】 計(jì)算多項(xiàng)式 的積分 。 ? 在命令窗口輸入多項(xiàng)式系數(shù) ， 求多項(xiàng)式的積分： ? p=[12 15 4 6]。 ? polyint(p) 641512 23 ???? xxx單擊此處編輯母版標(biāo)題樣式 ? 5. 多項(xiàng)式的部分分式展開(kāi) ? 對(duì)于多項(xiàng)式 b(x)和不含重根的 n階多項(xiàng)式 a(x)之比 ， 有如下的展開(kāi)： ? 式中 稱為極點(diǎn)（ Ploes）， 稱為留數(shù)（ Residues）， k(x)稱為直項(xiàng)（ Diret term）。假如 a(x)有 m重根p(j)=… .=p(j+m1), 則相應(yīng)的部分寫(xiě)成： )(...)( )(22`11 xkpxrpxrpxrxaxbnn ????????nppp ,..., 21 nrrr ,..., 21? ? ? ? mjmjjjjj pxrpxrpxr????????? 121 . . .單擊此處編輯母版標(biāo)題樣式 ? 在 MATLAB中 ， 計(jì)算部分分式展開(kāi)的函數(shù)為 residue， 其調(diào)用格式為： ? [r,p,k]=residue(d,a) 計(jì)算由多項(xiàng)式表達(dá)式 b(x) 和 a(x)分別為分子和分母時(shí)所得到的留數(shù) r、 極點(diǎn) p和直項(xiàng) k。 ? [b,a]=residue(r,p,k)將部分分式展開(kāi)轉(zhuǎn)換回多項(xiàng)式表達(dá)式 b(x)和 a(x)的系數(shù) 。 單擊此處編輯母版標(biāo)題樣式 ? 【 例 311】 兩多項(xiàng)式的比為 ，求部分分式展開(kāi)。 ? 在命令窗口輸入多項(xiàng)式系數(shù) ， 求部分分式展開(kāi)： ? a = [4 0 8 3]。 ? b = [ 5 3 2 7]。 ? [r, p, k] = residue(b,a) 3847235)()(323???????xxxxxxaxb單擊此處編輯母版標(biāo)題樣式 ? 【 例 312】 利用在例 311中求得的部分分式展開(kāi)結(jié)果轉(zhuǎn)換回多項(xiàng)式表達(dá)式 b(x)和 a(x)的系數(shù) 。 ? [b,a] = residue(r,p,k) ? 利用在例 311中求得的結(jié)果 ， 求多項(xiàng)式系數(shù)： ? 得到的結(jié)果為 ? 這是對(duì)于表達(dá)式分母的歸一形式 。 )()(323???????xxxxxxaxb單擊此處編輯母版標(biāo)題樣式 ? 6. 多項(xiàng)式擬合 ? 多項(xiàng)式擬合的函數(shù)為 polyfit,其調(diào)用格式為： ? p=polyfit(x,y,n) 應(yīng)用最小二乘法求出 n階擬合多項(xiàng)式 p(x)， 即用 p(x(i))擬合 y(i). ? [p,S]=polyfit(x,y,n) 返回 n階擬合多項(xiàng)式p(x)合包括誤差估計(jì)的結(jié)構(gòu) S。 ? [p,S， mu]=polyfit(x,y,n) 在這
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1