跳到主要内容

筑波大学 理工情報生命学術院 システム情報工学研究群 社会工学学位プログラム 2023年8月実施 微积分

Author

思齐塾, 祭音Myyura

Description

数学

問題IとIIの両方に答えよ。問題ごとに別々の解答用紙を使用せよ。 以下では、実数全体の集合を R\mathbb{R} とする。

I. 以下の問 [1] 及び [2] に答えよ。

[1] 平面 R2\mathbb{R}^2 上で定義された次の関数 f(x,y)f(x,y) について、以下の問 (1.1)-(1.4) に答えよ。

f(x,y)={xy3x2+y4((x,y)(0,0)のとき),0((x,y)=(0,0)のとき).f(x,y) = \begin{cases} \frac{xy^3}{x^2+y^4} & ((x,y) \neq (0,0) のとき), \\ 0 & ((x,y) = (0,0) のとき). \end{cases}

(1.1) 関数 f(x,y)f(x,y) は連続であることを示せ。

(1.2) 関数 f(x,y)f(x,y) は原点 (0,0)(0,0) で、両変数に関して偏微分可能であることを示せ。

(1.3) (x,y)(0,0)(x,y) \neq (0,0) のとき、関数 f(x,y)f(x,y) の偏導関数 fx(x,y),fy(x,y)f_x(x,y), f_y(x,y) を求めよ。

(1.4) 偏導関数 fx(x,y),fy(x,y)f_x(x,y), f_y(x,y) はともに原点 (0,0)(0,0) において連続でないことを示せ。

[2] ガンマ関数

Γ(s)=0xs1exdx,s>0\Gamma(s) = \int_0^{\infty} x^{s-1}e^{-x} dx, \quad s>0

について、以下の問 (2.1)-(2.3) に答えよ。

(2.1) 公式 0et2dt=π2\int_0^{\infty} e^{-t^2} dt = \frac{\sqrt{\pi}}{2} を利用して、 Γ(12)\Gamma(\frac{1}{2}) を求めよ。

(2.2) s>0s > 0 のとき、部分積分法を用いることにより、 Γ(s+1)=sΓ(s)\Gamma(s+1) = s\Gamma(s) を示せ。

(2.3) 各自然数 nn に対して、 Γ(n+12)\Gamma(n+\frac{1}{2}) を求めよ。

题目描述

原题开头要求问题 I、II 均作答并分别使用答题纸,但当前文件的 Description 与 Kai 只包含问题 I;问题 II 的题面缺失,无法由本文件唯一恢复。现有问题 I 包含以下两组题。

[1] 二元函数

在平面 R2\mathbb R^2 上定义

f(x,y)={xy3x2+y4,(x,y)(0,0),0,(x,y)=(0,0).f(x,y)= \begin{cases} \dfrac{xy^3}{x^2+y^4},&(x,y)\neq(0,0),\\[5pt] 0,&(x,y)=(0,0). \end{cases}
  1. 证明 f(x,y)f(x,y)R2\mathbb R^2 上连续。
  2. 证明 ff 在原点 (0,0)(0,0) 关于两个变量的偏导数都存在。
  3. (x,y)(0,0)(x,y)\neq(0,0) 时,求偏导函数 fx(x,y)f_x(x,y)fy(x,y)f_y(x,y)
  4. 证明 fxf_xfyf_y 在原点 (0,0)(0,0) 都不连续。

[2] Gamma 函数

s>0s>0,定义

Γ(s)=0xs1exdx.\Gamma(s)=\int_0^\infty x^{s-1}e^{-x}\,dx.
  1. 使用已知公式

    0et2dt=π2\int_0^\infty e^{-t^2}\,dt=\frac{\sqrt\pi}{2}

    Γ ⁣(12)\Gamma\!\left(\frac12\right)

  2. s>0s>0,使用分部积分证明递推公式

    Γ(s+1)=sΓ(s).\Gamma(s+1)=s\Gamma(s).
  3. 对每个自然数 nn,求

    Γ ⁣(n+12).\Gamma\!\left(n+\frac12\right).

Kai

解答:

[1] (1.1) 連続性: (x,y)(0,0)(x,y)\neq(0,0) では分母が正なので, ff は連続である. また, 2xy2x2+y42|x|y^2\leq x^2+y^4 より,

f(x,y)=yxy2x2+y4y2|f(x,y)| =|y|\frac{|x|y^2}{x^2+y^4} \leq\frac{|y|}{2}

が成り立つ. (x,y)(0,0)(x,y)\to(0,0) なら y/20|y|/2\to0 であるから, はさみうちの原理により f(x,y)0=f(0,0)f(x,y)\to0=f(0,0) となる. したがって, ff は原点でも連続であり,平面全体で連続である.

(1.2) 偏微分可能性:

fx(0,0)=limh0f(h,0)f(0,0)h=limh000h=0f_x(0,0) = \lim_{h\to 0} \frac{f(h,0) - f(0,0)}{h} = \lim_{h\to 0} \frac{0 - 0}{h} = 0
fy(0,0)=limk0f(0,k)f(0,0)k=limk000k=0f_y(0,0) = \lim_{k\to 0} \frac{f(0,k) - f(0,0)}{k} = \lim_{k\to 0} \frac{0 - 0}{k} = 0

したがって、 f(x,y)f(x,y) は原点 (0,0)(0,0) で偏微分可能です。

(1.3) 偏導関数の計算: (x,y)(0,0)(x,y) \neq (0,0) のとき

fx(x,y)=xxy3x2+y4=y3(x2+y4)xy3(2x)(x2+y4)2=y3(y4x2)(x2+y4)2f_x(x,y) = \frac{\partial}{\partial x} \frac{xy^3}{x^2+y^4} = \frac{y^3(x^2+y^4) - xy^3(2x)}{(x^2+y^4)^2} = \frac{y^3(y^4-x^2)}{(x^2+y^4)^2}
fy(x,y)=yxy3x2+y4=x(3y2(x2+y4)y3(4y3))(x2+y4)2=x(3x2y2y6)(x2+y4)2f_y(x,y) = \frac{\partial}{\partial y} \frac{xy^3}{x^2+y^4} = \frac{x(3y^2(x^2+y^4) - y^3(4y^3))}{(x^2+y^4)^2} = \frac{x(3x^2y^2 - y^6)}{(x^2+y^4)^2}

(1.4) 偏導関数の連続性:

fx(0,0)=0f_x(0,0) = 0
fx(x,y)=y3(y4x2)(x2+y4)2f_x(x,y) = \frac{y^3(y^4-x^2)}{(x^2+y^4)^2}

x=y2x = y^2 に沿って原点に近づくと

limy0fx(y2,y)=limy0y3(y4y4)(y4+y4)2=0\lim_{y\to 0} f_x(y^2, y) = \lim_{y\to 0} \frac{y^3(y^4 - y^4)}{(y^4 + y^4)^2} = 0

x=0x = 0 に沿って原点に近づくと

limy0fx(0,y)=limy0y7y8=limy01y\lim_{y\to 0} f_x(0,y) = \lim_{y\to 0} \frac{y^7}{y^8} = \lim_{y\to 0} \frac{1}{y}

これは存在しないので、 fxf_x は原点で連続ではありません。 また, x=y2x=y^2 に沿って原点に近づくと

fy(y2,y)=y2(3y6y6)(y4+y4)2=12f_y(y^2,y) =\frac{y^2(3y^6-y^6)}{(y^4+y^4)^2} =\frac12

となる一方, fy(0,0)=0f_y(0,0)=0 である.したがって, fyf_y も原点で連続ではない.

[2] (2.1) Γ(12)\Gamma(\frac{1}{2}) の計算:

Γ(12)=0x121exdx=0x12exdx\Gamma(\frac{1}{2}) = \int_0^{\infty} x^{\frac{1}{2}-1}e^{-x} dx = \int_0^{\infty} x^{-\frac{1}{2}}e^{-x} dx

x=t2x = t^2 を代入すると、 dx=2tdtdx = 2t dt

Γ(12)=0(t2)12et22tdt=20t1et2tdt=20et2dt=2π2=π\Gamma(\frac{1}{2}) = \int_0^{\infty} (t^2)^{-\frac{1}{2}}e^{-t^2} 2t dt = 2\int_0^{\infty} t^{-1}e^{-t^2} t dt = 2\int_0^{\infty} e^{-t^2} dt = 2 \cdot \frac{\sqrt{\pi}}{2} = \sqrt{\pi}

(2.2) Γ(s+1)=sΓ(s)\Gamma(s+1) = s\Gamma(s) の証明: 部分積分法を使う:

Γ(s+1)=0x(s+1)1exdx=0xsexdx\Gamma(s+1) = \int_0^{\infty} x^{(s+1)-1}e^{-x} dx = \int_0^{\infty} x^{s}e^{-x} dx

u=xs,dv=exdxu = x^s, dv = e^{-x} dx とすると, du=sxs1dx,v=exdu = sx^{s-1} dx, v = -e^{-x}

Γ(s+1)=[xsex]0+0exsxs1dx\Gamma(s+1) = [-x^se^{-x}]_0^{\infty} + \int_0^{\infty} e^{-x} sx^{s-1} dx
limxxsex=0\lim_{x\to \infty} x^s e^{-x} = 0
Γ(s+1)=0+s0xs1exdx=sΓ(s)\Gamma(s+1) = 0 + s\int_0^{\infty} x^{s-1} e^{-x} dx = s\Gamma(s)

(2.3) Γ(n+12)\Gamma(n+\frac{1}{2}) の計算: Γ(12)=π\Gamma(\frac{1}{2}) = \sqrt{\pi} を使用し, Γ(s+1)=sΓ(s)\Gamma(s+1) = s\Gamma(s) を繰り返し適用する。

Γ(n+12)=(n12)Γ(n12)=(n12)(n32)Γ(n32)==(n12)(n32)(12)Γ(12)\Gamma(n+\frac{1}{2}) = (n-\frac{1}{2})\Gamma(n-\frac{1}{2}) = (n-\frac{1}{2})(n-\frac{3}{2})\Gamma(n-\frac{3}{2}) = \cdots = (n-\frac{1}{2})(n-\frac{3}{2})\cdots(\frac{1}{2})\Gamma(\frac{1}{2})
Γ(n+12)=(2n1)(2n3)12nπ=(2n)!2n(2n)(2n2)2π=(2n)!2n2nn!π=(2n)!4nn!π\Gamma(n+\frac{1}{2}) = \frac{(2n-1)(2n-3)\cdots 1}{2^n} \sqrt{\pi} = \frac{(2n)!}{2^n(2n)(2n-2)\cdots 2} \sqrt{\pi} = \frac{(2n)!}{2^n 2^n n!} \sqrt{\pi} = \frac{(2n)!}{4^n n!} \sqrt{\pi}