JP4545683B2 - Interference wave removing apparatus and interference wave removing method - Google Patents

Interference wave removing apparatus and interference wave removing method Download PDF

Info

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
Application number
JP2005354522A
Other languages
Japanese (ja)
Other versions
JP2007158996A (en
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 JP2005354522A priority Critical patent/JP4545683B2/en
Publication of JP2007158996A publication Critical patent/JP2007158996A/en
Application granted granted Critical
Publication of JP4545683B2 publication Critical patent/JP4545683B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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 delay circuit 90 are multiplied by weights g and g * having appropriate complex values by multipliers 91 and 92, respectively. Subtractors 93 and 94 respectively subtract from e k−1 (t). These are obtained as mean square values by means of mean square value calculation devices 95 and 96 (steps 901 and 902 in FIG. 12). The load g * is a complex conjugate of the load g. The operation so far is performed so that the maximum search value d of the delay time is satisfied so as to satisfy the time axis of the impulse response of the unknown system, and the delay circuit 90 changes the d delay times to obtain two mean square values. These two mean square values are added by the adder 97 for each delay time and input to the mean square value evaluation circuit 98 (step 903 in FIG. 12).

加算した二乗平均値とその一次近似直線間の偏差が最大になる遅延時間を求め、それを遅延時間dとして出力する。以上の操作を行うことでインパルス応答値が局在する遅延時間dを推定する。 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)が遅延時間dを決定する原理となる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.

Figure 0004545683
Figure 0004545683

