

信号処理において、連続時間信号のパワースペクトルは、その信号を構成する周波数成分へのパワーの分布を表します。 [1]フーリエ解析 によれば、あらゆる物理信号は、いくつかの離散周波数、または連続範囲にわたる周波数スペクトルに分解できます。あらゆる種類の信号(ノイズを含む)の統計的平均をその周波数内容の観点から分析したものを、そのスペクトルと呼びます。
信号のエネルギーが有限の時間間隔に集中している場合、特にその総エネルギーが有限である場合は、エネルギースペクトル密度を計算できます。より一般的に使用されるのは、パワースペクトル密度(PSD、または単にパワースペクトル) です。これは、全時間にわたって存在する信号、または (特に測定の継続時間に関連して) 無限の時間間隔にわたって存在するのと同程度に長い時間間隔にわたって存在する信号に適用されます。PSD は、そのような信号の全時間にわたる総エネルギーが一般に無限であるため、単位時間あたりに見られるスペクトルエネルギー分布を指します。スペクトル成分の合計または積分により、総パワー (物理プロセスの場合) または分散 (統計プロセスの場合) が得られます。これは、パーセバルの定理で規定されているように、時間領域にわたって積分することによって得られるものと同じです。[1]
物理過程のスペクトルには、その性質に関する重要な情報が含まれていることがよくあります。たとえば、楽器のピッチと音色は、スペクトル解析からすぐに判断できます。光源の色は、極めて高い周波数で変動する電磁波の電場のスペクトルによって決まります。このような時系列からスペクトルを取得するには、フーリエ変換と、フーリエ解析に基づく一般化が必要です。多くの場合、時間領域は実際には特に使用されません。たとえば、分光器で光のスペクトルを取得するために分散プリズムが使用される場合や、特定の周波数に敏感な内耳の聴覚受容器への影響を通じて音が知覚される場合などです。
しかし、この記事では、時系列が既知である(少なくとも統計的な意味で)か、直接測定される(コンピュータでサンプリングされたマイクロフォンなど)状況に焦点を当てています。パワースペクトルは、統計信号処理や確率過程の統計的研究、および物理学や工学の多くの他の分野で重要です。通常、プロセスは時間の関数ですが、同様に空間周波数の観点から分解された空間領域のデータを議論することができます。[1]
ユニット
物理学では、信号は電磁波、音波、または機構の振動などの波である可能性があります。信号のパワースペクトル密度(PSD)は、信号に存在する電力を単位周波数あたりの周波数の関数として表します。パワースペクトル密度は通常、ワット/ヘルツ(略してW/Hz)のSI単位で表されます。[2]
たとえば、信号が電圧のみで定義されている場合、指定された振幅に関連付けられた一意の電力はありません。この場合、「電力」は単に信号の2乗で計算されます。これは、その信号が特定のインピーダンスに実際に供給する電力に常に比例するためです。したがって、PSDにはV 2 Hz −1の単位を使用できます。エネルギースペクトル密度(ESD)の単位はV 2 s Hz −1です。これは、エネルギーが電力に時間を掛けた単位(たとえば、ワット時)であるためです。[3]
一般的に、PSD の単位は周波数単位あたりの変動単位の比になります。したがって、たとえば、時間 (秒) に対する一連の変位値 (メートル単位) の PSD の単位は、メートル平方/ヘルツ (m 2 /Hz) になります。ランダム振動の解析では、加速度の PSD にg 2 Hz −1の単位がよく使用されます。ここで、g は重力加速度を表します。[4]
数学的には、信号または独立変数に物理的な次元を割り当てる必要はありません。以下の説明では、x ( t ) の意味は指定しませんが、独立変数は時間であると仮定します。
片面 vs 両面
PSDは、正の周波数のみの片側関数、または正と負の周波数の両方の振幅が半分の両側関数のいずれかになります。ノイズPSDは、一般的に工学では片側関数、物理学では両側関数です。[5]
意味
エネルギースペクトル密度
エネルギースペクトル密度は、信号または時系列のエネルギーが周波数とともにどのように分布しているかを表します。ここで、エネルギーという用語は、信号処理の一般的な意味で使用されています。[6]つまり、信号のエネルギーは次のとおりです。
エネルギースペクトル密度は、有限の総エネルギーを持つ過渡現象、つまりパルス状の信号に最も適しています。有限かどうかにかかわらず、パーセバルの定理(またはプランシュレルの定理)は、信号のエネルギーの別の表現を提供します。[7] ここで、は 周波数(Hz )での のフーリエ変換 の値です。この定理は離散時間の場合にも当てはまります。左辺の積分は信号のエネルギーであるため、 の値は、周波数間隔 の周波数での信号に含まれるエネルギーを表す、無限小周波数間隔で乗じられた密度関数として解釈できます。
したがって、のエネルギースペクトル密度は次のように定義される: [8]
関数と の自己相関はフーリエ変換のペアを形成し、その結果はウィーナー・ヒンチンの定理としても知られています(ピリオドグラムも参照)。
信号のエネルギー スペクトル密度を測定する物理的な例として、 がインピーダンスの伝送線路に沿って伝播する電気パルスの電位(ボルト 単位)を表し、伝送線路が整合抵抗器で終端されている (つまり、パルス エネルギーはすべて抵抗器に送られ、反射して戻ってくるものがない) とします。オームの法則により、時間 に抵抗器に送られる電力は に等しいため、総エネルギーはパルスの持続時間にわたって時間について積分することで求められます。周波数 におけるエネルギー スペクトル密度の値を求めるには、伝送線路と抵抗器の間に、対象周波数付近の狭い範囲の周波数 (たとえば ) のみを通過させるバンドパス フィルタを挿入し、抵抗器で消費される総エネルギーを測定します。この場合、 におけるエネルギー スペクトル密度の値は と推定されます。この例では、電力の単位は V 2 Ω −1、エネルギーの単位は V 2 s Ω −1 = Jであり、したがってエネルギー スペクトル密度の推定値は、必要に応じて J Hz −1の単位になります。多くの場合、除算の手順を忘れて、エネルギー スペクトル密度の単位が V 2 Hz −1になることがよくあります。
この定義は、離散時間にサンプリングされた信号など、可算無限個の値を持つ離散信号に簡単に一般化されます。 ここで、はの離散時間フーリエ変換です。 サンプリング間隔は、正しい物理単位を維持し、極限で連続ケースを回復するために必要です。しかし、数学科学では、間隔は 1 に設定されることが多く、一般性を犠牲にして結果を単純化します。(正規化された周波数 も参照してください)
パワースペクトル密度

