多项式的处理与分析.PDFVIP

多項式的處理與分析 本章介紹 MATLAB 處理多項式的相關指令，以便用於多項式的計算與分析， 包括：多項式的加、減、乘、除、求值、求根、微分、積分、矩陣的特徵多 項式、部份分式展開、多項式擬合等。 MATLAB 程式設計進階篇 7-1 多項式的加減乘除 一個 次的多項式可以表示成 n p (x ) a x n a x n1 a x a n n1 1 0 因此在 MATLAB 中，可以用一個長度為 n+1 的列向量來表示多項式p (x) 如 下： p [a ,a ,,a ,a ] n n1 1 0 舉例來說，我們可用 來表示一個三次多項式 3 2 p = [1, 2, 3, 1] p (x) = x + 2x + 3x + 1 。 多項式的加減，可直接由向量的加減而推出。若有兩多項式分別是p 1(x) = x 3 + x + 1 及p 2 (x) = x2 −x + 2 ，則其和與差可計算如下： Example 1: 07-多項式的處理與分析/polyPlus.m p1 = [1,0,1,1]; p2 = [0,1,-1,2]; p1 + p2 p1 - p2 ans = 1 1 0 3 ans = 1 -1 2 -1 必須注意的是：矩陣 與 的長度要一致，否則 就會產生運算 p1 p2 MATLAB 錯誤的訊息。 多項式的乘與除，可使用 conv 及 deconv 指令來達成，例如，欲求多項式p 1(x ) 與p 2 (x ) 的乘積，可輸入如下： 7-2 第 7 章 多項式的處理與分析 Example 2: 07-多項式的處理與分析/conv01.m p1 = [1, 0, 1, 1]; p2 = [0, 1, -1, 2]; p3 = conv(p1, p2) p3 = 0 1 -1 3 0 1 2 上例中的 與 多項式的乘積結果是 ，若改以多項式的表示法，即為 p1 p2 p3 p 3(x) 5 4 2 = x −x + 3x −x + 2 。 若欲求p 1(x ) 除以p 2 (x ) 所得的商式與餘式，可輸入： Example 3: 07-多項式的處理與分析/deconv01.m p1 = [1, 0, 1, 1]; p2 = [1, -1, 2]; [q, r] = deconv(p1, p2) q = 1 1 r = 0 0 0 -1 此即代表p 1(x) 除以p 2 (x) 後，得到的商式為q(x)