ただし、k=1のときは、誤差信号ek−1(t)は未知系の出力f(t)とする。ここで、dは遅延回路90において探索する遅延時間、e(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 delay circuit 90, e k (t) is an error signal, t is a time factor of the digital signal, g is a weight of the adaptive filter, and g * is a complex conjugate with the weight g of the adaptive filter. N ′ is the number of samples using the adaptive filter weight g, NN ′ is the number of samples using the adaptive filter weight g * , and x is the transmission signal.

次に、求められた誤差信号の二乗平均値の一次近似直線を求める。この直線と誤差信号の二乗平均値の偏差が最大値となる探索遅延時間dを遅延時間(^)d(なお、(^)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).

Figure 0004545683
Figure 0004545683

ここで、式(2)の右辺第1項は、次式の通りである。   Here, the first term on the right side of Equation (2) is as follows.

Figure 0004545683
Figure 0004545683

ここで、cは未知系のインパルス応答値、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=dになるとき、rxx(d−d)は最大値をもち、一次近似直線との偏差が最大になる。このときの探索遅延時間dを(^)dとすることでインパルス応答値が局在する遅延時間dを推定することができる(図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 (step 904 in FIG. 12).

特開2001−308751号公報JP 2001-308751 A

しかし、これまで説明してきた複素タップ選択適応フィルタを用いた干渉波除去装置では、送信信号が狭帯域で周期性をもつ場合、遅延時間評価関数も周期性をもち真値以外の遅延時間を推定値とすることがみられた。この誤った遅延時間推定値のため、未知伝播系のインパルス応答値の推定も不可能になり、抑圧性能が劣化することが見られた。   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.

送信信号が無変調の正弦波のように周期性をもつ信号の場合に、遅延時間推定の評価関数Φ’(d)の振幅が振動する問題点を説明する。 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=dとなるとき、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、未知系のインパルス応答もcのみ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.

Figure 0004545683
Figure 0004545683

ここで、送信信号が複素の正弦波e−j2πftとすると、式(4)は、以下のようになる。 Here, if the transmission signal is a complex sine wave e −j2πft , Equation (4) is as follows.

Figure 0004545683
Figure 0004545683

ここで、〈〉内は、期待値演算を表す。式(5)において、dの探索と無関係にe−j2πf(d−d1)により周期性が生じるため、遅延時間探索中にΦ’(d)の振幅が振動することがわかる。なお、e−j2πf(d−d1)における上付きの(d−d1)は、(d−d)である。これらのことから、送信信号が周期関数の場合、評価関数が周期性をもつため、遅延時間の推定を誤ってしまう恐れがあることが分かる。このときの遅延時間推定方法を、模式的に表したものを図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に係る干渉波除去装置の複素タップ選択適応フィルタの構成を示す図である。なお、各図中、同一符号は同一又は相当部分を示す。
Embodiment 1 FIG.
An interference wave canceling apparatus according to Embodiment 1 of the present invention will be described with reference to FIGS. 1 is a diagram showing a configuration of an interference wave canceling apparatus according to Embodiment 1 of the present invention. FIG. 2 is a diagram showing the configuration of the complex tap selection adaptive filter of the interference wave canceller according to Embodiment 1 of the present invention. In addition, in each figure, the same code | symbol shows the same or equivalent part.

図1において、この発明の実施の形態1に係る干渉波除去装置は、広帯域信号を発生する広帯域信号発生器1と、送信信号を発生する送信信号発生器2と、広帯域信号発生器1と送信信号発生器2を切り替えるスイッチ3と、広帯域信号や送信信号x(t)を未知系(unknown system)へ送信する送信アンテナ4と、未知系のインパルス応答が局在する場合にその遅延時間を推定し、その遅延時間に当たるタップのみ適応動作させ、複素数化を図った複素タップ選択適応フィルタ5と、推定した遅延時間などを記憶するメモリ6と、
所望波を受信する受信アンテナ7と、受信信号f(t)とプラントノイズn(t)を加算する加算器8と、加算器8の出力から適応フィルタ5の出力を減算して所望波e(t)を出力する減算器9(9,…,9,…,9)とが設けられている。なお、スイッチ3と送信アンテナ4の間には、変調器などの送信手段の図示が省略され、受信アンテナ7と加算器8の間には、復調器などの受信手段の図示が省略されている。また、送信アンテナ4から受信アンテナ7への矢印は、漏れこみ干渉波すなわち未知系を経由してきた送信信号を表す。
In FIG. 1, an interference wave canceling apparatus according to Embodiment 1 of the present invention includes a wideband signal generator 1 that generates a wideband signal, a transmission signal generator 2 that generates a transmission signal, a wideband signal generator 1, and a transmission. A switch 3 for switching the signal generator 2, a transmission antenna 4 for transmitting a broadband signal or a transmission signal x (t) to an unknown system, and an estimation of the delay time when the impulse response of the unknown system is localized A complex tap selection adaptive filter 5 adapted to adaptively operate only taps corresponding to the delay time, and a complex number, and a memory 6 for storing the estimated delay time, etc.
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 adaptive filter 5 from the output of the adder 8 Subtractors 9 (9 1 ,..., 9 k ,..., 9 K ) that output (t) are provided. Note that transmission means such as a modulator is not shown between the switch 3 and the transmission antenna 4, and reception means such as a demodulator is not shown between the reception antenna 7 and the adder 8. . An arrow from the transmission antenna 4 to the reception antenna 7 represents a leakage interference wave, that is, a transmission signal that has passed through an unknown system.

また、図2において、複素タップ選択適応フィルタ5は、遅延時間推定装置51,…,51,…,51と、M段の適応フィルタ52,…,52,…,52とが設けられている。なお、図1及び図2において、点線は制御信号を表す。 Further, in FIG. 2, the complex tap selection adaptive filter 5, the delay time estimating unit 51 1, ..., 51 k, ..., and 51 K, the adaptive filter 52 1 M stage, ..., 52 k, ..., and 52 K Is provided. 1 and 2, the dotted line represents a control signal.

つぎに、この実施の形態1に係る干渉波除去装置の動作について図面を参照しながら説明する。図3は、この発明の実施の形態1に係る干渉波除去装置の動作を示すフローチャートである。   Next, the operation of the interference wave canceller according to Embodiment 1 will be described with reference to the drawings. FIG. 3 is a flowchart showing the operation of the interference wave canceller according to Embodiment 1 of the present invention.

ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態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は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51、M段の適応フィルタ52、メモリ6が区別可能なように第1の制御信号(点線)を送る。複素タップ選択適応フィルタ5の動作を図3に示す。ここでは、図3中の「予め広帯域信号で行う処理」を行う。 Consider a case where a broadband signal (white noise or the like) is transmitted from the broadband signal generator 1 in advance. At this time, the switch 3 switches to the wideband signal generator 1 side. Further, the switch 3 uses the first control signal (1) so that the delay time estimation device 51 1 , the M-stage adaptive filter 52 1 , and the memory 6 can distinguish that the delay time is estimated using the wideband signal. Send a dotted line). The operation of the complex tap selection adaptive filter 5 is shown in FIG. Here, the “processing performed in advance with a broadband signal” in FIG. 3 is performed.

ステップ101、102、103において、広帯域信号を送信した場合、式(2)中の送信信号の自己相関関数rxxはd=dのみに相関値のピークが現れる。この場合、遅延時間評価関数Φ’(d)は唯一のピークをもつ。 When a wideband signal is transmitted in steps 101, 102, and 103, the autocorrelation function rxx of the transmission signal in equation (2) has a correlation value peak only at d = d i . In this case, the delay time evaluation function Φ ′ (d) has only one peak.

次に、ステップ104において、この唯一のピークをもつ探索パラメータdを遅延時間推定値(^)d1−preとする。(なお、下付きの1−preは1_preを表す。)この遅延時間推定値(^)d1−preをメモリ6に記憶する。ここまでで前処理は終了する。 Next, in step 104, and the only time delay estimate the search parameters d 1 having a peak (^) d 1-pre. (The subscript 1-pre represents 1_pre.) This delay time estimated value (^) d 1-pre is stored in the memory 6. Thus, the preprocessing ends.

実際に干渉波を抑圧処理する場合には、図3中の「干渉波抑圧時の処理」を行う。このとき、スイッチ3は送信信号発生器2側に切り替えて、遅延時間推定装置51、M段の適応フィルタ52、メモリ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 transmission signal generator 2 side and sends a second control signal (dotted line) to the delay time estimation device 51 1 , the M-stage adaptive filter 52 1 , and the memory 6.

ステップ151では、複素タップ選択適応フィルタ5の遅延時間推定の処理を行わず、予め広帯域信号を送信して求めた推定遅延時間(^)d1−preをメモリ6から呼び出し、直接前処理で求めた値を遅延時間推定値(^)d=(^)d1−preとする。この遅延時間(^)d=(^)d1−preにおいて、適応フィルタ52を動作させ、適応アルゴリズムを用いて、その荷重gを計算し、計算後のgを荷重推定値(^)gとする。適応アルゴリズムの種類については、正規化LMSアルゴリズムなど特に問わない。 In step 151, the delay time estimation processing of the complex tap selection adaptive filter 5 is not performed, but the estimated delay time (^) d 1-pre obtained by transmitting a broadband signal in advance is called from the memory 6 and directly obtained by preprocessing. This value is assumed to be a delay time estimated value (^) d 1 = (^) d 1-pre . In this delay time (^) d 1 = (^ ) d 1-pre, adaptive filter 52 1 is operated by using an adaptive algorithm to calculate the load g 1, load estimation value g 1 after calculation ( ^) and g 1. The type of adaptive algorithm is not particularly limited, such as a normalized LMS algorithm.

次に、ステップ152では、ステップ151で用いた遅延時間(^)d=(^)d1−preと求められた荷重gを用いて、誤差信号e(t)を求める。この行程は、次式で表すことができる。 Next, in step 152, the error signal e 1 (t) is obtained using the delay time (^) d 1 = (^) d 1-pre used in step 151 and the obtained load g 1 . This process can be expressed by the following equation.

Figure 0004545683
Figure 0004545683

ここで、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は、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、適応フィルタ52は、メモリ6から荷重g、gを呼び出して設定し、遅延時間推定装置51は、荷重g、gを用いてそれぞれ誤差信号の二乗平均値を求め、2つの二乗平均値から遅延時間評価関数を求め、この遅延時間評価関数から遅延時間を推定して推定遅延時間をメモリ6に記憶する。これらの前処理がステップ101〜104である。 That is, when the switch 3 is switched to the wideband signal generator 1 side, the switch 3 sends a first control signal to the delay time estimation device 51 1 , the M-stage adaptive filter 52 1 , and the memory 6. Based on this first control signal, the adaptive filter 52 1, a load g from the memory 6, and set by calling the g *, the delay time estimating unit 51 1, the load g, g * of each error signal by using A mean square value is obtained, a delay time evaluation function is obtained from the two mean square values, a delay time is estimated from the delay time evaluation function, and the estimated delay time is stored in the memory 6. These pre-processing are steps 101-104.

実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52は、適応アルゴリズムを用いて荷重gを求め、減算器9は、推定遅延時間と荷重gから誤差信号を求め、これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ151〜154である。 When the interference wave is actually suppressed, the switch 3 is switched to the transmission signal generator 2 side, and the second control signal is sent to the delay time estimation device 51 1 , the M-stage adaptive filter 52 1 , and the memory 6. Send. Based on this second control signal, the delay time estimating unit 51 1 functions as a delay circuit by setting calls the estimated delay time from the memory 6, the adaptive filter 52 1, a load g using an adaptive algorithm 1 determined, subtractor 9 1 obtains an error signal from the estimated delay time and the load g 1, which interference wave is suppressed by. These interference wave suppression processes are steps 151 to 154.

実施の形態2.
この発明の実施の形態2に係る干渉波除去装置について図4を参照しながら説明する。図4は、この発明の実施の形態2に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態2に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
Embodiment 2. FIG.
An interference wave removing apparatus according to Embodiment 2 of the present invention will be described with reference to FIG. FIG. 4 is a flowchart showing the operation of the interference wave canceller according to Embodiment 2 of the present invention. The configuration of the interference wave canceller according to Embodiment 2 of the present invention is the same as that of Embodiment 1 above.

ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態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は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51、タップ数M段の適応フィルタ52、メモリ6が区別可能なように第1の制御信号(点線)を送る。複素タップ選択適応フィルタ5の動作を図4に示す。ここでは、図4中の「予め広帯域信号で行う処理」を行う。 Consider a case where a broadband signal (white noise or the like) is transmitted from the broadband signal generator 1 in advance. At this time, the switch 3 switches to the wideband signal generator 1 side. Further, the switch 3 performs the first control so that the delay time estimation device 51 k , the M-tap adaptive filter 52 k , and the memory 6 can distinguish that the delay time is estimated using the wideband signal. Send a signal (dotted line). The operation of the complex tap selection adaptive filter 5 is shown in FIG. Here, the “processing to be performed in advance with a broadband signal” in FIG. 4 is performed.

ステップ201、202、203において、広帯域信号を送信した場合、式(2)中の送信信号の自己相関関数rxxはd=dのみに相関値のピークが現れる。この場合、遅延時間評価関数Φ’(d)は唯一のピークをもつ。 In step 201, 202 and 203, when transmitting broadband signals, the autocorrelation function r xx of the transmitted signal in the formula (2) is the peak of the correlation value appears only in d = d k. In this case, the delay time evaluation function Φ ′ (d) has only one peak.

次に、ステップ204において、この唯一のピークをもつ探索パラメータdを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。 Next, in step 204, the search parameter d i having the only peak is set as a delay time estimated value (^) d k-pre . This estimated delay time value (^) d k-pre is stored in the memory 6.

次に、ステップ205では、ステップ204で得られた遅延時間推定値で少タップ段数の適応フィルタ52を動作させ、適応アルゴリズムを用いて、その荷重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 step 204, using the adaptive algorithm, and calculates the load g k-pre, after calculation g load estimation value of the k-pre (^) and g k-pre.

次に、ステップ206では、遅延時間推定値(^)dk−preと荷重(^)gk−preより誤差信号ek−pre(t)を求める。 Next, in step 206, an error signal e k-pre (t) is obtained from the delay time estimated value (^) d k-pre and the load (^) g k-pre .

そして、ステップ207では、ステップ206で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、この繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまで、ステップ201〜206を繰り返す。ここまでで前処理は終了する。   In step 207, if the minimum value of the error signal obtained in step 206 is smaller than a preset threshold value, this repetitive processing is stopped. In other cases, k = k + 1 is set, and steps 201 to 206 are repeated until k = K. Thus, the preprocessing ends.

実際に干渉波を抑圧処理する場合には、図4中の「干渉波抑圧時の処理」を行う。このとき、スイッチ3は、送信信号発生器2側に切り替えて、遅延時間推定装置51、M段の適応フィルタ52、メモリ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 transmission signal generator 2 side and sends the second control signal (dotted line) to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6.

ステップ251では、複素タップ選択適応フィルタ5の遅延時間推定の処理を行わず、予め広帯域信号を送信して求めた推定遅延時間(^)dk−preをメモリ6から呼び出し、前処理で求めた(^)dk−preを第k段の遅延時間推定値(^)d=(^)dk−preとする。この遅延時間(^)d=(^)dk−preにおいて、タップ数M段の適応フィルタ52を動作させ、計算後の荷重値gを荷重推定値(^)gとし、この荷重値を未知系のインパルス応答推定値とする。適応アルゴリズムの種類については問わない。 In step 251, the delay time estimation process of the complex tap selection adaptive filter 5 is not performed, and the estimated delay time (^) d k-pre obtained by transmitting a broadband signal in advance is called from the memory 6 and obtained by preprocessing. Let (^) d k-pre be the k-th stage delay time estimate (^) d k = (^) d k-pre . In this delay time (^) d k = (^ ) d k-pre, operates the adaptive filter 52 k taps M stage, and the load value g k after calculation load estimation value (^) g k, this The load value is assumed to be an impulse response estimation value of an unknown system. It doesn't matter about the type of adaptive algorithm.

次に、ステップ252では、ステップ251で用いた遅延時間(^)d=(^)dk−preと求められた荷重(^)gを用いて、誤差信号e(t)を求める。この行程は次式で表すことができる。 Next, in step 252, the error signal e k (t) is obtained using the delay time (^) d k = (^) d k-pre used in step 251 and the obtained load (^) g k. . This process can be expressed by the following equation.

Figure 0004545683
Figure 0004545683

ただし、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 step 253, if the minimum value of the error signal obtained in step 252 is smaller than a preset threshold value, the iterative process is stopped. In other cases, k = k + 1 is set, and the processes of steps 251 and 252 are repeated until k = K.

これらの処理が終わった時点で全ての未知系のインパルス応答が局在する遅延時間と、その大きさの推定が行なえる。   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は、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、適応フィルタ52は、メモリ6から荷重g、gを呼び出して設定し、遅延時間推定装置51は、荷重g、gを用いてそれぞれ誤差信号の二乗平均値を求め、2つの二乗平均値から遅延時間評価関数を求め、この遅延時間評価関数から遅延時間を推定して推定遅延時間をメモリ6に記憶する。また、遅延時間推定装置51は、適応アルゴリズムを用いて推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の二乗平均値を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ201〜207である。 That is, when the switch 3 is switched to the wideband signal generator 1 side, the switch 3 sends a first control signal to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Based on the first control signal, the adaptive filter 52 k calls and sets the loads g and g * from the memory 6, and the delay time estimation device 51 k uses the loads g and g * , respectively. A mean square value is obtained, a delay time evaluation function is obtained from the two mean square values, a delay time is estimated from the delay time evaluation function, and the estimated delay time is stored in the memory 6. Further, the delay time estimation device 51 k estimates a load at the estimated delay time using an adaptive algorithm, and obtains an error signal from the estimated delay time and the estimated load. Until the minimum value of the error signal reaches the threshold value, the process from obtaining the first mean square value to the process obtaining the last error signal is repeated. These pre-processing are steps 201-207.

実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52は、適応アルゴリズムを用いて荷重(^)gを求めて、減算器9は、推定遅延時間と推定荷重(^)gから誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の荷重を求める処理から最後の誤差信号を求める処理までを繰り返す。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ251〜254である。 When the interference wave is actually suppressed, the switch 3 is switched to the transmission signal generator 2 side, and the second control signal is sent to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Send. Based on the second control signal, the delay time estimation device 51 k functions as a delay circuit by calling and setting the estimated delay time from the memory 6, and the adaptive filter 52 k uses a weight ( ^) G k is obtained, and the subtractor 9 k obtains an error signal from the estimated delay time and the estimated weight (^) g k . The process from obtaining the first load to the process obtaining the last error signal is repeated until the minimum value of the error signal reaches the threshold value. Thereby, the interference wave is suppressed. These interference wave suppression processes are steps 251 to 254.

実施の形態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 Embodiment 1 above.

ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態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は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51、タップ数M段の適応フィルタ52、メモリ6が区別可能なように第1の制御信号(点線)を送る。複素タップ選択適応フィルタ5の動作を図5に示す。ここでは、図5中の「予め広帯域信号で行う処理」を行う。 Consider a case where a broadband signal (white noise or the like) is transmitted from the broadband signal generator 1 in advance. At this time, the switch 3 switches to the wideband signal generator 1 side. Further, the switch 3 performs the first control so that the delay time estimation device 51 k , the M-tap adaptive filter 52 k , and the memory 6 can distinguish that the delay time is estimated using the wideband signal. Send a signal (dotted line). The operation of the complex tap selection adaptive filter 5 is shown in FIG. Here, the “processing performed in advance with a broadband signal” in FIG. 5 is performed.

ステップ301、302、303において、広帯域信号を送信した場合、式(2)中の送信信号の自己相関関数rxxはd=dのみに相関値のピークが現れる。この場合、遅延時間評価関数Φ’(d)は唯一のピークをもつ。 In Steps 301, 302, and 303, when a broadband signal is transmitted, a correlation value peak appears only in d k = d i in the autocorrelation function r xx of the transmission signal in Equation (2). In this case, the delay time evaluation function Φ ′ (d) has only one peak.

次に、ステップ304において、この唯一のピークをもつ探索パラメータdを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。 Next, in step 304, the search parameter d i having the only peak is set as a delay time estimated value (^) d k-pre . This estimated delay time value (^) d k-pre is stored in the memory 6.

次に、ステップ305では、ステップ304で得られた遅延時間推定値で少タップ段数の適応フィルタ52を動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとし、これをメモリ6に記憶する。 Next, in step 305, operates the adaptive filter 52 k of small number of tap stages in the delay time estimation value obtained in step 304, using the adaptive algorithm, and calculates the load g k-pre, after calculation g k-pre is assumed to be a load estimated value (^) g k-pre and is stored in the memory 6.

次に、ステップ306では、遅延時間(^)dk−preと荷重(^)gk−preより誤差信号ek−pre(t)を求める。 Next, in step 306, an error signal e k-pre (t) is obtained from the delay time (^) d k-pre and the load (^) g k-pre .

そして、ステップ307において、ステップ306で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまでステップ301〜306を繰り返す。ここまでで前処理は終了する。   In step 307, if the minimum value of the error signal obtained in step 306 is smaller than a preset threshold value, the iterative process is stopped. In other cases, k = k + 1 is set, and steps 301 to 306 are repeated until k = K. Thus, the preprocessing ends.

実際に干渉波を抑圧処理する場合には、図5中「干渉波抑圧時の処理」を行う。このとき、スイッチ3は送信信号発生器2側に切り替えて、遅延時間推定装置51、M段の適応フィルタ52、メモリ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 transmission signal generator 2 side and sends the second control signal (dotted line) to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6.

ステップ351では、複素タップ選択適応フィルタ5の遅延時間推定の処理を行わず、予め広帯域信号を送信して求めた推定遅延時間(^)dk−preと荷重(^)gk−preをメモリ6から呼び出し、前処理で求めた(^)dk−preを第k段の遅延時間推定値(^)d=(^)dk−preとする。この遅延時間(^)d=(^)dk−preにおいて、タップ数M段の適応フィルタ52を動作させ、その荷重gを計算し、計算後の荷重値を(^)gとし未知系のインパルス応答推定値とする。ここでは、荷重更新の初期値として推定荷重(^)gk−preを用いる。これにより、適応アルゴリズムの収束時間の短縮が可能になる。適応アルゴリズムの種類については問わない。 In step 351, the delay time estimation process of the complex tap selection adaptive filter 5 is not performed, and the estimated delay time (^) d k-pre and the weight (^) g k-pre obtained by transmitting a broadband signal in advance are stored in the memory. 6 and (^) d k-pre obtained in the pre-processing is set as the k-th stage delay time estimated value (^) d k = (^) d k-pre . In this delay time (^) dk = (^) dk-pre , the adaptive filter 52k with M taps is operated, the load gk is calculated, and the calculated load value is (^) gk. And an unknown impulse response estimate. Here, the estimated load (^) g k-pre is used as the initial value of the load update. Thereby, the convergence time of the adaptive algorithm can be shortened. It doesn't matter about the type of adaptive algorithm.

次に、ステップ352では、ステップ351で用いた遅延時間(^)d=(^)dk−preと求められた荷重(^)gを用いて、誤差信号e(t)を求める。この行程は次式で表すことができる。 Next, in step 352, the error signal e k (t) is obtained using the delay time (^) d k = (^) d k-pre used in step 351 and the obtained load (^) g k. . This process can be expressed by the following equation.

Figure 0004545683
Figure 0004545683

ただし、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 step 353, if the minimum value of the error signal obtained in step 352 is smaller than a preset threshold value, the iterative process is stopped. In other cases, k = k + 1 is set, and steps 351 and 352 are repeated until k = K.

これらの処理が終わった時点で全ての未知系のインパルス応答が局在する遅延時間とその大きさの推定が行なえる。   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は、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、適応フィルタ52は、メモリ6から荷重g、gを呼び出して設定し、遅延時間推定装置51は、荷重g、gを用いてそれぞれ誤差信号の二乗平均値を求め、2つの二乗平均値から遅延時間評価関数を求め、この遅延時間評価関数から遅延時間を推定して推定遅延時間をメモリ6に記憶する。また、遅延時間推定装置51は、適応アルゴリズムを用いて推定遅延時間での荷重を推定して推定荷重をメモリ6に記憶する。次に、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の二乗平均値を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ301〜307である。 That is, when the switch 3 is switched to the wideband signal generator 1 side, the switch 3 sends a first control signal to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Based on the first control signal, the adaptive filter 52 k calls and sets the loads g and g * from the memory 6, and the delay time estimation device 51 k uses the loads g and g * , respectively. A mean square value is obtained, a delay time evaluation function is obtained from the two mean square values, a delay time is estimated from the delay time evaluation function, and the estimated delay time is stored in the memory 6. Further, the delay time estimation device 51 k estimates the load at the estimated delay time using an adaptive algorithm and stores the estimated load in the memory 6. Next, an error signal is obtained from the estimated delay time and estimated load. Until the minimum value of the error signal reaches the threshold value, the process from obtaining the first mean square value to the process obtaining the last error signal is repeated. These preprocessing are steps 301 to 307.

実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52は、メモリ6から推定荷重を呼び出して初期値として設定し、適応アルゴリズムを用いて荷重gを推定し、減算器9は、推定遅延時間と推定荷重(^)gから誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の荷重を求める処理から最後の誤差信号を求める処理までを繰り返す。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ351〜354である。 When the interference wave is actually suppressed, the switch 3 is switched to the transmission signal generator 2 side, and the second control signal is sent to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Send. Based on the second control signal, the delay time estimation device 51 k functions as a delay circuit by calling and setting the estimated delay time from the memory 6, and the adaptive filter 52 k calls the estimated weight from the memory 6. As an initial value, the load g k is estimated using an adaptive algorithm, and the subtractor 9 k obtains an error signal from the estimated delay time and the estimated load (^) g k . The process from obtaining the first load to the process obtaining the last error signal is repeated until the minimum value of the error signal reaches the threshold value. Thereby, the interference wave is suppressed. These interference wave suppression processes are steps 351 to 354.

実施の形態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 Embodiment 1 above.

ここでは簡単のためにプラントノイズは十分小さい状況により説明する。上記の実施の形態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 memory 6 and used as an initial value at the time of calculating the k-th stage load when suppressing the interference wave, in the fourth embodiment, the load value obtained in advance during the processing of the broadband signal is used to suppress the interference. Sometimes used as a load value.

予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51、M段の適応フィルタ52、メモリ6が区別可能なように第1の制御信号(点線)を送る。ここでは、図6中の「予め広帯域信号で行う処理」を行う。 Consider a case where a broadband signal (white noise or the like) is transmitted from the broadband signal generator 1 in advance. At this time, the switch 3 switches to the wideband signal generator 1 side. In addition, the switch 3 uses the first control signal ( k ) so that the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6 can distinguish that the delay time is estimated using the wideband signal. Send a dotted line). Here, the “processing performed in advance with a broadband signal” in FIG. 6 is performed.

