The acoustic distance measurement (ADM) method based on phase interference between transmitted and reflected signal has been proposed as a method to accurately measure the distance. However, since its performance is deteriorated by ambient noise. In the previous study, complex ICA(Independent Component Analysis) was applied using adjacent 3ch microphone to reduce the noise. Although ICA achieves high performance for the noise reduction, it is not too practical method because of its high computational cost. This paper describes a practical method to reduce ambient noise in the ADM using adjacent 3ch observations. Concretely, the noise component is reduced by generating two difference signals among three signals of 3ch observations in the noisy environment. Then, the ADM using the crossspectrum method is applied to them. The difference signals can be obtained with a small computational cost. The validity of the proposed method is confirmed in the actual room under the condition where the SNR ratio is -10dB. As a result, the SNR is improved to 10.5 dB, demonstrating the validity of the proposed method.

Under ideal conditions with low noise, the acoustic distance measurement (ADM) method using the cross-spectrum for adjacent 2ch observations shows sufficient distance estimation performance. Increasing noise will deteriorate the distance estimation performance. To reduce noise, even if complex ICA is applied to adjacent 2ch noisy observations, sufficient separation performance for distance measurement cannot be exhibited, because the microphones are too close. Thus, if the ADM method using the cross-spectrum is applied to the separated signals of the adjacent 2ch noisy observations, an incorrect distance near 0 m is estimated because of the effect of measurement system. Therefore, by increasing the number of microphones (which increases information and can be expected to improve the separation performance), and applying complex ICA to improve the separation performance, ADM using the cross-spectrum for the new separation signals is proposed. This method enables us to obtain distance estimation that eliminates both environmental noise and the effect of the measurement system. The validity of the proposed method is also verified with the experiment of the actual distance measurement. As a result, it was found that the distance 1 m to the target can be measured under the noise of SNR (Signal to Noise Ratio) -10 dB.

The distance to a target is important information in various engineering fields. When measuring the distance to a target using audible sound, in general, an impulsive sound is transmitted toward the target and the time difference between the transmitted and reflected waves is obtained to determine the distance. However, this method cannot be applied to short distances where a reflected wave is buried in a transmitted wave. On the other hand, in the field of microwave radar, this distance measurement method has been proposed for plural targets, and can be applied to distances close to 0 m (i.e., a standing-wave radar). We have extended this distance measurement method so as to measure precisely the distance to the target in ranges close to 0 m by applying the principle of the standing wave radar to distance measurement using a sound signal. In this paper, first we will explain the principle of acoustical distance measurement methods based on phase interference using audible sound, and then by introducing an analytic signal, we describe in detail an acoustical distance measurement method, especially for a target that exists at a distance of around 0 m. Simulation results and experimental results are also shown.

我々はこれまで,可聴音域に拡張した送信波と反射波の位相干渉に基づく近距離推定法について検討してきた.本研究では,対象物のある環境とない環境での同期収録信号をスペクトル減算することで,バックグラウンド要素を除去する音響測距法を検討する.

特定の場所の音環境を立体的に再現する技術として,頭部伝達関数(Head Related Transfer Function:HRTF)を事前に計測して適用するバイノーラル方式と任意の位置にある複数の遠隔スピーカを利用するトランスオーラル方式が提案されている.バイノーラル方式はヘッドホンの着用が煩わしく,トランスオーラル方式は機器の個数や計算量の点で非常に複雑かつ大規模なシステムとなる問題がある.一方,我々はこれまで,スピーカとその近傍のマイクロホンで収録した直接音と反射音の位相干渉を利用した音響測距法を提案してきた.この音響測距法は,スピーカと受聴者の耳元間の距離計測に利用できる.この距離情報とスピーカ配置を考慮することで,トランスオーラル方式と比べて小規模な立体音場再現システムを実現できる可能性がある.そこで本研究では,受聴者の頭部の左右近傍に設置したスピーカ-マイクロホン対を利用した音響測距法に基づく立体音場再現手法を提案する.更に,理論の正当性を確認するための予備実験と,HRTFにより定位を付与した音楽による主観評価実験を行うことで,提案手法の有効性を確認する.

「怒り」，「喜び」，「悲しみ」などの感情を表現した自然音声の音響的特徴を「韻律（声の大きさ，高さ，発話速度の総称）」という視点から解析してきた。筆者らの研究の特徴は，より微妙なニュアンスの特徴まで調べることを目的に，それぞれの感情を，感情を含まない「平常」のほか，「弱い」，「中程度」，「強い」の4段階の度合に分類して解析を行い，それぞれの特徴を体系的に明らかにしてきたことである。
会話によるコミュニケーションにおいて，話者の心の内面にある感情がそのまま発声に反映され，更にそのまま聞き手に伝わるとは限らない。
本研究の目的は，(1)このような話者の内面の感情（話者の意図）がどの程度正しく聞き手に伝わるのか，(2)どのような表現のときに話者の意図した感情が正しく聞き手に伝わりやすいのか，を定量的に明らかにすることである。
そこでまず，発話の意図がどの程度聞き手の聴覚的印象と一致しているかを調