上記のエネルギー スペクトル密度の定義は、エネルギーが 1 つの時間ウィンドウに集中している過渡現象 (パルスのような信号) に適しています。この場合、信号のフーリエ変換が一般的に存在します。全時間にわたる連続信号の場合は、むしろ定常プロセス用に存在するパワー スペクトル密度(PSD)を定義する必要があります。これは、前述の簡単な例のように、信号または時系列のパワーが周波数にわたってどのように分布するかを説明します。ここで、パワーは実際の物理的パワーである場合もあれば、抽象的な信号での便宜上、信号の 2 乗値で単純に識別される場合もあります。たとえば、統計学者は、関数の時間的 (または別の独立変数) な変化を研究し、電気信号 (その他の物理的プロセスの中でも) との類似性を使用して、物理的パワーが関与していない場合でも、これをパワー スペクトルと呼ぶのが通例です。物理的な電圧源を作成し、それを 1オームの抵抗器の端子に適用すると、その抵抗器で消費される瞬間的な電力はワットで表されます。
したがって、信号の全時間にわたる平均電力は、任意の時間を中心にした周期を持つ次の時間平均で与えられます。
しかし、後続の計算を扱うには、積分の範囲で時間制限を扱うよりも、信号自体で時間制限を扱う方が便利です。そのため、平均電力の別の表現があり、ここで、およびは任意の期間内では 1 で、それ以外の場合は 0 です。 明らかに、上記の P の式が 0 でない場合は、T が無限に増加するのと同様に、積分は無限に増加する必要があります。これが、そのような場合に、発散積分である 信号のエネルギーを使用できない理由です。
信号の周波数成分を解析する際に、通常のフーリエ変換を計算したい場合があります。しかし、多くの関心信号に対して、フーリエ変換は正式には存在しません。[nb 1]とにかく、パーセバルの定理によれば、平均電力は次のように書き直すことができます。
すると、パワースペクトル密度は単純に上記の積分関数として定義される。[9] [10]
ここから、畳み込み定理により、を との時間畳み込みのフーリエ変換と見ることもできます。ここで、* は複素共役を表します。 と を作成することを考慮すると 、次が得られます。 ここで、畳み込み定理は3 行目から 4 行目に移るときに使用されています。
ここで、上記の時間畳み込みを周期で割り、極限を とすると、それはウィンドウ処理されていない信号 の自己相関関数となり、 がエルゴードである限り と表されます。これはほとんどの場合に当てはまりますが、すべての場合に当てはまるわけではありません。[nb 2]
ここから、再びのエルゴード性を仮定すると、パワースペクトル密度は自己相関関数のフーリエ変換として見つけられることがわかります(ウィーナー・ヒンチンの定理)。[11]
多くの著者は、この等式を実際にパワースペクトル密度を定義するために使用しています。 [12]
与えられた周波数帯域 (ただし )における信号のパワーは、周波数にわたって積分することで計算できます。 であるため、正の周波数帯域と負の周波数帯域に等しい量のパワーを帰属させることができ、次の形式で係数 2 を説明します(このような些細な係数は、使用する規則によって異なります)。 より一般的には、同様の手法を使用して、時間によって変化するスペクトル密度を推定できます。 この場合、時間間隔は無限大に近づくのではなく有限です。 未満の周波数はサンプリングされず、 の整数倍でない周波数の結果は独立していないため、スペクトル範囲と解像度が低下します。 このような時系列を 1 つだけ使用すると、推定されたパワー スペクトルは非常に「ノイズが多い」ものになります。 ただし、指定された時間ウィンドウで評価されたの実現の統計的アンサンブルに対応する多数の(または無限の)短期スペクトルを使用して、(上記の式で)期待値を評価できる場合は、このノイズを軽減できます。
エネルギースペクトル密度と同様に、パワースペクトル密度の定義は離散時間変数に一般化できます。前と同様に、離散時間にサンプリングされた信号と全測定期間のウィンドウを考察できます。PSD の単一の推定値は有限数のサンプリングによって得られることに注意してください。前と同様に、実際の PSD は(したがって) が無限大に近づき、期待値が正式に適用されるときに得られます。実際のアプリケーションでは、通常、有限測定 PSD を多数の試行にわたって平均化して、個々の測定の基礎となる物理プロセスの理論的な PSD のより正確な推定値を取得します。この計算された PSD は、ピリオドグラムと呼ばれることもあります。このピリオドグラムは、推定値の数と平均化時間間隔が無限大に近づくにつれて、真の PSD に収束します。[13]
2 つの信号の両方がパワースペクトル密度を持つ場合、クロススペクトル密度も同様に計算できます。PSD は自己相関に関連し、クロススペクトル密度は相互相関に関連します。
パワースペクトル密度の特性
PSDの特性としては次のようなものがある: [14]
- パワースペクトルは常に実数かつ非負であり、実数値プロセスのスペクトルも周波数の偶関数です。
- 連続確率過程x(t)に対して、自己相関関数Rxx ( t )は逆フーリエ変換を用いてそのパワースペクトルSxx(f)から再構成することができる。
- パーセバルの定理を使用すると、すべての周波数にわたってパワースペクトルを積分することで、プロセスの 分散(平均パワー)を計算できます。
- パワースペクトル密度を持つ実プロセスx ( t ) に対して、DCからfまでの周波数に含まれる平均帯域制限パワーを指定する積分スペクトルまたはパワースペクトル分布を次のように計算することができます。 [15]前述の全パワー(信号分散)の式は、 f → ∞ の特殊なケースであることに注意してください 。
クロスパワースペクトル密度
2 つの信号と があり、それぞれがパワー スペクトル密度とを持つ場合、クロス パワー スペクトル密度( CPSD ) またはクロス スペクトル密度( CSD )を定義することができます。まず、このような結合信号の平均パワーを考えてみましょう。
パワースペクトル密度の導出に使用したのと同じ表記法と方法を使用して、パーセバルの定理を利用して、次式を得ます 。 ここでも、 と の寄与はすでにわかっています。 であるため、クロスパワーへの完全な寄与は、一般に、個別のCPSDの実部の 2 倍から生じます。 前と同様に、ここからこれらの積を時間畳み込みのフーリエ変換として書き直します。これを周期で割って極限まで取ると、相互相関関数のフーリエ変換になります。[16] ここで、 はとの相互相関であり、は との相互相関です。 これを考慮すると、PSD は に対する CSD の特殊なケースであることがわかります。と が実信号 (電圧や電流など) の場合、それらのフーリエ変換と は通常、慣例により正の周波数に制限されます。したがって、一般的な信号処理では、完全なCPSD は、係数 2 でスケーリングされた CPSDの 1 つにすぎません。
離散信号x nとy nの場合、相互スペクトル密度と相互共分散の関係は
推定
スペクトル密度推定の目的は、一連の時間サンプルからランダム信号のスペクトル密度を推定することです。信号について何がわかっているかに応じて、推定手法にはパラメトリックまたはノンパラメトリックなアプローチが含まれ、時間領域または周波数領域分析に基づく場合があります。たとえば、一般的なパラメトリック手法では、観測値を自己回帰モデルに当てはめます。一般的なノンパラメトリック手法は、ピリオドグラムです。
スペクトル密度は通常、フーリエ変換法(ウェルチ法など)を使用して推定されますが、最大エントロピー法などの他の手法も使用できます。
関連概念
- 信号のスペクトル重心は、スペクトル密度関数の中点、つまり分布を 2 つの等しい部分に分割する周波数です。
- スペクトルエッジ周波数(SEF)は、通常「SEF x」と表され、与えられた信号の全電力のxパーセントがその下に位置する周波数を表します。通常、 xは75〜95の範囲です。これは特にEEGモニタリングで使用される一般的な指標であり、その場合SEFは麻酔の深さや睡眠段階を推定するためにさまざまな方法で使用されています。[17] [18]
- スペクトルエンベロープは、スペクトル密度のエンベロープ曲線です。これは、ある時点 (正確には、1 つのウィンドウ) を表します。たとえば、分光計を使用したリモート センシングでは、フィーチャのスペクトル エンベロープは、対象の各スペクトル バンドの明るさレベルの範囲によって定義される、そのスペクトル特性の境界です。
- スペクトル密度は、時間の関数ではなく、周波数の関数です。ただし、長い信号の小さなウィンドウのスペクトル密度を計算し、ウィンドウに関連付けられた時間に対してプロットすることができます。このようなグラフはスペクトログラムと呼ばれます。これは、短時間フーリエ変換やウェーブレットなどの多くのスペクトル分析手法の基礎となります。
- 「スペクトル」とは、一般的には、上で説明したように、周波数に対する信号内容の分布を表すパワー スペクトル密度を意味します。伝達関数(ボード線図、チャープなど) の場合、完全な周波数応答は、パワー対周波数と位相対周波数の 2 つの部分でグラフ化できます。つまり、位相スペクトル密度、位相スペクトル、またはスペクトル位相です。あまり一般的ではありませんが、2 つの部分は伝達関数の実部と虚部である場合があります。これは、周波数の関数として位相 (または同等の実部と虚部) も含む伝達関数の周波数応答と混同しないでください。時間領域のインパルス応答は、一般に、位相部分なしでパワー スペクトル密度のみから一意に復元することはできません。これらもフーリエ変換のペアですが、対称性 ( 自己相関の場合のように) がないため、フーリエ変換は実数値になります。超短パルス#スペクトル位相、位相ノイズ、群遅延を参照してください。
- 振幅スペクトル密度(ASD )に遭遇することがあります。これはPSDの平方根です。電圧信号のASDの単位はV Hz −1/2です。[19]これはスペクトルの形状がかなり一定である場合に役立ちます。なぜならASDの変動は信号の電圧レベル自体の変動に比例するからです。しかし、数学的にはPSDを使用するのが好ましいです。なぜならその場合にのみ、曲線の下の面積が全周波数または指定された帯域幅での実際の電力に関して意味があるからです。
アプリケーション
時間とともに変化する変数として表せる信号には、対応する周波数スペクトルがあります。これには、可視光(色として認識される)、音符 (ピッチとして認識される)、ラジオ/テレビ(周波数、または場合によっては波長 で指定される)、さらには地球の定期的な自転など、よく知られたエンティティが含まれます。これらの信号を周波数スペクトルの形で表示すると、受信信号の特定の側面や信号を生成する基礎となるプロセスが明らかになります。場合によっては、周波数スペクトルに正弦波成分に対応する明確なピークが含まれることがあります。また、基本ピークの高調波に対応するピークが存在する場合もあり、これは単純な正弦波ではない周期信号を示しています。または、連続スペクトルが、共鳴に対応して強く強調された狭い周波数間隔、またはノッチ フィルターによって生成されるようなほぼゼロの電力を含む周波数間隔を示すことがあります。
電気工学

