科学计算曲线拟合.pptxVIP

数据拟合方法;第四章数据拟合方法;例：考察某种纤维的强度与其拉伸倍数的关系,下表是实际测定的24个纤维样品的强度与相应的拉伸倍数的记录:;24个点大致分布在一条直线附近。;例已知实验数据如下,求拟合曲线。;线性函数拟合:假设拟合函数为

y=a+bx

需确定a和b.其几何意义是,确定一条直线使已有的数据点尽可能接近直线。;为了选取a、b,使得函数f(a,b)取极小值.令;xx1x2……xm

f(x)y1y2……ym;据实验数据分布特点选取，可选幂函数类、指数函数类、三角函数类等。;残差向量的各分量平方和记为：;由多元函数求极值的必要条件，有;上式为由n+1个方程组成的方程组，称正规方程组。;引入记号;将其表示成矩阵形式：;作为一种简单的情况，常使用多项式函数Pn(x)作为(xi,yi)(i=1,2,…,m)的拟合函数。;即正规方程组为;例.回到本节开始的实例，从散点图可以看出，;解得;若mn+1，则此方程组称超定方程组（方程个数未知数个数）;考虑正规方程组;故正规方程组矩阵形式为：;x12 3 4

y4101826;记系数矩阵为?，则;注：具体用几次多项式拟合，可据实际情况而定。可先画草图，将已知点描上去，看与什么函数相近，就以什么函数拟合。;指数函数拟合的例子;t;经计算N(2000)=64.1805。所以2000年的世界人口预测为64.1805亿。;MATLAB中的多项式拟合命令;x=-5:5;y=1./(1+x.^2);

t=-5:.05:5;y0=1./(1+t.^2);

p=polyfit(x,y,10);y1=polyval(p,t);

plot(t,y0,x,y,o,t,y1,.)