In many engineering fields, distance to target is fundamental and very important information. Acoustical signal, such as audible sound or ultrasonic, often plays an essential role in measurement of distance. To measure distance, the time delay of reflected wave measured with reference to transmitted pulse wave, i.e. time-of-flight (TOF) is used typically. This method, however, can not measure short distance because the transmitted wave, which has not attenuated enough as of reception of reflected waves, suppresses the reflected waves for short distance. Meanwhile, in the research field of microwave radar, to measure distances to multiple targets even if short, a novel method by using standing wave has been recently proposed. We apply the fundamental principle of this method to estimation of short distance for multiple targets using audible sound. In this paper, we focus our attention on audible band-limited noise. To verify the validity and effectiveness of our method, we perform computer simulations and experiments. As a result, it has been confirmed that short distances to multiple targets, such as the distances to two targets located at 0.1m and 1m, can be estimated properly by our method.

In order to establish an objective approach for evaluating human emotions, Emotional Spectrum Analysis Method (ESAM) has been proposed by Dr. T. Musha et al. It is possible to estimate emotion spectrums with four basic emotional states (i.e., anger/stress, joy, sadness and relaxation) by using ESAM and the measured data of brain waves. Now, we pay attention to the emotional evaluation of music sounds. More specifically, this paper describes the influence of a change of music sounds with different tones (from a minor scale to a major scale) on human emotions. The brain waves for six subjects are measured by electroencephalogram (EEG). Each emotion spectrum can be estimated by using the emotion spectrum analysis system developed by Brain Functions Laboratory, Inc. based on ESAM. The effectiveness of ESAM is experimentally confirmed by applying it to the emotion evaluation problem of the above two kinds of music sounds. We find out the possibility of ESAM for the objective evaluation in the field of music.

本研究では，風速から風雑音を推定し，その風雑音を風雑音下の観測信号から除去することにより信号のみを抽出する簡便な手法を提案する．具体的には，風速を入力とする線形・非線形のシステムモデルを用い風雑音の音圧レベルをまず推定する．ついで，その推定した風雑音を観測信号から音のエネルギースケールで減算し，風雑音成分を除去する．線形モデルとしてはMAフィルタ，非線形フィルタとしてニューラルネットワークを採用する．また，線形成分はMAフィルタで捉え，それからズレた非線形成分についてはニューラルネットワークで補正する組み合わせ手法も提案する．最後に，実際の風雑音混入下において音響信号を検出する問題に本手法を適用し，その有効性を検証した．(英文)

複数の混合音から複数の音源の位置を同定する手法を提案した。とくに今回は時間遅れを伴う混合を想定し、ブラインド分離により求まるパラメータを利用した。さらに逆フィルタ処理により信号のクリアな分離が可能となることも示した。（招待論文、英文）

近年脚光を浴びている独立成分分析（Independent Component Analysis; ICA）について、著者らの研究成果を中心に基礎から、とくに音響信号処理に関連して解説した。瞬時混合の場合だけでなく音響計測では無視できない畳み込み混合の場合にも焦点を当てた。

現実の音環境における不規則信号は、様々な環境要因により複雑な変動特性を示すのが常である。本論文では、このような実状況の環境騒音に対する統計的一評価手法を原理的に提案している。特に実用性に考察の主点を置き、ごく少数個のみのパラメータを持つ単純・標準化された結果を、まず理論的枠組みに採る。次に、単純から複雑への拡張原理に支えられた2種情報処理法として、ニューラルネットワークとファジー確率論理を採り上げ、現実の音環境における複雑さへも適用可能な形態へと方法論を拡張する。最後に音環境での実測データに適用することにより、方法論の有効性を3種の実際的場合へ適用して原理的に確認している。

騒音・振動環境分野では、自然要因の多様さや人間評価の非画一性に基づいて、様々な評価量が提案され、それぞれ固有な長短所を待ち合わせている。本研究では、L_x(x=5,50,95)評価量を一具体例にとり、各評価量間で説明変数から目的変数への一方向的な変換関係を求める回帰分析法とは異なり、互いに優劣なく双方向に変換可能な代表的一関連様式を見出す目的で主成分分析法にまず着目する。次いで、各L_x評価量間では本質的に非線形な関連様式を潜在させることから、線形関連様式を主とした公知の主成分分析法を非線形関連様式へ一般化できる一試みを、その原則論に力点をおいて提案する。更に、実L_xデータ間に適用しその実際的有効性をも実証する。

既発表の拡張型回帰分析法の解析基盤となるベイズ定理が原理的に情報量規範と等価な解を生むことを積極的に利用して、操作面上ベイズ定理より柔軟な極値探索型の情報量規範を適用することにより、複雑な住環境騒音を実用的に評価できる新たな一回帰分析法を提案した。情報量規範として入出力間の条件付エントロピーを最小にする評価規範を導入し、確率近似法を援用して入出力データを逐次採取することによってシステム同定ができる実用的手法を提案した。更にシステム同定時とは異なった変動形態を持つ他の任意入力がこの複雑なシステムに印加されたときの出力応答分布をきめ細かく予測できる一手法を見出した。本手法を実音響システムに対する応答出力分布の予測問題に適用し、その有効性の一端を実際的にも確認した。

