定義 2000Hz(ゼロを中心とした±1.5秒)の音声におけるsinc関数 sinc関数には、正規化された形式と正規化されていない形式の2種類があります。[ 1 ]
数学では、正規化されていない歴史的なsinc関数は、 x ≠ 0 に対して次のよう に定義されます。シン ( x ) = 罪 x x 。 {\displaystyle \operatorname {sinc} (x)={\frac {\sin x}{x}}.}
あるいは、正規化されていないsinc関数は、サンプリング関数 と呼ばれることが多く、Sa( x )と表記される。[ 2 ]
デジタル信号処理 と情報理論 では、正規化されたsinc関数は、 x ≠0 の 場合、一般的に次のように定義されます。シン ( x ) = 罪 ( π x ) π x 。 {\displaystyle \operatorname {sinc} (x)={\frac {\sin(\pi x)}{\pi x}}.}
いずれの場合も、 x = 0 の値は極限値として定義される。 シン ( 0 ) := リム x → 0 罪 ( 1 x ) 1 x = 1 {\displaystyle \operatorname {sinc} (0):=\lim _{x\to 0}{\frac {\sin(ax)}{ax}}=1} すべての実数a ≠ 0 に対して(極限は、はさみうちの定理 を使用して証明できます)。
正規化 により、実数全体にわたる関数の定積分は 1 になります (正規化されていない sinc 関数の同じ積分は π になります)。さらに便利な性質として、正規化された sinc 関数の零点は、x の非零整数値になります。
語源 この関数は、カーディナルサイン またはサインカーディナル 関数とも呼ばれています。[ 3 ] [ 4 ] 「sinc」という用語は、この関数の完全なラテン語名sinus cardinalis [ 5 ] の短縮形であり、 Philip M. Woodward と IL Davies が 1952 年の論文「情報理論と通信における逆確率 」で導入し、「この関数はフーリエ解析とその応用で非常に頻繁に現れるため、独自の表記法に値すると思われる」と述べています。[ 6 ] また、Woodward の 1953 年の著書「確率と情報理論、レーダーへの応用」 でも使用されています。[ 5 ] [ 7 ]
不動産 正規化されていない赤いsinc関数の局所的な極値(小さな白い点)は、青いcos関数との交点に対応します。 正規化されていないsinc関数のゼロ交差は π の非ゼロの整数倍で発生するが、正規化されたsinc関数のゼロ交差は非ゼロの整数で発生する。
正規化されていないsinc関数の局所極値は、 コサイン 関数との交点に対応します。つまり、sin ( ξ ) / ξ = cos( ξ ) がすべての点ξに対して成り立ち、 sin ( x ) / x の導関数がゼロになるため、 局所極値に達します。これはsinc関数の導関数から導かれます。 d d x シン ( x ) = { コス ( x ) − シン ( x ) x 、 x ≠ 0 0 、 x = 0 。 {\displaystyle {\frac {d}{dx}}\operatorname {sinc} (x)={\begin{cases}{\dfrac {\cos(x)-\operatorname {sinc} (x)}{x}},&x\neq 0\\0,&x=0\end{cases}}.}
正のx座標を持つ n 番目の極値のx 座標に関する無限級数の最初の数項は次のとおりである。x n = q − q − 1 − 2 3 q − 3 − 13 15 q − 5 − 146 105 q − 7 − ⋯ 、 {\displaystyle x_{n}=qq^{-1}-{\frac {2}{3}}q^{-3}-{\frac {13}{15}}q^{-5}-{\frac {146}{105}}q^{-7}-\cdots ,} どこ q = ( n + 1 2 ) π 、 {\displaystyle q=\left(n+{\frac {1}{2}}\right)\pi ,} また、奇数n の場合は局所最小値、偶数n の 場合は局所最大値となる。y 軸に関する対称性により、 x 座標− x n に極値が存在する。さらに、 ξ 0 = (0, 1) に絶対最大値が存在する。
正規化されたsinc関数は、無限積 として単純に表現できます。 罪 ( π x ) π x = ∏ n = 1 ∞ ( 1 − x 2 n 2 ) {\displaystyle {\frac {\sin(\pi x)}{\pi x}}=\prod _{n=1}^{\infty }\left(1-{\frac {x^{2}}{n^{2}}}\right)}
複素平面上に−2 − 2 i から2 + 2 i までプロットされた基本正弦関数sinc( z ) そして、オイラーの反射公式 を通して、ガンマ関数 Γ( x ) およびガウスのπ関数 と関連付けられる。 罪 ( π x ) π x = 1 Γ ( 1 + x ) Γ ( 1 − x ) = 1 Π ( x ) Π ( − x ) 。 {\displaystyle {\frac {\sin(\pi x)}{\pi x}}={\frac {1}{\Gamma (1+x)\Gamma (1-x)}}={\frac {1}{\Pi (x)\Pi (-x)}}.}
オイラー は[ 8 ] 罪 ( x ) x = ∏ n = 1 ∞ コス ( x 2 n ) 、 {\displaystyle {\frac {\sin(x)}{x}}=\prod _{n=1}^{\infty }\cos \left({\frac {x}{2^{n}}}\right),} そして、積和恒等式[ 9 ]により
sinc z = sin z / z の領域彩色 図∏ n = 1 k コス ( x 2 n ) = 1 2 k − 1 ∑ n = 1 2 k − 1 コス ( n − 1 / 2 2 k − 1 x ) 、 ∀ k ≥ 1 、 {\displaystyle \prod _{n=1}^{k}\cos \left({\frac {x}{2^{n}}}\right)={\frac {1}{2^{k-1}}}\sum _{n=1}^{2^{k-1}}\cos \left({\frac {n-1/2}{2^{k-1}}}x\right),\quad \forall k\geq 1,} オイラーの積は和として書き換えることができる 罪 ( x ) x = リム N → ∞ 1 N ∑ n = 1 N コス ( n − 1 / 2 N x ) 。 {\displaystyle {\frac {\sin(x)}{x}}=\lim _{N\to \infty }{\frac {1}{N}}\sum _{n=1}^{N}\cos \left({\frac {n-1/2}{N}}x\right).}
正規化されたsinc関数(通常の周波数)の連続フーリエ変換は rect ( f ) です。 ∫ − ∞ ∞ シン ( t ) e − 私 2 π f t d t = 長方形 ( f ) 、 {\displaystyle \int _{-\infty }^{\infty }\operatorname {sinc} (t)\,e^{-i2\pi ft}\,dt=\operatorname {rect} (f),} ここで、矩形関数は 引数 が− 1/2 から 1/2 の 間の場合に 1 となり、それ以外の場合は 0 と なり ます。これは 、 sinc フィルタ が理想的な (ブリックウォール 、つまり矩形周波数応答 )ローパスフィルタ であるという事実に対応します。
このフーリエ積分は、特別な場合を含めて ∫ − ∞ ∞ 罪 ( π x ) π x d x = 長方形 ( 0 ) = 1 {\displaystyle \int _{-\infty }^{\infty }{\frac {\sin(\pi x)}{\pi x}}\,dx=\operatorname {rect} (0)=1} これは広義積分 (ディリクレ積分を 参照) であり、収束するルベーグ積分 ではない。 ∫ − ∞ ∞ | 罪 ( π x ) π x | d x = + ∞ 。 {\displaystyle \int _{-\infty }^{\infty }\left|{\frac {\sin(\pi x)}{\pi x}}\right|\,dx=+\infty .}
正規化されたsinc関数は、サンプリングされ た 帯域制限 関数の補間において理想的な特性を備えています。
これは補間関数であり、sinc(0) = 1 、非ゼロ整数 k に対してはsinc( k ) = 0 となります 。 関数x k ( t ) = sinc( t − k ) ( k は 整数) は、 最高角周波数ω H = π (つまり、最高サイクル周波数f H = 1 / 2 ) を持つ関数 空間 L 2 ( R )の帯域制限関数の正規直交基底を 形成します。 これら2つのsinc関数のその他の特性は以下のとおりです。
正規化されていないsinc関数は、第1種零次球面ベッセル関数 j 0 ( x ) です。正規化されたsinc関数はj 0 (π x ) です。 ここでSi( x )は 正弦積分 であり、∫ 0 x 罪 ( θ ) θ d θ = Si ( x ) 。 {\displaystyle \int _{0}^{x}{\frac {\sin(\theta )}{\theta }}\,d\theta =\operatorname {Si} (x).} λ sinc( λx ) (正規化されていない) は、線形常微分方程式の 2 つの線形独立解のうちの 1 つです。 x d 2 y d x 2 + 2 d y d x + λ 2 x y = 0. {\displaystyle x{\frac {d^{2}y}{dx^{2}}}+2{\frac {dy}{dx}}+\lambda ^{2}xy=0.} もう一つは cos( λx ) / x で、これは対応する sinc 関数とは異なり、x = 0で有界ではありません。 正規化されたsinc関数を使用して、∫ − ∞ ∞ 罪 2 ( θ ) θ 2 d θ = π ⇒ ∫ − ∞ ∞ シン 2 ( x ) d x = 1 、 {\displaystyle \int _{-\infty }^{\infty }{\frac {\sin ^{2}(\theta )}{\theta ^{2}}}\,d\theta =\pi \quad \Rightarrow \quad \int _{-\infty }^{\infty }\operatorname {sinc} ^{2}(x)\,dx=1,} ∫ − ∞ ∞ 罪 ( θ ) θ d θ = ∫ − ∞ ∞ ( 罪 ( θ ) θ ) 2 d θ = π 。 {\displaystyle \int _{-\infty }^{\infty }{\frac {\sin(\theta )}{\theta }}\,d\theta =\int _{-\infty }^{\infty }\left({\frac {\sin(\theta )}{\theta }}\right)^{2}\,d\theta =\pi .} ∫ − ∞ ∞ 罪 3 ( θ ) θ 3 d θ = 3 π 4 。 {\displaystyle \int _{-\infty }^{\infty }{\frac {\sin ^{3}(\theta )}{\theta ^{3}}}\,d\theta ={\frac {3\pi }{4}}.} ∫ − ∞ ∞ 罪 4 ( θ ) θ 4 d θ = 2 π 3 。 {\displaystyle \int _{-\infty }^{\infty }{\frac {\sin ^{4}(\theta )}{\theta ^{4}}}\,d\theta ={\frac {2\pi }{3}}.} 以下の広義積分には、(正規化されていない)sinc関数が含まれています。∫ 0 ∞ d x x n + 1 = 1 + 2 ∑ k = 1 ∞ ( − 1 ) k + 1 ( k n ) 2 − 1 = 1 シン ( π n ) 。 {\displaystyle \int _{0}^{\infty }{\frac {dx}{x^{n}+1}}=1+2\sum _{k=1}^{\infty }{\frac {(-1)^{k+1}}{(kn)^{2}-1}}={\frac {1}{\operatorname {sinc} ({\frac {\pi }{n}})}}.}
ディラックのデルタ分布との関係 正規化されたsinc関数は、初期デルタ関数 として使用でき、次の弱極限 が成り立つことを意味します。
リム 1 → 0 罪 ( π x 1 ) π x = リム 1 → 0 1 1 シン ( x 1 ) = δ ( x ) 。 {\displaystyle \lim _{a\to 0}{\frac {\sin \left({\frac {\pi x}{a}}\right)}{\pi x}}=\lim _{a\to 0}{\frac {1}{a}}\operatorname {sinc} \left({\frac {x}{a}}\right)=\delta (x).}
これは通常の極限ではありません。なぜなら、左辺は収束しないからです。むしろ、それは次のことを意味します。
リム 1 → 0 ∫ − ∞ ∞ 1 1 シン ( x 1 ) φ ( x ) d x = φ ( 0 ) {\displaystyle \lim _{a\to 0}\int _{-\infty }^{\infty }{\frac {1}{a}}\operatorname {sinc} \left({\frac {x}{a}}\right)\varphi (x)\,dx=\varphi (0)}
フーリエ反転定理 からわかるように、すべてのシュワルツ関数 に対して、上記の式では、a → 0 のとき、sinc 関数の単位長さあたりの振動数は無限大に近づきます。しかしながら、式は常に± 1 / π x の包絡線内で振動し、 a の値に関係なく、その状態が続きます。
これは、 x = 0 以外のすべてのxに対して δ ( x ) がゼロであるという非公式なイメージを複雑にし、デルタ関数を分布ではなく関数として考えることの問題点を示しています。同様の状況はギブス現象 にも見られます。
また、標準的なディラック表現とすぐに関連付けることもできます。δ ( x ) {\displaystyle \delta (x)} 書くことによってb = 1 / 1 {\displaystyle b=1/a} そして
リム b → ∞ 罪 ( b π x ) π x = リム b → ∞ 1 2 π ∫ − b π b π e 私 k x d k = 1 2 π ∫ − ∞ ∞ e 私 k x d k = δ ( x ) 、 {\displaystyle \lim _{b\to \infty }{\frac {\sin \left(b\pi x\right)}{\pi x}}=\lim _{b\to \infty }{\frac {1}{2\pi }}\int _{-b\pi }^{b\pi }e^{ikx}dk={\frac {1}{2\pi }}\int _{-\infty }^{\infty }e^{ikx}dk=\delta (x),}
これは、デルタが積分の無限帯域幅極限として回復されることを明確に示している。
総和 このセクションにおけるすべての合計値は、正規化されていないsinc関数に関するものです。
1から∞ までの整数nに対する sinc( n ) の合計は、 π − 1 / 2 に等しい。
∑ n = 1 ∞ シン ( n ) = シン ( 1 ) + シン ( 2 ) + シン ( 3 ) + シン ( 4 ) + ⋯ = π − 1 2 。 {\displaystyle \sum _{n=1}^{\infty }\operatorname {sinc} (n)=\operatorname {sinc} (1)+\operatorname {sinc} (2)+\operatorname {sinc} (3)+\operatorname {sinc} (4)+\cdots ={\frac {\pi -1}{2}}.}
平方和も π − 1 / 2 に等しい: [ 10 ] [ 11 ]
∑ n = 1 ∞ シン 2 ( n ) = シン 2 ( 1 ) + シン 2 ( 2 ) + シン 2 ( 3 ) + シン 2 ( 4 ) + ⋯ = π − 1 2 。 {\displaystyle \sum _{n=1}^{\infty }\operatorname {sinc} ^{2}(n)=\operatorname {sinc} ^{2}(1)+\operatorname {sinc} ^{2}(2)+\operatorname {sinc} ^{2}(3)+\operatorname {sinc} ^{2}(4)+\cdots ={\frac {\pi -1}{2}}.}
加数 の 符号が交互になり、+ から始まる場合、 和は 1/2 になり ます 。 ∑ n = 1 ∞ ( − 1 ) n + 1 シン ( n ) = シン ( 1 ) − シン ( 2 ) + シン ( 3 ) − シン ( 4 ) + ⋯ = 1 2 。 {\displaystyle \sum _{n=1}^{\infty }(-1)^{n+1}\,\operatorname {sinc} (n)=\operatorname {sinc} (1)-\operatorname {sinc} (2)+\operatorname {sinc} (3)-\operatorname {sinc} (4)+\cdots ={\frac {1}{2}}.}
平方数と立方数 の 交互和も1/2 に 等しい : [ 12 ] ∑ n = 1 ∞ ( − 1 ) n + 1 シン 2 ( n ) = シン 2 ( 1 ) − シン 2 ( 2 ) + シン 2 ( 3 ) − シン 2 ( 4 ) + ⋯ = 1 2 、 {\displaystyle \sum _{n=1}^{\infty }(-1)^{n+1}\,\operatorname {sinc} ^{2}(n)=\operatorname {sinc} ^{2}(1)-\operatorname {sinc} ^{2}(2)+\operatorname {sinc} ^{2}(3)-\operatorname {sinc} ^{2}(4)+\cdots ={\frac {1}{2}},}
∑ n = 1 ∞ ( − 1 ) n + 1 シン 3 ( n ) = シン 3 ( 1 ) − シン 3 ( 2 ) + シン 3 ( 3 ) − シン 3 ( 4 ) + ⋯ = 1 2 。 {\displaystyle \sum _{n=1}^{\infty }(-1)^{n+1}\,\operatorname {sinc} ^{3}(n)=\operatorname {sinc} ^{3}(1)-\operatorname {sinc} ^{3}(2)+\operatorname {sinc} ^{3}(3)-\operatorname {sinc} ^{3}(4)+\cdots ={\frac {1}{2}}.}
シリーズ展開 正規化されていないsinc 関数のテイラー級数は、sin関数のテイラー級数から得ることができます(sin関数も x =0 で1の値をとります)。 罪 x x = ∑ n = 0 ∞ ( − 1 ) n x 2 n ( 2 n + 1 ) ! = 1 − x 2 3 ! + x 4 5 ! − x 6 7 ! + ⋯ {\displaystyle {\frac {\sin x}{x}}=\sum _{n=0}^{\infty }{\frac {(-1)^{n}x^{2n}}{(2n+1)!}}=1-{\frac {x^{2}}{3!}}+{\frac {x^{4}}{5!}}-{\frac {x^{6}}{7!}}+\cdots }
この級数はすべてのx に対して収束する。正規化されたバージョンは容易に導出できる。 罪 π x π x = 1 − π 2 x 2 3 ! + π 4 x 4 5 ! − π 6 x 6 7 ! + ⋯ {\displaystyle {\frac {\sin \pi x}{\pi x}}=1-{\frac {\pi ^{2}x^{2}}{3!}}+{\frac {\pi ^{4}x^{4}}{5!}}-{\frac {\pi ^{6}x^{6}}{7!}}+\cdots }
オイラーは、この級数を バーゼル問題を 解くための無限積形式の展開になぞらえたことで有名である。
高次元 1 次元 sinc 関数の積は、正方形のデカルト格子 (格子 )の多変数sinc 関数を容易に提供します。sinc C ( x , y ) = sinc( x ) sinc( y )で あり 、そのフーリエ変換は周波数空間の正方形 (つまり、2 次元空間で定義されたレンガの壁) の 指示関数 です。非デカルト格子 (たとえば、六角形格子 )の sinc 関数は、その格子のブリルアンゾーン の指示関数 であるフーリエ変換を 持つ関数です。たとえば、六角形格子の sinc 関数は、周波数空間の単位六角形の指示関数 であるフーリエ変換を持つ 関数です。非デカルト格子の場合、この関数は単純なテンソル積 では得られません。しかし、ブリルアンゾーンの幾何学的性質とゾノトープとの関連を利用することで、 六方格子 、体心立方格子 、面心立方格子 、その他の高次元格子のsinc関数の明示的な公式を明示的に導出することができる[ 13 ] 。
例えば、ベクトルの (整数)線形スパンによって 六角格子を 生成できます。u 1 = [ 1 2 3 2 ] そして u 2 = [ 1 2 − 3 2 ] 。 {\displaystyle \mathbf {u} _{1}={\begin{bmatrix}{\frac {1}{2}}\\{\frac {\sqrt {3}}{2}}\end{bmatrix}}\quad {\text{and}}\quad \mathbf {u} _{2}={\begin{bmatrix}{\frac {1}{2}}\\-{\frac {\sqrt {3}}{2}}\end{bmatrix}}.}
示す ξ 1 = 2 3 u 1 、 ξ 2 = 2 3 u 2 、 ξ 3 = − 2 3 ( u 1 + u 2 ) 、 x = [ x y ] 、 {\displaystyle {\boldsymbol {\xi }}_{1}={\tfrac {2}{3}}\mathbf {u} _{1},\quad {\boldsymbol {\xi }}_{2}={\tfrac {2}{3}}\mathbf {u} _{2},\quad {\boldsymbol {\xi }}_{3}=-{\tfrac {2}{3}}(\mathbf {u} _{1}+\mathbf {u} _{2}),\quad \mathbf {x} ={\begin{bmatrix}x\\y\end{bmatrix}},} この六角格子に対するsinc関数は [ 13 ] のよう に導出できる。シン H ( x ) = 1 3 ( コス ( π ξ 1 ⋅ x ) シン ( ξ 2 ⋅ x ) シン ( ξ 3 ⋅ x ) + コス ( π ξ 2 ⋅ x ) シン ( ξ 3 ⋅ x ) シン ( ξ 1 ⋅ x ) + コス ( π ξ 3 ⋅ x ) シン ( ξ 1 ⋅ x ) シン ( ξ 2 ⋅ x ) ) 。 {\displaystyle {\begin{aligned}\operatorname {sinc} _{\text{H}}(\mathbf {x} )={\tfrac {1}{3}}{\big (}&\cos \left(\pi {\boldsymbol {\xi }}_{1}\cdot \mathbf {x} \right)\operatorname {sinc} \left({\boldsymbol {\xi }}_{2}\cdot \mathbf {x} \right)\operatorname {sinc} \left({\boldsymbol {\xi }}_{3}\cdot \mathbf {x} \right)\\&{}+\cos \left(\pi {\boldsymbol {\xi }}_{2}\cdot \mathbf {x} \right)\operatorname {sinc} \left({\boldsymbol {\xi }}_{3}\cdot \mathbf {x} \right)\operatorname {sinc} \left({\boldsymbol {\xi }}_{1}\cdot \mathbf {x} \right)\\&{}+\cos \left(\pi {\boldsymbol {\xi }}_{3}\cdot \mathbf {x} \right)\operatorname {sinc} \left({\boldsymbol {\xi }}_{1}\cdot \mathbf {x} \right)\operatorname {sinc} \left({\boldsymbol {\xi }}_{2}\cdot \mathbf {x} \right){\big )}.\end{aligned}}}
この構成は、一般的な多次元格子に対するランチョス窓 の設計に使用できます。 [ 13 ]
シンチ 一部の著者は類推により双曲線 正弦基数関数 を定義している。[ 14 ] [ 15 ] [ 16 ]
s 私 n h c ( x ) = { シン ( x ) x 、 もし x ≠ 0 1 、 もし x = 0 {\displaystyle \mathrm {sinhc} (x)={\begin{cases}{\displaystyle {\frac {\sinh(x)}{x}},}&{\text{if }}x\neq 0\\{\displaystyle 1,}&{\text{if }}x=0\end{cases}}}
参考文献 ↑ Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. 編 (2010)、「数値計算法」、NIST 数学関数ハンドブック 、ケンブリッジ大学出版局、ISBN 978-0-521-19225-5 MR 2723248 . .↑ Singh, RP; Sapre, SD (2008). Communication Systems, 2E (図解 版). Tata McGraw-Hill Education. p. 15. ISBN 978-0-07-063454-1 。 15ページからの抜粋↑ Weisstein, Eric W. "Sinc関数" . mathworld.wolfram.com . 2023年6月7日 取得 . ↑ Merca, Mircea (2016-03-01). "The cardinal sine function and the Chebyshev–Stirling numbers" . Journal of Number Theory . 160 : 19–31 . doi : 10.1016/j.jnt.2015.08.018 . ISSN 0022-314X . S2CID 124388262 . 1 2 ポイントン、チャールズ A. (2003). デジタルビデオとHDTV . モーガン・カウフマン出版社. p. 147. ISBN 978-1-55860-792-7 。↑ Woodward, PM; Davies, IL (1952年3月)。 「電気通信における情報理論と逆確率」 (PDF) 。Proceedings of the IEE - Part III: Radio and Communication Engineering。99 ( 58): 37–44。doi: 10.1049 / pi - 3.1952.0011 。 ↑ ウッドワード、フィリップ M. (1953). 確率と情報理論、レーダーへの応用 . ロンドン: パーガモン・プレス. p. 29. ISBN 978-0-89006-103-9 . OCLC 488749777 . ↑ オイラー、レオンハルト (1735)。「逆数の級数の和について」。arXiv : math /0506415 。 ↑ Sanjar M. Abrarov; Brendan M. Quine (2015). "sinc 関数の不完全コサイン展開によるサンプリング: Voigt/複素誤差関数への応用" . Appl. Math. Comput . 258 : 425– 435. arXiv : 1407.0533 . doi : 10.1016/j.amc.2015.01.072 . ↑ 「高度な問題6241」。 アメリカ数学月報 。87 (6)。ワシントンDC: アメリカ数学協会 : 496–498 。 1980年6月–7月。doi : 10.1080 /00029890.1980.11995075 。 ↑ Robert Baillie; David Borwein ; Jonathan M. Borwein (2008 年12 月 )。 「 驚くべき Sinc 和と積分」。American Mathematical Monthly。115 ( 10 ) : 888–901。doi : 10.1080 / 00029890.2008.11920606。hdl : 1959.13/940062。JSTOR 27642636。S2CID 496934 。 ↑ Baillie, Robert (2008). "Fun with Fourier series". arXiv : 0806.0150v2 [ math.CA ]. 1 2 3 Ye, W.; Entezari, A. (2012 年 6 月). "多変数 Sinc 関数の幾何学的構成". IEEE Transactions on Image Processing . 21 (6): 2969–2979 . Bibcode : 2012ITIP...21.2969Y . doi : 10.1109/TIP.2011.2162421 . PMID 21775264. S2CID 15313688 . ↑エインズリー 、 マイケル(2010)。 ソナー性能モデリングの原理 。スプリンガー。p. 636。ISBN 9783540876625 。↑ Günter, Peter (2012). Nonlinear Optical Effects and Materials . Springer. p. 258. ISBN 9783540497134 。↑ Schächter, Levi (2013). Beam-Wave Interaction in Periodic and Quasi-Periodic Structures . Springer. p. 241. ISBN 9783662033982 .
Further reading Stenger, Frank (1993). Numerical Methods Based on Sinc and Analytic Functions . Springer Series on Computational Mathematics. Vol. 20. Springer-Verlag New York, Inc. doi :10.1007/978-1-4612-2706-9. ISBN 9781461276371 .