matlab在科学计算中的应用5+-+多项式、插值与数据拟合.pptVIP

5.3.1 多项式拟合 n次多项式： 曲线与数据点的残差为： 残差的平方和为： 为使其最小化，可令R关于 的偏导数为零，即： 或 或矩阵形式： 多项式拟合MATLAB命令：polyfit格式：p=polyfit(x,y,n) 选较密的插值点，用默认的线性插值算法进行插值 [x1,y1]=meshgrid(-3:.2:3,-2:.2:2); z1=interp2(x,y,z,x1,y1); surf(x1,y1,z1),axis([-3,3,-2,2,-0.7,1.5]) 立方和样条插值： z1=interp2(x,y,z,x1,y1,cubic); z2=interp2(x,y,z,x1,y1,spline); surf(x1,y1,z1),axis([-3,3,-2,2,-0.7,1.5]) figure;surf(x1,y1,z2),axis([-3,3,-2,2,-0.7,1.5]) 算法误差的比较 z=(x1.^2-2*x1).*exp(-x1.^2-y1.^2-x1.*y1); surf(x1,y1,abs(z-z1)),axis([-3,3,-2,2,0,0.08]) figure;surf(x1,y1,abs(z-z2)),axis([-3,3,-2,2,0,0.025]) 二维一般分布数据的插值 功能:可对