MATLAB程式设计：进阶篇多项式的处理与分析.pptVIP

不同次數的比較 (I) 不同的次數進行多項式擬合，「外插」（Extrapolation）時得到的差異相當大 範例7-17：polyfit02.m load census.mat cdate2 = min(cdate):(max(cdate)+30); curve = zeros(7, length(cdate2)); for i = 1:7 curve(i,:) = polyval(polyfit(cdate,pop,i), cdate2); end plot(cdate, pop,o, cdate2, curve); legend(實際資料, order=1, order=2,order=3, order=4, order=5, order=6, order=7); xlabel(年份); ylabel(人口（單位：百萬）); 不同次數的比較 (II) 次數越高，越能逼近給定的資料 不表示預測的準確度會提高 高次多項式對雜訊（Noise）敏感度較高，容易產生不準確的外插預測 多項式擬合的誤差及錯誤 上述範例中，常出現 MATLAB 警告訊息 表示實際資料的不足 或資料數值範圍的過大或過小 在進行多項式擬合時，最好先將資料進行必要的平移、放大或縮小範圍 MATLAB 程式設計進階篇：多項式的處理與分析 MATLAB 程式設計：進階篇多項式的處理與分析 張智星 (Roger Jang) jang@ /jang 台大資工系 多媒體檢索實驗室 7-1多項式的加減乘除 n 次多項式的表示法 MATLAB 中，用 n+1 的列向量表示 p(x) Ex. p = [1,2,3,1] 表示一個三次多項式 多項式的加減 (I) 多項式的加減，直接由向量的加減推出 兩多項式， 及 計算其和與差 範例7-1：polyPlus.m p1 = [1,0,1,1]; p2 = [0,1,-1,2]; p1 + p2 p1 - p2 多項式的加減 (II) 矩陣 p1 與 p2 的長度要一致，否則 MATLAB 會產生運算錯誤的訊息 ans = 1 1 0 3 ans = 1 -1 2 -1 多項式的乘除 (I) 多項式的乘除，使用 conv 及 deconv 指令達成 欲求多項式 與 的乘積 範例7-2：conv01.m 以多項式的表示法，為 p1 = [1, 0, 1, 1]; p2 = [0, 1, -1, 2]; p3 = conv(p1, p2) p3 = 0 1 -1 3 0 1 2 多項式的乘除 (II) 欲求 除以 的商式與餘式 範例7-3：deconv01.m 代表 除以 ，得到商式 ，餘式 p1 = [1, 0, 1, 1]; p2 = [1, -1, 2]; [q, r] = deconv(p1, p2) q = 1 1 r = 0 0 0 -1 多項式的加減乘除 函數 說明 p1 + p2 多項式 與 的和 p1 - p2 多項式 與 的差 conv(p1, p2) 多項式 與 的乘積 [q, r] = deconv(p1, p2) 多項式 除以 ，得到商式為 ，餘式為 7-2 多項式的求值、求根、微分與積分 計算多項式的值，使用 polyval 指令 範例7-4：polyval01.m p = [1 2 1]; x = 0:0.1:3; y = polyval(p, x); plot(x, y, -o); 計算多項式的值-1 x 和 y 是長度為 31 的向量，y(i) 的值為　　　　　 在 x = x(i) 的函數值 計算多項式的值-2 計算p(A)，A為方陣，可用polyvalm指令 範例7-5：polyvalm01.m 結果和 B = A^2 + 2*A + 1 相同 p = [1 2 1]; A = [1 2; 3 4]; B = polyvalm(p, A) B = 10 14 21 31 計算多項式的根 (I) 多項式的根，可用 roots 指令求得 多項式 的根 範例7-6