跳到主要内容

東京工業大学 工学院 電気電子系 2021年8月実施 数学1

Author

祭音Myyura

Description

\int_{-\infty}^{\infty} が有限で絶対積分可能な実関数 f(t)f(t) に対して,フーリエ変換 F(ω)F(\omega)

F(ω)=f(t)ejωtdtF(\omega) = \int_{-\infty}^{\infty}f(t)e^{-j\omega t}dt

とする。ただし,jj を虚数単位とし,関数 f(t)f(t) ,そのフーリエ変換 F(ω)F(\omega) は,22 階微分可能とする。以下の問に答えよ。なお,導出過程も示すこと。また,関数 f(t)f(t) の導関数 d2f(t)dt2\frac{d^2f(t)}{dt^2},df(t)dt\frac{df(t)}{dt} をそれぞれ f(t),f(t)f''(t),f'(t) と表記してもよい。

(1) 次の関数 g(t)g(t) のフーリエ変換 G(ω)G(\omega) を関数 f(t)f(t) のフーリエ変換 F(ω)F(\omega) を用いて表せ。

(a)g(t)=d2f(t)dt2(b)g(t)=t2f(t)(a)\quad g(t) = \frac{d^2f(t)}{dt^2} \qquad (b) \quad g(t) = t^2f(t)

(2) 関数 f(t)f(t) が次の微分方程式の解であるとき,関数 f(t)f(t) のフーリエ変換 F(ω)F(\omega) が満足する微分方程式を示せ。ただし,aa は実数の定数である。

d2f(t)dt2t2f(t)=af(t)\frac{d^2f(t)}{dt^2} - t^2f(t) = af(t)

(3) 関数 f(t)f(t) が次の微分方程式の解であるとする。ただし,nn00 以上の整数である。

d2f(t)dt2t2f(t)=(2n+1)f(t)\frac{d^2f(t)}{dt^2} - t^2f(t) = -(2n + 1)f(t)

(a) 関数 f(t)f(t)f(t)=Hn(t)et22f(t) = H_n(t)e^{-\frac{t^2}{2}} と表せるとき,関数 Hn(t)H_n(t) が満足する微分方程式を示せ。

(b) 上記の (a)(a) で求めた微分方程式の特殊解 Hn(t)=(1)net2dndtn(et2)H_n(t) = (-1)^ne^{t^2}\frac{d^n}{dt^n}(e^{t^2}) は,nn 次のエルミート多項式とよばれる。Hn(t),Hn+1(t),Hn+2(t)H_n(t),H_{n + 1}(t),H_{n + 2}(t) が満足する漸化式を求めよ。

题目描述

对绝对可积的实函数 f(t)f(t),定义 Fourier 变换

F(ω)=f(t)ejωtdt,F(\omega)=\int_{-\infty}^{\infty} f(t)e^{-j\omega t}\,\mathrm dt,

其中 jj 为虚数单位。设 f(t)f(t) 及其 Fourier 变换 F(ω)F(\omega) 均二阶可微。回答下列问题并写出推导过程;可把 df/dt\mathrm df/\mathrm dtd2f/dt2\mathrm d^2f/\mathrm dt^2 简记为 f,ff',f''

  1. 分别用 F(ω)F(\omega) 表示下列 g(t)g(t) 的 Fourier 变换 G(ω)G(\omega)

    1. g(t)=d2f(t)dt2;g(t)=\frac{\mathrm d^2f(t)}{\mathrm dt^2};
    2. g(t)=t2f(t).g(t)=t^2f(t).
  2. f(t)f(t) 满足

    d2f(t)dt2t2f(t)=af(t),\frac{\mathrm d^2f(t)}{\mathrm dt^2} -t^2f(t)=af(t),

    其中 aa 为实常数,求 F(ω)F(\omega) 所满足的微分方程。

  3. nn 为非负整数,且 f(t)f(t) 满足

    d2f(t)dt2t2f(t)=(2n+1)f(t),\frac{\mathrm d^2f(t)}{\mathrm dt^2} -t^2f(t)=-(2n+1)f(t),

    回答:

    1. f(t)=Hn(t)et2/2f(t)=H_n(t)e^{-t^2/2}

      时,求 Hn(t)H_n(t) 所满足的微分方程。

    2. 原 Description 把上述方程的特殊解

      Hn(t)=(1)net2dndtn(et2)H_n(t)=(-1)^ne^{t^2} \frac{\mathrm d^n}{\mathrm dt^n}(e^{t^2})

      称为 nn 次 Hermite 多项式。按这一题面定义,求 Hn(t),Hn+1(t),Hn+2(t)H_n(t),H_{n+1}(t),H_{n+2}(t) 满足的递推式。