前報告では、刻々の残響過程においてエネルギーの統計的な揺らぎ形態が特に個々のモーメント情報に中心をおいて時間と共にいかに推移してゆくかを考察した。前報告とは異なり、本報告ではまず、計測の実際と現実の揺らぎの複雑さに対応できる離散型の一般化Fokker-Planck方程式を導出し、残響室における残響曲線の揺らぎ分布全体に関する過渡的推移表現をSabineの残響則を反映した逐次形式の形で新たに導出する。具体化にあたっては、エネルギー変量における物理的対策の足場と評価におけるdB観測に整合できるよう、この確率分布の推移法則における初期分布の枠組として、特に対数正規分布を初項に持つ一般的な展開分布表示を採用している。最後に、残響室内における実測例に適用して本手法の有効性を実験的にも検証している。

任意周波数帯域幅の白色雑音入力を、拡散音場下の室に加えたとき、その不規則応答波が、任意K個の観測点において時間・空間的にどのような揺らぎ相関を持つか、その確率特性について理論的に考察する。具体的には、まずこれらの各種相関情報をすべて洩れなく導出できる任意多変量結合確率分布表示を一般表現の形で導出する。特に、2変量、3変量、4変量の場合について詳細に考察を行った。また、線形・非線形の各種相関情報を浮き彫りにするため、結合モーメント、条件付確率密度関数、条件付モーメントについての考察も行う。最後に、実際に残響室内において実験を行い、理論的に予測された曲線と実測値とを比較検討し、理論の有効性をも確認した。

前報告では、メモリある音・振動環境システムの簡便な1変量出力応答分布がメモリ効果とゼロメモリ型非線形効果を機能的に分離し近似表現ではあるが両者を反映した統一的形態で導出された。特にエネルギー変量のみに着目した前報告とは異なり、本報告ではまず、エネルギースケールにおける物理変量に足場を置きながらもdBスケールでの観測データへも整合できる揺らぎ分布表現の一般的枠組みとして、対数正規分布を表現の基幹にとり、更に揺らぎの多様さへ互いに冗長度なきよう階層的に対処できる直交展開型分布表現を新たに導入した。本手法の有効性は、ロック音楽と交通騒音の侵入した室内音響における実問題へ適用して実験的にも確認された。

任意帯域幅の白色雑音を入力した拡散音場下の2観測点間で、音の強さに関する互いの揺らぎが持つ空間相関と各観測点で個々に示す揺らぎ分布の間の関連性についてまず理論的に考察する。両者を関連付けた揺らぎの結合分布が、入力の任意周波数帯域幅と計器の自乗機能や任意時定数の影響を具体的に反映させた2変量統計的Laguerre展開型の分布表示として導出される。この表示を基盤にして2観測点間の相関効果を考察するとき、低次から高次へ線形相関のべき級数の形で階層的に個々に調査する分析的立場と、条件付き確率密度の分布関数全体に着目する総合的立場の、対照的な二つの見地をとる。最後に、残響室内の実測データに本理論を適用してその有効性を実験的にも検証する。

The additive model with external fluctuation factors and the multiplicative model with internal fluctuation factors are well-known as two fundamental models describing the stationary and non-stationary random phenomena in the actual engineering field. An actual stochastic systems described by an additive model in a short time interval must be very often described in a long time interval by a mixed type stochastic model containing the fluctuation of multiplicative factor, because the internal parameter of the system changes with the passage of time. In this paper, first, by introducing a mixed type stochastic model for the actual random process, a unified expression is proposed for the probability distribution function of the above non-stationary composite wave measured by an arbitrary non-linear observation system. Next, by use of the above theory, two types of new methods identifying the probability distribution function of a specific factor embedded in the composite wave are derived in a general form of series expansion expression. Finally, by applying the above theoretical results to various types of L_<eq> evaluation problems of road traffic noises and industry noises based on the actual noise observations corrupted by the background noise, the validity of proposed evaluation methods is also experimentally confirmed.

In general, the descriptive approach on the stochastic evaluation is typically effective for the complicated situation of actual sound insulation systems with an actual random noise excitation which is impossible to evaluate by use of the well-known structual approach. In this paper, after once introducing descriptively a linear regression type system model between input and output, the system parameters have been estimated not in a sense of the usual least square error evaluation but in a sense of laying stress on the edge or the whole shape of probability distribution function. More concretely, this estimation method has been originally proposed by use of a stochastic approximation method for the purpose of extracting more effective information latent in the remainder fluctuation around the above linear regression, under a new error criterion based on equalizing a noise evaluation index L_aα and a Kullbacks information quantity for modelled and actual systems. Furthermore, in an actual case with an existence of background noise, a simplified trial of removing this background noise is also discussed. Finally, the effectiveness of the proposed evaluation methods has been experimentally confirmed too by applying it to the actually complicated sound environmental system.

