.调查了北京市在1988—2006年期间每百户家庭平均拥有的照相机数,具体数据如下:年份 y 年份 y 1988 1 7.5 1998 11 59.6 1989 2 9.8 1999 12 62.2 1990 3 11.4 2000 13 66.5 1991 4 13.3 2001 14 72.7 1992 5 17.2 2002 15 77.2 1993 6 20.6 2003 16 82.4 1994 7 29.1 2004 17 85.4 1995 8 34.6 2005 18 86.8 1996 9 47.4 2006 19 87.2 1997 10 55.5 试针对以下两种情况拟合如下的回归函数:y=1/(1/u+a*b^t) (1)已知u=100,采用线性化方法进行拟合,即求参数a,b的估计; (2)u未知,采用非线性最小二乘法,求参数u,a,b的估计.
时间: 2024-04-01 07:33:35 浏览: 55
(1) 对回归方程进行线性化处理,即将原方程变形为 ln(1/y-1/u) = ln(a) + t*ln(b),令y* = 1/y-1/u,得到新的回归方程:ln(y*) = ln(a) + t*ln(b)
将数据带入新方程,得到如下表格:
| 年份 t | ln(y*) |
|--------|--------|
| 1 | 1.013 |
| 2 | 0.988 |
| 3 | 0.962 |
| 4 | 0.933 |
| 5 | 0.892 |
| 6 | 0.839 |
| 7 | 0.624 |
| 8 | 0.501 |
| 9 | 0.152 |
| 10 | -0.170 |
| 11 | -0.428 |
| 12 | -0.479 |
| 13 | -0.804 |
| 14 | -1.135 |
| 15 | -1.347 |
| 16 | -1.602 |
| 17 | -1.719 |
| 18 | -1.776 |
| 19 | -1.791 |
对上述数据进行最小二乘拟合,得到参数估计值:ln(a) = 0.617,ln(b) = 0.203,因此,a的估计值为exp(0.617) ≈ 1.854,b的估计值为exp(0.203) ≈ 1.225。
因此,拟合的回归方程为:y = 1/(1/100+1.854*1.225^t)。
(2) 对回归方程进行非线性最小二乘拟合,即将原方程变形为 y = 1/(1/u+a*b^t),将u、a、b看作参数,对样本数据进行拟合,使得残差平方和最小。
采用MATLAB软件中的nlinfit函数进行拟合,得到参数估计值:u的估计值为101.6,a的估计值为1.986,b的估计值为1.284。
因此,拟合的回归方程为:y = 1/(1/101.6+1.986*1.284^t)。
阅读全文