JP3183128B2 - Signal to noise power ratio improving apparatus and signal to noise power ratio improving method - Google Patents
Signal to noise power ratio improving apparatus and signal to noise power ratio improving methodInfo
- Publication number
- JP3183128B2 JP3183128B2 JP27563895A JP27563895A JP3183128B2 JP 3183128 B2 JP3183128 B2 JP 3183128B2 JP 27563895 A JP27563895 A JP 27563895A JP 27563895 A JP27563895 A JP 27563895A JP 3183128 B2 JP3183128 B2 JP 3183128B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- input signal
- power ratio
- cumulant
- parameters
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Noise Elimination (AREA)
- Filters That Use Time-Delay Elements (AREA)
Description
【0001】[0001]
【発明の属する技術分野】本発明は、振幅分布がガウス
分布である雑音(ガウス性雑音)を含む周期性信号の信
号対雑音電力比を改善する装置並びに方法に関するもの
である。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an apparatus and a method for improving a signal-to-noise power ratio of a periodic signal including noise whose amplitude distribution is Gaussian (Gaussian noise).
【0002】[0002]
【従来の技術】図10は例えば、United States Patent
5091890、■Method of Extracting Target Range And
Doppler Information From a Doppler Spread Signal■
に記載されている4次キュムラントを用いた信号対雑音
電力比改善装置の構成図である。図10において、1は
シフトレジスタ、2は4次キュムラント計算回路、3は
窓関数回路、4は離散フーリエ変換を行うものであり、
ここでは高速フーリエ変換回路を用いている。5は包絡
線検波を行うものであり、ここでは自乗検波回路を用い
ている。6は信号対雑音電力比計算回路である。2. Description of the Related Art FIG.
5091890, ■ Method of Extracting Target Range And
Doppler Information From a Doppler Spread Signal ■
FIG. 2 is a configuration diagram of a signal-to-noise power ratio improvement device using a fourth-order cumulant described in FIG. In FIG. 10, reference numeral 1 denotes a shift register, 2 denotes a fourth-order cumulant calculation circuit, 3 denotes a window function circuit, and 4 denotes a discrete Fourier transform.
Here, a fast Fourier transform circuit is used. Numeral 5 is for performing envelope detection. Here, a square detection circuit is used. Reference numeral 6 denotes a signal-to-noise power ratio calculation circuit.
【0003】4次キュムラントとは、高次統計量(High
er-Order Statistics)の一種であり、その定義およ
び、特性については、例えばC.L.Nikias and J.M.Mende
l著、■Signal Processing with Higher-Order Spectra
■、IEEE Signal ProcessingMagazine,pp10-37,July 19
93.に詳しく記載されている。そのため、ここでは本発
明に利用する4次キュムラントの定義と特性だけを以下
に簡単に述べる。A fourth-order cumulant is a high-order statistic (High
er-Order Statistics), whose definition and characteristics are described, for example, in CLNikias and JMMende
l, ■ Signal Processing with Higher-Order Spectra
■ 、 IEEE Signal ProcessingMagazine, pp10-37, July 19
93. Therefore, only the definition and characteristics of the quaternary cumulant used in the present invention will be briefly described below.
【0004】平均がゼロである信号に対する4次キュム
ラントc4 xは、式(1)に示すように、3つのタイムラグ
(t1、t2、t3)の関数で定義される。A fourth-order cumulant c 4 x for a signal having an average of zero is defined by a function of three time lags (t1, t2, t3) as shown in equation (1).
【0005】[0005]
【数1】 (Equation 1)
【0006】ただし、m4 x、m2 xはそれぞれ4次モーメ
ント関数、2次モーメント関数(自己相関関数)を表
し、以下の式で表される。Here, m 4 x and m 2 x represent a fourth moment function and a second moment function (autocorrelation function), respectively, and are represented by the following equations.
【0007】[0007]
【数2】 (Equation 2)
【0008】式(1)から式(8)において、共役複素数
(*)をとるのは、4つの入力信号の値(x(n)、x(n+t
1)、x(n+t2)、x(n+t3))の中で2カ所であれば、どこで
もかまわない。なお、入力が実信号である場合は、共役
複素数(*)をとる必要はない。In Equations (1) to (8), the conjugate complex number (*) is obtained because the values of four input signals (x (n), x (n + t)
1), x (n + t2), x (n + t3)) as long as there are two locations. When the input is a real signal, it is not necessary to take a conjugate complex number (*).
【0009】また、高次キュムラントには以下に示す特
長がある。 (1)振幅分布がガウス分布である信号に対する3次以
上のキュムラントが理論上ゼロになる。 (2)離散信号x(n)と離散信号y(n)が独立ならば、次式
に示す線形な関係がある。[0009] The high-order cumulant has the following features. (1) Cumulants of third or higher order for a signal whose amplitude distribution is Gaussian distribution are theoretically zero. (2) If the discrete signal x (n) and the discrete signal y (n) are independent, there is a linear relationship represented by the following equation.
【0010】[0010]
【数3】 (Equation 3)
【0011】(3)入力信号が式(10)に示すような複素
正弦波信号(Aは振幅、φは初期位相、ωは角周波数)
の場合の4次キュムラント関数は、式(11)に示すに、同
じ周波数成分ωを持つt1、t2、t3の関数の複素正弦波と
なる。(3) An input signal is a complex sine wave signal as shown in equation (10) (A is amplitude, φ is initial phase, ω is angular frequency)
Is a complex sine wave of a function of t1, t2, and t3 having the same frequency component ω as shown in Expression (11).
【0012】[0012]
【数4】 (Equation 4)
【0013】また、式(11)は、例えば、t1=t、t2=t3=0
とすることによって、式(12)に示すように、変数を1次
元化することができる。Equation (11) is, for example, t1 = t, t2 = t3 = 0
By doing so, the variables can be made one-dimensional as shown in Expression (12).
【0014】[0014]
【数5】 (Equation 5)
【0015】これらの特性により、ある所望信号に、所
望信号と独立なガウス性雑音が加わった信号に対して、
4次キュムラントを求めると、理論的には、ガウス性雑
音成分がゼロとなり、所望信号に対する4次キュムラン
トのみが残り、高い信号対雑音電力比改善効果(理論的
には無限大)が得られる。特に、所望信号が周期性信号
の場合には、その周波数成分の情報も、4次キュムラン
トに保存される。Due to these characteristics, a signal obtained by adding a Gaussian noise independent of a desired signal to a certain desired signal,
When the fourth-order cumulant is obtained, the Gaussian noise component becomes zero theoretically, only the fourth-order cumulant for the desired signal remains, and a high signal-to-noise power ratio improvement effect (theoretical infinity) is obtained. In particular, when the desired signal is a periodic signal, information on its frequency component is also stored in the fourth cumulant.
【0016】次に、上記従来例の信号対雑音電力比改善
装置の動作を図10に基づいて説明する。説明を簡単に
するために、ここでは入力信号x(n)を式(13)に示すよう
に、振幅A、角周波数ω、初期位相φの複素正弦波信号
に、その複素正弦波信号と独立で、平均ゼロ、分散σ2
のガウス性雑音g(n)が加わった信号をある時間間隔でサ
ンプリングした離散信号とする。Next, the operation of the conventional signal-to-noise power ratio improving apparatus will be described with reference to FIG. For the sake of simplicity, here, the input signal x (n) is converted into a complex sine wave signal having an amplitude A, an angular frequency ω, and an initial phase φ as shown in Expression (13). With zero mean and variance σ 2
The signal to which the Gaussian noise g (n) is added is a discrete signal sampled at a certain time interval.
【0017】[0017]
【数6】 (Equation 6)
【0018】また、式(13)で表される、入力信号の信号
対雑音電力比SNRinは、式(14)で表される。The signal-to-noise power ratio SNRin of the input signal, which is represented by the equation (13), is represented by the equation (14).
【0019】[0019]
【数7】 (Equation 7)
【0020】この入力信号x(n)は、シフトレジスタ1に
順次入力される。シフトレジスタ1では、Nサンプルの
入力信号を保持し、それらを同時に4次キュムラント計
算回路2に入力する。4次キュムラント計算回路2で
は、式(1)で示した3変数t1、t2、t3の4次キュムラン
ト関数において、例えば、t1=t、t2=t3=0(先に述べ
た、United States Patent 5091890では、t1=0、t2=t3=
tとしている。)とし、tを最大-(N-1)から(N-1)まで変
化させて求めた最大2N-1個の4次キュムラントを窓関数
回路3に出力する。ただし、式(11)にも示したように、
4次キュムラントの振幅は、入力信号の振幅の4乗とな
るために、窓関数回路3へは、振幅の4乗根をとった4
次キュムラントを出力する。窓関数回路3では、高速フ
ーリエ変換回路4における、4次キュムラント計算回路
2で求めた4次キュムラントの数が有限であることの影
響を低減するために、窓関数の重みを掛け、その結果
を、高速フーリエ変換回路4に出力する。高速フーリエ
変換回路4では、窓関数回路3の出力信号を、4次キュ
ムラント計算回路2で求めた4次キュムラントの数と入
力信号x(n)のサンプリング時間で決まる帯域ごとの成分
に分別し、自乗検波回路5に出力する。自乗検波回路5
では、各帯域の信号の電力を求め、信号対雑音電力比計
算回路6に出力する。信号対雑音電力比計算回路6で
は、入力信号x(n)の複素正弦波信号の周波数が分かって
いる、あるいは予測している場合、周波数の存在する帯
域の電力値から、周波数の存在する帯域以外の電力の平
均値を引いた値と、周波数の存在する帯域以外の電力の
平均値との比を求める。また、入力信号x(n)の複素正弦
波信号の周波数が分からない場合は、最も大きい電力値
からその他の帯域の電力の平均値を引いた値と、最も大
きい電力値を持つ帯域以外の帯域の電力の平均値との比
を求める。The input signal x (n) is sequentially input to the shift register 1. The shift register 1 holds the input signals of N samples and simultaneously inputs them to the fourth cumulant calculation circuit 2. In the fourth-order cumulant calculation circuit 2, in the fourth-order cumulant function of the three variables t1, t2, and t3 shown in the equation (1), for example, t1 = t, t2 = t3 = 0 (the aforementioned United States Patent 5091890) Then, t1 = 0, t2 = t3 =
and t. ), And outputs a maximum of 2N-1 fourth-order cumulants obtained by changing t from-(N-1) to (N-1) to the window function circuit 3. However, as shown in equation (11),
Since the amplitude of the fourth-order cumulant is the fourth power of the amplitude of the input signal, a fourth root of the amplitude is given to the window function circuit 3.
Output the next cumulant. The window function circuit 3 multiplies the weight of the window function in order to reduce the effect of the finite number of fourth-order cumulants obtained by the fourth-order cumulant calculation circuit 2 in the fast Fourier transform circuit 4, and multiplies the result by , To the fast Fourier transform circuit 4. The fast Fourier transform circuit 4 separates the output signal of the window function circuit 3 into components for each band determined by the number of fourth-order cumulants calculated by the fourth-order cumulant calculation circuit 2 and the sampling time of the input signal x (n), Output to the square detection circuit 5. Square detection circuit 5
Then, the power of the signal in each band is obtained and output to the signal-to-noise power ratio calculation circuit 6. In the signal-to-noise power ratio calculation circuit 6, when the frequency of the complex sine wave signal of the input signal x (n) is known or predicted, the band in which the frequency exists is calculated from the power value of the band in which the frequency exists. Then, a ratio of a value obtained by subtracting the average value of the powers other than the frequency and the average value of the power other than the frequency band is calculated. If the frequency of the complex sine wave signal of the input signal x (n) is not known, a value obtained by subtracting the average value of the power of the other bands from the largest power value and a band other than the band having the largest power value And the ratio of the average power to the average value.
【0021】[0021]
【発明が解決しようとする課題】上記のような信号対雑
音電力比改善装置では、入力信号のサンプル数が有限個
であるため、4次キュムラント計算回路2において、ガ
ウス性雑音に対する4次キュムラントが完全にゼロにな
らず、消え残り信号が生じる。その結果、信号対雑音電
力比改善効果が低くなるという問題点があった。In the signal-to-noise power ratio improving apparatus as described above, since the number of samples of the input signal is finite, the fourth-order cumulant for the Gaussian noise is reduced in the fourth-order cumulant calculation circuit 2. It does not become completely zero, and a remaining signal is generated. As a result, there is a problem that the signal-to-noise power ratio improvement effect is reduced.
【0022】本発明は、かかる問題を解決するためにな
されたものであり、ガウス性雑音に対する4次キュムラ
ントが完全にゼロにならず、消え残り信号があっても、
その消え残り信号をさらに小さくすることによって、従
来の装置よりも高い信号対雑音電力比改善効果を得るこ
とができる信号対雑音電力比改善装置および信号対雑音
電力比改善方法を提案する。The present invention has been made to solve such a problem, and the fourth-order cumulant for Gaussian noise is not completely zero, and even if there is a signal that remains undisturbed,
A signal-to-noise power ratio improvement device and a signal-to-noise power ratio improvement method capable of obtaining a higher signal-to-noise power ratio improvement effect than a conventional device by further reducing the unerased signal are proposed.
【0023】[0023]
【課題を解決するための手段】請求項1の発明に係わる
信号対雑音電力比改善装置は、入力信号のガウス性雑音
成分を抑圧する信号対雑音電力比改善装置において、入
力信号をサンプリングして一組の入力信号を生成する入
力信号取り込み手段と、上記入力信号取り込み手段から
の同一の一組の入力信号について、4次キュムラント関
数の3変数の中の1つを変数、他の2つをそれぞれパラ
メータとし、上記2つのパラメータを入力信号中の所望
信号に対する4次キュムラント関数の初期位相がすべて
同じになるようにした複数個の異なる2つのパラメータ
の組み合わせのそれぞれに対し、上記変数を変化させて
4次キュムラントの組を求める4次キュムラント計算処
理手段と、上記4次キュムラント計算処理手段からの4
次キュムラントの組のそれぞれについて、同じ変数ごと
に加算処理してその平均を求める加算処理手段とを備え
たものである。According to a first aspect of the present invention, there is provided a signal-to-noise power ratio improving apparatus for suppressing a Gaussian noise component of an input signal by sampling an input signal. For the input signal capturing means for generating a set of input signals, and for the same set of input signals from the input signal capturing means, one of the three variables of the fourth-order cumulant function is a variable and the other two are variables. The above parameters are varied for each of a plurality of different two parameter combinations in which the initial parameters of the fourth order cumulant function for the desired signal in the input signal are all the same. A fourth cumulant calculation processing means for obtaining a set of fourth cumulants by means of
An addition processing unit is provided for obtaining an average of the next cumulant set by adding the same variables.
【0024】請求項2の発明に係わる信号対雑音電力比
改善方法は、入力信号のガウス性雑音成分を抑圧する信
号対雑音電力比改善方法において、入力信号をサンプリ
ングして一組の入力信号を生成し、同一の上記一組の入
力信号について、4次キュムラント関数の3変数の中の
1つを変数、他の2つをそれぞれパラメータとし、上記
2つのパラメータを入力信号中の所望信号に対する4次
キュムラント関数の初期位相がすべて同じになるように
した複数個の異なる2つのパラメータの組み合わせのそ
れぞれに対し、上記変数を変化させて4次キュムラント
の組を求め、上記4次キュムラントの組のそれぞれにつ
いて、同じ変数ごとに加算処理してその平均を求めるも
のである。According to a second aspect of the present invention, there is provided a signal-to-noise power ratio improving method for suppressing a Gaussian noise component of an input signal. For the same set of input signals to be generated, one of the three variables of the fourth-order cumulant function is a variable, and the other two are parameters, respectively. For each of a plurality of different combinations of two parameters in which the initial phase of the second order cumulant function is all the same, the above variables are changed to obtain a set of fourth order cumulants. Are added to each of the same variables, and the average is obtained.
【0025】請求項3の発明に係わる信号対雑音電力比
改善装置は、入力信号のガウス性雑音成分を抑圧する信
号対雑音電力比改善装置において、入力信号をサンプリ
ングして一組の入力信号を生成する入力信号取り込み手
段と、上記入力信号取り込み手段からの同一の一組の入
力信号について、4次キュムラント関数の3変数の中の
1つを変数、他の2つをそれぞれパラメータとし、複数
個の異なる2つのパラメータの組み合わせのそれぞれに
対し、上記変数を変化させて4次キュムラントの組を求
める4次キュムラント計算処理手段と、上記4次キュム
ラント計算処理手段からの4次キュムラントの組のそれ
ぞれについて、各4次キュムラントに各組の上記2つの
パラメータに基づいて得られる量の関数シフトを与える
関数シフト処理手段と、上記関数シフト処理手段からの
出力を同じ変数ごとに加算処理してその平均を求める加
算処理手段とを備えたものである。According to a third aspect of the present invention, there is provided a signal to noise power ratio improving apparatus for suppressing a Gaussian noise component of an input signal. For the input signal capturing means to be generated and the same set of input signals from the input signal capturing means, one of the three variables of the fourth-order cumulant function is used as a variable and the other two are used as parameters, respectively. For each of the combinations of the two different parameters, the above variables are changed to obtain a fourth order cumulant set, and the fourth order cumulant set from the fourth order cumulant calculation processing means Function shift processing means for giving each fourth order cumulant a function shift of an amount obtained based on each set of the above two parameters. When, in which an adding processing means for obtaining an average by adding processing the output for each same variable from the function shift processing means.
【0026】請求項4の発明に係わる信号対雑音電力比
改善方法は、入力信号のガウス性雑音成分を抑圧する信
号対雑音電力比改善方法において、入力信号をサンプリ
ングして一組の入力信号を生成し、同一の上記一組の入
力信号について、4次キュムラント関数の3変数の中の
1つを変数、他の2つをそれぞれパラメータとし、複数
個の異なる2つのパラメータの組み合わせのそれぞれに
対し、上記変数を変化させて4次キュムラントの組を求
め、上記4次キュムラントの組のそれぞれについて、各
4次キュムラントに各組の上記2つのパラメータに基づ
いて得られる量の関数シフトを与え、同じ変数ごとに加
算処理してその平均を求めるものである。According to a fourth aspect of the present invention, there is provided a signal-to-noise power ratio improving method for suppressing a Gaussian noise component of an input signal. For the same set of input signals to be generated, one of the three variables of the fourth-order cumulant function is a variable, and the other two are parameters, respectively. For each of a plurality of different combinations of two parameters, The above variables are changed to obtain a set of fourth-order cumulants, and for each of the sets of fourth-order cumulants, each fourth-order cumulant is given a function shift of an amount obtained based on the two parameters of each set, and The average is obtained by performing addition processing for each variable.
【0027】請求項5の発明に係わる信号対雑音電力比
改善装置は、入力信号のガウス性雑音成分を抑圧する信
号対雑音電力比改善装置において、入力信号をサンプリ
ングして一組の入力信号を生成する入力信号取り込み手
段と、上記入力信号取り込み手段からの同一の一組の入
力信号について、4次キュムラント関数の3変数の中の
1つを変数、他の2つをそれぞれパラメータとし、複数
個の異なる2つのパラメータの組み合わせのそれぞれに
対し、上記変数を変化させて4次キュムラントの組を求
める4次キュムラント計算処理手段と、上記4次キュム
ラント計算処理手段からの4次キュムラントの組のそれ
ぞれについて、各4次キュムラントに各組の上記2つの
パラメータと上記入力信号中の所望信号の周波数情報に
基づいて位相補正を行う位相補償処理手段と、上記位相
処理手段からの出力を同じ変数ごとに加算処理してその
平均を求める加算処理手段とを備えたものである。According to a fifth aspect of the present invention, there is provided a signal to noise power ratio improving apparatus for suppressing a Gaussian noise component of an input signal. For the input signal capturing means to be generated and the same set of input signals from the input signal capturing means, one of the three variables of the fourth-order cumulant function is used as a variable and the other two are used as parameters, respectively. For each of the combinations of the two different parameters, the above variables are changed to obtain a fourth order cumulant set, and the fourth order cumulant set from the fourth order cumulant calculation processing means Phase correction for each fourth-order cumulant based on the two parameters of each set and the frequency information of the desired signal in the input signal A phase compensation processing means for, in which an adding processing means for obtaining an average by adding processing the output from said phase processing means for each same variable.
【0028】請求項6の発明に係わる信号対雑音電力比
改善方法は、入力信号のガウス性雑音成分を抑圧する信
号対雑音電力比改善方法において、入力信号をサンプリ
ングして一組の入力信号を生成し、同一の上記一組の入
力信号について、4次キュムラント関数の3変数の中の
1つを変数、他の2つをそれぞれパラメータとし、複数
個の異なる2つのパラメータの組み合わせのそれぞれに
対し、上記変数を変化させて4次キュムラントの組を求
め、上記4次キュムラントの組のそれぞれについて、各
4次キュムラントに各組の上記2つのパラメータと上記
入力信号中の所望信号の周波数情報に基づいて位相補正
を行い、同じ変数ごとに加算処理してその平均を求める
ものである。According to a sixth aspect of the present invention, there is provided a signal-to-noise power ratio improving method for suppressing a Gaussian noise component of an input signal. For the same set of input signals to be generated, one of the three variables of the fourth-order cumulant function is a variable, and the other two are parameters, respectively. For each of a plurality of different combinations of two parameters, By changing the variables, a set of fourth-order cumulants is obtained, and for each set of the fourth-order cumulants, each fourth-order cumulant is based on the two parameters of each set and frequency information of a desired signal in the input signal. In this case, phase correction is performed, and addition processing is performed for each of the same variables to obtain an average.
【0029】請求項7の発明に係わる信号対雑音電力比
改善装置は、請求項1、3、または5記載の信号対雑音
電力比改善装置において、加算処理手段に代えて、上記
加算処理手段へのそれぞれの組の入力を周波数成分に変
換して出力するフーリエ変換処理手段と、上記フーリエ
変換処理手段からの出力を同じ帯域ごとに加算処理して
その平均を求める周波数領域加算処理手段とを備えたも
のである。According to a seventh aspect of the present invention, there is provided the signal-to-noise power ratio improving apparatus according to the first, third or fifth aspect, wherein the addition processing means is replaced with the addition processing means. Fourier transform processing means for converting each set of inputs into a frequency component and outputting the same, and frequency domain addition processing means for adding the output from the Fourier transform processing means for each same band and obtaining the average thereof. It is a thing.
【0030】請求項8の発明に係わる信号対雑音電力比
改善方法は、請求項2、4、または6記載の信号対雑音
電力比改善方法において、同じ変数ごとに加算処理して
その平均を求める手順に代えて、フーリエ変換処理によ
り周波数成分に変換し、同じ帯域ごとに加算処理してそ
の平均を求めるものである。According to an eighth aspect of the present invention, there is provided the signal-to-noise power ratio improving method according to the second, fourth or sixth aspect of the present invention, wherein an addition process is performed for each same variable to obtain an average thereof. Instead of the procedure, the frequency component is converted into a frequency component by Fourier transform processing, the addition processing is performed for each same band, and the average is obtained.
【0031】請求項9の発明に係わる信号対雑音電力比
改善装置は、入力信号のガウス性雑音成分を抑圧する信
号対雑音電力比改善装置において、入力信号をサンプリ
ングして一組の入力信号を生成する入力信号取り込み手
段と、上記入力信号取り込み手段からの同一の一組の入
力信号について、4次キュムラント関数の3変数の中の
1つを変数、他の2つをそれぞれパラメータとし、複数
個の異なる2つのパラメータの組み合わせのそれぞれに
対し、上記変数を変化させて4次キュムラントの組を求
める4次キュムラント計算処理手段と、上記4次キュム
ラント計算処理手段からの4次キュムラントの組のそれ
ぞれについて、各4次キュムラントを周波数成分に変換
して出力するフーリエ変換処理手段と、上記フーリエ変
換処理手段からの出力を包絡線検波処理する包絡線検波
処理手段と、上記包絡線検波処理手段からの出力を同じ
帯域ごとに加算処理してその平均を求める加算処理手段
とを備えたものである。According to a ninth aspect of the present invention, there is provided a signal to noise power ratio improving apparatus for suppressing a Gaussian noise component of an input signal. For the input signal capturing means to be generated and the same set of input signals from the input signal capturing means, one of the three variables of the fourth-order cumulant function is used as a variable and the other two are used as parameters, respectively. For each of the combinations of the two different parameters, the above variables are changed to obtain a fourth order cumulant set, and the fourth order cumulant set from the fourth order cumulant calculation processing means A Fourier transform processing means for converting each fourth-order cumulant into a frequency component and outputting the frequency component; And envelope detection means for envelope detection processing power, in which an adding processing means for obtaining an average by adding processing the output from the envelope detection processing means for each same band.
【0032】請求項10の発明に係わる信号対雑音電力
比改善方法は、入力信号のガウス性雑音成分を抑圧する
信号対雑音電力比改善方法において、入力信号をサンプ
リングして一組の入力信号を生成し、同一の上記一組の
入力信号について、4次キュムラント関数の3変数の中
の1つを変数、他の2つをそれぞれパラメータとし、複
数個の異なる2つのパラメータの組み合わせのそれぞれ
に対し、上記変数を変化させて4次キュムラントの組を
求め、上記4次キュムラントの組のそれぞれについて、
各4次キュムラントをフーリエ変換処理により周波数成
分に変換し、包絡線検波処理を行い、同じ帯域ごとに加
算処理してその平均を求めるものである。According to a tenth aspect of the present invention, there is provided a signal-to-noise power ratio improving method for suppressing a Gaussian noise component of an input signal. For the same set of input signals to be generated, one of the three variables of the fourth-order cumulant function is a variable, and the other two are parameters, respectively. For each of a plurality of different combinations of two parameters, , By changing the above variables, to obtain a set of fourth-order cumulants, and for each of the set of fourth-order cumulants,
Each fourth-order cumulant is converted into a frequency component by a Fourier transform process, an envelope detection process is performed, an addition process is performed for each same band, and an average is obtained.
【0033】[0033]
実施の形態1.図1は、本発明の実施の形態1を示す信
号対雑音電力比改善装置の構成図であり、シフトレジス
タ1の出力をM個の4次キュムラント計算回路2に入力
し、その出力を加算回路7で加算するものである。加算
回路7以降の構成は図10の従来例と同等である。ま
た、信号対雑音電力比計算回路6は改善結果の評価用に
接続させてあるにすぎないものである。なお、ここでは
信号対雑音電力比改善方法の実施の形態についても併せ
て記す。Embodiment 1 FIG. FIG. 1 is a block diagram of a signal-to-noise power ratio improving apparatus according to a first embodiment of the present invention. The output of a shift register 1 is input to M fourth-order cumulant calculation circuits 2, and the output is added to an addition circuit. 7 is added. The configuration after the addition circuit 7 is the same as that of the conventional example of FIG. Further, the signal-to-noise power ratio calculation circuit 6 is merely connected for evaluating the improvement result. Here, an embodiment of the signal-to-noise power ratio improving method is also described.
【0034】上記の信号対雑音電力比改善装置の動作を
図1に基づいて説明する。ここで、入力信号x(n)は、従
来例と同様、式(13)に示したように、振幅A、角周波数
ω、初期位相φの複素正弦波信号に、その複素正弦波信
号と独立で、平均ゼロ、分散σ2のガウス性雑音g(n)が
加わったものとしている。図1において、シフトレジス
タ1までは、従来と同様である。すなわち、ある時間間
隔で、サンプリングされた入力信号x(n)は、シフトレジ
スタ1に順次入力される。シフトレジスタ1では、Nサ
ンプルの入力信号を保持する。図10の従来例では、1
つの4次キュムラント計算回路2に対して保持したNサ
ンプルの入力信号を出力していたのに対し、本実施例で
は、M個の4次キュムラント計算回路2に対して、Nサン
プルの同じ入力信号を出力する。各4次キュムラント計
算回路2では、4次キュムラント関数の3変数t1、t2、
t3の中の一つを変数tとし、他の2変数をパラメータと
して、tを最大-(N-1)から(N-1)まで変化させて、M組の
最大2N-1個の4次キュムラントを求める。すなわち、例
えばt1を変数tとした場合、他の2つの変数t2、t3は、
各4次キュムラント計算回路2ごとに、異なる組み合わ
せの定数として、4次キュムラントを求める。その際、
i(1≦i≦M)番目の4次キュムラント計算回路2でのt2
の値をai、t3の値をbiとすると、ai、biは各4次キュム
ラント計算回路2ごとに、Cを共通の任意の定数とし
て、式(15)で示す関係を満たすように設定する。なお、
式(15)の関係は、t1を変数tとし、4キュムラント関数
を式(1)に示したように、x(n+t2)とx(n+t3)を複素共役
として求めた場合に成り立つ関係であり、t2あるいはt3
を変数tとした場合や式(1)の他の位置を複素共役として
4次キュムラント関数を求めた場合は他の関係となる。The operation of the above signal-to-noise power ratio improving apparatus will be described with reference to FIG. Here, similarly to the conventional example, the input signal x (n) is converted into a complex sine wave signal having an amplitude A, an angular frequency ω, and an initial phase φ as shown in Expression (13), independently of the complex sine wave signal. It is assumed that Gaussian noise g (n) with zero mean and variance σ 2 is added. In FIG. 1, up to the shift register 1 is the same as the conventional one. That is, the input signal x (n) sampled at certain time intervals is sequentially input to the shift register 1. The shift register 1 holds an input signal of N samples. In the conventional example of FIG.
While the input signals of N samples held for the four fourth-order cumulant calculation circuits 2 are output, in the present embodiment, the same input signals of N samples are supplied to the M fourth-order cumulant calculation circuits 2. Is output. In each fourth-order cumulant calculation circuit 2, three variables t1, t2,
Using one of t3 as a variable t and the other two variables as parameters, t is changed from-(N-1) to (N-1), and M sets of up to 2N-1 quartic Seek cumulant. That is, for example, if t1 is a variable t, the other two variables t2 and t3 are
For each quaternary cumulant calculation circuit 2, a quaternary cumulant is obtained as a constant of a different combination. that time,
t2 in the i-th (1 ≦ i ≦ M) fourth order cumulant calculation circuit 2
Let ai be the value of t3 and bi be the value of t3. For each of the fourth cumulant calculation circuits 2, ai and bi are set so that C is a common arbitrary constant and satisfies the relationship shown in equation (15). In addition,
The relationship of equation (15) holds when x (n + t2) and x (n + t3) are obtained as complex conjugates, as shown in equation (1), where t1 is a variable t and 4 cumulant functions are shown in equation (1). Relationship, t2 or t3
Is a variable t, or when a fourth-order cumulant function is obtained by setting another position of the equation (1) as a complex conjugate, another relationship is obtained.
【0035】[0035]
【数8】 (Equation 8)
【0036】例えば、C=4とした場合は、ai、biは[a1=
1、b1=3]、[a2=0、b2=4]、[a3=2、b3=2]、[a4=-
1、b4=5]、……のように設定する。ただし、[a1=1、b
1=3]を用いた場合は、その後にa1とb1を入れ替えた
値、すなわち、[ai=3、bi=1]は用いない。For example, when C = 4, ai and bi are [a1 =
1, b1 = 3], [a2 = 0, b2 = 4], [a3 = 2, b3 = 2], [a4 =-
1, b4 = 5], ……. However, [a1 = 1, b
When [1 = 3] is used, a value obtained by exchanging a1 and b1 thereafter, that is, [ai = 3, bi = 1] is not used.
【0037】ai、biをこのように設定することによっ
て、複素正弦波信号のみを入力信号とした場合の4次キ
ュムラント関数(入力信号と同じ周波数の複素正弦波)
の初期位相を、式(16)に示すように、aiとbiの値によら
ず一定にすることができる。By setting ai and bi as described above, a fourth-order cumulant function when only a complex sine wave signal is used as an input signal (a complex sine wave having the same frequency as the input signal)
Can be made constant irrespective of the values of ai and bi, as shown in equation (16).
【0038】[0038]
【数9】 (Equation 9)
【0039】すなわち、複素正弦波信号のみを入力信号
とした場合の、M個の4次キュムラント計算回路2で求
まる、同じtに対する4次キュムラントの値をすべて同
じにすることができる。That is, when only the complex sine wave signal is used as the input signal, the values of the fourth order cumulants for the same t, which are obtained by the M fourth order cumulant calculation circuits 2, can be made the same.
【0040】一方、平均ゼロ、分散σ2のガウス性雑音g
(n)のみを入力信号とした場合、M個の4次キュムラント
計算回路2で求まるM組の4次キュムラント(先に述べ
たように、理論的には、ガウス性雑音に対する4次キュ
ムラントはゼロとなるが、入力信号のサンプル数が有限
個であるために消え残り信号が生じる。)は、aiとbiを
各4次キュムラント計算回路2ごとに、異なる組み合わ
せの定数とすると、同じtに対する4次キュムラントの
値は、各4次キュムラント計算回路2ごとに異なり、し
かも、それらの値の和の平均はゼロに近づく。On the other hand, Gaussian noise g with zero mean and variance σ 2
When only (n) is an input signal, M sets of fourth-order cumulants determined by the M fourth-order cumulant calculation circuits 2 (theoretically, the fourth-order cumulant for Gaussian noise is zero, as described above) However, the remaining signal is generated because the number of samples of the input signal is finite.) If ai and bi are constants of different combinations for each fourth-order cumulant calculation circuit 2, 4 for the same t The value of the next cumulant differs for each fourth-order cumulant calculation circuit 2, and the average of the sum of those values approaches zero.
【0041】このようにして求めたM組の4次キュムラ
ントは加算回路7に出力される。ただし、式(11)にも示
したように、4次キュムラントの振幅は、入力信号の振
幅の4乗となるために、加算回路7へは、振幅の4乗根
をとった4次キュムラントを出力する。加算回路7で
は、同じtごとにM組の4次キュムラントを足し、tごと
の平均をとる。このような処理を行うことによって、複
素正弦波信号成分に対する4次キュムラントの値を変化
させず、ガウス性雑音成分に対する4次キュムラントだ
けをゼロに近づけることができる。すなわち、信号対雑
音電力比が向上することができる。The M sets of fourth-order cumulants thus obtained are output to the adder circuit 7. However, since the amplitude of the fourth-order cumulant is the fourth power of the amplitude of the input signal as shown in the equation (11), the fourth-order cumulant having the fourth root of the amplitude is supplied to the adding circuit 7. Output. The adder circuit 7 adds M sets of fourth-order cumulants for each t, and takes an average for each t. By performing such processing, it is possible to make only the fourth-order cumulant for the Gaussian noise component close to zero without changing the value of the fourth-order cumulant for the complex sine wave signal component. That is, the signal-to-noise power ratio can be improved.
【0042】加算回路7の出力信号は、窓関数回路3に
入力される。窓関数回路3では、高速フーリエ変換回路
4における、加算回路7の出力信号の数が有限であるこ
との影響を低減するために、窓関数による重みを掛け、
その結果を、高速フーリエ変換回路4に出力する。高速
フーリエ変換回路4では、窓関数回路3の出力信号を、
加算回路7の出力信号の数と入力信号x(n)のサンプリン
グ時間で決まる帯域ごとの成分に分別し、自乗検波回路
5に出力する。以降は従来例と同様に、自乗検波回路5
において、各帯域の信号の電力を求め、信号対雑音電力
比計算回路6に出力する。信号対雑音電力比計算回路6
では、入力信号x(n)の複素正弦波信号の周波数が分かっ
ている、あるいは予測している場合は、周波数の存在す
る帯域の電力値から、周波数の存在する帯域以外の電力
の平均値を引いた値と、周波数の存在する帯域以外の電
力の平均値との比を求める。また、入力信号x(n)の複素
正弦波信号の周波数が分からない場合は、最も大きい電
力値からその他の帯域の電力の平均値を引いた値と、最
も大きい電力値を持つ帯域以外の帯域の電力の平均値と
の比を求める。The output signal of the adding circuit 7 is input to the window function circuit 3. In the window function circuit 3, in order to reduce the influence of the finite number of output signals of the adding circuit 7 in the fast Fourier transform circuit 4, weights are multiplied by a window function,
The result is output to the fast Fourier transform circuit 4. In the fast Fourier transform circuit 4, the output signal of the window function circuit 3 is
The signal is separated into components for each band determined by the number of output signals of the adder circuit 7 and the sampling time of the input signal x (n), and output to the square detection circuit 5. Thereafter, as in the conventional example, the square detection circuit 5
, The power of the signal in each band is obtained and output to the signal-to-noise power ratio calculation circuit 6. Signal to noise power ratio calculation circuit 6
Then, if the frequency of the complex sine wave signal of the input signal x (n) is known or predicted, the average value of the power outside the frequency band is calculated from the power value of the frequency band. The ratio between the subtracted value and the average value of the power other than the frequency band is calculated. If the frequency of the complex sine wave signal of the input signal x (n) is not known, a value obtained by subtracting the average value of the power of the other bands from the largest power value and a band other than the band having the largest power value And the ratio of the average power to the average value.
【0043】図2に、図10の従来の信号対雑音電力比
改善装置での自乗検波回路5の出力信号例を示す。ただ
し、入力信号の信号対雑音電力比SNRinは-5[dB]、入力
信号のサンプル数Nは32、入力する複素正弦波信号の規
格化周波数は9/32としている。また、図3に上記図2と
同じ入力信号を図1に示した実施の形態対1の信号対雑
音電力比改善装置に入力した場合の自乗検波回路5の出
力信号例を示す。ただし、1組の入力信号から10組の
4次キュムラントを求めている。なお、ここでai、bi
は、C=0として、[a1=0、b1=0]、[a2=-1、b2=1]、
[a3=-2、b3=2]、[a4=-3、b4=3]、……、[a10=-9、
b10=9]としている。FIG. 2 shows an example of an output signal of the square-law detection circuit 5 in the conventional signal-to-noise-power ratio improving apparatus of FIG. However, the signal-to-noise power ratio SNRin of the input signal is -5 [dB], the number N of samples of the input signal is 32, and the normalized frequency of the input complex sine wave signal is 9/32. FIG. 3 shows an example of the output signal of the square detection circuit 5 when the same input signal as in FIG. 2 is input to the signal-to-noise power ratio improving apparatus of the embodiment shown in FIG. However, ten sets of fourth-order cumulants are obtained from one set of input signals. Where ai, bi
Is [a1 = 0, b1 = 0], [a2 = -1, b2 = 1]
[A3 = -2, b3 = 2], [a4 = -3, b4 = 3],…, [a10 = -9,
b10 = 9].
【0044】図2の従来の信号対雑音電力比計算回路6
での信号対雑音電力比は、13.97[dB]であるのに対し、
図3の実施の形態1の信号対雑音電力比計算回路6での
信号対雑音電力比は、16.40[dB]となり、約2.5[dB]信号
対雑音電力比改善効果が向上していることが分かる。The conventional signal-to-noise power ratio calculation circuit 6 of FIG.
Signal-to-noise power ratio is 13.97 [dB],
The signal-to-noise power ratio in the signal-to-noise power ratio calculation circuit 6 of the first embodiment shown in FIG. 3 is 16.40 [dB], which means that the signal-to-noise power ratio improvement effect is improved by about 2.5 [dB]. I understand.
【0045】実施の形態2.図4は、本発明の実施の形
態2を示す信号対雑音電力比改善装置の構成図であり、
上記実施の形態1と異なる部分は、M個の4次キュムラ
ント計算回路2の後に、窓関数回路と、高速フーリエ変
換回路を4次キュムラント計算回路と同じように、M個
設け、加算回路7の後の窓関数回路3と、高速フーリエ
変換回路4をなくしたことである。なお、ここでは信号
対雑音電力比改善方法の実施の形態についても併せて記
す。Embodiment 2 FIG. 4 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 2 of the present invention.
The difference from the first embodiment is that a window function circuit and M fast Fourier transform circuits are provided after the M fourth-order cumulant calculation circuits 2 in the same manner as the fourth-order cumulant calculation circuit. The latter is that the window function circuit 3 and the fast Fourier transform circuit 4 are eliminated. Here, an embodiment of the signal-to-noise power ratio improving method is also described.
【0046】上記の信号対雑音電力比改善装置の動作を
図4に基づいて説明する。ここで、入力信号x(n)は、実
施の形態1と同様、式(13)に示したように、振幅A、角
周波数ω、初期位相φの複素正弦波信号に、その複素正
弦波信号と独立で、平均ゼロ、分散σ2のガウス性雑音g
(n)が加わったものとしている。図4において、4次キ
ュムラント計算回路2までは、上記実施の形態1と同様
である。すなわち、ある時間間隔で、サンプリングされ
た入力信号x(n)は、シフトレジスタ1に順次入力され
る。シフトレジスタ1では、Nサンプルの入力信号を保
持し、M個の4次キュムラント計算回路2に対して、Nサ
ンプルの同じ入力信号を出力する。各4次キュムラント
計算回路2では、4次キュムラント関数の3変数t1、t
2、t3の中の一つを変数tとし、他の2変数をパラメータ
として、tを最大-(N-1)から(N-1)まで変化させて、M組
の最大2N-1個の4次キュムラントを求める。すなわち、
例えばt1を変数tとした場合、他の2つの変数t2、t3
は、各4次キュムラント計算回路2ごとに、異なる組み
合わせの定数として、4次キュムラントを求める。その
際、i(1≦i≦M)番目の4次キュムラント計算回路2で
のt2の値をai、t3の値をbiとすると、ai、biは各4次キ
ュムラント計算回路2ごとに、Cを任意の定数として、
式(15)で示す関係を満たすように設定する。The operation of the above signal-to-noise power ratio improving apparatus will be described with reference to FIG. Here, similarly to the first embodiment, the input signal x (n) is converted into a complex sine wave signal having an amplitude A, an angular frequency ω, and an initial phase φ as shown in Expression (13). Gaussian noise g with zero mean and variance σ 2
(n) has been added. In FIG. 4, the configuration up to the fourth cumulant calculation circuit 2 is the same as that of the first embodiment. That is, the input signal x (n) sampled at certain time intervals is sequentially input to the shift register 1. The shift register 1 holds the input signals of N samples, and outputs the same input signals of N samples to the M fourth-order cumulant calculation circuits 2. In each fourth-order cumulant calculation circuit 2, three variables t1, t of the fourth-order cumulant function
2, one of t3 is a variable t, and the other two variables are parameters, and t is changed from-(N-1) to (N-1). Find the fourth cumulant. That is,
For example, if t1 is a variable t, the other two variables t2 and t3
Calculates a fourth-order cumulant as a constant of a different combination for each fourth-order cumulant calculation circuit 2. At this time, if the value of t2 in the i-th (1 ≦ i ≦ M) fourth-order cumulant calculation circuit 2 is ai and the value of t3 is bi, ai and bi are C Is an arbitrary constant,
The setting is made so as to satisfy the relationship shown in Expression (15).
【0047】実施の形態1では、このようにして求めた
M組の4次キュムラントを同じtに対して足していたのに
対し、本実施の形態2では、各4次キュムラントを周波
数成分に変換し、同じ帯域ごとに足す。すなわち、M個
の4次キュムラント計算回路2で求めたM組の4次キュ
ムラントは、それぞれに対応する窓関数回路3に入力さ
れる。ただし、式(11)にも示したように、4次キュムラ
ントの振幅は、入力信号の振幅の4乗となるために、各
窓関数回路3へは、振幅の4乗根をとった4次キュムラ
ントを出力する。各窓関数回路3では、各4次キュムラ
ント計算回路2の出力信号に対して、高速フーリエ変換
回路4における、4次キュムラントの数が有限であるこ
との影響を低減するために、窓関数の重みを掛け、その
結果を、それぞれに対応する高速フーリエ変換回路4に
出力する。各高速フーリエ変換回路4では、各窓関数回
路3の出力信号を、各4次キュムラントの数と入力信号
x(n)のサンプリング時間で決まる帯域ごとの成分に分別
し、加算回路7に出力する。加算回路7では、各帯域ご
とにM個の信号を加算し、各帯域ごとに平均値を求め
る。In the first embodiment, the value is obtained in this manner.
While M sets of fourth-order cumulants are added for the same t, in the second embodiment, each fourth-order cumulant is converted into a frequency component and added for each same band. That is, the M sets of quaternary cumulants calculated by the M quaternary cumulant calculation circuits 2 are input to the corresponding window function circuits 3. However, as shown in equation (11), the amplitude of the fourth-order cumulant is the fourth power of the amplitude of the input signal. Output cumulant. In each window function circuit 3, the weight of the window function is reduced in order to reduce the effect of the finite number of fourth-order cumulants in the fast Fourier transform circuit 4 on the output signal of each fourth-order cumulant calculation circuit 2. And outputs the result to the corresponding fast Fourier transform circuit 4. Each fast Fourier transform circuit 4 converts the output signal of each window function circuit 3 into the number of each fourth-order cumulant and the input signal.
It is separated into components for each band determined by the sampling time of x (n), and is output to the addition circuit 7. The adder circuit 7 adds M signals for each band and calculates an average value for each band.
【0048】実施の形態1でも示したように、ai、biの
値を式(15)で示す関係を満たすように設定することによ
って、入力信号を複素正弦波信号のみとした時の、M個
の4次キュムラント計算回路2で求まる、同じtに対す
る4次キュムラントの値をすべて同じにすることができ
る。よって、M個の4次キュムラント計算回路2で求め
た、M組の4次キュムラントをフーリエ変換した場合、
複素正弦波信号の周波数成分を含む帯域の値はすべて同
じになる。一方、平均ゼロ、分散σ2のガウス性雑音g
(n)のみを入力信号とした場合、aiとbiをM個の4次キュ
ムラント計算回路2ごとに、異なる組み合わせの定数と
して、M組の4次キュムラントを求めると、同じtに対す
る4次キュムラントの値は、各4次キュムラントごとに
異なり、しかも、それらの和の平均はゼロに近づく。そ
のため、それらM組の4次キュムラントをフーリエ変換
した時の各帯域の値も、各4次キュムラントごとに異な
り、各帯域ごとの値の和の平均もゼロに近づく。よっ
て、入力信号を、複素正弦波信号に、その複素正弦波信
号と独立で、平均ゼロ、分散σ2のガウス性雑音g(n)が
加わった信号として、M個の4次キュムラント計算回路
2で、式(15)のような条件のもとに求めたM組の4次キ
ュムラントを、それぞれフーリエ変換し、各帯域ごとに
足して平均すると、複素正弦波信号成分に対する信号を
変化させず、ガウス性雑音成分に対する信号だけをゼロ
に近づけることができ、結果として、信号対雑音電力比
を向上することができる。As described in the first embodiment, by setting the values of ai and bi so as to satisfy the relationship represented by the equation (15), the number of M signals when the input signal is only a complex sine wave signal is obtained. The values of the fourth-order cumulants for the same t obtained by the fourth-order cumulant calculation circuit 2 can be made the same. Therefore, when Fourier transform is performed on M sets of fourth-order cumulants obtained by the M fourth-order cumulant calculation circuits 2,
The values of the bands including the frequency components of the complex sine wave signal are all the same. On the other hand, Gaussian noise g with mean zero and variance σ 2
When only (n) is an input signal, ai and bi are set as constants of different combinations for each of the M quaternary cumulant calculation circuits 2 and M sets of quaternary cumulants are obtained. The values differ for each fourth-order cumulant, and the average of their sum approaches zero. Therefore, the value of each band when the M sets of fourth-order cumulants are subjected to Fourier transform also differs for each fourth-order cumulant, and the average of the sum of the values for each band approaches zero. Therefore, the input signal is a signal obtained by adding a Gaussian noise g (n) having a mean of zero and a variance of σ 2 to the complex sine wave signal independently of the complex sine wave signal. Then, the M sets of fourth-order cumulants obtained under the conditions as in equation (15) are Fourier-transformed, added and averaged for each band, and the signal for the complex sine wave signal component does not change. Only the signal for the Gaussian noise component can be made close to zero, and as a result, the signal-to-noise power ratio can be improved.
【0049】加算回路7の出力信号は、自乗検波回路5
に入力される。自乗検波回路5では、各帯域の信号の電
力を求め、信号対雑音電力比計算回路6に出力される。
信号対雑音電力比計算回路6では、入力信号x(n)の複素
正弦波信号の周波数が分かっている、あるいは予測して
いる場合は、周波数の存在する帯域の電力値から、周波
数の存在する帯域以外の電力の平均値を引いた値と、周
波数の存在する帯域以外の電力の平均値との比を求め
る。また、入力信号x(n)の複素正弦波信号の周波数が分
からない場合は、最も大きい電力値からその他の帯域の
電力の平均値を引いた値と、最も大きい電力値を持つ帯
域以外の帯域の電力の平均値との比を求める。The output signal of the adding circuit 7 is
Is input to The square detection circuit 5 obtains the power of the signal in each band and outputs it to the signal-to-noise power ratio calculation circuit 6.
In the signal-to-noise power ratio calculation circuit 6, when the frequency of the complex sine wave signal of the input signal x (n) is known or predicted, the frequency exists based on the power value of the frequency band. The ratio between the value obtained by subtracting the average value of the power outside the band and the average value of the power outside the band where the frequency exists is determined. If the frequency of the complex sine wave signal of the input signal x (n) is not known, a value obtained by subtracting the average value of the power of the other bands from the largest power value and a band other than the band having the largest power value And the ratio of the average power to the average value.
【0050】実施の形態3.図5は、本発明の実施の形
態3を示す信号対雑音電力比改善装置の構成図であり、
上記実施の形態1と異なる部分は、M個の4次キュムラ
ント計算回路の後に、M個の関数シフト回路8を設けた
ことである。なお、ここでは信号対雑音電力比改善方法
の実施の形態についても併せて記す。Embodiment 3 FIG. 5 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 3 of the present invention.
The difference from the first embodiment is that M function shift circuits 8 are provided after M fourth-order cumulant calculation circuits. Here, an embodiment of the signal-to-noise power ratio improving method is also described.
【0051】上記の信号対雑音電力比改善装置の動作を
図5に基づいて説明する。ここで、入力信号x(n)は、実
施の形態1と同様、式(13)に示したように、振幅A、角
周波数ω、初期位相φの複素正弦波信号に、その複素正
弦波信号と独立で、平均ゼロ、分散σ2のガウス性雑音g
(n)が加わったものとしている。図5において、M個の4
次キュムラント計算回路2の入力までは上記実施の形態
1と同じである。すなわち、ある時間間隔で、サンプリ
ングされた入力信号x(n)は、シフトレジスタ1に順次入
力される。シフトレジスタ1では、Nサンプルの入力信
号を保持し、M個の4次キュムラント計算回路2に対し
て、Nサンプルの同じ入力信号を出力する。各4次キュ
ムラント計算回路2では、4次キュムラント関数の3変
数t1、t2、t3の中の一つを変数tとし、他の2変数をパ
ラメータとして、tを最大-(N-1)から(N-1)まで変化させ
て、M組の最大2N-1個の4次キュムラントを求める。す
なわち、例えばt1を変数tとした場合、他の2つの変数t
2、t3は、各4次キュムラント計算回路2ごとに、異な
る組み合わせの定数として、4次キュムラントを求め
る。その際に、実施の形態1では、t2、t3の値ai、bi
(i=1,2,…,M)は、各4次キュムラント計算回路2ごと
に、式(15)で示す関係を満たす異なった組合せの定数と
するという制約があった。それに対し、本実施の形態で
は、ai、biは、各4次キュムラント計算回路2ごとに、
異なった組合せの定数であればよい。ただし、実施の形
態1と同様に、ai、biの値を入れ替えた値は同時には用
いない。The operation of the above signal-to-noise power ratio improving apparatus will be described with reference to FIG. Here, similarly to the first embodiment, the input signal x (n) is converted into a complex sine wave signal having an amplitude A, an angular frequency ω, and an initial phase φ as shown in Expression (13). Gaussian noise g with zero mean and variance σ 2
(n) has been added. In FIG. 5, M 4
The operation up to the input of the next cumulant calculation circuit 2 is the same as that of the first embodiment. That is, the input signal x (n) sampled at certain time intervals is sequentially input to the shift register 1. The shift register 1 holds the input signals of N samples, and outputs the same input signals of N samples to the M fourth-order cumulant calculation circuits 2. In each of the fourth-order cumulant calculation circuits 2, one of the three variables t1, t2, and t3 of the fourth-order cumulant function is set as a variable t, and the other two variables are set as parameters. N-1) to obtain M sets of up to 2N-1 quaternary cumulants. That is, for example, if t1 is a variable t, the other two variables t
2, t3 is to obtain the fourth cumulant as a constant of a different combination for each fourth cumulant calculation circuit 2. At this time, in the first embodiment, the values ai, bi of t2 and t3
(I = 1, 2,..., M) is restricted to a constant of a different combination that satisfies the relationship shown in Expression (15) for each fourth-order cumulant calculation circuit 2. On the other hand, in the present embodiment, ai and bi are calculated for each fourth-order cumulant calculation circuit 2,
What is necessary is just a constant of a different combination. However, as in the first embodiment, the values obtained by replacing the values of ai and bi are not used at the same time.
【0052】このようにして求めたM組の4次キュムラ
ントは、それぞれ対応する関数シフト回路8に入力され
る。ただし、式(11)にも示したように、4次キュムラン
トの振幅は、入力信号の振幅の4乗となるために、各関
数シフト回路8へは、振幅の4乗根をとった4次キュム
ラントを出力する。各関数シフト回路8では、各4次キ
ュムラント計算回路2で求めたtの変化に対する4次キ
ュムラントをt軸上の負の方向に(ai+bi)シフトさせ
る。例えば、[a1=1、b1=3]として4次キュムラントを
求めた場合、t=10に対して求めた4次キュムラントをt=
6に対して求めた4次キュムラントとする。なお、この
4次キュムラントをシフトさせる方向およびシフトさせ
る量は、t1を変数tとし、4キュムラント関数を式(1)に
示したように、x(n+t2)とx(n+t3)を複素共役として求め
た場合に成り立つ関係であり、t2あるいはt3を変数tと
した場合や式(1)の他の位置を複素共役として4次キュ
ムラント関数を求めた場合は異なるシフト方向、異なる
シフト量となる。このような処理を行うことによって、
複素正弦波信号のみを入力信号とした場合の4次キュム
ラント関数(入力信号と同じ周波数の複素正弦波)の初
期位相を、aiとbiの値によらず一定にすることができ
る。すなわち、複素正弦波信号のみを入力信号とした場
合の、M個の4次キュムラント計算回路2で求まった、
同じtに対する4次キュムラントの値をすべて同じにす
ることができる。一方、平均ゼロ、分散σ2のガウス性
雑音g(n)のみを入力信号とした場合は、同じtに対する
4次キュムラントの値は、各4次キュムラント計算回路
2ごとに異なり、しかも、それらの和の平均はゼロに近
づく。The M sets of fourth-order cumulants thus obtained are input to the corresponding function shift circuits 8. However, since the amplitude of the fourth-order cumulant is the fourth power of the amplitude of the input signal, as shown in the equation (11), the fourth-order root of the amplitude is given to each function shift circuit 8. Output cumulant. Each function shift circuit 8 shifts the fourth-order cumulant in the negative direction on the t-axis (ai + bi) with respect to the change of t obtained by each fourth-order cumulant calculation circuit 2. For example, when a fourth order cumulant is obtained as [a1 = 1, b1 = 3], the fourth order cumulant obtained for t = 10 is t = 10.
The fourth cumulant determined for 6. The direction and the amount of shift of the fourth cumulant are defined as x (n + t2) and x (n + t3), where t1 is a variable t and the 4 cumulant function is shown in equation (1). This relationship holds when calculated as a complex conjugate. When t2 or t3 is a variable t, or when a fourth-order cumulant function is calculated with another position of Equation (1) as a complex conjugate, the shift direction and shift amount are different. Becomes By performing such processing,
When only a complex sine wave signal is used as an input signal, the initial phase of a fourth-order cumulant function (a complex sine wave having the same frequency as the input signal) can be constant regardless of the values of ai and bi. In other words, when only the complex sine wave signal is used as the input signal, the number of the M fourth-order cumulant calculation circuits 2 is calculated.
All fourth order cumulant values for the same t can be the same. On the other hand, when only the Gaussian noise g (n) having zero mean and variance σ 2 is used as the input signal, the value of the fourth-order cumulant for the same t is different for each fourth-order cumulant calculation circuit 2. The average of the sums approaches zero.
【0053】このようにして求めたM組の関数シフト回
路8の出力信号は加算回路7に出力される。加算回路7
では、同じtごとにM組の関数シフト回路8の出力信号を
足し、tごとの平均をとる。このような処理を行うこと
によって、複素正弦波信号成分に対する4次キュムラン
トの値を変化させずに、ガウス性雑音成分に対する4次
キュムラントだけをゼロに近づけることができる。すな
わち、信号対雑音電力比を向上させることができる。The output signals of the M sets of function shift circuits 8 thus obtained are output to the adder circuit 7. Adder circuit 7
So the same t your capital adds the output signal of the M sets of function shift circuit 8, taking the average of each t. By performing such a process, it is possible to make only the fourth-order cumulant for the Gaussian noise component close to zero without changing the value of the fourth-order cumulant for the complex sine wave signal component. That is, it is possible to make improved signal-to-noise power ratio.
【0054】その後の処理は、上記実施の形態1と同様
に、加算回路7の出力信号は、窓関数回路3に入力され
る。窓関数回路3では、高速フーリエ変換回路4におけ
る、加算回路7の出力信号の数が有限であることの影響
を低減するために、窓関数による重みを掛け、その結果
を、高速フーリエ変換回路4に出力する。高速フーリエ
変換回路4では、窓関数回路3の出力信号を、加算回路
7の出力信号の数と入力信号x(n)のサンプリング時間で
決まる帯域ごとの成分に分別し、自乗検波回路5に出力
する。自乗検波回路5では、各帯域の信号の電力を求
め、信号対雑音電力比計算回路6に出力される。信号対
雑音電力比計算回路6では、入力信号x(n)の複素正弦波
信号の周波数が分かっている、あるいは予測している場
合は、周波数の存在する帯域の電力値から、周波数の存
在する帯域以外の電力の平均値を引いた値と、周波数の
存在する帯域以外の電力の平均値との比を求める。ま
た、入力信号x(n)の複素正弦波信号の周波数が分からな
い場合は、最も大きい電力値からその他の帯域の電力の
平均値を引いた値と、最も大きい電力値を持つ帯域以外
の帯域の電力の平均値との比を求める。In the subsequent processing, the output signal of the adding circuit 7 is input to the window function circuit 3 as in the first embodiment. The window function circuit 3 multiplies the weight by the window function in order to reduce the influence of the finite number of output signals of the adder circuit 7 in the fast Fourier transform circuit 4, and divides the result by the fast Fourier transform circuit 4. Output to The fast Fourier transform circuit 4 separates the output signal of the window function circuit 3 into components for each band determined by the number of output signals of the adder circuit 7 and the sampling time of the input signal x (n), and outputs it to the square detection circuit 5. I do. The square detection circuit 5 obtains the power of the signal in each band and outputs it to the signal-to-noise power ratio calculation circuit 6. In the signal-to-noise power ratio calculation circuit 6, when the frequency of the complex sine wave signal of the input signal x (n) is known or predicted, the frequency exists based on the power value of the frequency band. The ratio between the value obtained by subtracting the average value of the power outside the band and the average value of the power outside the band where the frequency exists is determined. If the frequency of the complex sine wave signal of the input signal x (n) is not known, a value obtained by subtracting the average value of the power of the other bands from the largest power value and a band other than the band having the largest power value And the ratio of the average power to the average value.
【0055】実施の形態4.図6は、本発明の実施の形
態4を示す信号対雑音電力比改善装置の構成図であり、
図4に基づいて説明した上記実施の形態2と異なる部分
は、M個の4次キュムラント計算回路2の後に、M個の関
数シフト回路8を設けたことである。なお、ここでは信
号対雑音電力比改善方法の実施の形態についても併せて
記す。Embodiment 4 FIG. 6 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 4 of the present invention.
The difference from the second embodiment described with reference to FIG. 4 is that M function shift circuits 8 are provided after M fourth-order cumulant calculation circuits 2. Here, an embodiment of the signal-to-noise power ratio improving method is also described.
【0056】上記の信号対雑音電力比改善装置の動作を
図6に基づいて説明する。ただし、入力信号x(n)は、実
施の形態2と同様、式(13)に示したように、振幅A、角
周波数ω、初期位相φの複素正弦波信号に、その複素正
弦波信号と独立で、平均ゼロ、分散σ2のガウス性雑音g
(n)が加わったものとしている。図6において、M個の4
次キュムラント計算回路2の入力までは上記実施の形態
2と同じである。すなわち、ある時間間隔で、サンプリ
ングされた入力信号x(n)は、シフトレジスタ1に順次入
力される。シフトレジスタ1では、Nサンプルの入力信
号を保持し、M個の4次キュムラント計算回路2に対し
て、Nサンプルの同じ入力信号を出力する。各4次キュ
ムラント計算回路2では、4次キュムラント関数の3変
数t1、t2、t3の中の一つを変数tとし、他の2変数をパ
ラメータとして、tを最大-(N-1)から(N-1)まで変化させ
て、M組の最大2N-1個の4次キュムラントを求める。す
なわち、例えばt1を変数tとした場合、他の2つの変数t
2、t3は、各4次キュムラント計算回路2ごとに、異な
る組み合わせの定数として、4次キュムラントを求め
る。その際に、実施の形態2では、t2、t3の値ai、bi
(i=1,2,…,M)は、各4次キュムラント計算回路2ごと
に、式(15)で示す関係を満たす異なった組合せの定数と
するという制約があった。それに対し、本実施の形態で
は、ai、biは、各4次キュムラント計算回路2ごとに、
異なった組合せの定数であればよい。ただし、実施の形
態2と同様に、ai、biの値を入れ替えた値は同時には用
いない。The operation of the signal to noise power ratio improving apparatus will be described with reference to FIG. However, as in the second embodiment, the input signal x (n) is a complex sine wave signal having an amplitude A, an angular frequency ω, and an initial phase φ as shown in Expression (13), and the complex sine wave signal Gaussian noise g with independent zero mean and variance σ 2
(n) has been added. In FIG. 6, M 4
The operation up to the input of the next cumulant calculation circuit 2 is the same as that of the second embodiment. That is, the input signal x (n) sampled at certain time intervals is sequentially input to the shift register 1. The shift register 1 holds the input signals of N samples, and outputs the same input signals of N samples to the M fourth-order cumulant calculation circuits 2. In each of the fourth-order cumulant calculation circuits 2, one of the three variables t1, t2, and t3 of the fourth-order cumulant function is set as a variable t, and the other two variables are set as parameters. N-1) to obtain M sets of up to 2N-1 quaternary cumulants. That is, for example, if t1 is a variable t, the other two variables t
2, t3 is to obtain the fourth cumulant as a constant of a different combination for each fourth cumulant calculation circuit 2. At this time, in the second embodiment, the values ai, bi of t2 and t3
(I = 1, 2,..., M) is restricted to a constant of a different combination that satisfies the relationship shown in Expression (15) for each fourth-order cumulant calculation circuit 2. On the other hand, in the present embodiment, ai and bi are calculated for each fourth-order cumulant calculation circuit 2,
What is necessary is just a constant of a different combination. However, as in the second embodiment, the values obtained by replacing the values of ai and bi are not used at the same time.
【0057】このようにして求めたM組の4次キュムラ
ンは、それぞれ対応する関数シフト回路8に入力され
る。ただし、式(11)にも示したように、4次キュムラン
トの振幅は、入力信号の振幅の4乗となるために、各関
数シフト回路8へは、振幅の4乗根をとった4次キュム
ラントを出力する。各関数シフト回路8では、各4次キ
ュムラント計算回路2で求めたtの変化に対する4次キ
ュムラントをt軸上の負の方向に(ai+bi)シフトさせ
る。なお、この4次キュムラントをシフトさせる方向お
よびシフトさせる量は、t1を変数tとし、4キュムラン
ト関数を式(1)に示したように、x(n+t2)とx(n+t3)を複
素共役として求めた場合に成り立つ関係であり、t2ある
いはt3を変数tとした場合や式(1)の他の位置を複素共役
として4次キュムラント関数を求めた場合は異なるシフ
ト方向、異なるシフト量となる。このような処理を行う
ことによって、複素正弦波信号のみを入力信号とした場
合の4次キュムラント関数(入力信号と同じ周波数の複
素正弦波)の初期位相を、aiとbiの値によらず一定にす
ることができる。すなわち、複素正弦波信号のみを入力
信号とした場合、M個の4次キュムラント計算回路2で
求まる、同じtに対する4次キュムラントの値をすべて
同じにすることができる。一方、平均ゼロ、分散σ2の
ガウス性雑音g(n)のみを入力信号とした場合は、同じt
に対する4次キュムラントの値は、各4次キュムラント
計算回路2ごとに異なり、しかも、それらの値の平均は
ゼロに近づく。The M sets of fourth-order cumulans determined as described above are input to the corresponding function shift circuits 8. However, since the amplitude of the fourth-order cumulant is the fourth power of the amplitude of the input signal, as shown in the equation (11), the fourth-order root of the amplitude is given to each function shift circuit 8. Output cumulant. Each function shift circuit 8 shifts the fourth-order cumulant in the negative direction on the t-axis (ai + bi) with respect to the change of t obtained by each fourth-order cumulant calculation circuit 2. The direction and the amount of shift of the fourth cumulant are defined as x (n + t2) and x (n + t3), where t1 is a variable t and the 4 cumulant function is shown in equation (1). This relationship holds when calculated as a complex conjugate. When t2 or t3 is a variable t, or when a fourth-order cumulant function is calculated with another position of Equation (1) as a complex conjugate, the shift direction and shift amount are different. Becomes By performing such processing, the initial phase of the fourth-order cumulant function (a complex sine wave having the same frequency as the input signal) when only the complex sine wave signal is used as the input signal is constant regardless of the values of ai and bi. Can be That is, when only the complex sine wave signal is used as the input signal, the values of the fourth order cumulants for the same t obtained by the M fourth order cumulant calculation circuits 2 can be all the same. On the other hand, when only the Gaussian noise g (n) having zero mean and variance σ 2 is used as the input signal, the same t
The value of the fourth-order cumulant differs for each fourth-order cumulant calculation circuit 2, and the average of those values approaches zero.
【0058】このようにして求めたM組の関数シフト回
路8の出力信号はそれぞれ対応する窓関数回路3に入力
される。各窓関数回路3では、各関数シフト回路8の出
力信号に対して、高速フーリエ変換回路4における、4
次キュムラントの数が有限であることの影響を低減する
ために、窓関数の重みを掛け、その結果を、それぞれに
対応する高速フーリエ変換回路4に出力する。その後の
処理は、上記実施の形態2と同様に、各高速フーリエ変
換回路4において、それぞれの窓関数回路3の出力信号
を、関数シフト回路8の出力信号の数と入力信号x(n)の
サンプリング時間で決まる帯域ごとの成分に分別し、加
算回路7に出力する。加算回路7では、各帯域ごとにM
個の信号を加算し、各帯域ごとに平均値を求める。この
ような処理を行うことによって、複素正弦波信号成分に
対する信号を変化させず、ガウス性雑音成分に対する信
号だけをゼロに近づけることができ、結果として、信号
対雑音電力比を向上することができる。加算回路7の出
力信号は、自乗検波回路5に入力される。自乗検波回路
5では、各帯域の信号の電力を求め、信号対雑音電力比
計算回路6に出力される。信号対雑音電力比計算回路6
では、入力信号x(n)の複素正弦波信号の周波数が分かっ
ている、あるいは予測している場合は、周波数の存在す
る帯域の電力値から、周波数の存在する帯域以外の電力
の平均値を引いた値と、周波数の存在する帯域以外の電
力の平均値との比を求める。また、入力信号x(n)の複素
正弦波信号の周波数が分からない場合は、最も大きい電
力値からその他の帯域の電力の平均値を引いた値と、最
も大きい電力値を持つ帯域以外の帯域の電力の平均値と
の比を求める。The output signals of the M sets of function shift circuits 8 thus obtained are input to the corresponding window function circuits 3. In each window function circuit 3, the output signal of each function shift circuit 8 is applied to the output signal of the fast Fourier transform circuit 4.
In order to reduce the influence of the finite number of next cumulants, the weight of the window function is multiplied, and the result is output to the corresponding fast Fourier transform circuit 4. In the subsequent processing, as in the second embodiment, in each fast Fourier transform circuit 4, the output signal of each window function circuit 3 is converted to the number of output signals of the function shift circuit 8 and the input signal x (n). The signal is separated into components for each band determined by the sampling time, and output to the adding circuit 7. In the addition circuit 7, M
The signals are added, and an average value is obtained for each band. By performing such processing, the signal for the Gaussian noise component can be made close to zero without changing the signal for the complex sine wave signal component, and as a result, the signal-to-noise power ratio can be improved. . The output signal of the addition circuit 7 is input to the square detection circuit 5. The square detection circuit 5 obtains the power of the signal in each band and outputs it to the signal-to-noise power ratio calculation circuit 6. Signal to noise power ratio calculation circuit 6
Then, if the frequency of the complex sine wave signal of the input signal x (n) is known or predicted, the average value of the power outside the frequency band is calculated from the power value of the frequency band. The ratio between the subtracted value and the average value of the power other than the frequency band is calculated. If the frequency of the complex sine wave signal of the input signal x (n) is not known, a value obtained by subtracting the average value of the power of the other bands from the largest power value and a band other than the band having the largest power value And the ratio of the average power to the average value.
【0059】実施の形態5.図7は、本発明の実施の形
態5を示す信号対雑音電力比改善装置の構成図であり、
上記実施の形態1と異なる部分は、M個の4次キュムラ
ント計算回路の後に、M個の位相補正回路9を設けたこ
とである。なお、ここでは信号対雑音電力比改善方法の
実施の形態についても併せて記す。Embodiment 5 FIG. FIG. 7 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 5 of the present invention.
The difference from the first embodiment is that M phase correction circuits 9 are provided after M fourth-order cumulant calculation circuits. Here, an embodiment of the signal-to-noise power ratio improving method is also described.
【0060】上記の信号対雑音電力比改善装置の動作を
図7に基づいて説明する。ただし、入力信号x(n)は、実
施の形態1と同様、式(13)に示したように、振幅A、角
周波数ω、初期位相φの複素正弦波信号に、その複素正
弦波信号と独立で、平均ゼロ、分散σ2のガウス性雑音g
(n)が加わったものとしている。図7において、M個の4
次キュムラント計算回路2の入力までは上記実施の形態
1と同じである。すなわち、ある時間間隔で、サンプリ
ングされた入力信号x(n)は、シフトレジスタ1に順次入
力される。シフトレジスタ1では、Nサンプルの入力信
号を保持し、M個の4次キュムラント計算回路2に対し
て、Nサンプルの同じ入力信号を出力する。各4次キュ
ムラント計算回路2では、4次キュムラント関数の3変
数t1、t2、t3の中の一つを変数tとし、他の2変数をパ
ラメータとして、tを最大-(N-1)から(N-1)まで変化させ
て、M組の最大2N-1個の4次キュムラントを求める。す
なわち、例えばt1を変数tとした場合、他の2つの変数t
2、t3は、各4次キュムラント計算回路2ごとに、異な
る組み合わせの定数として、4次キュムラントを求め
る。その際に、実施の形態1では、t2、t3の値ai、bi
(i=1,2,…,M)は、各4次キュムラント計算回路2ごと
に、式(15)で示す関係を満たす異なった組合せの定数と
するという制約があった。それに対し、本実施例では、
ai、biは、各4次キュムラント計算回路2ごとに、異な
った組合せの定数であればよい。ただし、実施の形態1
と同様に、ai、biの値を入れ替えた値は同時には用いな
い。The operation of the above signal-to-noise power ratio improving apparatus will be described with reference to FIG. However, as in the first embodiment, the input signal x (n) is converted into a complex sine wave signal having an amplitude A, an angular frequency ω, and an initial phase φ as shown in Expression (13). Gaussian noise g with independent zero mean and variance σ 2
(n) has been added. In FIG. 7, M 4
The operation up to the input of the next cumulant calculation circuit 2 is the same as that of the first embodiment. That is, the input signal x (n) sampled at certain time intervals is sequentially input to the shift register 1. The shift register 1 holds the input signals of N samples, and outputs the same input signals of N samples to the M fourth-order cumulant calculation circuits 2. In each of the fourth-order cumulant calculation circuits 2, one of the three variables t1, t2, and t3 of the fourth-order cumulant function is set as a variable t, and the other two variables are set as parameters. N-1) to obtain M sets of up to 2N-1 quaternary cumulants. That is, for example, if t1 is a variable t, the other two variables t
2, t3 is to obtain the fourth cumulant as a constant of a different combination for each fourth cumulant calculation circuit 2. At this time, in the first embodiment, the values ai, bi of t2 and t3
(I = 1, 2,..., M) is restricted to a constant of a different combination that satisfies the relationship shown in Expression (15) for each fourth-order cumulant calculation circuit 2. In contrast, in the present embodiment,
ai and bi may be constants of different combinations for each fourth-order cumulant calculation circuit 2. However, Embodiment 1
Similarly to the above, the values obtained by replacing the values of ai and bi are not used at the same time.
【0061】このようにして求めたM組の4次キュムラ
ントは、それぞれ対応する位相補正回路9に入力され
る。ただし、式(11)にも示したように、4次キュムラン
トの振幅は、入力信号の振幅の4乗となるために、各位
相補正回路9へは、振幅の4乗根をとった4次キュムラ
ントを出力する。各位相補正回路9では、図7に示した
所望信号の周波数情報から得られる複素正弦波信号の角
周波数ωと、ai、biを用いて式(17)で求まる値を各4次
キュムラントに掛け、位相の補正を行う。なお、式(17)
は、t1を変数tとし、4キュムラント関数を式(1)に示し
たように、x(n+t2)とx(n+t3)を複素共役として求めた場
合に成り立つ関係であり、t2あるいはt3を変数tとした
場合や式(1)の他の位置を複素共役として4次キュムラ
ント関数を求めた場合は他の式となる。The M sets of fourth-order cumulants thus obtained are input to the corresponding phase correction circuits 9. However, since the amplitude of the fourth-order cumulant is the fourth power of the amplitude of the input signal as shown in equation (11), the fourth-order cumulant having the fourth root of the amplitude is applied to each phase correction circuit 9. Output cumulant. Each phase correction circuit 9 multiplies each quadratic cumulant by the angular frequency ω of the complex sine wave signal obtained from the frequency information of the desired signal shown in FIG. 7 and the value obtained by equation (17) using ai and bi. And phase correction. Equation (17)
Is the relationship that holds when x (n + t2) and x (n + t3) are obtained as complex conjugates, as shown in equation (1), where t1 is the variable t and the four cumulant functions are shown in equation (1). When t3 is a variable t, or when a fourth-order cumulant function is obtained with another position of the expression (1) as a complex conjugate, another expression is obtained.
【0062】[0062]
【数10】 (Equation 10)
【0063】このような処理を行うことによって、複素
正弦波信号のみを入力信号とした場合の4次キュムラン
ト関数(入力信号と同じ周波数の複素正弦波)の初期位
相を、aiとbiの値によらず一定にすることができる。す
なわち、複素正弦波信号のみを入力信号とした場合、M
個の4次キュムラント計算回路2で求まった、同じtに
対する4次キュムラントの値をすべて同じにすることが
できる。一方、平均ゼロ、分散σ2のガウス性雑音g(n)
のみを入力信号とした場合は、同じtに対する4次キュ
ムラントの値は、各4次キュムラント計算回路2ごとに
異なり、しかも、それらの和の平均はゼロに近づく。By performing such processing, the initial phase of the fourth-order cumulant function (a complex sine wave having the same frequency as the input signal) when only the complex sine wave signal is used as the input signal is changed to the values of ai and bi. It can be kept constant. That is, when only a complex sine wave signal is used as an input signal, M
The values of the fourth-order cumulants for the same t obtained by the four fourth-order cumulant calculation circuits 2 can be all the same. On the other hand, Gaussian noise g (n) with mean zero and variance σ 2
If only the input signal is used, the value of the fourth-order cumulant for the same t differs for each fourth-order cumulant calculation circuit 2, and the average of the sum thereof approaches zero.
【0064】このようにして求めたM組の位相補正回路
9の出力信号は加算回路7に出力される。加算回路7で
は、同じtごとにM組の位相補正回路9の出力信号を足
し、tごとの平均をとる。このような処理を行うことに
よって、複素正弦波信号成分に対する4次キュムラント
の値を変化させずに、ガウス性雑音成分に対する4次キ
ュムラントだけをゼロに近づけることができる。すなわ
ち、信号対雑音電力比が向上することができる。The output signals of the M sets of phase correction circuits 9 thus obtained are output to the addition circuit 7. The adder circuit 7 adds the output signals of the M sets of the phase correction circuits 9 for each of the same t, and takes an average for each of the t. By performing such a process, it is possible to make only the fourth-order cumulant for the Gaussian noise component close to zero without changing the value of the fourth-order cumulant for the complex sine wave signal component. That is, the signal-to-noise power ratio can be improved.
【0065】その後の処理は、上記実施の形態1と同様
に、加算回路7の出力信号は、窓関数回路3に入力され
る。窓関数回路3では、高速フーリエ変換回路4におけ
る、加算回路7の出力信号の数が有限であることの影響
を低減するために、窓関数による重みを掛け、その結果
を、高速フーリエ変換回路4に出力する。高速フーリエ
変換回路4では、窓関数回路3の出力信号を、加算回路
7の出力信号の数と入力信号x(n)のサンプリング時間で
決まる帯域ごとの成分に分別し、自乗検波回路5に出力
する。以降は従来例と同様に、自乗検波回路5におい
て、各帯域の信号の電力を求め、信号対雑音電力比計算
回路6に出力される。信号対雑音電力比計算回路6で
は、入力信号x(n)の複素正弦波信号の周波数が分かって
いる、あるいは予測している場合は、周波数の存在する
帯域の電力値から、周波数の存在する帯域以外の電力の
平均値を引いた値と、周波数の存在する帯域以外の電力
の平均値との比を求める。また、入力信号x(n)の複素正
弦波信号の周波数が分からない場合は、最も大きい電力
値からその他の帯域の電力の平均値を引いた値と、最も
大きい電力値を持つ帯域以外の帯域の電力の平均値との
比を求める。In the subsequent processing, the output signal of the adding circuit 7 is input to the window function circuit 3 as in the first embodiment. The window function circuit 3 multiplies the weight by the window function in order to reduce the influence of the finite number of output signals of the adder circuit 7 in the fast Fourier transform circuit 4, and divides the result by the fast Fourier transform circuit 4. Output to The fast Fourier transform circuit 4 separates the output signal of the window function circuit 3 into components for each band determined by the number of output signals of the adder circuit 7 and the sampling time of the input signal x (n), and outputs it to the square detection circuit 5. I do. Thereafter, similarly to the conventional example, the power of the signal in each band is obtained in the square detection circuit 5 and output to the signal-to-noise power ratio calculation circuit 6. In the signal-to-noise power ratio calculation circuit 6, when the frequency of the complex sine wave signal of the input signal x (n) is known or predicted, the frequency exists based on the power value of the frequency band. The ratio between the value obtained by subtracting the average value of the power outside the band and the average value of the power outside the band where the frequency exists is determined. If the frequency of the complex sine wave signal of the input signal x (n) is not known, a value obtained by subtracting the average value of the power of the other bands from the largest power value and a band other than the band having the largest power value And the ratio of the average power to the average value.
【0066】実施の形態6.図8は、本発明の実施の形
態6を示す信号対雑音電力比改善装置の構成図であり、
上記実施の形態2と異なる部分は、M個の4次キュムラ
ント計算回路2の後に、M個の位相補正回路9を設けた
ことである。なお、ここでは信号対雑音電力比改善方法
の実施の形態についても併せて記す。Embodiment 6 FIG. FIG. 8 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 6 of the present invention.
The difference from the second embodiment is that M phase correction circuits 9 are provided after M fourth-order cumulant calculation circuits 2. Here, an embodiment of the signal-to-noise power ratio improving method is also described.
【0067】上記の信号対雑音電力比改善装置の動作を
図8に基づいて説明する。ただし、入力信号x(n)は、実
施の形態2と同様、式(13)に示したように、振幅A、角
周波数ω、初期位相φの複素正弦波信号に、その複素正
弦波信号と独立で、平均ゼロ、分散σ2のガウス性雑音g
(n)が加わったものとしている。図8において、M個の位
相補正回路9の入力までは、上記実施の形態2と同じで
ある。すなわち、ある時間間隔で、サンプリングされた
入力信号x(n)は、シフトレジスタ1に順次入力される。
シフトレジスタ1では、Nサンプルの入力信号を保持
し、M個の4次キュムラント計算回路2に対して、Nサン
プルの同じ入力信号を出力する。各4次キュムラント計
算回路2では、4次キュムラント関数の3変数t1、t2、
t3の中の一つを変数tとし、他の2変数をパラメータと
して、tを最大-(N-1)から(N-1)まで変化させて、M組の
最大2N-1個の4次キュムラントを求める。すなわち、例
えばt1を変数tとした場合、他の2つの変数t2、t3は、
各4次キュムラント計算回路2ごとに、異なる組み合わ
せの定数として、4次キュムラントを求める。その際
に、実施の形態2では、t2、t3の値ai、bi(i=1,2,…,
M)は、各4次キュムラント計算回路2ごとに、式(15)
で示す関係を満たす異なった組合せの定数とするという
制約があった。それに対し、本実施の形態では、ai、bi
は、各4次キュムラント計算回路2ごとに、異なった組
合せの定数であればよい。ただし、実施の形態2と同様
に、ai、biの値を入れ替えた値は同時には用いない。The operation of the above signal-to-noise power ratio improving apparatus will be described with reference to FIG. However, as in the second embodiment, the input signal x (n) is a complex sine wave signal having an amplitude A, an angular frequency ω, and an initial phase φ as shown in Expression (13), and the complex sine wave signal Gaussian noise g with independent zero mean and variance σ 2
(n) has been added. In FIG. 8, the operation up to the input of M phase correction circuits 9 is the same as that of the second embodiment. That is, the input signal x (n) sampled at certain time intervals is sequentially input to the shift register 1.
The shift register 1 holds the input signals of N samples, and outputs the same input signals of N samples to the M fourth-order cumulant calculation circuits 2. In each fourth-order cumulant calculation circuit 2, three variables t1, t2,
Using one of t3 as a variable t and the other two variables as parameters, t is changed from-(N-1) to (N-1), and M sets of up to 2N-1 quartic Seek cumulant. That is, for example, if t1 is a variable t, the other two variables t2 and t3 are
For each quaternary cumulant calculation circuit 2, a quaternary cumulant is obtained as a constant of a different combination. At that time, in the second embodiment, the values ai, bi (i = 1, 2,.
M) is given by equation (15) for each fourth-order cumulant calculation circuit 2.
There is a restriction that constants of different combinations satisfy the relationship shown by. On the other hand, in the present embodiment, ai, bi
May be a constant of a different combination for each fourth-order cumulant calculation circuit 2. However, as in the second embodiment, the values obtained by replacing the values of ai and bi are not used at the same time.
【0068】このようにして求めたM組の4次キュムラ
ントは、それぞれ対応する位相補正回路9に入力され
る。ただし、式(11)にも示したように、4次キュムラン
トの振幅は、入力信号の振幅の4乗となるために、各位
相補正回路9へは、振幅の4乗根をとった4次キュムラ
ントを出力する。各位相補正回路9では、複素正弦波信
号の角周波数ωと、ai、biを用いて式(17)で求まる値を
各4次キュムラントに掛け、位相の補正を行う。このよ
うな処理を行うことによって、複素正弦波信号のみを入
力信号とした場合の4次キュムラント関数(入力信号と
同じ周波数の複素正弦波)の初期位相を、aiとbiの値に
よらず一定にすることができる。すなわち、複素正弦波
信号のみを入力信号とした場合、M個の4次キュムラン
ト計算回路2で求まる、同じtに対する4次キュムラン
トの値をすべて同じにすることができる。一方、平均ゼ
ロ、分散σ2のガウス性雑音g(n)のみを入力信号とした
場合は、同じtに対する4次キュムラントの値は、各4
次キュムラント計算回路2ごとに異なり、しかも、それ
らの和の平均はゼロに近づく。The M sets of fourth-order cumulants thus obtained are input to the corresponding phase correction circuits 9. However, since the amplitude of the fourth-order cumulant is the fourth power of the amplitude of the input signal as shown in equation (11), the fourth-order cumulant having the fourth root of the amplitude is applied to each phase correction circuit 9. Output cumulant. Each phase correction circuit 9 corrects the phase by multiplying each fourth-order cumulant by the value obtained by Expression (17) using the angular frequency ω of the complex sine wave signal and ai and bi. By performing such processing, the initial phase of the fourth-order cumulant function (a complex sine wave having the same frequency as the input signal) when only the complex sine wave signal is used as the input signal is constant regardless of the values of ai and bi. Can be That is, when only the complex sine wave signal is used as the input signal, the values of the fourth order cumulants for the same t obtained by the M fourth order cumulant calculation circuits 2 can be all the same. On the other hand, when only the Gaussian noise g (n) having zero mean and variance σ 2 is used as the input signal, the value of the fourth order cumulant for the same t is 4
The next cumulant calculation circuit 2 differs, and the average of the sum thereof approaches zero.
【0069】このようにして求めたM組の位相補正回路
9の出力信号はそれぞれ対応する窓関数回路3に入力さ
れる。各窓関数回路3では、各位相補正回路9の出力信
号に対して、高速フーリエ変換回路4における、4次キ
ュムラントの数が有限であることの影響を低減するため
に、窓関数の重みを掛け、その結果を、それぞれに対応
する高速フーリエ変換回路4に出力する。その後の処理
は、上記実施の形態2と同様に、各高速フーリエ変換回
路4において、それぞれの窓関数回路3の出力信号を、
位相補正回路9の出力信号の数と入力信号x(n)のサンプ
リング時間で決まる帯域ごとの成分に分別し、加算回路
7に出力する。加算回路7では、各帯域ごとにM個の信
号を加算し、各帯域ごとに平均値を求める。このような
処理を行うことによって、複素正弦波信号成分に対する
信号を変化させず、ガウス性雑音成分に対する信号だけ
をゼロに近づけることができ、結果として、信号対雑音
電力比を向上することができる。加算回路7の出力信号
は、自乗検波回路5に入力される。自乗検波回路5で
は、各帯域の信号の電力を求め、信号対雑音電力比計算
回路6に出力される。信号対雑音電力比計算回路6で
は、入力信号x(n)の複素正弦波信号の周波数が分かって
いる、あるいは予測している場合は、周波数の存在する
帯域の電力値から、周波数の存在する帯域以外の電力の
平均値を引いた値と、周波数の存在する帯域以外の電力
の平均値との比を求める。また、入力信号x(n)の複素正
弦波信号の周波数が分からない場合は、最も大きい電力
値からその他の帯域の電力の平均値を引いた値と、最も
大きい電力値を持つ帯域以外の帯域の電力の平均値との
比を求める。The output signals of the M sets of phase correction circuits 9 obtained as described above are input to the corresponding window function circuits 3. Each window function circuit 3 multiplies the output signal of each phase correction circuit 9 by a weight of a window function in order to reduce the influence of the finite number of fourth-order cumulants in the fast Fourier transform circuit 4. , And outputs the result to the corresponding fast Fourier transform circuit 4. Subsequent processing is similar to that of the second embodiment described above. In each fast Fourier transform circuit 4, the output signal of each window function circuit 3 is
The signal is separated into components for each band determined by the number of output signals of the phase correction circuit 9 and the sampling time of the input signal x (n), and output to the addition circuit 7. The adder circuit 7 adds M signals for each band and calculates an average value for each band. By performing such processing, the signal for the Gaussian noise component can be made close to zero without changing the signal for the complex sine wave signal component, and as a result, the signal-to-noise power ratio can be improved. . The output signal of the addition circuit 7 is input to the square detection circuit 5. The square detection circuit 5 obtains the power of the signal in each band and outputs it to the signal-to-noise power ratio calculation circuit 6. In the signal-to-noise power ratio calculation circuit 6, when the frequency of the complex sine wave signal of the input signal x (n) is known or predicted, the frequency exists based on the power value of the frequency band. The ratio between the value obtained by subtracting the average value of the power outside the band and the average value of the power outside the band where the frequency exists is determined. If the frequency of the complex sine wave signal of the input signal x (n) is not known, a value obtained by subtracting the average value of the power of the other bands from the largest power value and a band other than the band having the largest power value And the ratio of the average power to the average value.
【0070】実施の形態7.図9は、本発明の実施の形
態7を示す信号対雑音電力比改善装置の構成図であり、
上記実施の形態2と異なる部分は、M個の高速フーリエ
変換回路4の後に設けた、M個の自乗検波回路5と、加
算回路7の後の自乗検波回路5をなくしたことである。
なお、ここでは信号対雑音電力比改善方法の実施の形態
についても併せて記す。Embodiment 7 FIG. FIG. 9 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 7 of the present invention.
The difference from the second embodiment is that the M square detection circuits 5 provided after the M fast Fourier transform circuits 4 and the square detection circuit 5 provided after the addition circuit 7 are eliminated.
Here, an embodiment of the signal-to-noise power ratio improving method is also described.
【0071】上記の信号対雑音電力比改善装置の動作を
図9に基づいて説明する。ただし、入力信号x(n)は、実
施の形態2と同様、式(13)に示したように、振幅A、角
周波数ω、初期位相φの複素正弦波信号に、その複素正
弦波信号と独立で、平均ゼロ、分散σ2のガウス性雑音g
(n)が加わったものとしている。図9において、M個の4
次キュムラント計算回路2の出力までは、上記実施の形
態2と同じである。すなわち、ある時間間隔で、サンプ
リングされた入力信号x(n)は、シフトレジスタ1に順次
入力される。シフトレジスタ1では、Nサンプルの入力
信号を保持し、M個の4次キュムラント計算回路2に対
して、Nサンプルの同じ入力信号を出力する。各4次キ
ュムラント計算回路2では、4次キュムラント関数の3
変数t1、t2、t3の中の一つを変数tとし、他の2変数を
パラメータとして、tを最大-(N-1)から(N-1)まで変化さ
せて、M組の最大2N-1個の4次キュムラントを求める。
すなわち、例えばt1を変数tとした場合、他の2つの変
数t2、t3は、各4次キュムラント計算回路2ごとに、異
なる組み合わせの定数として、4次キュムラントを求め
る。その際に、実施の形態2では、t2、t3の値ai、bi
(i=1,2,…,M)は、各4次キュムラント計算回路2ごと
に、式(15)で示す関係を満たす異なった組合せの定数と
するという制約があった。それに対し、本実施の形態で
は、ai、biは、各4次キュムラント計算回路2ごとに、
異なった組合せの定数であればよい。ただし、実施の形
態2と同様に、ai、biの値を入れ替えた値は同時には用
いない。The operation of the signal to noise power ratio improving apparatus will be described with reference to FIG. However, as in the second embodiment, the input signal x (n) is a complex sine wave signal having an amplitude A, an angular frequency ω, and an initial phase φ as shown in Expression (13), and the complex sine wave signal Gaussian noise g with independent zero mean and variance σ 2
(n) has been added. In FIG. 9, M 4
The operation up to the output of the next cumulant calculation circuit 2 is the same as that of the second embodiment. That is, the input signal x (n) sampled at certain time intervals is sequentially input to the shift register 1. The shift register 1 holds the input signals of N samples, and outputs the same input signals of N samples to the M fourth-order cumulant calculation circuits 2. In each fourth-order cumulant calculation circuit 2, the third-order cumulant function 3
One of the variables t1, t2, and t3 is set as a variable t, and the other two variables are set as parameters, and t is changed from-(N-1) to (N-1) to obtain a maximum of 2N- Find one quaternary cumulant.
That is, for example, when t1 is a variable t, the other two variables t2 and t3 determine the fourth cumulant as constants of different combinations for each fourth cumulant calculation circuit 2. At this time, in the second embodiment, the values ai, bi of t2 and t3
(I = 1, 2,..., M) is restricted to a constant of a different combination that satisfies the relationship shown in Expression (15) for each fourth-order cumulant calculation circuit 2. On the other hand, in the present embodiment, ai and bi are calculated for each fourth-order cumulant calculation circuit 2,
What is necessary is just a constant of a different combination. However, as in the second embodiment, the values obtained by replacing the values of ai and bi are not used at the same time.
【0072】このようにして求めたM組の4次キュムラ
ントは、それぞれに対応する窓関数回路3に出力され
る。ただし、式(11)にも示したように、4次キュムラン
トの振幅は、入力信号の振幅の4乗となるために、各窓
関数回路3へは、振幅の4乗根をとった4次キュムラン
トを出力する。各窓関数回路3では、各4次キュムラン
ト計算回路2で求まった4次キュムラントに対して、高
速フーリエ変換回路4における、4次キュムラントの数
が有限であることの影響を低減するために、窓関数の重
みを掛け、その結果を、それぞれに対応する高速フーリ
エ変換回路4に出力する。各高速フーリエ変換回路4で
は、それぞれの窓関数回路3の出力信号を、4次キュム
ラント関数値の数と入力信号x(n)のサンプリング時間で
決まる帯域ごとの成分に分別する。実施の形態2では、
各高速フーリエ変換回路4の出力信号は、加算回路7に
出力されていたのに対し、本実施の形態では、高速フー
リエ変換回路4の出力信号は、それぞれ対応する自乗検
波回路5に出力される。各自乗検波回路5では、各帯域
の信号を初期位相に無関係な電力値に変換し、その結果
を加算回路7に出力する。加算回路7では、各帯域ごと
にM個の信号を足し、平均をとる。このような処理を行
うことによって、複素正弦波信号成分に対する信号を変
化させず、ガウス性雑音成分に対する信号の分散のみを
小さくすることができ、前記式(14)から明らかなよう
に、結果として、信号対雑音電力比を向上することがで
きる。加算回路7の出力信号は、信号対雑音電力比計算
回路6に出力され、信号対雑音電力比計算回路6では、
入力信号x(n)の複素正弦波信号の周波数が分かってい
る、あるいは予測している場合は、周波数の存在する帯
域の電力値から、周波数の存在する帯域以外の電力の平
均値を引いた値と、周波数の存在する帯域以外の電力の
平均値との比を求める。また、入力信号x(n)の複素正弦
波信号の周波数が分からない場合は、最も大きい電力値
からその他の帯域の電力の平均値を引いた値と、最も大
きい電力値を持つ帯域以外の帯域の電力の平均値との比
を求める。The M sets of fourth-order cumulants determined in this way are output to the corresponding window function circuits 3. However, as shown in equation (11), the amplitude of the fourth-order cumulant is the fourth power of the amplitude of the input signal. Output cumulant. In each window function circuit 3, in order to reduce the influence of the finite number of fourth-order cumulants in the fast Fourier transform circuit 4 on the fourth-order cumulants obtained in each fourth-order cumulant calculation circuit 2, The weight of the function is multiplied, and the result is output to the corresponding fast Fourier transform circuit 4. Each fast Fourier transform circuit 4 classifies the output signal of each window function circuit 3 into components for each band determined by the number of fourth-order cumulant function values and the sampling time of the input signal x (n). In the second embodiment,
While the output signal of each fast Fourier transform circuit 4 has been output to the addition circuit 7, in the present embodiment, the output signal of the fast Fourier transform circuit 4 is output to the corresponding square detection circuit 5, respectively. . Each square detection circuit 5 converts the signal of each band into a power value irrelevant to the initial phase, and outputs the result to the addition circuit 7. The adding circuit 7 adds M signals for each band and takes an average. By performing such a process, it is possible to reduce only the variance of the signal with respect to the Gaussian noise component without changing the signal with respect to the complex sine wave signal component, and as a result, as is apparent from the equation (14), , The signal to noise power ratio can be improved. The output signal of the addition circuit 7 is output to the signal-to-noise power ratio calculation circuit 6, and the signal-to-noise power ratio calculation circuit 6
If the frequency of the complex sine wave signal of the input signal x (n) is known or predicted, the average value of the power outside the frequency band is subtracted from the power value of the frequency band. The ratio between the value and the average value of the power other than the band in which the frequency exists is determined. If the frequency of the complex sine wave signal of the input signal x (n) is not known, a value obtained by subtracting the average value of the power of the other bands from the largest power value and a band other than the band having the largest power value And the ratio of the average power to the average value.
【0073】なお、以上の実施の形態では、信号対雑音
電力比改善装置および信号対雑音電力比改善方法におい
て、4次キュムラント計算回路などをM個並列に設けた
装置例またはM組の並列的な処理例を示して説明した
が、本願発明では、ハード的にM個存在することや、並
列的な処理をすることは必須ではなく、ソフト的に順次
繰返し処理する装置および方法でも良いことは言うまで
もない。In the above embodiment, in the signal-to-noise power ratio improving apparatus and the signal-to-noise power ratio improving method, an example of an apparatus in which M fourth-order cumulant calculation circuits and the like are provided in parallel or M sets of parallel circuits Although the present invention has been described with reference to a simple processing example, in the present invention, it is not essential that there are M pieces of hardware, and that parallel processing is not essential, and that an apparatus and a method of sequentially and repeatedly processing software may be used. Needless to say.
【0074】[0074]
【発明の効果】請求項1又は2記載の発明によれば、入
力信号中の所望信号である周期性信号については、4次
キュムラント関数の3変数の中の1つを変数、他の2つ
をそれぞれパラメータとし、上記2つのパラメータを入
力信号中の所望信号に対する4次キュムラント関数の初
期位相がすべて同じになるようにした複数個の異なる2
つのパラメータの組み合わせのそれぞれに対し、上記変
数を変化させて4次キュムラントの組を求めると、4次
キュムラント関数の初期位相が一定となるため同じ変数
に対する4次キュムラントの値をすべて同じにすること
ができる。また、ガウス性雑音については同じ変数に対
する4次キュムラントの値は、複数個の異なる2つのパ
ラメータの組み合わせごとに異なり、しかも、それらの
値の平均はゼロに近づくので、4次キュムラントの組の
それぞれについて、同じ変数ごとに加算処理してその平
均を求めると信号対雑音電力比を改善させる効果があ
る。According to the first or second aspect of the present invention, for a periodic signal which is a desired signal in an input signal, one of three variables of the fourth-order cumulant function is a variable and the other two are variables. Are used as parameters, and the above two parameters are different from each other in such a manner that the initial phases of the fourth-order cumulant function for the desired signal in the input signal are all the same.
For each combination of the three parameters, when the above variables are changed to obtain a set of fourth-order cumulants, the initial phase of the fourth-order cumulant function becomes constant, so that the values of the fourth-order cumulants for the same variables are all the same. Can be. For Gaussian noise, the value of the fourth-order cumulant for the same variable differs for each combination of a plurality of different two parameters, and the average of those values approaches zero. Is obtained by adding each of the same variables and calculating the average, which has the effect of improving the signal-to-noise power ratio.
【0075】請求項3又は4記載の発明によれば、入力
信号中の所望信号である周期性信号については、4次キ
ュムラントの組のそれぞれについて、各4次キュムラン
トに各組の2つのパラメータに基づいて得られる量の関
数シフトを与えると、上記2つのパラメータを入力信号
中の所望信号に対する4次キュムラント関数の初期位相
がすべて同じになるように設定せずとも4次キュムラン
ト関数の初期位相が一定となるため、同じ変数に対する
4次キュムラントの値をすべて同じにすることができ
る。従って、請求項1又は2記載の発明と同様に4次キ
ュムラントの組のそれぞれについて、同じ変数ごとに加
算処理してその平均を求めると信号対雑音電力比を改善
させる効果がある。According to the third or fourth aspect of the present invention, for a periodic signal which is a desired signal in an input signal, for each set of fourth-order cumulants, each fourth-order cumulant has two parameters of each set. If the initial phase of the fourth-order cumulant function is not set so that the initial phases of the fourth-order cumulant function for the desired signal in the input signal are all the same, then the initial phase of the fourth-order cumulant function is given. Since it is constant, the values of the fourth-order cumulants for the same variable can all be the same. Therefore, similar to the first or second aspect of the present invention, if the addition processing is performed for each of the same variables for each set of the fourth-order cumulants and the average is obtained, there is an effect of improving the signal-to-noise power ratio.
【0076】請求項5又は6記載の発明によれば、入力
信号中の所望信号である周期性信号については、4次キ
ュムラントの組のそれぞれについて、各4次キュムラン
トに各組の2つのパラメータと上記入力信号中の所望信
号の周波数情報に基づいて位相補正を行うと、上記2つ
のパラメータを入力信号中の所望信号に対する4次キュ
ムラント関数の初期位相がすべて同じになるように設定
せずとも4次キュムラント関数の初期位相が一定となる
ため、同じ変数に対する4次キュムラントの値をすべて
同じにすることができる。従って、請求項1又は2記載
の発明と同様に4次キュムラントの組のそれぞれについ
て、同じ変数ごとに加算処理してその平均を求めると信
号対雑音電力比を改善させる効果がある。According to the fifth or sixth aspect of the present invention, for a periodic signal that is a desired signal in the input signal, for each set of fourth-order cumulants, each of the four-order cumulants has two parameters of each set. When the phase correction is performed based on the frequency information of the desired signal in the input signal, the four parameters need not be set so that the initial phases of the fourth-order cumulant functions for the desired signal in the input signal are all the same. Since the initial phase of the second-order cumulant function is constant, the values of the fourth-order cumulants for the same variable can be all the same. Therefore, similar to the first or second aspect of the present invention, if the addition processing is performed for each of the same variables for each set of the fourth-order cumulants and the average is obtained, there is an effect of improving the signal-to-noise power ratio.
【0077】請求項7又は8記載の発明によれば、同じ
変数ごとに加算処理してその平均を求める手段又は手順
に代えて、フーリエ変換処理により周波数成分に変換し
た後に、同じ帯域ごとに周波数領域で加算処理してその
平均を求めても、請求項1〜6記載の発明同様の信号対
雑音電力比を改善させる効果がある。According to the seventh or eighth aspect of the present invention, instead of the means or procedure for adding and processing the average for each of the same variables and obtaining the average thereof, the frequency is converted into frequency components by the Fourier transform processing, and then the frequency is changed for each of the same bands. Even if the addition processing is performed in the area and the average is obtained, there is an effect of improving the signal-to-noise power ratio as in the first to sixth aspects of the present invention.
【0078】請求項9又は10記載の発明によれば、4
次キュムラントの組のそれぞれについて、各4次キュム
ラントをフーリエ変換処理により周波数成分に変換し、
包絡線検波処理を行い、同じ帯域ごとに加算処理してそ
の平均を求めると、入力信号中の所望信号である周期性
信号については、2つのパラメータを入力信号中の所望
信号に対する4次キュムラント関数の初期位相がすべて
同じになるようにせずとも4次キュムラント関数の初期
位相に関係しない値で扱うことができ、また、ガウス性
雑音成分に対する信号の分散を小さくできるので、信号
対雑音電力比を改善させる効果がある。According to the ninth or tenth aspect of the present invention, 4
For each set of next cumulants, each fourth order cumulant is converted to a frequency component by Fourier transform processing,
When the envelope detection process is performed, the addition process is performed for each of the same bands, and the average is obtained. As for the periodic signal which is the desired signal in the input signal, two parameters are calculated by the fourth order cumulant function for the desired signal in the input signal. Can be treated as a value not related to the initial phase of the fourth-order cumulant function even if the initial phases are not all the same, and the variance of the signal with respect to the Gaussian noise component can be reduced. It has the effect of improving.
【図1】 この発明の実施の形態1の信号対雑音電力比
改善装置の構成図である。FIG. 1 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 1 of the present invention.
【図2】 従来の信号対雑音電力比改善装置による自乗
検波回路出力信号の一例を示す図である。FIG. 2 is a diagram illustrating an example of an output signal of a square detection circuit by a conventional signal-to-noise power ratio improving apparatus.
【図3】 この発明の実施の形態1の信号対雑音電力比
改善装置による自乗検波回路出力信号の一例を示す図で
ある。FIG. 3 is a diagram illustrating an example of an output signal of a square detection circuit by the signal-to-noise power ratio improving apparatus according to the first embodiment of the present invention;
【図4】 この発明の実施の形態2の信号対雑音電力比
改善装置の構成図である。FIG. 4 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 2 of the present invention.
【図5】 この発明の実施の形態3の信号対雑音電力比
改善装置の構成図である。FIG. 5 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 3 of the present invention.
【図6】 この発明の実施の形態4の信号対雑音電力比
改善装置の構成図である。FIG. 6 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 4 of the present invention.
【図7】 この発明の実施の形態5の信号対雑音電力比
改善装置の構成図である。FIG. 7 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to a fifth embodiment of the present invention.
【図8】 この発明の実施の形態6の信号対雑音電力比
改善装置の構成図である。FIG. 8 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 6 of the present invention.
【図9】 この発明の実施の形態7の信号対雑音電力比
改善装置の構成図である。FIG. 9 is a configuration diagram of a signal-to-noise power ratio improving apparatus according to Embodiment 7 of the present invention.
【図10】 従来の信号対雑音電力比改善装置の構成図
である。FIG. 10 is a configuration diagram of a conventional signal-to-noise power ratio improving apparatus.
1 シフトレジスタ、2 4次キュムラント計算回路、
3 窓関数回路、4 高速フーリエ変換回路、5 自乗
検波回路、6 信号対雑音電力比計算回路、7 加算回
路、8 関数シフト回路、9 位相補正回路。1 shift register, 4th order cumulant calculation circuit,
3 window function circuit, 4 fast Fourier transform circuit, 5 square detection circuit, 6 signal-to-noise power ratio calculation circuit, 7 addition circuit, 8 function shift circuit, 9 phase correction circuit.
フロントページの続き (56)参考文献 特開 平6−21838(JP,A) 特開 平5−14241(JP,A) 米国特許5337053(US,A) ”The general sour ce separation prob lem”Lacoume,J.L;Ga eta,M.Spectrum Est imation and Modeli ng,1990.,Fifth ASSP Workshop on,1990 Pag e(s)154−158 ”Cumulant techniq ues for the estima tion and cancellat ion of intersymbol interference”Hatz inakos,D.;Nikias, C.L.Spectrum Estim ation and Modelin g,1988.,Fourth ASSP Workshop on,1988 Pag e(s)204−207 (58)調査した分野(Int.Cl.7,DB名) G01S 7/00 - 7/42 H03H 21/00 H04B 15/00 Continuation of the front page (56) References JP-A-6-21838 (JP, A) JP-A-5-14241 (JP, A) U.S. Pat. No. 5,370,053 (US, A) "The general source separation probe lemme" Lacome, J. L; Ga eta, M .; Spectrum Estimation and Modeling, 1990. , Fifth ASSP Workshopon, 1990 Page (s) 154-158 "Cumulant techniques for the estima- tion and cancelation of intersymbol interference. Nikias, C .; L. Spectrum Estimation and Modeling, 1988. , Fourth ASPS Works on, 1988 Page (s) 204-207 (58) Fields investigated (Int. Cl. 7 , DB name) G01S 7 /00-7/42 H03H 21/00 H04B 15/00
Claims (10)
信号対雑音電力比改善装置において、入力信号をサンプ
リングして一組の入力信号を生成する入力信号取り込み
手段と、上記入力信号取り込み手段からの同一の一組の
入力信号について、4次キュムラント関数の3変数の中
の1つを変数、他の2つをそれぞれパラメータとし、上
記2つのパラメータを入力信号中の所望信号に対する4
次キュムラント関数の初期位相がすべて同じになるよう
にした複数個の異なる2つのパラメータの組み合わせの
それぞれに対し、上記変数を変化させて4次キュムラン
トの組を求める4次キュムラント計算処理手段と、上記
4次キュムラント計算処理手段からの4次キュムラント
の組のそれぞれについて、同じ変数ごとに加算処理して
その平均を求める加算処理手段とを備えたことを特徴と
する信号対雑音電力比改善装置。1. A signal-to-noise power ratio improving apparatus for suppressing a Gaussian noise component of an input signal, comprising: an input signal capturing means for sampling an input signal to generate a set of input signals; , One of the three variables of the fourth-order cumulant function as a variable and the other two as parameters, respectively, and the above two parameters are defined as 4 to the desired signal in the input signal.
Fourth-order cumulant calculation processing means for obtaining a set of fourth-order cumulants by changing the above variables for each of a plurality of different combinations of two parameters such that the initial phases of the second-order cumulant functions are all the same; A signal-to-noise-to-noise ratio improvement device, comprising: addition processing means for adding each of the sets of fourth-order cumulants from the fourth-order cumulant calculation processing means for each of the same variables to obtain an average thereof.
信号対雑音電力比改善方法において、入力信号をサンプ
リングして一組の入力信号を生成し、同一の上記一組の
入力信号について、4次キュムラント関数の3変数の中
の1つを変数、他の2つをそれぞれパラメータとし、上
記2つのパラメータを入力信号中の所望信号に対する4
次キュムラント関数の初期位相がすべて同じになるよう
にした複数個の異なる2つのパラメータの組み合わせの
それぞれに対し、上記変数を変化させて4次キュムラン
トの組を求め、上記4次キュムラントの組のそれぞれに
ついて、同じ変数ごとに加算処理してその平均を求める
ことを特徴とする信号対雑音電力比改善方法。2. A signal-to-noise power ratio improving method for suppressing a Gaussian noise component of an input signal, wherein the input signal is sampled to generate a set of input signals. One of the three variables of the next cumulant function is a variable, and the other two are parameters, respectively, and the above two parameters are set to 4 for the desired signal in the input signal.
For each of a plurality of different combinations of two parameters in which the initial phase of the second order cumulant function is all the same, the above variables are changed to obtain a set of fourth order cumulants. A signal-to-noise power ratio improving method, wherein an addition process is performed for each of the same variables and an average thereof is obtained.
信号対雑音電力比改善装置において、入力信号をサンプ
リングして一組の入力信号を生成する入力信号取り込み
手段と、上記入力信号取り込み手段からの同一の一組の
入力信号について、4次キュムラント関数の3変数の中
の1つを変数、他の2つをそれぞれパラメータとし、複
数個の異なる2つのパラメータの組み合わせのそれぞれ
に対し、上記変数を変化させて4次キュムラントの組を
求める4次キュムラント計算処理手段と、上記4次キュ
ムラント計算処理手段からの4次キュムラントの組のそ
れぞれについて、各4次キュムラントに各組の上記2つ
のパラメータに基づいて得られる量の関数シフトを与え
る関数シフト処理手段と、上記関数シフト処理手段から
の出力を同じ変数ごとに加算処理してその平均を求める
加算処理手段とを備えたことを特徴とする信号対雑音電
力比改善装置。3. A signal-to-noise power ratio improving apparatus for suppressing a Gaussian noise component of an input signal, comprising: an input signal capturing means for sampling an input signal to generate a set of input signals; , One of the three variables of the fourth-order cumulant function is a variable, and the other two are parameters, respectively. For each combination of a plurality of different two parameters, Is changed to obtain a set of quaternary cumulants, and for each set of quaternary cumulants from the quaternary cumulant calculation processing means, each quaternary cumulant is assigned to each of the above two parameters. Function shift processing means for giving the amount of function shift obtained based on the same variable and the output from the function shift processing means for the same variable. A signal-to-noise-power ratio improving apparatus, comprising: addition processing means for obtaining an average by performing addition processing on the signals.
信号対雑音電力比改善方法において、入力信号をサンプ
リングして一組の入力信号を生成し、同一の上記一組の
入力信号について、4次キュムラント関数の3変数の中
の1つを変数、他の2つをそれぞれパラメータとし、複
数個の異なる2つのパラメータの組み合わせのそれぞれ
に対し、上記変数を変化させて4次キュムラントの組を
求め、上記4次キュムラントの組のそれぞれについて、
各4次キュムラントに各組の上記2つのパラメータに基
づいて得られる量の関数シフトを与え、同じ変数ごとに
加算処理してその平均を求めることを特徴とする信号対
雑音電力比改善方法。4. A signal-to-noise power ratio improving method for suppressing a Gaussian noise component of an input signal, wherein the input signal is sampled to generate a set of input signals. One of the three variables of the next cumulant function is a variable, and the other two are parameters, respectively. For each of a plurality of different combinations of two parameters, the above variables are changed to obtain a fourth-order cumulant set. , For each of the above set of quaternary cumulants,
A signal-to-noise power ratio improvement method, characterized in that each fourth-order cumulant is provided with a function shift of an amount obtained based on each of the above two parameters of each set, and an addition process is performed for each same variable to obtain an average thereof.
信号対雑音電力比改善装置において、入力信号をサンプ
リングして一組の入力信号を生成する入力信号取り込み
手段と、上記入力信号取り込み手段からの同一の一組の
入力信号について、4次キュムラント関数の3変数の中
の1つを変数、他の2つをそれぞれパラメータとし、複
数個の異なる2つのパラメータの組み合わせのそれぞれ
に対し、上記変数を変化させて4次キュムラントの組を
求める4次キュムラント計算処理手段と、上記4次キュ
ムラント計算処理手段からの4次キュムラントの組のそ
れぞれについて、各4次キュムラントに各組の上記2つ
のパラメータと上記入力信号中の所望信号の周波数情報
に基づいて位相補正を行う位相補償処理手段と、上記位
相処理手段からの出力を同じ変数ごとに加算処理してそ
の平均を求める加算処理手段とを備えたことを特徴とす
る信号対雑音電力比改善装置。5. A signal-to-noise power ratio improving apparatus for suppressing a Gaussian noise component of an input signal, comprising: an input signal capturing means for sampling an input signal to generate a set of input signals; , One of the three variables of the fourth-order cumulant function is a variable, and the other two are parameters, respectively. For each combination of a plurality of different two parameters, And a fourth-order cumulant calculation processing means for obtaining a set of fourth-order cumulants by changing the above-mentioned parameters. For each of the fourth-order cumulant sets from the fourth-order cumulant calculation processing means, Phase compensation processing means for performing phase correction based on frequency information of a desired signal in the input signal, and output from the phase processing means. A signal-to-noise-to-noise-ratio improving apparatus, comprising: addition processing means for performing addition processing on the same variable for each variable and obtaining an average thereof.
信号対雑音電力比改善方法において、入力信号をサンプ
リングして一組の入力信号を生成し、同一の上記一組の
入力信号について、4次キュムラント関数の3変数の中
の1つを変数、他の2つをそれぞれパラメータとし、複
数個の異なる2つのパラメータの組み合わせのそれぞれ
に対し、上記変数を変化させて4次キュムラントの組を
求め、上記4次キュムラントの組のそれぞれについて、
各4次キュムラントに各組の上記2つのパラメータと上
記入力信号中の所望信号の周波数情報に基づいて位相補
正を行い、同じ変数ごとに加算処理してその平均を求め
ることを特徴とする信号対雑音電力比改善方法。6. A signal-to-noise power ratio improving method for suppressing a Gaussian noise component of an input signal, wherein the input signal is sampled to generate a set of input signals. One of the three variables of the next cumulant function is a variable, and the other two are parameters, respectively. For each of a plurality of different combinations of two parameters, the above variables are changed to obtain a fourth-order cumulant set. , For each of the above set of quaternary cumulants,
A signal pair comprising: performing a phase correction on each fourth-order cumulant based on the above two parameters of each set and frequency information of a desired signal in the input signal, and performing an addition process for each same variable to obtain an average thereof; Noise power ratio improvement method.
音電力比改善装置において、加算処理手段に代えて、上
記加算処理手段へのそれぞれの組の入力を周波数成分に
変換して出力するフーリエ変換処理手段と、上記フーリ
エ変換処理手段からの出力を同じ帯域ごとに加算処理し
てその平均を求める周波数領域加算処理手段とを備えた
ことを特徴とする信号対雑音電力比改善装置。7. A signal-to-noise power ratio improving apparatus according to claim 1, wherein each set of inputs to said addition processing means is converted into a frequency component and output instead of said addition processing means. A signal-to-noise power ratio improving apparatus, comprising: a Fourier transform processing means for performing the above processing; and a frequency domain addition processing means for adding an output from the Fourier transform processing means for each same band to obtain an average thereof.
音電力比改善方法において、同じ変数ごとに加算処理し
てその平均を求める手順に代えて、フーリエ変換処理に
より周波数成分に変換し、同じ帯域ごとに加算処理して
その平均を求めることを特徴とする信号対雑音電力比改
善方法。8. The method for improving a signal-to-noise ratio according to claim 2, 4 or 6, wherein a frequency component is converted to a frequency component by a Fourier transform instead of a procedure of adding each variable and calculating an average thereof. A signal-to-noise power ratio improving method, wherein an addition process is performed for each of the same bands to obtain an average thereof.
信号対雑音電力比改善装置において、入力信号をサンプ
リングして一組の入力信号を生成する入力信号取り込み
手段と、上記入力信号取り込み手段からの同一の一組の
入力信号について、4次キュムラント関数の3変数の中
の1つを変数、他の2つをそれぞれパラメータとし、複
数個の異なる2つのパラメータの組み合わせのそれぞれ
に対し、上記変数を変化させて4次キュムラントの組を
求める4次キュムラント計算処理手段と、上記4次キュ
ムラント計算処理手段からの4次キュムラントの組のそ
れぞれについて、各4次キュムラントを周波数成分に変
換して出力するフーリエ変換処理手段と、上記フーリエ
変換処理手段からの出力を包絡線検波処理する包絡線検
波処理手段と、上記包絡線検波処理手段からの出力を同
じ帯域ごとに加算処理してその平均を求める加算処理手
段とを備えたことを特徴とする信号対雑音電力比改善装
置。9. A signal-to-noise power ratio improving apparatus for suppressing a Gaussian noise component of an input signal, comprising: an input signal capturing means for sampling an input signal to generate a set of input signals; , One of the three variables of the fourth-order cumulant function is a variable, and the other two are parameters, respectively. For each combination of a plurality of different two parameters, Is changed to obtain a set of fourth-order cumulants, and for each set of fourth-order cumulants from the fourth-order cumulant calculation processing means, each fourth-order cumulant is converted into a frequency component and output. Fourier transform processing means, envelope detection processing means for performing envelope detection processing on the output from the Fourier transform processing means, A signal-to-noise power ratio improving apparatus, comprising: an addition processing means for adding an output from the envelope detection processing means for each of the same bands and obtaining an average thereof.
る信号対雑音電力比改善方法において、入力信号をサン
プリングして一組の入力信号を生成し、同一の上記一組
の入力信号について、4次キュムラント関数の3変数の
中の1つを変数、他の2つをそれぞれパラメータとし、
複数個の異なる2つのパラメータの組み合わせのそれぞ
れに対し、上記変数を変化させて4次キュムラントの組
を求め、上記4次キュムラントの組のそれぞれについ
て、各4次キュムラントをフーリエ変換処理により周波
数成分に変換し、包絡線検波処理を行い、同じ帯域ごと
に加算処理してその平均を求めることを特徴とする信号
対雑音電力比改善方法。10. A signal-to-noise power ratio improving method for suppressing a Gaussian noise component of an input signal, wherein the input signal is sampled to generate a set of input signals. One of the three variables of the next cumulant function is a variable, and the other two are parameters,
For each of a plurality of different combinations of two parameters, the above variables are changed to obtain a fourth-order cumulant set, and for each of the fourth-order cumulant sets, each fourth-order cumulant is converted into a frequency component by Fourier transform processing. A signal-to-noise power ratio improving method characterized by performing conversion, performing envelope detection processing, performing addition processing for each of the same bands, and calculating an average thereof.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP27563895A JP3183128B2 (en) | 1995-10-24 | 1995-10-24 | Signal to noise power ratio improving apparatus and signal to noise power ratio improving method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP27563895A JP3183128B2 (en) | 1995-10-24 | 1995-10-24 | Signal to noise power ratio improving apparatus and signal to noise power ratio improving method |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH09119980A JPH09119980A (en) | 1997-05-06 |
JP3183128B2 true JP3183128B2 (en) | 2001-07-03 |
Family
ID=17558256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP27563895A Expired - Fee Related JP3183128B2 (en) | 1995-10-24 | 1995-10-24 | Signal to noise power ratio improving apparatus and signal to noise power ratio improving method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3183128B2 (en) |
-
1995
- 1995-10-24 JP JP27563895A patent/JP3183128B2/en not_active Expired - Fee Related
Non-Patent Citations (2)
Title |
---|
"Cumulant techniques for the estimation and cancellation of intersymbol interference"Hatzinakos,D.;Nikias,C.L.Spectrum Estimation and Modeling,1988.,Fourth ASSP Workshop on,1988 Page(s)204−207 |
"The general source separation problem"Lacoume,J.L;Gaeta,M.Spectrum Estimation and Modeling,1990.,Fifth ASSP Workshop on,1990 Page(s)154−158 |
Also Published As
Publication number | Publication date |
---|---|
JPH09119980A (en) | 1997-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7590528B2 (en) | Method and apparatus for noise suppression | |
US7099417B2 (en) | Trace video filtering using wavelet de-noising techniques | |
US6878113B2 (en) | Ultrasonic diagnostic apparatus and method for processing ultrasonic signal | |
US5874916A (en) | Frequency selective TDOA/FDOA cross-correlation | |
US7170959B1 (en) | Tailored response cascaded integrator comb digital filter and methodology for parallel integrator processing | |
US5422912A (en) | Adaptive weak signal identification system | |
KR980012873A (en) | Spectrum direct spread signal receiver and acquisition circuit | |
CN116148802B (en) | Bistatic SAR phase synchronization anti-interference method | |
JP2797616B2 (en) | Noise suppression device | |
JP3183128B2 (en) | Signal to noise power ratio improving apparatus and signal to noise power ratio improving method | |
CN116930610A (en) | PCAL signal phase amplitude estimation method based on coherent accumulation | |
US20030149711A1 (en) | Filter bank approach to adaptive filtering method using independent component analysis | |
Narasimhan et al. | Improved Wigner–Ville distribution performance by signal decomposition and modified group delay | |
JP3134618B2 (en) | Ultrasonic signal processor | |
van der Byl et al. | Recursive sliding discrete Fourier transform with oversampled data | |
JP4344356B2 (en) | Detector, method, program, recording medium | |
JP3146093B2 (en) | Two-stage fast Fourier transform method | |
US20230393184A1 (en) | Device and methods for phase noise measurement | |
David et al. | Fast sequential LS estimation for sinusoidal modeling and decomposition of audio signals | |
WO2023234963A1 (en) | Device and methods for phase noise measurement | |
JP2789600B2 (en) | Frequency estimation method | |
JPH08328593A (en) | Spectrum analysis method | |
JP2837614B2 (en) | Unwanted signal suppression device | |
JP2006173920A (en) | Smart antenna assembly and compound smart antenna assembly using the same | |
RU2099720C1 (en) | Digital spectrum analyzer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080427 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090427 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100427 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100427 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110427 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120427 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120427 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130427 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130427 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140427 Year of fee payment: 13 |
|
LAPS | Cancellation because of no payment of annual fees |