As is well-known, the wave form of random noise and vibration appearing in our actual enviroment shows usually the variety and the complexity in its fluctuation pattern. The statistical methods of analyzing and evaluating the above random wave form were studied by many researchers, from the viewpoint of signal analysis. In this paper, after first setting a generalized method of statistically evaluating the acuity of fluctuating random wave based on the information on each order of differntial type state variables, especially for the stational random wave of arbitarary distribution type, a new trial toward the statistical evaluation on the locations of level crossing, peak and trough has been concretely considered as the special cases, in close connection with the above differential type state variables. That is, in the proposed theory, the effect of usual linear correlations and non-linear correlations among two and/or many sampled points of random wave on the above probability distribution form of level crossing, peak and trough has been fully discussed under introduction of multivariate joint probability distribution. Finally, the practical effectiveness of the present theory has been experimentally confirmed too by applying it to the acrual road traffic noise wave obeserved ina suburb of a large city.

スペクトル解析に関係する専門技術について、ハンドブックの一部を執筆した。内容は、脳波解析であり、「コンピュータの進歩と脳波解析」、「脳波の時間-周波数解析」、「独立成分分析」 に関し、それらの原理を説明するとともに、我々が行った脳波解析の結果も示した。

可聴音を用いて音源から物体までの距離を推定する手法について検討した。

感情音声の合成を目的として，声優が発声する自然音声の声質の特徴解析を行っている。今回は，音声信号のスペクトル傾斜に着目し，感情の種類との関係を調べた。その結果，「平常」に比べて「激怒」の方がスペクトル傾斜が大きいことがわかった。同様に，「喜び」におけるスペクトル傾斜も「平常」の場合より大きい。このことは，「激怒」および「喜び」の方が「平常」より高域が強調された声質に変化していることを意味している。２名の声優の音声５単語についてほぼ同様の結果が得られた。

音声に含まれる感情の種類・度合によって，音源に含まれる雑音のレベルがどのように変化するかを調べる研究の一環として，雑音量を定量化する指標である N/S 比について検討を行った。その中で，筆者らが提案した櫛形フィルタを用いた N/S 比の精度向上と計算の効率化を目的として，櫛の位置と幅を自動決定する手法を提案した。この手法では，残差パワーの周波数に関する二階微分値のピークを櫛の位置と幅の情報として用いている。この手法により，櫛形の自動調整が可能であるとの見通しがついた。

歌声には、人それぞれに声質がある。例えば、太い声、細い声、かすれ声、声域などが挙げられるが、このような個人の声質に合った歌曲を選曲するのはなかなか困難である。本研究では、歌声の声質の中でも特に声域に注目し、個人の声域により近い音域の歌曲の自動検索や指定曲を個人の声域の高さにより近くなるように移調するシステムを提案する。

音響分野に限らず信号の性質を調べるためにさまざまな方法が提案されており、なかでも時間-周波数分析法はしばしば用いられる。しかし、それらは一般に時間分解能と周波数分解能の間には不確定性の関係がある。近年、より精密な時間的変化，周波数的変化を同時に見るために時間-周波数でのパワー変化の分布（すなわち、正値時間-周波数分布）に着目する手法が提案されている。前報に引き続き本研究では、Copulaに基づく正値時間-周波数分布を導入し、音響信号の分析のために音声解析への応用を検討する。

対象物までの距離を推定するために超音波等がしばしば用いられ，その推定原理としては，送信パルスと対象物からの反射波との時間差が利用されている．しかし，この手法では対象物が数10cmに存在する状況では反射波が送信波に埋もれてしまい推定できなくなる．これに対してマイクロ波レーダの分野では，定在波を利用して近距離計測を行う手法が知られている．本研究では可聴音に着目し，スピーカから放射された音と対象物からの反射音の干渉によって生じる定在波を利用した距離の推定法を基礎から検討する．

計測自動制御学会関西支部からの依頼を受け講習会の講師を務めた。以下はその要旨である。 複数の独立信号が混合された観測信号から独立成分のみを取り出す独立成分分析 (Independent Component Analysis: ICA) が盛んに研究されている．これは単純な問題設定の割に解法がバラエティに富むためと，なにより潜在的な応用範囲が広いためと思われる．時間差のない信号の混合を Instantaneous mixture と呼び，時間差やシステムの影響を受けた場合を Convolutive mixture と呼んで区別している．このような状況において信号を分離することをブラインド分離(Blind Source Separation: BSS)問題と呼び，とくに Convolutive mixture の場合に信号源まで復元することをBlind Source Deconvolution問題と呼んでいる．ICAは原理的にはどちらの混合に対しても用いることができる． ICA はさまざまな分野で用いられるが，とくに音響信号については当初から特にカクテルパーティ効果などと関連して，音声分離のシミュレーシ

混合に時間差などの影響があるブラインド分離問題に対しては、周波数領域における複素ICA(独立成分分析) がよく用いられ、著者らも複素ICAの一解法を提案している。 本研究では、複素 ICA に基づいて時間差混合信号を分離する際得られる分離行列から遅延時間・振幅比の推定を行う。これを用いて混合信号の逆フィルタを構成し、実信号として帯域歪みのない分離・再生を目指す。

