• ベストアンサー

sech関数のフーリエ変換

aquarius_hiroの回答

  • ベストアンサー
回答No.1

sparklightさん、こんにちは。  F(ω) = ∫_{-∞}^{∞} dx exp[iωx] sech(ax) と定義すると、  sech(ax) = 1/cosh(ax) = 2/(exp[ax] + exp[-ax]) より、z=axとおいて、  F(ω) = (2/a) ∫_{-∞}^{∞} dz exp[i(ω/a)z] /(exp[z] + exp[-z]) … (1) ωを-ωに置換えても、積分変数zを-zに置き換えれば等しくなるので、これは偶関数です。 従って、ω/a > 0 と仮定して求めれば十分です。 (ω/a < 0 のケースは同じになります。) (1)を複素平面上の積分とみます。 ω/a > 0 と、exp[i(ω/a)z] の因子により、上半無限円上の同じ関数の積分は、0になります。 (簡単に言うと、z=x+iy とおくとき、上半無限円では、y→+∞としてよく、 指数関数 exp[i(ω/a)z] の肩は、i(ω/a)(x+iy) = i(ω/a) x - (ω/a)y となり、その実部が -(ω/a)y → -∞になるからです。) 従って、(1)は、  F(ω) = (2/a) ∫_C dz exp[i(ω/a)z] /(exp[z] + exp[-z]) … (2) という複素積分になります。ここで、Cは実軸に上半無限円を加えた閉曲線です。  1/(exp[z] + exp[-z]) = exp[z]/(exp[2z] + 1) は、この閉曲線の内部に1次の極を持ちます。すなわち、  exp[z] = - exp[-z]  exp[2z] = -1  2z = (2n+1)iπ  z = (2n+1)iπ/2 = inπ + iπ/2 ≡ z_n とおく が極です。n=0,1,2,… で、積分経路Cの内側には無限個の極があります。 この極以外では正則なので、経路を縮小すると、  F(ω) = (2/a) Σ_{n=0}^{∞} ∫_{C_n} dz exp[i(ω/a)z] /(exp[z] + exp[-z]) ここで、C_n は、z=z_n を囲む微小な閉曲線です。z=z_n の近傍では  1/(exp[z] + exp[-z]) = exp[z]/[exp(2z) + 1]   = exp[z]/[exp(2(z-z_n)+z_n) + 1]   = exp[z]/[exp(z_n) exp(2(z-z_n)) + 1]   = exp[z_n]/[- {1 + 2(z-z_n)} + 1]   = - i (-1)^n/2 × 1/(z-z_n) と展開できるので、これに留数定理を使うと、  F(ω) = (2/a) Σ_{n=0}^{∞} (-i) (-1)^n exp[i(ω/a)z_n]/2 × 2πi     = (2π/a) Σ_{n=0}^{∞} (-1)^n exp[-(ω/a)(2n+1)π/2]     = (2π/a) exp[-(ω/a)π/2]/[1 + exp[-(ω/a)π]]     = (2π/a) 1/[exp[(ω/a)π/2] + exp[-(ω/a)π/2]]     = (π/a) sech[πω/2a] と求まります。ω/a < 0 のほうも同じになることは上で述べたとおりです。

sparklight
質問者

お礼

大変細かく丁寧な解答をしていただき、ありがとうございました。 ただ、以下の部分からの計算が、自分の知識が足りないせいかわかりません。。 もしよろしければ以下の点について教えていただけますでしょうか? > = exp[z]/[exp(2(z-z_n)+z_n) + 1] > = exp[z]/[exp(z_n) exp(2(z-z_n)) + 1] > = exp[z_n]/[- {1 + 2(z-z_n)} + 1] 1行目のところですが、なぜexp[z]/[exp(2(z-z_n)+2z_n) + 1]にはならないのでしょうか?単純に式変形しただけのものではないのでしょうか? また、2行目から3行目への変形についても教えていただきたいです。 よろしくお願いします。

