跳到主要内容

九州大学 経済学府 経済工学専攻 2019年8月実施 経済数学

Author​

Miyake, 祭音Myyura

Description​

次の2問の両方について解答せよ。

問1​

(1)、(2) の中から1つを選んで解答せよ。

(1)​

3次元列ベクトル x1,x2,x3\boldsymbol x_1,\boldsymbol x_2,\boldsymbol x_3 および 03\boldsymbol0_3 を次のように定義する。

x1=(201),x2=(909),x3=(1t8),03=(000).\boldsymbol x_1=\begin{pmatrix}2\\0\\1\end{pmatrix},\quad \boldsymbol x_2=\begin{pmatrix}9\\0\\9\end{pmatrix},\quad \boldsymbol x_3=\begin{pmatrix}1\\t\\8\end{pmatrix},\quad \boldsymbol0_3=\begin{pmatrix}0\\0\\0\end{pmatrix}.

ただし、tt は実数である。また、3×33\times3 行列 XX を X=(x1,x2,x3)X=(\boldsymbol x_1,\boldsymbol x_2,\boldsymbol x_3) によって定義する。

(a) x1,x2,x3\boldsymbol x_1,\boldsymbol x_2,\boldsymbol x_3 が1次独立になるのは tt がどのような値のときか答えよ。

(b) c\boldsymbol c を3次元列ベクトルとする。XTX^{\mathsf T} が正則であるかどうかにかかわらず、XTXc=03X^{\mathsf T}X\boldsymbol c=\boldsymbol0_3 ならば、Xc=03X\boldsymbol c=\boldsymbol0_3 であることを示せ。ただし、右肩の添え字 T\mathsf T は転置を意味する。

(c) 行列 XTXX^{\mathsf T}X が正則になるのは tt がどのような値のときか答えよ。

(2)​

次の問いに答えよ。ただし、aa は正定数とする。

(a) 3次元ユークリッド空間 R3\mathbb R^3 内の xzxz 平面上の曲線 z=(x−a)2z=(x-a)^2(x≥0x\ge0)を zz 軸の周りに回転させた曲面が次のように書けることを示せ。

z=(x2+y2−a)2.z=\left(\sqrt{x^2+y^2}-a\right)^2.

(b) 次の関数 ff の極値を求めよ。

f(x,y)=(x2+y2−a)2.f(x,y)=\left(\sqrt{x^2+y^2}-a\right)^2.

問2​

(1)、(2) の中から1つを選んで解答せよ。

(1)​

三角形 ABCABC において、3つの辺 BC,CA,ABBC,CA,AB の長さをそれぞれ a,b,ca,b,c とし、三角形の内部の点 DD から、3本の直線 BC,CA,ABBC,CA,AB へ下した垂線の長さをそれぞれ x,y,zx,y,z とする。この三角形の面積を SS とおく。

(a) xx と yy の取りうる範囲を求めよ。

(b) x2+y2+z2x^2+y^2+z^2 の極小値を求めよ。

(c) (b) で求めた極小値が最小値になっていることを示せ。

(d) xyzxyz の最大値を求めよ。

(2)​

nn を自然数とする。また、p∈(0,1)p\in(0,1)、t∈Rt\in\mathbb R および λ>0\lambda>0 を定数とする。

(a) 確率変数 XX は2項分布 Bin⁡(n,p)\operatorname{Bin}(n,p) に従っているとする。確率関数は

pX(x)=nCxpx(1−p)n−x,x=0,1,…,np_X(x)={}_nC_xp^x(1-p)^{n-x},\qquad x=0,1,\ldots,n

である。期待値 E[exp⁡(tX)]E[\exp(tX)] を求めよ。

(b) 確率変数 YY はポアソン分布 Po⁡(λ)\operatorname{Po}(\lambda) に従っているとする。確率関数は

pY(y)=exp⁡(−λ)λyy!,y=0,1,…p_Y(y)=\exp(-\lambda)\frac{\lambda^y}{y!},\qquad y=0,1,\ldots

である。期待値 E[exp⁡(tY)]E[\exp(tY)] を求めよ。

(c) 確率変数 ZnZ_n は2項分布 Bin⁡(n,λ/n)\operatorname{Bin}(n,\lambda/n) に従っているとする。ただし、n>λn>\lambda とする。極限 lim⁡n→∞E[exp⁡(tZn)]\lim_{n\to\infty}E[\exp(tZ_n)] を求めよ。

题目描述​

以下两大题都必须作答。

问1 在 (1)、(2) 中选一题作答。

(1) 定义三维列向量

x1=(201),x2=(909),x3=(1t8),03=(000).\boldsymbol x_1=\begin{pmatrix}2\\0\\1\end{pmatrix},\quad \boldsymbol x_2=\begin{pmatrix}9\\0\\9\end{pmatrix},\quad \boldsymbol x_3=\begin{pmatrix}1\\t\\8\end{pmatrix},\quad \boldsymbol0_3=\begin{pmatrix}0\\0\\0\end{pmatrix}.

其中 tt 为实数,定义 3×33\times3 矩阵 X=(x1,x2,x3)X=(\boldsymbol x_1,\boldsymbol x_2,\boldsymbol x_3)。

(a) 求使 x1,x2,x3\boldsymbol x_1,\boldsymbol x_2,\boldsymbol x_3 线性无关的 tt 的取值。

(b) 设 c\boldsymbol c 为三维列向量。无论 XTX^{\mathsf T} 是否可逆,证明若 XTXc=03X^{\mathsf T}X\boldsymbol c=\boldsymbol0_3,则 Xc=03X\boldsymbol c=\boldsymbol0_3。右上标 T\mathsf T 表示转置。

(c) 求使矩阵 XTXX^{\mathsf T}X 可逆的 tt 的取值。

(2) 回答下列问题,aa 为正常数。

(a) 将三维欧氏空间 R3\mathbb R^3 中 xzxz 平面上的曲线 z=(x−a)2z=(x-a)^2(x≥0x\ge0)绕 zz 轴旋转,证明所得曲面可写为

z=(x2+y2−a)2.z=\left(\sqrt{x^2+y^2}-a\right)^2.

(b) 求函数 f(x,y)=(x2+y2−a)2f(x,y)=(\sqrt{x^2+y^2}-a)^2 的极值。

问2 在 (1)、(2) 中选一题作答。

(1) 三角形 ABCABC 的边 BC,CA,ABBC,CA,AB 的长度分别为 a,b,ca,b,c。从三角形内部的点 DD 向三条直线 BC,CA,ABBC,CA,AB 所作垂线的长度分别为 x,y,zx,y,z,三角形面积记为 SS。

(a) 求 x,yx,y 的取值范围。

(b) 求 x2+y2+z2x^2+y^2+z^2 的极小值。

(c) 证明 (b) 中求出的极小值也是最小值。

(d) 求 xyzxyz 的最大值。

(2) 设 nn 为自然数,p∈(0,1)p\in(0,1)、t∈Rt\in\mathbb R、λ>0\lambda>0 为常数。

(a) 随机变量 XX 服从二项分布 Bin⁡(n,p)\operatorname{Bin}(n,p),概率质量函数为

pX(x)=nCxpx(1−p)n−x,x=0,1,…,n.p_X(x)={}_nC_xp^x(1-p)^{n-x},\qquad x=0,1,\ldots,n.

求期望 E[exp⁡(tX)]E[\exp(tX)]。

(b) 随机变量 YY 服从泊松分布 Po⁡(λ)\operatorname{Po}(\lambda),概率质量函数为

pY(y)=exp⁡(−λ)λyy!,y=0,1,….p_Y(y)=\exp(-\lambda)\frac{\lambda^y}{y!},\qquad y=0,1,\ldots.

求期望 E[exp⁡(tY)]E[\exp(tY)]。

(c) 随机变量 ZnZ_n 服从二项分布 Bin⁡(n,λ/n)\operatorname{Bin}(n,\lambda/n),其中 n>λn>\lambda。求极限 lim⁡n→∞E[exp⁡(tZn)]\lim_{n\to\infty}E[\exp(tZ_n)]。