混合信号のブラインド分離については多くの研究が報告されており，著者らも以前に瞬時混合に対して Hermite モーメントと回転変換により，畳み込み混合に対して複素Hermiteモーメントとユニタリ変換(SU(2))により取り扱った． 本報告では，3つの複素信号が混合された場合に拡張する．具体的には，規格直交化された観測信号に対し，複素Hermite モーメントに基づく評価関数を最大化する3次元のユニタリ行列(SU(3))を推定する．最後に本手法の正当性をシミュレーションにより検証する．

音響分野に限らず信号の性質を調べるためにさまざまな方法が提案されている。なかでも時間-周波数分析法はしばしば用いられるが、一般に時間分解能と周波数分解能の間には不確定性の関係がある。近年、より精密な時間的変化，周波数的変化を同時に見るために時間-周波数でのパワー変化の分布（すなわち、正値時間-周波数分布）に着目する手法が提案されている。 本研究では、正値時間-周波数分布による音響信号の分析のために、Copulaに基づく正値時間-周波数分布の具体的な算出法を検討し、テンプレートとしてスペクトログラムを採用し、楽音への適用を行う。

混合信号のブラインド分離については多くの研究が報告されており，著者らも以前に複数の独立信号の同時混合を Hermite モーメントと回転変換により取り扱った．また，混合に時間差や残響などの影響がある場合には，混合信号を複素狭帯域化した後，分離する周波数（複素）領域での ICA も提案している．2信号に対しては複素HermiteモーメントとUnitary変換(SU(2))を利用するが，本報告では 3信号が畳み込み混合された場合に拡張する．具体的には，規格直交化された観測信号に対し，複素Hermite モーメントに基づく評価関数を最大化する3次元のUnitary行列(SU(3))を推定する．最後に本手法の正当性を検証するため，計算機シミュレーションを行う．

音響信号の性質を調べるためにさまざまな方法が提案されている。なかでも信号の時間的特徴と周波数領域での特徴を同時に調べることのできる時間-周波数分析法は強力なツールとしてしばしば用いられる。時間-周波数分析法には、スペクトログラム（短時間フーリエ変換）、Wigner分布、Choi-Williams分布などのいくつかの手法が知られている。一般に、時間-周波数分析ではエネルギー分布を用いるが、正値分布としての条件をすべては満足していない。 本報では、周辺分布などの条件を満足するCopulaに基づく正値時間-周波数分布に着目する。これにより、たとえばスペクトログラムなどに比べ、時間と周波数両方において分解能の向上が期待できる。最終的な目標を正値時間-周波数分布による音響信号の分析と定め、今回はCopulaに基づく正値時間-周波数分布の具体的な算出法を検討する。最後にテンプレートとしてスペクトログラムを採用し、簡単なシミュレーションや楽音への

電磁波の分野では、不規則に揺らいでいる媒質中を波が伝わるとき、波が小さな領域に閉じ込められる局在現象が知られている。よって、ランダム媒質あるいはランダムな装荷をもつ音場でも指数関数的な音の減衰があり、それを騒音抑圧に利用可能と思われる。 ランダム媒質による減衰効果を考察した既発表に引き続き、本研究では1次元音場の等価回路表現を用いて、ランダム装荷をもつダクトの減衰特性を検討し、シミュレーションにより正当性を検証する。

複数の信号が混合された観測信号から源信号を分離する問題をブラインド分離という。その解法の一つに源信号の統計的独立性に基づく独立成分分析(ICA)がある。本研究では、著者等の提案している Hermte モーメントに基づくICA 手法を、多チャンネルの脳波(EEG)データの分離に適用する。本手法により、ノイズ成分(アーチファクト)を取り除けることが分かった。 (英文)

源信号の統計的独立性に基づいて混合信号をブラインド的に分離する独立成分分析(ICA)が盛んに研究されている。著者等も、Hermite モーメントに基づく評価関数と回転変換により独自のICA アプローチを提案している。本研究では、対象信号以外の成分が存在する場合、評価関数にどのような影響を及ぼし、その結果、分離性能がどの程度劣化するかを理論的に考察した。シミュレーション実験でも良好な結果が得られている。

甲南大学の ORC(オープン・リサーチ・センター)の客員特別研究員として講演を依頼され、著者らの手法を中心に，ICA(独立成分分析) の基礎理論とその信号処理への応用について解説した．また，音声・画像でのシミュレーション結果や実際の生体信号への適用結果の一例も示した．

電磁波の分野では、不規則に揺らいでいる媒質中を波が伝わるとき、波が遠くまで伝わらず小さな領域に閉じ込められる現象が知られている(局在現象)。音場においても同様に指数関数的な減衰が期待でき、騒音抑圧に利用可能と思われる。本研究では既発表に引き続き、1次元音場の等価回路表現を用いて、ランダム性をもつダクトの減衰特性を検討し、シミュレーションにより正当性を検証する。

2信号の場合について複素Hermiteモーメントとユニタリ変換(SU(2))に基づく複素 ICA をすでに報告している。本報告では、3つの複素信号が混合された場合に拡張する。具体的には、規格直交化された観測信号に対し、複素Hermite モーメントに基づく評価関数を最大化する3次元のユニタリ行列(SU(3))を推定する。

