跳到主要内容

京都大学 情報学研究科 システム科学専攻 2008年8月実施 確率統計

Author

犬 (finalized by 祭音Myyura with assistance from GPT 6 Astra)

Description

問1

独立なランダム振幅 AA と位相 ϕ\phi をもつ X(t)=Asin(ωt+ϕ)X(t)=A\sin(\omega t+\phi) を考える。tt は時間、ω\omega は角周波数であり、X(t)=Ysinωt+ZcosωtX(t)=Y\sin\omega t+Z\cos\omega t によって Y,ZY,Z を定義する。

(i) Y,ZY,ZA,ϕA,\phi で表せ。

(ii) AA の確率密度が pA(x)=xex2/2p_A(x)=xe^{-x^2/2}x>0x>0)で、ϕ\phi(0,2π)(0,2\pi) の一様分布に従うとき、Y,ZY,Z の同時確率密度および X(t)X(t) の確率密度を求めよ。

(iii) AA(0,1)(0,1) の一様分布、ϕ\phi(0,2π)(0,2\pi) の一様分布に従うとき、E[X(t)],E[X(t)2]E[X(t)],E[X(t)^2] を求めよ。

問2

X1,,XnX_1,\ldots,X_n(0,θ)(0,\theta) の一様分布からの独立標本とし、Yn=max(X1,,Xn)Y_n=\max(X_1,\ldots,X_n) により未知母数 θ\theta を推定する。

(i) Pr(Ynx)=Pr(X1x,,Xnx)\Pr(Y_n\le x)=\Pr(X_1\le x,\ldots,X_n\le x) を示せ。

(ii) YnY_n の密度が 0<x<θ0<x<\thetapYn(x)=nxn1/θnp_{Y_n}(x)=nx^{n-1}/\theta^n、その他で 00 となることを示せ。

(iii) E[Yn]E[Y_n] を求め、不偏推定量かどうか述べよ。

(iv) E[Yn2]E[Y_n^2] と平均二乗推定誤差 E[(Ynθ)2]E[(Y_n-\theta)^2] を求めよ。

(v) cnYnc_nY_n の平均二乗推定誤差を最小にする定数 cnc_n を求め、YnY_n に比べた改善量を述べよ。

题目描述

问1

独立随机振幅 AA 和相位 ϕ\phi 产生信号 X(t)=Asin(ωt+ϕ)=Ysinωt+ZcosωtX(t)=A\sin(\omega t+\phi)=Y\sin\omega t+Z\cos\omega ttt 是时间,ω\omega 为角频率。

(i) 用 A,ϕA,\phi 表示 Y,ZY,Z

(ii) 若 pA(x)=xex2/2p_A(x)=xe^{-x^2/2}x>0x>0),ϕ\phi(0,2π)(0,2\pi) 上均匀分布,求 Y,ZY,Z 的联合密度和 X(t)X(t) 的密度。

(iii) 若 AA(0,1)(0,1) 上均匀分布,ϕ\phi(0,2π)(0,2\pi) 上均匀分布,求 E[X(t)]E[X(t)]E[X(t)2]E[X(t)^2]

问2

X1,,XnX_1,\ldots,X_n 独立服从 (0,θ)(0,\theta) 的均匀分布,以最大值 YnY_n 估计未知参数 θ\theta

(i) 证明 Pr(Ynx)=Pr(X1x,,Xnx)\Pr(Y_n\le x)=\Pr(X_1\le x,\ldots,X_n\le x)

(ii) 证明 YnY_n 的密度在 0<x<θ0<x<\theta 上为 nxn1/θnnx^{n-1}/\theta^n,其他为 00

(iii) 求 E[Yn]E[Y_n],判断是否无偏。

(iv) 求 E[Yn2]E[Y_n^2]E[(Ynθ)2]E[(Y_n-\theta)^2]

(v) 求使 cnYnc_nY_n 的均方误差最小的常数 cnc_n,并说明相对 YnY_n 的误差改善量。

Kai

問1

(i) 加法定理より Y=AcosϕY=A\cos\phiZ=AsinϕZ=A\sin\phi

(ii) 極座標変換のヤコビアンは (y,z)/(a,ϕ)=a|\partial(y,z)/\partial(a,\phi)|=a。よって

pY,Z(y,z)=12πe(y2+z2)/2.p_{Y,Z}(y,z)=\frac1{2\pi}e^{-(y^2+z^2)/2}.

Y,ZY,Z は独立な標準正規変数であり、sin2ωt+cos2ωt=1\sin^2\omega t+\cos^2\omega t=1 なので X(t)N(0,1)X(t)\sim N(0,1)。したがって pX(t)(x)=(2π)1/2ex2/2p_{X(t)}(x)=(2\pi)^{-1/2}e^{-x^2/2}

(iii) 位相について平均すると E[sin(ωt+ϕ)]=0E[\sin(\omega t+\phi)]=0E[sin2(ωt+ϕ)]=1/2E[\sin^2(\omega t+\phi)]=1/2E[A2]=1/3E[A^2]=1/3 より

E[X(t)]=0,E[X(t)2]=16.\boxed{E[X(t)]=0,\qquad E[X(t)^2]=\frac16}.

問2

(i) 最大値が xx 以下である事象は、全標本が xx 以下である事象に等しい。

(ii) 独立性より 0<x<θ0<x<\thetaPr(Ynx)=(x/θ)n\Pr(Y_n\le x)=(x/\theta)^n。微分すると所望の密度を得る。

(iii)、(iv) 密度から一般に E[Ynk]=nθk/(n+k)E[Y_n^k]=n\theta^k/(n+k)k>0k>0)なので

E[Yn]=nn+1θ,E[Yn2]=nn+2θ2,E[(Ynθ)2]=2θ2(n+1)(n+2).E[Y_n]=\frac n{n+1}\theta,\quad E[Y_n^2]=\frac n{n+2}\theta^2, \quad E[(Y_n-\theta)^2]=\frac{2\theta^2}{(n+1)(n+2)}.

特に YnY_n は不偏ではない。

(v) cc に関する二次式を平方完成すると

E[(cYnθ)2]=nθ2n+2(cn+2n+1)2+θ2(n+1)2.E[(cY_n-\theta)^2] =\frac{n\theta^2}{n+2}\left(c-\frac{n+2}{n+1}\right)^2 +\frac{\theta^2}{(n+1)^2}.

よって cn=(n+2)/(n+1)\boxed{c_n=(n+2)/(n+1)}。誤差の減少量は

2θ2(n+1)(n+2)θ2(n+1)2=nθ2(n+1)2(n+2).\frac{2\theta^2}{(n+1)(n+2)}-\frac{\theta^2}{(n+1)^2} =\boxed{\frac{n\theta^2}{(n+1)^2(n+2)}}.