Kai​

問 1​

(1)​

(a)​

det⁡X=−9t\det X = -9t なので、 x1,x2,x3\boldsymbol{x}_1, \boldsymbol{x}_2, \boldsymbol{x}_3 が1次独立になるのは、 t≠0t \ne 0 のときである。

(b)​
c=(abc)\begin{aligned} \boldsymbol{c} = \begin{pmatrix} a \\ b \\ c \end{pmatrix} \end{aligned}

とすると、

Xc=ax1+bx2+cx3\begin{aligned} X \boldsymbol{c} = a \boldsymbol{x}_1 + b \boldsymbol{x}_2 + c \boldsymbol{x}_3 \end{aligned}

なので、 XcX \boldsymbol{c} は x1,x2,x3\boldsymbol{x}_1, \boldsymbol{x}_2, \boldsymbol{x}_3 の1次結合で表される。

また、

XTXc=(x1T(Xc)x2T(Xc)x3T(Xc))\begin{aligned} X^T X \boldsymbol{c} = \begin{pmatrix} \boldsymbol{x}_1^T (X \boldsymbol{c}) \\ \boldsymbol{x}_2^T (X \boldsymbol{c}) \\ \boldsymbol{x}_3^T (X \boldsymbol{c}) \end{pmatrix} \end{aligned}

と表されるので、 XTXc=03X^T X \boldsymbol{c} = \boldsymbol{0}_3 は、 XcX \boldsymbol{c} が x1,x2,x3\boldsymbol{x}_1, \boldsymbol{x}_2, \boldsymbol{x}_3 のいずれとの内積も 00 であることを意味する。

したがって、 XTXc=03X^T X \boldsymbol{c} = \boldsymbol{0}_3 ならば、 Xc=03X \boldsymbol{c} = \boldsymbol{0}_3 である。

(cc)​
det⁡(XTX)=(det⁡XT)(det⁡X)=(det⁡X)2=81t2 \begin{aligned} \det \left( X^T X \right) &= \left( \det X^T \right) \left( \det X \right) \\ &= \left( \det X \right)^2 \\ &= 81t^2 \end{aligned}

なので、 XTXX^T X が正則なのは t≠0t \ne 0 のときである。

(2)​

(a)​

回転前の座標を ρ≥0\rho\ge0 とすると、回転後の点は

(x,y,z)=(ρcos⁡θ,ρsin⁡θ,(ρ−a)2),0≤θ<2π.\begin{gathered} (x,y,z)=(\rho\cos\theta,\rho\sin\theta,(\rho-a)^2),\\ 0\le\theta<2\pi. \end{gathered}

ρ=x2+y2\rho=\sqrt{x^2+y^2} より、曲面の式は z=(x2+y2−a)2z=(\sqrt{x^2+y^2}-a)^2 となる。

(b)​

r=x2+y2r=\sqrt{x^2+y^2} とおくと f=(r−a)2≥0f=(r-a)^2\ge0 である。 r=ar=a、すなわち円 x2+y2=a2x^2+y^2=a^2 上の各点で最小値 00 を取る。

原点では f(0,0)=a2f(0,0)=a^2 であり、0<r<2a0<r<2a ならば

f(x,y)−a2=r(r−2a)<0.f(x,y)-a^2=r(r-2a)<0.

したがって原点は厳密な極大点で、極大値は a2a^2 である。 r>0r>0 では ∇f=2(r−a)(x,y)/r\nabla f=2(r-a)(x,y)/r なので、他の極値はない。 なお r→∞r\to\infty で f→∞f\to\infty のため、大域的な最大値は存在しない。

問 2​

(1)​

(a)​

DD と3辺でできる三角形の面積を足すと、

ax+by+cz=2S.ax+by+cz=2S.

DD は内部の点なので x,y,z>0x,y,z>0 であり、

x>0,y>0,ax+by<2S.\boxed{x>0,\qquad y>0,\qquad ax+by<2S.}

逆に、この条件を満たす (x,y)(x,y) に対して z=(2S−ax−by)/c>0z=(2S-ax-by)/c>0 とおけば、正の重心座標 (ax,by,cz)/(2S)(ax,by,cz)/(2S) が和 11 を持ち、対応する内部の点 DD が存在する。

(b)​

制約 ax+by+cz=2Sax+by+cz=2S の下で Lagrange の未定乗数法を使うと 2x=μa, 2y=μb, 2z=μc2x=\mu a,\ 2y=\mu b,\ 2z=\mu c となる。 したがって

(x,y,z)=2Sa2+b2+c2(a,b,c),(x,y,z)=\frac{2S}{a^2+b^2+c^2}(a,b,c),

における極小値は

4S2a2+b2+c2.\boxed{\frac{4S^2}{a^2+b^2+c^2}}.
(c)​

Cauchy–Schwarz の不等式から

(2S)2=(ax+by+cz)2≤(a2+b2+c2)(x2+y2+z2).(2S)^2=(ax+by+cz)^2 \le(a^2+b^2+c^2)(x^2+y^2+z^2).

(b) の正の距離の組で等号が成り立つため、これは最小値である。

(d)​

ax,by,cz>0ax,by,cz>0 に相加相乗平均の不等式を適用すると、

abc xyz≤(ax+by+cz3)3=(2S3)3.abc\,xyz\le\left(\frac{ax+by+cz}{3}\right)^3 =\left(\frac{2S}{3}\right)^3.

ax=by=cz=2S/3ax=by=cz=2S/3 のとき等号が成り立つ。この点は三角形の重心であり、

max⁡xyz=8S327abc.\boxed{\max xyz=\frac{8S^3}{27abc}}.

(2)​

(a)​
E[exp⁡(tX)]=∑x=0nexp⁡(tx)⋅nCxpx(1−p)n−x=∑x=0n nCx(pet)x(1−p)n−x=(1−p+pet)n \begin{aligned} E \left[ \exp (tX) \right] &= \sum_{x=0}^n \exp(tx) \cdot {}_n C_x p^x (1-p)^{n-x} \\ &= \sum_{x=0}^n \ {}_n C_x \left( p e^t \right)^x (1-p)^{n-x} \\ &= \left( 1 - p + p e^t \right)^n \end{aligned}
(b)​
E[exp⁡(tY)]=∑y=0∞exp⁡(ty)⋅exp⁡(−λ)λyy!=exp⁡(−λ)∑y=0∞(λet)yy!=exp⁡(−λ)⋅exp⁡(λet)=exp⁡(λ(et−1)) \begin{aligned} E \left[ \exp (tY) \right] &= \sum_{y=0}^\infty \exp(ty) \cdot \exp(- \lambda) \frac{\lambda^y}{y!} \\ &= \exp(- \lambda) \sum_{y=0}^\infty \frac{\left( \lambda e^t \right)^y}{y!} \\ &= \exp(- \lambda) \cdot \exp \left( \lambda e^t \right) \\ &= \exp \left( \lambda \left( e^t - 1 \right) \right) \end{aligned}
(cc)​

(a) より、

E[exp⁡(tZn)]=(1+λ(et−1)n)n\begin{aligned} E \left[ \exp \left( tZ_n \right) \right] &= \left( 1 + \frac{\lambda \left( e^t - 1 \right)}{n} \right)^n \end{aligned}

なので、

lim⁡n→∞E[exp⁡(tZn)]=exp⁡(λ(et−1))\begin{aligned} \lim_{n \to \infty} E \left[ \exp \left( tZ_n \right) \right] &= \exp \left( \lambda \left( e^t - 1 \right) \right) \end{aligned}

である。右辺は Poisson⁡(λ)\operatorname{Poisson}(\lambda) の積率母関数であり、t=0t=0 の近傍で一致するので、ZnZ_n は Poisson⁡(λ)\operatorname{Poisson}(\lambda) に分布収束する。

Reference​