平成１６年電気関係学会関西支部連合大会の信号処理分野のシンポジウム「信号処理の現状と未来」に招かれ、独立成分分析(ICA) の原理と音響信号処理への応用について講演した。さらに、音源方位推定への本手法の適用例も示した。

電磁波が、不規則に揺らいでいる媒質があると局在現象が起こる。本研究では1次元音場の等価回路表現を用い、媒質がランダムな特性を示す場合、音でも同様の局在化が起こることを示し、減衰特性についてシミュレーションを行い、その正当性を検証する。

１次元ダクト内の音の振る舞いは電気工学の等価回路を用いて表すことができる。本研究では、マイクロホン２本の観測信号を用いて騒音を打ち消すために２段のフィルタを導入した。さらに帯域信号を用いたシミュレーションによりその正当性を検証した。

騒音を逆位相の音で打ち消す能動騒音制御(ANC)について、特に１次元ダクトに着目して、回路理論に基づき新たな制御理論を提案し、純音に関するシミュレーションによりその正当性を検証した。

近年盛んに研究されている独立成分分析(ICA)では、対象外成分が存在すると分離性能は劣化する。本研究では、著者らがすでに提案しているHermiteモーメントと回転変換に基づくICA 手法について，対象外成分によりどのように推定精度が劣化するか検討し，シミュレーションを通じて確認する．

分布定数系の等価回路により1 次元音場を書き直すと、主・副音源を用いた音場の制御が原理的に可能となる。これによると騒音制御用のフィルタ設計も容易である。本研究では、とくに1 次元ダクトを具体例とするシミュレーションを報告するが、結果は良好であった。

時間差混合信号に対する複素ブラインド分離の分離行列から遅延時間・振幅比の推定を行い、音源位置・方位を推定する手法を提案する。さらに推定パラメータを用いて混合信号の逆フィルタを構成することにより、帯域歪みのない実信号としての分離・再生を行う。

複素ICA（Independent Component Analysis）に基づいて時間差混合信号を分離する際得られる分離行列から遅延時間・振幅比の推定が行える。これにより音源位置の推定が原理的に可能となる。本研究では、音声によるシミュレーションを通じて理論の正当性を検証する。

すでに1 次元ダクトの伝送線路モデルを基本にして1 次元ダクトの能動騒音制御（Active Noise Control：ANC）を考察している。本研究では、既発表の理論に基づいて行ったシミュレーションを中心に報告した。結果は良好であり、理論の正当性の一部が示せた。

様々な観点から研究がなされている能動騒音制御（Active Noise Control：ANC）について、とくに1 次元音場に対し等価回路に基づき考察した。本手法は、既存の手法に比べ、ANC の原理ならびに、制御回路の周波数特性が単純で把握しやすいという特長がある。

既発表の1次元ダクトの等価回路に基づくANCは、既存の手法に比べ、ANCの原理ならびに、制御回路の周波数伝達特性が単純で把握しやすく、フィルタの構成も容易である。本報告では、簡単なシミュレーションを通じて理論の検証を行ったので報告する。

音源位置あるいは方向の推定のために、複素ICAに基づいて時間差混合信号の分離行列から遅延時間・振幅比の推定を行う。さらに、提案した逆フィルタによって、実際の音声信号で帯域歪みのない分離・再生ができることを、シミュレーションを通じて示す。

時間差混合信号に対し、複素ICAによる分離行列から遅延時間・振幅比の推定を行えば、分離と同時に音源位置あるいは方向の推定が原理的に可能である。さらに推定パラメータを用いた逆フィルタにより、帯域歪みのない分離や再生ができることを確かめた。

本報告では、前報告に引き続き、中型の防風スクリーンを装着した場合、また、風車発電機の騒音計測に用いられる計測システムを用いた場合の低周波雑音成分が、それぞれマイクロホン近傍の風速情報を用いて推定可能かどうかを実験的に検討している。

前報に引き続き、観測信号から風雑音を除去する手法を提案する。とくに実用的な観点から風速と風雑音の関係は線形モデルで記述し、その関係の時間変化も考慮する。さらに実際の風雑音混入下における音響信号検出の問題に本手法を適用し、その有効性を議論する。

The purpose of this research is to realize a distaince measurement system for detecting the distance to a human being using audible sound for devices such as nursing care and welfare robots that require a wide range of distance detection from near to far. However, the conventional method has a problem that the measurement accuracy is degraded by the influence of the measurement system between the loudspeaker and the microphone and the observation noise in the real environment. In this study, (1) a theory for the measurement is developed under high noise conditions, which is not possible with the conventional phase-interference-based acoustic distance measurement method, (2) a theory for eliminating the influence of the measurement system between the loudspeaker and the microphone is developed, and (3) the acoustic distance measurement method to a specific object is investigated. The results suggest the effectiveness of the proposed method.