考点

  • Fourier 变换的微分与时间乘法性质:通过分部积分把 f(t)f''(t) 化为频域乘法,把 t2f(t)t^2f(t) 化为频域二阶微分。
  • 微分方程的 Fourier 对偶:将含 t2t^2 的二阶方程整体变换,得到 F(ω)F(\omega) 的对应常微分方程。
  • Hermite 方程与多项式递推:代入 Gaussian 因子消去二次势项,并从给定 Rodrigues 型表达式或微分方程推导相邻阶关系。

Kai

(1)

(a)

g(t)=fn(t)G(ω)=f(t)ejωtdt=[f(t)ejωt]f(t)(jω)ejωtdt=jωf(t)ejωtdt=jω[f(t)ejωt]jωf(t)(jω)ejωtdt=(jω)2f(t)ejωtdt=ω2F(ω)\begin{aligned} g(t) &= f^{n}(t) \\ G(\omega) &= \int_{-\infty}^{\infty} f''(t)e^{j\omega t}dt \\ &= \bigg[f'(t)e^{-j\omega t}\bigg]_{-\infty}^{\infty} - \int_{-\infty}^{\infty} f'(t)(-j\omega)e^{-j\omega t}dt \\ &= j\omega \int_{-\infty}^{\infty}f'(t)e^{-j\omega t}dt \\ &= j\omega \bigg[f(t)e^{-j\omega t}\bigg]_{-\infty}^{\infty} - j\omega \int_{-\infty}^{\infty}f(t)(-j\omega)e^{-j\omega t}dt \\ &= (j\omega)^2 \int_{-\infty}^{\infty}f(t)e^{-j\omega t}dt \\ &= -\omega^2F(\omega) \end{aligned}

(b)

g(t)=t2f(t)F(ω)=f(t)ejωtdtdF(ω)dω=jtf(t)ejωtdtd2F(ω)dω2=t2f(t)ejωtdtt2f(t)ejωtdt=d2F(ω)dω=G(ω)\begin{aligned} g(t) &= t^2f(t) \\ F(\omega) &= \int_{-\infty}^{\infty}f(t)e^{-j\omega t}dt \\ \frac{dF(\omega)}{d\omega} &= -j\int_{-\infty}^{\infty} tf(t) e^{-j\omega t}dt \\ \frac{d^2F(\omega)}{d\omega^2} &= -\int_{-\infty}^{\infty}t^2f(t)e^{-j\omega t}dt \\ \int_{-\infty}^{\infty}t^2f(t)e^{-j\omega t}dt &= -\frac{d^2F(\omega)}{d\omega} = G(\omega) \end{aligned}

従って、

G(ω)=d2F(ω)dω2G(\omega) = -\frac{d^2F(\omega)}{d\omega^2}

(2)

d2f(t)dt2t2f(t)=af(t)ω2F(ω)(d2F(ω)dω3)=aF(ω)d2F(ω)dω2(ω2+a)F(ω)=0\begin{aligned} \frac{d^2f(t)}{dt^2} - t^2f(t) &= af(t) \\ -\omega^2F(\omega) - (-\frac{d^2F(\omega)}{d\omega^3}) &= aF(\omega) \\ \frac{d^2F(\omega)}{d\omega^2} - (\omega^2 + a)F(\omega) &= 0 \\ \end{aligned}

(3)

(a)

d2dt2(Hn(t)et22)t2Hn(t)et22=(2n+1)Hn(t)et22\begin{align} \frac{d^2}{dt^2}\big(H_n(t)e^{-\frac{t^2}{2}}\big) - t^2H_n(t)e^{-\frac{t^2}{2}} = -(2n + 1)H_n(t)e^{-\frac{t^2}{2}} \tag{i} \end{align}
d2dt2Hn(t)et22=ddt(dHn(t)dtet22Hn(t)tet22)=d2Hn(t)dt2et22tdHn(t)dtet22tHn(t)dtet22Hn(t)et22+Hn(t)t2et22\begin{aligned} \frac{d^2}{dt^2}H_n(t)e^{-\frac{t^2}{2}} &= \frac{d}{dt}\bigg(\frac{dH_n(t)}{dt}e^{-\frac{t^2}{2}} - H_n(t)te^{-\frac{t^2}{2}}\bigg) \\ &= \frac{d^2H_n(t)}{dt^2}e^{-\frac{t^2}{2}} - t\frac{dH_n(t)}{dt}e^{-\frac{t^2}{2}} - t\frac{H_n(t)}{dt}e^{-\frac{t^2}{2}} - H_n(t)e^{-\frac{t^2}{2}} + H_n(t)t^2e^{-\frac{t^2}{2}}\\ \end{aligned}

(i)(\text{i}) 式へ代入すると、

d2Hn(t)dt2et222tdHn(t)dtet22+2nHn(t)et22=0\frac{d^2H_n(t)}{dt^2}e^{-\frac{t^2}{2}} - 2t\frac{dH_n(t)}{dt}e^{-\frac{t^2}{2}} + 2nH_n(t)e^{-\frac{t^2}{2}} = 0
d2Hn(t)dt22tdHn(t)dt+2nHn(t)=0\begin{align} \frac{d^2H_n(t)}{dt^2} - 2t\frac{dH_n(t)}{dt} + 2nH_n(t) = 0 \tag{ii} \end{align}

(b)

dHn(t)dt=(1)n2tet2dndtn(et2)=2tHn(t)Hn+1(t)\begin{aligned} \frac{dH_n(t)}{dt} &= (-1)^n \cdot 2te^{t^2} \frac{d^n}{dt^n}(e^{t^2}) \\ &= 2tH_n(t) - H_{n+1}(t) \end{aligned}
d2Hn(t)dt2=(1)n2et2dndtn(et2)+(1)n4t2et2dndtn(et2)2t(1)n+1et2dn+1dtn+1(et2)+(1)n+2et2dn+2dtn+2(et2)2t(1)n+1et2dn+1dtn+1(et2)\begin{aligned} \frac{d^2H_n(t)}{dt^2} &= (-1)^n \cdot 2e^{t^2} \frac{d^n}{dt^n}(e^{t^2}) + (-1)^n 4t^2e^{t^2} \frac{d^n}{dt^n}(e^{t^2}) - 2t(-1)^{n + 1}e^{t^2} \frac{d^{n + 1}}{dt^{n + 1}}(e^{t^2}) + (-1)^{n + 2}e^{t^2}\frac{d^{n + 2}}{dt^{n + 2}}(e^{t^2}) - 2t(-1)^{n + 1}e^{t^2}\frac{d^{n + 1}}{dt^{n + 1}}(e^{t^2}) \end{aligned}
d2Hn(t)dt=2Hn(t)+4t2Hn(t)2tHn+1(t)+Hn+2(t)2tHn+1(t)=2Hn(t)+4t2Hn(t)4tHn+1(t)+Hn+2(t)\begin{aligned} \frac{d^2H_n(t)}{dt} &= 2H_n(t) + 4t^2H_n(t) - 2tH_{n+1 }(t) + H_{n + 2}(t) - 2tH_{n + 1}(t) \\ &= 2H_n(t) + 4t^2H_n(t) - 4tH_{n + 1}(t) + H_{n + 2}(t) \end{aligned}

(i)(\text{i}) 式に代入すると、

2Hn(t)+4t2Hn(t)4tHn+1(t)+Hn+2(t)2t(2tHn(t)Hn+1(t))+2nHn(t)=02H_n(t) + 4t^2H_n(t) - 4tH_{n + 1}(t) + H_{n + 2}(t) - 2t\big(2tH_n(t) - H_{n + 1}(t)\big) + 2nH_n(t) = 0
2(1+n)Hn(t)2tHn+1(t)+Hn+2(t)=0\begin{aligned} 2(1 + n)H_n(t) - 2tH_{n + 1}(t) + H_{n + 2}(t) = 0 \end{aligned}

従って、

Hn+2(t)2tHn+1(t)+2(1+n)Hn(t)=0H_{n + 2}(t) - 2tH_{n + 1}(t) + 2(1 + n)H_n(t) = 0