matlab在数值分析中的应用.pptxVIP

matlab在数值分析中的应用会计学第1页/共41页7.1 多项式插值原理7.2 分段线性插值7.3 离散数据的多项式拟合7.3.1线性最小二乘法拟合7.3.2非线性最小二乘法拟合第2页/共41页 插值与拟合是来源于实际、又广泛应用于实际的两种重要方法。随着计算机的不断发展及计算水平的不断提高，它们已在国民生产和科学研究等方面扮演着越来越重要的角色。 P10(t) f(t) f(x)第3页/共41页7.1 多项式插值原理已知 f(x) 在点 xi 上的函数值 yi=f(xi), (i=0,1,2,···,n)求多项式 P(x)=a0 + a1x +···+ anxn满足: P(xk)= yk (k = 0,1,…,n)第4页/共41页7.2 分段线性插值 所谓分段线性插值就是通过插值点用折线段连接起来逼近原曲线，这也是计算机绘制图形的基本原理。实现分段线性插值不需编制函数程序，MATLAB自身提供了内部函数interp1其主要用法如下：interp1(x,y,xi) 一维插值第5页/共41页◆ yi=interp1(x,y,xi) 对一组点(x,y) 进行插值，计算插值点xi的函数值。x为节点向量值，y为对应的节点函数值。如果y 为矩阵，则插值对y 的每一列进行，若y 的维数超出x 或 xi 的维数，则返回NaN。◆ yi=interp1(y,xi) 此格式默认x=1:n ，n为向量y的元素个数值，或等于矩阵y的size(y,1)。第6页/共41页◆ yi=interp1(x,y,xi,’method’) method用来指定插值的算法。默认为线性算法。其值常用的可以是如下的字符串。 ● nearest线性最近项插值。 ● linear线性插值。 ● spline三次样条插值。 ● cubic三次插值。 ● pchip分段三次Hermite插值。 所有的插值方法要求x是单调的。x 也可能并非连续等距的。第7页/共41页正弦曲线的插值示例：x = 0:10; y = x.*sin(x); xx = 0:.25:10; yy = interp1(x,y,xx); plot(x,y,kd,xx,yy)第8页/共41页例 产品产量拟合（离散数据的一维插值 ） year = 1900:10:2010; product = [75.995 91.972 105.711 123.203 131.669 150.697 179.323 203.212 226.505 249.633 256.344 267.893 ];p1995 = interp1(year,product,1995)x = 1900:1:2010;y = interp1(year,product,x,pchip);plot(year,product,o,x,y)第9页/共41页插值结果为：p1995 = 252.9885第10页/共41页 Matlab也能够完成二维插值的运算，相应的函数为interp2，使用方法与interpl基本相同，只是输入和输出的参数为矩阵，对应于二维平面上的数据点。二维插值zi=interp2(x, y, z, xi, yi, ‘method’)第11页/共41页 返回矩阵ZI，其元素包含对应于参量XI与YI（可以是向量、或同型矩阵）的元素，即Zi(i,j)←[Xi(i,j),yi(i,j)]。用户可以输入行向量和列向量Xi与Yi，此时，输出向量Zi与矩阵meshgrid(xi,yi)是同型的。同时取决于由输入矩阵X、Y与Z确定的二维函数Z=f(X,Y)。参量X与Y必须是单调的，且相同的划分格式，就像由命令meshgrid生成的一样。若Xi与Yi中有在X与Y范围之外的点，则相应地返回NAN（Not a Number）。第12页/共41页例：[X,Y] = meshgrid(-3:.25:3);Z = peaks(X,Y);[XI,YI] = meshgrid(-3:.125:3);ZZ = interp2(X,Y,Z,XI,YI);surfl(X,Y,Z);hold on;surfl(XI,YI,ZZ+15)axis([-3 3 -3 3 -5 20]);shading flathold off二维插值图第13页/共41页例years = 1950:10:1990;service = 10:10:30;wage = [150.697 199.592 187.625 179.323 195.072 250.287 203.212 179.092 322.767 226.505 153.706 426.730 249.633 120.281 598.243];w = interp2(service,y