関連するQ&A

  • 逆フーリエ変換

    フーリエ変換の問題を解いていて f(x)=1/2a(|x|<a) , 0(|x|>a)  をフーリエ変換したら sin(ax)/ax となったのですが、これを逆フーリエ変換したらf(x)になるはずですよね? 公式にあてはめて何度計算しても収束させることができず積分がうまくできないのですが、どうやればいいのでしょうか。。

  • フーリエ変換について

    フーリエ変換表などと載ってる「f(ax)→F(ν/a)/|a|」 を証明したいのですが、どうしても解けません。 f(x)の変換後をF(ν)、変換をZとして Z=∫f(ax)・exp(-2πiux)dx にy=ax、x=y/aとして代入して解いているのですが そもそもこのやり方が間違っているのか、 計算途中でミスしているのか、わかりません。 教えてください、お願いします。

  • cos(ax)をフーリエ変換する問題

    単なるcos(ax) でなく、 x<X0でf(x)=0、x≧x0 でf(x)=cos(ax)という関数ののフーリエ変換です。 (ただし ax0=π/2 です) ただのcos(ax) のフーリエ変換は、 ∫[-∞、+∞] {exp(iax) + exp(-iax) } exp(-ikx) dx =1/√2π {δ(k-a) + δ(k+a)} なので、 この半分に exp(ik何とか) を掛けたものかとなぁ と思うのですが、わかりません。 アドバイスをお願いします。

  • フーリエ変換の問題。助けてください!!

    フーリエ変換の問題です。 関数f(x)のフーリエ変換F(ω)および逆変換を次のように定義する。 F(ω)=∫[∞→-∞]f(x)e^(-iωx)dx f(x)=(1/2π)∫[∞→-∞]F(ω)e^iωx dω 関数f(x)=e^(-ax^2)のフーリエ変換はF(ω)=(√π/a)e^(-ω^2/(4a)) で与えられる。以下の問いに答えよ。 (1)∫[∞→-∞]e^((-x^2)/2)dx を求めよ。 (2)関数g(x)が次の方程式を満足する。 ∫[∞→-∞]∫[∞→-∞]g(x-y)g(y-z)g(z)dydz=(2π/√3)e^((-x^2)/6) g(x)を求めよ。 (1)はすぐに√2πであるとわかりました。(2)は朝からずっと考えていろいろと試行錯誤してみたのですが、まったく解けません。 どなたか、教えていただけませんか?私も引き続き考えますのでどうかわかる方、よろしくお願い致します。

  • フーリエ変換の問題

    f(x)のフーリエ変換F(x)は F(x)=∫[-∞,∞]f(x)exp(-iωx)dx で表される。次の関数のフーリエ変換を求めよ。 a>0として、 f(x)={0 (x>0 )            {-exp(ax) (x<=0) という問題があります。 私は F(x)=(iω-a)^(-1)[exp(a-iω)x](-∞→0)までやりました。 ここで、ちょっとわからないところがあります。 exp((a-iω)x)の値はx=0のときは1ですよね。 でも、x=-∞のときは、どうすればいいかわからなくなりました。 普通なら、exp(ax)=0ですよね。a>0(つまりaは0より大きければ),x=-∞なら ですが、a-iωは虚数であって0と比べられないですよね。ちょっとここでつまづいて... 問題の答えをみればexp((a-iω)x)=0 (x=-∞)って書いてありますけど、なぜそうなるか書いてないです。 ご指導お願いします!

  • フーリエ変換の問題

    f(x)=e^(-ax^2)   a(-∞≦x≦∞,a>0) のフーリエ変換が分かる方いましたら是非教えてください。よろしくお願いします。

  • フーリエ変換の問題について

    f(x)=e^(-ax^2)  (-∞≦x≦∞,a>0) のフーリエ変換が分かる方いましたら是非教えてください。よろしくお願いします。

  • フーリエ変換(スペクトル解析について)

    解き方を教えて下さい。よろしくお願いします。 g(x)={1-|x| (|x|≦1) 0 (|x|>1) のときスペクトルを求めよ。 答え・・・Sg(ω)= 2(1-cosω)/(πω^2) (ω>0) 1/π (ω = 0) g(x)のフーリエ変換は G(u)=2(1-cos u)/u^2 反転公式を適用し g(x)=1/π∫(0→∞) 2(1-cos ω/ω^2 ) cos ωx dω 角周波数にたいしてcos ωxの成分は2(1-cos ω/ω^2 ) であり求めるスペクトルは ω>0において Sg(ω) = 2(1-cosω)/(πω^2) というところまでは分かるのですが、ω=0の場合は どのように考えるのでしょうか?単純にω=0を代入すると 分母が0になってしまいます。 同様の問題で、 1 (0≦x<1)) f(x)= -1 (-1<x<0) 0 (|x|≧1) のフーリエ変換に反転公式を適用して連続スペクトルを出すという問題で、 これは奇関数なので正弦変換を用い、 F(u)=-2i∫(0→1) sin ux dx = 2(cos u -1) i/uとなり 反転公式f(x) = i/π∫(0→∞) F(u) sin ux duに代入し f(x) = i/π ∫(0→∞) (2(cos u -1) i/u) sin ux dx より、スペクトルはSf(ω)= (2i/π) * (cos ω -1) i/ω = - 2(cos ω -1)/(πω) (ω>1)? となりましたが、この答えは合っているでしょうか?また、ω=0の場合のスペクトルの式は どのようになるのでしょうか?

  • フーリエ変換

    f(x)=e^(-ax^2) のフーリエ変換を求めたいのですが、良く分かりません…。 お願いします。

  • フーリエ変換

    f(x)=e(-ax) =0 のフーリエ変換のやり方を教えてください。