跳到主要内容

京都大学 工学研究科 機械工学群 2022年8月実施 数学 問題1

Author​

祭音Myyura (Based on Miyake's answer refined with GPT 6 Astra)

Description​

  1. 実関数の内積を以下の式で定義する。
(f1(x),f2(x))=∫−ππf1(x)f2(x) dx.(1.1)(f_1(x),f_2(x))=\int_{-\pi}^{\pi}f_1(x)f_2(x)\,dx.\tag{1.1}

関数 f(x),g(x)f(x),g(x) を

f(x)=a012π+a11πcos⁡x+a21πsin⁡x,(1.2)f(x)=a_0\frac1{\sqrt{2\pi}}+a_1\frac1{\sqrt\pi}\cos x+a_2\frac1{\sqrt\pi}\sin x,\tag{1.2}
g(x)=b012π+b11πcos⁡x+b21πsin⁡x(1.3)g(x)=b_0\frac1{\sqrt{2\pi}}+b_1\frac1{\sqrt\pi}\cos x+b_2\frac1{\sqrt\pi}\sin x\tag{1.3}

とする。ここで、a0,a1,a2,b0,b1,b2a_0,a_1,a_2,b_0,b_1,b_2 は実定数である。以下の問いに答えよ。

1-1​

(f(x),g(x))(f(x),g(x)) を計算せよ。

1-2​

演算子 O^1=2+ddx\hat O_1=2+\dfrac{d}{dx} を考えると

(f(x),O^1g(x))=(a0a1a2)A(b0b1b2)(1.4)(f(x),\hat O_1g(x))= \begin{pmatrix}a_0&a_1&a_2\end{pmatrix}A \begin{pmatrix}b_0\\b_1\\b_2\end{pmatrix}\tag{1.4}

をみたす3行3列の実行列 AA が存在する。AA を具体的な値で示せ。

1-3​

演算子 O^2=1+2ddx+2d2dx2\hat O_2=1+2\dfrac{d}{dx}+2\dfrac{d^2}{dx^2} について、(f(x),(O^2)3f(x))(f(x),(\hat O_2)^3f(x)) を計算せよ。

1-4​

次の微分方程式について考える。

(2d4dx4+d2dx2)f(x)=λf(x).(1.5)\left(2\frac{d^4}{dx^4}+\frac{d^2}{dx^2}\right)f(x)=\lambda f(x).\tag{1.5}

境界条件は f(π)=−1,f′(−π)=0f(\pi)=-1,f'(-\pi)=0 とし、λ\lambda は実定数とする。ここで f′(x)f'(x) は f(x)f(x) の導関数である。この微分方程式が成立する全ての λ\lambda の値とそれぞれの λ\lambda に対応した a0,a1,a2a_0,a_1,a_2 の値を答えよ。

1-5​

実定数 α0,α1,…,αn,β1,β2,…,βn\alpha_0,\alpha_1,\ldots,\alpha_n,\beta_1,\beta_2,\ldots,\beta_n を用いて関数 gn(x)g_n(x) を

gn(x)=12α0+∑k=1n(αkcos⁡kx+βksin⁡kx)(1.6)g_n(x)=\frac12\alpha_0+\sum_{k=1}^n(\alpha_k\cos kx+\beta_k\sin kx)\tag{1.6}

と定義する。関数 gn(x)g_n(x) と 2π2\pi を周期とする次の関数

h(x)=∣x∣(−π≤x<π),h(x+2π)=h(x)(1.7)h(x)=|x|\quad(-\pi\le x<\pi),\qquad h(x+2\pi)=h(x)\tag{1.7}

との自乗平均誤差

12π∫−ππ(h(x)−gn(x))2 dx(1.8)\frac1{2\pi}\int_{-\pi}^{\pi}(h(x)-g_n(x))^2\,dx\tag{1.8}

が最小となるように α0,αk,βk (k=1,…,n)\alpha_0,\alpha_k,\beta_k\ (k=1,\ldots,n) を決定する。このとき α0,αk,βk (k=1,…,n)\alpha_0,\alpha_k,\beta_k\ (k=1,\ldots,n) はフーリエ係数となっている。α0,αk,βk (k=1,…,n)\alpha_0,\alpha_k,\beta_k\ (k=1,\ldots,n) を計算するとともに、lim⁡n→∞(gn(x),gn(x))\displaystyle\lim_{n\to\infty}(g_n(x),g_n(x)) の値がいくらになるかを計算し、その過程を簡潔に説明せよ。

题目描述​

  1. 实函数的内积定义为
(f1(x),f2(x))=∫−ππf1(x)f2(x) dx.(1.1)(f_1(x),f_2(x))=\int_{-\pi}^{\pi}f_1(x)f_2(x)\,dx.\qquad(1.1)

设函数 f(x),g(x)f(x),g(x) 为

f(x)=a012π+a11πcos⁡x+a21πsin⁡x,(1.2)f(x)=a_0\frac1{\sqrt{2\pi}}+a_1\frac1{\sqrt\pi}\cos x+a_2\frac1{\sqrt\pi}\sin x,\qquad(1.2)
g(x)=b012π+b11πcos⁡x+b21πsin⁡x,(1.3)g(x)=b_0\frac1{\sqrt{2\pi}}+b_1\frac1{\sqrt\pi}\cos x+b_2\frac1{\sqrt\pi}\sin x,\qquad(1.3)

其中 a0,a1,a2,b0,b1,b2a_0,a_1,a_2,b_0,b_1,b_2 均为实常数。回答以下问题。

1-1

计算 (f(x),g(x))(f(x),g(x))。

1-2

考虑算子 O^1=2+ddx\hat O_1=2+\dfrac{d}{dx},存在一个 3×33\times3 实矩阵 AA 满足

(f(x),O^1g(x))=(a0a1a2)A(b0b1b2).(1.4)(f(x),\hat O_1g(x))= \begin{pmatrix}a_0&a_1&a_2\end{pmatrix}A \begin{pmatrix}b_0\\b_1\\b_2\end{pmatrix}.\qquad(1.4)

求出 AA 的具体数值。

1-3

对于算子 O^2=1+2ddx+2d2dx2\hat O_2=1+2\dfrac{d}{dx}+2\dfrac{d^2}{dx^2},计算 (f(x),(O^2)3f(x))(f(x),(\hat O_2)^3f(x))。

1-4

考虑微分方程

(2d4dx4+d2dx2)f(x)=λf(x).(1.5)\left(2\frac{d^4}{dx^4}+\frac{d^2}{dx^2}\right)f(x)=\lambda f(x).\qquad(1.5)

边界条件为 f(π)=−1,f′(−π)=0f(\pi)=-1,f'(-\pi)=0,λ\lambda 为实常数,f′(x)f'(x) 为 f(x)f(x) 的导函数。求使该微分方程成立的所有 λ\lambda 值,以及每个 λ\lambda 对应的 a0,a1,a2a_0,a_1,a_2。

1-5

用实常数 α0,α1,…,αn,β1,β2,…,βn\alpha_0,\alpha_1,\ldots,\alpha_n,\beta_1,\beta_2,\ldots,\beta_n 定义函数

gn(x)=12α0+∑k=1n(αkcos⁡kx+βksin⁡kx).(1.6)g_n(x)=\frac12\alpha_0+\sum_{k=1}^n(\alpha_k\cos kx+\beta_k\sin kx).\qquad(1.6)

将函数 gn(x)g_n(x) 与以下周期为 2π2\pi 的函数进行比较:

h(x)=∣x∣(−π≤x<π),h(x+2π)=h(x).(1.7)h(x)=|x|\quad(-\pi\le x<\pi),\qquad h(x+2\pi)=h(x).\qquad(1.7)

选择 α0,αk,βk (k=1,…,n)\alpha_0,\alpha_k,\beta_k\ (k=1,\ldots,n),使均方误差

12π∫−ππ(h(x)−gn(x))2 dx(1.8)\frac1{2\pi}\int_{-\pi}^{\pi}(h(x)-g_n(x))^2\,dx\qquad(1.8)

最小。这时 α0,αk,βk (k=1,…,n)\alpha_0,\alpha_k,\beta_k\ (k=1,\ldots,n) 就是傅里叶系数。求这些系数,同时计算 lim⁡n→∞(gn(x),gn(x))\displaystyle\lim_{n\to\infty}(g_n(x),g_n(x)) 的值,并简要说明计算过程。

Kai​

1-1​

e0=1/2π,e1=cos⁡x/π,e2=sin⁡x/πe_0=1/\sqrt{2\pi},e_1=\cos x/\sqrt\pi,e_2=\sin x/\sqrt\pi は (ei,ej)=δij(e_i,e_j)=\delta_{ij} を満たすので

(f,g)=a0b0+a1b1+a2b2.\boxed{(f,g)=a_0b_0+a_1b_1+a_2b_2}.

1-2​

De0=0,De1=−e2,De2=e1De_0=0,De_1=-e_2,De_2=e_1 だから、O^1=2+D\hat O_1=2+D の行列表現は

A=(2000210−12).\boxed{A=\begin{pmatrix}2&0&0\\0&2&1\\0&-1&2\end{pmatrix}}.

1-3​

O^2=1+2D+2D2\hat O_2=1+2D+2D^2 の行列表現とその三乗は

B=(1000−120−2−1),B3=(100011−20211).B=\begin{pmatrix}1&0&0\\0&-1&2\\0&-2&-1\end{pmatrix},\qquad B^3=\begin{pmatrix}1&0&0\\0&11&-2\\0&2&11\end{pmatrix}.

反対称な非対角成分は二次形式で打ち消し合うので

(f,O^23f)=a02+11a12+11a22.\boxed{(f,\hat O_2^3f)=a_0^2+11a_1^2+11a_2^2}.

1-4​

2D4+D22D^4+D^2 の行列表現は diag⁡(0,1,1)\operatorname{diag}(0,1,1) である。境界条件から f≠0f\ne0 なので λ=0\lambda=0 または 11。

λ=0\lambda=0 では a1=a2=0a_1=a_2=0、f(π)=−1f(\pi)=-1 から a0=−2πa_0=-\sqrt{2\pi}。

λ=1\lambda=1 では a0=0a_0=0 であり、f(π)=−a1/π=−1f(\pi)=-a_1/\sqrt\pi=-1、f′(−π)=−a2/π=0f'(-\pi)=-a_2/\sqrt\pi=0 から a1=π,a2=0a_1=\sqrt\pi,a_2=0。

したがって

(λ,a0,a1,a2)=(0,−2π,0,0),(1,0,π,0).\boxed{(\lambda,a_0,a_1,a_2)=(0,-\sqrt{2\pi},0,0),\quad(1,0,\sqrt\pi,0)}.

1-5​

誤差を各係数で微分すると、最小点では h−gnh-g_n が 1,cos⁡kx,sin⁡kx1,\cos kx,\sin kx に直交する。三角関数系の直交性より

α0=1π∫−ππ∣x∣ dx=π,βk=0.\alpha_0=\frac1\pi\int_{-\pi}^{\pi}|x|\,dx=\boxed{\pi}, \qquad\boxed{\beta_k=0}.

さらに部分積分により

αk=2π∫0πxcos⁡kx dx=2((−1)k−1)πk2={−4/(πk2)k が奇数,0k が偶数.\begin{gathered} \alpha_k=\frac2\pi\int_0^\pi x\cos kx\,dx\\ =\boxed{\frac{2((-1)^k-1)}{\pi k^2}}\\ =\begin{cases}-4/(\pi k^2)&k\text{ が奇数},\\0&k\text{ が偶数}.\end{cases} \end{gathered}

直交性から誤差は各係数について正の二次式であり、この停留点が一意な最小点である。gng_n は hh のフーリエ部分和で、h∈L2[−π,π]h\in L^2[-\pi,\pi] より gn→hg_n\to h が L2L^2 の意味で成立する。Parseval の等式によって

lim⁡n→∞(gn,gn)=(h,h)=2∫0πx2 dx=2π33.\boxed{\lim_{n\to\infty}(g_n,g_n)=(h,h) =2\int_0^\pi x^2\,dx=\frac{2\pi^3}{3}}.

Reference​