ステップ401、402、403において、広帯域信号を送信した場合、式(2)中の送信信号の自己相関関数rxxはd=dのみに相関値のピークが現れる。この場合、遅延時間評価関数Φ’(d)は唯一のピークをもつ。 In step 401, when transmitting broadband signals, the autocorrelation function r xx of the transmitted signal in the formula (2) is the peak of the correlation value appears only in d = d k. In this case, the delay time evaluation function Φ ′ (d) has only one peak.

次に、ステップ404において、この唯一のピークをもつ探索パラメータdを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。 Next, in step 404, the search parameter d i having the only peak is set as a delay time estimated value (^) d k-pre . This estimated delay time value (^) d k-pre is stored in the memory 6.

次に、ステップ405では、ステップ404で得られた遅延時間推定値で少タップ段数の適応フィルタ52を動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとし、これをメモリ6に記憶する。 Next, in step 405, operates the adaptive filter 52 k of small number of tap stages in the delay time estimation value obtained in step 404, using the adaptive algorithm, and calculates the load g k-pre, after calculation g k-pre is assumed to be a load estimated value (^) g k-pre and is stored in the memory 6.

次に、ステップ406では、遅延時間(^)dk−preと荷重(^)gk−preより誤差信号ek−pre(t)を求める。 Next, in step 406, an error signal e k-pre (t) is obtained from the delay time (^) d k-pre and the load (^) g k-pre .

