跳到主要内容

京都大学 工学研究科 機械工学群 2023年度 数学 問題1

Author​

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

Description​

実関数の内積を

(f1,f2)=∫−ππf1(x)f2(x) dx(f_1,f_2)=\int_{-\pi}^{\pi}f_1(x)f_2(x)\,dx

で定義する。実定数 a0,a1,a2,b0,b1,b2a_0,a_1,a_2,b_0,b_1,b_2 を用いて

f(x)=a02π+a1πcos⁡x+a2πsin⁡x,f(x)=\frac{a_0}{\sqrt{2\pi}}+\frac{a_1}{\sqrt\pi}\cos x+\frac{a_2}{\sqrt\pi}\sin x,
g(x)=b02π+b1πcos⁡x+b2πsin⁡xg(x)=\frac{b_0}{\sqrt{2\pi}}+\frac{b_1}{\sqrt\pi}\cos x+\frac{b_2}{\sqrt\pi}\sin x

とする。以下の問いに答えよ。

1-1 (f,g)(f,g) を計算せよ。

1-2 演算子 O^1=2+d/dx\hat O_1=2+d/dx に対して

(f,O^1g)=(a0 a1 a2)A(b0b1b2)(f,\hat O_1g)=(a_0\ a_1\ a_2)A \begin{pmatrix}b_0\\b_1\\b_2\end{pmatrix}

を満たす 3×33\times3 実行列 AA を具体的な値で示せ。

1-3 O^2=1+2d/dx+2d2/dx2\hat O_2=1+2d/dx+2d^2/dx^2 について (f,O^23f)(f,\hat O_2^3f) を計算せよ。

1-4 次の微分方程式を考える。

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

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

1-5 実定数 α0,…,αn,β1,…,βn\alpha_0,\ldots,\alpha_n,\beta_1,\ldots,\beta_n により

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

を定義する。周期 2π2\pi の関数 h(x)=∣x∣ (−π≤x<π)h(x)=|x|\ (-\pi\le x<\pi)、h(x+2π)=h(x)h(x+2\pi)=h(x) との自乗平均誤差

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

を最小にする α0,αk,βk (k=1,…,n)\alpha_0,\alpha_k,\beta_k\ (k=1,\ldots,n) を求めよ。これらはフーリエ係数となる。また lim⁡n→∞(gn,gn)\lim_{n\to\infty}(g_n,g_n) を計算し、その過程を簡潔に説明せよ。

题目描述​

实函数内积定义为 (f1,f2)=∫−ππf1f2 dx(f_1,f_2)=\int_{-\pi}^{\pi}f_1f_2\,dx。函数 f,gf,g 在正交归一基 1/2π,cos⁡x/π,sin⁡x/π1/\sqrt{2\pi},\cos x/\sqrt\pi,\sin x/\sqrt\pi 下的系数分别为 (a0,a1,a2)(a_0,a_1,a_2) 和 (b0,b1,b2)(b_0,b_1,b_2),系数均为实数。

1-1:求 (f,g)(f,g)。

1-2:求微分算符 O^1=2+d/dx\hat O_1=2+d/dx 的矩阵 AA,使 (f,O^1g)=aTAb(f,\hat O_1g)=a^{\mathsf T}Ab。

1-3:令 O^2=1+2d/dx+2d2/dx2\hat O_2=1+2d/dx+2d^2/dx^2,求 (f,O^23f)(f,\hat O_2^3f)。

1-4:在上述三维函数空间内,求满足 (2d4/dx4+d2/dx2)f=λf(2d^4/dx^4+d^2/dx^2)f=\lambda f 及 f(π)=−1,f′(−π)=0f(\pi)=-1,f'(-\pi)=0 的所有 λ,a0,a1,a2\lambda,a_0,a_1,a_2。

1-5:用三角多项式 gn=α0/2+∑k=1n(αkcos⁡kx+βksin⁡kx)g_n=\alpha_0/2+\sum_{k=1}^n(\alpha_k\cos kx+\beta_k\sin kx) 逼近周期延拓的 h(x)=∣x∣h(x)=|x|。求使均方误差 (2π)−1∫−ππ(h−gn)2 dx(2\pi)^{-1}\int_{-\pi}^{\pi}(h-g_n)^2\,dx 最小的全部系数,并求 lim⁡n→∞(gn,gn)\lim_{n\to\infty}(g_n,g_n),简述推导。

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 が偶数.\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}

直交性から誤差は各係数について正の二次式であり、この停留点が一意な最小点である。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​