(1) We introduced L and U filters having not been well-known as morphological filters and the digital pulse transform (DPT) constructed of L and U filters. (2) We proposed a new method for detecting randomly generated signals buried in noise, in particular, a new method to detect saccades in human eye movement data. We attained error-free detection of the saccades by this method which has not been achieved by traditional linear filters. This method can be widely applied to detection problems in biomedical environmental data processing. (3) Agaist the commomly accepted theory that median filtering outperforms traditional linear filtering in preserving signal edges while reducing noise, we theoretically proved that moving average filter is a better smoother preserving signal edges than the median filter under a practical wide range condition.

The purpose of this study is to realize the distance measuring system for detecting the distance by using the audible sound, especially in a close range (e.g., where a nursing-care robot care for human). (1) Based on the previously proposed acoustic distance measurement method (ADM) using a standing wave (or phaseinterference), we improved the original ADM method, and we considered the properties and potentials of the improved ADM method. (2) We considered ADM method to measure the distance shorter than the minimum measurable distance which is defined by the bandwidth of transmitted signal. Furthermore, (3) we considered the principle method for the moving target though the original ADM method is not applicable to the moving target. Thus, it becomes possible to measure both the velocity and distance for the moving target. This implies that the distance between microphone and target is measurable eve in a case when both microphone and target are moving.

The purpose of this study is to realize a rear-sonar of automobile for detecting obstacles(or targets) behind the automobile, using the exhaust sound of the automobile emitted from the automobile muffler. (1) We considered the fundamental properties and potentials of the distance estimation method based on the phase interference of audible sound. (2) We developed several methods, so as to reduce the error due to geometrical positions of the loudspeaker, target and microphone and the error due to the effect of the observation noise and measurement system. Furthermore, (3) we proposed a distance measurement method to measure the distance shorter than the minimum measurable distance.

(1) We proposed a new mutual information criterion for deciding an optimal mixture distribution and an optimal threshold. (2) We proposed the usage of Renyi's alpha-entropy as a measure of width of marginal distribution and obtained very intuitive new and simple ICA methods. (3) We proposed an adaptive method for smoothing the periodogram for positive time-frequency representation. For impulse train time-frequency analysis, we proposed an inverse domain demodulation method of IPFM and an objective evaluation method for demodulators. Further, we proposed distance estimation method exploiting time-frequency feature of reflected signal from a moving object.

It is known that a wave cannot propagate long distance in random media and is confined in a small region when the field is excited. In physics it is called the "localization", the phenomenon which is closely related to the cut-off (non-propagating) state of waves in a certain broad frequency band. It is expected that this property can be utilized as a muffler to reduce the sound noise because the sound propagating in random media or random loadings would exhibit strong exponential decay along distance. On the other hand, a one-dimensional acoustic duct can be expressed by a transmission line model commonly used in the study of microwave waveguide. In this research, I focus on the sound decaying property of one-dimensional duct (random duct) and demonstrate the validity of the proposed theory through a computer simulation. The results I obtained in 2005 and 2006 are listed below.[2005 year] : As the first step of research, I expressed the sound field in duct by a transmission line model, on which I theoretically analyzed the ordinary one-dimensional duct and carried out computer simulation. Based on the circuit representation, I theoretically proposed the sound reduction property in a case when there exist random media in duct. To confirm the validity of the analysis, I performed the computer simulation and showed the existence of the strong exponential decay. These results have been published in a scientific journal and presented in domestic conference.[2006 year] : Following to the fundamental consideration on the sound reduction problem with randomness in 2005, I theoretically proposed the sound reduction property in a case when there exist random loadings in duct and the extensive simulation makes it clear that the proposed method also enables us to cause the strong exponential decay. These results have been presented in international conferences.

Sound source localization problem has been important in various actual fields of acoustics and many approaches have been developed e.g., MUSIC, etc. Most of them utilize only the lower order statistics of observations such as a covariance matrix and the orthogonality (or uncorrelation) between the objective signal and noise. Recently , Independent Component Analysis (ICA) becomes popular to solve a blind signal separation problem where independent sources are mixed and then measured by multiple sensors. The ICA of complex signals has been also proposed for a case when observations have time differences between source signals or are affected by system (room etc.).In this research, some trials of estimating the sound source positions or directions are proposed from various viewpoints. In particular, most of them are based on the statistical independency between source signals or between the source signal and an additional noise. The results we obtained in 2001 and 2002 are listed below.[2001 year]: As the first step of research, the correlation method is introduced to calculate the delay time between observations at each microphone and then to estimate the sound direction. In the second approach, MUSIC method and Kalman filter are employed, respectively, to estimate the direction of source and to detect the amplitude information for sound source localization.The third method is a trial of estimating the time delays or the distances between receivers and sound sources based on the complex ICA, especially for the mixed observations with time delays. Further, we can easily construct an inverse filter to separate the (over-all) real mixed signals using the parameters estimated by ICA, while the usual complex ICA required synthesizing all frequency components. Finally, through fundamental simulation, we demonstrate the high-quality blind separation as well as the validity for the estimation of the source parameters.[2002 year]: Following to the fundamental consideration on the sound source localization problem in 2001, the extensive simulation makes it clear that the proposed method also enables us to find the directions of arrival. These results have been presented in domestic or international conferences.