そして、ステップ407において、ステップ406で求められた誤差信号の最小値が予め設定しておいたしきい値よりも小さくなる場合には、繰り返しの処理を止める。それ以外の場合にはk=k+1として、k=Kとなるまでステップ401〜406を繰り返す。ここまでで前処理は終了する。   In step 407, if the minimum value of the error signal obtained in step 406 is smaller than a preset threshold value, the iterative process is stopped. In other cases, k = k + 1 is set, and steps 401 to 406 are repeated until k = K. Thus, the preprocessing ends.

実際に干渉波を抑圧処理する場合には、図6中「干渉波抑圧時の処理」を行う。このとき、スイッチ3は送信信号発生器2側に切り替えて、遅延時間推定装置51、M段の適応フィルタ52、メモリ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 transmission signal generator 2 side and sends the second control signal (dotted line) to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6.

ステップ451では、複素タップ選択適応フィルタ5の遅延時間推定の処理と、タップ数M段の適応フィルタ52の荷重演算を行わず、予め広帯域信号を送信して求めた推定遅延時間(^)dk−preと推定荷重(^)gk−preをメモリ6から呼び出し、この遅延時間推定値(^)dk−preと荷重値(^)gk−preを用いて、誤差信号e(t)を求める。これにより、適応アルゴリズムの演算が不要になる。この行程は次式で表すことができる。 In step 451, the processing of the delay time estimation of the complex tap selection adaptive filter 5, without load calculation of the adaptive filter 52 k taps M stage, estimated delay time obtained by sending a pre-wideband signal (^) d k-pre and the estimated load (^) g k-pre are called from the memory 6, and the error signal e k (by using the delay time estimated value (^) d k-pre and the load value (^) g k-pre t). This eliminates the need for adaptive algorithm operations. This process can be expressed by the following equation.

