第二章-参数估计
数理统计: 参数估计方法
引例
设总体 X 服从参数为 的指数分布, 未知,
X1 , X 2 , , X n 是来自X的样本, x1 , x2 , , xn 是
相应的样本值,求 的矩估计量和矩估计值.
解 因为 E( X ) 所以 用样本矩替换总体矩, 得 的矩估计量
ˆ
1 n
n i 1
Xi
X
(
x)
1
e
x
,
x0
0,
其他.
但参数 未知。已知参数的取值范围,记为 。
给出样本的一组观察值,如何推断总体的分布?
【思路】给出 的估计,则得到对总体分布的推断。
【方法】根据一定的原则,从 中找到一个值(点) 作为的 估计。
点估计
2
点估计定义
设总体 X 的分布函数 F ( x; ) 的形式为已知,
的估计量.
4
二、估计量的评选标准 1. 无偏性
定义 若 X1, X 2 ,, X n 为总体 X 的一个样本,
是包含在总体 X 的分布中的待估参数, 若估计量ˆ ˆ( X1 , X 2 ,, X n )的数学期望 E(ˆ) 存在, 且对于任意 有
E(ˆ) 则称ˆ 是 的无偏估计量,否则称为有偏的.
(2) lim S 2 2 a.s. (强大数定律) n
即样本方差是总体方差2的强相合估计, 也是相合估计.
12
C. 样本标准差
其观察值:
S
S2
1 n1
n i 1
Xi
X
2
;
s
1 n1
n i 1
( xi
第二章参数估计(作业)
3 . 7 0 3 . 3 0
3 . 2 8 3 . 0 5
3 . 3 5 3 . 3 3
3 . 2 0 3 . 2 7
3 . 1 2 3 . 2 8
3 . 2 5 3 . 2 5
2 。构造两个总体方差比 1
2 的 95%的置信区间。 2
2 答案: 已知, x1 =3.33, =0.006, 根据自由度 n1=21-1=20 和 n2=21-1=20, x 2 =3.27, s12 =0.06, s2
z 2
s =3.31± 0.53,则该校大学生平均上网时间 n
的置信区间为(2.78,3.84) 。 当置信水平为 99%时,z/2=2.58 , x 的置信区间为(2.62,0.69,则该校大学生平均上网时间 n
3、在一项家电市场调查中,随机抽取了 200 个居民户,调查他们是否拥有某一品牌的电视 机。其中拥有该品牌电视机的家庭占 23%。求总体比例的置信区间,置信水平分别为 90% 和 95%。 答案:已知 n=200,P=23%,则
第二章参数估计
1、某快餐店想要估计每位顾客午餐的平均花费金额,在为期 3 周的时间里选取 49 名顾客
组成了一个简单随机样本。 (1) 假定总体标准差为 15 元,求样本均值的抽样标准误差; (2) 在 95%的置信水平下,求边际误差; (3) 如果样本均值为 120 元,求总体均值 的 95%的置信区间。
6、生产工序的方差是工序质量的一个重要度量。当方差较大时,需要对工序进行改进以减 小方差。两部机器生产的袋茶重量(单位:g)的数据如下:
机 3 3 器 . . 1 4 2 5 0 机 3 3 器 . . 2 2 3 2 8
3 . 2 2 3 . 3 0
第二章 参数估计.pdf
22、设总体 X 在区间 [, +1] 上服从均匀分布,则 的矩估计 ˆ =
;
3
D(ˆ) =
。
23、设总体 X ~ N(, 2 ) ,若 和 2 均未知, n 为样本容量,总体均值 的置 信水平为1 − 的置信区间为 (X − , X + ) ,则 的值为________;
24、在实际问题中求某参数的置信区间时,总是希望置信水平愈 愈好,而置
解: E(ˆ1) = E(ˆ2), D(ˆ1) D(ˆ2) . 12、设ˆ1 和ˆ2 均是未知参数 的无偏估计量,且 E(ˆ12 ) E(ˆ22 ) ,则其中的统计
量 更有效。
13、在参数的区间估计 (1,2 ) 中,当样本容量 n 固定时,精度2 −1 提高时,置
信度1 −
。
14、设 X1, X 2 ,, X n 是来自总体 X ~ N(,1) 的样本,则 的置信度为 0.95 的置
9、什么是最优无偏估计量? 10、什么是一致最小方差无偏估计量? 11、有效估计量和最优无偏估计量的关系是什么? 12、什么叫均方误差最小估计量? 13、叙述一致估计量的概念。 14、试述评价一个置信区间好坏的标准。 15、描述区间估计中样本容量、精度、置信度的关系。
三、单选题 1、设总体未知参数 的估计量 满足 E( ) = ,则 一定是 的( )
的关系为
。
6 、 称 统 计 量 T = T ( X1, X 2 ,, X n ) 为 可 估 函 数 g() 的 ( 弱 ) 一 致 估 计 量 是
指
。
7、判断对错:设总体 X ~ N(, 2 ) ,且 与 2 都未知,设 X1, X 2 ,..., X n 是来自
1
该总体的一个样本,设用矩法求得 的估计量为 ˆ1 、用极大似然法求得 的
第二章参数估计
第二章 参数估计【学习目标】1、掌握矩估计的替代原则;会求已知分布中未知参数的矩估计(值)2、熟练掌握极大似然估计的思想及求法3、估计量的评价标准:无偏性、有效性、相合性的定义4、统计量的无偏性的判断;两个无偏估计的有效性判断;会用Fisher 信息量及c-R 下界进行统计量的UMVUE 充分性判断5、掌握区间估计的定义6、单个正态总体均值的区间估计(包括方差已知、方差未知);单个正态总体方差的区间估计(包括均值已知、均值未知)7、两个正态总体均值差的区间估计(方差未知);两个正态总体方差比的区间估计 8、单侧置信区间的求法 【典型例题讲解】例1、设1,,n X X 是来自均匀分布(,1)U θθ+的总体的容量为n 的样本,其中θ-∞<<+∞为未知参数,试证:θ的极大似然估计量不止一个,例如1(1)ˆXθ=,2()ˆ1n X θ=-,3(1)()11ˆ()22n XXθ=+-都是θ的极大似然估计。
解:(,1)U θθ+分布的密度函数为11()0x f x θθ≤≤+⎧=⎨⎩其他似然函数(1)()11()0n x x L θθθ≤≤≤+⎧=⎨⎩其他由于在(1)()1n x x θθ≤≤≤+上()L θ为常数,所以凡是满足:(1)()ˆˆ1n x x θθ≤≤≤+的ˆθ均为θ的极大似然估计。
从而(1)1(1)ˆX θ=满足此条件,故1(1)ˆX θ=是θ的极大似然估计;(2)由于()(1)1n X X -≤,故2()(1)()2ˆˆ11n n X X X θθ=-≤≤=+,所以2()ˆ1n Xθ=-为θ的极大似然估计;(3)由于()(1)1n X X -≤,故(1)()(1)12n X X X +-≤,(1)()()12n n X X X ++≥,从而有3(1)()(1)()(1)()31111ˆˆ()()12222n n n XXXXXXθθ=+-≤≤≤++=+,故3ˆθ也为θ的极大似然估计。
第二章-参数估计性质的证明
一、参数估计1ˆβ与2ˆβ是Y 的线性函数仅说明2ˆβ是Y 的线性函数,1ˆβ的情况类似可以得到。
由2ˆβ的估计式, 即2ˆβ由随机变量Y 线性表出,从这一关系式还可理解到2ˆβ的随机性是由Y 带来的。
参数估计线性性质的重要性,是可以基于Y 的统计分布建立参数估计1ˆβ和2ˆβ统计分布,这对利用1ˆβ和2ˆβ对真实参数1β和2β的统计推断带来了极大的方便。
对于i k 有如下一些性质,1、0i k =∑。
2、221()i i k X X =-∑∑。
3、1i i k X =∑。
4、()1i i k X X -=∑。
二、最小二乘估计1ˆβ与2ˆβ的无偏性质 仅说明2ˆβ是2β的无偏估计,1ˆβ的无偏性类似可证。
由2ˆβ的关于Y 的线性表出式, 对2ˆβ求数学期望并考虑零均值假定(()0i E u =) 所以这就证明了2ˆβ具有无偏性。
同理有11ˆ()E ββ=。
三、最小二乘估计1ˆβ与2ˆβ的最小方差性质 对于OLS 估计式1ˆβ和2ˆβ,已知其方差为 这里只证明2ˆ()Var β最小,1ˆ()Var β最小的证明可以类似得出。
任设2β的另一个线性无偏估计为*2β,即 其中因为*2β也是2β的无偏估计,即*22()E ββ=,必须有0i w =∑,1i i w X =∑同时22i w σ=∑ [因为2()i Var Y σ=]上式最后一项中0= (因为0i w =∑,1i i w X =∑) 所以而20σ≥,因为i i w k ≠,则有2()0i i w k -≥,为此只有i i w k =时,*22ˆ()()Var Var ββ=,由于*2β是任意设定的2β的线性无偏估计式,这表明2β的OLS 估计式具有最小方差性。
一般地,将具有最小方差性的无偏估计量,称为该估计量满足有效性。
四、参数估计1ˆβ与2ˆβ的一致性 若样本容量n 趋于无穷时,有21()ni i X X =-→∞∑,则下面只证明2ˆβ具有一致性,1ˆβ的一致性类似可证。
数值分析答案第二章参数估计习题
f(x)= () { > − ex λ ) λ 0λ ( x解: λe , x ≥ 0
第二章 参数估计 1.设母体X具有负指数分布,它的分布密度 −λ x 为 λe , x ≥ 0 f(x)= 0, x < 0 其中 λ > 0 。试用矩法求的估计量。 解:x e(λ ) f(x)=
0
1
θ −1
dx =
θ θ +1
X 估计EX
X ∴θ = 1− X
1 e 5.设母体X的密度为 f ( x) = 2σ
−
x
σ
, −∞ < x < ∞
试求 σ 的最大似然估计;并问所得估计量是 否的无偏估计. ∑x x n 解: n 1 −σ 1 n − σ
i
L = ∏ f ( xi ) = ∏
i =1 i =1
ln L = n ln θ + (θ − 1)∑ ln xi
i
0, 其他 n
i =1
( θ >0 )
n i =1
d ln L n ^= − n = + ∑ ln xi = 0,∴θ θ i dθ ∑ ln xi
i
2矩法估计
EX =
用
X 用估计EX
+∞
−∞
∫ x ⋅ f ( x)dx = ∫ x ⋅θ ⋅ x
2
给定置信概率1−α 即
P ( x − uα
2
σ/ n
,有 uα ,使
2
P{ u ≤ uα } = 1 − α
第二章 参数估计
0
x 2de
x
2xe
x
dx
2
xde
x
0
x
0
0
2 e dx 2 2
0
9
例4:设X1, … , Xn为取自 N ( , 2 ) 总体的
样本,求参数 , 2 的矩估计。
: E( X ) D( X ) 2 E( X 2 ) [E( X )]2
极大似然法是由德国数学家G.F.Gauss在1821年提 出的.然而这个方法通常归于英国统计学家 R.A.Fisher,因为他在1912年里发现了这一方法,并 且首先研究了这种方法的性质.
设总体的密度函数为f(x,θ), θ为待估参数,θ∈Θ,Θ
为参数空间.当给定样本观察值 x (x1, x2 , xn )后,f(x,
以随便给的,所以根据统计思想建立各种点估计方法
和评价点估计的好坏标准便是估计问题的研究中心.
这里先介绍三个常用的标准:无偏性、有效性和一致
性.
1
有效性
^
^
设 i i ( X1,, X n ), i 1, 2分别是参数 的两个无偏估计,
^
^
^
^
若D 1 D 2 至少有一个n使 成立 , 则称 1比 2 有效.
总体k阶矩 样本k阶矩
k E(Xk )
Ak
1 n
n i 1
X
k i
的矩估计量是
约定:若
是未知参数的矩估计,则u()的矩
估计为u(
),
6
例2、:设X1, … , Xn为取自参数为的指数分布 总体的样本,求的矩估计。
参数估计2
n
e n
i
x !
i 1 n i 1
ii ) ln L( x1 , x 2 ,..., x n ; ) xi ln n ln xi !
i 1
xi ln L( x1 , x2 ,...,xn ; ) i 1 n 0 iii)令 : 1 n iv)解之得 : xi x为 的极大似然估计值 , n i 1 1 n X i X 为 的极大似然估计量 . n i 1
(1)正态分布N (u, 2 ) (2)指数分布Z ( ) (3)均匀分布U (a, b) (4)二项分布B(n, p) (3)泊松分布 ( ) 试求其中未知参数的矩 估计. 解 : (1)
因为X ~ N ( , 2 ), E ( X ) , D( X ) 2 故有 X ,
注2
若 为 的矩估计量, g ( )为 的连续函数, 亦称g ( )为g ( )
2 2 例如S n 为总体方差D( X )的矩估计量, 则S n S n 为标准差 D( X )
的矩估计量. 的矩估计量.
例1.1
设X 1 , X 2 ,..., X n为来自正态总体 X 的样本, X的分布为
i 1 n n
( X为连续型)
(1.4) (1.5)
或
L( x1 , x2 ,..., xn ) PX i xi ;
i 1
( X为离散型)
达到最大值
L( x1 , x2 ,..., xn ; ) max L( x1 , x2 ,..., xn ; )
(1) 利用求导法求极大然估 计步骤 i )建立似然函数: L( x1 , x 2 ,..., x n ; 1 , 2 ,..., r ) f ( xi ; 1 , 2 ,..., r )
第二章 参数估计2-3 区间估计
I=0.814
上页 下页 返回
钢厂铁水含碳量X 例3. 钢厂铁水含碳量 ~ N(µ,0.1082), 现在随机测定 该厂9炉铁水得 炉铁水得X=4.484,求在置信度为 求在置信度为0.95 的条件 该厂 炉铁水得 求在置信度为 下铁水平均含碳量的置信区间。 下铁水平均含碳量的置信区间。 解
置信区间为
上页
下页
返回
联合方差
上页
下页
返回
1、 µ1 - µ2的1-α置信区间 、 α (1)、 σ12 、σ22已知 、
由于 X −Y ~ N(µ1 − µ2 ,
选取
2 2 σ1 σ2
n1
+
n2
)
因此置信度为1-α 因此置信度为 α的µ1 - µ2置信区间可为
上页
下页
返回
(2)、σ12 、σ22未知,且n1,n2较大 如大于 、 未知, 较大(如大于 如大于50)
=27.5, ,
=6.26, ,
上页
下页
返回
测量一批铅锭的比重,设铅锭的比重X 例6. 测量一批铅锭的比重,设铅锭的比重 ~ N(µ, 现进行16次检测得铅锭的比重有 σ2),现进行 次检测得铅锭的比重有 现进行 次检测得铅锭的比重有X=2.705, , S2=0.0292,试求总体 的均值µ和方差 σ2置信度为 求总体X的均值 0.95 的置信区间。 的置信区间。 解 (1)求µ的置信区间 σ2未知 n=16,α=0.05. 求 的置信区间, 未知, α 选取 查表得 置信区间为
(二)、总体X数学期望 (二)、总体X数学期望µ未知 数学期望µ 样本X 的无偏估计. 样本 1,X2, • • • , Xn, 且S2是σ2的无偏估计
选取样本函数
《应用数理统计》吴翊李永乐第二章 参数估计课后习题参考答案
第二章 参数估计课后习题参考答案2.1 设总体X 服从二项分布()n X X X p p N B ,,,,11,,21 <<为其子样,求N 及p 的矩法估计。
解:()()()p Np X D Np X E -==1,令()⎪⎩⎪⎨⎧-==p Np S Np X 12解上述关于N 、p 的方程得:2.2 对容量为n 的子样,对密度函数22(),0(;)0,0x x f x x x ααααα⎧-⎪=⎨⎪≤≥⎩其中参数α的矩法估计。
解:122()()a E x xx dx ααα==-⎰22022()x x dx ααα=-⎰2321221333ααααααα=-=-= 所以 133a x α∧== 其中121,21(),,,n n x x x x x x x n=+++为n 个样本的观察值。
2.3 使用一测量仪器对同一值进行了12次独立测量,其结果为(单位:mm) 232.50,232.48,232.15,232.52,232.53,232.30 232.48,232.05,232.45,232.60,232.47,232.30 试用矩法估计测量的真值和方差(设仪器无系统差)。
⎪⎪⎩⎪⎪⎨⎧-=-==X S p S X X p X N 2221ˆˆˆ解:()()()∑∑====-====ni i ni i S X X n X D X X n X E 12210255.014025.23212.4 设子样1.3,0.6,1.7,2.2,0.3,1.1是来自具有密度函数()10,1,<<=βββx f 的总体,试用矩法估计总体均值、总体方差及参数β。
解:()()()()4.22ˆ2,1,407.012.1101221========-===⎰⎰∑∑==X Xdx xdx x xf X E x f XX n S X n X ni i ni i ββββββββ参数:总体方差:总体均值:2.5 设n X X X ,,,21 为()1N ,μ的一个字样,求参数μ的MLE ;又若总体为()21N σ,的MLE 。
第二章多元正态分布的参数估计
第二章多元正态分布的参数估计多元正态分布是在多个随机变量之间存在相互依赖关系时使用的一种概率分布。
它在许多统计分析和机器学习领域中都有广泛的应用。
在实际应用中,我们通常需要使用样本数据对多元正态分布的参数进行估计。
多元正态分布由均值向量和协方差矩阵两个参数来描述。
均值向量表示各个随机变量的平均值,而协方差矩阵表示各个随机变量之间的协方差。
参数估计的目标就是通过样本数据来估计这两个参数。
首先,我们需要收集一个具有充分样本量的数据集。
对于一个具有n个样本的多元正态分布,我们可以将样本数据表示为一个n行d列的矩阵X,其中每一行是一个d维的样本向量。
其中n表示样本数量,d表示随机变量的个数。
接下来,我们可以根据样本数据来估计多元正态分布的均值向量和协方差矩阵。
1.均值向量的估计:多元正态分布的均值向量可以通过样本均值向量来估计。
样本均值向量的计算公式如下:μ = (1/n) * Σxi其中μ是估计得到的均值向量,xi表示样本矩阵X的第i行。
2.协方差矩阵的估计:多元正态分布的协方差矩阵可以通过样本协方差矩阵来估计。
Σ=(1/(n-1))*(X-μ)'*(X-μ)其中Σ是估计得到的协方差矩阵,X是样本矩阵,μ是估计得到的均值向量。
需要注意的是,在计算协方差矩阵时,我们使用的是样本协方差矩阵而不是总体协方差矩阵。
这是因为样本协方差矩阵能更好地反映样本数据的真实情况。
以上就是多元正态分布的参数估计方法。
通过样本数据,我们可以使用样本均值向量和样本协方差矩阵来估计多元正态分布的参数。
这些参数估计能为我们提供关于多元正态分布的统计属性和特征,进而用于进一步的分析和应用。
应用数理统计(武汉理工大)2-参数估计
1
D(S 2 )nI (
2)
n 1 n
1,
n
故S 2是渐进有效的。
第二章 参数估计
例: 设总体X (), X1, X 2 , , X n是X的一个样本, 讨论的无偏估计X的有效性。
解:lnp( X
,)
ln
X e
X!
X
ln
ln( X
!)
区间估计的关键: 用合适的方法确定两个统计量
1(X1, X2 , , Xn), 2(X1, X2 , , Xn)
第二章 参数估计
1.区间估计的定义及计算步骤
3) 区间估计的例子
例1 设总体X~N(μ , σ2), σ2已知,μ未知,设X1,…,Xn是X的样本, 求μ的置信度为1-α的置信区间。
)
2
n
,
D(ˆ2 )
D(nZ )
n2D(Z )
n2
n
2
2
当n 1时,显然D(ˆ1) D(ˆ2 ),故ˆ1比ˆ2有效。
第二章 参数估计
最小方差无偏估计问题 设 若 及T对 任(g意X(1, , X)的2都,任有一 , XD无n()T是 偏) g估(D计()T的量')一, T '个 ( X无1, X偏2估 , 计, X量n ), 则 无称 偏T估(计X1,, X或2 ,者,称X为n )是最g优(无)的偏一估致计最。小方差
其它类型的估计,如 贝叶斯估计…
第二章 参数估计
2.1参数的点估计
1. 矩估计 2. 极大似然估计 3. 点估计量的评价
第二章 参数估计
1 n k Xi 的一个样本,则 Ak n i 1
是 的无偏估
k
例 2.3:设总体 X 的数学期望 与方差 存在,
75 90 54 22 6 2 1
1.3 极大似然估计 1.3.1 极大似然法所依据的基本原理是: 已经 实现的事件是概率最大的事件,或者说,概 率最大的事件最可能出现。
例如:1. 有两位同学一起进行实弹射击, 两人共同射击一个目 标,事先并不知道谁的技术较好,让每人各打一发,有一人 击中目标,认为击中的同学的技术比击不中的技术要好 2.一个老猎人常领一个新手进山打猎,遇见一只飞奔的野兔, 他们各发一弹,野兔被打中了,但身上只有一个弹孔,到底 是谁打中的呢?一般,绝大多数人认为是老猎人打中。 3.医生看病,在问明症状后(包括必要的检查) ,做诊断的总 是对那些可能直接引起这些症状的疾病多加考虑的; 4.当机器发生故障,有经验的修理工,首先总从易损部件,薄 弱环节查起; 5.公安人员在破一起他杀案时, 谁是凶手?首先把与被害者密 切来往又有作案可能性的人列为重点嫌疑对象。 6.有一事件,知道它发生的概率 p 只可能是 0.01 或 0.09,再 一次观察中这事件发生了,试问这事件发生的概率是什么,
例 1.5:设总体 X 服从 —分布 , ,其密 度函数为,
1 f x x e ,x 0 1
x
求 , 的矩法估计量。 例1. 6:求相关系数 的矩法估计量。 例 1.7: 某炸药制造厂, 一天中发生的着火现 象的次数 X~ P , 未知:试估计 着火次数 着火 k 次的天数 0 1 2 3 456
1 n
a,b 的极大似然估计。
例 1.15 : 设 总 体 X 的 概 率 密 度 为
第二章 参数估计
ˆ = q ( X , K , X ) , q k k 1 n
k = 1, 2, L , m
(2.2)
ˆ 为 q 的矩估计, g ( x 若 q ) 为连续函数,则也称 g (qˆ k k k ) 为 g (q k ) 的矩估计.
【例 2.1】 设总体 X 服从参数为 l 的泊松分布,X 1 , K , X n 为来自总体的样本, 求l 的 矩估计. 解: a1 = EX = l
i =1
定义 2.1:设总体 X 的概率函数为 f ( x;q ) , x1 ,L , x n 是来自总体的样本,则称
n
L(q ) = Õ f ( xi ;q )
i =1
(2.4)
为总体 X 对应样本 x1 ,L , x n 的似然函数.
L(q ) 越大,越有利于样本 x1 ,K , x n 被观察到.
-l ì l x e ï f ( x 0,1, 2, L 其它
或简写为
f ( x) =
-l l x e
x !
x = 0,1, 2, L
§2.1 点估计
我们经常会遇到这样的问题: 总体 X 的分布函数 F ( x,q ) 的形式已知, 但其中的参数q 未知, 希望利用 X 的样本 x1 ,K , x 这类问题称为参数的点估计 (point n 对 q 的值进行估计, estimation)问题. 比如,已知某种电子元件的寿命 X ~ N ( m , s ) ,即 X 的分布密度
P( X = xi ) = p( xi ,q ), i = 1, 2,L ,
其中q 为未知参数,q Î Q . 设 X 1 , K , X n 是来自总体 X 的一组样本, 观察值为 x1 ,K , x n .我们把观察到的样本看成 结果,而需要判断的是未知参数q 的取值,根据最大似然原理,应该选取一个最有利于结 果的发生的q 值作为 qˆ .
应用数理统计吴翊李永乐第二章-参数估计课后习题参考答案
《应用数理统计》吴翊李永乐第二章-参数估计课后习题参考答案(总19页)--本页仅作为文档封面,使用时请直接删除即可----内页可以根据需求调整合适字体及大小--第二章 参数估计课后习题参考答案设总体X 服从二项分布()n X X X p p N B ,,,,11,,21 <<为其子样,求N 及p 的矩法估计。
解:()()()p Np X D Np X E -==1,令()⎪⎩⎪⎨⎧-==p Np S Np X 12 解上述关于N 、p 的方程得:对容量为n 的子样,对密度函数22(),0(;)0,0x x f x x x ααααα⎧-⎪=⎨⎪≤≥⎩其中参数α的矩法估计。
解:122()()a E x xx dx ααα==-⎰2222()x x dx ααα=-⎰2321221333ααααααα=-=-= 所以 133a x α∧== 其中121,21(),,,n n x x x x x x x n=+++为n 个样本的观察值。
使用一测量仪器对同一值进行了12次独立测量,其结果为(单位:mm) ,,,,,⎪⎪⎩⎪⎪⎨⎧-=-==X S p S X X p X N 2221ˆˆˆ,,,,,试用矩法估计测量的真值和方差(设仪器无系统差)。
解:()()()∑∑====-====ni ini i S XX nX D X X n X E 12210255.014025.2321设子样,,,,,是来自具有密度函数()10,1,<<=βββx f 的总体,试用矩法估计总体均值、总体方差及参数β。
解:()()()()4.22ˆ2,1,407.012.1101221========-===⎰⎰∑∑==X Xdx xdx x xf X E x f XX n S X n X ni i ni i ββββββββ参数:总体方差:总体均值:设n X X X ,,,21 为()1N ,μ的一个字样,求参数μ的MLE ;又若总体为()21N σ,的MLE 。
第二章参数估计理论_3
ˆ 估计准则为使均方误差E{(θ − θ ) 2 }最小。即 ˆ arg min E{(θ − θ ) } = E{(∑ wi xi − θ ) 2 }=E{e 2 }
2 wi i =1 N
求偏导数
∂E{e 2 } ∂e 2 ∂e = E{ } = 2 E{e } = 2 E{exi } = 0 ∂wi ∂wi ∂wi i = 1,..., N
物理含义:在上述观测模型中,误差向量e各分量间不相关 并且具有相同的方差时,最小二乘估计才是无偏的和具有最 小方差的估计。
且系数矩阵列满秩时
最小二乘(LS)估计
加权最小二乘估计(为了克服误差分量具有相关性或方差不同问题)
ˆ 目标:使损失函数Q(θ )=e H We最小 ˆ ˆ θˆ = arg min e H We = e H e = ( Aθ − b) H W (Aθ − b)
最小二乘(LS)估计 如何求加权矩阵W?
假设误差向量的自协方差 cov(e ) = σ 2V ,V 是一个已知正定矩阵。 由于V 正定 ⇒ V = PP H , 其中P非奇异。 用P -1左乘观测模型b = Aθ + e得
P -1b = P -1 Aθ + P -1e
⇒ x = Bθ + ε 其中x = P -1b, B = P -1 A, ε = P -1e 而 cov(ε ) = cov( P -1e ) = σ 2 I 即新的误差向量ε = P -1e各分量间不相关,且具有相同方差。 ˆ ∴ θ = ( B H B )-1 B H x
第二章 参数估计理论 (3) 线性均方估计(LMMSE) 最小二乘估计(LS)
2.5 线性均方估计
贝叶斯估计需要知道待估计量的PDF或者后验分布函数f(θ|x) 最大似然估计会导致非线性问题的求解。
第二章-多元正态分布的参数估计
11 Σ 21
31
12 22 32
13 23 33
Σ11
Σ
21
Σ12
22
则
X (1)
X1
X
2
~
N2 ( μ(1) ,
Σ11)
其中
μ (1)
1
2
Σ11
11 21
12
22
在此我们应该注意到,如果 X ( X1, X 2 ,L , X p ) 服从 p
X
X1 X2
,
μ
1 2
,
Σ
2 1
1 2
1 2
2 2
易见,ρ是X1和 X2的相关系数。当|ρ|<1时,可得X的
概率密度函数为:
f
x1,
x2
1
21 2
1
2
exp 2
1
1 2
x1 1 1
2
2
x1 1 1
x2 2 2
x2 2 2
2
二元正态分布的密度曲面图
X3
1
aμ
(0,1,
0)
2
2
3
11 12 aΣa (0,1, 0) 21 22
31 32
13 0
23
1
22
33 0
(2) 其中
AX
1
0
0 0
0 1
X X X
1 2 3
X1
X
3
~
N
(Aμ
,AΣA
)
Aμ
1 0
0 0
0 1
1 2 3
5 1
11
则X2和X3不独立,X1和(X2,X3)独立。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
arg min wi W
E{e2}
e
令 E{e2} 0,得 wk
E{exi} 0,i 1,..., N ,即
E
N
wk xk
xi
0,
i
1,
...,
N
k1
令gi E{ xi}, Rij E{xi x j }, 则有
试确定均值和方差的最大似然估计。
解:由题设可得
N
f (x1,..., xN | , 2 ) i 1
1
e( xi )2 /(2 2 )
2
(2
2
)
N
/
2
exp
(
1
2
2
N i 1
( xi
)2
)
从而有
L( ) ln[ f (x1,..., xN | , )]
R( , ) E{C(, )}
使二次风险函数最小的估计称为最小均方误差(MMSE)估计:
RMMSE E{( )2}
...
(ˆ
)2
f
( x1 , ...,
xN
, )dx1...dxN d
...
[
(
)2
评价两个估计子有效性的量度:均值和方差。
均值:b( ) E{ } E{}
方差: var( ) E{( E{})2}
均值和方差:M 2 ( ) E{( )2} var( ) b2 ( )
2.2 Fisher信息与Cramer-Rao不等式
损失函数 :
C(, ) | |,绝对损失函数
C( , ) | |2,二次型损失函数
C( ,
)
0,|
|
均匀损失函数
1,| |
损失函数是随机变量x的函数,因而本身也是随机的。
损失函数的数学期望称为风险函数:
N
Rik wk gi , i 1,..., N
k 1
采用向量表示可得
Rw g
线性均方估计假定随机变量的相关函数已知,其本质上 就是线性滤波器。
2.6 最小二乘估计
2.6.1 最小二乘估计
问题描述:b A e
N
LS
arg min
i 1
ei2
arg min( A b)T ( A b)
如何判断一个估计子是否是最优的。
随机变量x的品质函数
V (x)
ห้องสมุดไป่ตู้
ln f (x | )
f (x |)
f (x |)
品质函数的性质: E{V (x)} 0, var(V (x)) E{V 2 (x)}
品质函数的方差越大越好,还是越小越好?
品质函数的方差称为Fisher信息,表征样本中包含被估计 量的信息量:
f (
| x1,..., xN )d ]
令 RMM SfE(x10..,., x得N )dx1...dxN
ˆMMSE
f
(
|
x1,...,
xN )d
E{
|
x1,...xN }
2.4 最大似然估计
问似题然描函述数:的定M义L :aLr(gm)axlnf (fx1(,x..1.,,.x..N,
J ( )
E
ln
f
(x
|
)
2
Fisher信息值越大,意味着相邻估计子的区分度越高。
Cramer-Rao不等式:令x=(x1,…,xN)为样本向量。若参数
估计
是真实参数
的无偏估计,且
f
(x |)和
2
f
(x |)
2
存在,则 的均方误差所能达到的下界等于Fisher信息的倒
第二章 参数估计理论
估计理论的分类:参数化与非参数化(基于模型与无模型)。 参数化方法假定数据服从已知结构的概率模型,但是模型的 某些参数未知。例如线性均方估计,ARMA谱估计。 非参数化方法不假定数据服从某种特定的概率模型。例如基 于离散傅立叶变换的功率谱估计和高阶谱估计。 参数化方法的效率高,但是适应范围窄。
数,即
var( ) E{( E{})2} E{( )2}
1
J ( )
2.3 Bayes估计
Bayes公式:
f ( | x) f (x | ) f ( ) , p( | x) p(x | ) p( )
f (x)
p(x)
意义:结果到原因(推断)原因到结果(机理)
2.1 估计子的性能
2.1.1 无偏估计与渐进无偏估计
参数估计:由N个样本数据推出p个参数的值。
估计子:真实参数的近似。
参数估计的偏差:b( ) E{ } E{}
无偏估计:b( ) 0
渐进无偏估计:lim b( N ) 0
N
2.1.2 估计子的有效性
N 2
ln(2 )
N 2
ln(
2
)
1
2
2
N
( xi
i 1
)2
分别求偏导,可得
L
1
2
N
(xi ) 0
i 1
L
2
N
2 2
1
2 4
N
(xi )2
i 1
0
解以上二方程,可得
ML
1 N
N
xi
i 1
_
x
N
代价函数J ei2 ( A b)T ( A b)
i 1
令
dJ
2 AT A 2 AT b 0,即
|)
xN |
)
特别地,当x1,…,xN为独立的观测样本,N似然函数为:
L( ) ln[ f (x1 | )... f (xN | )] ln f (xi | )
令
L(
)
0,可求出
ML
i 1
例2.4.1 x1,…,xN为独立的观测样本,且服从正态分布
f ( x, , 2 ) 1 e( x )2 /(2 2 ) 2
与 样ML本均N1值和iN1方(差xi 的无x_ )偏2 估计比较,可知均值的估计是无
偏的,方差的估计是渐进无偏的。
2.5 线性均方估计
N
问题描述: ˆLMS
wi xi
i 1
wi
arg min wi W
E
N i 1
wi xi
2