第六章计算方法.pptVIP

* 第六章 计算方法 非线性方程求解 多项式插值与曲线拟合 数值微分与数值积分 求常微分方程数值解命令 ? ? ? ? x y (k = 1,2,·······) 已知函数 f(x),求 x 使 f(x)=0 —— 非线性方程 例 x = tan x f=inline(x-tan(x)) x=fzero(f,3.5) x = 4.4934 求一元函数零点 fzero(fun, x0) X = fzero(fun, x0)寻找fun在x0附近的零点 例 x 3 – 30 x 2 + 2552 =0 p=[1 -30 0 2552]; roots(p) ans = 26.3146 11.8615 -8.1761 r x 半径r =10 cm的球体,密度? =0.638.浸入水深度 x=? 15/18 求多项式零点: roots(p); 例1.求方程 x3 – 4x + 5 = 0 的解 P=[1 0 -4 5]; roots(P) ans = -2.4567 1.2283 + 0.7256i 1.2283 - 0.7256i 例2 求 x3 –8x2 +6x– 30=0的解 P=[1 -8 6 -30]; r=roots(P) r = 7.7260 0.1370 + 1.9658i 0.1370 - 1.9658i 多项式插值原理 已知 f(x) 在点 xi 上的函数值 yi=f(xi), (i=0,1,2,···,n) 求多项式 P(x)=a0 + a1x +···+ anxn 满足: P(xk)= yk (k = 0,1,…,n) P10(t) f(t) f(x) 一维插值: yi = interp1(x, y, xi, ‘method ’) 立方插值 cubic 样条插值 spline 线性插值 linear 最近点插值 nearest method x = 0:10; y = sin(x); xi = 0:.25:10; yi = interp1(x,y,xi); plot(x,y,o,xi,yi) 二维插值 zi=interp2(x, y, z, xi, yi, ‘method’) 三维插值 vi = interp3(x,y,z,v,xi,yi,zi, ‘method’) 17/18 例1.正弦曲线数据插值试验 已知数据表 x x1 x2 ·········· xm f(x) y1 y2 ·········· ym 求拟合函数: ? (x) = a0 + a1x + ······+ anxn 使得 离散数据的多项式拟合 达到最小 多项式拟合命令 P=polyfit(X,Y,n) 求出(最小二乘意义下)n次拟合多项式 P(x)=a0xn + a1xn-1 + ······ + an-1x + an 计算结果为系数: P=[ a0, a1, ······ , an-1, an ] 多项式求值命令: y1=polyval(P, x) 其中,P是n次多项式的系数,x是自变量的值,y1是多项式在x处的值 17/18 4/16 定积分与积分和式 右矩形和 h 1 0.5 0.2 ······ Sn 5.2908 5.1044 4.9835 ······ 4.8999 5/16 左矩形 梯形 右矩形 4.4429 4.8669 5.2908 4.6804 4.8924 5.1044 4.8139 4.8987 4.9835 4.8572 4.8996 4.9420 梯形数