Figure 0004545683
Figure 0004545683

ただし、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は、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、適応フィルタ52は、メモリ6から荷重g、gを呼び出して設定し、遅延時間推定装置51は、荷重g、gを用いてそれぞれ誤差信号の二乗平均値を求め、2つの二乗平均値から遅延時間評価関数を求め、この遅延時間評価関数から遅延時間を推定して推定遅延時間をメモリ6に記憶する。また、遅延時間推定装置51は、適応アルゴリズムを用いて推定遅延時間での荷重を推定して推定荷重をメモリ6に記憶する。次に、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の二乗平均値を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ401〜407である。 That is, when the switch 3 is switched to the wideband signal generator 1 side, the switch 3 sends a first control signal to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Based on the first control signal, the adaptive filter 52 k calls and sets the loads g and g * from the memory 6, and the delay time estimation device 51 k uses the loads g and g * , respectively. A mean square value is obtained, a delay time evaluation function is obtained from the two mean square values, a delay time is estimated from the delay time evaluation function, and the estimated delay time is stored in the memory 6. Further, the delay time estimation device 51 k estimates the load at the estimated delay time using an adaptive algorithm and stores the estimated load in the memory 6. Next, an error signal is obtained from the estimated delay time and estimated load. Until the minimum value of the error signal reaches the threshold value, the process from obtaining the first mean square value to the process obtaining the last error signal is repeated. These pre-processing are steps 401-407.

実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52は、メモリ6から推定荷重を呼び出して設定する。減算器9は、推定遅延時間と推定荷重(^)gから誤差信号を求める。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ451〜452である。 When the interference wave is actually suppressed, the switch 3 is switched to the transmission signal generator 2 side, and the second control signal is sent to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Send. Based on the second control signal, the delay time estimation device 51 k functions as a delay circuit by calling and setting the estimated delay time from the memory 6, and the adaptive filter 52 k calls the estimated weight from the memory 6. To set. The subtractor 9 k obtains an error signal from the estimated delay time and the estimated load (^) g k . Thereby, the interference wave is suppressed. These interference wave suppression processes are steps 451 to 452.

実施の形態5.
この発明の実施の形態5に係る干渉波除去装置について図7を参照しながら説明する。図7は、この発明の実施の形態5に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態5に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
Embodiment 5 FIG.
An interference wave removing apparatus according to Embodiment 5 of the present invention will be described with reference to FIG. FIG. 7 is a flowchart showing the operation of the interference wave canceller according to Embodiment 5 of the present invention. Note that the configuration of the interference wave canceller according to Embodiment 5 of the present invention is the same as that of Embodiment 1 above.

ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態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は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51、M段の適応フィルタ52、メモリ6が区別可能なように第1の制御信号(点線)を送る。ここでは、図7中の「予め広帯域信号で行う処理」を行う。 Consider a case where a broadband signal (white noise or the like) is transmitted from the broadband signal generator 1 in advance. At this time, the switch 3 switches to the wideband signal generator 1 side. Further, the switch 3 uses the first control signal (1) so that the delay time estimation device 51 1 , the M-stage adaptive filter 52 1 , and the memory 6 can distinguish that the delay time is estimated using the wideband signal. Send a dotted line). Here, the “processing performed in advance with a broadband signal” in FIG. 7 is performed.