信号のパワースペクトルの概念と使用は、電気工学、特に無線通信、レーダー、関連システム、および受動リモートセンシング技術を含む電子通信システムにおいて基本的なものです。信号のパワースペクトルを観察および測定するには、スペクトルアナライザーと呼ばれる電子機器が使用されます。
スペクトル アナライザは、入力信号の短時間フーリエ変換(STFT)の振幅を測定します。分析対象の信号が定常プロセスであると考えられる場合、STFT は、そのパワー スペクトル密度の平滑化された推定値として適しています。
宇宙論
原始的変動、つまり初期宇宙における密度の変動は、空間スケールの関数として変動のパワーを示すパワースペクトルによって定量化されます。
参照
注記
- ^ 例えば、(Risken & Frank 1996、p. 30) などの著者は、正規化されていないフーリエ変換を正式な方法で使用して、パワースペクトル密度の定義を定式化しています 。 ここで、はディラックのデルタ関数です。このような正式な記述は、直感を導くのに役立つ場合がありますが、常に細心の注意を払って使用する必要があります。
- ^ ウィーナー・ヒンチンの定理は、より弱い仮説の下では、任意の広義の定常過程に対してこの式の意味を成す。 は絶対積分可能である必要はなく、存在するだけでよい。しかし、積分は通常どおりには解釈できない。この式は、関数の代わりに分布(統計的累積分布関数の意味ではなく、ローラン・シュワルツの意味で)を含むと解釈しても意味を成す。 が連続である場合、ボホナーの定理を使用して、そのフーリエ変換が分布関数が F である(ただし、必ずしも関数としてではなく、必ずしも確率密度を持つとは限らない)正の測度 として存在することを証明できる。
- ^ abc P Stoica & R Moses (2005). 「信号のスペクトル分析」(PDF)。
- ^ マラル 2004年。
- ^ ノートン&カルチュブ 2003.
- ^ ビロリーニ 2007、83ページ。
- ^ Paschotta, Rüdiger. 「パワースペクトル密度」. rp-photonics.com . 2024年4月15日時点のオリジナルよりアーカイブ。2024年6月26日閲覧。
- ^ オッペンハイム&ヴェルゲーズ 2016年、12ページ。
- ^ スタイン2000、108、115頁。
- ^ オッペンハイム&ヴェルゲーズ 2016年、14ページ。
- ^ オッペンハイム&ヴェルゲーズ 2016年、422-423頁。
- ^ ミラー&チルダーズ 2012年、429-431頁。
- ^ ミラー&チルダーズ 2012年、433ページ。
- ^ Dennis Ward Ricker (2003). エコー信号処理. Springer. ISBN 978-1-4020-7395-3。
- ^ ブラウン&ファン 1997.
- ^ ミラー&チルダーズ 2012年、431ページ。
- ^ ダベンポート&ルート 1987年。
- ^ William D Penny (2009). 「信号処理コース、第 7 章」
- ^ イランマネシュ & ロドリゲス・ビジェガス 2017.
- ^ イムティアス & ロドリゲス・ビジェガス 2014.
- ^ Michael Cerna & Audrey F. Harvey (2000). 「FFT ベースの信号解析と測定の基礎」(PDF)。
参考文献
- ビロリーニ、アレッサンドロ (2007)。信頼性工学。ベルリン、ニューヨーク:シュプリンガーサイエンス&ビジネスメディア。ISBN 978-3-540-49388-4。
- Brown, Robert Grover; Hwang, Patrick YC (1997)。ランダム信号と応用カルマンフィルタリング入門、MATLAB演習と解答付き。ニューヨーク: Wiley- Liss。ISBN 978-0-471-12839-7。
- ダベンポート、ウィルバー B. (Jr); ルート、ウィリアム L. (1987)。ランダム信号とノイズの理論入門。ニューヨーク: Wiley-IEEE 出版。ISBN 978-0-87942-235-6。
- Imtiaz, Syed Anas; Rodriguez-Villegas, Esther (2014)。「シングル チャネル EEG を使用した REM 睡眠検出のための低計算コスト アルゴリズム」。Annals of Biomedical Engineering。42 ( 11 ): 2344–59。doi : 10.1007 /s10439-014-1085-6。PMC 4204008。PMID 25113231。
- Iranmanesh, Saam; Rodriguez-Villegas, Esther (2017). 「チップ上の超低電力睡眠スピンドル検出システム」。IEEE Transactions on Biomedical Circuits and Systems . 11 (4): 858–866. doi :10.1109/TBCAS.2017.2690908. hdl : 10044/1/46059 . PMID 28541914. S2CID 206608057.
- Maral, Gerard (2004)。VSATネットワーク。ウェスト・サセックス、イングランド、ホーボーケン、ニュージャージー:Wiley。ISBN 978-0-470-86684-9。
- ミラー、スコット、チルダーズ、ドナルド (2012)。確率とランダムプロセス。ボストン、マサチューセッツ州: アカデミック プレス。ISBN 978-0-12-386981-4. OCLC 696092052.
- Norton, MP; Karczub, DG (2003)。エンジニアのための騒音・振動解析の基礎。ケンブリッジ: ケンブリッジ大学出版局。ISBN 978-0-521-49913-2。
- オッペンハイム、アラン V.; ヴェルゲーズ、ジョージ C. (2016)。信号、システム、推論。ボストン: ピアソン。ISBN 978-0-13-394328-3。
- リスクン、ハンネス、フランク、ティル (1996)。フォッカー・プランク方程式。ニューヨーク: シュプリンガー・サイエンス & ビジネス・メディア。ISBN 978-3-540-61530-9。
- スタイン、ジョナサン Y. (2000)。デジタル信号処理。ニューヨーク・ワインハイム:ワイリー・インターサイエンス。ISBN 978-0-471-29546-4。
外部リンク
- パワースペクトル密度 Matlab スクリプト
