第四次课拟合与插值.pptxVIP

拟合与插值 问题的提出 在一些实际问题中，有时只知道未知函数y=f(x)在有限个点x0,x1,…,xn上的值f(x0),f(x1),…,f(xn)，这相当于已知一个数表 然后根据表中的数据来构造一个简单的函数p(x)作为未知函数f(x)的近似函数，去参与有关f(x)的运算，这类问题称为数据逼近问题. 拟合 插值 目录 CONTENTS 曲线拟合 已知平面上若干个点 ， 求一个较简单的函数（曲线） , 使 尽可能的靠近数据点，在某种意义下达到最优。 称为拟合函数。 曲线拟合的一般提法 使残差的平方和最小，即 最小二乘法 多项式拟合 Matlab多项式拟合命令---拟合函数为多项式 将多项式表示成习惯的形式。 是多项式系数，字符 为自变量 按数组规则计算 处多项式的值 功能: 功能: 功能: X y=a1xn+a2xn-1+……...+anx+an+1 p=[a1,a2,...,an,an+1] x=1:0.1:2; y=[2.1,3.2,2.1,2.5,3.2,3.5,3.4,4.1,4.7,5.0,4.8]; p2=polyfit(x,y,2) p3=polyfit(x,y,3); p7=polyfit(x,y,7); disp(‘二次拟合函数),f2=poly2str(p2,x) disp(‘三次拟合函数),f3=poly2str(p3,x); disp(‘七次拟合函数),f7=poly2str(p7,x); x1=1:0.01:2; y2=polyval(p2,x1); y3=polyval(p3,x1); y7=polyval(p7,x1); plot(x,y,rp,x1,y2,--,x1,y3,k-.,x1,y7); 例 多项式拟合 p2= 1.3869 -1.2608 2.141 二次拟合函数 f2 = 1.3869 x^2 - 1.2608 x + 2.141 拟合应用——人口预测 下表是1971年到1990年我国总人口的统计数字，试根据1971年到1985年这15年人口的统计数字用多种方法预测未来20年的人口数字，并比较1986年到1990年间预测人口数字与实际统计数字的差异，在你所使用的几种预测方法中找出一种较为合理的预测方法。 年份 人口统计数字 年份 人口统计数字 1971 8.5229 1981 10.0072 1972 8.7177 1982 10.1654 1973 8.9211 1983 10.3008 1974 9.0859 1984 10.4357 1975 9.2420 1985 10.5851 1976 9.3717 1986 10.7507 1977 9.4974 1987 10.9300 1978 9.6259 1988 11.1026 1979 9.7542 1989 11.2704 1980 9.8705 1990 11.4333 拟合应用——人口预测 年份 统计数字 3次拟合 4次拟合 　 年份 统计数字 3次拟合 4次拟合 1971 8.5229 8.5263 8.5143 　 1989 11.2704 11.4329 11.1129 1972 8.7177 8.7278 8.7329 　 1990 11.4333 11.7003 11.2058 1973 8.9211 8.9105 8.9209 　 1991 　 11.9965 11.2704 1974 9.0859 9.0771 9.0855 　 1992 　 12.3242 11.2991 1975 9.242 9.2301 9.2331 　 1993 　 12.6859 11.283 1976 9.3717 9.3722 9.3692 　 1994 　 13.0844 11.2128 1977 9.4974 9.5061 9.4987 　 1995 　 13.5222 11.078 1978 9.6259 9.6344 9.6254 　 1996 　 14.0021 10.8676 1979 9.7542 9.7598 9.7524 　 1997 　 14.5266 10.5696 1980 9.8705 9.8848 9.88 　 1998 　 15.0984 10.171 1981 10.0072 10.0122 10.0152 　 1999 　 15.7201 9.6583 1982 10.1654 10.1446 10.153 　 2000 　 16.3944 9.017 1983 10.3008 1