ステップ501において、広帯域信号を送信した場合、広帯域信号x(t)と、この広帯域信号が未知系を通過した受信信号f(t)との相互相関関数rxfを計算する。これは次式で表すことができる。dは遅延時間探索パラメータとし、未知系のインパルス応答の時間軸を満たすように与える。 In step 501, when a broadband signal is transmitted, a cross-correlation function r xf between the broadband signal x (t) and the received signal f (t) through which the broadband signal has passed through the unknown system is calculated. This can be expressed as: d is a delay time search parameter and is given so as to satisfy the time axis of the impulse response of the unknown system.

Figure 0004545683
Figure 0004545683

ここで、〈〉内は期待値演算、*は複素共役をそれぞれ表す。   Here, <> represents an expected value calculation, and * represents a complex conjugate.

次に、ステップ502において、式(10)はd=dのみに相関値のピークが現れる。この唯一のピークを与える探索パラメータdを遅延時間推定値(^)d1−preとする。この遅延時間推定値(^)d1−preをメモリ6に記憶する。ここまでで前処理は終了する。 Next, in step 502, the peak of the correlation value appears only in d 1 = d i in equation (10). Let the search parameter d i giving this unique peak be the delay time estimate (^) d 1 -pre . This estimated delay time value (^) d 1-pre is stored in the memory 6. Thus, the preprocessing ends.

以下、ステップ551〜554の処理は、上記の実施の形態1のステップ151〜154とそれぞれ同様であるので説明は省略する。   Hereinafter, the processing in steps 551 to 554 is the same as that in steps 151 to 154 in the first embodiment, and a description thereof will be omitted.

すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、遅延時間推定装置51は、広帯域信号x(t)と、この広帯域信号が未知系を通過してきた受信信号f(t)の相互相関関数を計算し、相互相関値が最大となる遅延時間を推定遅延時間としてメモリ6に記憶する。これらの前処理がステップ501〜502である。 That is, when the switch 3 is switched to the wideband signal generator 1 side, the switch 3 sends a first control signal to the delay time estimation device 51 1 , the M-stage adaptive filter 52 1 , and the memory 6. Based on this first control signal, the delay time estimating unit 51 1 calculates the wideband signal x (t), the cross-correlation function of the broadband signal is a received signal f that has passed through the unknown system (t), The delay time that maximizes the cross-correlation value is stored in the memory 6 as the estimated delay time. These pre-processing are steps 501 to 502.

実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52は、適応アルゴリズムを用いて荷重gを求め、減算器9は、推定遅延時間と荷重gから誤差信号を求め、これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ551〜554である。 When the interference wave is actually suppressed, the switch 3 is switched to the transmission signal generator 2 side, and the second control signal is sent to the delay time estimation device 51 1 , the M-stage adaptive filter 52 1 , and the memory 6. Send. Based on this second control signal, the delay time estimating unit 51 1 functions as a delay circuit by setting calls the estimated delay time from the memory 6, the adaptive filter 52 1, a load g using an adaptive algorithm 1 determined, subtractor 9 1 obtains an error signal from the estimated delay time and the load g 1, which interference wave is suppressed by. These interference wave suppression processes are steps 551 to 554.

実施の形態6.
この発明の実施の形態6に係る干渉波除去装置について図8を参照しながら説明する。図8は、この発明の実施の形態6に係る干渉波除去装置の動作を示すフローチャートである。なお、この発明の実施の形態6に係る干渉波除去装置の構成は、上記の実施の形態1と同様である。
Embodiment 6 FIG.
An interference wave canceller according to Embodiment 6 of the present invention will be described with reference to FIG. FIG. 8 is a flowchart showing the operation of the interference wave canceller according to Embodiment 6 of the present invention. Note that the configuration of the interference wave canceller according to Embodiment 6 of the present invention is the same as that of Embodiment 1 above.

ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態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は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51、M段の適応フィルタ52、メモリ6が区別可能なように第1の制御信号(点線)を送る。この前処理では、図8中の「予め広帯域信号で行う処理」を行う。 Consider a case where a broadband signal (white noise or the like) is transmitted from the broadband signal generator 1 in advance. At this time, the switch 3 switches to the wideband signal generator 1 side. In addition, the switch 3 uses the first control signal ( k ) so that the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6 can distinguish that the delay time is estimated using the wideband signal. Send a dotted line). In this pre-processing, “processing performed in advance with a broadband signal” in FIG. 8 is performed.

ステップ601において、広帯域信号を送信し、広帯域信号と、この広帯域信号が未知系を通過した受信信号との相互相関関数rxfを計算する。 In step 601, a wideband signal is transmitted, and a cross-correlation function r xf between the wideband signal and a received signal that has passed through the unknown system is calculated.

次に、ステップ602において、相互相関関数rxfはd=dのみに相関値のピークが現れる。この唯一のピークをもつ探索パラメータdを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。 Next, in step 602, the cross-correlation function r xf has a correlation value peak only at d = d k . The search parameter d k having the only peak is assumed to be a delay time estimated value (^) d k-pre . This estimated delay time value (^) d k-pre is stored in the memory 6.

次に、ステップ603では、ステップ602で得られた遅延時間推定値(^)dk−preで少タップ段数の適応フィルタ52を動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとする。 Next, in step 603, the adaptive filter 52 k having a small number of taps is operated with the delay time estimated value (^) d k-pre obtained in step 602, and the weight g k-pre is calculated using an adaptive algorithm. The calculated g k-pre is defined as a load estimated value (^) g k-pre .

次に、ステップ604では、遅延時間推定値(^)dk−preと荷重(^)gk−preを用いて誤差信号ek−pre(t)を次式により計算する。 Next, in step 604, 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 .

Figure 0004545683
Figure 0004545683

ただし、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 step 605, if the minimum value of the error signal obtained in step 604 is smaller than a preset threshold value, the iterative process is stopped. In other cases, k = k + 1 is set, and steps 601 to 604 are repeated until k = K. Thus, the preprocessing ends.

以下、ステップ651〜654の処理は、上記の実施の形態2のステップ251〜254とそれぞれ同様であるので説明は省略する。   Hereinafter, the processing in steps 651 to 654 is the same as that in steps 251 to 254 in the second embodiment, and a description thereof will be omitted.

すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、遅延時間推定装置51は、広帯域信号x(t)と、この広帯域信号が未知系を通過してきた受信信号f(t)の相互相関関数を計算し、相互相関値が最大となる遅延時間を推定遅延時間としてメモリ6に記憶する。また、遅延時間推定装置51は、適応アルゴリズムを用いて推定遅延時間での荷重を推定し、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の相互相関関数を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ601〜605である。 That is, when the switch 3 is switched to the wideband signal generator 1 side, the switch 3 sends a first control signal to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Based on the first control signal, the delay time estimation device 51 k calculates a cross-correlation function between the wideband signal x (t) and the received signal f (t) through which the wideband signal has passed through the unknown system, The delay time that maximizes the cross-correlation value is stored in the memory 6 as the estimated delay time. Further, the delay time estimation device 51 k estimates a load at the estimated delay time using an adaptive algorithm, and obtains an error signal from the estimated delay time and the estimated load. The process from obtaining the first cross-correlation function to the process obtaining the last error signal is repeated until the minimum value of the error signal reaches the threshold value. These pre-processing are steps 601 to 605.

