数学 および信号処理 において、解析信号とは 負の周波数 成分を持たない複素数値関数 である。[ 1 ] 解析信号の実部と虚部は、ヒルベルト変換 によって互いに関連付けられた実数値関数である。
実数値 関数の解析的表現は 、元の関数とそのヒルベルト変換からなる解析信号です。この表現は、多くの数学的操作を容易にします。基本的な考え方は、 実数値関数の フーリエ変換 (またはスペクトル )の負の周波数成分は、スペクトルのエルミート対称性 により不要であるということです。これらの負の周波数成分は、代わりに複素数値関数を扱うことを厭わなければ、情報損失なく破棄できます。これにより、関数の特定の属性がより容易に把握できるようになり、シングルサイドバンドなどの変調および復調技術の導出が容易になります。
操作対象の関数に負の周波数成分がない限り(つまり、解析的 である限り)、複素数から実数への変換は虚数部を捨てるだけで済みます。解析的表現はフェーザ 概念の一般化です。[ 2 ] フェーザは時間不変の振幅、位相、周波数に制限されますが、解析信号は時間可変のパラメータを許容します。
意味 解析信号を生成するための伝達関数 もしs ( t ) {\displaystyle s(t)} はフーリエ変換を持つ実数値 関数ですS ( f ) {\displaystyle S(f)} (どこf {\displaystyle f} は周波数を表す実数値である)ならば、変換はに関してエルミート対称性を持つ。 f = 0 {\displaystyle f=0} 軸:
S ( − f ) = S ( f ) * 、 {\displaystyle S(-f)=S(f)^{*},} どこS ( f ) * {\displaystyle S(f)^{*}} は、 S ( f ) {\displaystyle S(f)} 機能:
S 1 ( f ) ≜ { 2 S ( f ) 、 のために f > 0 、 S ( f ) 、 のために f = 0 、 0 、 のために f < 0 = 2 u ( f ) ⏟ 1 + サイン ( f ) S ( f ) = S ( f ) + サイン ( f ) S ( f ) 、 {\displaystyle {\begin{aligned}S_{\mathrm {a} }(f)&\triangleq {\begin{cases}2S(f),&{\text{for}}\ f>0,\\S(f),&{\text{for}}\ f=0,\\0,&{\text{for}}\ f<0\end{cases}}\\&=\underbrace {2\operatorname {u} (f)} _{1+\operatorname {sgn}(f)}S(f)=S(f)+\operatorname {sgn}(f)S(f),\end{aligned}}} どこ
u ( f ) {\displaystyle \operatorname {u} (f)} はヘヴィサイド階段関数 です。サイン ( f ) {\displaystyle \operatorname {sgn}(f)} は符号関数 です。非負の周波数 成分のみを含むS ( f ) {\displaystyle S(f)} 。そして、エルミート対称性により、この操作は可逆である。S ( f ) {\displaystyle S(f)} :
S ( f ) = { 1 2 S 1 ( f ) 、 のために f > 0 、 S 1 ( f ) 、 のために f = 0 、 1 2 S 1 ( − f ) * 、 のために f < 0 (エルミート対称性) = 1 2 [ S 1 ( f ) + S 1 ( − f ) * ] 。 {\displaystyle {\begin{aligned}S(f)&={\begin{cases}{\frac {1}{2}}S_{\mathrm {a} }(f),&{\text{for}}\ f>0,\\S_{\mathrm {a} }(f),&{\text{for}}\ f=0,\\{\frac {1}{2}}S_{\mathrm {a} }(-f)^{*},&{\text{for}}\ f<0\ {\text{(エルミート対称性)}}\end{cases}}\\&={\frac {1}{2}}[S_{\mathrm {a} }(f)+S_{\mathrm {a} }(-f)^{*}].\end{aligned}}} 分析信号 s ( t ) {\displaystyle s(t)} は、の逆フーリエ変換である。S 1 ( f ) {\displaystyle S_{\mathrm {a} }(f)} :
s 1 ( t ) ≜ F − 1 [ S 1 ( f ) ] = F − 1 [ S ( f ) + サイン ( f ) ⋅ S ( f ) ] = F − 1 { S ( f ) } ⏟ s ( t ) + F − 1 { サイン ( f ) } ⏟ j 1 π t * F − 1 { S ( f ) } ⏟ s ( t ) ⏞ 畳み込み = s ( t ) + j [ 1 π t * s ( t ) ] ⏟ H [ s ( t ) ] = s ( t ) + j s ^ ( t ) 、 {\displaystyle {\begin{aligned}s_{\mathrm {a} }(t)&\triangleq {\mathcal {F}}^{-1}[S_{\mathrm {a} }(f)]\\&={\mathcal {F}}^{-1}[S(f)+\operatorname {sgn}(f)\cdot S(f)]\\&=\underbrace {{\mathcal {F}}^{-1}\{S(f)\}} _{s(t)}+\overbrace {\underbrace {{\mathcal {F}}^{-1}\{\operatorname {sgn}(f)\}} _{j{\frac {1}{\pi t}}}*\underbrace {{\mathcal {F}}^{-1}\{S(f)\}} _{s(t)}} ^{\text{畳み込み}}\\&=s(t)+j\underbrace {\left[{1 \over \pi t}*s(t)\right]} _{\operatorname {\mathcal {H}} [s(t)]}\\&=s(t)+j{\hat {s}}(t),\end{aligned}}} どこ
s ^ ( t ) ≜ H [ s ( t ) ] {\displaystyle {\hat {s}}(t)\triangleq \operatorname {\mathcal {H}} [s(t)]} はヒルベルト変換 であるs ( t ) {\displaystyle s(t)} ;* {\displaystyle *} はバイナリ畳み込み 演算子です。j {\displaystyle j} は虚数単位 です。留意すべき点として、s ( t ) = s ( t ) * δ ( t ) 、 {\displaystyle s(t)=s(t)*\delta (t),} これは、負の周波数成分を直接除去するフィルタリング操作として表現することもできます。
s 1 ( t ) = s ( t ) * [ δ ( t ) + j 1 π t ] ⏟ F − 1 { 2 u ( f ) } 。 {\displaystyle s_{\mathrm {a} }(t)=s(t)*\underbrace {\left[\delta (t)+j{1 \over \pi t}\right]} _{{\mathcal {F}}^{-1}\{2u(f)\}}.}
例
例1 s ( t ) = コス ( ω t ) 、 {\displaystyle s(t)=\cos(\omega t),} どこ ω > 0. {\displaystyle \omega >0.} それから:
s ^ ( t ) = コス ( ω t − π 2 ) = 罪 ( ω t ) 、 s 1 ( t ) = s ( t ) + j s ^ ( t ) = コス ( ω t ) + j 罪 ( ω t ) = e j ω t 。 {\displaystyle {\begin{aligned}{\hat {s}}(t)&=\cos \left(\omega t-{\frac {\pi }{2}}\right)=\sin(\omega t),\\s_{\mathrm {a} }(t)&=s(t)+j{\hat {s}}(t)=\cos(\omega t)+j\sin(\omega t)=e^{j\omega t}.\end{aligned}}} 最後の等式はオイラーの公式 であり、その系 はコス ( ω t ) = 1 2 ( e j ω t + e j ( − ω ) t ) 。 {\textstyle \cos(\omega t)={\frac {1}{2}}\left(e^{j\omega t}+e^{j(-\omega )t}\right).} 一般に、単純な正弦波の解析的表現は、複素指数関数を用いて表し、負の周波数 成分を破棄し、正の周波数成分を2倍することによって得られます。また、複数の正弦波の和の解析的表現は、個々の正弦波の解析的表現の和となります。
例2 ここでは、オイラーの公式を用いて負の周波数を特定し、除外します。
s ( t ) = コス ( ω t + θ ) = 1 2 ( e j ( ω t + θ ) + e − j ( ω t + θ ) ) {\displaystyle s(t)=\cos(\omega t+\theta )={\frac {1}{2}}\left(e^{j(\omega t+\theta )}+e^{-j(\omega t+\theta )}\right)} それから:
s 1 ( t ) = { e j ( ω t + θ ) = e j | ω | t ⋅ e j θ 、 もし ω > 0 、 e − j ( ω t + θ ) = e j | ω | t ⋅ e − j θ 、 もし ω < 0. {\displaystyle s_{\mathrm {a} }(t)={\begin{cases}e^{j(\omega t+\theta )}\ \ =\ e^{j|\omega |t}\cdot e^{j\theta },&{\text{if}}\ \omega >0,\\e^{-j(\omega t+\theta )}=\ e^{j|\omega |t}\cdot e^{-j\theta },&{\text{if}}\ \omega <0.\end{cases}}}
物件
瞬時振幅と位相 青色で示された関数と、赤色で示されたその解析的表現の大きさ。包絡線効果を示している。 解析信号は極座標 でも表現できる。
s 1 ( t ) = s m ( t ) e j ϕ ( t ) 、 {\displaystyle s_{\mathrm {a} }(t)=s_{\mathrm {m} }(t)e^{j\phi (t)},} ここで、以下の時間変動量が導入される。
s m ( t ) ≜ | s 1 ( t ) | {\displaystyle s_{\mathrm {m} }(t)\triangleq |s_{\mathrm {a} }(t)|} これは瞬間振幅 またはエンベロープ と呼ばれます。ϕ ( t ) ≜ 引数 [ s 1 ( t ) ] {\displaystyle \phi (t)\triangleq \arg \!\left[s_{\mathrm {a} }(t)\right]} これは瞬時位相 または位相角 と呼ばれます。添付の図では、青い曲線はs ( t ) {\displaystyle s(t)} そして赤い曲線は対応するs m ( t ) {\displaystyle s_{\mathrm {m} }(t)} 。
展開された瞬時 位相 の時間微分はラジアン/秒 の単位を持ち、瞬時角周波数 と呼ばれます。
ω ( t ) ≜ d ϕ d t ( t ) 。 {\displaystyle \omega (t)\triangleq {\frac {d\phi }{dt}}(t).} したがって、瞬時周波数 (ヘルツ単位 )は次のようになります。
f ( t ) ≜ 1 2 π ω ( t ) 。 {\displaystyle f(t)\triangleq {\frac {1}{2\pi }}\omega (t).} [ 3 ] 瞬時振幅、瞬時位相、瞬時周波数は、信号の局所的な特徴を測定および検出するために、いくつかの用途で使用されます。信号の解析的表現のもう1つの用途は、変調信号 の復調に関するものです。極座標は、振幅変調 と位相(または周波数)変調の影響を都合よく分離し、特定の種類の信号を効果的に復調します。
解析信号を多変数信号に拡張する 解析信号の概念は、通常は時間である単一変数の信号については明確に定義されています。2つ以上の変数の信号の場合、解析信号はさまざまな方法で定義できますが、以下に2つのアプローチを示します。
単一遺伝子シグナル 解析信号の実部と虚部は、1変数信号に対して定義されるベクトル値単一遺伝子信号 の2つの要素に対応します。ただし、単一遺伝子信号は任意の数の変数に容易に拡張でき、n 変数信号の場合には( n +1) 次元のベクトル値関数 が生成されます。
注記 ↑ 「複素包絡線(または複素振幅)」 [ 6 ] ↑ 「複素包絡線(または複素振幅)」、p. 586 [ 7 ] ↑ 「複素エンベロープは、時間の関数としての複素振幅の拡張解釈である。」p. 85 [ 8 ]
参考文献 ↑ Smith, JO「解析信号とヒルベルト変換フィルタ」、『離散フーリエ変換(DFT)の数学とオーディオアプリケーション』第2版、 https://ccrma.stanford.edu/~jos/r320/Analytic_Signals_Hilbert_Transform.html、またはhttps://www.dsprelated.com/freebooks/mdft/Analytic_Signals_Hilbert_Transform.html、オンライン書籍、2007年版、2021年4月29日アクセス。 1 2 ブレイスウェル、ロン。『フーリエ変換とその応用 』マグロウヒル、2000年、361-362頁 ↑ B. Boashash、「信号の瞬時周波数の推定と解釈 - パートI:基礎」、Proceedings of the IEEE、第80巻、第4号、519~538ページ、1992年4月 ↑ Justice, J. (1979-12-01). "音楽計算における解析的信号処理". IEEE Transactions on Acoustics, Speech, and Signal Processing . 27 (6): 670– 684. doi : 10.1109/TASSP.1979.1163321 . ISSN 0096-3518 . ↑ B. Boashash、「時間周波数信号解析におけるウィグナー分布の使用に関する注記」、IEEE Trans. on Acoustics, Speech, and Signal Processing、第26巻、第9号、1987年 ↑ Hlawatsch, Franz; Auger, François (2013-03-01). 時間周波数解析 . John Wiley & Sons. ISBN 9781118623831 。↑ ドリガーズ、ロナルド・G. (2003年1月1日). 光学工学百科事典:アベラス、1-1024ページ . CRC Press. ISBN 9780824742508 。↑ 岡本健一(2001年1月1日) 『地球環境リモートセンシング 』IOS Press. ISBN 9781586031015 。
さらに読む レオン・コーエン著、『時間周波数解析』 、プレンティス・ホール、アッパー・サドル・リバー、1995年。 フレデリック・W・キング著『ヒルベルト変換』 第2巻、ケンブリッジ大学出版局、ケンブリッジ、2009年。 B. Boashash、『時間周波数信号解析と処理:包括的な参考文献』 、Elsevier Science、オックスフォード、2003年。