数据拟合的非线性模型多项式拟合的正交化方法超定方程组最小二乘解.pptVIP

* 数据拟合的非线性模型 多项式拟合的正交化方法 超定方程组最小二乘解 温度数据的平滑处理 《数值分析》 18 ? ? ? ? 数据拟合的非线性模型 求拟合函数 f(x, c0, c1, ···, cn )满足 观测数据 2/16 例1.已知人口统计数据 利用最小二乘法求指数拟合 y = c e a x x x1 x2 ·········· xm f y1 y2 ·········· ym 年 1991 1992 1993 1994 1995 1996 数量 11.58 11.72 11.85 11.98 12.11 12.24 指数函数拟合人口统计数据(单位：亿) (1)计算对数值 yk = ln Nk ( k = 1,2,…,6) (2) 列出未知数a 、b的超定方程组 a + b t k = y k ( k = 1，2，……，6) (3)求超定方程组的最小二乘解,得 N= ea+bx , 并预测 N(2000), N(2008) N2000=12.7971, N2008=13.9783 设 N= ea+bt ? lnN=a+bt 令 y = lnN, 有 y(t)=a+bt t 1991 1992 1993 1994 1995 1996 ln N 2.45 2.46 2.47 2.48 2.49 2.50 3/16 例2.利用极坐标观察值确定慧星轨道 其中，p 为参数，e为偏心率 令k=1/r r 2.70 2.00 1.61 1.20 1.02 ? 48? 67? 83? 108? 126? k 0.3704 0.5000 0.6211 0.8333 0.9804 ? 0.8378 1.1694 1.4486 1.8850 2.1991 线性方程组 ( j= 1,2,3,4,5 ) ? ? p=?, e = ? 4/16 r=[2.70 2.00 1.61 1.20 1.02]; fai=[48 67 83 108 126]*pi/180; n=length(r);b=ones(n,1); G=[1./r cos(fai)]; X=G\b; p=X(1);e=X(2); p = 1.4521, e = 0.7027 残差平方和: 0.0022 5/16 fa=pi/4:.1:3.5*pi/2; ro=p./(1-e*cos(fa)); polar(fai,r,ro),hold on polar(fa,ro) R=r-p./(1-e*cos(fai)) R=sum(R.^2) 例3. 温度数据的拟合问题: 中午1:00到第二天中午12:00温度记录数据如下 66, 66, 65, 64, 63, 63, 62, 61 60, 60, 59, 58, 58, 58, 58, 58 57, 57, 57, 58, 60, 64, 67, 68 用最小二乘法求拟合函数 f(x)=E + A cos(C x) + B sin( D x) E=?, A=？, B=? , C=? , D=? 6/16 定义 设 ?j(x) (j=0,1,···,n )是多项式系列, 是一点列, 多项式?i (x), ?j (x)在该 点列上的值分别为 则称 ?k(x) (k=0,1,···,n )为点列 上的 正交多项式系. (j ≠i) 如果 ?i = [?i(x1), ?i(x2), ······, ?i(xm)]T ?j = [?j(x1), ?j(x2), ······, ?j(xm)]T 7/16 ····························· p0(x) = 1, pj(x) = x pj-1(x) ( j=1,2,····,n ) ( j=2,3,···,n ) ? 8/16 给定数据表 x x1 x2 ··········