実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52は、適応アルゴリズムを用いて荷重gを求め、減算器9は、推定遅延時間と推定荷重(^)gから誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の荷重を求める処理から最後の誤差信号を求める処理までを繰り返す。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ651〜654である。 When the interference wave is actually suppressed, the switch 3 is switched to the transmission signal generator 2 side, and the second control signal is sent to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Send. Based on the second control signal, the delay time estimation device 51 k functions as a delay circuit by calling and setting the estimated delay time from the memory 6, and the adaptive filter 52 k uses the adaptive algorithm to load g k is obtained, and the subtractor 9 k obtains an error signal from the estimated delay time and the estimated weight (^) g k . The process from obtaining the first load to the process obtaining the last error signal is repeated until the minimum value of the error signal reaches the threshold value. Thereby, the interference wave is suppressed. These interference wave suppression processes are steps 651 to 654.

実施の形態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 Embodiment 1 above.

ここでは簡単のためにプラントノイズは十分小さい状況により説明する。この実施の形態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は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51、M段の適応フィルタ52、メモリ6が区別可能なように第1の制御信号(点線)を送る。前処理では、図9中の「予め広帯域信号で行う処理」を行う。 Consider a case where a broadband signal (white noise or the like) is transmitted from the broadband signal generator 1 in advance. At this time, the switch 3 switches to the wideband signal generator 1 side. In addition, the switch 3 uses the first control signal ( k ) so that the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6 can distinguish that the delay time is estimated using the wideband signal. Send a dotted line). In the pre-processing, “processing performed in advance with a broadband signal” in FIG. 9 is performed.

ステップ701において、広帯域信号を送信し、広帯域信号と、この広帯域信号が未知系を通過した受信信号の相互相関関数rxfを計算する。 In step 701, a wideband signal is transmitted, and a cross-correlation function r xf of the wideband signal and a received signal that has passed through the unknown system is calculated.

次に、ステップ702において、相互相関関数rxfはd=dのみに相関値のピークが現れる。この唯一のピークをもつ探索パラメータdを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。 Next, in step 702, the cross-correlation function r xf has a correlation value peak only at d = d k . The search parameter d having the only peak is assumed to be a delay time estimated value (^) d k-pre . This estimated delay time value (^) d k-pre is stored in the memory 6.

次に、ステップ703では、ステップ702で得られた遅延時間推定値(^)dk−preで少タップ段数の適応フィルタを動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとし、これをメモリ6に記憶する。 Next, in step 703, an adaptive filter with a small number of taps is operated with the delay time estimated value (^) d k-pre obtained in step 702, and the load g k-pre is calculated using an adaptive algorithm. Then, g k-pre after the calculation is set as a load estimated value (^) g k-pre and stored in the memory 6.

次に、ステップ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 .

Figure 0004545683
Figure 0004545683

ただし、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 step 705, if the minimum value of the error signal obtained in step 704 is smaller than a preset threshold value, the iterative process is stopped. In other cases, k = k + 1 is set, and steps 701 to 704 are repeated until k = K. Thus, the preprocessing ends.

以下、ステップ751〜754の処理は、上記の実施の形態3のステップ351〜354とそれぞれ同様であるので説明は省略する。   Hereinafter, the processing in steps 751 to 754 is the same as that in steps 351 to 354 in the third embodiment, and a description thereof will be omitted.

すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、遅延時間推定装置51は、広帯域信号x(t)と、この広帯域信号が未知系を通過してきた受信信号f(t)の相互相関関数を計算し、相互相関値が最大となる遅延時間を推定遅延時間としてメモリ6に記憶する。また、遅延時間推定装置51は、適応アルゴリズムを用いて推定遅延時間での荷重を推定して推定荷重をメモリ6に記憶する。次に、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の相互相関関数を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ701〜705である。 That is, when the switch 3 is switched to the wideband signal generator 1 side, the switch 3 sends a first control signal to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Based on the first control signal, the delay time estimation device 51 k calculates a cross-correlation function between the wideband signal x (t) and the received signal f (t) through which the wideband signal has passed through the unknown system, The delay time that maximizes the cross-correlation value is stored in the memory 6 as the estimated delay time. Further, the delay time estimation device 51 k estimates the load at the estimated delay time using an adaptive algorithm and stores the estimated load in the memory 6. Next, an error signal is obtained from the estimated delay time and estimated load. The process from obtaining the first cross-correlation function to the process obtaining the last error signal is repeated until the minimum value of the error signal reaches the threshold value. These preprocessing are steps 701 to 705.

実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52は、メモリ6から推定荷重を呼び出して初期値として設定し、適応アルゴリズムを用いて荷重gを推定し、減算器9は、推定遅延時間と推定荷重(^)gから誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の荷重を求める処理から最後の誤差信号を求める処理までを繰り返す。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ751〜754である。 When the interference wave is actually suppressed, the switch 3 is switched to the transmission signal generator 2 side, and the second control signal is sent to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Send. Based on the second control signal, the delay time estimation device 51 k functions as a delay circuit by calling and setting the estimated delay time from the memory 6, and the adaptive filter 52 k calls the estimated weight from the memory 6. As an initial value, the load g k is estimated using an adaptive algorithm, and the subtractor 9 k obtains an error signal from the estimated delay time and the estimated load (^) g k . The process from obtaining the first load to the process obtaining the last error signal is repeated until the minimum value of the error signal reaches the threshold value. Thereby, the interference wave is suppressed. These interference wave suppression processes are steps 751 to 754.

実施の形態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 Embodiment 1 above.

ここでは簡単のためにプラントノイズは十分小さい状況により説明する。上記の実施の形態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 memory 6 and used as an initial value at the time of calculating the k-th stage load when suppressing the interference wave, in the eighth embodiment, the load value obtained in advance during the processing of the wideband signal is used for interference suppression. Sometimes used as a load value.

予め広帯域信号発生器1より広帯域信号(白色ノイズ等)を送信する場合を考える。このとき、スイッチ3は広帯域信号発生器1側に切り替える。また、スイッチ3は、広帯域信号を用いて遅延時間の推定を行っていることを遅延時間推定装置51、M段の適応フィルタ52、メモリ6が区別可能なように第1の制御信号(点線)を送る。前処理では図10中の「予め広帯域信号で行う処理」を行う。 Consider a case where a broadband signal (white noise or the like) is transmitted from the broadband signal generator 1 in advance. At this time, the switch 3 switches to the wideband signal generator 1 side. In addition, the switch 3 uses the first control signal ( k ) so that the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6 can distinguish that the delay time is estimated using the wideband signal. Send a dotted line). In the pre-processing, “processing performed in advance with a broadband signal” in FIG. 10 is performed.

