WO2007126054A1 - 信号分離装置、信号分離方法、情報記録媒体、ならびに、プログラム - Google Patents

信号分離装置、信号分離方法、情報記録媒体、ならびに、プログラム Download PDF

Info

Publication number
WO2007126054A1
WO2007126054A1 PCT/JP2007/059172 JP2007059172W WO2007126054A1 WO 2007126054 A1 WO2007126054 A1 WO 2007126054A1 JP 2007059172 W JP2007059172 W JP 2007059172W WO 2007126054 A1 WO2007126054 A1 WO 2007126054A1
Authority
WO
WIPO (PCT)
Prior art keywords
matrix
columns
rows
row
column
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.)
Ceased
Application number
PCT/JP2007/059172
Other languages
English (en)
French (fr)
Inventor
Andrzej Cichocki
Rafal Zdunek
Shunichi Amari
Gen Hori
Ken Umeno
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
RIKEN
Original Assignee
RIKEN
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by RIKEN filed Critical RIKEN
Priority to US12/226,651 priority Critical patent/US8285773B2/en
Priority to EP07742607A priority patent/EP2012555A4/en
Publication of WO2007126054A1 publication Critical patent/WO2007126054A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S3/00—Systems employing more than two channels, e.g. quadraphonic
    • H04S3/02—Systems employing more than two channels, e.g. quadraphonic of the matrix type, i.e. in which input signals are combined algebraically, e.g. after having been phase shifted with respect to each other
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2134—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0272—Voice signal separating

