跳到主要内容

九州大学 工学府 土木工学専攻 2021年8月実施 数学 (ALL)

Author​

Miyake, 祭音Myyura

Description​

試験問題は【問題1】から【問題6】の計6問です。すべての問題を解答してください。

【問題1】(15点)​

関数 y=xlog⁡xy=x\log x の概形を、解答用紙に与えられた範囲内で示せ。この際、以下の2点に注意して作図すること。

  • 図中の範囲内に極値あるいは変曲点があるかを調べた上で、存在する場合にはその座標の値を図中と解答欄に示せ。ない場合には「なし」と示すこと。
  • xx 軸との交点は必ず求めること。また、この点での接線を、破線の直線で示し、この傾きの値を図中と解答欄に示せ。

問題1の作図範囲

【問題2】(15点)​

双曲面 z=f(x,y)=xyz=f(x,y)=xy が円柱面 x2+y2=4x^2+y^2=4 により切り取られる部分の表面積 SS を、以下の手順により求めよ。

(1) 求めるべき図形は原点に対して対称であることを使えば、x≥0,y≥0x\ge0,y\ge0 の部分の表面積の4倍として評価できる。この部分領域を D={(x,y)≤4, x≥0, y≥0}D=\{(x,y)\le4,\ x\ge0,\ y\ge0\} として与えることにし、求めるべき表面積 SS を与える式を2重積分の形として x,yx,y を使って表せ(答えは S=○○S=\text{○○} の形で解答すること)。このとき、2重積分の領域は部分領域 DD を使って与えること。

(2) (1) で与えた式を領域積分し、表面積 SS の値を求めよ。((1) の結果を使わずに、そのまま求めても良い。)

【問題3】(24点)​

1階の常微分方程式①を考える。以下の (1)~(3) の問いに答えよ。

X(x,y) dx+Y(x,y) dy=0.(①)X(x,y)\,dx+Y(x,y)\,dy=0.\tag{①}

(1) 式①が ∂X/∂y=∂Y/∂x\partial X/\partial y=\partial Y/\partial x を満たすとき、完全微分方程式(または完全微分形)という。完全微分方程式②の一般解を求めよ。

(y−x3) dx+(x−sin⁡y) dy=0.(②)(y-x^3)\,dx+(x-\sin y)\,dy=0.\tag{②}

(2) 式①が完全微分方程式でない場合、式③のように両辺に x,yx,y の関数 P(x,y)≠0P(x,y)\ne0 を乗じると完全微分方程式になることがある。このとき、P(x,y)P(x,y) を積分因子(または積分因数)という。もし、積分因子が式④に示すように yy のみの関数 Q(y)Q(y) であるとき、Q(y)Q(y) は、X,Y,x,yX,Y,x,y を用いてどのように表されるか、示せ。

なお、例えば、X(x,y)X(x,y) の yy による偏微分は、∂X/∂y\partial X/\partial y のように表してよい。

P(x,y)X(x,y) dx+P(x,y)Y(x,y) dy=0.(③)P(x,y)X(x,y)\,dx+P(x,y)Y(x,y)\,dy=0.\tag{③}
Q(y)X(x,y) dx+Q(y)Y(x,y) dy=0.(④)Q(y)X(x,y)\,dx+Q(y)Y(x,y)\,dy=0.\tag{④}

(3) 微分方程式⑤は、yy のみの関数を積分因子として持つことがあらかじめ分かっている。(2) の結果を用いて、式⑤の一般解を求めよ。

(xy2−y3) dx+(1−xy2) dy=0.(⑤)(xy^2-y^3)\,dx+(1-xy^2)\,dy=0.\tag{⑤}

【問題4】(30点)​

以下の (1)~(3) の問いに答えよ。

(1) 次の行列 AA について、nn 乗 AnA^n を求めよ。nn は2以上の整数とする。

A=(0100001000010000).A=\begin{pmatrix}0&1&0&0\\0&0&1&0\\0&0&0&1\\0&0&0&0\end{pmatrix}.

