最小二乘法解决人口预测问题假设人口满足指数增长类型.x(t)=xo×e^rt,r为人口增长率,x.为

网友回答

最小二乘法与数据拟合
一、问题某公交公司1路车过去20个季度内的客流量(单位：百万)如下表：
季度序号 1 2 3 4 5 6 7 8 9 10
客流量 1.85 2.18 1.6 2.31 1.93 2.35 1.638 2.51 1.92 2.49
季度序号 11 12 13 14 15 16 17 18 19 20
客流量 2.01 2.68 2.05 2.77 2.23 2.95 2.46 2.87 2.35 3.07
试确定客流量与季度序号之间的函数关系,并预测未来八个季度的客流量.
二、实验目的
掌握最小二乘法的原理并会用于解决实际问题；学会用mathematica中的算符进行曲线拟合.
三、预备知识
1 最小二乘法
在许多实际问题中,往往需要根据实验测得两个变量x与y的若干组实验数据(x1,y1),…(xn,yn)来建立这两个变量的函数关系的近似式,这样得到的函数近似式称为经验公式.
通过对实验数据的处理,能够判断x、y大体上满足某种类型的函数关系y＝f(x,a1,a2,…,as),但是其中s个参数a1,a2,…,as的值需要通过n组实验数据来确定,通常可以这样来确定参数：选择参数a1,a2,…,as,使得f(x,a1,a2,…,as)在x1,x2 …xn处的函数值与实验数据 y1,y2 …yn 的偏差的平方和为最小,就是使
d＝ (1)为最小,这种方法称为最小二乘法.当f(x,a1,a2,…,as)是s个参数的线性函数时,利用求极值与解线性方程组的方法可以解决.
例如,若x、y大体上满足线性关系即f(x,a,b)＝ax+b ,则
d(a,b)＝ (2)
由多元极值的求法有
(3)解上述关于a、b的二元一次方程组得
a＝ , b＝ 从而求得经验公式y＝ax＋b .
d＝ 的大小是衡量经验公式精度的一种尺度.
线性函数是最简单最常用的经验公式,有一些实际问题,它们的经验公式可能不是线性函数,我们可以把它化为线性函数来讨论,例如y=kemx,两边取对数得lny=mx+lnk ,令z=lny,b=lnk,即可化为z＝mx＋b.
2 mathematica中数据拟合算符的用法
在数据处理中常常设法用一个函数按照某种法则去描述一组数据,这就是数据拟合.上面介绍的最小二乘法就是一种最常用的数据拟合方法.mathematica中最基本的数据拟合算符是 fit[ ] ,语法为
fit[数据,拟合函数的基函数列表,变量]
线性函数拟合的基函数为1,x ,n阶多项式拟合的基函数是1,x,x2,…xn .
例 一册书的成本费y与印刷的册数x有关,统计数据如下：
xi(千册) 1 2 3 4 5 6 7 8 9 10
yi（元） 10.15 5.52 4.08 2.85 2.11 1.62 1.41 1.30 1.21 1.15
试用y=a+ 去拟合上述数据.
mathematica程序及运行结果如下：
data={10.15,5.52,4.08,2.85,2.11,1.62,1.41,1.30,1.21,1.15};
fit[data,{1,1/x},x]
四、实验内容与要求
1 画出实验问题的数据图,并粗略估计这些数据与什么类型的函数比较吻合?
2 取经验公式为线性函数y＝ax＋b 按照最小二乘法的原理用mathematica编程解实验问题.
3 取经验公式为y＝ax＋b ＋c sin[ x]+d cos[ x] ,用mathematica中算符fit[]来求解实验问题,并与内容2的精度比较,对比实际情况,你能得出什么?
五、操作提示
12 拟合程序及运行结果如下： 预测程序及运行结果如下：3 程序及运行结果如下：计算两种经验公式的精度可以看出第二种较好,这与客流量呈季节被动变动的实际情况吻合.怎样用 mathematica 拟合二元函数?数据拟合 由一组已知数据（xk,yk）（k=1,2,…,n）,求函数的近似解析式y=f（x）,就是数据拟合问题,当然函数还可以是多元的. Mathematica提供了进行数据拟合的函数： Fit[data,funs,vars] 对数据data用最小二乘法求函数表funs中各函数的一个线性组合作为所求的近似解析式,其中vars是自变量或自变量的表. 例如： Fit[data,{1,x},x] 求形为y=a+bx的近似函数式. Fit[data,{1,x,x2},x] 求形为y=a+b