Definitions

  • the present invention relates to a signal separation device, a signal separation method, and a computer or digital method suitable for improving accuracy as much as possible and making convergence as fast as possible in signal separation techniques such as independent component analysis.
  • the present invention relates to a computer-readable information recording medium in which a program realized on a signal processor is recorded, and the program. Background art
  • Patent Document 1 Japanese Unexamined Patent Publication No. 2003-141102
  • Patent Document 1 there is a technology that uses principal component analysis or independent component analysis to separate changes in the amount of a chemical substance into a plurality of components and group the causes of the generation of the chemical substance by the components. Proposed.
  • observation signals are received from m channels, consider an observation signal matrix X of m rows and T columns in which the observed values in the time direction are arranged in the column direction and the rows are arranged in the channel order.
  • Matrix operations c (',') include matrix products, matrix convolution, and nonlinear link functions of various matrices.
  • Matrix function J ( ⁇ , ⁇ , ⁇ ) that accepts an m-by-T matrix, an m-by-n matrix, and an n-by-T matrix and returns a scalar value
  • X-AS maximum absolute value of the element
  • the mean square of each element of, the sum of squares of each element of the matrix (X-AS), etc. can be used. It may also impose constraints such as non-neg ativity, sparseness, and statistical independence as needed.
  • the present invention is for solving the above-described problems, and in a signal separation technique such as independent component analysis, a signal separation device and a signal suitable for improving accuracy as much as possible and making convergence as fast as possible Separation methods, as well as computer or digital signal processing It is an object to provide a computer-readable information recording medium in which a program realized on a sessa is recorded, and the program.
  • the signal separation device includes an m-row T-column observation signal matrix X, an m-row n-column mixing ratio matrix A, an n-row T-column source signal matrix S, and an m-row row.
  • the iterative estimator accepts an m-by-T matrix, an m-by-n matrix, and an n-by-T matrix and returns a scalar value J (', ⁇ , ⁇ )
  • the matrix function J (', ⁇ ) is repeatedly estimated at a predetermined number of iterations.
  • the iterative calculation unit includes L matrices A 1, A 1, A 2,.
  • the result output unit mixes m rows and n columns from the obtained L matrices A, A, A, ..., A.
  • the signal separation device of the present invention may further include a repetition control unit and may be configured as follows.
  • the repetition control unit performs the obtained L matrixes A, A, A, ..., A and n rows and T columns.
  • the result output unit calculates and outputs the matrices A and S when the predetermined convergence condition is satisfied.
  • the result output unit outputs a noise matrix N of m rows and T columns,
  • N X-c (A, S)
  • the matrix N can be further output by further calculation according to
  • the predetermined matrix operation c (′, ′) can be configured to be a matrix product.
  • the predetermined matrix operation c (',') can be configured to be a convolution operation.
  • the predetermined matrix operation c (',') can be configured to be a nonlinear link function.
  • a signal separation method includes an m-row T-column observation signal matrix X, an m-row n-column mixing ratio matrix A, an n-row T-column source signal matrix S, m Row T column noise matrix N for a given matrix operation c (', ⁇ )
  • a signal separation device having an iterative estimation unit, an iterative calculation unit, and a result output unit, and includes an iterative estimation process, an iterative calculation process, and a result output process.
  • the configuration is as follows.
  • the iterative estimation unit accepts an m-row T-column matrix, an m-by-n matrix, and an n-by-T matrix, and returns a scalar value J (', ⁇ Given an m-by-T matrix, an m-by-n matrix, and an n-by-T matrix, the matrix function J (', ⁇ ).
  • the combination of an m-by-n matrix and an n-by-T matrix of arguments that are extremal values of the scalar values returned by ⁇ , ⁇ ) are repeatedly estimated at a predetermined number of iterations.
  • the iterative calculation unit includes L matrices A 1, A 1, A 2,.
  • the result output unit includes the obtained L matrices A 1, A 2, A 3,.
  • the signal separation device further includes a repetition control unit, further includes a repetition control step, and can be configured as follows.
  • the repetition control unit performs the obtained L matrices A 1, A 2, A 3,
  • the iteration is repeated until A and the n-row T-column matrix S satisfy the predetermined convergence condition.
  • a program according to another aspect of the present invention is configured to cause a computer or a digital signal processor to function as each unit of the signal separation device or to execute the signal separation method.
  • the program of the present invention can be recorded on a computer-readable information storage medium such as a compact disk, flexible disk, hard disk, magneto-optical disk, digital video disk, magnetic tape, and semiconductor memory.
  • a computer-readable information storage medium such as a compact disk, flexible disk, hard disk, magneto-optical disk, digital video disk, magnetic tape, and semiconductor memory.
  • the above program can be distributed and sold via a computer communication network independently of a computer or a digital signal processor on which the program is executed.
  • the information storage medium can be distributed and sold independently of a computer or a digital signal processor.
  • a signal separation apparatus a signal separation method, and a signal separation technique suitable for improving accuracy as much as possible and making convergence as fast as possible for signal separation techniques such as independent component analysis.
  • a computer-readable information recording medium in which a program for realizing these on a computer or a digital signal processor is recorded, and the program.
  • FIG. 1 is a schematic diagram showing a schematic configuration of a transmission device and a signal separation device according to one embodiment of the present invention.
  • FIG. 2 is a schematic diagram showing a configuration of a signal separation device according to the present embodiment.
  • FIG. 3 is a flowchart showing a control flow of signal separation processing executed by the signal separation device.
  • FIG. 4 is a graph showing the waveform of the source signal in the first experiment.
  • FIG. 5 is a graph showing the waveform of the observation signal in the first experiment.
  • FIG. 8 is a graph showing the waveform of the source signal in the second experiment.
  • FIG. 9 is a graph showing the waveform of the observation signal in the second experiment.
  • FIG. 11 is a graph showing the SIR value (Singnal to Interference Ratio) in the second experiment.
  • FIG. 12 Waveform obtained in the second experiment by the conventional method called NMF Lee-Seung algorithm.
  • FIG. 13 is a diagram showing an image serving as a source signal in a third experiment.
  • FIG. 14 is a diagram showing an image corresponding to an observation signal in a third experiment.
  • FIG. 17 is a diagram showing an image of a source signal obtained by a conventional technique called NMF Lee-Seung algorithm in the third experiment.
  • FIG. 18 is a graph showing the SIR value of the signal in the method of the present embodiment in the third experiment.
  • component analysis method is changed to the principal component analysis (component analysis by sphering is performed by the same method as this). And can be replaced by non-negative matrix factorization, and these embodiments are also within the scope of the present invention.
  • the following matrix X is used as an input for signal separation.
  • the cost function J (X, A, S) is minimum or minimum (depending on the type of cost function, the maximum or maximum is In the following, it shall be such that the cost function is minimized in order to facilitate understanding. That is, the following holds.
  • the above-mentioned cost function uses the absolute value of the element with the maximum absolute value (maximum absolute value of the element) of the matrix (X-AS) and the matrix (X-AS) Root mean square of elements A power of the sum of squares of each element of the matrix (X-AS) Amari's a divergence, Kullback Leibl er divergence ⁇ Trobenius norm, Jenssen Shannon divergence, etc.
  • an iterative method (gradient method, conjugate gradient method, Newton method, etc.) is employed. In general, this iteration is repeated until A and S satisfy a predetermined iteration termination condition. For example, the repetition end condition is satisfied when the specific power of the sum of squares of the difference between the elements before and after the iteration of the matrix with respect to the sum of squares of the elements of A and S becomes smaller than a predetermined value ⁇ .
  • the calculation may be performed by imposing constraint conditions such as non-negativeness, looseness, spatiotemporal decorrelation, and smoothness independence based on the known properties of the signal to be solved.
  • such a function of performing independent component analysis is regarded as one calculation module called an iterative estimation unit.
  • an iterative estimation unit For the processing executed by the iterative estimator, various independent component analysis techniques can be applied in addition to the existing techniques.
  • FIG. 1 is an explanatory diagram showing the state of input and output to the iterative estimation unit. This will be described below with reference to this figure.
  • the iterative estimation unit 101 is realized by a general computer CPU (Central Processing Unit) working together with a storage device such as a RAM (Random Access Memory).
  • a general computer CPU Central Processing Unit
  • RAM Random Access Memory
  • matrices X, A, and S are given as inputs.
  • the matrices X, A, and S are stored in a RAM or the like.
  • the library (program) function that implements the iterative estimation unit 101 is stored in the memory in which the matrices X, A, and S are secured. The address is passed as an argument.
  • the matrix X can be considered as a constant matrix that defines conditions in the iterative method when considered broadly as the force corresponding to the observed signal matrix as described above.
  • the matrixes A and S given here are a mixture ratio matrix and a source signal matrix, but can be considered as initial values that are approximate values of true values in the iterative method. Therefore, it is stored in RAM etc.
  • the matrices A and S to be initialized need to be initialized with some value.
  • each element of matrix A is given a random number or a constant (such as 0, 1 or — 1) so that the element distribution is sparse, and each element of matrix S is a random number.
  • a constant such as 0, 1 or — 1
  • constants (0, 1, — 1 etc.) or various techniques are considered (“Initialization” in the figure). Further, as will be described later, this is preferable in many cases where the result output by the iterative estimation unit 101 in the previous process may be set as the initial value.
  • the cost function and the iteration end condition are also given to the iteration estimation unit 101. These may be directly specified in the code (program) for realizing the iterative inference unit 101, or may be by specifying the address of a library function.
  • the iterative estimation unit 101 receives various pieces of information as described above, repeats the iteration until the iteration end condition is satisfied, and further approximates the matrices A and S given as inputs to true values.
  • Matrix A ', S' is output.
  • the output of A ', S' means that the iteration end condition has been satisfied.
  • the following conditions can be considered as the iteration termination condition.
  • the cost function J (X, A, S) is treated as a combination of two cost functions K (X, A, S), H (X, A, S).
  • K X, A, S
  • H X, A, S
  • 1 (, and, may select the same cost function (that is, K, H, and J are the same cost function), and select different cost functions as appropriate. It's okay.
  • the iterative estimation unit 101 as described above is used at least L (L ⁇ 2) times by a subroutine call, but the arguments given at that time are different. Therefore, in order to facilitate understanding, in the following description, it is assumed that L iteration estimation units 101 are individually prepared.
  • FIG. 2 is a schematic diagram showing a configuration of the signal separation device according to the present embodiment. This will be described below with reference to this figure.
  • the signal separation device 100 includes L iteration estimation units 101, an iteration calculation unit 102, and a result output unit 103.
  • a storage unit 105 is provided.
  • L corresponds to “the number of steps” and “the number of layers”.
  • the storage unit 105 an area for storing the following information is secured.
  • the area in 05 is indicated by the same symbol.
  • A is an n-by-n matrix.
  • FIG. 3 is a flowchart showing a control flow of signal separation processing executed by the signal separation device.
  • FIG. 3 is a flowchart showing a control flow of signal separation processing executed by the signal separation device.
  • the matrix X (that is, the matrix S) is entered.
  • Column S can be identified with matrix S.
  • step S302 and step S303 For initialization of values in step S302 and step S303, various techniques can be applied in addition to the initialization using the initial values described above. In the present embodiment, these initialization processes are realized by the CPU cooperating with the RAM or the like.
  • step S304 the repetition control unit 104 does not satisfy the convergence condition! /, Between (step S308), The processing from step S304 to step S308 is repeated (step S304).
  • the convergence condition various forces as well as the iteration end condition can be adopted.When the change of the matrix S, ⁇ , S and the matrix A, ⁇ , A in the iteration becomes extremely small, The convergence condition is
  • the convergence condition determined at the end of the iteration may always be true. In this case, the repetition is performed only once, and in reality, only sequential processing is performed.
  • the repetition control unit 104 is realized by the CPU cooperating with the RAM or the like.
  • the iteration calculation unit 102 repeats the processing of the following steps S305 to S307 while incrementing the counter variable i by 1 until the initial value 1 reaches the end value L. (Step S305).
  • the iterative calculation unit 102 is realized by the CPU cooperating with the RAM or the like.
  • the i-th iterative estimation unit 101 is given matrices S 1, A 2 and S as inputs, and the above
  • step S306 The calculation process of the iterative method is performed (step S306).
  • step S306 In the calculation process of the iterative method in step S306, as described above, the iteration is repeated until the iteration end condition is satisfied. Which iteration end condition is to be selected is appropriately selected according to the application field and purpose. can do.
  • the matrices that approach the "true value" by the iterative method in the process of step S306 are A and S.
  • the current estimation result in the iterative estimation unit 101 is stored in the storage areas A and S.
  • step S308 the process returns to step S304 and repeats (step S308).
  • the result output unit 103 sets the mixture ratio matrix A as
  • the noise matrix N is obtained by multiplying the matrix as
  • N X-AS
  • step 309 corresponds to obtaining the matrix A as described above.
  • the observation signal matrix X and the source signal matrix S are two-dimensional matrices as described above, but the mixing ratio matrix A is a three-dimensional matrix.
  • the elements of the three-dimensional matrix are expressed in the form of [ ⁇ , ⁇ , ⁇ ] as above, and the value of the range where the elements are not defined is set to 0 according to the convention of tensor operation. It is assumed that the sum of all the scans of the range that the subscript of ⁇ can take as an index in the tensor is taken.
  • FIG. 4 is a graph showing the waveform of the source signal
  • FIG. 5 is a graph showing the waveform of the observation signal.
  • FIG. 8 is a graph showing the waveform of the source signal
  • FIG. 9 is a graph showing the waveform of the observation signal.
  • Figure 12 shows the waveform obtained by the conventional method called NMF Lee-Seung algorithm.
  • the SIR value is less than 8 dB, which indicates the high performance of the method of the present embodiment.
  • the third experiment is an application in image processing.
  • FIG. 13 is a diagram showing an image that is a source signal
  • FIG. 14 is a diagram showing an image corresponding to an observation signal
  • Fig. 17 shows the conventional NMF Lee-Seung algorithm.
  • FIG. 18 is a graph showing the value of the SIR value of the signal in the method of the present embodiment.
  • a signal separation technique such as independent component analysis, a signal separation device, a signal separation method, and a method suitable for improving accuracy as much as possible and making convergence as fast as possible.
  • a computer-readable information recording medium in which a program for realizing these on a computer or a digital signal processor is recorded, and the program.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computational Linguistics (AREA)
  • Quality & Reliability (AREA)
  • Artificial Intelligence (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Human Computer Interaction (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Multimedia (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

　信号分離装置（１００）において、繰返計算部（１０２）は、反復推測部（１０１）に、観測信号行列Xを反復法により独立成分分析させ、その結果得られた源信号行列をさらに反復法により独立成分分析させることを繰り返させ、結果出力部（１０３）は、繰返しの過程で得られた各混合比行列の積を、観測信号行列Xに対する混合比行列Aとして出力し、繰り返しの最後に得られた源信号行列を観測信号行列Xに対する源信号行列Sとして出力する。繰返制御部（１０４）は、繰返計算部（１０２）に、各混合比行列や各源信号行列が収束条件を満たすまで、計算制御を繰り返させる。反復推測部（１０１）における反復法による反復回数は、固定回数でも良いし、収束がするまで反復を行うこととしても良い。

Description

明 細 書
信号分離装置、信号分離方法、情報記録媒体、ならびに、プログラム 技術分野
[0001] 本発明は、独立成分分析等の信号分離技術において、精度をできるだけ向上し、 収束をできるだけ速くするのに好適な、信号分離装置、信号分離方法、ならびに、こ れらをコンピュータもしくはディジタル信号プロセッサ上にて実現するプログラムを記 録したコンピュータ読取可能な情報記録媒体、ならびに、当該プログラムに関する。 背景技術
[0002] 従来から、主成分分析 (スフイアリングによる成分分析を含む。）、独立成分分析、 非負行列因子化（Non-negative Matfix Factorization； NMF)による成分分析など、 あるデータ列が与えられたときに、当該データ列を複数の成分 (主成分、独立成分、 因子等）に分離する手法が提案されている。たとえば、以下のような文献で、このよう な技術が開示されている。
特許文献 1：特開 2003— 141102号公報
[0003] 特許文献 1では、主成分分析や独立成分分析を用いて、化学物質の量の変化を 複数の成分に分離し、当該成分により、当該化学物質の生成の原因をグループ化す る技術が提案されている。
[0004] 一般に、独立成分分析においては、
(1) m個のチャンネルから観測信号を受け付けたときに、その時間方向の観測値を 列方向に並べた行を、チャンネルの順序に並べた m行 T列の観測信号行列 X を考える。
[0005] そして、観測信号行列 Xを、
(2) n個のチャンネルの源信号の時間方向の源信号値を列方向に並べた行を、チ ヤンネルの順序に並べた n行 T列の源信号行列 Sと、
(3)源信号の各チャンネル力 観測信号の各チャンネルまでの経路の様子を表す m行 n列の混合比行列 Aと、
(4) m行 T列の雑音行列 Nと、 に、分離する。
[0006] その際に、所定の行列演算 c(' , ')について
(5) X = c(A,S) + N
を満たすように分離をする。行列演算 c(' , ')としては、行列の積、行列のコンボリュー シヨンのほか、各種の行列の非線型リンク関数が用いられる。
[0007] また、この際に、勾配法、共役勾配法、ニュートン法などの反復法を用いるのが一 般的であるが、
(6) m行 T列の行列と m行 n列の行列と n行 T列の行列とを受け付けてスカラー値を返 す行列関数 J( · , · , ·)
をコスト関数として採用する。
[0008] 具体的には、 J(X,A,S)に対して、 Xを固定して、 A, Sを変化させたときに、 J(X,A,S)の 値が最小 (極小。一般には、極大もしくは極小、すなわち、「極値」。）となるような A, S の組合せを計算する。
[0009] コスト関数としては、たとえば、 J(X,A,S)として、行列 (X - AS)の絶対値最大の要素の 絶対値 (要素の最大絶対値)や、行列 (X - AS)の各要素の自乗平均、行列 (X - AS) の各要素の総自乗和等を用いることができる。また、必要に応じて、非負性 (non-neg ativity)、疎性（sparseness)、統計的独立性（statistical independence)などの制約を 課すことちある。
発明の開示
発明が解決しょうとする課題
[0010] し力しながら、単純に反復法を適用するのみでは、源信号同士の強弱の差が大き い場合や観測信号同士の差が小さい場合 (観測信号同士が似ている場合)には、収 束に時間がかかり、しかも信号の分離性能が落ちるという問題があった。
[0011] したがって、このような場合であっても、高速に収束し、信号の分離性能を向上させ るような信号分離技術に対する要望は強 、。
[0012] 本発明は、上記の課題を解決するためのもので、独立成分分析等の信号分離技術 において、精度をできるだけ向上し、収束をできるだけ速くするのに好適な、信号分 離装置、信号分離方法、ならびに、これらをコンピュータもしくはディジタル信号プロ セッサ上にて実現するプログラムを記録したコンピュータ読取可能な情報記録媒体、 ならびに、当該プログラムを提供することを目的とする。
課題を解決するための手段
[0013] 本発明の第 1の観点に係る信号分離装置は、 m行 T列の観測信号行列 Xを、 m行 n 列の混合比行列 Aと n行 T列の源信号行列 Sと m行 T列の雑音行列 Nとであって、所定 の行列演算 c(' , ')について
X = c(A,S) + N
を満たす行列 A, Sに分離し、反復推測部、繰返計算部、結果出力部を備え、以下の ように構成する。
[0014] すなわち、反復推測部は、 m行 T列の行列と m行 n列の行列と n行 T列の行列とを受 け付けてスカラー値を返す行列関数 J(' , · , ·)に対して、 m行 T列の行列と m行 n列の行 列と n行 T列の行列とが与えられると、当該 m行 T列の行列を固定したときに当該行列 関数 J(' , · , ·)が返すスカラー値を極値ィ匕する引き数の m行 n列の行列と n行 T列の行列 との組合せを、所定の繰返し回数で反復推測する。
[0015] 一方、繰返計算部は、 L個の行列 A , A , A , · ··, Aであって、
1 2 3 L
c(A ,c(A ,c(A ,c(- --,c(A , A )· ··))))
1 2 3 L-l L
が m行 n列の行列となる L個の行列を繰返しにより計算し、
S = X
0
として、 i (= 1, 2, 3, · ··, L)回目の繰返しにおいて、反復推測部に、 i回目の繰返しに おける行列関数 J (· , · , ·)であって、 m行 T列の行列と m行 n列の行列と n行 T列の行列と を受け付けてスカラー値を返す行列関数 J (· , · , に対して、 m行 T列の行列 S と m行 n
i i-1 列の所定の初期値行列と n行 T列の所定の初期値行列とを与えて、反復推測させ、 得られた m行 η列の行列を Αとし、 n行 T列の行列を Sとする。
[0016] さらに、結果出力部は、求められた L個の行列 A , A , A , · ··, Aから、 m行 n列の混
1 2 3 L
合比行列 Aを
A = c(A ,c(A ,c(A ,c(- --,c(A ,A )· ··))))
1 2 3 L-l L
により計算し、求められた n行 T列の行列 Sを、 n行 T列の源信号行列 Sとして、行列 A, し
Sを出力する。 [0017] また、本発明の信号分離装置は、繰返制御部をさらに備え、以下のように構成する ことができる。
[0018] すなわち、繰返制御部は、求められた L個の行列 A , A , A , · ··, Aと n行 T列の行
1 2 3 L
列 Sとが所定の収束条件を満たすまで、繰返計算部に当該繰返しをやり直させ、 i回 し
目の繰返しにおける m行 n列の所定の初期値行列および n行 T列の所定の初期値行 列として、やり直しの前回における Aおよびやり直しの前回における Sを与えさせる。
[0019] 一方、結果出力部は、当該所定の収束条件が満たされると、行列 A, Sを計算して 出力する。
[0020] また、本発明の信号分離装置において、結果出力部は、 m行 T列の雑音行列 Nを、
N = X - c(A,S)
によりさらに計算して、行列 Nをさらに出力するように構成することができる。
[0021] また、本発明の信号分離装置において、当該所定の行列演算 c(' , ')は、行列の積 であるように構成することができる。
[0022] また、本発明の信号分離装置において、当該所定の行列演算 c(' , ')は、コンボリュ ーシヨン演算であるように構成することができる。
[0023] また、本発明の信号分離装置において、当該所定の行列演算 c(' , ')は、非線型リン ク関数であるように構成することができる。
[0024] また、本発明の信号分離装置において、反復推測部における反復回数は、 i (= 1,
2, 3, · ··, L)回目のいずれの繰返しにおいても 1回であるように構成することができる
[0025] また、本発明の信号分離装置において、反復推測部に対する i (= 1, 2, 3, · ··, L)回 目の行列関数 J (· , · , ·)は、 、ずれも等 U、ように構成することができる。
[0026] 本発明のその他の観点に係る信号分離方法は、 m行 T列の観測信号行列 Xを、 m行 n列の混合比行列 Aと、 n行 T列の源信号行列 Sと、 m行 T列の雑音行列 Nであって、所 定の行列演算 c(' , ·)について
X = c(A,S) + N
を満たす行列 A, S, Nに分離し、反復推測部、繰返計算部、結果出力部を有する信 号分離装置にて実行され、反復推測工程、繰返計算工程、結果出力工程を備え、 以下のように構成する。
[0027] すなわち、反復推測工程では、反復推測部が、 m行 T列の行列と m行 n列の行列と n 行 T列の行列とを受け付けてスカラー値を返す行列関数 J(' , · , に対して、 m行 T列の 行列と m行 n列の行列と n行 T列の行列とが与えられると、当該 m行 T列の行列を固定 したときに当該行列関数 J(' , · , ·)が返すスカラー値を極値ィ匕する引き数の m行 n列の 行列と n行 T列の行列との組合せを、所定の繰返し回数で反復推測する。
[0028] 一方、繰返計算工程では、繰返計算部が、 L個の行列 A , A , A , · ··, Aであって、
1 2 3 し c(A ,c(A ,c(A ,c("*,c(A ，A )· ··》》
1 2 3 L-l し
が m行 n列の行列となる L個の行列を繰返しにより計算し、
S = X
0
として、 i (= 1, 2, 3, · ··, L)回目の繰返しにおいて、反復推測部に、 i回目の繰返しに おける行列関数 J (· , · , ·)であって、 m行 T列の行列と m行 n列の行列と n行 T列の行列と を受け付けてスカラー値を返す行列関数 J
i (· , · , に対して、 m行 T列の行列 S
i-1と m行 n 列の所定の初期値行列と n行 T列の所定の初期値行列とを与えて、反復推測させ、 得られた m行 η列の行列を Αとし、 n行 T列の行列を Sとする。
[0029] さらに、結果出力工程では、結果出力部が、求められた L個の行列 A , A , A , · ··,
1 2 3
Aから、 m行 n列の混合比行列 Aを
し
A = c(A ，c(A ，c(A ,c("*,c(A ，A )· ··》))
1 2 3 L-l し
により計算し、求められた n行 T列の行列 Sを、 n行 T列の源信号行列 Sとして、行列 A, し
Sを出力する。
[0030] また、本発明の信号分離方法において、当該信号分離装置は、繰返制御部をさら に有し、繰返制御工程をさらに備え、以下のように構成することができる。
[0031] すなわち、繰返制御工程では、繰返制御部が、求められた L個の行列 A , A , A ,
1 2 3
· ··, Aと n行 T列の行列 Sとが所定の収束条件を満たすまで、繰返計算部に当該繰 し し
返しをやり直させ、 i回目の繰返しにおける m行 n列の所定の初期値行列および n行 T 列の所定の初期値行列として、やり直しの前回におけるへおよびやり直しの前回にお ける Sを与えさせる。
[0032] 一方、結果出力工程では、当該所定の収束条件が満たされると、行列 A, Sを計算 して出力する。
[0033] 本発明のその他の観点に係るプログラムは、コンピュータもしくはディジタル信号プ 口セッサを上記の信号分離装置の各部として機能させ、もしくは、上記の信号分離方 法を実行させるように構成する。
[0034] また、本発明のプログラムは、コンパクトディスク、フレキシブルディスク、ハードディ スク、光磁気ディスク、ディジタルビデオディスク、磁気テープ、半導体メモリ等のコン ピュータ読取可能な情報記憶媒体に記録することができる。
[0035] 上記プログラムは、プログラムが実行されるコンピュータやディジタル信号プロセッ サとは独立して、コンピュータ通信網を介して配布 '販売することができる。また、上記 情報記憶媒体は、コンピュータやディジタル信号プロセッサとは独立して配布'販売 することができる。
発明の効果
[0036] 本発明によれば、独立成分分析等の信号分離技術にお!、て、精度をできるだけ向 上し、収束をできるだけ速くするのに好適な、信号分離装置、信号分離方法、ならび に、これらをコンピュータもしくはディジタル信号プロセッサ上にて実現するプログラム を記録したコンピュータ読取可能な情報記録媒体、ならびに、当該プログラムを提供 することができる。
図面の簡単な説明
[0037] [図 1]本発明の実施形態の一つに係る送信装置および信号分離装置の概要構成を 示す模式図である。
[図 2]本実施形態に係る信号分離装置の構成を示す模式図である。
[図 3]本信号分離装置にて実行される信号分離処理の制御の流れを示すフローチヤ ートである。
[図 4]第 1の実験における源信号の波形の様子を示すグラフである。
[図 5]第 1の実験における観測信号の波形の様子を示すグラフである。
[図 6]第 1の実験において、 L = 10として共役勾配法を用いて本実施形態の手法によ り 4つの信号に独立成分分析を行った場合の波形のグラフである。
[図 7]第 1の実験において、従来の手法により 4つ信号に独立成分分析を行った場合 の波形のグラフである。
[図 8]第 2の実験における源信号の波形の様子を示すグラフである。
[図 9]第 2の実験における観測信号の波形の様子を示すグラフである。
[図 10]第 2の実験において L = 5として甘利の a - divergenceNMFアルゴリズムを反 復推測部にぉ 、て用いた本実施形態の手法により得られた波形のグラフである。
[図 11]第 2の実験における SIR値（Singnal to Interference Ratio)を示すグラフである
[図 12]第 2の実験において NMF Lee-Seungアルゴリズムという従来の手法により得ら れた波形である。
[図 13]第 3の実験における源信号となる画像を示す図である。
[図 14]第 3の実験における観測信号に相当する画像を示す図である。
[図 15]第 3の実験において L = 5として本実施形態の手法を採用した場合に得られる 源信号の様子を示す図である。
[図 16]第 3の実験において L = 2として本実施形態の手法を採用した場合に得られる 源信号の様子を示す図である。
[図 17]第 3の実験において、 NMF Lee-Seungアルゴリズムという従来の手法により得 られた源信号の画像を示す図である。
[図 18]第 3の実験における本実施形態の手法における信号の SIR値の値を示すダラ フである。
符号の説明
[0038] 100 信号分離装置
101 反復推測部
102 繰返計算部
103 結果出力部
104 繰返制御部
105 記憶部
発明を実施するための最良の形態
[0039] 以下に本発明の実施形態を説明する。なお、以下に説明する実施形態は説明の ためのものであり、本発明の範囲を制限するものではない。したがって、当業者であ ればこれらの各要素もしくは全要素をこれと均等なものに置換した実施形態を採用 することが可能である力 これらの実施形態も本発明の範囲に含まれる。
[0040] なお、以下の説明では、成分分析の手法として独立成分分析を例としてとりあげて 説明するが、これと同様の手法によって、成分分析の手法を主成分分析 (スフィアリン グによる成分分析を含む。）非負行列因子化に置き換えることができ、これらの実施 形態も本発明の範囲に含まれる。
[0041] また、以下の説明では、理解を容易にするため、行列の演算 c(' , ')として、行列の 積を例にあげて説明する力 行列のコンボリューシヨンや各種の行列の非線型リンク 関数を利用することも可能であり、これらを採用した場合も本発明の範囲に含まれる。 実施例 1
[0042] まず、本実施例における計算処理の理解を容易にするため、各種記号について、 説明する。
[0043] 本実施例では、以下の行列 Xを、信号分離の入力として用いる。
(1) m行 T列の観測信号行列 X
[0044] そして、以下の行列 A, S, Nを信号分離の結果として得る。
(2) m行 n列の混合比行列 A
(3) n行 T列の源信号行列 S
(4) m行 T列の雑音行列 N
[0045] ただし、これらの行列の間には、以下の関係が成立しなければならない。
(5) X = AS + N
[0046] また、 Xを固定して A, Sを微小に変化させたとしたときに、コスト関数 J(X,A,S)最小あ るいは極小 (コスト関数の種類によっては、最大あるいは極大を含む「極値」。以下で は、理解を容易にするため、コスト関数を極小にするものを求めるものとする。）となる ものでなければならない。すなわち、以下が成立する。
(6) (A,S) = argmin J(X,A,S)
(A,S)
[0047] 一般的な独立成分分析の手法では、上記のコスト関数として行列 (X - AS)の絶対 値最大の要素の絶対値 (要素の最大絶対値)や、行列 (X - AS)の各要素の自乗平均 、行列 (X - AS)の各要素の総自乗和等のほ力 甘利の a divergence, Kullback Leibl er divergence ^ト robenius norm、 Jenssen Shannon divergence等 用 ヽること力 eさ
[0048] また、上記の行列 A, Sを求める計算手法としては、反復法 (勾配法、共役勾配法、 ニュートン法等)を採用する。通常、この反復は、 A, Sが所定の反復終了条件を満た すようになるまで繰り返されるのが一般的である。たとえば、 A, Sの要素の総自乗和 に対する当該行列の反復の前後の要素の差の総自乗和の比力 所定の値 εよりも 小さくなつた場合に反復終了条件が満たされる等である。
[0049] もっとも、粗い分析で十分な場合等には、反復回数を数回に固定することもありうる 。後者の場合は、反復終了条件として「反復回数〇〇回」を指定したことと等価である
[0050] このほか、解くべき信号の既知の性質に基づいて、非負性、疎性、時空間無相関 性、 smoothness独立性などの束縛条件を課して、計算を行うこともある。
[0051] 本実施形態では、このような独立成分分析を行う機能を、反復推測部という 1つの 計算モジュールとしてとらえるものとする。反復推測部が実行する処理については、 既存の技術の他、各種の独立成分分析の技術を適用することができる。
[0052] 図 1は、反復推測部に対する入力と出力の様子を示す説明図である。以下、本図 を参照して説明する。
[0053] 反復推測部 101は、一般的なコンピュータの CPU (Central Processing Unit)が RA M (Random Access Memory)などの記憶装置と共働して実現される。
[0054] 反復推測部 101に対しては、入力として行列 X, A, Sが与えられる。一般には、行列 X, A, Sは RAM等に記憶されており、本実施形態では、反復推測部 101を実現する ライブラリ（プログラム）関数に、行列 X, A, Sが確保されているメモリ内のアドレスが引 き数として渡される構成をとる。
[0055] ここで、行列 Xは、上記のように観測信号行列に相当するものである力 広く考える と、反復法における条件を定める定数行列と考えることができる。一方、ここで与えら れる行列 A, Sは、混合比行列と源信号行列であるが、広く考えると、反復法において 、真の値の近似値となる初期値と考えることができる。したがって、 RAM等に記憶さ れる行列 A, Sは、何らかの値で初期化される必要がある。
[0056] 何ら情報がない場合には、行列 Aの各要素には、要素分布が疎であるように乱数や 定数 (0や 1や— 1等)を与え、行列 Sの各要素には乱数や定数 (0, 1 , — 1等)を与え たり、種々の技術が考えられる（図中の「初期化」）。また、後述するように、反復推測 部 101が前回の処理において出力した結果を初期値とすることとしても良ぐ多くの 場合、これが好ましい。
[0057] このほか、コスト関数と反復終了条件も、反復推測部 101に対して与えられる。これ らは、反復推測部 101を実現するためのコード (プログラム）に直接指定されることも あるし、ライブラリ関数のアドレスを指定することによる場合もありうる。
[0058] さて、反復推測部 101は、上記のような諸情報を受け付けて、反復終了条件が満た されるまで反復を繰り返し、入力として与えられた行列 A, Sを、さらに真の値に近付け た行列 A', S'を出力する。 A', S'が出力された、ということは、反復終了条件が満たさ れた、ということである。ここで、反復終了条件としては、以下のようなものが考えられ る。
( 1)反復回数が 1回になったとき
(2)反復回数が所定の η回になったとき (η≥2)
(3) (∑ m ∑ n |A'[i,k]-A[i,k]|2 + ∑ n ∑ T |S'[i,k]-S[i,k]|2) ≤ ε
i=l k=l k=l t=l
(4) (max m max n |A，[i，k]— A[i，k]|2 + max n max T |S'[i,k]-S[i,k]|2) ≤ ε i=l k=l k=l t=l
(5)∑ m ∑ n |A'[i,k]-A[i,k]|2≤ ε
i=l k=l
(6) max m max n |A'[i,k]-A[i,k]|2 ≤ ε
i=l k=l
(7)∑ n ∑ T |S'[i,k]— S[i,k]|2 ≤ ε
k=l t=l
(8) max n max T |S'[i,k]-S[i,k]|2 ≤ ε
k=l t=l
[0059] ただし、出力方法としては、入力として与えられた行列 A, Sの RAM等内における領 域に、 A' , S'の内容を上書き（図中の「上書き」。）することとして、出力結果としても良 い。本実施形態では、主にそのような態様を採用する。
[0060] さて、上記の事項に合わせて説明すると、反復法においては、 Xを固定して J(X,A,S) 力 S小さくなるような方向に A, Sを移動させた結果を A', S'とするのである力 2つの行 列 A, Sをまとめて勾配法等に適用することはのは難しい場合もある。 [0061] そこで、コスト関数 J(X,A,S)を、 2つのコスト関数 K(X,A,S), H(X,A,S)の組み合わせと して取り扱う。ここで、 1 ( ， と ， は、同じコスト関数を選択しても良いし (すな わち、 K, H, Jの 3つは同じコスト関数である。）、適宜異なるコスト関数を選択しても良 い。
[0062] そして、
(1)まず、 Χ,Αを固定して、 K(X,A,S)を小さくする方向に Sを移動させた結果を S'とし
(2)次に、 X,S'を固定して、 H(X,A,S')を小さくする方向に Aを移動させた結果を A'と する
ことを、反復の一回分と考えるのが典型的である。
[0063] あるいはその逆に、
(1)まず、 X,Sを固定して、 H(X,A,S)を小さくする方向に Aを移動させた結果を A'とし
(2)次に、 Χ,Α'を固定して、 K(X,A',S)を小さくする方向に Sを移動させた結果を S'と する
ことを、反復の一回分としても良い。
[0064] さて、本実施形態においては、上記のような反復推測部 101を、サブルーチン呼び 出しによって少なくとも L (L≥2)回用いるが、その際に与える引き数等が異なる。そこ で、理解を容易にするため、以下では、個別に反復推測部 101を L個用意したものと して説明する。
[0065] 図 2は、本実施形態に係る信号分離装置の構成を示す模式図である。以下、本図 を参照して説明する。
[0066] 信号分離装置 100は、 L個の反復推測部 101、繰返計算部 102、結果出力部 103
、繰返制御部 104のほか、記憶部 105を備える。
[0067] 本実施形態では、反復推測の多段構成によって独立成分分析を行うため、 Lは、「 段 (step)数」「層（layer)数」に相当するものである。
[0068] 記憶部 105には、以下の情報を記憶するための領域が確保される。
(l) m行 T列の観測信号行列 X。これは、信号分離の象となる観測値である。 (2) L個の行列 A , · ··, A。これは、一時的な計算のために用いられる。
1 し
(3) L個の n行 T列の行列 S ， …， S。これは、一時的な計算のために用いられる。
1 し
(4) m行 n列の混合比行列 A。これは、信号分離の結果の一つである。
(5) m行 T列の雑音行列 N。これは、信号分離の結果の一つである。
[0069] なお、信号分離の結果の一つである n行 T列の源信号行列 Sは、 Sと同一視すること し
ができる。
[0070] また、漸化式による表現の理解を容易にするため、仮想的に、行列 Xを Sと同一視
0
することとする。
[0071] さらに、以下では、混乱のな!、限り、行列名と行列の要素の値を記憶する記憶部 1
05内での領域とを同じ記号で表記するものとする。
[0072] なお、 L個の行列 A , · ··, Aは、所定の演算 c(' , ')
1 し
c(A ,c(A ,c(A ,c("*,c(A ，A )· ··》》
1 2 3 L-l し
が m行 n列の行列となるようなものであり、本実施形態では、 c(' , ')として、行列の積を 用いている。
[0073] したがって、本実施形態では、 L個の行列の行数と列数を定める数列
n = m, η , η , · ··, n , n = n
0 1 2 L-l L
がー意に定まり、 Aは、 n 行 n列の行列であることになる。
i i-1 i
[0074] 図 3は、本信号分離装置にて実行される信号分離処理の制御の流れを示すフロー チャートである。以下、本図を参照して説明する。
[0075] 信号分離装置 100における処理が開始されると、行列 X (すなわち、行列 S )に、入
0 力となる観測信号の値を設定する (ステップ S301)。
[0076] 次に、行列 S , S S , Sを適切な値で初期化する (ステップ S 302)。ここで、行
1 2,…， L-l L
列 Sは行列 Sと同一視することができる。
し
[0077] さらに、行列 A , A , · ··, A , Aを適切な値で初期化する（ステップ S303)。
1 2 L-l L
[0078] ステップ S302およびステップ S303における値の初期化は、上述した初期値による 初期化のほか、種々の技術を適用することができる。本実施形態においては、これら の初期化処理は、 CPUが RAM等と共働することによって実現される。
[0079] そして、繰返制御部 104は、収束条件が満たされな!/、間 (ステップ S308)、以下の ステップ S 304〜ステップ S 308の処理を繰り返す (ステップ S 304)。収束条件につい ては、反復終了条件同様種々のものを採用することができる力 当該繰返しにおける 行列 S , · ··, Sおよび行列 A , · ··, Aの変化が極めて小さくなつたときに、収束条件が
1 し 1 し
満たされたものと考えるのが典型的である。
[0080] もっとも、適用分野によっては、繰返しの終端で判断される収束条件は常に真であ ることとしても良い。この場合、当該繰返しは 1回だけ行われ、実際には単なる逐次処 理を行うのみとなる。
[0081] したがって、繰返制御部 104は、 CPUが RAM等と共働することによって、実現され る。
[0082] さて、当該繰返しの中で、繰返計算部 102は、カウンタ変数 iを、初期値 1から終了 値 Lになるまで、 1ずつ増やしながら、以下のステップ S305〜ステップ S307の処理 を繰り返す (ステップ S305)。
[0083] したがって、繰返計算部 102は、 CPUが RAM等と共働することによって、実現され る。
[0084] すなわち、 i番目の反復推測部 101に、入力として、行列 S , A , Sを与えて、上記
i-1 i i
のような反復法の計算処理を行わせる (ステップ S306)。
[0085] ステップ S306における反復法の計算処理では、上述の通り、反復終了条件が満た されるまで反復を繰り返すが、いずれの反復終了条件を選択するかは、適用分野や 目的に応じて適宜選択することができる。
[0086] また、潘目の反復推測部 101に与えられるコスト関数 Jおよび反復終了条件 (0は、 i
= 1, · ··, Lのすべてにおいて同じものとしても良いし、異なるものとしても良い。
[0087] ステップ S306の処理における反復法により「真の値」に近付く行列は、 Aおよび Sで ある。そして、上記のように、記憶領域 Aおよび Sに、反復推測部 101における今回の 推測結果が格納されることになる。
[0088] このようにして、反復推測部 101における反復推測を i = 1, · ··, Lの間繰り返すと (ス テツプ S307)、ステップ S302、ステップ S303【こお!/、て行歹 US , · ··, Sおよび行歹 [JA
1 L 1
, · ··, Aに設定された初期値よりも、精度の高い値が、同じ領域に格納されることにな し
る。 [0089] さて、この値を吟味することによって、収束条件が満たされたた力否かを判定し、収 束条件が満たされていなければ、ステップ S304に戻って繰り返しを続ける (ステップ S308)。
[0090] 一方、収束条件が満たされていれば、結果出力部 103は、混合比行列 Aを、
Α = Α Α · ··Α A
1 2 L-l L
のように行列の乗算をすることによって求め（ステップ S309)、雑音行列 Nを、
N = X - AS
のように行列の演算をすることによって求め（ステップ S310)、行列 A, S (= S ), Nを し 結果として RAM等内に格納することによって出力して (ステップ S 311)、本処理を終 了する。したがって、 CPUが RAM等と共働して、結果出力部として機能することとな る。
[0091] さて、上記のような繰返しは、信号分離を以下のように行うことに相当する。なお、以 下の式変形では、理解を容易にするため、雑音行列については表記を省略している
X =
S = A S
0 1
S = A S
1 2：
S = A S
2 3：
S = A S ；
L-2 L-l L-l
S = A S
L-l L L
= A S
L
[0092] これらをまとめると、以下のように表記することができる。
X = A A A - --A S
1 2 3 L
= AS
[0093] ステップ 309における乗算は、行列 Aを上記のように求めていることに相当するので ある。
[0094] (所定の行列の演算） 上記実施形態では、所定の行列演算 c(?)として、行列の積を用いていた。すなわ ち、行列 Zの i行 t列の要素を Z[i,t]とし、行列 Zの行数を row(Z)、行列 Zの列数を col(Z)と 書くとき、行列 Z, Wであって、 col(Z) = row(W)が成立する行列 Z, Wについて、所定の 演算は、
c(Z,W)[i,t] =∑ col(z) Z[i,k]W[k,t]
k=l
と定義される。
[0095] c(' , ·)として、行列のコンボリューシヨン演算を考える場合は、 2次元の行列（2階の テンソル）と 3次元の行列（3階のテンソル）とを考慮した、いわゆるテンソル演算を行う ことになる。このコンボリューシヨン演算される行列を、混合作用素（mixing operator) と呼ぶ。
[0096] 観測信号行列 Xと源信号行列 Sは、上記のように 2次元の行列であるが、混合比行 列 Aは 3次元の行列となる。 3次元の行列の要素を、上記と同様に [· , · , ·]の形式で表 現することとし、また、テンソル演算の慣習にしたがって、要素が定義されていない範 囲の値を 0とおくこととし、∑の添字がテンソル内のインデックスとしてとりうる範囲をす ベて走査した総和をとることとする。
[0097] すると、 3階のテンソル Aと 2階のテンソル Sとのコンボリューシヨンは、 [数 1]のように 定義される。
[0098] [数 1] c(AyS) [i, t] = y Y A[ , k, s]S[k, t - s
[0099] 一方、 3階のテンソル Aと 3階のテンソル Bとのコンボリューシヨンは、 [数 2]のように 定義される。
[0100] [数 2]
Figure imgf000017_0001
[0101] また、非線型リンク関数を用いる場合、
c(A,S)[i,t] = K∑ Z[i,k]W[k,t]) のように演算を定義する。関数 ·)としては、典型的には、ニューラルネットの分野でよ く用いられるステップ関数の原点付近を滑らかにした関数を利用することができる。た とえば、
ί χ) = 2 arctan(k χ〃 π
において、 kを十分に大きくしたような関数である。ただし一般には、任意の非線型関 数を適用することが可能である。
[0102] (実験結果）
以下では、上記の実施形態と従来の手法とを比較する各種の実験結果につ!、て説 明する。なお、以下の実験は、すべて計算機シミュレーションによっている。
[0103] 第 1の実験は、複数の音源から発せられる音波（一般には、信号源から発せられる 源信号の波動)を、近い位置に配置されたマイクで集めたとき (一般には、受信機や センサーで波動を観測して観測信号を得るとき）に相当する信号分離の様子を対比 するものである。以下では、 NMF (Non-negative Matrix Factorization)による信号分 離の分野で用いられるベンチマーク用データを用いて実験する。
[0104] 図 4は、源信号の波形の様子を示すグラフであり、図 5は、観測信号の波形の様子 を示すグラフである。
[0105] 図 4に示す源信号の波形に対して、条件数が 15,000程度のヒルベルト行列を用い て混合を行うと、図 5に示すような、互いに波形の類似した観測信号が得られる。
[0106] 図 6は、 L = 10として共役勾配法を用いて本実施形態の手法により 4つの信号に独 立成分分析を行った場合の波形のグラフであり、図 7は、これに対して、従来の手法 により 4つの信号に独立成分分析を行った場合の波形のグラフである。
[0107] 従来の手法（図 7)では、元の波形（図 4)とはまったく異なる成分に分解がされてし まっているが、本実施形態の手法（図 6)によると、元の波形（図 4)に極めて近いもの が得られて 、ることがわ力る。
[0108] 第 2の実験は、雑音が大き!/、場合の信号分離の様子を示すものである。
[0109] 図 8は、源信号の波形の様子を示すグラフであり、図 9は、観測信号の波形の様子 を示すグラフである。
[0110] 図 8中の上段 alに示すように、雑音レベルがかなり大きいため、図 9に示すように、 観測信号にも雑音が乗って 、ることがわかる。
[0111] 図 10は、 L = 5として甘利の a - divergenceNMFアルゴリズムを反復推測部におい て用いた本実施形態の手法により得られた波形のグラフであり、図 11は、 SIR値 (Sin gnal to Interference Ratioノ 不すグフフで &)る。
[0112] 図 8と対比すると、極めて類似したものが得られていることがわかり、また、 SIR値も 4
0dB〜53dB程度と、性能が良いことがわかる。
[0113] 図 12は、 NMF Lee-Seungアルゴリズムという従来の手法により得られた波形である
[0114] 従来手法は、 SIR値が 8dB未満となり、本実施形態の手法の性能の高さがわかる。
[0115] 第 3の実験は、画像処理における適用である。
[0116] 図 13は、源信号となる画像を示す図であり、図 14は、観測信号に相当する画像を 示す図であり、図 15は、 L = 5として本実施形態の手法を採用した場合に得られる源 信号の様子であり、図 16は、 L = 2として本実施形態の手法を採用した場合に得られ る源信号の様子であり、図 17は、 NMF Lee-Seungアルゴリズムという従来の手法によ り得られた源信号の画像であり、図 18は、本実施形態の手法における信号の SIR値 の値を示すグラフである。
[0117] 図 13に示すように、源信号は 4つであり、図 14に示すように、これらを 9個のチャン ネルで観測する。
[0118] L = 5として本発明の手法により独立成分分析すると、図 15に示すように、源信号が ほぼ完璧に復元される。 SIR値は、図 18に示すように、 46dB程度である。
[0119] 一方、 L = 2として本発明の手法により独立成分分析すると、図 16に示すように、分 離は不十分であり、 SIR値も 15dB未満と低い。 Lの値を増やすことの効果がよくあら われている。
[0120] 従来手法による分離の結果は、図 17に示すようになり、やはり分離は不十分で、 SI R値は 10db未満である。
産業上の利用可能性
[0121] 本発明によれば、独立成分分析等の信号分離技術において、精度をできるだけ向 上し、収束をできるだけ速くするのに好適な、信号分離装置、信号分離方法、ならび に、これらをコンピュータもしくはディジタル信号プロセッサ上にて実現するプログラム を記録したコンピュータ読取可能な情報記録媒体、ならびに、当該プログラムを提供 することができる。

Claims

請求の範囲
[1] m行 T列の観測信号行列 Xを、 m行 n列の混合比行列 Aと n行 T列の源信号行列 Sと m 行 T列の雑音行列 Nとであって、所定の行列演算 c(' , ·)につ 、て
X = c(A,S) + N
を満たす行列 A, Sに分離する信号分離装置であって、
m行 T列の行列と m行 n列の行列と n行 T列の行列とを受け付けてスカラー値を返す 行列関数 J(' , · , ·)に対して、 m行 T列の行列と m行 n列の行列と n行 T列の行列とが与え られると、当該 m行 T列の行列を固定したときに当該行列関数 ， ·， ·)が返すスカラー 値を極値化する弓 Iき数の m行 n列の行列と n行 T列の行列との組合せを、所定の繰返 し回数で反復推測する反復推測部、
L個の行列 A , A , A , · ··, Aであって、
1 2 3 L
c(A ,c(A ,c(A ,c(- --,c(A , A )· ··))))
1 2 3 L-l L
が m行 n列の行列となる L個の行列を繰返しにより計算する繰返計算部であって、 S = X
0
として、 i (= 1, 2, 3, · ··, L)回目の繰返しにおいて、前記反復推測部に、 i回目の繰返 しにおける行列関数 J (· , · , ·)であって、 m行 T列の行列と m行 n列の行列と n行 T列の行 列とを受け付けてスカラー値を返す行列関数 J (· , · , ·)に対して、 m行 T列の行列 S と
i i-1 m行 n列の所定の初期値行列と n行 T列の所定の初期値行列とを与えて、反復推測さ せ、得られた m行 η列の行列を Αとし、 n行 T列の行列を Sとする繰返計算部、
前記求められた L個の行列 A , A , A , · ··, Aから、 m行 n列の混合比行列 Aを
1 2 3 L
A = c(A ,c(A ,c(A ,c(- --,c(A ,A )· ··))))
1 2 3 L-l L
により計算し、
前記求められた n行 T列の行列 Sを、 n行 T列の源信号行列 Sとして、行列 A, Sを出 し
力する結果出力部
を備えることを特徴とする信号分離装置。
[2] 請求項 1に記載の信号分離装置であって、
前記求められた L個の行列 A , A , A , · ··, Aと n行 T列の行列 Sとが所定の収束条
1 2 3 し し
件を満たすまで、前記繰返計算部に当該繰返しをやり直させ、 i回目の繰返しにおけ る m行 n列の所定の初期値行列および n行 T列の所定の初期値行列として、やり直し の前回におけるへおよびやり直しの前回における を与えさせる繰返制御部
をさらに備え、
前記結果出力部は、当該所定の収束条件が満たされると、行列 A, Sを計算して出 力する
ことを特徴とする信号分離装置。
[3] 請求項 1に記載の信号分離装置であって、
前記結果出力部は、 m行 T列の雑音行列 Nを、
N = X - c(A,S)
によりさらに計算して、行列 Nをさらに出力する
ことを特徴とする信号分離装置。
[4] 請求項 1に記載の信号分離装置であって、当該所定の行列演算 , ·)は、行列の 禾 M "ある
ことを特徴とする信号分離装置。
[5] 請求項 1に記載の信号分離装置であって、
当該所定の行列演算 c(' , ·)は、コンボリューシヨン演算である
ことを特徴とする信号分離装置。
[6] 請求項 1に記載の信号分離装置であって、
当該所定の行列演算 c(' , ·)は、非線型リンク関数である
ことを特徴とする信号分離装置。
[7] 請求項 1に記載の信号分離装置であって、
前記反復推測部における反復回数は、 i (= 1, 2, 3, · ··, L)回目のいずれの繰返し においても 1回である
ことを特徴とする信号分離装置。
[8] 請求項 1に記載の信号分離装置であって、
前記反復推測部に対する i (= 1, 2, 3, · ··, L)回目の行列関数 J (' , ' , ')は、いずれも 等しい
ことを特徴とする信号分離装置。 [9] m行 T列の観測信号行列 Xを、 m行 η列の混合比行列 Αと、 n行 T列の源信号行列 Sと
、 m行 T列の雑音行列 Nであって、所定の行列演算 c(' , ')について
X = c(A,S) + N
を満たす行列 A, S, Nに分離し、反復推測部、繰返計算部、結果出力部を有する信 号分離装置にて実行される信号分離方法であって、
前記反復推測部が、 m行 T列の行列と m行 n列の行列と n行 T列の行列とを受け付け てスカラー値を返す行列関数 J(' , · , に対して、 m行 T列の行列と m行 n列の行列と n行 T列の行列とが与えられると、当該 m行 T列の行列を固定したときに当該行列関数 J(' , • , が返すスカラー値を極値ィ匕する引き数の m行 n列の行列と n行 T列の行列との組 合せを、所定の繰返し回数で反復推測する反復推測工程、
前記繰返計算部が、 L個の行列 A , A , A , · ··, Aであって、
1 2 3 し
c(A ,c(A ,c(A ,c("*,c(A ，A )· ··》》
1 2 3 L-l し
が m行 n列の行列となる L個の行列を繰返しにより計算し、
S = X
0
として、 i (= 1, 2, 3, · ··, L)回目の繰返しにおいて、前記反復推測部に、 i回目の繰返 しにおける行列関数 J (· , · , ·)であって、 m行 T列の行列と m行 n列の行列と n行 T列の行 列とを受け付けてスカラー値を返す行列関数 J (· , · , ·)に対して、 m行 T列の行列 S と
i i-1 m行 n列の所定の初期値行列と n行 T列の所定の初期値行列とを与えて、反復推測さ せ、得られた m行 η列の行列を Αとし、 n行 T列の行列を Sとする繰返計算工程、 前記結果出力部が、前記求められた L個の行列 A , A , A , · ··, Aから、 m行 n列の
1 2 3 し
混合比行列 Aを
A = c(A ,c(A ,c(A ,c(- --,c(A ,A )· ··))))
1 2 3 L-l L
により計算し、前記求められた n行 T列の行列 Sを、 n行 T列の源信号行列 Sとして、行 し
列 A, Sを出力する結果出力工程
を備えることを特徴とする信号分離方法。
[10] 請求項 9に記載の信号分離方法であって、当該信号分離装置は、繰返制御部をさ らに有し、
前記繰返制御部が、前記求められた L個の行列 A , A , A , · ··, Aと n行 T列の行列 sとが所定の収束条件を満たすまで、前記繰返計算部に当該繰返しをやり直させ、 i し
回目の繰返しにおける m行 n列の所定の初期値行列および n行 T列の所定の初期値 行列として、やり直しの前回における Aおよびやり直しの前回における Sを与えさせる 繰返制御工程
をさらに備え、
前記結果出力工程では、当該所定の収束条件が満たされると、行列 A, Sを計算し て出力する
ことを特徴とする信号分離方法。
コンピュータに、
m行 T列の観測信号行列 Xを、 m行 n列の混合比行列 Aと、 n行 T列の源信号行列 Sと 、 m行 T列の雑音行列 Nであって、所定の行列演算 c(' , ')について
X = c(A,S) + N
を満たす行列 A, S, Nに分離させるプログラムであって、
m行 T列の行列と m行 n列の行列と n行 T列の行列とを受け付けてスカラー値を返す 行列関数 J(' , · , ·)に対して、 m行 T列の行列と m行 n列の行列と n行 T列の行列とが与え られると、当該 m行 T列の行列を固定したときに当該行列関数 ， ·， ·)が返すスカラー 値を極値化する弓 Iき数の m行 n列の行列と n行 T列の行列との組合せを、所定の繰返 し回数で反復推測する反復推測部、
L個の行列 A , A , A , · ··, Aであって、
1 2 3 L
c(A ,c(A ,c(A ,c(- --,c(A , A )· ··))))
1 2 3 L-l L
が m行 n列の行列となる L個の行列を繰返しにより計算する繰返計算部であって、 S = X
0
として、 i (= 1, 2, 3, · ··, L)回目の繰返しにおいて、前記反復推測部に、 i回目の繰返 しにおける行列関数 J (· , · , ·)であって、 m行 T列の行列と m行 n列の行列と n行 T列の行 列とを受け付けてスカラー値を返す行列関数 J (· , · , ·)に対して、 m行 T列の行列 S と
i i-1 m行 n列の所定の初期値行列と n行 T列の所定の初期値行列とを与えて、反復推測さ せ、得られた m行 η列の行列を Αとし、 n行 T列の行列を Sとする繰返計算部、
前記求められた L個の行列 A , A , A , · ··, Aから、 m行 n列の混合比行列 Aを A = c(A ,c(A ,c(A ,c(- --,c(A ,A )· ··))))
1 2 3 L-l L
により計算し、前記求められた n行 T列の行列 Sを、 n行 T列の源信号行列 Sとして、行 し
列 A, Sを出力する結果出力部
として機能させることを特徴とするプログラムを記録したコンピュータ読取可能な情 報記録媒体。
[12] 請求項 11に記載のプログラムを記録したコンピュータ読取可能な情報記録媒体で あって、当該コンピュータを、
前記求められた L個の行列 A , A , A , · ··, Aと n行 T列の行列 Sとが所定の収束条
1 2 3 し し
件を満たすまで、前記繰返計算部に当該繰返しをやり直させ、 i回目の繰返しにおけ る m行 n列の所定の初期値行列および n行 T列の所定の初期値行列として、やり直し の前回における Aおよびやり直しの前回における Sを与えさせる繰返制御部
としてさらに機能させ、
前記結果出力部は、当該所定の収束条件が満たされると、行列 A, Sを計算して出 力する
ように機能させることを特徴とするプログラムを記録したコンピュータ読取可能な情 報記録媒体。
[13] コンピュータに、
m行 T列の観測信号行列 Xを、 m行 n列の混合比行列 Aと、 n行 T列の源信号行列 Sと 、 m行 T列の雑音行列 Nであって、所定の行列演算 c(' , ')について
X = c(A,S) + N
を満たす行列 A, S, Nに分離させるプログラムであって、
m行 T列の行列と m行 n列の行列と n行 T列の行列とを受け付けてスカラー値を返す 行列関数 J(' , · , ·)に対して、 m行 T列の行列と m行 n列の行列と n行 T列の行列とが与え られると、当該 m行 T列の行列を固定したときに当該行列関数 ， ·， ·)が返すスカラー 値を極値化する弓 Iき数の m行 n列の行列と n行 T列の行列との組合せを、所定の繰返 し回数で反復推測する反復推測部、
L個の行列 A , A , A , · ··, Aであって、
1 2 3 L
c(A ,c(A ,c(A ,c(- --,c(A , A )· ··)))) が m行 n列の行列となる L個の行列を繰返しにより計算する繰返計算部であって、 S = X
0
として、 i (= 1, 2, 3, · ··, L)回目の繰返しにおいて、前記反復推測部に、 i回目の繰返 しにおける行列関数 J (· , · , ·)であって、 m行 T列の行列と m行 n列の行列と n行 T列の行 列とを受け付けてスカラー値を返す行列関数 J (· , · , ·)に対して、 m行 T列の行列 S と
i i-1 m行 n列の所定の初期値行列と n行 T列の所定の初期値行列とを与えて、反復推測さ せ、得られた m行 η列の行列を Αとし、 n行 T列の行列を Sとする繰返計算部、 前記求められた L個の行列 A , A , A , · ··, Aから、 m行 n列の混合比行列 Aを
1 2 3 L
A = c(A ,c(A ,c(A ,c(- --,c(A ,A )· ··))))
1 2 3 L-l L
により計算し、前記求められた n行 T列の行列 Sを、 n行 T列の源信号行列 Sとして、行 し
列 A, Sを出力する結果出力部
として機能させることを特徴とするプログラム。
請求項 13に記載のプログラムであって、当該コンピュータを、
前記求められた L個の行列 A , A , A , · ··, Aと n行 T列の行列 Sとが所定の収束条
1 2 3 し し
件を満たすまで、前記繰返計算部に当該繰返しをやり直させ、 i回目の繰返しにおけ る m行 n列の所定の初期値行列および n行 T列の所定の初期値行列として、やり直し の前回における Aおよびやり直しの前回における Sを与えさせる繰返制御部
としてさらに機能させ、
前記結果出力部は、当該所定の収束条件が満たされると、行列 A, Sを計算して出 力する
ように機能させることを特徴とするプログラム。
PCT/JP2007/059172 2006-04-27 2007-04-27 信号分離装置、信号分離方法、情報記録媒体、ならびに、プログラム Ceased WO2007126054A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US12/226,651 US8285773B2 (en) 2006-04-27 2007-04-27 Signal separating device, signal separating method, information recording medium, and program
EP07742607A EP2012555A4 (en) 2006-04-27 2007-04-27 SIGNAL DISCONNECTION, SIGNAL SEPARATION, INFORMATION RECORDING MEDIUM AND PROGRAM

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2006124167A JP4952979B2 (ja) 2006-04-27 2006-04-27 信号分離装置、信号分離方法、ならびに、プログラム
JP2006-124167 2006-04-27

Publications (1)

Publication Number Publication Date
WO2007126054A1 true WO2007126054A1 (ja) 2007-11-08

Family

ID=38655575

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2007/059172 Ceased WO2007126054A1 (ja) 2006-04-27 2007-04-27 信号分離装置、信号分離方法、情報記録媒体、ならびに、プログラム

Country Status (4)

Country Link
US (1) US8285773B2 (ja)
EP (1) EP2012555A4 (ja)
JP (1) JP4952979B2 (ja)
WO (1) WO2007126054A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109192220A (zh) * 2018-08-02 2019-01-11 福建师范大学福清分校 一种比例微分控制的混合信号自适应快速分离方法
CN110969248A (zh) * 2018-09-28 2020-04-07 中国石油化工股份有限公司 基于rbf-a2c的强化学习方法及系统

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4958241B2 (ja) * 2008-08-05 2012-06-20 日本電信電話株式会社 信号処理装置、信号処理方法、信号処理プログラムおよび記録媒体
US8340943B2 (en) * 2009-08-28 2012-12-25 Electronics And Telecommunications Research Institute Method and system for separating musical sound source
KR101225932B1 (ko) 2009-08-28 2013-01-24 포항공과대학교 산학협력단 음악 음원 분리 방법 및 장치
KR101272972B1 (ko) 2009-09-14 2013-06-10 한국전자통신연구원 음원 데이터베이스를 사용하지 않는 음악 음원 분리 방법 및 장치
KR20120031854A (ko) * 2010-09-27 2012-04-04 한국전자통신연구원 시간 및 주파수 특징을 이용하는 음악 음원 분리 장치 및 방법
JP6054142B2 (ja) * 2012-10-31 2016-12-27 株式会社東芝 信号処理装置、方法およびプログラム
CN105989851B (zh) * 2015-02-15 2021-05-07 杜比实验室特许公司 音频源分离
US9668066B1 (en) * 2015-04-03 2017-05-30 Cedar Audio Ltd. Blind source separation systems
CN110269609B (zh) * 2019-07-08 2021-09-28 江苏理工学院 基于单通道的脑电信号中眼电伪迹分离方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003141102A (ja) * 2001-11-05 2003-05-16 Inst Of Physical & Chemical Res 化学物質分類装置、化学物質分類方法、ならびに、プログラム
JP2003271166A (ja) * 2002-03-14 2003-09-25 Nissan Motor Co Ltd 入力信号処理方法および入力信号処理装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6535666B1 (en) * 1995-06-02 2003-03-18 Trw Inc. Method and apparatus for separating signals transmitted over a waveguide
US6026217A (en) * 1996-06-21 2000-02-15 Digital Equipment Corporation Method and apparatus for eliminating the transpose buffer during a decomposed forward or inverse 2-dimensional discrete cosine transform through operand decomposition storage and retrieval
US7158559B2 (en) * 2002-01-15 2007-01-02 Tensor Comm, Inc. Serial cancellation receiver design for a coded signal processing engine
GB0229473D0 (en) * 2002-12-18 2003-01-22 Qinetiq Ltd Signal separation system and method
US7039546B2 (en) * 2003-03-04 2006-05-02 Nippon Telegraph And Telephone Corporation Position information estimation device, method thereof, and program
FI20031079A0 (fi) * 2003-07-16 2003-07-16 Nokia Corp Menetelmä tiedonsiirtoresurssien kontrolloimiseksi, sekä kontrolleri
KR100580843B1 (ko) * 2003-12-22 2006-05-16 한국전자통신연구원 Ｖ―ｂｌａｓｔ에서 채널전달함수행렬 처리장치 및 그의처리방법

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003141102A (ja) * 2001-11-05 2003-05-16 Inst Of Physical & Chemical Res 化学物質分類装置、化学物質分類方法、ならびに、プログラム
JP2003271166A (ja) * 2002-03-14 2003-09-25 Nissan Motor Co Ltd 入力信号処理方法および入力信号処理装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP2012555A4 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109192220A (zh) * 2018-08-02 2019-01-11 福建师范大学福清分校 一种比例微分控制的混合信号自适应快速分离方法
CN109192220B (zh) * 2018-08-02 2022-11-11 福建师范大学福清分校 一种比例微分控制的混合信号自适应快速分离方法
CN110969248A (zh) * 2018-09-28 2020-04-07 中国石油化工股份有限公司 基于rbf-a2c的强化学习方法及系统
CN110969248B (zh) * 2018-09-28 2023-04-07 中国石油化工股份有限公司 基于rbf-a2c的强化学习方法及系统

Also Published As

Publication number Publication date
US8285773B2 (en) 2012-10-09
JP4952979B2 (ja) 2012-06-13
EP2012555A1 (en) 2009-01-07
US20090234901A1 (en) 2009-09-17
EP2012555A4 (en) 2011-06-08
JP2007298564A (ja) 2007-11-15

Similar Documents

Publication Publication Date Title
WO2007126054A1 (ja) 信号分離装置、信号分離方法、情報記録媒体、ならびに、プログラム
Oymak et al. On the role of attention in prompt-tuning
US11238881B2 (en) Weight matrix initialization method to improve signal decomposition
Hu et al. Stein neural sampler
Low et al. Complexity of implementing trotter steps
Vershynin Spectral norm of products of random and deterministic matrices
Zhang et al. Optimized sparse fractional Fourier transform: Principle and performance analysis
EP1290816B1 (en) Method of extracting a signal
EP2312576A2 (en) Method and system for reducing dimensionality of the spectrogram of a signal produced by a number of independent processes
Rössler Runge–Kutta methods for Itô stochastic differential equations with scalar noise
US20220130407A1 (en) Method for isolating sound, electronic equipment, and storage medium
Ito et al. Grouped generalized estimating equations for longitudinal data analysis
Zhao et al. Penalized pairwise pseudo likelihood for variable selection with nonignorable missing data
JP6099032B2 (ja) 信号処理装置、信号処理方法及びコンピュータプログラム
Zhou et al. Gaussian processes with errors in variables: Theory and computation
Phandoidaen et al. Empirical process theory for locally stationary processes
Chen et al. Online active regression
Moliner et al. Unsupervised estimation of nonlinear audio effects: Comparing diffusion-based and adversarial approaches
Hadar High-order hidden Markov models-estimation and implementation
Hara et al. Design of graph signal sampling matrices for arbitrary signal subspaces
Wu Convergence of the randomized block Gauss-Seidel method
Zhou et al. Gaussian processes with errors in variables: Theory and computation
Yakovlev et al. Implicit score matching meets denoising score matching: improved rates of convergence and log-density Hessian estimation
Dmitriev et al. Efficient sampling with discrete diffusion models: Sharp and adaptive guarantees
US20210117793A1 (en) Data processing system and data processing method

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 07742607

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 12226651

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2007742607

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE