参数的点估计

  • 格式:docx
  • 大小:13.45 KB
  • 文档页数:3

实验一 参数的点估计

实验目的:

通过本实验,使学生以 Matlab 为工具掌握参数点估计的计算方法的计算机实现;对常见分布,掌握生成点估计量值的模拟方法,通过观察不同样本量下估计量的值在真实参数周围的分布情况,获得估计量的值在真实参数周围分布情况及其随样本量增加所发生变化的数值经验.

实验过程:

本实验选择使用服从正态分布的数据进行设计,通过函数生成一组服从正太分布的样本,进行点估计。

通过matlab的normrnd函数生成均值为15,方差为2.52,容量为10的一组随机正太分布数据。

>> w = normrnd(15,2.5,1,10)

w =

12.2059 18.1516 16.6504 14.8303 14.5119 14.4560 14.2422

15.0576 15.1282 17.0652

总体w均值的矩估计:

>> mu_ju = mean(w)

mu_ju =

15.2299

总体w方差的矩估计:

>> sigma2_ju = moment(w,2)

sigma2_ju =

2.5357

点估计的最大似然法,使用matlab的mle函数,进行置信度为97%的参数估计,得到参数的均值的点估计和方差的点估计值。

paramhat中有两个数值,第一个是正态分布均值(mu)的点估计,后一个是正态分布的标准差(sigma)。paramint中第1,2两列分布对应均值和标准差的区间估计。从结果来看,mu和sigma的区间估计都很窄,而且很接近我们生成随机数时使用的参数值。

>> alpha=0.03;

>> [paramhat,paramint] = mle('norm',w,alpha)

paramhat =

15.2299 1.5924

paramint =

13.8638 1.1118

16.5961 3.2955