JP4545683B2 - Interference wave removing apparatus and interference wave removing method - Google Patents
Interference wave removing apparatus and interference wave removing method Download PDFInfo
- Publication number
- JP4545683B2 JP4545683B2 JP2005354522A JP2005354522A JP4545683B2 JP 4545683 B2 JP4545683 B2 JP 4545683B2 JP 2005354522 A JP2005354522 A JP 2005354522A JP 2005354522 A JP2005354522 A JP 2005354522A JP 4545683 B2 JP4545683 B2 JP 4545683B2
- Authority
- JP
- Japan
- Prior art keywords
- estimated
- delay time
- signal
- memory
- unknown system
- 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
Links
Images
Description
この発明は、予め広帯域信号を送信し、未知系(未知伝播系)のインパルス応答値が局在する遅延時間を推定する複素タップ選択適応フィルタによる干渉波除去装置及び干渉波除去方法に関するものである。 The present invention relates to an interference wave removal apparatus and an interference wave removal method using a complex tap selection adaptive filter that transmits a broadband signal in advance and estimates a delay time in which an impulse response value of an unknown system (an unknown propagation system) is localized. .
適応フィルタによる未知系のシステム同定を応用した例として、干渉波除去装置、エコーキャンセラ、適応等化器などが知られている。ここでは干渉波除去装置の一手法として用いられるタップ選択適応フィルタについて説明する(例えば、特許文献1参照)。タップ選択適応フィルタとは、時間軸上に局在するインパルス応答をもつような未知系の同定を行う場合、インパルス応答が局在する遅延時間のみにおいてタップ数M段のサブ適応フィルタをそれぞれ動作させるものである。 As examples of applying unknown system identification by an adaptive filter, an interference wave canceller, an echo canceller, an adaptive equalizer, and the like are known. Here, a tap selection adaptive filter used as one method of an interference wave removing apparatus will be described (for example, see Patent Document 1). The tap selection adaptive filter, when identifying an unknown system having an impulse response localized on the time axis, operates the sub-adaptive filter having M taps only in the delay time where the impulse response is localized. Is.
このタップ選択適応フィルタは、少タップ段数のサブ適応フィルタを複数個カスケードに配置したもので構成される。 This tap selection adaptive filter is composed of a plurality of sub-adaptive filters having a small number of tap stages arranged in cascade.
未知系の同定の処理では、まず、インパルス応答が局在する時間を推定する。このときには荷重をある複素数値とその位相を変化させた荷重の2種類を用いて遅延時間評価関数を計算して遅延時間を求める。ここでの遅延時間はサンプル間隔で規格化したものである。次いで、インパルス応答値に相当する荷重を求める。荷重はそれぞれ接続された遅延時間推定装置において推定された遅延時間において、タップ数M段の適応フィルタを動作させて求める。上記の遅延時間推定と荷重推定の処理を第1段のサブ適応フィルタ群から第K段のサブ適応フィルタ群まで順次行う。このような構成により未知系のシステム同定を行う。 In the unknown system identification process, first, the time for which the impulse response is localized is estimated. At this time, a delay time evaluation function is calculated using two types of loads, ie, a complex value and a load whose phase is changed, to obtain a delay time. The delay time here is normalized by the sample interval. Next, a load corresponding to the impulse response value is obtained. The load is obtained by operating an adaptive filter with M taps in the delay time estimated by the connected delay time estimation devices. The delay time estimation and load estimation processes are sequentially performed from the first-stage sub-adaptive filter group to the K-th-stage sub-adaptive filter group. The system identification of an unknown system is performed with such a configuration.
未知系の出力が時間軸において局在したインパルス応答である場合、従来の適応フィルタのタップ段数は時間軸においてインパルス応答波形に見合うだけ必要となる。干渉波が広帯域化したり時間軸の広がりが大きくなるなどして適応フィルタのタップ段数が増えてくると演算負荷が大きくなり処理に時間がかかる。そこで、インパルス応答が局在する遅延時間においてのみ適応フィルタを動作させるタップ選択適応フィルタが演算負荷低減の面などから有利であると考えられている。 When the output of the unknown system is an impulse response localized on the time axis, the number of tap stages of the conventional adaptive filter is required only to match the impulse response waveform on the time axis. If the number of tap stages of the adaptive filter increases due to the wide band of interference waves or the spread of the time axis, etc., the computation load increases and processing takes time. Therefore, it is considered that a tap selection adaptive filter that operates the adaptive filter only in the delay time in which the impulse response is localized is advantageous from the viewpoint of reducing the calculation load.
送信信号が狭帯域で周期性をもつ場合、タップ選択適応フィルタの遅延時間評価関数も周期性をもつために、遅延時間の推定を誤り、未知系のシステム同定の精度が劣化してしまう現象が見られた。この原因を遅延時間評価関数の理論にしたがい、図11、図12に基づいて説明する。説明を簡単にするためプラントノイズは十分小さく、かつ未知系のインパルス応答値が局在する遅延時間の数とサブ適応フィルタ群数は一致している状況を想定する。 When the transmission signal has a narrow band and periodicity, the delay time evaluation function of the tap selection adaptive filter also has periodicity, so that the delay time estimation error may occur and the accuracy of system identification of unknown systems may deteriorate. It was seen. This cause will be described with reference to FIGS. 11 and 12 according to the theory of the delay time evaluation function. In order to simplify the explanation, it is assumed that the plant noise is sufficiently small and the number of delay times in which the impulse response values of the unknown system are localized matches the number of sub-adaptive filter groups.
遅延時間評価関数では以下の処理を行う。ここでは、第k段の遅延時間を推定する場合について述べる。図11において、送信信号x(t)を遅延回路90で遅延させた信号に適当な複素数値を持つ荷重g、g*を乗算器91、92でそれぞれ乗算し、第k−1段の誤差信号ek−1(t)から減算器93、94でそれぞれ減算する。これらを二乗平均値算出装置95、96において二乗平均値をそれぞれ求める(図12のステップ901、902)。なお、荷重g*は、荷重gの複素共役をとったものである。ここまでの操作を遅延時間の最大探索値dを未知系のインパルス応答の時間軸を満たすように与え、遅延回路90において遅延時間をd種類変化させ、2つの二乗平均値を求める。遅延時間毎にこれら2つの二乗平均値を加算器97で加算し二乗平均値評価回路98に入力する(図12のステップ903)。
The delay time evaluation function performs the following processing. Here, a case where the delay time of the k-th stage is estimated will be described. In FIG. 11, the signals obtained by delaying the transmission signal x (t) by the
加算した二乗平均値とその一次近似直線間の偏差が最大になる遅延時間を求め、それを遅延時間dkとして出力する。以上の操作を行うことでインパルス応答値が局在する遅延時間dkを推定する。 Deviation between adding the mean square value and its primary approximate line calculated delay time becomes maximum, and outputs it as the delay time d k. By performing the above operation, the delay time d k where the impulse response value is localized is estimated.
この二乗平均誤差を求める式を式(1)に示す。この式(1)が遅延時間dkを決定する原理となる2つの誤差信号の二乗平均値を加算した遅延時間評価関数を示す式である。 An equation for obtaining this mean square error is shown in Equation (1). The equation (1) is an expression showing the delay time evaluation function obtained by adding the mean square value of the two error signal as a principle for determining the delay time d k.
ただし、k=1のときは、誤差信号ek−1(t)は未知系の出力f(t)とする。ここで、dは遅延回路90において探索する遅延時間、ek(t)は誤差信号、tはディジタル信号の時間因子、gは適応フィルタの荷重、g*は適応フィルタの荷重gと複素共役な値を持つ適応フィルタの荷重、N’は適応フィルタの荷重gを用いるサンプル数、N−N’は適応フィルタの荷重g*を用いるサンプル数、xは送信信号を表す。
However, when k = 1, the error signal e k−1 (t) is the output f (t) of the unknown system. Here, d is a delay time searched in the
次に、求められた誤差信号の二乗平均値の一次近似直線を求める。この直線と誤差信号の二乗平均値の偏差が最大値となる探索遅延時間dを遅延時間(^)dk(なお、(^)dはdの上に^があることを表す。)とする。この遅延時間推定方法を模式的に表したものを図13に示す。この図13において、実線は誤差信号の二乗平均値を、点線はその一次近似直線をそれぞれ表している。横軸には遅延時間をとったものである。図中矢印で示した遅延時間において二乗平均値とその一次近似した直線間の偏差が最大値をとり、この遅延時間を第k段の遅延時間とする。偏差の最大値によって遅延時間が推定できる原理は、式(1)を変形し次式で表すことができる。 Next, a first-order approximate line of the mean square value of the obtained error signal is obtained. The search delay time d at which the deviation between the mean square value of the straight line and the error signal becomes the maximum value is the delay time (^) d k (where (^) d indicates that ^ is on d). . FIG. 13 schematically shows this delay time estimation method. In FIG. 13, the solid line represents the mean square value of the error signal, and the dotted line represents the linear approximation line. The horizontal axis shows the delay time. In the delay time indicated by the arrow in the figure, the deviation between the mean square value and the linearly approximated straight line takes the maximum value, and this delay time is set as the k-th delay time. The principle that the delay time can be estimated by the maximum value of the deviation can be expressed by the following equation by modifying the equation (1).
ここで、式(2)の右辺第1項は、次式の通りである。 Here, the first term on the right side of Equation (2) is as follows.
ここで、ckは未知系のインパルス応答値、rxxはサンプル数N’個の送信信号の自己相関関数、r’xxはサンプル数N−N’個の送信信号の自己相関関数、rnx、rxnは送信信号とプラントノイズの相互相関関数、Kはサブ適応フィルタ群数、*は複素共役を表す。 Here, ck is an impulse response value of an unknown system, r xx is an autocorrelation function of a transmission signal having N ′ samples, r ′ xx is an autocorrelation function of a transmission signal having N−N ′ samples, and r nx , R xn is a cross-correlation function between a transmission signal and plant noise, K is the number of sub-adaptive filter groups, and * is a complex conjugate.
式(2)において、右辺第1項は、定数値となる。なぜならば、式(3)より、送信信号に白色雑音を考えているために送信信号の自己相関値は0となり、送信信号とプラントノイズの相互相関値も0となり、残りの項は定数値をとるからである。 In Expression (2), the first term on the right side is a constant value. This is because, from equation (3), since white noise is considered in the transmission signal, the autocorrelation value of the transmission signal is 0, the cross-correlation value of the transmission signal and plant noise is also 0, and the remaining terms are constant values. Because it takes.
式(2)の右辺第2項と第3項は、dをパラメータとして遅延時間探索中に遅延時間がd=dkになるとき、rxx(d−dk)は最大値をもち、一次近似直線との偏差が最大になる。このときの探索遅延時間dを(^)dkとすることでインパルス応答値が局在する遅延時間dkを推定することができる(図12のステップ904)。
The second term and the third term on the right side of Equation (2) are such that when the delay time becomes d = d k during the delay time search using d as a parameter, r xx (d−d k ) has a maximum value, The deviation from the approximate line is maximized. Search delay time d at this time (^) impulse response value by the d k can be used to estimate the delay time d k localized (
しかし、これまで説明してきた複素タップ選択適応フィルタを用いた干渉波除去装置では、送信信号が狭帯域で周期性をもつ場合、遅延時間評価関数も周期性をもち真値以外の遅延時間を推定値とすることがみられた。この誤った遅延時間推定値のため、未知伝播系のインパルス応答値の推定も不可能になり、抑圧性能が劣化することが見られた。 However, in the interference wave canceller using the complex tap selection adaptive filter described so far, when the transmission signal has a narrow band and periodicity, the delay time evaluation function also has periodicity and estimates a non-true delay time. Value was seen. Due to this incorrect delay time estimate, it was impossible to estimate the impulse response value of the unknown propagation system, and the suppression performance was degraded.
送信信号が無変調の正弦波のように周期性をもつ信号の場合に、遅延時間推定の評価関数Φ’(dk)の振幅が振動する問題点を説明する。 The problem that the amplitude of the evaluation function Φ ′ (d k ) for estimating the delay time oscillates when the transmission signal is a signal having periodicity such as an unmodulated sine wave will be described.
式(2)に示す遅延時間評価関数Φ’(d)の右辺第2項と第3項に注目すると、d=dkとなるとき、rxxとr’xxは最大値をとる。 When attention is paid to the second term and the third term on the right side of the delay time evaluation function Φ ′ (d) shown in Expression (2), r xx and r ′ xx take the maximum values when d = d k .
以下、ここで送信信号が複素の正弦波である場合を考える。簡単のためサブ適応フィルタ群数はK=1、未知系のインパルス応答もc1のみ1種類とする。また、自己相関関数rxxとr’xxはサンプル数が異なるが同一と考えられる。これらを考慮して式(2)の右辺第2項と第3項を変形すると、以下のようになる。 Hereinafter, a case where the transmission signal is a complex sine wave will be considered. Sub adaptive filter number set for the sake of simplicity K = 1, the impulse response also c 1 only one type of the unknown system. The autocorrelation functions r xx and r ′ xx are considered to be the same although the number of samples is different. Taking these into account, the second term and the third term on the right side of Equation (2) are transformed as follows.
ここで、送信信号が複素の正弦波e−j2πftとすると、式(4)は、以下のようになる。 Here, if the transmission signal is a complex sine wave e −j2πft , Equation (4) is as follows.
ここで、〈〉内は、期待値演算を表す。式(5)において、dの探索と無関係にe−j2πf(d−d1)により周期性が生じるため、遅延時間探索中にΦ’(d)の振幅が振動することがわかる。なお、e−j2πf(d−d1)における上付きの(d−d1)は、(d−d1)である。これらのことから、送信信号が周期関数の場合、評価関数が周期性をもつため、遅延時間の推定を誤ってしまう恐れがあることが分かる。このときの遅延時間推定方法を、模式的に表したものを図14に示す。 Here, the inside of <> represents an expected value calculation. In formula (5), since periodicity is generated by e −j2πf (d−d1) regardless of the search for d, it can be seen that the amplitude of Φ ′ (d) oscillates during the delay time search. The superscript (d−d1) in e −j2πf (d− d1) is (d−d 1 ). From these facts, it can be understood that when the transmission signal is a periodic function, the evaluation function has periodicity, and therefore the estimation of the delay time may be erroneous. FIG. 14 schematically shows the delay time estimation method at this time.
この発明は、上述のような課題を解決するためになされたもので、その目的は、送信信号が狭帯域で周期性をもつ場合に、低演算量である複素タップ選択適応フィルタを用いるために、狭帯域信号の信号処理に対応させ、未知系からの干渉波を除去することができる干渉波除去装置及び干渉波除去方法を得るものである。 The present invention has been made to solve the above-described problems, and its purpose is to use a complex tap selection adaptive filter having a low amount of computation when a transmission signal has a narrow band and periodicity. Thus, an interference wave removing apparatus and an interference wave removing method capable of removing an interference wave from an unknown system in correspondence with signal processing of a narrowband signal are obtained.
この発明に係る干渉波除去装置は、広帯域信号を発生する広帯域信号発生器と、送信信号を発生する送信信号発生器と、前記広帯域信号を未知系へ送信するために前記広帯域信号発生器側に切り替えられると第1の制御信号を出力し、前記送信信号を前記未知系へ送信するために前記送信信号発生器側に切り替えられると第2の制御信号を出力するスイッチと、推定される遅延時間を記憶するためのメモリと、前記第1の制御信号が入力されると、遅延時間評価関数に基づき遅延時間を推定して推定遅延時間を前記メモリに記憶するとともに、前記第2の制御信号が入力されると、前記メモリから推定遅延時間を呼び出し適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求める複素タップ選択適応フィルタと、前記送信信号が前記未知系を通過した受信信号から前記複素タップ選択適応フィルタの出力を減算することにより干渉波を抑圧する減算器とを設けたものである。 An interference wave canceller according to the present invention includes a wideband signal generator that generates a wideband signal, a transmission signal generator that generates a transmission signal, and a broadband signal generator that transmits the wideband signal to an unknown system. A switch that outputs a first control signal when switched, and outputs a second control signal when switched to the transmission signal generator in order to transmit the transmission signal to the unknown system, and an estimated delay time When the first control signal is input to the memory for storing the delay time, the delay time is estimated based on the delay time evaluation function, the estimated delay time is stored in the memory, and the second control signal is When input, the estimated delay time is called from the memory, a weight at the estimated delay time is estimated using an adaptive algorithm, and an error signal is obtained from the estimated delay time and the estimated weight. And flop selection adaptive filter, in which the transmission signal is provided a subtracter reduce interference signals by subtracting the output of the complex tap selection adaptive filter from the received signal that has passed through the unknown system.
この発明に係る干渉波除去方法は、広帯域信号を未知系へ送信し、前記広帯域信号の遅延時間評価関数に基づき遅延時間を推定して推定遅延時間をメモリに記憶するステップと、送信信号を前記未知系へ送信し、前記メモリから推定遅延時間を呼び出し適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求めるステップと、前記送信信号が前記未知系を通過した受信信号から前記誤差信号を減算することにより干渉波を抑圧するステップとを設けたものである。 The interference wave elimination method according to the present invention includes a step of transmitting a wideband signal to an unknown system, estimating a delay time based on a delay time evaluation function of the wideband signal, and storing the estimated delay time in a memory; Transmitting to the unknown system, calling the estimated delay time from the memory, estimating a load at the estimated delay time using an adaptive algorithm, obtaining an error signal from the estimated delay time and the estimated weight, and transmitting the transmitted signal to the unknown And a step of subtracting the error signal from the received signal that has passed through the system to suppress the interference wave.
この発明に係る干渉波除去装置及び方法は、送信信号が狭帯域で周期性をもつ場合に、低演算量である複素タップ選択適応フィルタを用いるために、狭帯域信号の信号処理に対応させ、未知系からの干渉波を除去することができるという効果を奏する。 The interference wave canceling apparatus and method according to the present invention is adapted to signal processing of a narrowband signal in order to use a complex tap selection adaptive filter having a low calculation amount when the transmission signal has a periodicity in a narrowband, There is an effect that interference waves from unknown systems can be removed.
実施の形態1.
この発明の実施の形態1に係る干渉波除去装置について図1から図3までを参照しながら説明する。図1は、この発明の実施の形態1に係る干渉波除去装置の構成を示す図である。また、図2は、この発明の実施の形態1に係る干渉波除去装置の複素タップ選択適応フィルタの構成を示す図である。なお、各図中、同一符号は同一又は相当部分を示す。
An interference wave canceling apparatus according to
図1において、この発明の実施の形態1に係る干渉波除去装置は、広帯域信号を発生する広帯域信号発生器1と、送信信号を発生する送信信号発生器2と、広帯域信号発生器1と送信信号発生器2を切り替えるスイッチ3と、広帯域信号や送信信号x(t)を未知系(unknown system)へ送信する送信アンテナ4と、未知系のインパルス応答が局在する場合にその遅延時間を推定し、その遅延時間に当たるタップのみ適応動作させ、複素数化を図った複素タップ選択適応フィルタ5と、推定した遅延時間などを記憶するメモリ6と、
所望波を受信する受信アンテナ7と、受信信号f(t)とプラントノイズn(t)を加算する加算器8と、加算器8の出力から適応フィルタ5の出力を減算して所望波ek(t)を出力する減算器9(91,…,9k,…,9K)とが設けられている。なお、スイッチ3と送信アンテナ4の間には、変調器などの送信手段の図示が省略され、受信アンテナ7と加算器8の間には、復調器などの受信手段の図示が省略されている。また、送信アンテナ4から受信アンテナ7への矢印は、漏れこみ干渉波すなわち未知系を経由してきた送信信号を表す。
In FIG. 1, an interference wave canceling apparatus according to
A receiving antenna 7 for receiving a desired wave, an adder 8 for adding the received signal f (t) and plant noise n (t), desired wave e k by subtracting the output of the
また、図2において、複素タップ選択適応フィルタ5は、遅延時間推定装置511,…,51k,…,51Kと、M段の適応フィルタ521,…,52k,…,52Kとが設けられている。なお、図1及び図2において、点線は制御信号を表す。
Further, in FIG. 2, the complex tap selection
つぎに、この実施の形態1に係る干渉波除去装置の動作について図面を参照しながら説明する。図3は、この発明の実施の形態1に係る干渉波除去装置の動作を示すフローチャートである。
Next, the operation of the interference wave canceller according to
ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態1に係る干渉波除去装置は、サブ適応フィルタ群数がK=1である場合で、予め広帯域信号で遅延時間を推定しておき、荷重はそれぞれの干渉波が入射するときに求める。 Here, for the sake of simplicity, the description will be made based on a situation where the plant noise is sufficiently small. In the interference wave canceling apparatus according to the first embodiment, when the number of sub-adaptive filter groups is K = 1, the delay time is estimated in advance with a broadband signal, and the load is applied when each interference wave is incident. Ask.
予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置511、M段の適応フィルタ521、メモリ6が区別可能なように第1の制御信号(点線)を送る。複素タップ選択適応フィルタ5の動作を図3に示す。ここでは、図3中の「予め広帯域信号で行う処理」を行う。
Consider a case where a broadband signal (white noise or the like) is transmitted from the
ステップ101、102、103において、広帯域信号を送信した場合、式(2)中の送信信号の自己相関関数rxxはd=diのみに相関値のピークが現れる。この場合、遅延時間評価関数Φ’(d)は唯一のピークをもつ。
When a wideband signal is transmitted in
次に、ステップ104において、この唯一のピークをもつ探索パラメータd1を遅延時間推定値(^)d1−preとする。(なお、下付きの1−preは1_preを表す。)この遅延時間推定値(^)d1−preをメモリ6に記憶する。ここまでで前処理は終了する。
Next, in
実際に干渉波を抑圧処理する場合には、図3中の「干渉波抑圧時の処理」を行う。このとき、スイッチ3は送信信号発生器2側に切り替えて、遅延時間推定装置511、M段の適応フィルタ521、メモリ6へ第2の制御信号(点線)を送る。
When the interference wave is actually suppressed, the “process when suppressing interference wave” in FIG. 3 is performed. At this time, the switch 3 switches to the
ステップ151では、複素タップ選択適応フィルタ5の遅延時間推定の処理を行わず、予め広帯域信号を送信して求めた推定遅延時間(^)d1−preをメモリ6から呼び出し、直接前処理で求めた値を遅延時間推定値(^)d1=(^)d1−preとする。この遅延時間(^)d1=(^)d1−preにおいて、適応フィルタ521を動作させ、適応アルゴリズムを用いて、その荷重g1を計算し、計算後のg1を荷重推定値(^)g1とする。適応アルゴリズムの種類については、正規化LMSアルゴリズムなど特に問わない。
In
次に、ステップ152では、ステップ151で用いた遅延時間(^)d1=(^)d1−preと求められた荷重g1を用いて、誤差信号e1(t)を求める。この行程は、次式で表すことができる。
Next, in
ここで、Hは複素共役転置を示す。 Here, H indicates complex conjugate transpose.
これらの処理が終わった時点で、未知系のインパルス応答が局在する遅延時間と、その大きさの推定が行なえる(ステップ153)。 When these processes are completed, the delay time in which the impulse response of the unknown system is localized and the magnitude thereof can be estimated (step 153).
推定された未知系のインパルス応答の遅延時間とその大きさから未知系を模擬し、受信信号から模擬した未知系出力信号を減じて干渉波の抑圧を行う(ステップ154)。 The unknown system is simulated from the estimated delay time of the impulse response of the unknown system and its magnitude, and the unknown system output signal is subtracted from the received signal to suppress the interference wave (step 154).
すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置511、M段の適応フィルタ521、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、適応フィルタ521は、メモリ6から荷重g、g*を呼び出して設定し、遅延時間推定装置511は、荷重g、g*を用いてそれぞれ誤差信号の二乗平均値を求め、2つの二乗平均値から遅延時間評価関数を求め、この遅延時間評価関数から遅延時間を推定して推定遅延時間をメモリ6に記憶する。これらの前処理がステップ101〜104である。
That is, when the switch 3 is switched to the
実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置511、M段の適応フィルタ521、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置511は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ521は、適応アルゴリズムを用いて荷重g1を求め、減算器91は、推定遅延時間と荷重g1から誤差信号を求め、これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ151〜154である。
When the interference wave is actually suppressed, the switch 3 is switched to the
実施の形態2.
この発明の実施の形態2に係る干渉波除去装置について図4を参照しながら説明する。図4は、この発明の実施の形態2に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態2に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
An interference wave removing apparatus according to
ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態2では、サブ適応フィルタ群数がK≧2である場合で、予め広帯域信号を送信し、1,…,K段の遅延時間推定値とそれぞれのサブ適応フィルタの荷重求め、これらのうち1,…,K段の遅延時間推定値を実際の干渉波抑圧時に用いるものである。 Here, for the sake of simplicity, the description will be made based on a situation where the plant noise is sufficiently small. In the second embodiment, when the number of sub-adaptive filter groups is K ≧ 2, a wideband signal is transmitted in advance, 1,..., K-stage delay time estimation values and respective sub-adaptive filter weights are obtained. Of these, the estimated delay time values of 1,..., K stages are used for actual interference wave suppression.
予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は、広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51k、タップ数M段の適応フィルタ52k、メモリ6が区別可能なように第1の制御信号(点線)を送る。複素タップ選択適応フィルタ5の動作を図4に示す。ここでは、図4中の「予め広帯域信号で行う処理」を行う。
Consider a case where a broadband signal (white noise or the like) is transmitted from the
ステップ201、202、203において、広帯域信号を送信した場合、式(2)中の送信信号の自己相関関数rxxはd=dkのみに相関値のピークが現れる。この場合、遅延時間評価関数Φ’(d)は唯一のピークをもつ。
In
次に、ステップ204において、この唯一のピークをもつ探索パラメータdiを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。
Next, in
次に、ステップ205では、ステップ204で得られた遅延時間推定値で少タップ段数の適応フィルタ52kを動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとする。
Next, in step 205, operates the adaptive filter 52 k of small number of tap stages in the delay time estimation value obtained in
次に、ステップ206では、遅延時間推定値(^)dk−preと荷重(^)gk−preより誤差信号ek−pre(t)を求める。
Next, in
そして、ステップ207では、ステップ206で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、この繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまで、ステップ201〜206を繰り返す。ここまでで前処理は終了する。
In
実際に干渉波を抑圧処理する場合には、図4中の「干渉波抑圧時の処理」を行う。このとき、スイッチ3は、送信信号発生器2側に切り替えて、遅延時間推定装置51k、M段の適応フィルタ52k、メモリ6へ第2の制御信号(点線)を送る。
When the interference wave is actually suppressed, the “process when suppressing interference wave” in FIG. 4 is performed. At this time, the switch 3 switches to the
ステップ251では、複素タップ選択適応フィルタ5の遅延時間推定の処理を行わず、予め広帯域信号を送信して求めた推定遅延時間(^)dk−preをメモリ6から呼び出し、前処理で求めた(^)dk−preを第k段の遅延時間推定値(^)dk=(^)dk−preとする。この遅延時間(^)dk=(^)dk−preにおいて、タップ数M段の適応フィルタ52kを動作させ、計算後の荷重値gkを荷重推定値(^)gkとし、この荷重値を未知系のインパルス応答推定値とする。適応アルゴリズムの種類については問わない。
In
次に、ステップ252では、ステップ251で用いた遅延時間(^)dk=(^)dk−preと求められた荷重(^)gkを用いて、誤差信号ek(t)を求める。この行程は次式で表すことができる。
Next, in
ただし、k=1のときは、ek−1(t)は未知系の出力f(t)とする。 However, when k = 1, e k−1 (t) is the output f (t) of the unknown system.
そして、ステップ253において、ステップ252で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまで、ステップ251、252の処理を繰り返す。
In
これらの処理が終わった時点で全ての未知系のインパルス応答が局在する遅延時間と、その大きさの推定が行なえる。 When these processes are completed, it is possible to estimate the delay time in which the impulse responses of all unknown systems are localized and the magnitude thereof.
推定された全ての未知系のインパルス応答の遅延時間とその大きさから未知系を模擬し、干渉波の除去を行う(ステップ254)。 The unknown system is simulated from the estimated delay times and the magnitudes of the impulse responses of all unknown systems, and interference waves are removed (step 254).
すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、適応フィルタ52kは、メモリ6から荷重g、g*を呼び出して設定し、遅延時間推定装置51kは、荷重g、g*を用いてそれぞれ誤差信号の二乗平均値を求め、2つの二乗平均値から遅延時間評価関数を求め、この遅延時間評価関数から遅延時間を推定して推定遅延時間をメモリ6に記憶する。また、遅延時間推定装置51kは、適応アルゴリズムを用いて推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の二乗平均値を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ201〜207である。
That is, when the switch 3 is switched to the
実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51kは、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52kは、適応アルゴリズムを用いて荷重(^)gkを求めて、減算器9kは、推定遅延時間と推定荷重(^)gkから誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の荷重を求める処理から最後の誤差信号を求める処理までを繰り返す。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ251〜254である。
When the interference wave is actually suppressed, the switch 3 is switched to the
実施の形態3.
この発明の実施の形態3に係る干渉波除去装置について図5を参照しながら説明する。図5は、この発明の実施の形態3に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態3に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
Embodiment 3 FIG.
An interference wave removing apparatus according to Embodiment 3 of the present invention will be described with reference to FIG. FIG. 5 is a flowchart showing the operation of the interference wave canceller according to Embodiment 3 of the present invention. Note that the configuration of the interference wave canceller according to Embodiment 3 of the present invention is the same as that of
ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態3では、サブ適応フィルタ群数がK≧1である場合で、予め広帯域信号を送信し、1,…,K段の遅延時間推定値とそれぞれのサブ適応フィルタの荷重を求め、これらを実際の干渉波処理時に用いるものである。 Here, for the sake of simplicity, the description will be made based on a situation where the plant noise is sufficiently small. In the third embodiment, when the number of sub-adaptive filter groups is K ≧ 1, a wideband signal is transmitted in advance, and a delay time estimated value of 1,..., K stages and a weight of each sub-adaptive filter are obtained. These are used during actual interference wave processing.
予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は、広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51k、タップ数M段の適応フィルタ52k、メモリ6が区別可能なように第1の制御信号(点線)を送る。複素タップ選択適応フィルタ5の動作を図5に示す。ここでは、図5中の「予め広帯域信号で行う処理」を行う。
Consider a case where a broadband signal (white noise or the like) is transmitted from the
ステップ301、302、303において、広帯域信号を送信した場合、式(2)中の送信信号の自己相関関数rxxはdk=diのみに相関値のピークが現れる。この場合、遅延時間評価関数Φ’(d)は唯一のピークをもつ。
In
次に、ステップ304において、この唯一のピークをもつ探索パラメータdiを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。
Next, in
次に、ステップ305では、ステップ304で得られた遅延時間推定値で少タップ段数の適応フィルタ52kを動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとし、これをメモリ6に記憶する。
Next, in
次に、ステップ306では、遅延時間(^)dk−preと荷重(^)gk−preより誤差信号ek−pre(t)を求める。
Next, in
そして、ステップ307において、ステップ306で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまでステップ301〜306を繰り返す。ここまでで前処理は終了する。
In
実際に干渉波を抑圧処理する場合には、図5中「干渉波抑圧時の処理」を行う。このとき、スイッチ3は送信信号発生器2側に切り替えて、遅延時間推定装置51k、M段の適応フィルタ52k、メモリ6へ第2の制御信号(点線)を送る。
When the interference wave is actually suppressed, the “processing when suppressing the interference wave” in FIG. 5 is performed. At this time, the switch 3 switches to the
ステップ351では、複素タップ選択適応フィルタ5の遅延時間推定の処理を行わず、予め広帯域信号を送信して求めた推定遅延時間(^)dk−preと荷重(^)gk−preをメモリ6から呼び出し、前処理で求めた(^)dk−preを第k段の遅延時間推定値(^)dk=(^)dk−preとする。この遅延時間(^)dk=(^)dk−preにおいて、タップ数M段の適応フィルタ52kを動作させ、その荷重gkを計算し、計算後の荷重値を(^)gkとし未知系のインパルス応答推定値とする。ここでは、荷重更新の初期値として推定荷重(^)gk−preを用いる。これにより、適応アルゴリズムの収束時間の短縮が可能になる。適応アルゴリズムの種類については問わない。
In
次に、ステップ352では、ステップ351で用いた遅延時間(^)dk=(^)dk−preと求められた荷重(^)gkを用いて、誤差信号ek(t)を求める。この行程は次式で表すことができる。
Next, in
ただし、k=1のときは、ek−1(t)は未知系の出力f(t)とする。 However, when k = 1, e k−1 (t) is the output f (t) of the unknown system.
そして、ステップ353において、ステップ352で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまでステップ351、352の処理を繰り返す。
In
これらの処理が終わった時点で全ての未知系のインパルス応答が局在する遅延時間とその大きさの推定が行なえる。 When these processes are completed, it is possible to estimate the delay time and the size of the impulse responses of all unknown systems.
推定された全ての未知系のインパルス応答の遅延時間とその大きさから未知系を模擬し、干渉波の除去を行う(ステップ354)。 The unknown system is simulated from the estimated delay times and the magnitudes of the impulse responses of all unknown systems, and interference waves are removed (step 354).
すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、適応フィルタ52kは、メモリ6から荷重g、g*を呼び出して設定し、遅延時間推定装置51kは、荷重g、g*を用いてそれぞれ誤差信号の二乗平均値を求め、2つの二乗平均値から遅延時間評価関数を求め、この遅延時間評価関数から遅延時間を推定して推定遅延時間をメモリ6に記憶する。また、遅延時間推定装置51kは、適応アルゴリズムを用いて推定遅延時間での荷重を推定して推定荷重をメモリ6に記憶する。次に、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の二乗平均値を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ301〜307である。
That is, when the switch 3 is switched to the
実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51kは、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52kは、メモリ6から推定荷重を呼び出して初期値として設定し、適応アルゴリズムを用いて荷重gkを推定し、減算器9kは、推定遅延時間と推定荷重(^)gkから誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の荷重を求める処理から最後の誤差信号を求める処理までを繰り返す。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ351〜354である。
When the interference wave is actually suppressed, the switch 3 is switched to the
実施の形態4.
この発明の実施の形態4に係る干渉波除去装置について図6を参照しながら説明する。図6は、この発明の実施の形態4に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態4に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
Embodiment 4 FIG.
An interference wave removing apparatus according to Embodiment 4 of the present invention will be described with reference to FIG. FIG. 6 is a flowchart showing the operation of the interference wave canceller according to Embodiment 4 of the present invention. Note that the configuration of the interference wave canceller according to Embodiment 4 of the present invention is the same as that of
ここでは簡単のためにプラントノイズは十分小さい状況により説明する。上記の実施の形態3では、予め広帯域信号について処理し、1,…,K段の遅延時間推定値を求め、これらの遅延時間推定値における、1,…,K段の適応フィルタの荷重もそれぞれメモリ6に記憶し、干渉波抑圧時の第k段の荷重計算時の初期値に使用するものであるが、この実施の形態4では、予め広帯域信号の処理時に得られた荷重値を干渉抑圧時にも荷重値として用いるものである。
Here, for the sake of simplicity, the description will be made based on a situation where the plant noise is sufficiently small. In the third embodiment, a wideband signal is processed in advance to obtain 1,..., K stages of delay time estimates, and the load of the 1,. Although stored in the
予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51k、M段の適応フィルタ52k、メモリ6が区別可能なように第1の制御信号(点線)を送る。ここでは、図6中の「予め広帯域信号で行う処理」を行う。
Consider a case where a broadband signal (white noise or the like) is transmitted from the
ステップ401、402、403において、広帯域信号を送信した場合、式(2)中の送信信号の自己相関関数rxxはd=dkのみに相関値のピークが現れる。この場合、遅延時間評価関数Φ’(d)は唯一のピークをもつ。
In
次に、ステップ404において、この唯一のピークをもつ探索パラメータdiを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。
Next, in
次に、ステップ405では、ステップ404で得られた遅延時間推定値で少タップ段数の適応フィルタ52kを動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとし、これをメモリ6に記憶する。
Next, in
次に、ステップ406では、遅延時間(^)dk−preと荷重(^)gk−preより誤差信号ek−pre(t)を求める。
Next, in
そして、ステップ407において、ステップ406で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまでステップ401〜406を繰り返す。ここまでで前処理は終了する。
In
実際に干渉波を抑圧処理する場合には、図6中「干渉波抑圧時の処理」を行う。このとき、スイッチ3は送信信号発生器2側に切り替えて、遅延時間推定装置51k、M段の適応フィルタ52k、メモリ6へ第2の制御信号(点線)を送る。
When the interference wave is actually suppressed, the “processing when suppressing the interference wave” in FIG. 6 is performed. At this time, the switch 3 switches to the
ステップ451では、複素タップ選択適応フィルタ5の遅延時間推定の処理と、タップ数M段の適応フィルタ52kの荷重演算を行わず、予め広帯域信号を送信して求めた推定遅延時間(^)dk−preと推定荷重(^)gk−preをメモリ6から呼び出し、この遅延時間推定値(^)dk−preと荷重値(^)gk−preを用いて、誤差信号ek(t)を求める。これにより、適応アルゴリズムの演算が不要になる。この行程は次式で表すことができる。
In
ただし、k=1のときは、ek−1(t)は未知系の出力f(t)とする。 However, when k = 1, e k−1 (t) is the output f (t) of the unknown system.
これらの処理が終わった時点で全ての未知系のインパルス応答が局在する遅延時間とその大きさの推定が行なえる。 When these processes are completed, it is possible to estimate the delay time and the size of the impulse responses of all unknown systems.
推定された全ての未知系のインパルス応答の遅延時間とその大きさから未知系を模擬し、干渉波の除去を行う(ステップ452)。 The unknown system is simulated from the estimated delay times and the magnitudes of the impulse responses of all unknown systems, and interference waves are removed (step 452).
すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、適応フィルタ52kは、メモリ6から荷重g、g*を呼び出して設定し、遅延時間推定装置51kは、荷重g、g*を用いてそれぞれ誤差信号の二乗平均値を求め、2つの二乗平均値から遅延時間評価関数を求め、この遅延時間評価関数から遅延時間を推定して推定遅延時間をメモリ6に記憶する。また、遅延時間推定装置51kは、適応アルゴリズムを用いて推定遅延時間での荷重を推定して推定荷重をメモリ6に記憶する。次に、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の二乗平均値を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ401〜407である。
That is, when the switch 3 is switched to the
実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51kは、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52kは、メモリ6から推定荷重を呼び出して設定する。減算器9kは、推定遅延時間と推定荷重(^)gkから誤差信号を求める。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ451〜452である。
When the interference wave is actually suppressed, the switch 3 is switched to the
実施の形態5.
この発明の実施の形態5に係る干渉波除去装置について図7を参照しながら説明する。図7は、この発明の実施の形態5に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態5に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
An interference wave removing apparatus according to
ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態5は、サブ適応フィルタ群数がK=1である場合で、予め広帯域信号で遅延時間を推定しておき、荷重はそれぞれの干渉波が入射するときに求める。 Here, for the sake of simplicity, the description will be made based on a situation where the plant noise is sufficiently small. In the fifth embodiment, when the number of sub-adaptive filter groups is K = 1, the delay time is estimated in advance with a wideband signal, and the load is obtained when each interference wave is incident.
予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置511、M段の適応フィルタ521、メモリ6が区別可能なように第1の制御信号(点線)を送る。ここでは、図7中の「予め広帯域信号で行う処理」を行う。
Consider a case where a broadband signal (white noise or the like) is transmitted from the
ステップ501において、広帯域信号を送信した場合、広帯域信号x(t)と、この広帯域信号が未知系を通過した受信信号f(t)との相互相関関数rxfを計算する。これは次式で表すことができる。dは遅延時間探索パラメータとし、未知系のインパルス応答の時間軸を満たすように与える。
In
ここで、〈〉内は期待値演算、*は複素共役をそれぞれ表す。 Here, <> represents an expected value calculation, and * represents a complex conjugate.
次に、ステップ502において、式(10)はd1=diのみに相関値のピークが現れる。この唯一のピークを与える探索パラメータdiを遅延時間推定値(^)d1−preとする。この遅延時間推定値(^)d1−preをメモリ6に記憶する。ここまでで前処理は終了する。
Next, in
以下、ステップ551〜554の処理は、上記の実施の形態1のステップ151〜154とそれぞれ同様であるので説明は省略する。
Hereinafter, the processing in
すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置511、M段の適応フィルタ521、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、遅延時間推定装置511は、広帯域信号x(t)と、この広帯域信号が未知系を通過してきた受信信号f(t)の相互相関関数を計算し、相互相関値が最大となる遅延時間を推定遅延時間としてメモリ6に記憶する。これらの前処理がステップ501〜502である。
That is, when the switch 3 is switched to the
実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置511、M段の適応フィルタ521、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置511は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ521は、適応アルゴリズムを用いて荷重g1を求め、減算器91は、推定遅延時間と荷重g1から誤差信号を求め、これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ551〜554である。
When the interference wave is actually suppressed, the switch 3 is switched to the
実施の形態6.
この発明の実施の形態6に係る干渉波除去装置について図8を参照しながら説明する。図8は、この発明の実施の形態6に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態6に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
An interference wave canceller according to
ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態6では、サブ適応フィルタ群数がK≧2である場合で、予め広帯域信号について処理し、1,…,K段の遅延時間推定値とそれぞれのサブ適応フィルタの荷重求め、これらのうち、1,…,K段の遅延時間推定値を実際の干渉波処理時に用いるものである。 Here, for the sake of simplicity, the description will be made based on a situation where the plant noise is sufficiently small. In the sixth embodiment, when the number of sub-adaptive filter groups is K ≧ 2, wideband signals are processed in advance, 1,..., K-stage delay time estimation values and respective sub-adaptive filter weights are obtained, Of these, the estimated delay time values of 1,..., K stages are used during actual interference wave processing.
予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51k、M段の適応フィルタ52k、メモリ6が区別可能なように第1の制御信号(点線)を送る。この前処理では、図8中の「予め広帯域信号で行う処理」を行う。
Consider a case where a broadband signal (white noise or the like) is transmitted from the
ステップ601において、広帯域信号を送信し、広帯域信号と、この広帯域信号が未知系を通過した受信信号との相互相関関数rxfを計算する。
In
次に、ステップ602において、相互相関関数rxfはd=dkのみに相関値のピークが現れる。この唯一のピークをもつ探索パラメータdkを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。
Next, in
次に、ステップ603では、ステップ602で得られた遅延時間推定値(^)dk−preで少タップ段数の適応フィルタ52kを動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとする。
Next, in
次に、ステップ604では、遅延時間推定値(^)dk−preと荷重(^)gk−preを用いて誤差信号ek−pre(t)を次式により計算する。
Next, in
ただし、k=1のときは、ek−1(t)は未知系の出力f(t)とする。 However, when k = 1, e k−1 (t) is the output f (t) of the unknown system.
そして、ステップ605において、ステップ604で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまでステップ601〜604を繰り返す。ここまでで前処理は終了する。
In
以下、ステップ651〜654の処理は、上記の実施の形態2のステップ251〜254とそれぞれ同様であるので説明は省略する。
Hereinafter, the processing in
すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、遅延時間推定装置51kは、広帯域信号x(t)と、この広帯域信号が未知系を通過してきた受信信号f(t)の相互相関関数を計算し、相互相関値が最大となる遅延時間を推定遅延時間としてメモリ6に記憶する。また、遅延時間推定装置51kは、適応アルゴリズムを用いて推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の相互相関関数を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ601〜605である。
That is, when the switch 3 is switched to the
実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51kは、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52kは、適応アルゴリズムを用いて荷重gkを求め、減算器9kは、推定遅延時間と推定荷重(^)gkから誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の荷重を求める処理から最後の誤差信号を求める処理までを繰り返す。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ651〜654である。
When the interference wave is actually suppressed, the switch 3 is switched to the
実施の形態7.
この発明の実施の形態7に係る干渉波除去装置について図9を参照しながら説明する。図9は、この発明の実施の形態7に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態7に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
Embodiment 7 FIG.
An interference wave canceling apparatus according to Embodiment 7 of the present invention will be described with reference to FIG. FIG. 9 is a flowchart showing the operation of the interference wave canceller according to Embodiment 7 of the present invention. Note that the configuration of the interference wave canceller according to Embodiment 7 of the present invention is the same as that of
ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態7では、サブ適応フィルタ群数がK≧1である場合で、予め広帯域信号について処理し、1,…,K段の遅延時間推定値とそれぞれのサブ適応フィルタの荷重求め、これらを実際の干渉波抑圧時に用いるものである。 Here, for the sake of simplicity, the description will be made based on a situation where the plant noise is sufficiently small. In the seventh embodiment, when the number of sub-adaptive filter groups is K ≧ 1, wideband signals are processed in advance, 1,..., K stages of delay time estimation values and respective sub-adaptive filter weights are obtained, Is used when actual interference waves are suppressed.
予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51k、M段の適応フィルタ52k、メモリ6が区別可能なように第1の制御信号(点線)を送る。前処理では、図9中の「予め広帯域信号で行う処理」を行う。
Consider a case where a broadband signal (white noise or the like) is transmitted from the
ステップ701において、広帯域信号を送信し、広帯域信号と、この広帯域信号が未知系を通過した受信信号の相互相関関数rxfを計算する。
In
次に、ステップ702において、相互相関関数rxfはd=dkのみに相関値のピークが現れる。この唯一のピークをもつ探索パラメータdを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。
Next, in
次に、ステップ703では、ステップ702で得られた遅延時間推定値(^)dk−preで少タップ段数の適応フィルタを動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとし、これをメモリ6に記憶する。
Next, in
次に、ステップ704では、遅延時間推定値(^)dk−preと荷重(^)gk−preを用いて誤差信号ek−pre(t)を次式により計算する。 Next, in step 704, the error signal e k-pre (t) is calculated by the following equation using the delay time estimated value (^) d k-pre and the load (^) g k-pre .
ただし、k=1のとき、ek−1−pre(t)は未知系の出力f(t)とする。 However, when k = 1, ek-1-pre (t) is an unknown output f (t).
そして、ステップ705において、ステップ704で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまでステップ701〜704を繰り返す。ここまでで前処理は終了する。
In
以下、ステップ751〜754の処理は、上記の実施の形態3のステップ351〜354とそれぞれ同様であるので説明は省略する。
Hereinafter, the processing in
すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、遅延時間推定装置51kは、広帯域信号x(t)と、この広帯域信号が未知系を通過してきた受信信号f(t)の相互相関関数を計算し、相互相関値が最大となる遅延時間を推定遅延時間としてメモリ6に記憶する。また、遅延時間推定装置51kは、適応アルゴリズムを用いて推定遅延時間での荷重を推定して推定荷重をメモリ6に記憶する。次に、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の相互相関関数を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ701〜705である。
That is, when the switch 3 is switched to the
実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51kは、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52kは、メモリ6から推定荷重を呼び出して初期値として設定し、適応アルゴリズムを用いて荷重gkを推定し、減算器9kは、推定遅延時間と推定荷重(^)gkから誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の荷重を求める処理から最後の誤差信号を求める処理までを繰り返す。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ751〜754である。
When the interference wave is actually suppressed, the switch 3 is switched to the
実施の形態8.
この発明の実施の形態8に係る干渉波除去装置について図10を参照しながら説明する。図10は、この発明の実施の形態8に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態8に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
Embodiment 8 FIG.
An interference wave canceller according to Embodiment 8 of the present invention will be described with reference to FIG. FIG. 10 is a flowchart showing the operation of the interference wave canceller according to Embodiment 8 of the present invention. Note that the configuration of the interference wave canceller according to Embodiment 8 of the present invention is the same as that of
ここでは簡単のためにプラントノイズは十分小さい状況により説明する。上記の実施の形態7では、予め広帯域信号について処理し、1,…,K段の遅延時間推定値を求め、それぞれの遅延時間推定値における、1,…,K段の適応フィルタの荷重もそれぞれメモリ6に記憶し、干渉波抑圧時の第k段の荷重計算時の初期値に使用するものであるが、この実施の形態8では、予め広帯域信号の処理時に得られた荷重値を干渉抑圧時にも荷重値として用いるものである。
Here, for the sake of simplicity, the description will be made based on a situation where the plant noise is sufficiently small. In the seventh embodiment, a wideband signal is processed in advance to obtain 1,..., K-stage delay time estimates, and the 1,. Although stored in the
予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51k、M段の適応フィルタ52k、メモリ6が区別可能なように第1の制御信号(点線)を送る。前処理では図10中の「予め広帯域信号で行う処理」を行う。
Consider a case where a broadband signal (white noise or the like) is transmitted from the
ステップ801において、広帯域信号を送信し、広帯域信号と、この広帯域信号が未知系を通過した受信信号との相互相関関数rxfを計算する。
In
次に、ステップ802において、相互相関関数rxfはd=dkのみに相関値のピークが現れる。この唯一のピークをもつ探索パラメータdを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。
Next, in
次に、ステップ803では、ステップ802で得られた遅延時間推定値(^)dk−preで少タップ段数の適応フィルタ52kを動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとし、これをメモリ6に記憶する。
Next, in
次に、ステップ804では、遅延時間推定値(^)dk−preと荷重(^)gk−preを用いて誤差信号ek−pre(t)を次式により計算する。 Next, in step 804, the error signal e k-pre (t) is calculated by the following equation using the delay time estimated value (^) d k-pre and the load (^) g k-pre .
ただし、k=1のとき、ek−1−pre(t)は未知系の出力f(t)とする。 However, when k = 1, ek-1-pre (t) is an unknown output f (t).
そして、ステップ805において、ステップ804で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまでステップ801〜804を繰り返す。ここまでで前処理は終了する。
In
以下、ステップ851〜852の処理は、上記の実施の形態4のステップ451〜452とそれぞれ同様であるので説明は省略する。
Hereinafter, the processing in
すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、遅延時間推定装置51kは、広帯域信号x(t)と、この広帯域信号が未知系を通過してきた受信信号f(t)の相互相関関数を計算し、相互相関値が最大となる遅延時間を推定遅延時間としてメモリ6に記憶する。また、遅延時間推定装置51kは、適応アルゴリズムを用いて推定遅延時間での荷重を推定して推定荷重をメモリ6に記憶する。次に、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の相互相関関数を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ801〜805である。
That is, when the switch 3 is switched to the
実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51k、M段の適応フィルタ52k、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51kは、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52kは、メモリ6から推定荷重を呼び出して設定する。減算器9kは、推定遅延時間と推定荷重(^)gkから誤差信号を求める。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ851〜852である。
When the interference wave is actually suppressed, the switch 3 is switched to the
1 広帯域信号発生器、2 送信信号発生器、3 スイッチ、4 送信アンテナ、5 複素タップ選択適応フィルタ、6 メモリ、7 受信アンテナ、8 加算器、9、91、9k 減算器、511、51k 遅延時間推定装置、521、52k 適応フィルタ。
1-wideband signal generator, second transmission signal generator, 3 switches, 4 transmit antennas, 5 complex tap selection adaptive filter, 6 memory, 7 receiving antenna, 8 adders, 9, 9 1, 9 k
Claims (12)
送信信号を発生する送信信号発生器と、
前記広帯域信号を未知系へ送信するために前記広帯域信号発生器側に切り替えられると第1の制御信号を出力し、前記送信信号を前記未知系へ送信するために前記送信信号発生器側に切り替えられると第2の制御信号を出力するスイッチと、
推定される遅延時間を記憶するためのメモリと、
前記第1の制御信号が入力されると、遅延時間評価関数に基づき遅延時間を推定して推定遅延時間を前記メモリに記憶するとともに、前記第2の制御信号が入力されると、前記メモリから推定遅延時間を呼び出し適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求める複素タップ選択適応フィルタと、
前記送信信号が前記未知系を通過した受信信号から前記複素タップ選択適応フィルタの出力を減算することにより干渉波を抑圧する減算器と
を備えたことを特徴とする干渉波除去装置。 A broadband signal generator for generating a broadband signal;
A transmission signal generator for generating a transmission signal;
When switched to the broadband signal generator side to transmit the broadband signal to an unknown system, a first control signal is output, and switched to the transmission signal generator side to transmit the transmission signal to the unknown system A switch that outputs a second control signal when received,
A memory for storing the estimated delay time;
When the first control signal is input, a delay time is estimated based on a delay time evaluation function and the estimated delay time is stored in the memory. When the second control signal is input, the memory A complex tap selection adaptive filter that calls an estimated delay time, estimates a weight at the estimated delay time using an adaptive algorithm, and obtains an error signal from the estimated delay time and the estimated weight;
An interference wave canceling apparatus comprising: a subtractor that suppresses an interference wave by subtracting an output of the complex tap selection adaptive filter from a reception signal in which the transmission signal passes through the unknown system.
送信信号を発生する送信信号発生器と、
前記広帯域信号を未知系へ送信するために前記広帯域信号発生器側に切り替えられると第1の制御信号を出力し、前記送信信号を前記未知系へ送信するために前記送信信号発生器側に切り替えられると第2の制御信号を出力するスイッチと、
推定される遅延時間及び荷重を記憶するためのメモリと、
前記第1の制御信号が入力されると、遅延時間評価関数に基づき遅延時間を推定して推定遅延時間を前記メモリに記憶し、適応アルゴリズムを用いて前記推定遅延時間での荷重を推定して推定荷重を前記メモリに記憶するとともに、前記第2の制御信号が入力されると、前記メモリから推定遅延時間及び推定荷重を呼び出し前記推定荷重を初期値として適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び前記第2の制御信号の入力時の推定荷重より誤差信号を求める複素タップ選択適応フィルタと、
前記送信信号が前記未知系を通過した受信信号から前記複素タップ選択適応フィルタの出力を減算することにより干渉波を抑圧する減算器と
を備えたことを特徴とする干渉波除去装置。 A broadband signal generator for generating a broadband signal;
A transmission signal generator for generating a transmission signal;
When switched to the broadband signal generator side to transmit the broadband signal to an unknown system, a first control signal is output, and switched to the transmission signal generator side to transmit the transmission signal to the unknown system A switch that outputs a second control signal when received,
A memory for storing estimated delay times and weights;
When the first control signal is input, the delay time is estimated based on the delay time evaluation function, the estimated delay time is stored in the memory, and the load at the estimated delay time is estimated using an adaptive algorithm. When the estimated load is stored in the memory and the second control signal is input, the estimated delay time and the estimated load are called from the memory, and the estimated load is used as an initial value to determine the estimated delay time using the adaptive algorithm. A complex tap selection adaptive filter that obtains an error signal from an estimated delay time and an estimated weight at the time of input of the second control signal;
An interference wave canceling apparatus comprising: a subtractor that suppresses an interference wave by subtracting an output of the complex tap selection adaptive filter from a reception signal in which the transmission signal passes through the unknown system.
送信信号を発生する送信信号発生器と、
前記広帯域信号を未知系へ送信するために前記広帯域信号発生器側に切り替えられると第1の制御信号を出力し、前記送信信号を前記未知系へ送信するために前記送信信号発生器側に切り替えられると第2の制御信号を出力するスイッチと、
推定される遅延時間及び荷重を記憶するためのメモリと、
前記第1の制御信号が入力されると、遅延時間評価関数に基づき遅延時間を推定して推定遅延時間を前記メモリに記憶し、適応アルゴリズムを用いて前記推定遅延時間での荷重を推定して推定荷重を前記メモリに記憶するとともに、前記第2の制御信号が入力されると、前記メモリから推定遅延時間及び推定荷重を呼び出し、推定遅延時間及び推定荷重より誤差信号を求める複素タップ選択適応フィルタと、
前記送信信号が前記未知系を通過した受信信号から前記複素タップ選択適応フィルタの出力を減算することにより干渉波を抑圧する減算器と
を備えたことを特徴とする干渉波除去装置。 A broadband signal generator for generating a broadband signal;
A transmission signal generator for generating a transmission signal;
When switched to the broadband signal generator side to transmit the broadband signal to an unknown system, a first control signal is output, and switched to the transmission signal generator side to transmit the transmission signal to the unknown system A switch that outputs a second control signal when received,
A memory for storing estimated delay times and weights;
When the first control signal is input, the delay time is estimated based on the delay time evaluation function, the estimated delay time is stored in the memory, and the load at the estimated delay time is estimated using an adaptive algorithm. A complex tap selection adaptive filter that stores an estimated weight in the memory and, when the second control signal is input, calls the estimated delay time and estimated weight from the memory and obtains an error signal from the estimated delay time and estimated weight. When,
An interference wave canceling apparatus comprising: a subtractor that suppresses an interference wave by subtracting an output of the complex tap selection adaptive filter from a reception signal in which the transmission signal passes through the unknown system.
送信信号を発生する送信信号発生器と、
前記広帯域信号を未知系へ送信するために前記広帯域信号発生器側に切り替えられると第1の制御信号を出力し、前記送信信号を前記未知系へ送信するために前記送信信号発生器側に切り替えられると第2の制御信号を出力するスイッチと、
推定される遅延時間を記憶するためのメモリと、
前記第1の制御信号が入力されると、前記広帯域信号と前記広帯域信号が前記未知系を通過してきた受信信号との相互相関関数に基づき遅延時間を推定して推定遅延時間を前記メモリに記憶するとともに、前記第2の制御信号が入力されると、前記メモリから推定遅延時間を呼び出し適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求める複素タップ選択適応フィルタと、
前記送信信号が前記未知系を通過した受信信号から前記複素タップ選択適応フィルタの出力を減算することにより干渉波を抑圧する減算器と
を備えたことを特徴とする干渉波除去装置。 A broadband signal generator for generating a broadband signal;
A transmission signal generator for generating a transmission signal;
When switched to the broadband signal generator side to transmit the broadband signal to an unknown system, a first control signal is output, and switched to the transmission signal generator side to transmit the transmission signal to the unknown system A switch that outputs a second control signal when received,
A memory for storing the estimated delay time;
When the first control signal is input, a delay time is estimated based on a cross-correlation function between the broadband signal and the received signal through which the broadband signal has passed through the unknown system, and the estimated delay time is stored in the memory. In addition, when the second control signal is input, the estimated delay time is called from the memory, the load at the estimated delay time is estimated using an adaptive algorithm, and an error signal is obtained from the estimated delay time and the estimated load. A complex tap selection adaptive filter;
An interference wave canceling apparatus comprising: a subtractor that suppresses an interference wave by subtracting an output of the complex tap selection adaptive filter from a reception signal in which the transmission signal passes through the unknown system.
送信信号を発生する送信信号発生器と、
前記広帯域信号を未知系へ送信するために前記広帯域信号発生器側に切り替えられると第1の制御信号を出力し、前記送信信号を前記未知系へ送信するために前記送信信号発生器側に切り替えられると第2の制御信号を出力するスイッチと、
推定される遅延時間及び荷重を記憶するためのメモリと、
前記第1の制御信号が入力されると、前記広帯域信号と前記広帯域信号が前記未知系を通過してきた受信信号との相互相関関数に基づき遅延時間を推定して推定遅延時間を前記メモリに記憶し、適応アルゴリズムを用いて前記推定遅延時間での荷重を推定して推定荷重を前記メモリに記憶するとともに、前記第2の制御信号が入力されると、前記メモリから推定遅延時間及び推定荷重を呼び出し前記推定荷重を初期値として適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び前記第2の制御信号の入力時の推定荷重より誤差信号を求める複素タップ選択適応フィルタと、
前記送信信号が前記未知系を通過した受信信号から前記複素タップ選択適応フィルタの出力を減算することにより干渉波を抑圧する減算器と
を備えたことを特徴とする干渉波除去装置。 A broadband signal generator for generating a broadband signal;
A transmission signal generator for generating a transmission signal;
When switched to the broadband signal generator side to transmit the broadband signal to an unknown system, a first control signal is output, and switched to the transmission signal generator side to transmit the transmission signal to the unknown system A switch that outputs a second control signal when received,
A memory for storing estimated delay times and weights;
When the first control signal is input, a delay time is estimated based on a cross-correlation function between the broadband signal and the received signal through which the broadband signal has passed through the unknown system, and the estimated delay time is stored in the memory. Then, the load at the estimated delay time is estimated using an adaptive algorithm and the estimated load is stored in the memory. When the second control signal is input, the estimated delay time and the estimated load are read from the memory. A complex tap selection adaptive filter that estimates a load at the estimated delay time using an adaptive algorithm with the estimated weight as an initial value and obtains an error signal from the estimated delay time and the estimated weight at the time of input of the second control signal When,
An interference wave canceling apparatus comprising: a subtractor that suppresses an interference wave by subtracting an output of the complex tap selection adaptive filter from a reception signal in which the transmission signal passes through the unknown system.
送信信号を発生する送信信号発生器と、
前記広帯域信号を未知系へ送信するために前記広帯域信号発生器側に切り替えられると第1の制御信号を出力し、前記送信信号を前記未知系へ送信するために前記送信信号発生器側に切り替えられると第2の制御信号を出力するスイッチと、
推定される遅延時間及び荷重を記憶するためのメモリと、
前記第1の制御信号が入力されると、前記広帯域信号と前記広帯域信号が前記未知系を通過してきた受信信号との相互相関関数に基づき遅延時間を推定して推定遅延時間を前記メモリに記憶し、適応アルゴリズムを用いて前記推定遅延時間での荷重を推定して推定荷重を前記メモリに記憶するとともに、前記第2の制御信号が入力されると、前記メモリから推定遅延時間及び推定荷重を呼び出し、推定遅延時間及び推定荷重より誤差信号を求める複素タップ選択適応フィルタと、
前記送信信号が前記未知系を通過した受信信号から前記複素タップ選択適応フィルタの出力を減算することにより干渉波を抑圧する減算器と
を備えたことを特徴とする干渉波除去装置。 A broadband signal generator for generating a broadband signal;
A transmission signal generator for generating a transmission signal;
When switched to the broadband signal generator side to transmit the broadband signal to an unknown system, a first control signal is output, and switched to the transmission signal generator side to transmit the transmission signal to the unknown system A switch that outputs a second control signal when received,
A memory for storing estimated delay times and weights;
When the first control signal is input, a delay time is estimated based on a cross-correlation function between the broadband signal and the received signal through which the broadband signal has passed through the unknown system, and the estimated delay time is stored in the memory. Then, the load at the estimated delay time is estimated using an adaptive algorithm and the estimated load is stored in the memory. When the second control signal is input, the estimated delay time and the estimated load are read from the memory. A complex tap selection adaptive filter that obtains an error signal from the call, estimated delay time and estimated weight;
An interference wave canceling apparatus comprising: a subtractor that suppresses an interference wave by subtracting an output of the complex tap selection adaptive filter from a reception signal in which the transmission signal passes through the unknown system.
送信信号を前記未知系へ送信し、前記メモリから推定遅延時間を呼び出し適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求めるステップと、
前記送信信号が前記未知系を通過した受信信号から前記誤差信号を減算することにより干渉波を抑圧するステップと
を含むことを特徴とする干渉波除去方法。 Transmitting a broadband signal to an unknown system, estimating a delay time based on a delay time evaluation function of the broadband signal, and storing the estimated delay time in a memory;
Transmitting a transmission signal to the unknown system, calling an estimated delay time from the memory, estimating a load at the estimated delay time using an adaptive algorithm, and obtaining an error signal from the estimated delay time and the estimated weight;
A method of suppressing interference waves by subtracting the error signal from a reception signal in which the transmission signal passes through the unknown system.
送信信号を前記未知系へ送信し、前記メモリから推定遅延時間及び推定荷重を呼び出し前記推定荷重を初期値として適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び前記送信信号の送信時の推定荷重より誤差信号を求めるステップと、
前記送信信号が前記未知系を通過した受信信号から前記誤差信号を減算することにより干渉波を抑圧するステップと
を含むことを特徴とする干渉波除去方法。 A wideband signal is transmitted to an unknown system, a delay time is estimated based on a delay time evaluation function of the wideband signal, and the estimated delay time is stored in a memory, and a weight at the estimated delay time is estimated using an adaptive algorithm. Storing the estimated load in the memory;
A transmission signal is transmitted to the unknown system, an estimated delay time and an estimated weight are called from the memory, an estimated algorithm is used to estimate a load at the estimated delay time using an estimated algorithm, and the estimated delay time and the transmission Obtaining an error signal from an estimated load at the time of signal transmission;
A method of suppressing interference waves by subtracting the error signal from a reception signal in which the transmission signal passes through the unknown system.
送信信号を前記未知系へ送信し、前記メモリから推定遅延時間及び推定荷重を呼び出し、推定遅延時間及び推定荷重より誤差信号を求めるステップと、
前記送信信号が前記未知系を通過した受信信号から前記誤差信号を減算することにより干渉波を抑圧するステップと
を含むことを特徴とする干渉波除去方法。 A wideband signal is transmitted to an unknown system, a delay time is estimated based on a delay time evaluation function of the wideband signal, and the estimated delay time is stored in a memory, and a weight at the estimated delay time is estimated using an adaptive algorithm. Storing the estimated load in the memory;
Transmitting a transmission signal to the unknown system, calling an estimated delay time and estimated load from the memory, and obtaining an error signal from the estimated delay time and estimated load;
A method of suppressing interference waves by subtracting the error signal from a reception signal in which the transmission signal passes through the unknown system.
送信信号を前記未知系へ送信し、前記メモリから推定遅延時間を呼び出し適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求めるステップと、
前記送信信号が前記未知系を通過した受信信号から前記誤差信号を減算することにより干渉波を抑圧するステップと
を含むことを特徴とする干渉波除去方法。 Transmitting a broadband signal to an unknown system, estimating a delay time based on a cross-correlation function between the broadband signal and a received signal through which the broadband signal has passed through the unknown system, and storing the estimated delay time in a memory;
Transmitting a transmission signal to the unknown system, calling an estimated delay time from the memory, estimating a load at the estimated delay time using an adaptive algorithm, and obtaining an error signal from the estimated delay time and the estimated weight;
A method of suppressing interference waves by subtracting the error signal from a reception signal in which the transmission signal passes through the unknown system.
送信信号を前記未知系へ送信し、前記メモリから推定遅延時間及び推定荷重を呼び出し前記推定荷重を初期値として適応アルゴリズムを用いて前記推定遅延時間での荷重を推定し、推定遅延時間及び前記送信信号の送信時の推定荷重より誤差信号を求めるステップと、
前記送信信号が前記未知系を通過した受信信号から前記誤差信号を減算することにより干渉波を抑圧するステップと
を含むことを特徴とする干渉波除去方法。 A wideband signal is transmitted to an unknown system, a delay time is estimated based on a cross-correlation function between the wideband signal and a received signal that has passed through the unknown system, and the estimated delay time is stored in a memory and adaptively Estimating a load at the estimated delay time using an algorithm and storing the estimated load in the memory;
A transmission signal is transmitted to the unknown system, an estimated delay time and an estimated weight are called from the memory, an estimated algorithm is used to estimate a load at the estimated delay time using an estimated algorithm, and the estimated delay time and the transmission Obtaining an error signal from an estimated load at the time of signal transmission;
A method of suppressing interference waves by subtracting the error signal from a reception signal in which the transmission signal passes through the unknown system.
送信信号を前記未知系へ送信し、前記メモリから推定遅延時間及び推定荷重を呼び出し、推定遅延時間及び推定荷重より誤差信号を求めるステップと、
前記送信信号が前記未知系を通過した受信信号から前記誤差信号を減算することにより干渉波を抑圧するステップと
を含むことを特徴とする干渉波除去方法。 A wideband signal is transmitted to an unknown system, a delay time is estimated based on a cross-correlation function between the wideband signal and a received signal that has passed through the unknown system, and the estimated delay time is stored in a memory and adaptively Estimating a load at the estimated delay time using an algorithm and storing the estimated load in the memory;
Transmitting a transmission signal to the unknown system, calling an estimated delay time and estimated load from the memory, and obtaining an error signal from the estimated delay time and estimated load;
A method of suppressing interference waves by subtracting the error signal from a reception signal in which the transmission signal passes through the unknown system.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005354522A JP4545683B2 (en) | 2005-12-08 | 2005-12-08 | Interference wave removing apparatus and interference wave removing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005354522A JP4545683B2 (en) | 2005-12-08 | 2005-12-08 | Interference wave removing apparatus and interference wave removing method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2007158996A JP2007158996A (en) | 2007-06-21 |
JP4545683B2 true JP4545683B2 (en) | 2010-09-15 |
Family
ID=38242733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2005354522A Expired - Fee Related JP4545683B2 (en) | 2005-12-08 | 2005-12-08 | Interference wave removing apparatus and interference wave removing method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4545683B2 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102281044B (en) * | 2010-06-12 | 2016-07-20 | 澜起科技(上海)有限公司 | Blind adaptive filtering device and the application thereof of arrowband interference can be eliminated |
JP2012182694A (en) * | 2011-03-02 | 2012-09-20 | Panasonic Corp | Waveform equalization device |
US10567013B2 (en) | 2014-10-17 | 2020-02-18 | Tohoku University | Rotary system mounted piezoelectric generator, self-powered wireless communication terminal and wireless communication system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001308751A (en) * | 2000-02-15 | 2001-11-02 | Mitsubishi Electric Corp | Echo canceller and echo canceling method |
JP2003087217A (en) * | 2001-09-17 | 2003-03-20 | Sumitomo Electric Ind Ltd | Sneak canceler |
JP2003273831A (en) * | 2002-03-19 | 2003-09-26 | Sumitomo Electric Ind Ltd | Sneak canceller |
JP2004064681A (en) * | 2002-07-31 | 2004-02-26 | Japan Science & Technology Corp | Adaptive filter |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0746783B2 (en) * | 1985-09-18 | 1995-05-17 | 日本電気株式会社 | Interference signal elimination method |
JPH02312326A (en) * | 1989-05-26 | 1990-12-27 | Mitsubishi Electric Corp | Radio transmission and reception equipment |
JP2957289B2 (en) * | 1991-02-04 | 1999-10-04 | 防衛庁技術研究本部長 | Transceiver |
JPH0556009A (en) * | 1991-08-22 | 1993-03-05 | Mitsubishi Electric Corp | Transmitter-receiver |
JPH05218903A (en) * | 1992-01-07 | 1993-08-27 | Nec Corp | Echo canceller |
JP3584603B2 (en) * | 1996-01-05 | 2004-11-04 | 三菱電機株式会社 | Transceiver |
-
2005
- 2005-12-08 JP JP2005354522A patent/JP4545683B2/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001308751A (en) * | 2000-02-15 | 2001-11-02 | Mitsubishi Electric Corp | Echo canceller and echo canceling method |
JP2003087217A (en) * | 2001-09-17 | 2003-03-20 | Sumitomo Electric Ind Ltd | Sneak canceler |
JP2003273831A (en) * | 2002-03-19 | 2003-09-26 | Sumitomo Electric Ind Ltd | Sneak canceller |
JP2004064681A (en) * | 2002-07-31 | 2004-02-26 | Japan Science & Technology Corp | Adaptive filter |
Also Published As
Publication number | Publication date |
---|---|
JP2007158996A (en) | 2007-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4697465B2 (en) | Signal processing method, signal processing apparatus, and signal processing program | |
JP6363324B2 (en) | Signal processing apparatus, signal processing method, and signal processing program | |
JP3177562B2 (en) | Low delay subband adaptive filter device | |
US8243941B2 (en) | Active noise reduction system and active noise reduction method | |
US9830900B2 (en) | Adaptive equalizer, acoustic echo canceller device, and active noise control device | |
US9286881B2 (en) | Echo canceler and echo detector | |
JP4632047B2 (en) | Signal processing method and apparatus | |
JP2503747B2 (en) | FIR type eco-canceller | |
EP0288577B1 (en) | Echo canceller with short processing delay and decreased multiplication number and method for controlling an echo signal | |
KR100299290B1 (en) | Echo cancellation methods and echo cancellers that implement such processes | |
US9536539B2 (en) | Nonlinear acoustic echo signal suppression system and method using volterra filter | |
JP2850814B2 (en) | Method and apparatus for identifying unknown system using adaptive filter | |
JP4545683B2 (en) | Interference wave removing apparatus and interference wave removing method | |
JP4591685B2 (en) | Double talk state determination method, echo cancellation method, double talk state determination device, echo cancellation device, and program | |
JP2002135171A (en) | System identification method | |
US5737409A (en) | Echo removing apparatus | |
US20100260352A1 (en) | System identification device and system identification method | |
KR100738713B1 (en) | Echo cancellers for sparse channels | |
TW200412728A (en) | Noise matching for echo cancellers | |
EP3667662B1 (en) | Acoustic echo cancellation device, acoustic echo cancellation method and acoustic echo cancellation program | |
JP5167871B2 (en) | Propagation delay time estimator, program and method, and echo canceller | |
JP2017098861A (en) | Echo canceller and echo cancellation method | |
US20060288067A1 (en) | Reduced complexity recursive least square lattice structure adaptive filter by means of approximating the forward error prediction squares using the backward error prediction squares | |
JP5193618B2 (en) | Interference wave canceller | |
JP4367243B2 (en) | Adaptive phasing device, its program and adaptive phasing system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080606 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100621 |
|
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: 20100629 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100630 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130709 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |