JP3788117B2 - レーダ信号処理装置 - Google Patents

レーダ信号処理装置 Download PDF

Info

Publication number
JP3788117B2
JP3788117B2 JP21371699A JP21371699A JP3788117B2 JP 3788117 B2 JP3788117 B2 JP 3788117B2 JP 21371699 A JP21371699 A JP 21371699A JP 21371699 A JP21371699 A JP 21371699A JP 3788117 B2 JP3788117 B2 JP 3788117B2
Authority
JP
Japan
Prior art keywords
frequency
circuit
phase correction
target
reception signal
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.)
Expired - Fee Related
Application number
JP21371699A
Other languages
English (en)
Other versions
JP2001042031A (ja
Inventor
直樹 橋田
隆三郎 臼井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP21371699A priority Critical patent/JP3788117B2/ja
Publication of JP2001042031A publication Critical patent/JP2001042031A/ja
Application granted granted Critical
Publication of JP3788117B2 publication Critical patent/JP3788117B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Description

【0001】
【発明の属する技術分野】
この発明は、例えば、高分解能レーダ装置のレーダ信号処理装置における目標受信信号の位相ずれの補正に関するものである。
【0002】
【従来の技術】
図12は高分解能レーダ装置のレーダ信号処理装置の構成図であり、図中、1はレーダ装置から入力された目標受信信号SMを内部で処理できるデータ形式に変換するデータインタフェース部、2はデータインタフェース部1で変換された目標受信信号をパルス圧縮するパルス圧縮部、3はパルス圧縮部2でパルス圧縮された目標受信信号の時間による距離ずれを補正する距離補正部、4は距離補正部3で補正された目標受信信号の時間による位相ずれを補正する位相補正部、5は位相補正部4で補正された目標受信信号のドップラ周波数を分離する周波数分析部、6は周波数分析部5で周波数分析された目標受信信号の周波数スペクトルを画像データに変換する検波部、7は検波部6で得られた画像データと表示器のインタフェースを調整し表示画像データを生成する表示器インタフェース部、SMはレーダ装置から入力された目標受信信号、RSは距離補正部3で時間による距離ずれが補正された目標受信信号、RGはレーダと目標重心との初期距離、RDは位相補正部4で時間による位相ずれが補正された目標受信信号、Dは表示画像データである。
【0003】
図13は図12のレーダ信号処理装置における従来の位相補正部4の構成図であり、図中、RS,RG,RD及び4は図12と同じであり、8は距離補正部3より出力された時間による距離ずれが補正された目標受信信号RS及びレーダと目標重心との初期距離RGを格納するバッファ回路、9はバッファ回路8より出力されたレーダと目標重心との初期距離RGにおける目標受信信号を時間方向に小区間で周波数分析する区分周波数分析回路、10は区分周波数分析回路9で得られた周波数と振幅の波形に対して振幅値が最大となる周波数を基準点周波数として検出する振幅値最大検出回路、11は振幅値最大検出回路10で検出された基準点周波数の時間方向に対する軌跡を平滑化する平滑化回路、12は平滑化回路11で平滑化された基準点周波数の時間方向に対する軌跡から位相補正量を算出する位相補正量算出回路、13は位相補正量算出回路12で算出された位相補正量を用いてバッファ回路8より出力された目標受信信号の位相を補正する位相補正回路、GSはレーダと目標重心との初期距離RGにおける目標受信信号である。
【0004】
次に、動作について説明する。レーダ装置から入力された目標受信信号SMはデータインタフェース部1で内部で処理できるデータ形式に変換され、パルス圧縮部2でパルス圧縮された後、距離補正部3で時間による距離ずれを補正し、目標受信信号RSとして位相補正部4に出力する。また、距離補正部3ではレーダと目標重心との初期距離RGを算出し、位相補正部4に出力する。位相補正部4では時間による距離ずれが補正された目標受信信号RS及びレーダと目標重心との初期距離RGを用いて目標受信信号RSの時間による位相ずれを補正し、目標受信信号RDとして周波数分析部5に出力する。この目標受信信号RDは周波数分析部5で周波数分析されることによって周波数スペクトルに変換され、検波部6で画像データに変換された後、表示器インタフェース部7で表示器とのインタフェースを調整し、表示画像データDとして出力される。
【0005】
次に、位相補正部4の動作について説明する。距離補正部3から入力された時間による距離ずれが補正された目標受信信号RS及びレーダと目標重心との初期距離RGはバッファ回路8に格納され、目標受信信号RS及びレーダと目標重心との初期距離RGにおける目標受信信号GSとして出力される。このレーダと目標重心との初期距離RGにおける目標受信信号GSは区分周波数分析回路9で時間方向に小区間で周波数分析され、得られた周波数と振幅の波形に対して振幅値最大検出回路10で振幅値が最大となる周波数を基準点周波数として検出した後、平滑化回路11に出力される。平滑化回路11では基準点周波数の時間方向に対する軌跡を平滑化し、平滑化された軌跡から位相補正量算出回路12で位相補正量を算出する。位相補正回路13は位相補正量算出回路12で算出した位相補正量を用いてバッファ回路8から出力された時間による距離ずれが補正された目標受信信号RSの位相を補正し、時間による位相ずれが補正された目標受信信号RDとして周波数分析部5へ出力する。
【0006】
更に、位相補正部4を図14を用いて説明する。図14は位相補正部4の処理方法を示した図である。時間による距離ずれが補正された目標受信信号RSをSi,j (ここで、iはレンジビン番号、jはパルスヒット番号、i,jは自然数である。)、レーダと目標重心との初期距離RGの存在するレンジビン番号をrと定義するとレーダと目標重心との初期距離RGにおける目標受信信号GSはSr,jと表され、図14(a)のような波形が得られる。Sr,jに対し、区分周波数分析回路9で時間方向(パルスヒット方向)に小区間で周波数分析を行うと図14(b)のような波形が得られ、周波数fmと振幅Am kの関係(ここで、kは区分周波数分析番号、mは周波数ビン番号、k,mは自然数である。)は“数1”で表される。
【0007】
【数1】
Figure 0003788117
【0008】
振幅値最大検出回路10で各区分周波数分析番号kに対し、振幅Am kが最大値をとる時の周波数を検出し、それを基準点周波数fkとすると、時間tkと基準点周波数fkの関係は図14(c)のプロットのようになる。平滑化回路11で図14(c)のプロットに対し、平滑化を行うと図14(c)の実線のような波形が得られ、時間tkと周波数f' kの関係は“数2”で表される。
【0009】
【数2】
Figure 0003788117
【0010】
位相補正量算出回路12では、位相補正量Wjを“数3”で算出する。
【0011】
【数3】
Figure 0003788117
【0012】
位相補正回路13では、Si,j の位相を位相補正量Wjを用いて“数4”で補正する。但し、時間による位相ずれが補正された目標受信信号RDをS'i,j と定義する。
【0013】
【数4】
Figure 0003788117
【0014】
【発明が解決しようとする課題】
上記のような従来のレーダ信号処理装置では、区分周波数分析後の周波数と振幅の波形が多峰で、かつ振幅値が最大となる峰の位置が時間毎に大きく変動するような場合、検出した基準点周波数がふらついて正確な位相補正量を算出することができず、画像がぼけたり、にじんだりするという課題があった。
【0015】
この発明はかかる課題を解決するためになされたもので、検出した基準点周波数がふらついて正確な位相補正量を算出することができず、画像がぼけたり、にじんだりするということを防止するレーダ信号処理装置を得ることを目的とする。
【0016】
【課題を解決するための手段】
第1の発明に係るレーダ信号処理装置は、上記位相補正部を上記距離補正部で補正された目標受信信号及びレーダと目標重心との初期距離を格納するバッファ回路と、上記バッファ回路より出力されたレーダと目標重心との初期距離における目標受信信号を時間方向に小区間で周波数分析する区分周波数分析回路と、上記区分周波数分析回路で得られた周波数と振幅の波形について高周波成分を除去する低域通過フィルタと、上記低域通過フィルタにより高周波成分が除去された周波数と振幅の波形についてデータを切り出すデータ切り出し回路と、上記データ切り出し回路でデータを切り出す際に必要な閾値を設定する閾値設定回路と、上記データ切り出し回路で切り出された周波数と振幅の波形について極大となる点を検出しその時の周波数を算出する極大点周波数検出回路と、上記極大点周波数検出回路で検出された周波数の平均値を基準点周波数とする平均周波数検出回路と、上記平均周波数検出回路で検出された基準点周波数の時間方向に対する軌跡を平滑化する平滑化回路と、上記平滑化回路で平滑化された基準点周波数の時間方向に対する軌跡から位相補正量を算出する位相補正量算出回路と、上記位相補正量算出回路で算出された位相補正量を用いて上記バッファ回路から出力された目標受信信号の位相を補正する位相補正回路とによって構成したものである。
【0017】
また、第2の発明に係るレーダ信号処理装置は上記位相補正部の基準点周波数検出手段として、上記平均周波数検出回路を上記極大点周波数検出回路で検出された周波数の中心値を基準点周波数とする中心周波数検出回路に置き換えたものである。
【0018】
第3の発明に係るレーダ信号処理装置は上記位相補正部の基準点周波数検出手段として、上記平均周波数検出回路を上記極大点周波数検出回路で検出された周波数の最小値を基準点周波数とする最小周波数検出回路に置き換えたものである。
【0019】
また、第4の発明に係るレーダ信号処理装置は上記位相補正部の基準点周波数検出手段として、上記平均周波数検出回路を上記極大点周波数検出回路で検出された周波数の最大値を基準点周波数とする最大周波数検出回路に置き換えたものである。
【0020】
第5の発明に係るレーダ信号処理装置は上記位相補正部の基準点周波数検出手段として、上記平均周波数検出回路を上記極大点周波数検出回路で検出された周波数の中央値を基準点周波数とする中央周波数検出回路に置き換えたものである。
【0021】
【発明の実施の形態】
実施の形態1.
図1は、図12に示すレーダ信号処理装置におけるこの発明の位相補正部の実施の一形態を示すものである。図において、4,RS,RG及びRDは図12及び図13と同じである。また、図において、8,9,11,12,13及びGSは図13と同じである。22は区分周波数分析回路9で得られた周波数と振幅の波形について高周波成分を除去する低域通過フィルタ、14は低域通過フィルタ22により高周波成分が除去された周波数と振幅の波形についてデータを切り出すデータ切り出し回路、15はデータ切り出し回路14でデータを切り出す際に必要な閾値を設定する閾値設定回路、16はデータ切り出し回路14で切り出された周波数と振幅の波形について極大となる点を検出しその時の周波数を算出する極大点周波数検出回路、17は極大点周波数検出回路16で検出された周波数の平均値を基準点周波数とする平均周波数検出回路である。
【0022】
次に、上記図1のように構成された位相補正部4の動作について説明する。距離補正部3から入力された時間による距離ずれが補正された目標受信信号RS及びレーダと目標重心との初期距離RGはバッファ回路8に格納され、目標受信信号RS及びレーダと目標重心との初期距離RGにおける目標受信信号GSとして出力される。このレーダと目標重心との初期距離RGにおける目標受信信号GSは区分周波数分析回路9で時間方向に小区間で周波数分析され、低域通過フィルタ22で高周波成分を除去した後、得られた周波数と振幅の波形に対してデータ切り出し回路14で各振幅値が閾値設定回路15で設定した閾値を越える領域についてデータが切り出される。閾値設定回路15では、例えば、メインローブレベルとサイドローブレベルの間に一定の閾値を設定する固定スレッショルドやアダプティブに閾値を設定するCFAR(Constant False Alarm Rate)等を用いて、各区分周波数分析後の波形毎に閾値を設定する。極大点周波数検出回路16ではデータ切り出し回路14で切り出された周波数と振幅の波形について全ての極大点を検出し、その時の周波数を算出する。平均周波数検出回路17では極大点周波数検出回路16で算出された全極大点における周波数に対し、それらの平均周波数を基準点周波数として検出した後、平滑化回路11に出力する。平滑化回路11では基準点周波数の時間方向に対する軌跡を平滑化し、平滑化された軌跡から位相補正量算出回路12で位相補正量を算出する。位相補正回路13は位相補正量算出回路12で算出した位相補正量を用いてバッファ回路8から出力された時間による距離ずれが補正された目標受信信号RSの位相を補正し、時間による位相ずれが補正された目標受信信号RDとして周波数分析部5へ出力する。
【0023】
次に、上記図1のように構成された位相補正部4を図14,図6及び図7を用いて説明する。図14は位相補正部4の処理方法、図6は区分周波数分析後の波形、低域通過フィルタ通過後の波形及び極大点周波数の検出方法を示した図、図7は平均周波数検出回路17における基準点周波数の検出方法を示した図である。時間による距離ずれが補正された目標受信信号RSをSi,j (ここで、iはレンジビン番号、jはパルスヒット番号、i,jは自然数である。)、レーダと目標重心との初期距離RGの存在するレンジビン番号をrと定義するとレーダと目標重心との初期距離RGにおける目標受信信号GSはSr,jと表され、図14(a)のような波形が得られる。Sr,jに対し、区分周波数分析回路9で時間方向(パルスヒット方向)に小区間で周波数分析を行うと図14(b)のような波形が得られ、周波数fmと振幅Am kの関係(ここで、kは区分周波数分析番号、mは周波数ビン番号、k,mは自然数である。)は“数1”で表される。
【0024】
各区分周波数分析番号kにおいて周波数fmと振幅Am kの波形は図6(a)のように表される。図6(b)は、図6(a)で表される波形に対して、低域通過フィルタ22により高周波成分を除去したものである。図6(c)では、図6(b)で表される波形に対し、各振幅値が閾値設定回路15で設定した閾値(例えば、図6(c)中の一点破線u)を越える周波数領域についてデータ切り出し回路14でデータ切り出しを行い、次に切り出された周波数領域について極大点周波数検出回路16で全ての極大点を検出し、その時の周波数f'n k(ここで、n=1,・・・,Nであり、nは極大点番号、Nは極大点数、n,Nは自然数である。)を算出する。
【0025】
平均周波数検出回路17では、図7のように極大点周波数検出回路16で検出された極大点周波数f'n kを用いて、“数5”で基準点周波数fkを算出する。
【0026】
【数5】
Figure 0003788117
【0027】
時間tkと基準点周波数fkの関係は図14(c)のプロットのようになる。平滑化回路11で図14(c)のプロットに対し、平滑化を行うと図14(c)の実線のような波形が得られ、時間tkと周波数f' kの関係は“数2”で表される。
【0028】
位相補正量算出回路12では、位相補正量Wjを“数3”で算出する。
【0029】
位相補正回路13では、Si,j の位相を位相補正量Wjを用いて“数4”で補正する。但し、時間による位相ずれが補正された目標受信信号RDをS'i,j と定義する。
【0030】
実施の形態2.
図2に示される実施の形態では、上記実施の形態1における位相補正部4の基準点周波数検出手段として、平均周波数検出回路17を極大点周波数検出回路16で検出された周波数の中心値を基準点周波数とする中心周波数検出回路18に置き換えている。このような実施態様によれば、区分周波数分析後の周波数と振幅の波形に対し、低域通過フィルタ22で高周波成分を除去した後、極大点の周波数のみに着目して基準点周波数を決定するため、高周波成分による極大点の誤検出を低減し、波形の振幅値に依存しない安定した基準点周波数を検出することができる。更に、極大点の周波数から中心周波数を検出するので、計算量の削減が可能である。
【0031】
次に、上記図2のように構成された位相補正部4の動作について説明する。距離補正部3から入力された時間による距離ずれが補正された目標受信信号RS及びレーダと目標重心との初期距離RGはバッファ回路8に格納され、目標受信信号RS及びレーダと目標重心との初期距離RGにおける目標受信信号GSとして出力される。このレーダと目標重心との初期距離RGにおける目標受信信号GSは区分周波数分析回路9で時間方向に小区間で周波数分析され、低域通過フィルタ22で高周波成分を除去した後、得られた周波数と振幅の波形に対してデータ切り出し回路14で各振幅値が閾値設定回路15で設定した閾値を越える領域についてデータが切り出される。閾値設定回路15では、例えば、メインローブレベルとサイドローブレベルの間に一定の閾値を設定する固定スレッショルドやアダプティブに閾値を設定するCFAR(Constant False Alarm Rate)等を用いて、各区分周波数分析後の波形毎に閾値を設定する。極大点周波数検出回路16ではデータ切り出し回路14で切り出された周波数と振幅の波形について全ての極大点を検出し、その時の周波数を算出する。中心周波数検出回路18では極大点周波数検出回路16で算出された全極大点における周波数に対し、それらの中心周波数を基準点周波数として検出した後、平滑化回路11に出力する。平滑化回路11では基準点周波数の時間方向に対する軌跡を平滑化し、平滑化された軌跡から位相補正量算出回路12で位相補正量を算出する。位相補正回路13は位相補正量算出回路12で算出した位相補正量を用いてバッファ回路8から出力された時間による距離ずれが補正された目標受信信号RSの位相を補正し、時間による位相ずれが補正された目標受信信号RDとして周波数分析部5へ出力する。
【0032】
次に、上記図2のように構成された位相補正部4を図6及び図8を用いて説明する。図8は中心周波数検出回路18における基準点周波数の検出方法を示した図である。各区分周波数分析番号kにおいて周波数fmと振幅Am kの波形は図6(a)のように表される。図6(b)は、図6(a)で表される波形に対して、低域通過フィルタ22により高周波成分を除去したものである。図6(c)では、図6(b)で表される波形に対し、各振幅値が閾値設定回路15で設定した閾値(例えば、図6(c)中の一点破線u)を越える周波数領域についてデータ切り出し回路14でデータ切り出しを行い、次に切り出された周波数領域について極大点周波数検出回路16で全ての極大点を検出し、その時の周波数f'n k(ここで、n=1,・・・,Nであり、nは極大点番号、Nは極大点数、n,Nは自然数である。)を算出する。
【0033】
中心周波数検出回路18では、図8のように極大点周波数検出回路16で検出された極大点周波数f'n kを用いて、“数6”で基準点周波数fkを算出する。
【0034】
【数6】
Figure 0003788117
【0035】
実施の形態3.
図3に示される実施の形態では、上記実施の形態1における位相補正部4の基準点周波数検出手段として、平均周波数検出回路17を極大点周波数検出回路16で検出された周波数の最小値を基準点周波数とする最小周波数検出回路19に置き換えている。このような実施態様によれば、区分周波数分析後の周波数と振幅の波形に対し、低域通過フィルタ22で高周波成分を除去した後、極大点の周波数のみに着目して基準点周波数を決定するため、高周波成分による極大点の誤検出を低減し、波形の振幅値に依存しない安定した基準点周波数を検出することができる。更に、極大点の周波数から最小周波数を検出するので、計算量の削減が可能である。
【0036】
次に、上記図3のように構成された位相補正部4の動作について説明する。距離補正部3から入力された時間による距離ずれが補正された目標受信信号RS及びレーダと目標重心との初期距離RGはバッファ回路8に格納され、目標受信信号RS及びレーダと目標重心との初期距離RGにおける目標受信信号GSとして出力される。このレーダと目標重心との初期距離RGにおける目標受信信号GSは区分周波数分析回路9で時間方向に小区間で周波数分析され、低域通過フィルタ22で高周波成分を除去した後、得られた周波数と振幅の波形に対してデータ切り出し回路14で各振幅値が閾値設定回路15で設定した閾値を越える領域についてデータが切り出される。閾値設定回路15では、例えば、メインローブレベルとサイドローブレベルの間に一定の閾値を設定する固定スレッショルドやアダプティブに閾値を設定するCFAR(Constant False Alarm Rate)等を用いて、各区分周波数分析後の波形毎に閾値を設定する。極大点周波数検出回路16ではデータ切り出し回路14で切り出された周波数と振幅の波形について全ての極大点を検出し、その時の周波数を算出する。最小周波数検出回路19では極大点周波数検出回路16で算出された全極大点における周波数に対し、それらの最小周波数を基準点周波数として検出した後、平滑化回路11に出力する。平滑化回路11では基準点周波数の時間方向に対する軌跡を平滑化し、平滑化された軌跡から位相補正量算出回路12で位相補正量を算出する。位相補正回路13は位相補正量算出回路12で算出した位相補正量を用いてバッファ回路8から出力された時間による距離ずれが補正された目標受信信号RSの位相を補正し、時間による位相ずれが補正された目標受信信号RDとして周波数分析部5へ出力する。
【0037】
次に、上記図3のように構成された位相補正部4を図6及び図9を用いて説明する。図9は最小周波数検出回路19における基準点周波数の検出方法を示した図である。各区分周波数分析番号kにおいて周波数fmと振幅Am kの波形は図6(a)のように表される。図6(b)は、図6(a)で表される波形に対して、低域通過フィルタ22により高周波成分を除去したものである。図6(c)では、図6(b)で表される波形に対し、各振幅値が閾値設定回路15で設定した閾値(例えば、図6(c)中の一点破線u)を越える周波数領域についてデータ切り出し回路14でデータ切り出しを行い、次に切り出された周波数領域について極大点周波数検出回路16で全ての極大点を検出し、その時の周波数f'n k(ここで、n=1,・・・,Nであり、nは極大点番号、Nは極大点数、n,Nは自然数である。)を算出する。
【0038】
最小周波数検出回路19では、図9のように極大点周波数検出回路16で検出された極大点周波数f'n kを用いて、“数7”で基準点周波数fkを算出する。
【0039】
【数7】
Figure 0003788117
【0040】
実施の形態4.
図4に示される実施の形態では、上記実施の形態1における位相補正部4の基準点周波数検出手段として、平均周波数検出回路17を極大点周波数検出回路16で検出された周波数の最大値を基準点周波数とする最大周波数検出回路20に置き換えている。このような実施態様によれば、区分周波数分析後の周波数と振幅の波形に対し、低域通過フィルタ22で高周波成分を除去した後、極大点の周波数のみに着目して基準点周波数を決定するため、高周波成分による極大点の誤検出を低減し、波形の振幅値に依存しない安定した基準点周波数を検出することができる。更に、極大点の周波数から最大周波数を検出するので、計算量の削減が可能である。
【0041】
次に、上記図4のように構成された位相補正部4の動作について説明する。距離補正部3から入力された時間による距離ずれが補正された目標受信信号RS及びレーダと目標重心との初期距離RGはバッファ回路8に格納され、目標受信信号RS及びレーダと目標重心との初期距離RGにおける目標受信信号GSとして出力される。このレーダと目標重心との初期距離RGにおける目標受信信号GSは区分周波数分析回路9で時間方向に小区間で周波数分析され、低域通過フィルタ22で高周波成分を除去した後、得られた周波数と振幅の波形に対してデータ切り出し回路14で各振幅値が閾値設定回路15で設定した閾値を越える領域についてデータが切り出される。閾値設定回路15では、例えば、メインローブレベルとサイドローブレベルの間に一定の閾値を設定する固定スレッショルドやアダプティブに閾値を設定するCFAR(Constant False Alarm Rate)等を用いて、各区分周波数分析後の波形毎に閾値を設定する。極大点周波数検出回路16ではデータ切り出し回路14で切り出された周波数と振幅の波形について全ての極大点を検出し、その時の周波数を算出する。最大周波数検出回路20では極大点周波数検出回路16で算出された全極大点における周波数に対し、それらの最大周波数を基準点周波数として検出した後、平滑化回路11に出力する。平滑化回路11では基準点周波数の時間方向に対する軌跡を平滑化し、平滑化された軌跡から位相補正量算出回路12で位相補正量を算出する。位相補正回路13は位相補正量算出回路12で算出した位相補正量を用いてバッファ回路8から出力された時間による距離ずれが補正された目標受信信号RSの位相を補正し、時間による位相ずれが補正された目標受信信号RDとして周波数分析部5へ出力する。
【0042】
次に、上記図4のように構成された位相補正部4を図6及び図10を用いて説明する。図10は最大周波数検出回路20における基準点周波数の検出方法を示した図である。各区分周波数分析番号kにおいて周波数fmと振幅Am kの波形は図6(a)のように表される。図6(b)は、図6(a)で表される波形に対して、低域通過フィルタ22により高周波成分を除去したものである。図6(c)では、図6(b)で表される波形に対し、各振幅値が閾値設定回路15で設定した閾値(例えば、図6(c)中の一点破線u)を越える周波数領域についてデータ切り出し回路14でデータ切り出しを行い、次に切り出された周波数領域について極大点周波数検出回路16で全ての極大点を検出し、その時の周波数f'n k(ここで、n=1,・・・,Nであり、nは極大点番号、Nは極大点数、n,Nは自然数である。)を算出する。
【0043】
最大周波数検出回路20では、図10のように極大点周波数検出回路16で検出された極大点周波数f'n kを用いて、“数8”で基準点周波数fkを算出する。
【0044】
【数8】
Figure 0003788117
【0045】
実施の形態5.
図5に示される実施の形態では、上記実施の形態1における位相補正部4の基準点周波数検出手段として、平均周波数検出回路17を極大点周波数検出回路16で検出された周波数の中央値を基準点周波数とする中央周波数検出回路21に置き換えている。このような実施態様によれば、区分周波数分析後の周波数と振幅の波形に対し、低域通過フィルタ22で高周波成分を除去した後、極大点の周波数のみに着目して基準点周波数を決定するため、高周波成分による極大点の誤検出を低減し、波形の振幅値に依存しない安定した基準点周波数を検出することができる。更に、極大点の周波数から中央周波数を検出するので、計算量の削減が可能である。
【0046】
次に、上記図5のように構成された位相補正部4の動作について説明する。距離補正部3から入力された時間による距離ずれが補正された目標受信信号RS及びレーダと目標重心との初期距離RGはバッファ回路8に格納され、目標受信信号RS及びレーダと目標重心との初期距離RGにおける目標受信信号GSとして出力される。このレーダと目標重心との初期距離RGにおける目標受信信号GSは区分周波数分析回路9で時間方向に小区間で周波数分析され、低域通過フィルタ22で高周波成分を除去した後、得られた周波数と振幅の波形に対してデータ切り出し回路14で各振幅値が閾値設定回路15で設定した閾値を越える領域についてデータが切り出される。閾値設定回路15では、例えば、メインローブレベルとサイドローブレベルの間に一定の閾値を設定する固定スレッショルドやアダプティブに閾値を設定するCFAR(Constant False Alarm Rate)等を用いて、各区分周波数分析後の波形毎に閾値を設定する。極大点周波数検出回路16ではデータ切り出し回路14で切り出された周波数と振幅の波形について全ての極大点を検出し、その時の周波数を算出する。中央周波数検出回路21では極大点周波数検出回路16で算出された全極大点における周波数に対し、それらの中央周波数を基準点周波数として検出した後、平滑化回路11に出力する。平滑化回路11では基準点周波数の時間方向に対する軌跡を平滑化し、平滑化された軌跡から位相補正量算出回路12で位相補正量を算出する。位相補正回路13は位相補正量算出回路12で算出した位相補正量を用いてバッファ回路8から出力された時間による距離ずれが補正された目標受信信号RSの位相を補正し、時間による位相ずれが補正された目標受信信号RDとして周波数分析部5へ出力する。
【0047】
次に、上記図5のように構成された位相補正部4を図6及び図11を用いて説明する。図11は中央周波数検出回路21における基準点周波数の検出方法を示した図である。各区分周波数分析番号kにおいて周波数fmと振幅Am kの波形は図6(a)のように表される。図6(b)は、図6(a)で表される波形に対して、低域通過フィルタ22により高周波成分を除去したものである。図6(c)では、図6(b)で表される波形に対し、各振幅値が閾値設定回路15で設定した閾値(例えば、図6(c)中の一点破線u)を越える周波数領域についてデータ切り出し回路14でデータ切り出しを行い、次に切り出された周波数領域について極大点周波数検出回路16で全ての極大点を検出し、その時の周波数f'n k(ここで、n=1,・・・,Nであり、nは極大点番号、Nは極大点数、n,Nは自然数である。)を算出する。
【0048】
中央周波数検出回路21では、極大点周波数検出回路16で検出された極大点周波数f'n kを用いて、“数9”で基準点周波数fkを算出する。
【0049】
【数9】
Figure 0003788117
【0050】
“数9”の右辺において極大点数Nが奇数である場合、(1+N)/2=nは整数となるため、該当する極大点周波数f'n kは存在する。従って、容易に基準点周波数fkを算出することができる。しかし、極大点数Nが偶数である場合には、(1+N)/2=nは整数とならない(2で割り切れず、実数となる)ため、該当する極大点周波数f'n kは存在しない。そこで、極大点数Nが偶数の場合には、例えば以下に示す方法のうち、いずれかの方法で基準点周波数fkを算出するようにすればよい。
【0051】
“数10”に(1+N)/2の小数点以下を切り捨てる方法、“数11”に(1+N)/2の小数点以下を切り上げる方法、“数12”に“数10”と“数11”の平均をとる方法を示す。なお、図11では“数12”を用いて基準点周波数fkを算出している。
【0052】
【数10】
Figure 0003788117
【0053】
【数11】
Figure 0003788117
【0054】
【数12】
Figure 0003788117
【0055】
【発明の効果】
第1の発明は以上説明した通り、区分周波数分析後の周波数と振幅の波形から低域通過フィルタにより高周波成分を除去した後、設定した閾値に基づいてデータを切り出し、切り出した波形の極大点に着目して、全極大点における周波数の平均値を基準点周波数として検出するので、まず最初に高周波成分による極大点の誤検出を低減することができる。また、周波数と振幅の波形が多峰で、かつ振幅値が最大となる峰の位置が時間毎に大きく変動するような場合でも変動の影響を受けにくく、安定して基準点周波数を検出することができ、画像のぼけやにじみを除去することができるという効果がある。
【0056】
また、第2の発明は区分周波数分析後の周波数と振幅の波形から低域通過フィルタにより高周波成分を除去した後、設定した閾値に基づいてデータを切り出し、切り出した波形の極大点に着目して、全極大点における周波数の中心値を基準点周波数として検出するので、まず最初に高周波成分による極大点の誤検出を低減することができる。また、計算量を削減できるとともに、周波数と振幅の波形が多峰で、かつ振幅値が最大となる峰の位置が時間毎に大きく変動するような場合でも変動の影響を受けにくく、安定して基準点周波数を検出することができ、画像のぼけやにじみを除去することができるという効果がある。
【0057】
第3の発明は区分周波数分析後の周波数と振幅の波形から低域通過フィルタにより高周波成分を除去した後、設定した閾値に基づいてデータを切り出し、切り出した波形の極大点に着目して、全極大点における周波数の最小値を基準点周波数として検出するので、まず最初に高周波成分による極大点の誤検出を低減することができる。また、計算量を削減できるとともに、周波数と振幅の波形が多峰で、かつ振幅値が最大となる峰の位置が時間毎に大きく変動するような場合でも変動の影響を受けにくく、安定して基準点周波数を検出することができ、画像のぼけやにじみを除去することができるという効果がある。
【0058】
また、第4の発明は区分周波数分析後の周波数と振幅の波形から低域通過フィルタにより高周波成分を除去した後、設定した閾値に基づいてデータを切り出し、切り出した波形の極大点に着目して、全極大点における周波数の最大値を基準点周波数として検出するので、まず最初に高周波成分による極大点の誤検出を低減することができる。また、計算量を削減できるとともに、周波数と振幅の波形が多峰で、かつ振幅値が最大となる峰の位置が時間毎に大きく変動するような場合でも変動の影響を受けにくく、安定して基準点周波数を検出することができ、画像のぼけやにじみを除去することができるという効果がある。
【0059】
第5の発明は区分周波数分析後の周波数と振幅の波形から低域通過フィルタにより高周波成分を除去した後、設定した閾値に基づいてデータを切り出し、切り出した波形の極大点に着目して、全極大点における周波数の中央値を基準点周波数として検出するので、まず最初に高周波成分による極大点の誤検出を低減することができる。また、計算量を削減できるとともに、周波数と振幅の波形が多峰で、かつ振幅値が最大となる峰の位置が時間毎に大きく変動するような場合でも変動の影響を受けにくく、安定して基準点周波数を検出することができ、画像のぼけやにじみを除去することができるという効果がある。
【図面の簡単な説明】
【図1】 この発明の実施の形態1を示す位相補正部の構成図である。
【図2】 この発明の実施の形態2を示す位相補正部の構成図である。
【図3】 この発明の実施の形態3を示す位相補正部の構成図である。
【図4】 この発明の実施の形態4を示す位相補正部の構成図である。
【図5】 この発明の実施の形態5を示す位相補正部の構成図である。
【図6】 区分周波数分析後の波形、低域通過フィルタ通過後の波形及び極大点周波数の検出方法を示す図である。
【図7】 平均周波数検出回路における基準点周波数の検出方法を示す図である。
【図8】 中心周波数検出回路における基準点周波数の検出方法を示す図である。
【図9】 最小周波数検出回路における基準点周波数の検出方法を示す図である。
【図10】 最大周波数検出回路における基準点周波数の検出方法を示す図である。
【図11】 中央周波数検出回路における基準点周波数の検出方法を示す図である。
【図12】 高分解能レーダ装置におけるレーダ信号処理装置の構成図である。
【図13】 従来の位相補正部の構成図である。
【図14】 位相補正部の処理方法を示す図である。
【符号の説明】
1 データインタフェース部、2 パルス圧縮部、3 距離補正部、 4 位相補正部、5 周波数分析部、6 検波部、7 表示器インタフェース部、8 バッファ回路、9 区分周波数分析回路、10 振幅値最大検出回路、11 平滑化回路、12 位相補正量算出回路、13 位相補正回路、14 データ切り出し回路、15 閾値設定回路、16 極大点周波数検出回路、17 平均周波数検出回路、18 中心周波数検出回路、19 最小周波数検出回路、20 最大周波数検出回路、21 中央周波数検出回路、22 低域通過フィルタ。

Claims (5)

  1. レーダ装置から入力された目標受信信号をパルス圧縮するパルス圧縮部と、上記パルス圧縮部でパルス圧縮された目標受信信号の時間による
    距離ずれを補正する距離補正部と、上記距離補正部で補正された目標受信信号の時間による位相ずれを補正する位相補正部と、上記位相補正部で補正された目標受信信号のドップラ周波数を分離する周波数分析部と、上記周波数分析部で周波数分析された目標受信信号の周波数スペクトルを画像データに変換する検波部とを備えたレーダ信号処理装置において、上記位相補正部を上記距離補正部で補正された目標受信信号及びレーダと目標重心との初期距離を格納するバッファ回路と、上記バッファ回路より出力されたレーダと目標重心との初期距離における目標受信信号を時間方向に小区間で周波数分析する区分周波数分析回路と、上記区分周波数分析回路で得られた周波数と振幅の波形について周波数方向における信号の振幅変動のうち変動周期が小さい成分を除去する低域通過フィルタと、上記低域通過フィルタにより上記変動周期が小さい成分が除去された周波数と振幅の波形について閾値を越える周波数領域のデータを切り出すデータ切り出し回路と、上記データ切り出し回路でデータを切り出す際に必要な上記閾値を設定する閾値設定回路と、上記データ切り出し回路で切り出された周波数と振幅の波形について全ての極大となる点を検出しその時の周波数を算出する極大点周波数検出回路と、上記極大点周波数検出回路で検出された周波数の平均値を基準点周波数とする平均周波数検出回路と、上記平均周波数検出回路で検出された基準点周波数の時間方向に対する軌跡を平滑化する平滑化回路と、上記平滑化回路で平滑化された基準点周波数の時間方向に対する軌跡から位相補正量を算出する位相補正量算出回路と、上記位相補正量算出回路で算出された位相補正量を用いて上記バッファ回路から出力された目標受信信号の位相を補正する位相補正回路とによって構成したことを特徴とするレーダ信号処理装置。
  2. レーダ装置から入力された目標受信信号をパルス圧縮するパルス圧縮部と、上記パルス圧縮部でパルス圧縮された目標受信信号の時間による距離ずれを補正する距離補正部と、上記距離補正部で補正された目標受信信号の時間による位相ずれを補正する位相補正部と、上記位相補正部で補正された目標受信信号のドップラ周波数を分離する周波数分析部と、上記周波数分析部で周波数分析された目標受信信号の周波数スペクトルを画像データに変換する検波部とを備えたレーダ信号処理装置において、上記位相補正部を上記距離補正部で補正された目標受信信号及びレーダと目標重心との初期距離を格納するバッファ回路と、上記バッファ回路より出力されたレーダと目標重心との初期距離における目標受信信号を時間方向に小区間で周波数分析する区分周波数分析回路と、上記区分周波数分析回路で得られた周波数と振幅の波形について周波数方向における信号の振幅変動のうち変動周期が小さい成分を除去する低域通過フィルタと、上記低域通過フィルタにより上記変動周期が小さい成分が除去された周波数と振幅の波形について閾値を越える周波数領域のデータを切り出すデータ切り出し回路と、上記データ切り出し回路でデータを切り出す際に必要な上記閾値を設定する閾値設定回路と、上記データ切り出し回路で切り出された周波数と振幅の波形について全ての極大となる点を検出しその時の周波数を算出する極大点周波数検出回路と、上記極大点周波数検出回路で検出された周波数の中心値を基準点周波数とする中心周波数検出回路と、上記中心周波数検出回路で検出された基準点周波数の時間方向に対する軌跡を平滑化する平滑化回路と、上記平滑化回路で平滑化された基準点周波数の時間方向に対する軌跡から位相補正量を算出する位相補正量算出回路と、上記位相補正量算出回路で算出された位相補正量を用いて上記バッファ回路から出力された目標受信信号の位相を補正する位相補正回路とによって構成したことを特徴とするレーダ信号処理装置。
  3. レーダ装置から入力された目標受信信号をパルス圧縮するパルス圧縮部と、上記パルス圧縮部でパルス圧縮された目標受信信号の時間による距離ずれを補正する距離補正部と、上記距離補正部で補正された目標受信信号の時間による位相ずれを補正する位相補正部と、上記位相補正部で補正された目標受信信号のドップラ周波数を分離する周波数分析部と、上記周波数分析部で周波数分析された目標受信信号の周波数スペクトルを画像データに変換する検波部とを備えたレーダ信号処理装置において、上記位相補正部を上記距離補正部で補正された目標受信信号及びレーダと目標重心との初期距離を格納するバッファ回路と、上記バッファ回路より出力されたレーダと目標重心との初期距離における目標受信信号を時間方向に小区間で周波数分析する区分周波数分析回路と、上記区分周波数分析回路で得られた周波数と振幅の波形について周波数方向における信号の振幅変動のうち変動周期が小さい成分を除去する低域通過フィルタと、上記低域通過フィルタにより上記変動周期が小さい成分が除去された周波数と振幅の波形について閾値を越える周波数領域のデータを切り出すデータ切り出し回路と、上記データ切り出し回路でデータを切り出す際に必要な上記閾値を設定する閾値設定回路と、上記データ切り出し回路で切り出された周波数と振幅の波形について全ての極大となる点を検出しその時の周波数を算出する極大点周波数検出回路と、上記極大点周波数検出回路で検出された周波数の最小値を基準点周波数とする最小周波数検出回路と、上記最小周波数検出回路で検出された基準点周波数の時間方向に対する軌跡を平滑化する平滑化回路と、上記平滑化回路で平滑化された基準点周波数の時間方向に対する軌跡から位相補正量を算出する位相補正量算出回路と、上記位相補正量算出回路で算出された位相補正量を用いて上記バッファ回路から出力された目標受信信号の位相を補正する位相補正回路とによって構成したことを特徴とするレーダ信号処理装置。
  4. レーダ装置から入力された目標受信信号をパルス圧縮するパルス圧縮部と、上記パルス圧縮部でパルス圧縮された目標受信信号の時間による距離ずれを補正する距離補正部と、上記距離補正部で補正された目標受信信号の時間による位相ずれを補正する位相補正部と、上記位相補正部で補正された目標受信信号のドップラ周波数を分離する周波数分析部と、上記周波数分析部で周波数分析された目標受信信号の周波数スペクトルを画像データに変換する検波部とを備えたレーダ信号処理装置において、上記位相補正部を上記距離補正部で補正された目標受信信号及びレーダと目標重心との初期距離を格納するバッファ回路と、上記バッファ回路より出力されたレーダと目標重心との初期距離における目標受信信号を時間方向に小区間で周波数分析する区分周波数分析回路と、上記区分周波数分析回路で得られた周波数と振幅の波形について周波数方向における信号の振幅変動のうち変動周期が小さい成分を除去する低域通過フィルタと、上記低域通過フィルタにより上記変動周期が小さい成分が除去された周波数と振幅の波形について閾値を越える周波数領域のデータを切り出すデータ切り出し回路と、上記データ切り出し回路でデータを切り出す際に必要な上記閾値を設定する閾値設定回路と、上記データ切り出し回路で切り出された周波数と振幅の波形について全ての極大となる点を検出しその時の周波数を算出する極大点周波数検出回路と、上記極大点周波数検出回路で検出された周波数の最大値を基準点周波数とする最大周波数検出回路と、上記最大周波数検出回路で検出された基準点周波数の時間方向に対する軌跡を平滑化する平滑化回路と、上記平滑化回路で平滑化された基準点周波数の時間方向に対する軌跡から位相補正量を算出する位相補正量算出回路と、上記位相補正量算出回路で算出された位相補正量を用いて上記バッファ回路から出力された目標受信信号の位相を補正する位相補正回路とによって構成したことを特徴とするレーダ信号処理装置。
  5. レーダ装置から入力された目標受信信号をパルス圧縮するパルス圧縮部と、上記パルス圧縮部でパルス圧縮された目標受信信号の時間による距離ずれを補正する距離補正部と、上記距離補正部で補正された目標受信信号の時間による位相ずれを補正する位相補正部と、上記位相補正部で補正された目標受信信号のドップラ周波数を分離する周波数分析部と、上記周波数分析部で周波数分析された目標受信信号の周波数スペクトルを画像データに変換する検波部とを備えたレーダ信号処理装置において、上記位相補正部を上記距離補正部で補正された目標受信信号及びレーダと目標重心との初期距離を格納するバッファ回路と、上記バッファ回路より出力されたレーダと目標重心との初期距離における目標受信信号を時間方向に小区間で周波数分析する区分周波数分析回路と、上記区分周波数分析回路で得られた周波数と振幅の波形について周波数方向における信号の振幅変動のうち変動周期が小さい成分を除去する低域通過フィルタと、上記低域通過フィルタにより上記変動周期が小さい成分が除去された周波数と振幅の波形について閾値を越える周波数領域のデータを切り出すデータ切り出し回路と、上記データ切り出し回路でデータを切り出す際に必要な上記閾値を設定する閾値設定回路と、上記データ切り出し回路で切り出された周波数と振幅の波形について全ての極大となる点を検出しその時の周波数を算出する極大点周波数検出回路と、上記極大点周波数検出回路で検出された周波数の中央値を基準点周波数とする中央周波数検出回路と、上記中央周波数検出回路で検出された基準点周波数の時間方向に対する軌跡を平滑化する平滑化回路と、上記平滑化回路で平滑化された基準点周波数の時間方向に対する軌跡から位相補正量を算出する位相補正量算出回路と、上記位相補正量算出回路で算出された位相補正量を用いて上記バッファ回路から出力された目標受信信号の位相を補正する位相補正回路とによって構成したことを特徴とするレーダ信号処理装置。
JP21371699A 1999-07-28 1999-07-28 レーダ信号処理装置 Expired - Fee Related JP3788117B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP21371699A JP3788117B2 (ja) 1999-07-28 1999-07-28 レーダ信号処理装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP21371699A JP3788117B2 (ja) 1999-07-28 1999-07-28 レーダ信号処理装置

Publications (2)

Publication Number Publication Date
JP2001042031A JP2001042031A (ja) 2001-02-16
JP3788117B2 true JP3788117B2 (ja) 2006-06-21

Family

ID=16643820

Family Applications (1)

Application Number Title Priority Date Filing Date
JP21371699A Expired - Fee Related JP3788117B2 (ja) 1999-07-28 1999-07-28 レーダ信号処理装置

Country Status (1)

Country Link
JP (1) JP3788117B2 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100950301B1 (ko) * 2006-06-02 2010-03-31 (주)연엔지니어링 레이더형 파랑측정 시스템에서의 유의파고 보정방법
JP2009156807A (ja) * 2007-12-27 2009-07-16 Toshiba Corp 測角装置
CN119780841B (zh) * 2024-12-26 2025-09-23 扬州宇安电子科技股份有限公司 一种基于fpga的雷达脉冲信号采样存储方法及装置

Also Published As

Publication number Publication date
JP2001042031A (ja) 2001-02-16

Similar Documents

Publication Publication Date Title
JP6031267B2 (ja) 干渉検出装置、干渉除去器、レーダ装置、干渉検出方法および干渉検出用プログラム
JP2006510236A (ja) ディジタル濾波方法
CN1812488B (zh) 估计视频信号中的水平同步的数字视频处理系统和方法
JP3788117B2 (ja) レーダ信号処理装置
US8433739B2 (en) Methods and systems for detecting repetitive synchronized signal events
JP3788079B2 (ja) レーダ信号処理装置
JP3747856B2 (ja) レーダ信号処理装置
JP4635628B2 (ja) レーダ信号処理装置及びレーダ信号処理方法
JP3788157B2 (ja) レーダ信号処理装置
US8438203B2 (en) Methods and systems for processing and displaying data
US20210059562A1 (en) Dynamic vital-sign detection system and method
EP3053134B1 (en) Method for noise reduction in an image sequence
JP2822827B2 (ja) レーダ信号処理装置
JP2002116251A (ja) レーダ信号処理装置
US7324157B2 (en) Repeat field detecting apparatus, video progressive conversion reproducing apparatus, repeat field detecting method, program, and recording medium
JP3223697B2 (ja) レーダ信号処理装置
JP4063069B2 (ja) レーダ信号処理装置
JP3296234B2 (ja) クラッタ抑圧装置およびクラッタ抑圧方法
JP2020085591A (ja) レーダ信号処理装置及びレーダ信号処理プログラム
JP3791442B2 (ja) レーダ信号処理装置
JP2004061332A (ja) レーダ信号処理装置
JPH10206524A (ja) 目標検出装置
JP3770176B2 (ja) レーダ信号処理装置
JP2004101364A (ja) レーダ装置
JPH11258342A (ja) レーダ信号処理器

Legal Events

Date Code Title Description
RD01 Notification of change of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7421

Effective date: 20040624

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040825

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20051125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20051129

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060124

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060307

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060320

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100407

Year of fee payment: 4

LAPS Cancellation because of no payment of annual fees