(2) 次の行列 BB について、nn 乗 BnB^n を求めよ。nn は2以上の整数、b≠0b\ne0 とする。

B=(b1000b1000b1000b).B=\begin{pmatrix}b&1&0&0\\0&b&1&0\\0&0&b&1\\0&0&0&b\end{pmatrix}.

(3) 次の行列 CC を対角化せよ。

C=(13−305−12).C=\begin{pmatrix}13&-30\\5&-12\end{pmatrix}.

【問題5】(8点)​

確率密度が以下の f(x)f(x) で与えられる連続分布について、平均と分散を求めよ。

f(x)={2(1−x)(0≤x≤1),0(x<0, x>1).f(x)=\begin{cases}2(1-x)&(0\le x\le1),\\0&(x<0,\ x>1).\end{cases}

【問題6】(8点)​

A工場で製造された鉄筋から無作為に9本を選び、引張強度を測定したところ、標本平均 Xˉ\bar X は 500 N/mm2500\,\mathrm{N/mm^2} であった。A工場で製造される鉄筋の引張強度は正規分布に従い、引張強度の母標準偏差 σ\sigma が 75 N/mm275\,\mathrm{N/mm^2} であることが分かっているとき、A工場で製造される鉄筋の引張強度の母平均 μ\mu の95%信頼区間を求めよ。解答にあたり、下の附表を参照してよい。

標準正規分布の両側の陰影部分

附表:陰影部の面積の和が α\alpha となる zz の値

α\alphazz
0.012.576
0.022.326
0.051.960
0.101.645
0.201.282

题目描述​

试卷共6题,全部作答。

问题1(15分) 在答题纸给出的范围内画出 y=xlog⁡xy=x\log x 的概形,注意以下两点。

  • 判断图中范围内是否有极值点或拐点,若存在,在图中及答题栏标明坐标;不存在时写“无”。
  • 必须求出与 xx 轴的交点,用虚线画出该点处的切线,并在图中及答题栏标明斜率。

作图范围见原图,两轴均为 −0.5-0.5 到 2.52.5。

问题2(15分) 按以下步骤,求双曲面 z=f(x,y)=xyz=f(x,y)=xy 被圆柱面 x2+y2=4x^2+y^2=4 截取部分的表面积 SS。

(1) 利用图形关于原点的对称性,可将面积算作 x≥0,y≥0x\ge0,y\ge0 部分的4倍。原文将这部分区域记为 D={(x,y)≤4, x≥0, y≥0}D=\{(x,y)\le4,\ x\ge0,\ y\ge0\}。用 x,yx,y 将表面积 SS 表示为二重积分,答案写成 S=○○S=\text{○○} 的形式,积分区域使用 DD 表示。

(2) 计算 (1) 中的区域积分,求 SS 的值。也可以不使用 (1) 的结果直接求解。

问题3(24分) 考虑一阶常微分方程①,回答 (1)~(3)。

X(x,y) dx+Y(x,y) dy=0.X(x,y)\,dx+Y(x,y)\,dy=0.

(1) 当式①满足 ∂X/∂y=∂Y/∂x\partial X/\partial y=\partial Y/\partial x 时,称为恰当微分方程(全微分形式)。求式②的通解:

(y−x3) dx+(x−sin⁡y) dy=0.(y-x^3)\,dx+(x-\sin y)\,dy=0.

(2) 若式①不是恰当方程,有时像式③那样在两边乘以 x,yx,y 的非零函数 P(x,y)P(x,y) 后可变为恰当方程,称 P(x,y)P(x,y) 为积分因子。若积分因子像式④那样只依赖 yy,记为 Q(y)Q(y),请用 X,Y,x,yX,Y,x,y 表示 Q(y)Q(y)。

例如,X(x,y)X(x,y) 对 yy 的偏导数可以写作 ∂X/∂y\partial X/\partial y。

式③、④分别为

P(x,y)X(x,y) dx+P(x,y)Y(x,y) dy=0,P(x,y)X(x,y)\,dx+P(x,y)Y(x,y)\,dy=0,
Q(y)X(x,y) dx+Q(y)Y(x,y) dy=0.Q(y)X(x,y)\,dx+Q(y)Y(x,y)\,dy=0.

(3) 已知微分方程⑤具有只依赖 yy 的积分因子,利用 (2) 的结果求其通解:

(xy2−y3) dx+(1−xy2) dy=0.(xy^2-y^3)\,dx+(1-xy^2)\,dy=0.

问题4(30分) 回答 (1)~(3)。

(1) 对以下矩阵 AA 求 AnA^n,其中 nn 为不小于2的整数。

A=(0100001000010000).A=\begin{pmatrix}0&1&0&0\\0&0&1&0\\0&0&0&1\\0&0&0&0\end{pmatrix}.

(2) 对以下矩阵 BB 求 BnB^n,其中 nn 为不小于2的整数,b≠0b\ne0。

B=(b1000b1000b1000b).B=\begin{pmatrix}b&1&0&0\\0&b&1&0\\0&0&b&1\\0&0&0&b\end{pmatrix}.

(3) 对角化矩阵

C=(13−305−12).C=\begin{pmatrix}13&-30\\5&-12\end{pmatrix}.

问题5(8分) 求概率密度为下式的连续分布的均值和方差。

f(x)={2(1−x)(0≤x≤1),0(x<0, x>1).f(x)=\begin{cases}2(1-x)&(0\le x\le1),\\0&(x<0,\ x>1).\end{cases}

问题6(8分) 从A工厂生产的钢筋中随机抽取9根,测得抗拉强度样本均值 Xˉ=500 N/mm2\bar X=500\,\mathrm{N/mm^2}。已知该厂钢筋抗拉强度服从正态分布,总体标准差 σ=75 N/mm2\sigma=75\,\mathrm{N/mm^2},求总体均值 μ\mu 的95%置信区间。可以参考上方附图及附表;表中 zz 对应两侧阴影面积之和为 α\alpha。

Kai​

【問題 1】​

定義域は x>0x>0 であり、

y′=log⁡x+1,y′′=1x>0.y'=\log x+1,\qquad y''=\frac1x>0.

したがって、0<x<e−10<x<e^{-1} で減少、x>e−1x>e^{-1} で増加し、極小点は (e−1,−e−1)(e^{-1},-e^{-1})、極大点と変曲点はない。xx 軸との交点は (1,0)(1,0)、そこでの接線の傾きは 11 で、接線は y=x−1y=x-1 となる。また lim⁡x→0+xlog⁡x=0\lim_{x\to0+}x\log x=0 だが、原点はグラフに含まれない。

y=x log x と交点での接線

【問題 2】​

原文の領域の表記 D={(x,y)≤4, x≥0, y≥0}D=\{(x,y)\le4,\ x\ge0,\ y\ge0\} は、円柱面の条件から D={(x,y):x2+y2≤4, x≥0, y≥0}D=\{(x,y):x^2+y^2\le4,\ x\ge0,\ y\ge0\} と解釈する。

(1)​

曲面 z=f(x,y)z=f(x,y) 上の点 P=(x,y,f(x,y))P=(x,y,f(x,y)), Q=(x+Δx,y,f(x+Δx,y))Q=(x+\Delta x,y,f(x+\Delta x,y)), R=(x,y+Δy,f(x,y+Δy))R=(x,y+\Delta y,f(x,y+\Delta y)) を考えると、 Δx,Δy\Delta x, \Delta y の1次までで

PQ→≃(Δx,0,∂f∂xΔx)=(1,0,∂f∂x)ΔxPR→≃(0,Δy,∂f∂yΔy)=(0,1,∂f∂y)Δy\begin{aligned} \overrightarrow{PQ} &\simeq \left( \Delta x, 0, \frac{\partial f}{\partial x} \Delta x \right) = \left( 1, 0, \frac{\partial f}{\partial x} \right) \Delta x \\ \overrightarrow{PR} &\simeq \left( 0, \Delta y, \frac{\partial f}{\partial y} \Delta y \right) = \left( 0, 1, \frac{\partial f}{\partial y} \right) \Delta y \end{aligned}

であり、これらのベクトル積とその大きさは

PQ→×PR→=(−∂f∂x,−∂f∂y,1)ΔxΔy∣PQ→×PR→∣=∣ΔxΔy∣(∂f∂x)2+(∂f∂y)2+1=∣ΔxΔy∣x2+y2+1\begin{aligned} \overrightarrow{PQ} \times \overrightarrow{PR} &= \left( - \frac{\partial f}{\partial x}, - \frac{\partial f}{\partial y}, 1 \right) \Delta x \Delta y \\ \left| \overrightarrow{PQ} \times \overrightarrow{PR} \right| &= \left| \Delta x \Delta y \right| \sqrt{ \left( \frac{\partial f}{\partial x} \right)^2 + \left( \frac{\partial f}{\partial y} \right)^2 + 1 } \\ &= \left| \Delta x \Delta y \right| \sqrt{ x^2 + y^2 + 1 } \end{aligned}

である。 よって、

S=4∬Ddxdyx2+y2+1\begin{aligned} S = 4 \iint_D dx dy \sqrt{x^2+y^2+1} \end{aligned}

がわかる。

(2)​

x=rcos⁡θ, y=rsin⁡θx = r \cos \theta, \ y = r \sin \theta によって2次元極座標 (r,θ)(r, \theta) を導入すると、 dxdy=rdrdθdxdy = r dr d \theta であり、次のように計算できる:

S=4∫02drrr2+1∫0π/2dθ=2π3[(r2+1)3/2]02=2(55−1)π3\begin{aligned} S &= 4 \int_0^2 dr r \sqrt{r^2 + 1} \int_0^{\pi/2} d \theta \\ &= \frac{2 \pi}{3} \left[ \left( r^2 + 1 \right)^{3/2} \right]_0^2 \\ &= \frac{2 \left( 5 \sqrt{5} - 1 \right) \pi}{3} \end{aligned}

【問題 3】​

(1)​

∂f∂x=y−x3∂f∂y=x−sin⁡y\begin{align} \frac{\partial f}{\partial x} &= y - x^3 \tag{a} \\ \frac{\partial f}{\partial y} &= x - \sin y \tag{b} \end{align}

を満たす関数 f(x,y)f(x,y) を求める。 式 (aa) より、

f(x,y)=xy−14x4+g(y)\begin{align} f(x,y) = xy - \frac{1}{4} x^4 + g(y) \tag{c} \end{align}

なる関数 g(y)g(y) が存在することがわかる。 (cc) を (bb) に代入して整理すると、

dgdy=−sin⁡y∴  g(y)=cos⁡y+C         ( C は任意定数 )∴  f(x,y)=xy−14x4+cos⁡y+C         ( C は任意定数 )\begin{aligned} \frac{dg}{dy} &= - \sin y \\ \therefore \ \ g(y) &= \cos y + C \ \ \ \ \ \ \ \ \text{ ( $C$ は任意定数 )} \\ \therefore \ \ f(x,y) &= xy - \frac{1}{4} x^4 + \cos y + C \ \ \ \ \ \ \ \ \text{ ( $C$ は任意定数 )} \end{aligned}

を得る。 よって、求める一般解は

xy−14x4+cos⁡y+C=0         ( C は任意定数 )\begin{aligned} xy - \frac{1}{4} x^4 + \cos y + C = 0 \ \ \ \ \ \ \ \ \text{ ( $C$ は任意定数 )} \end{aligned}

である。

(2)​

式④が完全微分方程式であるための条件は、

∂∂yQ(y)X(x,y)=∂∂xQ(y)Y(x,y)dQdyX+Q∂X∂y=Q∂Y∂x1QdQdy=1X(∂Y∂x−∂X∂y)\begin{aligned} \frac{\partial}{\partial y} Q(y) X(x,y) &= \frac{\partial}{\partial x} Q(y) Y(x,y) \\ \frac{dQ}{dy} X + Q \frac{\partial X}{\partial y} &= Q \frac{\partial Y}{\partial x} \\ \frac{1}{Q} \frac{dQ}{dy} &= \frac{1}{X} \left( \frac{\partial Y}{\partial x} - \frac{\partial X}{\partial y} \right) \end{aligned}

なので、

Q(y)=exp⁡(∫dy1X(x,y)(∂Y∂x−∂X∂y))\begin{aligned} Q(y) &= \exp \left( \int dy \frac{1}{X(x,y)} \left( \frac{\partial Y}{\partial x} - \frac{\partial X}{\partial y} \right) \right) \end{aligned}

を得る。ただし右辺の被積分関数が yy のみの関数となることが必要である。

(3)​

(2) の X(x,y),Y(x,y)X(x,y),Y(x,y) は今の場合、

X(x,y)=xy2−y3,  Y(x,y)=1−xy2\begin{aligned} X(x,y) = xy^2 - y^3 , \ \ Y(x,y) = 1 - xy^2 \end{aligned}

であるから、

1X(x,y)(∂Y∂x−∂X∂y)=−y2−(2xy−3y2)xy2−y3=−2y\begin{aligned} \frac{1}{X(x,y)} \left( \frac{\partial Y}{\partial x} - \frac{\partial X}{\partial y} \right) &= \frac{-y^2 - (2xy-3y^2) }{xy^2 - y^3} \\ &= - \frac{2}{y} \end{aligned}

であり、積分因子として

Q(y)=1y2\begin{aligned} Q(y) = \frac{1}{y^2} \end{aligned}

を考えればよい。ただし y≠0y\ne0 とする。 このとき与えられた微分方程式は、

(x−y)dx+(1y2−x)dy=0\begin{aligned} (x-y)dx + \left( \frac{1}{y^2} - x \right) dy = 0 \end{aligned}

となり、 (1) と同じ方法を使って、一般解

12x2−xy−1y+C=0         ( C は任意定数 )\begin{aligned} \frac{1}{2} x^2 - xy - \frac{1}{y} + C = 0 \ \ \ \ \ \ \ \ \text{ ( $C$ は任意定数 )} \end{aligned}

を得る。なお、積分因子を掛ける前の方程式には特異解 y=0y=0 もある。

【問題 4】​

(1)​

A2=(0010000100000000)A3=(0001000000000000)An=(0000000000000000)        (n=4,5,6,⋯ ) \begin{aligned} A^2 &= \begin{pmatrix} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix} \\ A^3 &= \begin{pmatrix} 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix} \\ A^n &= \begin{pmatrix} 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix} \ \ \ \ \ \ \ \ (n = 4, 5, 6, \cdots) \end{aligned}

(2)​

4次の単位行列を EE とすると、

B=bE+A\begin{aligned} B = bE+A \end{aligned}

であり、まず、

B2=b2E+2bA+A2=(b22b100b22b100b22b000b2)\begin{aligned} B^2 &= b^2 E + 2bA + A^2 \\ &= \begin{pmatrix} b^2 & 2b & 1 & 0 \\ 0 & b^2 & 2b & 1 \\ 0 & 0 & b^2 & 2b \\ 0 & 0 & 0 & b^2 \end{pmatrix} \end{aligned}

がわかる。 また、

B3=b3E+3b2A+3bA2+A3=(b33b23b10b33b23b00b33b2000b3)B4=b4E+4b3A+6b2A2+4bA3+A4=b4E+4b3A+6b2A2+4bA3=(b44b36b24b0b44b36b200b44b3000b4)\begin{aligned} B^3 &= b^3 E + 3b^2A + 3bA^2 + A^3 \\ &= \begin{pmatrix} b^3 & 3b^2 & 3b & 1 \\ 0 & b^3 & 3b^2 & 3b \\ 0 & 0 & b^3 & 3b^2 \\ 0 & 0 & 0 & b^3 \end{pmatrix} \\ B^4 &= b^4 E + 4b^3A + 6b^2A^2 + 4bA^3 + A^4 \\ &= b^4 E + 4b^3A + 6b^2A^2 + 4bA^3 \\ &= \begin{pmatrix} b^4 & 4b^3 & 6b^2 & 4b \\ 0 & b^4 & 4b^3 & 6b^2 \\ 0 & 0 & b^4 & 4b^3 \\ 0 & 0 & 0 & b^4 \end{pmatrix} \end{aligned}

