JPH0526919A - スペクトラム・アナライザ - Google Patents
スペクトラム・アナライザInfo
- Publication number
- JPH0526919A JPH0526919A JP17676591A JP17676591A JPH0526919A JP H0526919 A JPH0526919 A JP H0526919A JP 17676591 A JP17676591 A JP 17676591A JP 17676591 A JP17676591 A JP 17676591A JP H0526919 A JPH0526919 A JP H0526919A
- Authority
- JP
- Japan
- Prior art keywords
- fast fourier
- fourier transform
- analog
- digital
- digital data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Abstract
(57)【要約】 (修正有)
【目的】高速フーリエ変換に固有のスカロップ・ロスを
補正して測定レベル確度を向上したスペクトラム・アナ
ライザを提供する。 【構成】入力信号をアナログ−ディジタル変換するアナ
ログ−ディジタル変換器9を具備し、アナログ−ディジ
タル変換されたディジタル・データを一時的に記憶する
メモリ11を具備し、上記メモリ11に記憶されたディ
ジタル・データに対し窓関数としてガウス関数を掛ける
掛算器を具備し、ディジタル・データに対して窓関数を
掛けた結果のものに対して高速フーリエ変換処理を施す
高速フーリエ変換器を具備し、高速フーリエ変換処理結
果に対して3次関数近似による補間処理を施す補間処理
器を具備し、補間処理結果を表示する表示器13を具備
したスペクトラム・アナライザ。
補正して測定レベル確度を向上したスペクトラム・アナ
ライザを提供する。 【構成】入力信号をアナログ−ディジタル変換するアナ
ログ−ディジタル変換器9を具備し、アナログ−ディジ
タル変換されたディジタル・データを一時的に記憶する
メモリ11を具備し、上記メモリ11に記憶されたディ
ジタル・データに対し窓関数としてガウス関数を掛ける
掛算器を具備し、ディジタル・データに対して窓関数を
掛けた結果のものに対して高速フーリエ変換処理を施す
高速フーリエ変換器を具備し、高速フーリエ変換処理結
果に対して3次関数近似による補間処理を施す補間処理
器を具備し、補間処理結果を表示する表示器13を具備
したスペクトラム・アナライザ。
Description
【0001】
【産業上の利用分野】この発明は、スペクトラム・アナ
ライザに関し、特に高速フーリエ変換(FFT)に固有
のスカロップ・ロス(Scallop Loss)を補
正して測定レベル確度を向上したスペクトラム・アナラ
イザに関する。
ライザに関し、特に高速フーリエ変換(FFT)に固有
のスカロップ・ロス(Scallop Loss)を補
正して測定レベル確度を向上したスペクトラム・アナラ
イザに関する。
【0002】
【従来の技術】この発明の従来例を図1を参照して説明
する。図1において、分析されるべき信号fS1は入力端
子1に入力され、第1の低域フィルタ2を介して第1の
ミキサ3の一方の入力に供給される。4は第1のミキサ
3の他方の入力に接続された第1の局部発振器であり、
周波数fL01の信号を出力するものである。ミキサ3の
周波数fS2(=fL02 ±fS1)の出力信号は中心周波数
fIF1の帯域フィルタタ5を介して第2のミキサ6の一
方の入力端子に印加されている。7は第2のミサ6の他
方の入力に接続された第2の局部発振器であり、周波数
fLO2の信号を出力するものである。ミキサ6の周波数
fS3(=fLO2±fS2)の出力信号は第2の低域フィル
タ8を介してA/D変換器9に印加されている。10は
サンプリング周波数(fSmP)発振器である。A/D変
換器9においてサンプリング周波数(fSmP)でA/D
変換されたディジタルIF信号fIFdは、メモリ11に
一時的に記憶され、次いで内臓コンピュータ処理部12
においてFFT処理されて周波数領域データに変換後、
表示器13に表示される。なお、FFTにおいては、一
般に、信号にFFT処理を施すに先だってこれに対して
窓関数が掛けられる。
する。図1において、分析されるべき信号fS1は入力端
子1に入力され、第1の低域フィルタ2を介して第1の
ミキサ3の一方の入力に供給される。4は第1のミキサ
3の他方の入力に接続された第1の局部発振器であり、
周波数fL01の信号を出力するものである。ミキサ3の
周波数fS2(=fL02 ±fS1)の出力信号は中心周波数
fIF1の帯域フィルタタ5を介して第2のミキサ6の一
方の入力端子に印加されている。7は第2のミサ6の他
方の入力に接続された第2の局部発振器であり、周波数
fLO2の信号を出力するものである。ミキサ6の周波数
fS3(=fLO2±fS2)の出力信号は第2の低域フィル
タ8を介してA/D変換器9に印加されている。10は
サンプリング周波数(fSmP)発振器である。A/D変
換器9においてサンプリング周波数(fSmP)でA/D
変換されたディジタルIF信号fIFdは、メモリ11に
一時的に記憶され、次いで内臓コンピュータ処理部12
においてFFT処理されて周波数領域データに変換後、
表示器13に表示される。なお、FFTにおいては、一
般に、信号にFFT処理を施すに先だってこれに対して
窓関数が掛けられる。
【0003】
【発明が解決しようとする課題】上述の如くに高速フー
リエ変換された周波数領域データには、スカロップ・ロ
スと称するレベル変動(低下)が必然的に付随する。こ
こで、スカロップ・ロスについて図2を参照して簡単に
説明する。高速フーリエ変換においては、信号に窓関数
を掛けることにより周波数スペクトルの漏れが少なくな
るようにしている。ところが、この場合、周波数は量子
化されているので、観測される信号の周波数foが量子
化された周波数とたまたま共通している場合に限って正
しいレベルが観測されるに過ぎず、量子化された周波数
の間の周波数における信号は如何に補間しようとも必ず
真のレベルより低いレベルとして観測されることとな
る。周波数foのところのレベルと量子化された周波数
1単位の1/2だけずれたところのレベルとの間の差を
スカロープ・ロスと称している。このロスは概ね1dB
程度である。レベル確度の向上のみに着目してその他の
特性の劣化は度外視した特殊な窓関数でも採用しない限
り、FFTポイント数を如何に増大しても、このスカロ
ップ・ロスは小さくはならない性質のものである。
リエ変換された周波数領域データには、スカロップ・ロ
スと称するレベル変動(低下)が必然的に付随する。こ
こで、スカロップ・ロスについて図2を参照して簡単に
説明する。高速フーリエ変換においては、信号に窓関数
を掛けることにより周波数スペクトルの漏れが少なくな
るようにしている。ところが、この場合、周波数は量子
化されているので、観測される信号の周波数foが量子
化された周波数とたまたま共通している場合に限って正
しいレベルが観測されるに過ぎず、量子化された周波数
の間の周波数における信号は如何に補間しようとも必ず
真のレベルより低いレベルとして観測されることとな
る。周波数foのところのレベルと量子化された周波数
1単位の1/2だけずれたところのレベルとの間の差を
スカロープ・ロスと称している。このロスは概ね1dB
程度である。レベル確度の向上のみに着目してその他の
特性の劣化は度外視した特殊な窓関数でも採用しない限
り、FFTポイント数を如何に増大しても、このスカロ
ップ・ロスは小さくはならない性質のものである。
【0004】この発明は、上述の通りの欠点、問題を解
消、除去したスペクトラム・アナライザを提供しようと
するものである。
消、除去したスペクトラム・アナライザを提供しようと
するものである。
【0005】
【課題を解決するための手段】スペクトラム・アナライ
ザにおいて、先ず分析されるべき時間領域データを入力
し、次いで時間領域データ入力に窓関数を掛け、その結
果のものにFFT演算処理を施して時間領域データを周
波数領域データに変換し、この周波数領域データに補間
処理或いはピーク値補正処理を施し、補正処理を施され
たデータを表示用メモリに転送し、表示する。
ザにおいて、先ず分析されるべき時間領域データを入力
し、次いで時間領域データ入力に窓関数を掛け、その結
果のものにFFT演算処理を施して時間領域データを周
波数領域データに変換し、この周波数領域データに補間
処理或いはピーク値補正処理を施し、補正処理を施され
たデータを表示用メモリに転送し、表示する。
【0006】
【実施例】この発明の実施例を図1を用いて説明をす
る。図自体はこの発明と従来例との間に相違するところ
はない。メモリ部11に記憶された内容についての内臓
コンピュータ12による処理の流れの概要は以下の通り
である。 第1 時間領域データの入力 第2 時間領域データ入力に窓関数を掛ける 第3 FFT演算処理を施し、時間領域データを周波数
領域データに変換する 第4 補間処理或いはピーク値補正処理を施す 第5 補正処理を施されたデータを表示用メモリに転送
し、表示する 高速フーリエ変換においては、信号に窓関数を掛けるこ
とにより周波数スペクトルの漏れを少なくなるようにし
ているが、この場合、周波数が量子化されているので、
観測される信号の周波数foが量子化された周波数とた
またま共通している場合に限って正しいレベルが観測さ
れるに過ぎず、量子化された周波数の間の周波数におけ
る信号は如何に補間しようとも必ず真のレベルより低い
レベルとして観測されるに到る、ということについては
上述した。図2には、FFT処理された周波数領域デー
タのピークの近傍が示されている。周波数foのところ
のレベルと量子化された周波数1単位の1/2だけずれ
たところ即ち周波数f-0.5およびf0.5のところのレベ
ルとの間の差をスカロープ・ロスと言うのであるが、こ
の発明はこの周波数領域データのピークの近傍を関数近
似により補間することにより真のピーク値を近似しよう
とするものである。
る。図自体はこの発明と従来例との間に相違するところ
はない。メモリ部11に記憶された内容についての内臓
コンピュータ12による処理の流れの概要は以下の通り
である。 第1 時間領域データの入力 第2 時間領域データ入力に窓関数を掛ける 第3 FFT演算処理を施し、時間領域データを周波数
領域データに変換する 第4 補間処理或いはピーク値補正処理を施す 第5 補正処理を施されたデータを表示用メモリに転送
し、表示する 高速フーリエ変換においては、信号に窓関数を掛けるこ
とにより周波数スペクトルの漏れを少なくなるようにし
ているが、この場合、周波数が量子化されているので、
観測される信号の周波数foが量子化された周波数とた
またま共通している場合に限って正しいレベルが観測さ
れるに過ぎず、量子化された周波数の間の周波数におけ
る信号は如何に補間しようとも必ず真のレベルより低い
レベルとして観測されるに到る、ということについては
上述した。図2には、FFT処理された周波数領域デー
タのピークの近傍が示されている。周波数foのところ
のレベルと量子化された周波数1単位の1/2だけずれ
たところ即ち周波数f-0.5およびf0.5のところのレベ
ルとの間の差をスカロープ・ロスと言うのであるが、こ
の発明はこの周波数領域データのピークの近傍を関数近
似により補間することにより真のピーク値を近似しよう
とするものである。
【0007】この発明を更に詳しく説明する。入力信号
fS1は、上述の通り、入力端子1に入力され、第1の低
域フィルタ2から第2の低域フィルタ8に到る周波数変
換処理を施されて信号fS3とされる。この信号fS3はア
ナログ−ディジタル変換器9に印加されてディジタルI
F信号fIFDに変換される。このアナログ−ディジタル
変換された時間領域データであるディジタル・データI
F信号fIFDは一時的にメモリ11に記憶される。上記
メモリ11に記憶されたディジタル・データは、先ず、
窓関数としてガウス関数を掛ける掛算器に印加される。
ところで、ガウス関数は下記のとおりのものであり、そ
して時間領域データが1周期の正数倍の場合、正数倍で
はない場合の周波数領域データは図3に示される通りで
ある。
fS1は、上述の通り、入力端子1に入力され、第1の低
域フィルタ2から第2の低域フィルタ8に到る周波数変
換処理を施されて信号fS3とされる。この信号fS3はア
ナログ−ディジタル変換器9に印加されてディジタルI
F信号fIFDに変換される。このアナログ−ディジタル
変換された時間領域データであるディジタル・データI
F信号fIFDは一時的にメモリ11に記憶される。上記
メモリ11に記憶されたディジタル・データは、先ず、
窓関数としてガウス関数を掛ける掛算器に印加される。
ところで、ガウス関数は下記のとおりのものであり、そ
して時間領域データが1周期の正数倍の場合、正数倍で
はない場合の周波数領域データは図3に示される通りで
ある。
【0008】
【数1】 ただし、N=FFTポイント数、α=3.5、θ=ラジ
アン、*=畳み込み、−N/2≦n≦N/2である。 次いで、高速フーリエ変換処理器において、掛算器にお
いてディジタル・データに対して窓関数を掛けた結果に
対して高速フーリエ変換を施す。最後に、補間処理器に
おいて、高速フーリエ変換処理結果に対して3次関数近
似による補間処理を施し、この補間処理結果を表示部1
3に表示する。
アン、*=畳み込み、−N/2≦n≦N/2である。 次いで、高速フーリエ変換処理器において、掛算器にお
いてディジタル・データに対して窓関数を掛けた結果に
対して高速フーリエ変換を施す。最後に、補間処理器に
おいて、高速フーリエ変換処理結果に対して3次関数近
似による補間処理を施し、この補間処理結果を表示部1
3に表示する。
【0009】ここで、近似関数の次数についてである
が、
が、
【0010】
【数2】 となり、結局2次関数となる。以上の結果から、dB単
位に変換されたデータにより近似する場合は、関数の次
数は2次以上であればよいことがわかる。
位に変換されたデータにより近似する場合は、関数の次
数は2次以上であればよいことがわかる。
【0011】
【発明の効果】従来、スカロップ・ロスは概ね1dB程
度であったところ、この発明によれば、図4に示される
程度のロスとすることができた。
度であったところ、この発明によれば、図4に示される
程度のロスとすることができた。
【図1】スペクトラム・アナライザのブロック図。
【図2】FFT処理された周波数領域データのピークの
近傍を示す図。
近傍を示す図。
【図3】時間領域データが1周期の正数倍の場合、正数
倍ではない場合の周波数領域データを示す図。
倍ではない場合の周波数領域データを示す図。
【図4】スカロップ・ロスを示す図。
9 アナログ−ディジタル変換器 11 メモリ 12 内臓コンピュータ処理部 13 表示器
Claims (1)
- 【特許請求の範囲】 【請求項1】 入力信号をアナログ−ディジタル変換す
るアナログ−ディジタル変換器を具備し、アナログ−デ
ィジタル変換されたディジタル・データを一時的に記憶
するメモリを具備し、上記メモリに記憶されたディジタ
ル・データに対し窓関数としてガウス関数を掛ける掛算
器を具備し、ディジタル・データに対して窓関数を掛け
た結果のものに対して高速フーリエ変換処理を施す高速
フーリエ変換器を具備し、高速フーリエ変換処理結果に
対して3次関数近似による補間処理を施す補間処理器を
具備し、補間処理結果を表示する表示器を具備すること
を特徴とするスペクトラム・アナライザ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP17676591A JP2904425B2 (ja) | 1991-07-17 | 1991-07-17 | スペクトラム・アナライザ |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP17676591A JP2904425B2 (ja) | 1991-07-17 | 1991-07-17 | スペクトラム・アナライザ |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH0526919A true JPH0526919A (ja) | 1993-02-05 |
JP2904425B2 JP2904425B2 (ja) | 1999-06-14 |
Family
ID=16019430
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP17676591A Expired - Fee Related JP2904425B2 (ja) | 1991-07-17 | 1991-07-17 | スペクトラム・アナライザ |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2904425B2 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003222646A (ja) * | 2002-01-31 | 2003-08-08 | Mitsubishi Electric Corp | 周波数分析装置 |
JP2006516722A (ja) * | 2003-01-10 | 2006-07-06 | テラビュー リミテッド | イメージング方法及びその関連装置 |
JP2007199030A (ja) * | 2006-01-30 | 2007-08-09 | Hioki Ee Corp | 表示制御装置 |
JPWO2006019196A1 (ja) * | 2004-08-20 | 2008-05-08 | 日本電気株式会社 | 電磁妨害低減量の計算方法、電磁妨害低減量計算装置、計算プログラム及び電子回路 |
-
1991
- 1991-07-17 JP JP17676591A patent/JP2904425B2/ja not_active Expired - Fee Related
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003222646A (ja) * | 2002-01-31 | 2003-08-08 | Mitsubishi Electric Corp | 周波数分析装置 |
JP2006516722A (ja) * | 2003-01-10 | 2006-07-06 | テラビュー リミテッド | イメージング方法及びその関連装置 |
US9829433B2 (en) | 2003-01-10 | 2017-11-28 | Teraview Limited | Imaging techniques and associated apparatus |
JPWO2006019196A1 (ja) * | 2004-08-20 | 2008-05-08 | 日本電気株式会社 | 電磁妨害低減量の計算方法、電磁妨害低減量計算装置、計算プログラム及び電子回路 |
JP2007199030A (ja) * | 2006-01-30 | 2007-08-09 | Hioki Ee Corp | 表示制御装置 |
Also Published As
Publication number | Publication date |
---|---|
JP2904425B2 (ja) | 1999-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Grandke | Interpolation algorithms for discrete Fourier transforms of weighted signals | |
JP3338370B2 (ja) | 周波数分析方法及びこの方法を用いた掃引型スペクトラム・アナライザ | |
US4896285A (en) | Calculation of filter factors for digital filter | |
US5576978A (en) | High resolution frequency analyzer and vector spectrum analyzer | |
US7161511B2 (en) | Linearization system and method | |
US20040083069A1 (en) | Method for optimum spectrum analysis | |
US5117179A (en) | Swept signal analysis instrument and method | |
EP0945731B1 (en) | Real-time signal analyzer | |
TW438981B (en) | Method for separating a desired signal from an interfering signal | |
JPH07294566A (ja) | 高分解能周波数分析装置を用いたドップラ補償装置 | |
CN113361331A (zh) | 基于加窗插值fft的工频干扰消除方法、系统和介质 | |
US6014620A (en) | Power spectral density estimation method and apparatus using LPC analysis | |
EP1248395B1 (en) | Reference frequency spur cancellation in synthesized measurement receivers | |
US20210231789A1 (en) | Doppler signal processing device and method thereof for interference suppression | |
JPH0526919A (ja) | スペクトラム・アナライザ | |
EP0841569B1 (en) | A spectrum analyser | |
JPH07235861A (ja) | 加重平均を使用したサンプリング周波数変換方法 | |
JP2867769B2 (ja) | 音響測定方法およびその装置 | |
JP2004150825A (ja) | スペクトル分析装置およびスペクトル分析方法 | |
CN109541309B (zh) | 一种频谱分析仪及其信号处理方法 | |
JP4344356B2 (ja) | 検波装置、方法、プログラム、記録媒体 | |
US5903483A (en) | Frequency domain kernel processor | |
CN114866099B (zh) | 滤波方法、装置和电子设备 | |
JPH0526921A (ja) | スペクトラム・アナライザ | |
US6690977B1 (en) | Data processor and data processing method for wavemeter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19990223 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |