跳到主要内容

京都大学 情報学研究科 システム科学専攻 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,Z を A,ϕA,\phi で表せ。

(ii) AA の確率密度が pA(x)=xe−x2/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⁡(Yn≤x)=Pr⁡(X1≤x,…,Xn≤x)\Pr(Y_n\le x)=\Pr(X_1\le x,\ldots,X_n\le x) を示せ。

(ii) YnY_n の密度が 0<x<θ0<x<\theta で pYn(x)=nxn−1/θ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 t。tt 是时间,ω\omega 为角频率。

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

(ii) 若 pA(x)=xe−x2/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⁡(Yn≤x)=Pr⁡(X1≤x,…,Xn≤x)\Pr(Y_n\le x)=\Pr(X_1\le x,\ldots,X_n\le x)。

(ii) 证明 YnY_n 的密度在 0<x<θ0<x<\theta 上为 nxn−1/θ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\phi、Z=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 は独立な標準正規変数であり、sin⁡2ωt+cos⁡2ω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/2e−x2/2p_{X(t)}(x)=(2\pi)^{-1/2}e^{-x^2/2}。

(iii) 位相について平均すると E[sin⁡(ωt+ϕ)]=0E[\sin(\omega t+\phi)]=0、E[sin⁡2(ωt+ϕ)]=1/2E[\sin^2(\omega t+\phi)]=1/2。E[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<\theta で Pr⁡(Yn≤x)=(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(c−n+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)}}.