から、 n=3,4,5,⋯n = 3, 4, 5, \cdots について

Bn=bnE+nbn−1A+n(n−1)2bn−2A2+n(n−1)(n−2)6bn−3A3=(bnnbn−1n(n−1)2bn−2n(n−1)(n−2)6bn−30bnnbn−1n(n−1)2bn−200bnnbn−1000bn)\begin{gathered} B^n=b^nE+nb^{n-1}A+\frac{n(n-1)}2b^{n-2}A^2 \\ {}+\frac{n(n-1)(n-2)}6b^{n-3}A^3 \\ =\begin{pmatrix} b^n&nb^{n-1}&\frac{n(n-1)}2b^{n-2}&\frac{n(n-1)(n-2)}6b^{n-3}\\ 0&b^n&nb^{n-1}&\frac{n(n-1)}2b^{n-2}\\ 0&0&b^n&nb^{n-1}\\0&0&0&b^n \end{pmatrix} \end{gathered}

がわかる。

(3)​

det⁡(λI−C)=(λ+2)(λ−3)\det(\lambda I-C)=(\lambda+2)(\lambda-3) である。固有値 −2,3-2,3 に対する固有ベクトルとして、それぞれ (2,1)T,(3,1)T(2,1)^{\mathsf T},(3,1)^{\mathsf T} を取れる。したがって、

