跳到主要内容

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

Author​

Miyake, 祭音Myyura

Description​

題意の要約​

出典:九州大学 2020年度 経済数学。 問1・問2の両方に解答し、それぞれで (1)、(2) の一方を選ぶ。

問1​

(1)​

t∈Rt\in\mathbb R とし、

x1=(201),x2=(909),x3=(1t8),03=(000),X=(x1,x2,x3).\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},\quad X=(\boldsymbol x_1,\boldsymbol x_2,\boldsymbol x_3).
  • (a) 3列が線形独立となる tt の条件を求める。
  • (b) 任意の3次元ベクトル c\boldsymbol c について、X⊤Xc=03X^\top X\boldsymbol c=\boldsymbol0_3 から Xc=03X\boldsymbol c=\boldsymbol0_3 が従うことを、X⊤XX^\top X の正則性を仮定せず示す。⊤\top は転置を表す。
  • (c) X⊤XX^\top X が正則となる tt の条件を求める。

(2)​

a>0a>0 とする。

  • (a) xzxz 平面内の曲線 z=(x−a)2z=(x-a)^2(x≥0x\ge0)を zz 軸のまわりに回転した曲面の式が z=(x2+y2−a)2z=(\sqrt{x^2+y^2}-a)^2 となることを示す。
  • (b) f(x,y)=(x2+y2−a)2f(x,y)=(\sqrt{x^2+y^2}-a)^2 の極値を求める。

問2​

(1)​

三角形 ABCABC の辺 BC,CA,ABBC,CA,AB の長さを a,b,ca,b,c、面積を SS とする。 内部の点 DD から直線 BC,CA,ABBC,CA,AB までの距離を順に x,y,zx,y,z とする。

  • (a) (x,y)(x,y) の取り得る範囲を求める。
  • (b) x2+y2+z2x^2+y^2+z^2 の極小値を求める。
  • (c) (b) の値が最小値でもあることを示す。
  • (d) xyzxyz の最大値を求める。

(2)​

n∈Nn\in\mathbb N、0<p<10<p<1、t∈Rt\in\mathbb R、λ>0\lambda>0 とする。

  • (a) X∼Bin⁡(n,p)X\sim\operatorname{Bin}(n,p) の確率関数は pX(x)=(nx)px(1−p)n−xp_X(x)=\binom nx p^x(1-p)^{n-x}(x=0,…,nx=0,\ldots,n)。E[etX]E[e^{tX}] を求める。
  • (b) Y∼Poisson⁡(λ)Y\sim\operatorname{Poisson}(\lambda) の確率関数は pY(y)=e−λλy/y!p_Y(y)=e^{-\lambda}\lambda^y/y!(y=0,1,…y=0,1,\ldots)。E[etY]E[e^{tY}] を求める。
  • (c) n>λn>\lambda とし、Zn∼Bin⁡(n,λ/n)Z_n\sim\operatorname{Bin}(n,\lambda/n) について lim⁡n→∞E[etZn]\lim_{n\to\infty}E[e^{tZ_n}] を求める。

题目描述​

第 1、2 题都要作答;每题在 (1)、(2) 中任选一项。

第 1 题 (1):令 t∈Rt\in\mathbb R,并以

x1=(2,0,1)⊤,x2=(9,0,9)⊤,x3=(1,t,8)⊤\boldsymbol x_1=(2,0,1)^\top,\quad \boldsymbol x_2=(9,0,9)^\top,\quad \boldsymbol x_3=(1,t,8)^\top

为列构成矩阵 XX。

  1. 求三个向量线性无关时 tt 的条件。
  2. 不假设 X⊤XX^\top X 可逆,证明对任意三维向量 c\boldsymbol c,若 X⊤Xc=0X^\top X\boldsymbol c=\boldsymbol0,则 Xc=0X\boldsymbol c=\boldsymbol0。
  3. 求 Gram 矩阵 X⊤XX^\top X 可逆时 tt 的条件。

第 1 题 (2):设 a>0a>0。

  1. 将 xzxz 平面上的曲线 z=(x−a)2z=(x-a)^2(x≥0x\ge0)绕 zz 轴旋转,证明所得曲面为 z=(x2+y2−a)2z=(\sqrt{x^2+y^2}-a)^2。
  2. 求函数 f(x,y)=(x2+y2−a)2f(x,y)=(\sqrt{x^2+y^2}-a)^2 的极值。

第 2 题 (1):三角形 ABCABC 的三边 BC,CA,ABBC,CA,AB 长度依次为 a,b,ca,b,c,面积为 SS。内部点 DD 到这三条边所在直线的距离依次为 x,y,zx,y,z。

  1. 求 (x,y)(x,y) 的取值范围。
  2. 求 x2+y2+z2x^2+y^2+z^2 的极小值。
  3. 证明该极小值也是最小值。
  4. 求 xyzxyz 的最大值。

第 2 题 (2):设 n∈Nn\in\mathbb N、0<p<10<p<1、λ>0\lambda>0、t∈Rt\in\mathbb R。

  1. 对 X∼Bin⁡(n,p)X\sim\operatorname{Bin}(n,p),求矩母函数 E(etX)E(e^{tX})。
  2. 对 Y∼Poisson⁡(λ)Y\sim\operatorname{Poisson}(\lambda),求矩母函数 E(etY)E(e^{tY})。
  3. 当 n>λn>\lambda、Zn∼Bin⁡(n,λ/n)Z_n\sim\operatorname{Bin}(n,\lambda/n) 时,求 lim⁡n→∞E(etZn)\lim_{n\to\infty}E(e^{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π.(x,y,z)=(\rho\cos\theta,\rho\sin\theta,(\rho-a)^2),\qquad 0\le\theta<2\pi.

ρ=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) に分布収束する。