The main objective of this research is to make a plan and investigation for a joint research project on "Agricultural Resource Management taking into account Eco-system and Environment Problems"with IIASA (International Institute for Applied Systems Analysis) and Prof. R.H.Richardson of Texas University. To this end, 1. We co-organized a workshop on "Agricultural Resource Management taking into account Environment Problems" with the Society of Animal Production Systems, and observed Mt. Sanbe national park in which grass ground was recovered by pasturage. 2. We observed Saito Farm and others which succeeded in creating excellent environment with beautiful view by reclamation by pasturaged cattles, and co-organized a symposium on "Agriculture in Harmonious with Natural Eco-systems" with the Society of Animal Production Systems, and Open Network for Nature Conservation, In the symposium, we discussed and exchanged ideas on the basis of lectures i) Holistic Management in Agriculture in Harmonious with Eco-systems by Prof. R.H.Richardson of Texas University, and ii) Land Use Problems taking into account Earth Scale Environment by Dr. G.Fischer of IIASA. 3. We investigated the activities of Lake Biwa Research Institute, and had a meeting with a fisherman in Biwa Lake and a large scale farmer to discuss fishing and farming from the viewpoint of environmental management. Through the above activities, we clarified the theme and direction of our planning international joint research in the future.

本研究では、我々の生活と切り放すことができない音環境システムに特に着目し、その信号処理を目指して、多変量解析の分野で公知の主成分分析法と因子分析法の拡張を試みた。まず、特に音環境における確率評価量(例えば、Lx、Leq)間での双方向変換関係を見いだす目的で、公知の主成分分析法を理論的に拡張した。具体的には、線形関係および自乗誤差評価(すなわち、線形相関を利用)のみに基づく従来の手法とは異なり、非線形モデルの導入や、高次相関情報の誤差評価への利用などを行なった。さらに、公知の回帰分析法と対比させながら提案手法を実測データに適用し、良好な結果を得た。この成分分析の拡張手法には次のような特長がある。すなわち、1)実現象ごもつ複雑さや各変数間の非線形性に対応できる、2)多変数間の線形相関情報のみでなく非線形相関をも利用できる、3)従来の主成分分析法をスペシャルケースで含んでいる、4)2変数のみでなく、3変数以上の場合のもそのまま本手法が拡張できる、などである。ついで、音環境における複雑な多数遮音システムを、暗騒音に埋もれた出力観測のみから同時に同定したり、その出力応答を(騒音評価量とも関連し)揺らぎ分布全体において予測する目的で、特に実用的な立場から因子分析法の考え方を拡張した。具体的には、共通因子と独自因子に実体的メカニズムをまず反映させて、それぞれ騒音入力インテンシティ、各観測点での暗騒音インテンシティとして捉え、インテンシティスケールでの線形モデルに基づきシステムパラメータとして因子負荷量を推定した。この結果を用いて、暗騒音の影響がない場合、すなわち入力騒音のみに対する出力騒音の分布予測を行うことができる。さらに、本手法を実際の音環境データへ適用することによりその有効性を実験的にも確認できた。この因子分析の拡張手法は次のような特長を持っている。すなわち、1)複数システムを同時に扱える、2)入力および暗騒音が未知でも適用できる、3)計測において入出力間の同期をとる必要がない、などである。

本研究では、もともとの変動が有限な範囲内に留まっていたり、(飽和型非線形系の)計測器のダイナミック・レンジなどの存在により有限な範囲内の値を示す実際的な場合において、非定常入力を伴う系の応答解析を最終目標としている。今回は研究の初期的段階にあることから、特に飽和型非線形系における非定常変動観測時の雑音対策を、システムの物理メカニズムと直結するパワースケールとの関連で考究した。すなわち、有限レベル変動範囲のダイナミック・レンジをもつ観測データを用いて、任意分布型の外来雑音に汚された未知非定常信号(特に、パワースケールのような正の物理量)を、動的に推定してゆく推定アルゴリズムをベイズ定理に基づき開発した。具体的には、飽和の影響を受ける前の任意変動の観測値と任意非ガウス型変動を示す未知信号間の線形・非線形の各種相情報を階層的に反映した新たな信号復元法を、広義ディジタルフィルタの形で見い出した。更に本研究で得られた理論的結果を、シミュレーションデータや残響室内における暗騒音混入下の実音響データに適用し、その有効性を検証した。本研究で得られた理論の特長を列挙すると以下のとおりである。1)本手法は、外来雑音の混入とダイナミックレンジの存在に整合している。2)実システムが本質的にもつ非ガウス性、非線形性に対応できる。3)推定アルゴリズムが実用的である。すなわち、観測レンジ内では従来のベイズフィルタを形式的に採用し、観測レンジの上限、下限では確率密度の集中を簡易的に考慮している。4)スペシャルケースとして、ダイナミック・レンジが十分広い場合には、従来のベイズフィルタを理論的に包含している。本研究をもとに、今後、非定常な確率現象と計測における有限性を伴うあらゆる実分野への適用とその成果が期待できる(たとえば、機械振動、地盤振動、道路交通騒音などにおける非定常揺らぎの評価や解析など)。