P=(2311),P−1=(−131−2),P−1CP=(−2003).\begin{gathered} P=\begin{pmatrix}2&3\\1&1\end{pmatrix},\\ P^{-1}=\begin{pmatrix}-1&3\\1&-2\end{pmatrix},\\ P^{-1}CP=\begin{pmatrix}-2&0\\0&3\end{pmatrix}. \end{gathered}

【問題 5】​

平均は

μ=∫−∞∞xf(x)dx=2∫01(x−x2)dx=2[x22−x33]01=13\begin{aligned} \mu &= \int_{- \infty}^\infty x f(x) dx \\ &= 2 \int_0^1 \left( x - x^2 \right) dx \\ &= 2 \left[ \frac{x^2}{2} - \frac{x^3}{3} \right]_0^1 \\ &= \frac{1}{3} \end{aligned}

であり、分散は

σ2=∫−∞∞(x−μ)2f(x)dx=∫−∞∞x2f(x)dx−μ2=2∫01(x2−x3)dx−μ2=2[x33−x44]01−19=118\begin{aligned} \sigma^2 &= \int_{- \infty}^\infty (x - \mu)^2 f(x) dx \\ &= \int_{- \infty}^\infty x^2 f(x) dx - \mu^2 \\ &= 2 \int_0^1 \left( x^2 - x^3 \right) dx - \mu^2 \\ &= 2 \left[ \frac{x^3}{3} - \frac{x^4}{4} \right]_0^1 - \frac{1}{9} \\ &= \frac{1}{18} \end{aligned}

である。

【問題 6】​

1.96⋅759=49 \begin{aligned} 1.96 \cdot \frac{75}{\sqrt{9}} = 49 \end{aligned}

なので、求める95%信頼区間の下端は 500−49=451 N/mm2500-49=451 \ \mathrm{N/mm^2} であり、 上端は 500+49=549 N/mm2500+49=549 \ \mathrm{N/mm^2} である。

Reference​