ステップ801において、広帯域信号を送信し、広帯域信号と、この広帯域信号が未知系を通過した受信信号との相互相関関数rxfを計算する。 In step 801, a wideband signal is transmitted, and a cross-correlation function r xf between the wideband signal and a received signal that has passed through the unknown system is calculated.

次に、ステップ802において、相互相関関数rxfはd=dのみに相関値のピークが現れる。この唯一のピークをもつ探索パラメータdを遅延時間推定値(^)dk−preとする。この遅延時間推定値(^)dk−preをメモリ6に記憶する。 Next, in step 802, the cross-correlation function r xf has a correlation value peak only at d = d k . The search parameter d having the only peak is assumed to be a delay time estimated value (^) d k-pre . This estimated delay time value (^) d k-pre is stored in the memory 6.

次に、ステップ803では、ステップ802で得られた遅延時間推定値(^)dk−preで少タップ段数の適応フィルタ52を動作させ、適応アルゴリズムを用いて、その荷重gk−preを計算し、計算後のgk−preを荷重推定値(^)gk−preとし、これをメモリ6に記憶する。 Next, in step 803, the adaptive filter 52 k having a small number of tap stages is operated with the delay time estimation value (^) d k-pre obtained in step 802, and the weight g k-pre is calculated using an adaptive algorithm. The calculated g k-pre is set as a load estimated value (^) g k-pre and stored in the memory 6.

次に、ステップ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 .

Figure 0004545683
Figure 0004545683

ただし、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 step 805, if the minimum value of the error signal obtained in step 804 is smaller than a preset threshold value, the iterative process is stopped. In other cases, k = k + 1 is set, and steps 801 to 804 are repeated until k = K. Thus, the preprocessing ends.

以下、ステップ851〜852の処理は、上記の実施の形態4のステップ451〜452とそれぞれ同様であるので説明は省略する。   Hereinafter, the processing in steps 851 to 852 is the same as that in steps 451 to 452 of the above-described fourth embodiment, and thus description thereof is omitted.

すなわち、スイッチ3が広帯域信号発生器1側に切り替えられると、スイッチ3は、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6に第1の制御信号を送る。この第1の制御信号に基づいて、遅延時間推定装置51は、広帯域信号x(t)と、この広帯域信号が未知系を通過してきた受信信号f(t)の相互相関関数を計算し、相互相関値が最大となる遅延時間を推定遅延時間としてメモリ6に記憶する。また、遅延時間推定装置51は、適応アルゴリズムを用いて推定遅延時間での荷重を推定して推定荷重をメモリ6に記憶する。次に、推定遅延時間及び推定荷重より誤差信号を求める。この誤差信号の最小値がしきい値に到達するまで最初の相互相関関数を求める処理から最後の誤差信号を求める処理までを繰り返す。これらの前処理がステップ801〜805である。 That is, when the switch 3 is switched to the wideband signal generator 1 side, the switch 3 sends a first control signal to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Based on the first control signal, the delay time estimation device 51 k calculates a cross-correlation function between the wideband signal x (t) and the received signal f (t) through which the wideband signal has passed through the unknown system, The delay time that maximizes the cross-correlation value is stored in the memory 6 as the estimated delay time. Further, the delay time estimation device 51 k estimates the load at the estimated delay time using an adaptive algorithm and stores the estimated load in the memory 6. Next, an error signal is obtained from the estimated delay time and estimated load. The process from obtaining the first cross-correlation function to the process obtaining the last error signal is repeated until the minimum value of the error signal reaches the threshold value. These pre-processing are steps 801 to 805.

実際に干渉波を抑圧処理する場合には、スイッチ3が送信信号発生器2側に切り替えられて、遅延時間推定装置51、M段の適応フィルタ52、及びメモリ6へ第2の制御信号を送る。この第2の制御信号に基づいて、遅延時間推定装置51は、メモリ6から推定遅延時間を呼び出して設定することにより遅延回路として機能し、適応フィルタ52は、メモリ6から推定荷重を呼び出して設定する。減算器9は、推定遅延時間と推定荷重(^)gから誤差信号を求める。これにより干渉波が抑圧される。これらの干渉波抑圧処理がステップ851〜852である。 When the interference wave is actually suppressed, the switch 3 is switched to the transmission signal generator 2 side, and the second control signal is sent to the delay time estimation device 51 k , the M-stage adaptive filter 52 k , and the memory 6. Send. Based on the second control signal, the delay time estimation device 51 k functions as a delay circuit by calling and setting the estimated delay time from the memory 6, and the adaptive filter 52 k calls the estimated weight from the memory 6. To set. The subtractor 9 k obtains an error signal from the estimated delay time and the estimated load (^) g k . Thereby, the interference wave is suppressed. These interference wave suppression processes are steps 851 to 852.

この発明の実施の形態1に係る干渉波除去装置の構成を示す図である。It is a figure which shows the structure of the interference wave removal apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る干渉波除去装置の複素タップ選択適応フィルタの構成を示す図である。It is a figure which shows the structure of the complex tap selection adaptive filter of the interference wave removal apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る干渉波除去装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the interference wave removal apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態2に係る干渉波除去装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the interference wave removal apparatus which concerns on Embodiment 2 of this invention. この発明の実施の形態3に係る干渉波除去装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the interference wave removal apparatus which concerns on Embodiment 3 of this invention. この発明の実施の形態4に係る干渉波除去装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the interference wave removal apparatus which concerns on Embodiment 4 of this invention. この発明の実施の形態5に係る干渉波除去装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the interference wave removal apparatus which concerns on Embodiment 5 of this invention. この発明の実施の形態6に係る干渉波除去装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the interference wave removal apparatus which concerns on Embodiment 6 of this invention. この発明の実施の形態7に係る干渉波除去装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the interference wave removal apparatus which concerns on Embodiment 7 of this invention. この発明の実施の形態8に係る干渉波除去装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the interference wave removal apparatus which concerns on Embodiment 8 of this invention. 従来の遅延時間推定装置の原理的構成を示す概念図である。It is a conceptual diagram which shows the fundamental structure of the conventional delay time estimation apparatus. 従来のタップ選択適応フィルタの動作を示すフローチャートである。It is a flowchart which shows operation | movement of the conventional tap selection adaptive filter. 広帯域信号における遅延時間推定の手法を示す概念図である。It is a conceptual diagram which shows the method of delay time estimation in a wideband signal. 周期性をもつ狭帯域信号における遅延時間推定の手法を示す概念図である。It is a conceptual diagram which shows the method of delay time estimation in the narrow-band signal with periodicity.

符号の説明Explanation of symbols

1 広帯域信号発生器、2 送信信号発生器、3 スイッチ、4 送信アンテナ、5 複素タップ選択適応フィルタ、6 メモリ、7 受信アンテナ、8 加算器、9、9、9 減算器、51、51 遅延時間推定装置、52、52 適応フィルタ。 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 subtractor 51 1, 51 k delay time estimation device, 52 1 , 52 k adaptive filter.

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.
JP2005354522A 2005-12-08 2005-12-08 Interference wave removing apparatus and interference wave removing method Expired - Fee Related JP4545683B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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