数值分析13线性插值和二次插值公式课件.ppt

插值方法的应用 代数插值问题 线性插值与二次插值公式 拉格朗日插值公式;引例1.正弦函数 sin x 的计算问题;(1)复杂函数的计算; (2)函数表中非表格点计算 (3)光滑曲线的绘制; (4)提高照片分辩率算法 (5)定积分的离散化处理; (6)微分方程的离散化处理; (7)积分方程的离散化处理;;引例2. 误差函数;已知f(x)在点xi上的函数值 yi=f(xi), (i=0,1,2,···,n);点,则满足插值条件 P(xk)= yk (k = 0,1,…,n) 的n次插值多项式 P(x)=a0 + a1x +……+ anxn 存在而且是唯一的。;方程组系数矩阵取行列式;构造3次多项式P(x) 逼近 Erf(x);x=0:.6:1.8; y=erf(x); x=x;y=y; A=[ones(4,1) x x.^2 x.^3]; p=A\y; a0=p(1);a1=p(2); a2=p(3);a3=p(4); t=0:.2:2; u=a0+a1*t+a2*t.^2+a3*t.^3; plot(x,y,o,t,u);过两点直线方程;记;二次插值问题;二次插值函数: L(x)=l0(x)y0+l1(x)y1+l2(x)y2，;二次插值基函数图形;二次插值的一个应用——极值点近似计??;拉格朗日插值公式;Runge反例: , (-5≤x≤5);x=-5:5;y=1./(1+x.^2); t=-5:0.05:5;y1=1./(1+t.^2);n=length(t); for i=1:n z=t(i);s=0; for k=1:11 Lk=1;u=x(k); for j=1:11 if j~=k,Lk=Lk*(z-x(j))/(u-x(j)); end end s=s+Lk*y(k); end y2(i)=s; end plot(x,y,ko,t,y1,t,y2,r)