MATLAB應用(擬合).ppt

三數(shù)據(jù)擬合 用插值的方法對一函數(shù)進行近似 要求所得到的插值多項式經(jīng)過已知插值節(jié)點 在n比較大的情況下 插值多項式往往是高次多項式 這也就容易出現(xiàn)振蕩現(xiàn)象 龍格現(xiàn)象 即雖然在插值節(jié)點上沒有誤差 但在插值節(jié)點之外插值誤差變得很大 從 整體 上看 插值逼近效果將變得 很差 所謂數(shù)據(jù)擬合是求一個簡單的函數(shù) 例如是一個低次多項式 不要求通過已知的這些點 而是要求在整體上 盡量好 的逼近原函數(shù) 這時 在每個已知點上就會有誤差 數(shù)據(jù)擬合就是從整體上使誤差 盡量的小一些 3 1多項式擬合 n次多項式 曲線與數(shù)據(jù)點的殘差為 殘差的平方和為 為使其最小化 可令R關于的偏導數(shù)為零 即 或或矩陣形式 多項式擬合MATLAB命令 polyfit格式 p polyfit x y n x0 0 1 1 y0 x0 2 3 x0 5 exp 5 x0 sin x0 p3 polyfit x0 y0 3 vpa poly2sym p3 10 可以如下顯示多項式ans 2 839962923 x 3 4 789842696 x 2 1 943211631 x 5975248921e 1 例 繪制擬合曲線 x 0 01 1 ya x 2 3 x 5 exp 5 x sin x y1 polyval p3 x plot x y1 x ya x0 y0 o 就不同的次數(shù)進行擬合 p4 polyfit x0 y0 4 y2 polyval p4 x p5 polyfit x0 y0 5 y3 polyval p5 x p8 polyfit x0 y0 8 y4 polyval p8 x plot x ya x0 y0 o x y2 x y3 x y4 擬合最高次數(shù)為8的多項式 vpa poly2sym p8 5 ans 8 2586 x 8 43 566 x 7 101 98 x 6 140 22 x 5 125 29 x 4 74 450 x 3 27 672 x 2 4 9869 x 42037e 6Taylor冪級數(shù)展開 symsx y x 2 3 x 5 exp 5 x sin x vpa taylor y 9 5 ans 5 x 28 x 2 77 667 x 3 142 x 4 192 17 x 5 204 96 x 6 179 13 x 7 131 67 x 8多項式表示數(shù)據(jù)模型是不唯一的 即是兩個多項式函數(shù)完全不同 在某一區(qū)域內(nèi)其曲線將特別近似 多項式擬合的效果并不一定總是很精確的 x0 1 2 0 10 10 y0 1 1 25 x0 2 x 1 01 1 ya 1 1 25 x 2 p3 polyfit x0 y0 3 y1 polyval p3 x p5 polyfit x0 y0 5 y2 polyval p5 x p8 polyfit x0 y0 8 y3 polyval p8 x p10 polyfit x0 y0 10 y4 polyval p10 x plot x ya x y1 x y2 x y3 x y4 例 用Taylor冪級數(shù)展開效果將更差 symsx y 1 1 25 x 2 p taylor y x 10 p 1 25 x 2 625 x 4 15625 x 6 390625 x 8多項式擬合效果 x1 1 0 01 1 ya 1 1 25 x1 2 y1 subs p x x1 plot x1 ya x1 y1 3 2函數(shù)線性組合的曲線擬合方法 該方程的最小二乘解為 其中 例 x 0 0 2 0 4 0 7 0 9 0 92 0 99 1 2 1 4 1 48 1 5 y 2 88 2 2576 1 9683 1 9258 2 0862 2 109 2 1979 2 5409 2 9627 3 155 3 2052 A ones size x exp 3 x cos 2 x exp 4 x x 2 c A y c1 c c1 1 22002 3397 0 67970 8700 圖形顯示 x0 0 0 01 1 5 A1 ones size x0 exp 3 x0 cos 2 x0 exp 4 x0 x0 2 y1 A1 c plot x0 y1 x y x 數(shù)據(jù)分析 x 1 1052 1 2214 1 3499 1 4918 1 6487 1 8221 2 0138 2 2255 2 4596 2 7183 3 6693 y 0 6795 0 6006 0 5309 0 4693 0 4148 0 3666 0 3241 0 2864 0 2532 0 2238 0 1546 plot x y x y 例 分別對x y進行對數(shù)變換 x1 log x y1 log y plot x1 y1 A x1 ones size x1 c A y1 c 1 2339 0 2630 exp c 2 ans 0 7687 x 0 0 1 1 y x 2 3 x 5 exp 5 x sin x n 8 A fori 1 n 1 A i x n 1 i end c A y vpa poly2sym c 5 ans 8 2586 x 8 43 566 x 7 101 98 x 6 140 22 x 5 125 29 x 4 74 450 x 3 27 672 x 2 4 9869 x 42037e 6 例 3 3最小二乘曲線擬合 格式 a jm lsqcurvefit Fun a0 x y 例 x 0 1 10 y 0 12 exp 0 213 x 0 54 exp 0 17 x sin 1 23 x f inline a 1 exp a 2 x a 3 exp a 4 x sin a 5 x a x xx res lsqcurvefit f 1 1 1 1 1 x y xx resOptimizationterminatedsuccessfully RelativefunctionvaluechangingbylessthanOPTIONS TolFunans 0 11970 21250 54040 17021 2300res 7 1637e 007 修改最優(yōu)化選項 ff optimset ff TolFun 1e 20 ff TolX 1e 15 修改精度限制 xx res lsqcurvefit f 1 1 1 1 1 x y ff xx resOptimizationterminatedsuccessfully RelativefunctionvaluechangingbylessthanOPTIONS TolFunans 0 12000 21300 54000 17001 2300res 9 5035e 021 繪制曲線 x1 0 0 01 10 y1 f xx x1 plot x1 y1 x y o 例 x 0 1 0 1 1 y 2 3201 2 6470 2 9707 3 2885 3 6008 3 9090 4 2147 4 5191 4 8232 5 1275 functiony c8f3 a x y a 1 x a 2 x 2 exp a 3 x a 4 a lsqcurvefit c8f3 1 2 2 3 x y a Maximumnumberoffunctionevaluationsexceeded increaseoptions MaxFunEvalsans 2 45752 45571 44372 0720 繪制曲線 y1 c8f3 a x plot x y x y1 3 4B樣條函數(shù)及其MATLAB表示 格式S spapi k x y 例 x0 0 0 4 1 2 pi y0 sin x0 ezplot sin t 0 pi holdon sp1 csapi x0 y0 fnplt sp1 三次分段多項式樣條插