跳到主要内容

大阪大学 情報科学研究科 情報工学 2019年度 数学解析と信号処理

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

(1-1) Γ(z)=0xz1exdx\Gamma(z)=\int_0^\infty x^{z-1}e^{-x}\,dxz>0\Re z>0)について Γ(1/2)=π\Gamma(1/2)=\sqrt\pi を示せ。

(1-2) t>0t>0ss を複素変数とする。f(t)=1/tf(t)=1/\sqrt t のLaplace変換が π/s\sqrt{\pi/s} であることを示せ。

(1-3) G(s)=1/(s(s1))G(s)=1/(\sqrt{s}(s-1)) の逆Laplace変換を求めよ。

(2) nn を整数、zz を複素変数とし、離散時間の線形時不変システムを考える。

(2-1) x(n)=δ(n)δ(n1)+2δ(n2)x(n)=\delta(n)-\delta(n-1)+2\delta(n-2)zz 変換を求めよ。δ(n)\delta(n) は単位インパルスを表す。

(2-2) y(n)=x(n)+2x(n1)3x(n2)y(n)=x(n)+2x(n-1)-3x(n-2) となる線形時不変システムの伝達関数を求めよ。

(2-3) H(z)=1+2z1+z2H(z)=1+2z^{-1}+z^{-2} の振幅特性と位相特性を求めよ。

Kai

(1)

(1-1) x=u2x=u^2 とすると Γ(1/2)=20eu2du=I\Gamma(1/2)=2\int_0^\infty e^{-u^2}\,du=I。2重積分を極座標に変換して

I2=R2e(u2+v2)dudv=2π0rer2dr=π.I^2=\int_{\mathbb R^2}e^{-(u^2+v^2)}\,du\,dv=2\pi\int_0^\infty re^{-r^2}\,dr=\pi.

I>0I>0 より Γ(1/2)=π\boxed{\Gamma(1/2)=\sqrt\pi}

(1-2) 実数 s>0s>0 では x=stx=st と置き

L[t1/2](s)=s1/2Γ(1/2)=π/s.\mathcal L[t^{-1/2}](s)=s^{-1/2}\Gamma(1/2)=\boxed{\sqrt{\pi/s}}.

両辺は s>0\Re s>0 で正則なので、恒等定理により同じ半平面で成立する(平方根は主値)。

(1-3) 畳込み定理により

L1[G](t)=1π0tetuudu=2etπ0tev2dv=eterf(t).\boxed{\mathcal L^{-1}[G](t)=\frac1{\sqrt\pi}\int_0^t\frac{e^{t-u}}{\sqrt u}\,du =\frac{2e^t}{\sqrt\pi}\int_0^{\sqrt t}e^{-v^2}\,dv=e^t\operatorname{erf}(\sqrt t)}.

(2)

(2-1) X(z)=1z1+2z2\boxed{X(z)=1-z^{-1}+2z^{-2}}

(2-2) Y(z)=(1+2z13z2)X(z)Y(z)=(1+2z^{-1}-3z^{-2})X(z) より H(z)=1+2z13z2\boxed{H(z)=1+2z^{-1}-3z^{-2}}

(2-3) z=eiωz=e^{i\omega} とすると

H(eiω)=(1+eiω)2=4cos2(ω/2)eiω.H(e^{i\omega})=(1+e^{-i\omega})^2=4\cos^2(\omega/2)e^{-i\omega}.

したがって

H(eiω)=4cos2(ω/2),argH(eiω)=ω(mod2π).\boxed{|H(e^{i\omega})|=4\cos^2(\omega/2),\qquad\arg H(e^{i\omega})=-\omega\pmod{2\pi}}.

ω=(2k+1)π\omega=(2k+1)\pi では振幅0のため位相は定義されない。