JP7032961B2 - Controller design equipment and programs for acoustic signals - Google Patents

Description

本発明は、音響信号用の制御器を設計する技術に関し、特に、制御対象内の制御点に、スピーカ再生により所望の音響信号を提示するための制御器を設計する制御器設計装置及びプログラムに関する。 The present invention relates to a technique for designing a controller for an acoustic signal, and more particularly to a controller design device and a program for designing a controller for presenting a desired acoustic signal by speaker reproduction at a control point in a controlled object. ..

従来、音響システムにおいて、ある制御位置に所望の音響信号を提示する方法が知られている。図１９は、音響システムの例を説明する概略図である。この音響システムは、音場である制御対象Ｇｐの逆システムＨinv、５台のスピーカ１０１－１～１０１－５及び２台のマイクロホン１０２－１，１０２－２により構成される。 Conventionally, in an acoustic system, a method of presenting a desired acoustic signal at a certain control position is known. FIG. 19 is a schematic diagram illustrating an example of an acoustic system. This acoustic system is composed of a reverse system Hinv of the controlled object Gp, which is a sound field, five speakers 101-1 to 101-5, and two microphones 102-1 and 102-2.

マイクロホン１０２－１，１０２－２の位置が制御位置であり、スピーカ１０１－１～１０１－５は、２つの制御位置に所望の音響信号を提示するための音響信号を再生する音源である。マイクロホン１０２－１，１０２－２の位置を制御点といい、スピーカ１０１－１～１０１－５を２次音源という。 The positions of the microphones 102-1 and 102-2 are the control positions, and the speakers 101-1 to 101-5 are sound sources for reproducing the acoustic signals for presenting the desired acoustic signals to the two control positions. The positions of the microphones 102-1 and 102-2 are called control points, and the speakers 101-1 to 101-5 are called secondary sound sources.

制御対象Ｇｐは、制御点に所望の音響信号を提示するための制御が行われる音場であり、２次音源から制御点までの音響伝達関数によってモデル化される。図１９の例では、制御対象Ｇｐの音響伝達関数は、制御点の数（２）×２次音源の数（５）の音響伝達関数行列として表される。制御対象Ｇｐは、２次音源を入力端、制御点を出力端とする状態空間モデルとして表されることもある。 The controlled target Gp is a sound field in which control is performed to present a desired acoustic signal to the control point, and is modeled by an acoustic transfer function from the secondary sound source to the control point. In the example of FIG. 19, the acoustic transfer function of the controlled object Gp is represented as an acoustic transfer function matrix of the number of control points (2) × the number of secondary sound sources (5). The controlled object Gp may be represented as a state space model in which the secondary sound source is the input end and the control point is the output end.

状態空間モデルの一般形は、入力ベクトル（入力信号）をｕ、出力ベクトル（出力信号）をｙとして、以下の式にて表される。

｛Ａ，Ｂ，Ｃ，Ｄ｝は、それぞれシステムの係数行列（システム行列）、入力側の係数行列（入力行列）、出力側の係数行列（出力行列）、直達項の係数行列である。 The general form of the state space model is expressed by the following equation, where u is the input vector (input signal) and y is the output vector (output signal).

{A, B, C, D} are a system coefficient matrix (system matrix), an input side coefficient matrix (input matrix), an output side coefficient matrix (output matrix), and a direct term coefficient matrix, respectively.

図２０は、前記式（１）で表す状態空間モデルをブロック線図で表した図である。図２０において、Ａ，Ｂ，Ｃ，Ｄで表した各ブロックは、入力信号と係数行列との積演算が行われる乗算器であり、Ｚ^-1で表したブロックは、１サンプル分の遅延（サンプルシフト）が行われる遅延器である。 FIG. 20 is a block diagram representation of the state space model represented by the equation (1). In FIG. 20, each block represented by A, B, C, and D is a multiplier in which the product operation of the input signal and the coefficient matrix is performed, and the block represented by Z ^-1 is a delay for one sample (delay for one sample). It is a delayer that performs sample shift).

また、このブロック線図には、Ｂで表したブロックの出力信号とＡで表したブロックの出力信号との加算演算が行われる加算器が含まれる。さらに、Ｃで表したブロックの出力信号とＤで表したブロックの出力信号との加算演算が行われる加算器が含まれる。 Further, this block diagram includes an adder in which an addition operation of the output signal of the block represented by B and the output signal of the block represented by A is performed. Further, an adder for performing an addition operation between the output signal of the block represented by C and the output signal of the block represented by D is included.

図１９に戻って、制御点に所望の音響信号を提示するためには、制御対象Ｇｐの逆システムＨinvを適切に設計した上で、逆システムＨinvの入力端に所望の音響信号を入力すればよい。逆システムＨinvは、制御対象Ｇｐと逆の特性を持つシステムである。制御対象Ｇｐ及び逆システムＨinvを音響伝達関数行列及び状態空間モデルで記述する場合、その関係は以下の式にて表される。

Ｉは単位行列である。 Returning to FIG. 19, in order to present the desired acoustic signal to the control point, after appropriately designing the inverse system Hinv of the controlled object Gp, the desired acoustic signal is input to the input end of the inverse system Hinv. good. The reverse system Hinv is a system having characteristics opposite to those of the controlled target Gp. When the controlled object Gp and the inverse system Hinv are described by the acoustic transfer function matrix and the state space model, the relationship is expressed by the following equation.

I is an identity matrix.

図２１は、仮想音源から制御点に所望の音響信号を提示する音響システムの例を説明する概略図である。この音響システムは、２次音源とは別の位置にある音源であるスピーカ１０３から、所望の音響信号を再生することを模擬するシステムである。この音源を仮想音源という。仮想音源から制御点までの音響伝達関数及び状態空間モデルによってモデル化される音場を目標システムＧｔとする。 FIG. 21 is a schematic diagram illustrating an example of an acoustic system that presents a desired acoustic signal from a virtual sound source to a control point. This acoustic system is a system that simulates the reproduction of a desired acoustic signal from the speaker 103, which is a sound source located at a position different from the secondary sound source. This sound source is called a virtual sound source. The target system Gt is a sound field modeled by an acoustic transfer function from a virtual sound source to a control point and a state space model.

図２２は、目標システムＧｔの入力端に音響信号を入力する音響システムの例を説明する概略図である。図２２に示すように、制御対象Ｇｐの逆システムＨinvと目標システムＧｔとを直列に接続し、目標システムＧｔの入力端に音響信号を入力することにより、仮想音源からの音響信号の再生を模擬することができる。 FIG. 22 is a schematic diagram illustrating an example of an acoustic system in which an acoustic signal is input to an input end of the target system Gt. As shown in FIG. 22, the inverse system Hinv of the controlled target Gp and the target system Gt are connected in series, and the acoustic signal is input to the input end of the target system Gt to simulate the reproduction of the acoustic signal from the virtual sound source. can do.

図２３は、制御器Ｈを備えた音響システムの例を説明する概略図である。図２３に示すように、この音響システムは、図２２に示した逆システムＨinvの特性と目標システムＧｔの特性とを合わせた制御器Ｈを備えている。この音響システムでは、制御器Ｈに対し、逆システムＨinv及び目標システムＧｔの両者の処理を行わせることが可能である。 FIG. 23 is a schematic diagram illustrating an example of an acoustic system including the controller H. As shown in FIG. 23, this acoustic system includes a controller H that combines the characteristics of the inverse system Hinv shown in FIG. 22 with the characteristics of the target system Gt. In this acoustic system, it is possible to have the controller H process both the reverse system Hinv and the target system Gt.

制御器Ｈは、以下の式を満たすように設計される。

前記式（３）は、以下の式に変換される。これにより、仮想音源からの音響信号の再生が模擬される。

The controller H is designed to satisfy the following equation.

The above equation (3) is converted into the following equation. As a result, the reproduction of the acoustic signal from the virtual sound source is simulated.

このような音響システムにおいて、制御点を聴取者の両耳鼓膜位置に定めた再生法をトランスオーラル再生法という。トランスオーラル再生法は、ある音響信号に所定の信号処理を施し、スピーカから音響信号を再生することにより、そのスピーカが置かれていない別の位置から音響信号が再生されているかのように、聴取者に知覚させる技術である。 In such an acoustic system, a reproduction method in which a control point is set at the position of both eardrums of the listener is called a transoral reproduction method. In the transoral reproduction method, a certain acoustic signal is subjected to predetermined signal processing, and the acoustic signal is reproduced from a speaker, so that the acoustic signal is reproduced as if the acoustic signal is reproduced from another position where the speaker is not placed. It is a technique that makes people perceive.

図２４は、トランスオーラル再生法を実現する音響システムの例を説明する図である。聴取者１０４の両方の耳１０５－１，１０５－２における鼓膜位置が所望の音響信号を提示する制御位置であり、両方の耳１０５－１，１０５－２の位置は制御点である。音響信号を再生するスピーカ１０１－１～１０１－５の音源は２次音源である。音響伝達関数は、スピーカ１０１－１～１０１－５から両方の耳１０５－１，１０５－２の鼓膜位置までの伝達関数であり、特に、頭部伝達関数と呼ばれる。 FIG. 24 is a diagram illustrating an example of an acoustic system that realizes a transoral reproduction method. The eardrum positions in both ears 105-1 and 105-2 of the listener 104 are control positions that present the desired acoustic signal, and the positions of both ears 105-1 and 105-2 are control points. The sound sources of the speakers 101-1 to 101-5 that reproduce the acoustic signal are secondary sound sources. The acoustic transfer function is a transfer function from the speakers 101-1 to 101-5 to the eardrum positions of both ears 105-1 and 105-2, and is particularly called a head-related transfer function.

図２４の例では、音場である制御対象Ｇｐは、制御点の数（２）×２次音源の数（５）の頭部伝達関数行列、または５入力２出力の状態空間モデルとして表現される。目標システムＧｔも同様に、頭部伝達関数行列または状態空間モデルとして表現される。 In the example of FIG. 24, the controlled object Gp, which is a sound field, is represented as a head-related transfer function matrix of the number of control points (2) × the number of secondary sound sources (5), or a state space model with 5 inputs and 2 outputs. To. The target system Gt is also expressed as a head related transfer function matrix or a state space model.

図２２に示した音響システムと同様に、図２４に示す音響システムでは、制御対象Ｇｐの逆システムＨinvを適切に設計した上で、逆システムＨinvと目標システムＧｔとを直列に接続し、目標システムＧｔの入力端に音響信号を入力する。 Similar to the acoustic system shown in FIG. 22, in the acoustic system shown in FIG. 24, after appropriately designing the inverse system Hinv of the controlled target Gp, the inverse system Hinv and the target system Gt are connected in series to form the target system. An acoustic signal is input to the input end of Gt.

これにより、仮想音源からの音響信号の再生を模擬することができ、聴取者１０４に対し、仮想音源から音響信号が再生されているかのように、知覚させることができる。 As a result, the reproduction of the acoustic signal from the virtual sound source can be simulated, and the listener 104 can be perceived as if the acoustic signal is being reproduced from the virtual sound source.

図２５は、制御器Ｈを備えた、トランスオーラル再生法を実現する音響システムの例を説明する概略図である。図２５に示すように、この音響システムは、図２４に示した逆システムＨinvの特性と目標システムＧｔの特性とを合わせた制御器Ｈを備えている。 FIG. 25 is a schematic diagram illustrating an example of an acoustic system that realizes a transoral reproduction method and includes a controller H. As shown in FIG. 25, this acoustic system includes a controller H that combines the characteristics of the inverse system Hinv shown in FIG. 24 with the characteristics of the target system Gt.

トランスオーラル再生法においても、図２３に示した音響システムの場合と同様に、制御器Ｈに対し、逆システムＨinv及び目標システムＧｔの両者の処理を行わせることが可能である。 Also in the transoral reproduction method, as in the case of the acoustic system shown in FIG. 23, it is possible to cause the controller H to process both the reverse system Hinv and the target system Gt.

制御器Ｈは、前記式（３）を満たすように設計される。前述のとおり、前記式（３）は前記式（４）に変換されるから、仮想音源からの音響信号の再生を模擬することができる。 The controller H is designed to satisfy the above equation (3). As described above, since the equation (3) is converted into the equation (4), it is possible to simulate the reproduction of the acoustic signal from the virtual sound source.

従来、制御器Ｈを設計する手法として、制御対象Ｇｐ及び目標システムＧｔを状態空間でモデル化し、前記式（３）を満たす制御器Ｈも同様にモデル化する手法が提案されている（例えば、非特許文献１を参照）。 Conventionally, as a method for designing a controller H, a method has been proposed in which a controlled object Gp and a target system Gt are modeled in a state space, and a controller H satisfying the above equation (3) is also modeled in the same manner (for example). See Non-Patent Document 1).

この手法は、制御器ＨのＨ∞ノルムを評価関数として制御器Ｈのゲインを定量化し、その値が最小となるように、制御器Ｈを設計するものである。詳細については、非特許文献２，３、本件特許出願の出願人及び発明者によりなされた、本件特許出願時に未公開の特願２０１７－３５５１７号、及び、本件特許出願の出願人等及び発明者等によりなされた、本件特許出願時に未公開の特願２０１７－１７２０３５号に記載されている。 In this method, the gain of the controller H is quantified using the H∞ norm of the controller H as an evaluation function, and the controller H is designed so that the value is minimized. For details, see Non-Patent Documents 2 and 3, Japanese Patent Application No. 2017-35517, which was made by the applicant and inventor of the Patent Application, which was unpublished at the time of filing the Patent Application, and the applicant and the inventor of the Patent Application. It is described in Japanese Patent Application No. 2017-172035, which was not published at the time of filing the patent application.

制御器Ｈのゲインは、制御器Ｈの入力信号及び出力信号の大きさの比であり、前述の音響システムのような系の摂動及び外乱は、このゲインに応じて増幅される。したがって、このゲインを小さくすることができれば、摂動及び外乱の影響を小さく抑えることができる。 The gain of the controller H is the ratio of the magnitudes of the input signal and the output signal of the controller H, and the perturbation and disturbance of the system as in the above-mentioned acoustic system are amplified according to this gain. Therefore, if this gain can be reduced, the effects of perturbations and disturbances can be suppressed to a small level.

音響システムの摂動及び外乱は、例えば、スピーカ１０１－１～１０１－５の向きまたは位置がずれる状況、マイクロホン１０２－１，１０２－２の向きまたは位置がずれる状況により生じる。また、聴取者１０４が別の者に代わる状況、外部からノイズが混入する状況等によっても生じる。 Perturbations and disturbances in the acoustic system are caused, for example, by situations where the speakers 101-1 to 101-5 are misaligned or misaligned, and microphones 102-1 and 102-2 are misaligned or misaligned. It also occurs due to a situation in which the listener 104 is replaced by another person, a situation in which noise is mixed in from the outside, and the like.

松井他、“Design Method for Transaural Reproduction Controller Based on H-Infinity Norm Constraint”、Proceedings of the SICE Annual Conference 2016 Tsukuba Japan ISBN978-4-907764-50-0、公益社団法人計測自動制御学会、September 20-23,2016Matsui et al., “Design Method for Transaural Reproduction Controller Based on H-Infinity Norm Constraint”, Proceedings of the SICE Annual Conference 2016 Tsukuba Japan ISBN978-4-907764-50-0, Society of Instrument and Control Engineers, September 20-23 , 2016 松井他、“出力追従制御を応用したトランスオーラル再生制御器の設計”、日本音響学会誌73巻5号(2017）、p.281-290Matsui et al., "Design of Transoral Regeneration Controller Applying Output Tracking Control", Journal of Acoustical Society of Japan, Vol. 73, No. 5 (2017), p.281-290 松井他、“出力追従制御を応用したトランスオーラル再生制御器の緩和処理法”、2017年秋季日本音響学会研究発表会講演論文集CD-ROM、1-P-31、2017、p.729-730Matsui et al., "Mitigation processing method for transoral reproduction controller applying output tracking control", Proceedings of the 2017 Autumn Meeting of the Acoustical Society of Japan CD-ROM, 1-P-31, 2017, p.729-730

前述のとおり、制御対象Ｇｐ、目標システムＧｔ、及び前記式（３）を満たす制御器Ｈを、前記式（１）に示した状態空間モデルで表すことにより、制御器Ｈを設計することができる。この手法は、前記式（２）～（４）の関係が成立することを前提としている。 As described above, the controller H can be designed by expressing the controlled object Gp, the target system Gt, and the controller H satisfying the equation (3) by the state space model shown in the equation (1). .. This method is based on the premise that the relationships of the above equations (2) to (4) are established.

しかしながら、一般的な音響システムにおいては、様々な要因から、前記式（２）～（４）の関係が成立しなくなることが多い。例えば、スピーカ１０１－１～１０１－５の向きまたは位置、マイクロホン１０２－１，１０２－２の向きまたは位置が僅かにずれるだけでも、その間の音響伝達関数が変わるからである。 However, in a general acoustic system, the relationships of the above equations (2) to (4) often do not hold due to various factors. For example, even if the orientation or position of the speakers 101-1 to 101-5 and the orientation or position of the microphones 102-1 and 102-2 are slightly deviated, the acoustic transfer function between them changes.

また、図２５において、聴取者１０４が別の者に代わると、頭部及び耳１０５－１，１０５－２の形状が変わるため、頭部伝達関数も変わる。また、音場での反射または残響、２次音源以外の音源から発せられる音（雑音）は外乱となる。 Further, in FIG. 25, when the listener 104 is replaced by another person, the shapes of the head and the ears 105-1 and 105-2 change, so that the head-related transfer function also changes. Further, the reflection or reverberation in the sound field, and the sound (noise) emitted from a sound source other than the secondary sound source becomes a disturbance.

前記式（２）～（４）に従って厳密に設計された逆システムＨinv及び制御器Ｈは、これらの系の摂動及び外乱が勘案されておらず、変化に対して脆弱である。音響システムの場合、例えば図１９からわかるように、制御点では、２次音源から再生される音響信号の重ね合わせにより、観察される音響信号が合成される。 The inverse system Hinv and the controller H, which are strictly designed according to the above equations (2) to (4), are vulnerable to changes without taking into account the perturbations and disturbances of these systems. In the case of an acoustic system, for example, as can be seen from FIG. 19, at the control point, the observed acoustic signal is synthesized by superimposing the acoustic signal reproduced from the secondary sound source.

系の摂動または外乱によって音響伝達関数が変わり、その位相が反転する場合、加算されるはずの音響信号が相殺され、または相殺されるはずの音響信号が加算増幅される等の問題が生じることになる。特に、振幅の大きな音響信号を逆相で相殺している場合、系の摂動または外乱によって位相が反転した結果、大きな振幅同士が加算され、大きな誤差とノイズが発生することになる。 When the acoustic transfer function changes due to the perturbation or disturbance of the system and its phase is inverted, problems such as the acoustic signals that should be added are canceled out or the acoustic signals that should be added up are added and amplified. Become. In particular, when an acoustic signal having a large amplitude is canceled by the opposite phase, the phase is inverted due to the perturbation or disturbance of the system, and as a result, the large amplitudes are added to each other, resulting in a large error and noise.

一般に制御器Ｈは、そのゲインが小さいほど、系の摂動及び外乱に対して頑健となる。しかし、前記式（２）～（４）の関係が成立しなくなる場合があり、前述の非特許文献１の手法では、ゲインを十分に小さくすることができない場合があり得る。 In general, the smaller the gain of the controller H, the more robust it is to the perturbations and disturbances of the system. However, the relationship of the above equations (2) to (4) may not be established, and the above-mentioned method of Non-Patent Document 1 may not be able to sufficiently reduce the gain.

前述の非特許文献１の手法は、前記式（３）が満たされるように、仮想音源からの再生が制御点にて完全に再現されることを条件として、その条件の下でゲインを最小化している。この場合、前記式（３）を条件とすると、ゲインを十分に小さくする解が得られないことがあるという問題があった。 The method of Non-Patent Document 1 described above minimizes the gain under the condition that the reproduction from the virtual sound source is completely reproduced at the control point so that the equation (3) is satisfied. ing. In this case, if the above equation (3) is a condition, there is a problem that a solution that sufficiently reduces the gain may not be obtained.

そこで、本発明は前記課題を解決するためになされたものであり、その目的は、制御性能の低下をある程度許容することで、制御器のゲインを小さくすることが可能な制御器設計装置及びプログラムを提供することにある。 Therefore, the present invention has been made to solve the above problems, and an object thereof is a controller design device and a program capable of reducing the gain of the controller by allowing a certain degree of deterioration in control performance. Is to provide.

前記目的を達成するために、請求項１の制御器設計装置は、信号が観察される所定数の制御点と前記所定数の制御点へ信号をそれぞれ提示する所定数の２次音源とを含む制御対象についての逆システムであって、前記所定数の２次音源へ信号を出力する前記逆システムの特性と、前記所定数の制御点へ信号をそれぞれ提示する所定数の仮想音源を含む目標システムの特性とを合わせた制御器を設計する音響信号用の制御器設計装置において、前記制御器が、第１遅延器、仮制御器及び第２遅延器を備え、前記制御対象を、前記所定数の２次音源から提示される第１信号と前記所定数の制御点にて観察される第２信号との間の関係について、直達項を持たない状態空間モデルで表し、前記目標システムを、前記所定数の仮想音源から提示される第３信号と前記所定数の制御点にて観察される第４信号との間の関係について、直達項を持たない状態空間モデルで表した場合に、前記第１信号、前記第２信号、前記第３信号及び前記第４信号に基づいて、前記制御対象及び前記目標システムを同定する同定部と、前記第１信号及び前記第２信号に基づいて、前記制御対象から遅延時間を分離し、当該遅延時間を分離した分離制御対象を求めると共に、前記第３信号及び前記第４信号に基づいて、前記目標システムから遅延時間を分離し、当該遅延時間を分離した分離目標システムを求める遅延時間分離部と、入力信号を前記第３信号とし、出力信号を、前記仮制御器のゲイン、及び前記分離制御対象と前記分離目標システムとの間の出力追従誤差とし、前記ゲインと前記出力追従誤差との間のバランスを定めるパラメータをバランス係数として、前記入力信号を入力する前記仮制御器、前記仮制御器により出力された信号を入力する前記分離制御対象、前記仮制御器により出力された信号に前記バランス係数を乗算し、乗算結果を前記ゲインとして出力する演算器、前記入力信号を入力する前記分離目標システム、及び、前記分離制御対象により出力された信号から前記分離目標システムにより出力された信号を減算し、減算結果を前記出力追従誤差として出力する減算器からなる全体システムを構成するシステム構成部と、前記システム構成部により構成された前記全体システムのＨ∞ノルムが最小となるように、前記仮制御器のパラメータを決定するパラメータ決定部と、前記遅延時間分離部により分離された前記遅延時間に基づいて、前記第３信号を入力する前記第１遅延器の遅延時間を算出すると共に、前記第１信号を出力する前記第２遅延器の遅延時間を算出する遅延時間算出部と、を備え、前記全体システムを、直達項を持つ状態空間モデルで表し、前記全体システムの伝達関数をＨall（ｚ）、前記全体システムの係数行列をＡall、入力側の係数行列をＢall、出力側の係数行列をＣall及び直達項の係数行列をＤallとし、前記仮制御器を、直達項を持つ状態空間モデルで表し、前記仮制御器の係数行列をＡ0、入力側の係数行列をＢ0、出力側の係数行列をＣ0及び直達項の係数行列をＤ0とし、遅延時間を分離した前記制御対象である遅延分離制御対象を、直達項を持たない状態空間モデルで表し、前記遅延分離制御対象の係数行列をＡｗ＃、入力側の係数行列をＢｐ＃、出力側の係数行列をＣｗ＃とし、遅延時間を分離した前記目標システムである遅延分離目標システムを、直達項を持たない状態空間モデルで表し、前記遅延分離目標システムの係数行列を前記Ａｗ＃、入力側の係数行列をＢｔ＃、出力側の係数行列を前記Ｃｗ＃とし、前記バランス係数をαとして、以下の数式：

で表される場合に、前記パラメータ決定部が、前記全体システムＨallのＨ∞ノルムを最小化する条件の数式：

を線形行列不等式で表したときに、パラメータγが最小となるように、前記仮制御器の係数行列Ａ0、前記入力側の係数行列Ｂ0、前記出力側の係数行列Ｃ0及び前記直達項の係数行列Ｄ0を、前記パラメータとして決定する、ことを特徴とする。 In order to achieve the above object, the controller design device according to claim 1 includes a predetermined number of control points in which a signal is observed and a predetermined number of secondary sound sources that present signals to the predetermined number of control points. A target system that is a reverse system for a controlled object and includes the characteristics of the reverse system that outputs signals to the predetermined number of secondary sound sources and a predetermined number of virtual sound sources that present signals to the predetermined number of control points. In a controller design device for an acoustic signal that designs a controller that matches the characteristics of the above, the controller includes a first delay device, a temporary controller, and a second delay device, and the control target is the predetermined number. The relationship between the first signal presented from the secondary sound source of No. 1 and the second signal observed at the predetermined number of control points is represented by a state space model having no direct term, and the target system is described as described above. The relationship between the third signal presented by a predetermined number of virtual sound sources and the fourth signal observed at the predetermined number of control points is represented by a state space model having no direct term. The control unit based on the identification unit that identifies the control target and the target system based on the one signal, the second signal, the third signal, and the fourth signal, and the control based on the first signal and the second signal. The delay time was separated from the target, the separation control target with the delay time separated was obtained, and the delay time was separated from the target system based on the third signal and the fourth signal, and the delay time was separated. The delay time separation unit for obtaining the separation target system, the input signal is the third signal, and the output signal is the gain of the temporary controller and the output tracking error between the separation control target and the separation target system. The temporary controller for inputting the input signal, the separation control target for inputting the signal output by the temporary controller, and the temporary controller with a parameter for determining the balance between the gain and the output tracking error as a balance coefficient. From the arithmetic unit that multiplies the signal output by the controller and outputs the multiplication result as the gain, the separation target system that inputs the input signal, and the signal output by the separation control target. A system configuration unit consisting of a subtractor that subtracts the signal output by the separation target system and outputs the subtraction result as the output tracking error, and H∞ of the entire system configured by the system configuration unit. The parameter determination unit that determines the parameters of the temporary controller and the delay time separated by the delay time separation unit so that the norm is minimized. Based on this, a delay time calculation unit for calculating the delay time of the first delay device for inputting the third signal and calculating the delay time of the second delay device for outputting the first signal is provided. The whole system is represented by a state space model having a direct term, the transfer function of the whole system is Hall (z), the coefficient matrix of the whole system is Aall, the coefficient matrix of the input side is Ball, and the coefficient matrix of the output side is. Let the coefficient matrix of the Call and the direct term be Dall, the temporary controller is represented by a state space model having the direct term, the coefficient matrix of the temporary controller is A0, the coefficient matrix of the input side is B0, and the coefficient matrix of the output side. Is C0 and the coefficient matrix of the direct term is D0, the delay separation control target which is the control target separated by the delay time is represented by a state space model having no direct term, and the coefficient matrix of the delay separation control target is Aw #. , The coefficient matrix on the input side is Bp #, the coefficient matrix on the output side is Cw #, and the delay separation target system, which is the target system in which the delay time is separated, is represented by a state space model having no direct term, and the delay separation is performed. Let Aw # be the coefficient matrix of the target system, Bt # be the coefficient matrix on the input side, Cw # be the coefficient matrix on the output side, and α be the balance coefficient.

When represented by, the formula of the condition that the parameter determination unit minimizes the H∞ norm of the whole system Hall:

Is expressed by a linear matrix inequality, the coefficient matrix A0 of the temporary controller, the coefficient matrix B0 on the input side, the coefficient matrix C0 on the output side, and the coefficient matrix of the direct term so that the parameter γ is minimized. It is characterized in that D0 is determined as the parameter.

さらに、請求項２のプログラムは、コンピュータを、請求項１に記載の制御器設計装置として機能させることを特徴とする。 Further, the program of claim 2 is characterized in that the computer functions as the controller design device according to claim 1.

以上のように、本発明によれば、制御性能の低下をある程度許容することで、制御器のゲインを小さくすることが可能となる。したがって、従来技術では制御器のゲインを小さくすることができない場合に、本発明では、制御器の設計に多少の誤差は生じるが、ゲインを小さくすることができ、結果として、系の摂動及び外乱に対して一層頑健な制御器を設計することが可能となる。 As described above, according to the present invention, it is possible to reduce the gain of the controller by allowing a certain degree of deterioration in control performance. Therefore, when the gain of the controller cannot be reduced by the prior art, the present invention can reduce the gain, although some errors occur in the design of the controller, and as a result, the perturbation and disturbance of the system. It is possible to design a more robust controller for the above.

本発明の実施形態による制御器設計装置のハードウェア構成を示す概略図である。It is a schematic diagram which shows the hardware composition of the controller design apparatus by embodiment of this invention. 制御器を含む音響システムの構成例を示す図である。It is a figure which shows the structural example of the acoustic system including a controller. 制御対象Ｇｐを、直達項を持たない状態空間モデルで表した場合のブロック線図である。It is a block diagram when the controlled object Gp is represented by the state space model which does not have a direct term. 制御部の機能構成例を示すブロック図である。It is a block diagram which shows the functional structure example of a control part. 制御器Ｈの構成例を示すブロック図である。It is a block diagram which shows the structural example of the controller H. 制御器Ｈ0を、直達項を持つ状態空間モデルで表した場合のブロック線図である。It is a block diagram when the controller H0 is represented by the state space model having a direct term. 制御部の処理例を示すフローチャートである。It is a flowchart which shows the processing example of a control part. 制御対象Ｇｐ及び目標システムＧｔを結合した結合システムＧｗを説明する図である。It is a figure explaining the coupling system Gw which combined the control target Gp and the target system Gt. 遅延時間を分離した制御対象Ｇｐ＃を説明する図である。It is a figure explaining the controlled object Gp # which separated the delay time. 遅延時間を分離した目標システムＧｔ＃を説明する図である。It is a figure explaining the target system Gt # which separated the delay time. 全体システムＨallの構成例を示すブロック線図である。It is a block diagram which shows the configuration example of the whole system Hall. 図１１を簡単化したブロック線図である。FIG. 11 is a simplified block diagram of FIG. 制御部の他の例（第１構成例）を示すブロック図である。It is a block diagram which shows the other example (first configuration example) of a control part. 第１構成例の制御部により設計された制御器の適用例１を示すブロック図である。It is a block diagram which shows the application example 1 of the controller designed by the control part of 1st configuration example. 第１構成例の制御部により設計された制御器の適用例２を示すブロック図である。It is a block diagram which shows the application example 2 of the controller designed by the control part of 1st configuration example. 制御部の他の例（第２構成例）を示すブロック図である。It is a block diagram which shows another example (second composition example) of a control part. 第２構成例の制御部により設計された制御器の適用例１を示すブロック図である。It is a block diagram which shows the application example 1 of the controller designed by the control part of the 2nd configuration example. 第２構成例の制御部により設計された制御器の適用例２を示すブロック図である。It is a block diagram which shows the application example 2 of the controller designed by the control part of the 2nd configuration example. 音響システムの例を説明する概略図である。It is a schematic diagram explaining an example of an acoustic system. 式（１）で表す状態空間モデルをブロック線図で表した図であるIt is a figure which represented the state space model represented by the equation (1) by a block diagram. 仮想音源から制御点に所望の音響信号を提示する音響システムの例を説明する概略図である。It is a schematic diagram explaining an example of an acoustic system which presents a desired acoustic signal from a virtual sound source to a control point. 目標システムＧｔの入力端に音響信号を入力する音響システムの例を説明する概略図である。It is the schematic explaining the example of the acoustic system which inputs an acoustic signal to the input end of the target system Gt. 制御器Ｈを備えた音響システムの例を説明する概略図である。It is a schematic diagram explaining the example of the acoustic system provided with the controller H. トランスオーラル再生法を実現する音響システムの例を説明する図である。It is a figure explaining the example of the acoustic system which realizes a transoral reproduction method. 制御器Ｈを備えた、トランスオーラル再生法を実現する音響システムの例を説明する概略図である。It is a schematic diagram explaining the example of the acoustic system which realizes the transoral reproduction method equipped with a controller H.

以下、本発明を実施するための形態について図面を用いて詳細に説明する。本発明は、所定の全体システムのＨ∞ノルムを評価関数として、制御器Ｈのゲイン及び誤差を定量化し、制御器Ｈの誤差をある程度許容するように、Ｈ∞ノルムを最小化する制御器Ｈのパラメータを決定する。 Hereinafter, embodiments for carrying out the present invention will be described in detail with reference to the drawings. The present invention uses the H∞ norm of a predetermined overall system as an evaluation function to quantify the gain and error of the controller H, and minimizes the H∞ norm so as to allow the error of the controller H to some extent. Determine the parameters of.

これにより、従来技術では制御器Ｈのゲインを十分に小さくすることができない場合に、制御器Ｈの制御性能の低下をある程度許容することにより、制御器Ｈの設計に多少の誤差は生じるが、ゲインを小さくすることができる。 As a result, when the gain of the controller H cannot be sufficiently reduced by the prior art, some errors occur in the design of the controller H by allowing the control performance of the controller H to deteriorate to some extent. The gain can be reduced.

〔ハードウェア構成〕
まず、本発明の実施形態による制御器設計装置のハードウェア構成について説明する。図１は、本発明の実施形態による制御器設計装置のハードウェア構成を示す概略図である。この制御器設計装置１は、制御器Ｈを設計する装置である。制御器設計装置１は、ＣＰＵ１１と、プログラム及びテーブル等を記憶するＲＯＭ及びＲＡＭからなる記憶部１２と、アプリケーションのプログラム、テーブル及びデータ等を記憶する記憶装置（例えばハードディスク装置）１３と、当該制御器設計装置１のオペレータによるキーボード及びマウス等の操作に伴い、所定のデータを入力制御する操作／入力部１４と、オペレータに対しデータ入力操作等を促すための画面情報を表示器に出力する表示出力インタフェース部１５と、インターネット等のネットワークを介してプログラム及びデータの送受信を行う通信部１６と、を備えて構成され、これらの構成部はシステムバス１７を介して相互に接続される。 [Hardware configuration]
First, the hardware configuration of the controller design device according to the embodiment of the present invention will be described. FIG. 1 is a schematic diagram showing a hardware configuration of a controller design device according to an embodiment of the present invention. The controller design device 1 is a device for designing the controller H. The controller design device 1 includes a CPU 11, a storage unit 12 including a ROM and a RAM for storing programs, tables, and the like, a storage device (for example, a hard disk device) 13 for storing application programs, tables, data, and the like, and the control thereof. An operation / input unit 14 for inputting and controlling predetermined data and a display for outputting screen information for prompting the operator to perform a data input operation, etc., in accordance with an operation of a keyboard, a mouse, etc. by the operator of the device design device 1. It is configured to include an output interface unit 15 and a communication unit 16 that transmits / receives programs and data via a network such as the Internet, and these components are connected to each other via a system bus 17.

記憶装置１３には、制御器設計装置１の基本的な機能を提供するＯＳ（オペレーティングシステム）プログラム、制御器設計プログラム、及び、制御器設計プログラムにて使用する各種テーブル及びデータ等が記憶されている。 The storage device 13 stores an OS (operating system) program that provides the basic functions of the controller design device 1, a controller design program, and various tables and data used in the controller design program. There is.

制御器設計プログラムは、制御対象Ｇｐの逆システムＨinvの特性と、所望の特性を付与する目標システムＧｔの特性とを合わせた制御器Ｈを設計するためのプログラムである。具体的には、制御器設計プログラムは、制御対象Ｇｐを、直達項を持たない状態空間モデルで表現することを前提に、所定の全体システムのＨ∞ノルムを評価関数として、制御器Ｈのゲイン及び誤差を定量化し、制御器Ｈの誤差をある程度許容するように、Ｈ∞ノルムを最小化する制御器Ｈのパラメータを決定する。 The controller design program is a program for designing a controller H that combines the characteristics of the inverse system Hinv of the controlled object Gp and the characteristics of the target system Gt that imparts desired characteristics. Specifically, the controller design program assumes that the controlled object Gp is represented by a state-space model having no direct term, and uses the H∞ norm of a predetermined overall system as an evaluation function to obtain the gain of the controller H. And the parameters of the controller H that minimize the H∞ norm are determined so that the error is quantified and the error of the controller H is tolerated to some extent.

尚、制御器設計プログラムは、当該制御器設計装置１が処理を行うときに、ＣＰＵ１１により記憶装置１３から記憶部１２のＲＡＭに読み出されて実行される。また、各種テーブル及びデータは、制御器設計プログラムの実行に伴い生成され、ＣＰＵ１１によって記憶部１２のＲＡＭから記憶装置１３へ書き込まれ、また、制御器設計プログラムの実行に伴い、ＣＰＵ１１によって記憶装置１３から記憶部１２のＲＡＭに読み出される。 The controller design program is read from the storage device 13 to the RAM of the storage unit 12 by the CPU 11 and executed when the controller design device 1 performs processing. Further, various tables and data are generated by the execution of the controller design program, are written from the RAM of the storage unit 12 to the storage device 13 by the CPU 11, and are written by the CPU 11 by the CPU 11 as the controller design program is executed. Is read into the RAM of the storage unit 12.

ここで、ＯＳプログラムは、ＣＰＵ１１により実行され、制御器設計装置１の基本的な機能として、記憶部１２、記憶装置１３、操作／入力部１４、表示出力インタフェース部１５及び通信部１６を管理する。そして、このＯＳプログラムがＣＰＵ１１によって実行された状態で、前述の制御器設計プログラムが実行される。 Here, the OS program is executed by the CPU 11 and manages the storage unit 12, the storage device 13, the operation / input unit 14, the display output interface unit 15, and the communication unit 16 as basic functions of the controller design device 1. .. Then, the controller design program described above is executed while this OS program is executed by the CPU 11.

制御部１０は、ＣＰＵ１１及び記憶部１２により構成され、ＣＰＵ１１が記憶装置１３に記憶された制御器設計プログラムを記憶部１２に読み出して実行することにより、制御器設計装置１全体を統括制御する。図１は、制御器設計プログラムが記憶装置１３から記憶部１２に読み出された状態を示している。このように、制御器設計装置１は、図１に示したハードウェア構成により、制御部１０が制御器設計プログラムに従って各種処理を行う。 The control unit 10 is composed of a CPU 11 and a storage unit 12, and the CPU 11 collectively controls the entire controller design device 1 by reading the controller design program stored in the storage device 13 into the storage unit 12 and executing the control unit 10. FIG. 1 shows a state in which the controller design program is read from the storage device 13 to the storage unit 12. As described above, in the controller design device 1, the control unit 10 performs various processes according to the controller design program according to the hardware configuration shown in FIG.

〔制御器Ｈ〕
次に、図１に示した制御器設計プログラムにより設計される制御器Ｈについて説明する。図２は、制御器Ｈを含む音響システムの構成例を示す図である。この音響システムは、制御器Ｈ、及び音場である制御対象Ｇｐにより構成される。 [Control H]
Next, the controller H designed by the controller design program shown in FIG. 1 will be described. FIG. 2 is a diagram showing a configuration example of an acoustic system including a controller H. This acoustic system is composed of a controller H and a controlled object Gp which is a sound field.

制御器Ｈは、制御対象Ｇｐの逆システムＨinvの特性と、所望の特性を付与する目標システムＧｔの特性とを備えている。制御器Ｈは、ｍ２個の音響信号（仮想音源信号）u2-1～u2-m2を入力し、ｍ１個の音響信号（２次音源信号）u1-1～u1-m1を生成し、ｍ１個の音響信号u1-1～u1-m1を２次音源であるスピーカ１０１－１～１０１－ｍ１へ出力する。ｍ２個の音響信号u2-1～u2-m2は、後述する仮想音源であるスピーカ１０３－１～１０３－ｍ２が出力する信号に相当する。 The controller H has the characteristics of the inverse system Hinv of the controlled target Gp and the characteristics of the target system Gt that imparts desired characteristics. The controller H inputs m2 acoustic signals (virtual sound source signals) u2-1 to u2-m2, generates m1 acoustic signals (secondary sound source signals) u1-1 to u1-m1, and m1 pieces. The acoustic signals u1-1 to u1-m1 are output to the speakers 101-1 to 101-m1 which are secondary sound sources. The m2 acoustic signals u2-1 to u2-m2 correspond to the signals output by the speakers 103-1 to 103-m2, which are virtual sound sources described later.

制御対象Ｇｐには、ｍ１個のスピーカ１０１－１～１０１－ｍ１が配置され、聴取者１０４の耳１０５－１，１０５－２の位置が制御点である。この制御点において音響信号が観察される。制御対象Ｇｐの入力点であるスピーカ１０１－１～１０１－ｍ１の位置が提示点であり、当該提示点から制御点へ音響信号が提示される。制御対象Ｇｐは、制御器Ｈからｍ１個の音響信号u1-1～u1-m1を入力し、スピーカ１０１－１～１０１－ｍ１から耳１０５－１，１０５－２へ音響信号u1-1～u1-m1を提示する。これにより、耳１０５－１，１０５－２にて音響信号y1,y2が観察される。 M1 speakers 101-1 to 101-m1 are arranged in the controlled target Gp, and the positions of the ears 105-1 and 105-2 of the listener 104 are control points. An acoustic signal is observed at this control point. The positions of the speakers 101-1 to 101-m1 which are the input points of the control target Gp are the presentation points, and the acoustic signal is presented from the presentation points to the control points. The control target Gp inputs m1 acoustic signals u1-1 to u1-m1 from the controller H, and the acoustic signals u1-1 to u1 are input from the speakers 101-1 to 101-m1 to the ears 105-1 and 105-2. Present -m1. As a result, the acoustic signals y1 and y2 are observed in the ears 105-1 and 105-2.

目標システムＧｔには、仮想音源であるｍ２個のスピーカ１０３－１～１０３－ｍ２が配置される。スピーカ１０３－１～１０３－ｍ２からのｍ２個の音響信号u2-1～u2-m2により、聴取者１０４の耳１０５－１，１０５－２の位置である制御点において音響信号が観察される。目標システムＧｔは、所望の音響信号を再生することを模擬するシステムであり、ｍ２個の音響信号u2-1～u2-m2を入力し、音響信号u2-1～u2-m2に対し、耳１０５－１，１０５－２の位置である制御点にて実現したい所望の音響特性が観察される。 In the target system Gt, m2 speakers 103-1 to 103-m2, which are virtual sound sources, are arranged. The m2 acoustic signals u2-1 to u2-m2 from the speakers 103-1 to 103-m2 observe the acoustic signals at the control points located at the ears 105-1 and 105-2 of the listener 104. The target system Gt is a system that simulates the reproduction of a desired acoustic signal, inputs m2 acoustic signals u2-1 to u2-m2, and ears 105 for the acoustic signals u2-1 to u2-m2. The desired acoustic characteristics desired to be achieved are observed at the control points at positions -1 and 105-2.

このように、制御器Ｈは、制御対象Ｇｐの逆システムＨinvの特性と、所望の特性を付与する目標システムＧｔの特性とを備え、逆システムＨinvが請け負う処理と、目標システムＧｔが請け負う処理とを一括して行う。 As described above, the controller H has the characteristics of the inverse system Hinv of the controlled target Gp and the characteristics of the target system Gt that imparts the desired characteristics, and the processing undertaken by the inverse system Hinv and the processing undertaken by the target system Gt. Do all at once.

〔制御対象Ｇｐ〕
本発明の実施形態では、制御対象Ｇｐは、直達項を持たない状態空間モデルとして表現する。 [Control target Gp]
In the embodiment of the present invention, the controlled object Gp is expressed as a state-space model having no direct term.

図３は、制御対象Ｇｐを、直達項を持たない状態空間モデルで表した場合のブロック線図であり、本発明の実施形態にて想定するモデルである。図２０に示した直達項を持つ状態空間モデルとこの図３の状態空間モデルとを比較すると、図２０では直達項の係数行列Ｄが存在するのに対し、図３では直達項の係数行列Ｄが存在しない点で相違する。その他は同じである。 FIG. 3 is a block diagram when the controlled object Gp is represented by a state-space model having no direct term, and is a model assumed in the embodiment of the present invention. Comparing the state-space model with the direct term shown in FIG. 20 with the state-space model of FIG. 3, the coefficient matrix D of the direct term exists in FIG. 20, whereas the coefficient matrix D of the direct term exists in FIG. The difference is that there is no. Others are the same.

直達項を持たない状態空間モデルの一般形は、以下の式にて表される。

前記式（５）は、直達項を持つ状態空間モデルの前記式（１）と異なり、直達項の係数行列Ｄと入力信号ｕの乗算結果Ｄｕ（ｋ）が存在しない。 The general form of the state-space model without a direct term is expressed by the following equation.

Unlike the equation (1) of the state space model having the direct term, the equation (5) does not have the multiplication result Du (k) of the coefficient matrix D of the direct term and the input signal u.

図３に示すように、入力信号ｕと出力信号ｙとの間に、１サンプル時間の遅れを生じさせる遅延器（Ｚ^-1）が存在する。このため、入力信号ｕのサンプル時間と出力信号ｙのサンプル時間とが同じではなく、この状態空間モデルは、１サンプル時間の遅延を有することとなる。 As shown in FIG. 3, there is a delay device (Z ^-1 ) that causes a delay of one sample time between the input signal u and the output signal y. Therefore, the sample time of the input signal u and the sample time of the output signal y are not the same, and this state-space model has a delay of one sample time.

〔制御器設計装置１の制御部１０〕
次に、図１に示した制御器設計装置１の制御部１０について説明する。図４は、制御部１０の機能構成例を示すブロック図であり、当該制御部１０が制御器設計プログラムの処理を実行する際の機能構成を示している。この制御部１０は、同定部２０、遅延時間分離部２１、システム構成部２２、パラメータ決定部２３及び遅延時間補償部（遅延時間算出部）２４を備えている。 [Control unit 10 of controller design device 1]
Next, the control unit 10 of the controller design device 1 shown in FIG. 1 will be described. FIG. 4 is a block diagram showing an example of the functional configuration of the control unit 10, and shows the functional configuration when the control unit 10 executes the processing of the controller design program. The control unit 10 includes an identification unit 20, a delay time separation unit 21, a system configuration unit 22, a parameter determination unit 23, and a delay time compensation unit (delay time calculation unit) 24.

図２に示した二次音源であるｍ１個のスピーカ１０１－１～１０１－ｍ１から出力される音響信号u1-1～u1-m1、仮想音源であるｍ２個のスピーカ１０３－１～１０３－ｍ２から出力される音響信号u2-1～u2-m2、及び制御点の耳１０５－１，１０５－２に提示される音響信号y1,y2が測定され、制御部１０は、これらの信号を測定信号として入力する。 The acoustic signals u1-1 to u1-m1 output from m1 speakers 101-1 to 101-m1, which are secondary sound sources shown in FIG. 2, and m2 speakers 103-1 to 103-m2, which are virtual sound sources. The acoustic signals u2-1 to u2-m2 output from the speaker and the acoustic signals y1 and y2 presented to the ear 105-1 and 105-2 of the control point are measured, and the control unit 10 measures these signals. Enter as.

制御部１０は、これらの測定信号である音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に基づいて、制御器設計プログラムの処理にて制御器Ｈを設計することで、制御器Ｈのパラメータを決定して出力する。制御器Ｈのパラメータは、制御器（仮制御器）Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0、遅延器（第１遅延器）Ｒｄの遅延時間Ｒd1～Ｒdm2、及び遅延器（第２遅延器）Ｒｃの遅延時間Ｒc1～Ｒcm1である。これらのパラメータの詳細については後述する。 The control unit 10 designs the controller H by processing the controller design program based on the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 which are the measurement signals. As a result, the parameters of the controller H are determined and output. The parameters of the controller H are the coefficient matrices A0, B0, C0, D0 of the controller (temporary controller) H0, the delay times Rd1 to Rdm2 of the delay device (first delay device) Rd, and the delay device (second delay device). ) The delay time of Rc is Rc1 to Rcm1. Details of these parameters will be described later.

制御部１０により決定されたこれらのパラメータを用いて構成される制御器Ｈは、図２に示したとおり、音響信号u2-1～u2-m2を入力し、音響信号u1-1～u1-m1を出力する。つまり、これらのパラメータは、入力信号である音響信号u2-1～u2-m2から出力信号である音響信号u1-1～u1-m1を生成するためのデータである。 As shown in FIG. 2, the controller H configured by using these parameters determined by the control unit 10 inputs the acoustic signals u2-1 to u2-m2 and inputs the acoustic signals u1-1 to u1-m1. Is output. That is, these parameters are data for generating acoustic signals u1-1 to u1-m1 which are output signals from the acoustic signals u2-1 to u2-m2 which are input signals.

図５は、制御器Ｈの構成例を示すブロック図であり、図６は、制御器Ｈ0を、直達項を持つ状態空間モデルで表した場合のブロック線図である。図５を参照して、この制御器Ｈは、遅延器Ｒｄ、制御器Ｈ0及び遅延器Ｒｃを備えている。 FIG. 5 is a block diagram showing a configuration example of the controller H, and FIG. 6 is a block diagram when the controller H0 is represented by a state space model having a direct term. With reference to FIG. 5, the controller H includes a delayer Rd, a controller H0 and a delayer Rc.

遅延器Ｒｄは、ｍ２個の遅延器３０－１，３０－２，・・・，３０－ｍ２により構成される。遅延器３０－１は、音響信号u2-1を入力し、音響信号u2-1を遅延時間Ｒd1だけ遅延させ、遅延後の音響信号u2-1を制御器Ｈ0に出力する。遅延器３０－２は、音響信号u2-2を入力し、音響信号u2-2を遅延時間Ｒd2だけ遅延させ、遅延後の音響信号u2-2を制御器Ｈ0に出力する。同様に、遅延器３０－ｍ２は、音響信号u2-m2を入力し、音響信号u2-m2を遅延時間Ｒdm2だけ遅延させ、遅延後の音響信号u2-m2を制御器Ｈ0に出力する。 The delay device Rd is composed of m2 delay devices 30-1, 30-2, ..., 30-m2. The delay device 30-1 inputs the acoustic signal u2-1, delays the acoustic signal u2-1 by the delay time Rd1, and outputs the delayed acoustic signal u2-1 to the controller H0. The delay device 30-2 inputs the acoustic signal u2-2, delays the acoustic signal u2-2 by the delay time Rd2, and outputs the delayed acoustic signal u2-2 to the controller H0. Similarly, the delay device 30-m2 inputs the acoustic signal u2-m2, delays the acoustic signal u2-m2 by the delay time Rdm2, and outputs the delayed acoustic signal u2-m2 to the controller H0.

制御器Ｈ0は、遅延器Ｒｄから遅延後の音響信号u2-1，u2-2～u2-m2を入力し、図６のブロック線図に示す処理（前記式（１）の処理に相当、後述する式（１１）の処理）を行う。そして、制御器Ｈ0は、ｍ１個の遅延前の音響信号u1-1,u1-2～u1-m1を生成して遅延器Ｒｃに出力する。 The controller H0 inputs the delayed acoustic signals u2-1, u2-2 to u2-m2 from the delay device Rd, and processes shown in the block diagram of FIG. 6 (corresponding to the process of the above equation (1), which will be described later. The processing of the equation (11) is performed. Then, the controller H0 generates m1 pre-delayed acoustic signals u1-1, u1-2 to u1-m1 and outputs them to the delay device Rc.

図６を参照して、制御器Ｈ0は、乗算器４０，４３，４４，４５、遅延器４２及び加算器４１，４６を備えている。乗算器４０は、入力信号ｕである音響信号u2-1～u2-m2（u2）を入力し、音響信号u2に係数行列Ｂ0を乗算し、乗算結果を加算器４１に出力する。加算器４１は、乗算器４０から乗算結果を入力すると共に、乗算器４３から乗算結果を入力し、両方の乗算結果を加算し、加算結果を遅延器４２に出力する。 With reference to FIG. 6, controller H0 includes multipliers 40, 43, 44, 45, delayers 42 and adders 41, 46. The multiplier 40 inputs acoustic signals u2-1 to u2-m2 (u2), which are input signals u, multiplies the acoustic signal u2 by the coefficient matrix B0, and outputs the multiplication result to the adder 41. The adder 41 inputs the multiplication result from the multiplier 40, inputs the multiplication result from the multiplier 43, adds both multiplication results, and outputs the addition result to the delayer 42.

遅延器４２は、加算器４１から加算結果を入力し、加算結果を１サンプル分遅延させ、１サンプル遅延後の加算結果ｘを乗算器４３，４４に出力する。乗算器４３は、遅延器４２から加算結果ｘを入力し、加算結果ｘに係数行列Ａ0を乗算し、乗算結果を加算器４１に出力する。 The delay device 42 inputs the addition result from the adder 41, delays the addition result by one sample, and outputs the addition result x after the delay of one sample to the multipliers 43 and 44. The multiplier 43 inputs the addition result x from the delay device 42, multiplies the addition result x by the coefficient matrix A0, and outputs the multiplication result to the adder 41.

乗算器４４は、遅延器４２から加算結果ｘを入力し、加算結果ｘに係数行列Ｃ0を乗算し、乗算結果を加算器４６に出力する。乗算器４５は、音響信号u2を入力し、音響信号u2に係数行列Ｄ0を乗算し、乗算結果を加算器４６に出力する。 The multiplier 44 inputs the addition result x from the delay device 42, multiplies the addition result x by the coefficient matrix C0, and outputs the multiplication result to the adder 46. The multiplier 45 inputs the acoustic signal u2, multiplies the acoustic signal u2 by the coefficient matrix D0, and outputs the multiplication result to the adder 46.

加算器４６は、乗算器４４から乗算結果を入力すると共に、乗算器４５から乗算結果を入力し、両方の乗算結果を加算し、加算結果を遅延前の音響信号u1-1,u1-2～u1-m1（u1）である出力信号ｙとして遅延器Ｒｃに出力する。 The adder 46 inputs the multiplication result from the multiplier 44, inputs the multiplication result from the multiplier 45, adds both multiplication results, and outputs the addition result to the acoustic signals u1-1, u1-2 to before delay. It is output to the delay device Rc as an output signal y which is u1-m1 (u1).

図５に戻って、遅延器Ｒｃは、ｍ１個の遅延器３１－１，３１－２，・・・，３１－ｍ１により構成される。遅延器３１－１は、制御器Ｈ0から遅延前の音響信号u1-1を入力し、遅延前の音響信号u1-1を遅延時間Ｒc1だけ遅延させ、音響信号u1-1をスピーカ１０１－１へ出力する。遅延器３１－２は、制御器Ｈ0から遅延前の音響信号u1-2を入力し、遅延前の音響信号u1-2を遅延時間Ｒc2だけ遅延させ、音響信号u1-2をスピーカ１０１－２へ出力する。同様に、遅延器３１－ｍ１は、制御器Ｈ0から遅延前の音響信号u1-m1を入力し、遅延前の音響信号u1-m1を遅延時間Ｒcm1だけ遅延させ、音響信号u1-m1をスピーカ１０１－ｍ１へ出力する。 Returning to FIG. 5, the delay device Rc is composed of m1 delay devices 31-1, 31-2, ..., 31-m1. The delayr 31-1 inputs the acoustic signal u1-1 before the delay from the controller H0, delays the acoustic signal u1-1 before the delay by the delay time Rc1, and sends the acoustic signal u1-1 to the speaker 101-1. Output. The delay device 31-2 inputs the acoustic signal u1-2 before the delay from the controller H0, delays the acoustic signal u1-2 before the delay by the delay time Rc2, and sends the acoustic signal u1-2 to the speaker 101-2. Output. Similarly, the delay device 31-m1 inputs the acoustic signal u1-m1 before the delay from the controller H0, delays the acoustic signal u1-m1 before the delay by the delay time Rcm1, and makes the acoustic signal u1-m1 the speaker 101. Output to -m1.

図６に示した制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0は、制御部１０のパラメータ決定部２３により決定され、図５に示した遅延器Ｒｄの遅延時間Ｒd1～Ｒdm2及び遅延器Ｒｃの遅延時間Ｒc1～Ｒcm1は、制御部１０の遅延時間補償部２４により決定される。 The coefficient matrices A0, B0, C0, D0 of the controller H0 shown in FIG. 6 are determined by the parameter determination unit 23 of the control unit 10, and the delay times Rd1 to Rdm2 and the delayers Rc of the delay devices Rd shown in FIG. 5 are determined. The delay times Rc1 to Rcm1 are determined by the delay time compensating unit 24 of the control unit 10.

〔制御部１０の処理〕
図７は、制御部１０の処理例を示すフローチャートである。以下、制御部１０の処理例について詳細に説明する。前述のとおり、制御部１０は、音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に基づいて、制御器設計プログラムの処理にて制御器Ｈを設計することで、制御器Ｈのパラメータ（制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0、遅延器Ｒｄの遅延時間Ｒd1～Ｒdm2及び遅延器Ｒｃの遅延時間Ｒc1～Ｒcm1）を決定する。 [Processing of control unit 10]
FIG. 7 is a flowchart showing a processing example of the control unit 10. Hereinafter, a processing example of the control unit 10 will be described in detail. As described above, the control unit 10 designs the controller H by processing the controller design program based on the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2. , The parameters of the controller H (coefficient matrix A0, B0, C0, D0 of the controller H0, the delay times Rd1 to Rdm2 of the delay device Rd, and the delay times Rc1 to Rcm1 of the delay device Rc) are determined.

（制御対象Ｇｐ及び目標システムＧｔを同定：ステップＳ７０１）
制御部１０の同定部２０は、音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に基づいて、制御対象Ｇｐ及び目標システムＧｔを同定する（ステップＳ７０１）。 (Identification of controlled target Gp and target system Gt: step S701)
The identification unit 20 of the control unit 10 identifies the control target Gp and the target system Gt based on the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 (step S701).

これにより、制御対象Ｇｐ及び目標システムＧｔを状態空間モデルで表した場合の係数行列が得られる。 As a result, a coefficient matrix when the controlled object Gp and the target system Gt are represented by a state space model can be obtained.

図８は、制御対象Ｇｐ及び目標システムＧｔを結合した結合システムＧｗを説明する図である。制御対象Ｇｐにおける２次音源であるスピーカ１０１－１～１０１－ｍ１の数はｍ１であり、目標システムＧｔにおける仮想音源の数はｍ２であり、制御点である耳１０５－１，１０５－２の数はｐ＝２である。 FIG. 8 is a diagram illustrating a coupling system Gw in which a controlled target Gp and a target system Gt are coupled. The number of speakers 101-1 to 101-m1 which are secondary sound sources in the controlled target Gp is m1, the number of virtual sound sources in the target system Gt is m2, and the ears 105-1 and 105-2 which are control points. The number is p = 2.

制御対象Ｇｐ及び目標システムＧｔを結合したシステムを結合システムＧｗと定義すると、結合システムＧｗは、ｍ１＋ｍ２入力及びｐ出力のシステムとなる。結合システムＧｗにおいて、第１番目～第ｍ１番目の入力を制御対象Ｇｐの入力とし、第ｍ１＋１番目～第ｍ１＋ｍ２番目の入力を目標システムＧｔの入力とする。 If the system in which the controlled target Gp and the target system Gt are combined is defined as the combined system Gw, the combined system Gw becomes a system of m1 + m2 input and p output. In the combined system Gw, the first to m1st inputs are the inputs of the control target Gp, and the m1 + 1st to the m1 + m2th inputs are the inputs of the target system Gt.

具体的には、まず、同定部２０は、直達項を持たない状態空間モデルで表した結合システムＧｗに対し、既存の手法を用いて、結合システムＧｗの係数行列Ａｗ，Ｂｗ，Ｃｗを求める。 Specifically, first, the identification unit 20 obtains the coefficient matrices Aw, Bw, and Cw of the coupling system Gw by using the existing method for the coupling system Gw represented by the state space model having no direct term.

例えば、２次音源の音響信号u1-1～u1-m1、仮想音源の音響信号u2-1～u2-m2、及びその応答である制御点の音響信号y1,y2が測定（収音）されるものとする。同定部２０は、これらの測定信号である音響信号を入力し、２次音源及び仮想音源と制御点との間のインパルス応答を算出する。そして、同定部２０は、インパルス応答に基づいて、特異値分解法により、係数行列Ａｗ，Ｂｗ，Ｃｗを求める。尚、同定部２０は、これらの測定信号である音響信号を入力し、部分空間法により、係数行列Ａｗ，Ｂｗ，Ｃｗを求めるようにしてもよい。 For example, the acoustic signals u1-1 to u1-m1 of the secondary sound source, the acoustic signals u2-1 to u2-m2 of the virtual sound source, and the acoustic signals y1 and y2 of the control points that are the responses thereof are measured (collected). It shall be. The identification unit 20 inputs these measurement signals, acoustic signals, and calculates an impulse response between the secondary sound source and the virtual sound source and the control point. Then, the identification unit 20 obtains the coefficient matrices Aw, Bw, and Cw by the singular value decomposition method based on the impulse response. The identification unit 20 may input the acoustic signals which are the measurement signals and obtain the coefficient matrices Aw, Bw, and Cw by the subspace method.

尚、特異値分解法及び部分空間法は既知であるから、ここでは詳細な説明を省略する。特異値分解法及び部分空間法は一例であり、他の手法を用いるようにしてもよい。音響信号u1-1～u1-m1,u2-1～u2-m2の伝搬にはむだ時間が含まれるため、直達項の係数行列Ｄｗはゼロとする。 Since the singular value decomposition method and the subspace method are known, detailed description thereof will be omitted here. The singular value decomposition method and the subspace method are examples, and other methods may be used. Since the propagation of the acoustic signals u1-1 to u1-m1 and u2-1 to u2-m2 includes dead time, the coefficient matrix Dw of the direct term is set to zero.

直達項を持たない状態空間モデルで表した結合システムＧｗは、入力信号をｕ、出力信号をｙ、状態変数をｘとして、以下の式にて表される。

｛Ａｗ，Ｂｗ，Ｃｗ｝は、それぞれ結合システムＧｗの係数行列（システム行列）、入力側の係数行列（入力行列）、出力側の係数行列（出力行列）である。ここで、係数行列Ｂｗの１～ｍ１列をＢｐとし、ｍ１＋１～ｍ２列をＢｔとする。 The coupling system Gw represented by a state-space model having no direct term is expressed by the following equation, where u is an input signal, y is an output signal, and x is a state variable.

{Aw, Bw, Cw} are a coefficient matrix (system matrix) of the coupling system Gw, a coefficient matrix on the input side (input matrix), and a coefficient matrix on the output side (output matrix), respectively. Here, the 1st to m1 columns of the coefficient matrix Bw are referred to as Bp, and the m1 + 1 to m2 columns are referred to as Bt.

直達項を持たない状態空間モデルで表した制御対象Ｇｐは、入力信号をｕ、出力信号をｙ、状態変数をｘとして、以下の式にて表される。

｛Ａｗ，Ｂｐ，Ｃｗ｝は、それぞれ制御対象Ｇｐの係数行列（システム行列）、入力側の係数行列（入力行列）、出力側の係数行列（出力行列）である。 The controlled object Gp represented by the state space model having no direct term is expressed by the following equation, where u is an input signal, y is an output signal, and x is a state variable.

{Aw, Bp, Cw} are a coefficient matrix (system matrix) of the controlled object Gp, a coefficient matrix on the input side (input matrix), and a coefficient matrix on the output side (output matrix), respectively.

直達項を持たない状態空間モデルで表した目標システムＧｔは、入力信号をｕ、出力信号をｙ、状態変数をｘとして、以下の式にて表される。

｛Ａｗ，Ｂｔ，Ｃｗ｝は、それぞれ目標システムＧｔの係数行列（システム行列）、入力側の係数行列（入力行列）、出力側の係数行列（出力行列）である。 The target system Gt represented by a state-space model having no direct term is expressed by the following equation, where u is an input signal, y is an output signal, and x is a state variable.

{Aw, Bt, Cw} are a coefficient matrix (system matrix) of the target system Gt, a coefficient matrix on the input side (input matrix), and a coefficient matrix on the output side (output matrix), respectively.

このように、同定部２０により、制御対象Ｇｐについて、直達項を持たない状態空間モデルで表した場合の係数行列Ａｗ，Ｂｐ，Ｃｗが得られ、目標システムＧｔについて、直達項を持たない状態空間モデルで表した場合の係数行列Ａｗ，Ｂｔ，Ｃｗが得られる。 In this way, the identification unit 20 obtains the coefficient matrices Aw, Bp, and Cw when the controlled object Gp is represented by a state space model having no direct term, and the target system Gt is a state space having no direct term. The coefficient matrices Aw, Bt, and Cw expressed by the model can be obtained.

（制御対象Ｇｐ及び目標システムＧｔから遅延時間を分離：ステップＳ７０２）
制御部１０の遅延時間分離部２１は、音響信号u1-1～u1-m1,y1,y2に基づいて、制御対象Ｇｐの遅延時間を求め、制御対象Ｇｐから遅延時間を分離し、制御対象（分離制御対象）Ｇｐ＃を求める。また、遅延時間分離部２１は、音響信号u2-1～u2-m2,y1,y2に基づいて、目標システムＧｔの遅延時間を求め、目標システムＧｔから遅延時間を分離し、目標システム（分離目標システム）Ｇｔ＃を求める（ステップＳ７０２）。 (Separate the delay time from the controlled target Gp and the target system Gt: step S702)
The delay time separation unit 21 of the control unit 10 obtains the delay time of the control target Gp based on the acoustic signals u1-1 to u1-m1, y1, y2, separates the delay time from the control target Gp, and controls the control target ( Separation control target) Gp # is obtained. Further, the delay time separation unit 21 obtains the delay time of the target system Gt based on the acoustic signals u2-1 to u2-m2, y1, y2, separates the delay time from the target system Gt, and separates the delay time from the target system Gt. System) Gt # is obtained (step S702).

これにより、制御対象Ｇｐ＃及び目標システムＧｔ＃を状態空間モデルで表した場合の係数行列が得られる。 As a result, a coefficient matrix when the controlled object Gp # and the target system Gt # are represented by a state space model can be obtained.

具体的には、遅延時間分離部２１は、制御対象Ｇｐに対し、第ｉ（ｉ＝１，・・・，ｍ１）番目の音響信号u1-iを入力してから、出力信号である音響信号y1,y2のうちのいずれかから非零の値が出力されるまでの間の時間を求める。そして、遅延時間分離部２１は、当該時間から１サンプル時間減算し、減算結果の時間を第ｉ番目の音響信号u1-iの遅延時間とする。 Specifically, the delay time separation unit 21 inputs the i (i = 1, ..., M1) th acoustic signal u1-i to the controlled target Gp, and then outputs the acoustic signal, which is an output signal. Find the time between y1 and y2 until a non-zero value is output. Then, the delay time separation unit 21 subtracts one sample time from the time, and sets the time of the subtraction result as the delay time of the i-th acoustic signal u1-i.

図９は、遅延時間を分離した制御対象Ｇｐ＃を説明する図である。制御対象Ｇｐは、入力した信号を遅延時間分遅延させる遅延器Ｒｐ、及び遅延時間を分離した制御対象Ｇｐ＃を備えて構成される。 FIG. 9 is a diagram illustrating a controlled target Gp # in which the delay time is separated. The control target Gp is configured to include a delay device Rp that delays the input signal by the delay time, and a control target Gp # in which the delay time is separated.

遅延器Ｒｐは、音響信号u1-1～u1-m1のそれぞれに対応した遅延時間（遅延要素）を対角成分に持つ伝達関数行列にて表される。つまり、遅延器Ｒｐのｉ行ｉ列の要素は、制御対象Ｇｐのｉ番目の入力に対応する遅延器である。 The delay device Rp is represented by a transfer function matrix having a delay time (delay element) corresponding to each of the acoustic signals u1-1 to u1-m1 as diagonal components. That is, the element of the i-th row and i-column of the delay device Rp is the delay device corresponding to the i-th input of the controlled target Gp.

遅延時間を分離した制御対象Ｇｐ＃は、直達項を持たない状態空間モデルで表され、入力信号をｕ、出力信号をｙ、状態変数をｘとして、以下の式にて表される。

｛Ａｗ＃，Ｂｐ＃，Ｃｗ＃｝は、それぞれ制御対象Ｇｐ＃の係数行列（システム行列）、入力側の係数行列（入力行列）、出力側の係数行列（出力行列）である。 The controlled object Gp # with the delay time separated is represented by a state space model having no direct term, and is represented by the following equation, where u is an input signal, y is an output signal, and x is a state variable.

{Aw #, Bp #, Cw #} are a coefficient matrix (system matrix) of the controlled object Gp #, a coefficient matrix on the input side (input matrix), and a coefficient matrix on the output side (output matrix), respectively.

また、遅延時間分離部２１は、目標システムＧｔに対し、第ｉ（ｉ＝１，・・・，ｍ２）番目の音響信号u2-iを入力してから、出力信号である音響信号y1,y2のうちのいずれかから非零の値が出力されるまでの間の時間を求める。そして、遅延時間分離部２１は、当該時間から１サンプル時間減算し、減算結果の時間を第ｉ番目の音響信号u2-iの遅延時間とする。 Further, the delay time separation unit 21 inputs the i-th (i = 1, ..., M2) th acoustic signal u2-i to the target system Gt, and then outputs the acoustic signals y1, y2. Find the time between any of these and the output of a non-zero value. Then, the delay time separation unit 21 subtracts one sample time from the time, and sets the time of the subtraction result as the delay time of the i-th acoustic signal u2-i.

図１０は、遅延時間を分離した目標システムＧｔ＃を説明する図である。目標システムＧｔは、入力した信号を遅延時間分遅延させる遅延器Ｒｔ、及び遅延時間を分離した目標システムＧｔ＃を備えて構成される。 FIG. 10 is a diagram illustrating a target system Gt # in which delay times are separated. The target system Gt includes a delay device Rt that delays the input signal by the delay time, and a target system Gt # that separates the delay times.

遅延器Ｒｔは、音響信号u2-1～u2-m2のそれぞれに対応した遅延時間（遅延要素）を対角成分に持つ伝達関数行列にて表される。つまり、遅延器Ｒｔのｉ行ｉ列の要素は、目標システムＧｔのｉ番目の入力に対応する遅延器である。 The delay device Rt is represented by a transfer function matrix having a delay time (delay element) corresponding to each of the acoustic signals u2-1 to u2-m2 as diagonal components. That is, the element of the i-th row and i-th column of the delay device Rt is the delay device corresponding to the i-th input of the target system Gt.

遅延時間を分離した目標システムＧｔ＃は、直達項を持たない状態空間モデルで表され、入力信号をｕ、出力信号をｙ、状態変数をｘとして、以下の式にて表される。

｛Ａｗ＃，Ｂｔ＃，Ｃｗ＃｝は、それぞれ目標システムＧｔ＃の係数行列（システム行列）、入力側の係数行列（入力行列）、出力側の係数行列（出力行列）である。 The target system Gt # in which the delay time is separated is represented by a state space model having no direct term, and is represented by the following equation, where u is an input signal, y is an output signal, and x is a state variable.

{Aw #, Bt #, Cw #} are the coefficient matrix (system matrix) of the target system Gt #, the coefficient matrix on the input side (input matrix), and the coefficient matrix on the output side (output matrix), respectively.

このように、遅延時間分離部２１により、遅延時間を分離した制御対象Ｇｐ＃について、直達項を持たない状態空間モデルで表した場合の係数行列Ａｗ＃，Ｂｐ＃，Ｃｗ＃が得られる。また、遅延時間を分離した目標システムＧｔ＃について、直達項を持たない状態空間モデルで表した場合の係数行列Ａｗ＃，Ｂｔ＃，Ｃｗ＃が得られる。 In this way, the delay time separation unit 21 obtains coefficient matrices Aw #, Bp #, and Cw # when the controlled object Gp # with the delay time separated is represented by a state-space model having no direct term. Further, for the target system Gt # in which the delay time is separated, the coefficient matrices Aw #, Bt #, and Cw # when expressed by a state space model having no direct term can be obtained.

（全体システムＨallを構成：ステップＳ７０３）
制御部１０のシステム構成部２２は、制御器Ｈ0、ステップＳ７０２にて遅延時間が分離された制御対象Ｇｐ＃及び目標システムＧｔ＃、減算器並びにバランス係数αを用いて、全体システムＨallを構成する（ステップＳ７０３）。 (Constructing the entire system Hall: step S703)
The system configuration unit 22 of the control unit 10 configures the entire system Hall by using the controller H0, the control target Gp # whose delay time is separated in step S702, the target system Gt #, the subtractor, and the balance coefficient α. (Step S703).

図１１は、全体システムＨallの構成例を示すブロック線図である。この全体システムＨallは、制御器Ｈ0、制御対象Ｇｐ＃、目標システムＧｔ＃、演算器αＩ及び減算器４７を備えている。全体システムＨallの入力信号ｕは音響信号u2-1～u2-m2であり、出力信号y1は制御器Ｈ0のゲインにより増幅される信号であり、制御器Ｈ0のゲインに比例する。出力信号y2は、制御対象Ｇｐ＃と目標システムＧｔ＃との間の出力追従誤差である。バランス係数αは、２つの出力である制御器Ｈ0のゲインにより増幅される信号と制御対象Ｇｐ＃及び目標システムＧｔ＃の出力追従誤差との間のバランスを決定するためのパラメータである。 FIG. 11 is a block diagram showing a configuration example of the entire system Hall. This whole system Hall includes a controller H0, a controlled object Gp #, a target system Gt #, an arithmetic unit αI, and a subtractor 47. The input signals u of the entire system Hall are acoustic signals u2-1 to u2-m2, and the output signal y1 is a signal amplified by the gain of the controller H0 and is proportional to the gain of the controller H0. The output signal y2 is an output tracking error between the controlled target Gp # and the target system Gt #. The balance coefficient α is a parameter for determining the balance between the signal amplified by the gain of the controller H0, which is two outputs, and the output tracking error of the controlled object Gp # and the target system Gt #.

制御器Ｈ0及び目標システムＧｔ＃は、入力信号ｕを入力する。演算器αＩは、制御器Ｈ0により出力された信号にαＩを乗算し、乗算結果である制御器Ｈ0のゲインにより増幅される信号を出力信号y1として出力する。Ｉは単位行列である。制御対象Ｇｐ＃は、制御器Ｈ0により出力された信号を入力する。減算器４７は、制御対象Ｇｐ＃により出力された信号から目標システムＧｔ＃により出力された信号を減算し、減算結果である制御対象Ｇｐ＃と目標システムＧｔ＃との間の出力追従誤差を出力信号y2として出力する。 The controller H0 and the target system Gt # input the input signal u. The arithmetic unit αI multiplies the signal output by the controller H0 by αI, and outputs the signal amplified by the gain of the controller H0 which is the multiplication result as the output signal y1. I is an identity matrix. The control target Gp # inputs the signal output by the controller H0. The subtractor 47 subtracts the signal output by the target system Gt # from the signal output by the control target Gp #, and outputs an output tracking error between the control target Gp # and the target system Gt #, which is the subtraction result. Output as signal y2.

制御器Ｈ0は、直達項を持つ状態空間モデルで表され、入力信号をｕ、出力信号をｙ、状態変数をｘとして、以下の式にて表される。

｛Ａ0，Ｂ0，Ｃ0，Ｄ0｝は、それぞれ制御器Ｈ0の係数行列（システム行列）、入力側の係数行列（入力行列）、出力側の係数行列（出力行列）、直達項の係数行列である。 The controller H0 is represented by a state space model having a direct term, and is represented by the following equation, where u is an input signal, y is an output signal, and x is a state variable.

{A0, B0, C0, D0} are the coefficient matrix (system matrix) of the controller H0, the coefficient matrix on the input side (input matrix), the coefficient matrix on the output side (output matrix), and the coefficient matrix of the direct term, respectively. ..

図１１に示した全体システムＨallは、出力信号y1,y2（制御器Ｈ0のゲインにより増幅される信号、及び制御対象Ｇｐ＃と目標システムＧｔ＃との間の出力追従誤差）を共に小さくすることが目的となる。 The overall system Hall shown in FIG. 11 reduces both the output signals y1 and y2 (the signal amplified by the gain of the controller H0 and the output tracking error between the controlled target Gp # and the target system Gt #). Is the purpose.

このように、システム構成部２２により、入力信号ｕを音響信号u2-1～u2-m2とし、出力信号y1,y2を、制御器Ｈ0のゲインにより増幅される信号、及び制御対象Ｇｐ＃と目標システムＧｔ＃との間の出力追従誤差とし、係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0で表される制御器Ｈ0、係数行列Ａｗ＃，Ｂｐ＃，Ｃｗ＃で表される制御対象Ｇｐ＃、係数行列Ａｗ＃，Ｂｔ＃，Ｃｗ＃で表される目標システムＧｔ＃、減算器４７及びバランス係数αを用いた全体システムＨallが構成される。 In this way, the system configuration unit 22 sets the input signal u as the acoustic signal u2-1 to u2-m2, and the output signals y1 and y2 as the signal amplified by the gain of the controller H0, the control target Gp #, and the target. The output tracking error with the system Gt # is the controller H0 represented by the coefficient matrix A0, B0, C0, D0, the controlled object Gp # represented by the coefficient matrix Aw #, Bp #, Cw #, and the coefficient matrix. The whole system Hall using the target system Gt # represented by Aw #, Bt #, and Cw #, the subtractor 47, and the balance coefficient α is configured.

（バランス係数αの初期値を設定：ステップＳ７０４）
制御部１０のパラメータ決定部２３は、制御器Ｈ0のゲインにより増幅される信号と制御対象Ｇｐ＃及び目標システムＧｔ＃の出力追従誤差との間のバランスを決定するためのバランス係数αの初期値として、予め設定された十分に小さい任意の値を設定する（ステップＳ７０４）。 (Set the initial value of the balance coefficient α: step S704)
The parameter determination unit 23 of the control unit 10 determines the initial value of the balance coefficient α for determining the balance between the signal amplified by the gain of the controller H0 and the output tracking error of the control target Gp # and the target system Gt #. , A preset and sufficiently small arbitrary value is set (step S704).

バランス係数αが大きいほど、制御器Ｈ0のゲインにより増幅される信号に比重が置かれ、バランス係数αが小さいほど、制御対象Ｇｐ＃と目標システムＧｔ＃との間の出力追従誤差に比重が置かれる。 The larger the balance coefficient α, the more weight is placed on the signal amplified by the gain of the controller H0, and the smaller the balance coefficient α, the more weight is placed on the output tracking error between the controlled target Gp # and the target system Gt #. It will be taken.

（パラメータ（係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0）を決定：ステップＳ７０５～Ｓ７０９）
パラメータ決定部２３は、バランス係数αを用い、ステップＳ７０３にて構成した全体システムＨallのＨ∞ノルムを評価関数として、制御器Ｈ0のゲイン、及び制御対象Ｇｐ＃と目標システムＧｔ＃との間の出力追従誤差を定量化し、全体システムＨallのＨ∞ノルムを最小化する制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定する（ステップＳ７０５）。 (Determining the parameters (coefficient matrix A0, B0, C0, D0): steps S705 to S709)
The parameter determination unit 23 uses the balance coefficient α and uses the H∞ norm of the entire system Hall configured in step S703 as an evaluation function to obtain the gain of the controller H0 and between the control target Gp # and the target system Gt #. The coefficient matrix A0, B0, C0, D0 of the controller H0 that quantifies the output tracking error and minimizes the H∞ norm of the entire system Hall is determined (step S705).

ここで、Ｈ∞ノルムとゲインとは、Ｈ∞ノルムが大きいほどゲインの上限が大きくなり、Ｈ∞ノルムが小さいほどゲインの上限が小さくなる関係にある。制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定する処理の詳細については後述する。 Here, the H∞ norm and the gain have a relationship that the larger the H∞ norm, the larger the upper limit of the gain, and the smaller the H∞ norm, the smaller the upper limit of the gain. The details of the process for determining the coefficient matrices A0, B0, C0, and D0 of the controller H0 will be described later.

ステップＳ７０５及び後述するステップＳ７０６～Ｓ７０９により、出力追従誤差をある程度許容しながら（大きくしながら）、ゲインが小さくなるように、制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0が決定される。 In steps S705 and steps S706 to S709 described later, the coefficient matrices A0, B0, C0, and D0 of the controller H0 are determined so that the gain becomes small while allowing the output tracking error to some extent (while increasing).

パラメータ決定部２３は、ステップＳ７０５から移行して、制御器Ｈ0のＨ∞ノルムであるゲインを算出する。例えば、パラメータ決定部２３は、制御器Ｈ0の全周波数について特異値を算出し、周波数毎の特異値の中から最大値を特定し、当該最大値を制御器Ｈ0のＨ∞ノルムとする。制御器Ｈ0のＨ∞ノルムを算出する処理は既知であるから、ここでは詳細な説明を省略する。そして、パラメータ決定部２３は、制御器Ｈ0のＨ∞ノルムが予め設定された閾値以下であるか否かを判定する（ステップＳ７０６）。 The parameter determination unit 23 shifts from step S705 to calculate the gain which is the H∞ norm of the controller H0. For example, the parameter determination unit 23 calculates a singular value for all frequencies of the controller H0, specifies the maximum value from the singular values for each frequency, and sets the maximum value as the H∞ norm of the controller H0. Since the process of calculating the H∞ norm of the controller H0 is known, detailed description thereof will be omitted here. Then, the parameter determination unit 23 determines whether or not the H∞ norm of the controller H0 is equal to or less than a preset threshold value (step S706).

パラメータ決定部２３は、ステップＳ７０６において、制御器Ｈ0のＨ∞ノルムが予め設定された閾値以下でないと判定した場合（ステップＳ７０６：Ｎ）、バランス係数αに予め設定された定数（正の実数）を加算することで、バランス係数αが少し大きくなるように設定する（ステップＳ７０７）。そして、パラメータ決定部２３は、ステップＳ７０５へ移行し、新たなバランス係数αのときの係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定する。 When the parameter determination unit 23 determines in step S706 that the H∞ norm of the controller H0 is not equal to or less than the preset threshold value (step S706: N), the constant (positive real number) preset in the balance coefficient α. Is added so that the balance coefficient α becomes slightly larger (step S707). Then, the parameter determination unit 23 proceeds to step S705 and determines the coefficient matrices A0, B0, C0, and D0 when the new balance coefficient α is used.

パラメータ決定部２３は、ステップＳ７０６において、制御器Ｈ0のＨ∞ノルムが予め設定された閾値以下であると判定した場合（ステップＳ７０６：Ｙ）、ステップＳ７０８へ移行する。 When the parameter determination unit 23 determines in step S706 that the H∞ norm of the controller H0 is equal to or less than a preset threshold value (step S706: Y), the parameter determination unit 23 proceeds to step S708.

パラメータ決定部２３は、ステップＳ７０６から移行して、全体システムＨallのＨ∞ノルムであるゲインを算出する。例えば、パラメータ決定部２３は、全体システムＨallの全周波数について特異値を算出し、周波数毎の特異値の中から最大値を特定し、当該最大値を全体システムＨallのＨ∞ノルムとする。全体システムＨallのＨ∞ノルムを算出する処理は既知であるから、ここでは詳細な説明を省略する。そして、パラメータ決定部２３は、全体システムＨallのＨ∞ノルムが予め設定された閾値以下であるか否かを判定する（ステップＳ７０８）。 The parameter determination unit 23 shifts from step S706 to calculate the gain which is the H∞ norm of the entire system Hall. For example, the parameter determination unit 23 calculates a singular value for all frequencies of the entire system Hall, specifies the maximum value from the singular values for each frequency, and sets the maximum value as the H∞ norm of the entire system Hall. Since the process of calculating the H∞ norm of the entire system Hall is known, detailed description thereof will be omitted here. Then, the parameter determination unit 23 determines whether or not the H∞ norm of the entire system Hall is equal to or less than a preset threshold value (step S708).

パラメータ決定部２３は、ステップＳ７０８において、全体システムＨallのＨ∞ノルムが予め設定された閾値以下でないと判定した場合（ステップＳ７０８：Ｎ）、バランス係数αから予め設定された定数（正の実数）を減算することで、バランス係数αが少し小さくなるように設定する（ステップＳ７０９）。そして、パラメータ決定部２３は、ステップＳ７０５へ移行し、新たなバランス係数αのときの係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定する。 When the parameter determination unit 23 determines in step S708 that the H∞ norm of the entire system Hall is not equal to or less than a preset threshold value (step S708: N), a constant (positive real number) preset from the balance coefficient α. Is set so that the balance coefficient α becomes slightly smaller by subtracting (step S709). Then, the parameter determination unit 23 proceeds to step S705 and determines the coefficient matrices A0, B0, C0, and D0 when the new balance coefficient α is used.

パラメータ決定部２３は、ステップＳ７０８において、全体システムＨallのＨ∞ノルムが予め設定された閾値以下であると判定した場合（ステップＳ７０８：Ｙ）、ステップＳ７１０へ移行する。 When the parameter determination unit 23 determines in step S708 that the H∞ norm of the entire system Hall is equal to or less than a preset threshold value (step S708: Y), the parameter determination unit 23 proceeds to step S710.

ここで、ステップＳ７０７にてバランス係数αが大きくなるように設定されることにより、ステップＳ７０５にて、新たなバランス係数αを用いて制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0が決定される。そして、ステップＳ７０６にて算出される制御器Ｈ0のＨ∞ノルムであるゲインは、以前よりも小さくなる。 Here, by setting the balance coefficient α to be large in step S707, the coefficient matrices A0, B0, C0, and D0 of the controller H0 are determined in step S705 using the new balance coefficient α. To. Then, the gain, which is the H∞ norm of the controller H0 calculated in step S706, becomes smaller than before.

これは、図１１を参照して、バランス係数αが大きくなると、出力信号y1である制御器Ｈ0のゲインにより増幅された信号は、出力追従誤差y2と比較して全体システムＨallのＨ∞ノルムの大きさへの寄与が大きくなる。新たな制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0は、全体システムＨallのＨ∞ノルムを最小化する方向で決定され、結果として全体システムＨallのＨ∞ノルムへの寄与の大きい新たな制御器Ｈ0のゲインが抑制されるからである。つまり、バランス係数αが大きくなるほど、制御器Ｈ0のゲインは小さくなる。その一方で、出力追従誤差y2は相対的に大きくなり、全体システムＨallのＨ∞ノルムも大きくなる。 This is because, referring to FIG. 11, when the balance coefficient α becomes large, the signal amplified by the gain of the controller H0 which is the output signal y1 becomes the H∞ norm of the whole system Hall as compared with the output tracking error y2. The contribution to size increases. The coefficient matrix A0, B0, C0, D0 of the new controller H0 is determined in the direction of minimizing the H∞ norm of the whole system Hall, and as a result, the new control that greatly contributes to the H∞ norm of the whole system Hall. This is because the gain of the device H0 is suppressed. That is, the larger the balance coefficient α, the smaller the gain of the controller H0. On the other hand, the output tracking error y2 becomes relatively large, and the H∞ norm of the entire system Hall also becomes large.

一方、ステップＳ７０９にてバランス係数αが小さくなるように設定されることにより、ステップＳ７０５にて、新たなバランス係数αを用いて制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0が決定される。そして、ステップＳ７０６にて算出される制御器Ｈ0のＨ∞ノルムであるゲインは、以前よりも大きくなる。つまり、バランス係数αが小さくなるほど、制御器Ｈ0のゲインは大きくなる。その一方で、出力追従誤差y2は相対的に小さくなり、全体システムＨallのＨ∞ノルムも小さくなる。 On the other hand, by setting the balance coefficient α to be small in step S709, the coefficient matrices A0, B0, C0, D0 of the controller H0 are determined in step S705 using the new balance coefficient α. .. Then, the gain, which is the H∞ norm of the controller H0 calculated in step S706, becomes larger than before. That is, the smaller the balance coefficient α, the larger the gain of the controller H0. On the other hand, the output tracking error y2 becomes relatively small, and the H∞ norm of the entire system Hall also becomes small.

このように、パラメータ決定部２３により、制御対象Ｇｐ＃と目標システムＧｔ＃との出力追従誤差をある程度許容しながら（大きくしながら）、制御器Ｈ0ゲインが小さくなるように、制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0が決定される。これにより、図５に示した制御器Ｈのパラメータのうち、制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0が決定される。 In this way, the parameter determination unit 23 allows the output tracking error between the controlled target Gp # and the target system Gt # to some extent (while increasing), and the coefficient of the controller H0 is reduced so that the controller H0 gain becomes smaller. The matrices A0, B0, C0, D0 are determined. As a result, among the parameters of the controller H shown in FIG. 5, the coefficient matrices A0, B0, C0, and D0 of the controller H0 are determined.

（遅延時間を補償：ステップＳ７１０）
遅延時間補償部２４は、ステップＳ７０２にて制御対象Ｇｐ及び目標システムＧｔから分離した遅延時間を補償するため、当該遅延時間に基づいて、図５に示した制御器Ｈに含まれる遅延器Ｒｄの遅延時間Ｒd1～Ｒdm2及び遅延器Ｒｃの遅延時間Ｒc1～Ｒcm1を算出する（ステップＳ７１０）。 (Compensation for delay time: step S710)
In order to compensate for the delay time separated from the controlled target Gp and the target system Gt in step S702, the delay time compensating unit 24 of the delay device Rd included in the controller H shown in FIG. 5 based on the delay time. The delay times Rd1 to Rdm2 and the delay times Rc1 to Rcm1 of the delay device Rc are calculated (step S710).

これにより、ステップＳ７０２にて制御対象Ｇｐ及び目標システムＧｔから分離した遅延時間を補償するための遅延器Ｒｄ，遅延器Ｒｃが構成される。 As a result, the delay device Rd and the delay device Rc for compensating for the delay time separated from the controlled target Gp and the target system Gt in step S702 are configured.

具体的には、遅延時間補償部２４は、ステップＳ７０２にて目標システムＧｔから分離したそれぞれの遅延時間（音響信号u2-1～u2-m2に対応するそれぞれの遅延時間）のうち、最小時間をβとする。遅延時間補償部２４は、最小時間βから、ステップＳ７０２にて制御対象Ｇｐから分離したそれぞれの遅延時間を減算し、遅延器Ｒｃを構成する遅延器３１－１～３１－ｍ１のそれぞれの遅延時間を求める。 Specifically, the delay time compensation unit 24 sets the minimum time among the delay times (each delay time corresponding to the acoustic signals u2-1 to u2-m2) separated from the target system Gt in step S702. Let it be β. The delay time compensation unit 24 subtracts each delay time separated from the controlled target Gp in step S702 from the minimum time β, and delay times of the delay devices 31-1 to 31-m1 constituting the delay device Rc. Ask for.

遅延時間補償部２４は、ステップＳ７０２にて目標システムＧｔから分離したそれぞれの遅延時間（音響信号u2-1～u2-m2に対応するそれぞれの遅延時間）から、最小時間βを減算し、遅延器Ｒｄを構成する遅延器３０－１～３０－ｍ２のそれぞれの遅延時間を求める。 The delay time compensation unit 24 subtracts the minimum time β from each delay time (each delay time corresponding to the acoustic signals u2-1 to u2-m2) separated from the target system Gt in step S702, and the delay device 24. The delay times of the delay devices 30-1 to 30-m2 constituting Rd are obtained.

ただし、遅延時間補償部２４は、遅延器Ｒｃを構成する遅延器３１－１～３１－ｍ１の遅延時間のうち、いずれかの遅延時間が負値ｒとなった場合、当該遅延時間が０となるように、それぞれの遅延時間に絶対値｜ｒ｜を加算する。これにより、遅延器Ｒｃを構成する遅延器３１－１～３１－ｍ１の全ての遅延時間が０以上の値となる。 However, the delay time compensation unit 24 sets the delay time to 0 when any of the delay times of the delay devices 31-1 to 31-m1 constituting the delay device Rc has a negative value r. The absolute value | r | is added to each delay time so as to be. As a result, all the delay times of the delay devices 31-1 to 31-m1 constituting the delay device Rc become a value of 0 or more.

このように、遅延時間補償部２４により、ステップＳ７０２にて制御対象Ｇｐ及び目標システムＧｔから分離した遅延時間を補償するための遅延器Ｒｄ及び遅延器Ｒｃが構成される。これにより、図５に示した制御器Ｈのパラメータのうち、遅延器Ｒｄの遅延時間Ｒd1～Ｒdm2及び遅延器Ｒｃの遅延時間Ｒc1～Ｒcm1が決定される。 As described above, the delay time compensating unit 24 configures the delay device Rd and the delay device Rc for compensating the delay time separated from the controlled target Gp and the target system Gt in step S702. As a result, among the parameters of the controller H shown in FIG. 5, the delay times Rd1 to Rdm2 of the delay device Rd and the delay times Rc1 to Rcm1 of the delay device Rc are determined.

以上のように、制御器設計装置１の制御部１０により、音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に基づいて、制御器設計プログラムの処理にて制御器Ｈが設計され、制御器Ｈのパラメータ（制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0、遅延器Ｒｄの遅延時間Ｒd1～Ｒdm2及び遅延器Ｒｃの遅延時間Ｒc1～Ｒcm1）が決定される。 As described above, the control unit 10 of the controller design device 1 controls by processing the controller design program based on the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2. The device H is designed, and the parameters of the controller H (coefficient matrix A0, B0, C0, D0 of the controller H0, the delay time Rd1 to Rdm2 of the delay device Rd, and the delay time Rc1 to Rcm1 of the delay device Rc) are determined. ..

〔制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0の決定処理（ステップＳ７０５）〕
次に、図７のステップＳ７０５に示した、パラメータ決定部２３による制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0の決定処理について詳細に説明する。 [Determining process of coefficient matrix A0, B0, C0, D0 of controller H0 (step S705)]
Next, the determination process of the coefficient matrices A0, B0, C0, D0 of the controller H0 by the parameter determination unit 23 shown in step S705 of FIG. 7 will be described in detail.

この手法は、全体システムＨallのＨ∞ノルムを最小化する条件式を、線形行列不等式（ＬＭＩ）の制約のもとでの最適化問題で表し、当該最適化問題を求解することで、全体システムＨallのＨ∞ノルムを最小化する制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定するものである。 In this method, the conditional expression that minimizes the H∞ norm of the whole system Hall is expressed by an optimization problem under the constraint of the linear matrix inequality (LMI), and the whole system is solved by solving the optimization problem. It determines the coefficient matrices A0, B0, C0, D0 of the controller H0 that minimizes the H∞ norm of Hall.

図１１に示した全体システムＨallは、直達項を持つ状態空間モデルで表され、入力信号をｕ、出力信号をｙ、状態変数をｘとして、以下の式にて表される。

さらに、y1及びy2を以下の式で表す。

The whole system Hall shown in FIG. 11 is represented by a state space model having a direct term, and is represented by the following equation, where u is an input signal, y is an output signal, and x is a state variable.

Furthermore, y1 and y2 are expressed by the following equations.

そうすると、図１１は図１２のように簡単化される。図１２は、図１１を簡単化したブロック線図である。図１２に示したブロック線図から、全体システムＨallのＨ∞ノルムを最小化する制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定する。 Then, FIG. 11 is simplified as shown in FIG. FIG. 12 is a block diagram obtained by simplifying FIG. 11. From the block diagram shown in FIG. 12, the coefficient matrices A0, B0, C0, D0 of the controller H0 that minimizes the H∞ norm of the entire system Hall are determined.

この制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0の決定には、線形行列不等式の制約のもとでの最適化問題を用いる。全体システムＨallについて、Ｈ∞性能を特徴づける線形行列不等式を記述することができれば、既存の計算アルゴリズムを適用することにより、係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0の数値解を得ることができる。 An optimization problem under the constraint of the linear matrix inequality is used to determine the coefficient matrices A0, B0, C0, D0 of the controller H0. If the linear matrix inequality that characterizes the H∞ performance can be described for the entire system Hall, the numerical solutions of the coefficient matrices A0, B0, C0, and D0 can be obtained by applying the existing calculation algorithm.

そのために、全体システムＨallを表す前記式（１２）（すなわち後述する式（１４））を、後述する式（１６）に代入し、当該式（１６）が、係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0及び後述する正定行列Ｐについて線形行列不等式となるように、当該式（１６）を式変換する。 Therefore, the above equation (12) representing the whole system Hall (that is, the equation (14) described later) is substituted into the equation (16) described later, and the equation (16) is the coefficient matrix A0, B0, C0, D0. And the equation (16) is transformed so that the linear matrix inequality is obtained for the regular matrix P described later.

これにより、当該式（１６）は、全体システムＨallのＨ∞性能を特徴づける線形行列不等式として記述することができるから、全体システムＨallのＨ∞ノルムを最小化する制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定することができる。 As a result, the equation (16) can be described as a linear matrix inequality that characterizes the H∞ performance of the overall system Hall. Therefore, the coefficient matrix A0 of the controller H0 that minimizes the H∞ norm of the overall system Hall, B0, C0, D0 can be determined.

以下、具体的に説明する。前記式（１２）を参照して、全体システムＨallの伝達関数Ｈall（ｚ）を以下の式にて表記することにする（ドイルの記法）。

｛Ａall，Ｂall，Ｃall，Ｄall｝は、全体システムＨallの係数行列（システム行列）、入力側の係数行列（入力行列）、出力側の係数行列（出力行列）、直達項の係数行列である。 Hereinafter, a specific description will be given. With reference to the above equation (12), the transfer function Hall (z) of the whole system Hall will be expressed by the following equation (Doyle notation).

{Aall, Ball, Hall, Dall} is a coefficient matrix (system matrix) of the whole system Hall, a coefficient matrix on the input side (input matrix), a coefficient matrix on the output side (output matrix), and a coefficient matrix of direct terms.

ここで、全体システムＨallのＨ∞ノルムを評価関数として、その最小化問題は、以下のように定式化することができる。

前記式（１５）は、全体システムＨallのＨ∞ノルムがパラメータγよりも小さいという条件の下で、パラメータγを最小化する制御器Ｈ0を求める、すなわち、全体システムＨallのＨ∞ノルムを最小化する制御器Ｈ0を求めることを示している。 Here, the minimization problem can be formulated as follows, using the H∞ norm of the entire system Hall as an evaluation function.

The equation (15) finds the controller H0 that minimizes the parameter γ under the condition that the H∞ norm of the whole system Hall is smaller than the parameter γ, that is, the H∞ norm of the whole system Hall is minimized. It shows that the controller H0 to be used is obtained.

しかしながら、前記式（１４）のままでは、前記式（１５）に基づいて全体システムＨallのＨ∞ノルムを最小化する制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定することができない。そこで、以下に示すように、等価変換と変数変換により前記式（１４）を線形行列式に変換し、全体システムＨallのＨ∞ノルムを最小化する条件式である前記式（１５）を半正定値計画問題に置き換えることで、制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定する。 However, with the equation (14) as it is, the coefficient matrices A0, B0, C0, D0 of the controller H0 that minimizes the H∞ norm of the entire system Hall cannot be determined based on the equation (15). Therefore, as shown below, the equation (14), which is a conditional equation that minimizes the H∞ norm of the entire system Hall, is semidefinitely defined by converting the equation (14) into a linear determinant by equivalent transformation and change of variables. By substituting for the value planning problem, the coefficient matrices A0, B0, C0, D0 of the controller H0 are determined.

既存の定理より、前記式（１５）の小なり条件に係る制約は、以下の式を満たす正定行列Ｐが存在することと等価である。

From the existing theorem, the constraint relating to the lesser condition of the above equation (15) is equivalent to the existence of a positive-definite matrix P satisfying the following equation.

しかしながら、前記式（１４）を代入した前記式（１６）は、｛Ｐ，Ａ0，Ｃ0｝について双線形項を含む双線形行列不等式であるため、数値的な求解が困難である。そこで、前記式（１６）を、線形行列不等式に変換する。 However, since the equation (16) to which the equation (14) is substituted is a bilinear matrix inequality including a bilinear term for {P, A0, C0}, it is difficult to numerically solve it. Therefore, the above equation (16) is converted into a linear matrix inequality.

パラメータ決定部２３は、前記式（１６）が変換された線形行列不等式を用いて、正定行列Ｐが存在し、かつ｛Ａall，Ｂall，Ｃall，Ｄall｝を算出することを前提に、パラメータγを最小化することができる。すなわち、パラメータ決定部２３は、前記式（１６）が変換された線形行列不等式を用いて、パラメータγが最小となるように、係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を算出することができる。 The parameter determination unit 23 sets the parameter γ on the premise that the positive-definite matrix P exists and {Aall, Ball, Hall, Dall} is calculated using the linear matrix inequality obtained by converting the equation (16). Can be minimized. That is, the parameter determination unit 23 can calculate the coefficient matrices A0, B0, C0, and D0 so that the parameter γ is minimized by using the linear matrix inequality obtained by converting the equation (16).

以下、前記式（１６）が、係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0及び正定行列Ｐについて線形行列不等式となるように、式変換する。 Hereinafter, the equation (16) is transformed so that the coefficient matrices A0, B0, C0, D0 and the canonical matrix P become linear matrix inequalities.

まず、正定行列Ｐを、正定行列Ｙ，Ｖ及び正則行列Ｎを用いて、以下の式のように分割する。

First, the positive-definite matrix P is divided by using the positive-definite matrices Y and V and the regular matrix N as shown in the following equation.

次に、前記式（１７）の逆行列Ｐ^－１を、正定行列Ｘ，Ｗ及び正則行列Ｍを用いて、以下の式のように表す。

Next, the inverse matrix P ^-1 of the above equation (17) is expressed as the following equation using the positive-definite matrices X and W and the regular matrix M.

さらに、以下の式のとおり、行列Π₁，Π₂，Πを定義する。

Furthermore, the matrices Π ₁ , Π ₂ , and Π are defined as shown in the following equation.

前記式（１９）を使用して表す以下の式は、前記式（１６）と等価である。

The following equation expressed using the equation (19) is equivalent to the equation (16).

次に、以下の式のとおり、新たな係数行列Ａ0＄，Ｂ0＄，Ｃ0＄，Ｄ0＄を定義する。

Next, a new coefficient matrix A0 $, B0 $, C0 $, D0 $ is defined as shown in the following equation.

これにより、前記式（２０）の要素であるΠ₁ ^TＰΠ₁，Π₁ ^TＰＡallΠ₁，ＣallΠ₁，Π₁ ^TＰＢallは、以下の式にて表される。

As a result, the elements of the above equation (20), Π ₁ ^T P Π ₁ , Π ₁ ^T P Aall Π ₁ , Call Π ₁ , Π ₁ ^T P Ball, are expressed by the following equations.

これにより、前記式（２２）のΠ₁ ^TＰΠ₁，Π₁ ^TＰＡallΠ₁，ＣallΠ₁，Π₁ ^TＰＢallを要素とする前記式（２０）は、以下の式で表される設計パラメータに関して線形となる。

尚、前記式（２０）の要素であるΠ₁ ^TＡall^TＰΠ₁，Π₁ ^TＣall^T，Ｂall^TＰΠ₁は、それぞれΠ₁ ^TＰＡallΠ₁，ＣallΠ₁，Π₁ ^TＰＢallの転置であり、これも同様の変換によって、前記式（２３）で表される設計パラメータに関して線形となる。 As a result, the equation (20) having Π ₁ ^T P Π ₁ , Π ₁ ^T P Aall Π ₁ , Call Π ₁ , Π ₁ ^T P Ball of the equation (22) as elements is linear with respect to the design parameters expressed by the following equations. It becomes.

The elements of the above equation (20), Π ₁ ^T Aall ^T P Π ₁ , Π _{1 T Hall T, and Ball T P Π 1, are transposes of Π 1} ^{T PAall Π} _{1 ,} ^Ball _{Π 1} , and Π ₁ ^T P Ball, respectively. This is also made linear with respect to the design parameter represented by the above equation (23) by the same transformation.

つまり、パラメータ決定部２３は、前記式（２０）の要素であるΠ₁ ^TＰΠ₁，Π₁ ^TＰＡallΠ₁，ＣallΠ₁，Π₁ ^TＰＢall及びその転置であるΠ₁ ^TＡall^TＰΠ₁，Π₁ ^TＣall^T，Ｂall^TＰΠ₁を前記式（２１）及び前記式（２２）の変数行列及びそれらの転置行列で置き換えることにより、線形行列不等式として記述した前記式（２０）について、正定行列Ｖ，Ｗが存在する条件の下で、パラメータγが最小となるように、係数行列Ａ0＄，Ｂ0＄，Ｃ0＄，Ｄ0＄を決定する。 That is, the parameter determination unit 23 is the elements of the above equation (20), Π ₁ ^T P Π ₁ , Π ₁ ^T P Aall Π ₁ , Call Π ₁ , Π ₁ ^T P Ball, and their transposes Π ₁ ^T Aall ^T P Π ₁ , Π. For the equation (20) described as a linear matrix inequality by replacing ₁ ^T Call ^T and Ball ^T PΠ ₁ with the variable matrices of the equations (21) and (22) and their transposed matrices, the canonical matrix V Under the condition that W is present, the coefficient matrices A0 $, B0 $, C0 $, and D0 $ are determined so that the parameter γ is minimized.

決定された係数行列Ａ0＄，Ｂ0＄，Ｃ0＄，Ｄ0＄から、制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0は、以下の式で計算される。

From the determined coefficient matrix A0 $, B0 $, C0 $, D0 $, the coefficient matrix A0, B0, C0, D0 of the controller H0 is calculated by the following equation.

ただし、Ｍ，Ｎは、以下の式を満たす任意の正則行列である。

However, M and N are arbitrary regular matrices that satisfy the following equations.

このように、前記式（１５）と等価である前記式（１６）は、全体システムＨallのＨ∞性能を特徴づける線形行列不等式として記述することができることから、全体システムＨallのＨ∞ノルムを最小化する制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定することができる。 Thus, the equation (16), which is equivalent to the equation (15), can be described as a linear matrix inequality that characterizes the H∞ performance of the overall system Hall, and thus minimizes the H∞ norm of the overall system Hall. The coefficient matrix A0, B0, C0, D0 of the controller H0 to be converted can be determined.

以上のように、本発明の実施形態の制御器設計装置１によれば、制御対象Ｇｐを、直達項を持たない状態空間モデルで表現することを前提に、制御器Ｈのパラメータを、全体システムＨallのＨ∞ノルムを評価関数として、制御器Ｈのゲイン、及び制御対象Ｇｐと目標システムＧｔとの間の出力追従誤差を定量化し、出力追従誤差をある程度許容して、ゲインを小さく抑えるように、制御器Ｈのパラメータを決定する。 As described above, according to the controller design device 1 of the embodiment of the present invention, the parameters of the controller H are set to the entire system on the premise that the controlled object Gp is represented by a state space model having no direct term. Using Hall's H∞ norm as an evaluation function, quantify the gain of the controller H and the output tracking error between the controlled target Gp and the target system Gt, allow the output tracking error to some extent, and keep the gain small. , Determine the parameters of the controller H.

具体的には、同定部２０は、音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に基づいて、制御対象Ｇｐ及び目標システムＧｔを結合した結合システムＧｗの係数行列Ａｗ，Ｂｗ，Ｃｗを求め、結合システムＧｗの係数行列Ａｗ，Ｂｗ，Ｃｗから、制御対象Ｇｐの係数行列Ａｗ，Ｂｐ，Ｃｗを求め、目標システムＧｔの係数行列Ａｗ，Ｂｔ，Ｃｗを求めることで、制御対象Ｇｐ及び目標システムＧｔを同定する。 Specifically, the identification unit 20 is a coefficient of the coupling system Gw in which the controlled target Gp and the target system Gt are coupled based on the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2. Find the matrices Aw, Bw, Cw, find the coefficient matrices Aw, Bp, Cw of the control target Gp from the coefficient matrices Aw, Bw, Cw of the coupling system Gw, and find the coefficient matrices Aw, Bt, Cw of the target system Gt. In, the controlled target Gp and the target system Gt are identified.

遅延時間分離部２１は、音響信号u1-1～u1-m1,y1,y2に基づいて、制御対象Ｇｐの遅延時間を求め、制御対象Ｇｐから遅延時間を分離し、制御対象Ｇｐ＃を求める。また、遅延時間分離部２１は、音響信号u2-1～u2-m2,y1,y2に基づいて、目標システムＧｔの遅延時間を求め、目標システムＧｔから遅延時間を分離し、目標システムＧｔ＃を求める The delay time separation unit 21 obtains the delay time of the control target Gp based on the acoustic signals u1-1 to u1-m1, y1, y2, separates the delay time from the control target Gp, and obtains the control target Gp #. Further, the delay time separation unit 21 obtains the delay time of the target system Gt based on the acoustic signals u2-1 to u2-m2, y1, y2, separates the delay time from the target system Gt, and sets the target system Gt #. Ask

システム構成部２２は、制御器Ｈ0、制御対象Ｇｐ＃、目標システムＧｔ＃、減算器４７及びバランス係数αを用いて、入力信号ｕを音響信号u2-1～u2-m2とし、出力信号y1,y2を、制御器Ｈ0のゲインにより増幅される信号、及び制御対象Ｇｐ＃と目標システムＧｔ＃との間の出力追従誤差とする全体システムＨallを構成する。 The system configuration unit 22 uses the controller H0, the controlled object Gp #, the target system Gt #, the subtractor 47, and the balance coefficient α to set the input signal u as the acoustic signal u2-1 to u2-m2, and the output signal y1, It constitutes an overall system Hall in which y2 is a signal amplified by the gain of the controller H0 and an output tracking error between the controlled target Gp # and the target system Gt #.

パラメータ決定部２３は、バランス係数αを用い、全体システムＨallのＨ∞ノルムを評価関数として、制御器Ｈのゲイン、及び制御対象Ｇｐと目標システムＧｔとの間の出力追従誤差を定量化し、全体システムＨallのＨ∞ノルムを最小化する制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0を決定する。 The parameter determination unit 23 uses the balance coefficient α to quantify the gain of the controller H and the output tracking error between the control target Gp and the target system Gt using the H∞ norm of the entire system Hall as an evaluation function, and the whole is Determine the coefficient matrices A0, B0, C0, D0 of the controller H0 that minimizes the H∞ norm of the system Hall.

遅延時間補償部２４は、制御対象Ｇｐ及び目標システムＧｔから分離した遅延時間に基づいて、遅延器Ｒｄの遅延時間Ｒd1～Ｒdm2及び遅延器Ｒｃの遅延時間Ｒc1～Ｒcm1を求める。 The delay time compensation unit 24 obtains the delay times Rd1 to Rdm2 of the delay device Rd and the delay times Rc1 to Rcm1 of the delay device Rc based on the delay time separated from the controlled target Gp and the target system Gt.

これにより、制御器Ｈのパラメータを、制御器Ｈのゲインを指標とする最適化問題として決定することができ、制御器Ｈの制御性能の低下をある程度許容することで、制御器Ｈのゲインを小さくすることができる。特に、従来技術では制御器Ｈのゲインを十分に小さくすることができない場合に、本発明の実施形態では、制御器Ｈの設計に多少の誤差は生じるが、ゲインを小さくすることができる。したがって、制御器Ｈのゲインを陽に抑制するため、系の外乱及び制御器Ｈ自身の摂動に対して高い頑健性を持つ制御器Ｈを設計することが可能となる。 As a result, the parameter of the controller H can be determined as an optimization problem using the gain of the controller H as an index, and the gain of the controller H can be determined by allowing a certain degree of deterioration in the control performance of the controller H. It can be made smaller. In particular, when the gain of the controller H cannot be sufficiently reduced by the prior art, the gain can be reduced in the embodiment of the present invention, although some errors occur in the design of the controller H. Therefore, in order to positively suppress the gain of the controller H, it is possible to design the controller H having high robustness against the disturbance of the system and the perturbation of the controller H itself.

〔制御部１０の他の構成例〕
次に、制御器設計装置１に備えた制御部１０の他の構成例について説明する。第１構成例は、音響信号を複数の周波数帯域に分割した場合の帯域毎に、制御器設計プログラムの処理にて制御器Ｈを設計することで、制御器Ｈのパラメータを決定する例である。第２構成例は、音響信号を複数の周波数帯域に分割し、さらに音響信号のサンプルの間引きを行った場合の帯域毎に、制御器設計プログラムの処理にて制御器Ｈを設計することで、制御器Ｈのパラメータを決定する例である。 [Other Configuration Examples of Control Unit 10]
Next, another configuration example of the control unit 10 provided in the controller design device 1 will be described. The first configuration example is an example in which the parameters of the controller H are determined by designing the controller H by the processing of the controller design program for each band when the acoustic signal is divided into a plurality of frequency bands. .. In the second configuration example, the acoustic signal is divided into a plurality of frequency bands, and the controller H is designed by processing the controller design program for each band when the sample of the acoustic signal is thinned out. This is an example of determining the parameters of the controller H.

（第１構成例：帯域分割後の音響信号にて制御器Ｈを設計）
まず、制御部１０の他の例として第１構成例について説明する。前述のとおり、第１構成例は、音響信号を複数の周波数帯域に分割した場合の帯域毎に、制御器設計プログラムの処理にて制御器Ｈを設計することで、制御器Ｈのパラメータを決定する例である。 (First configuration example: Controller H is designed with the acoustic signal after band division)
First, a first configuration example will be described as another example of the control unit 10. As described above, in the first configuration example, the parameters of the controller H are determined by designing the controller H by the processing of the controller design program for each band when the acoustic signal is divided into a plurality of frequency bands. This is an example.

図１３は、制御部１０の他の例（第１構成例）を示すブロック図である。この制御部１０－１は、帯域通過フィルタ５０、及び図４に示した制御部１０を備えている。 FIG. 13 is a block diagram showing another example (first configuration example) of the control unit 10. The control unit 10-1 includes a bandpass filter 50 and a control unit 10 shown in FIG.

帯域通過フィルタ５０は、音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を入力し、音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に対し、所定のカットオフ周波数にて帯域通過フィルタ処理を行う。そして、帯域通過フィルタ５０は、帯域通過フィルタ処理後の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を制御部１０に出力する。 The band passage filter 50 inputs acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2, and acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1. Band-passing filter processing is performed on, y2 at a predetermined cutoff frequency. Then, the bandpass filter 50 outputs the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 after the bandpass filter processing to the control unit 10.

制御部１０は、帯域通過フィルタ５０から帯域通過フィルタ処理後の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を入力する。そして、制御部１０は、帯域通過フィルタ処理後の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に基づいて、制御器設計プログラムの処理にて制御器Ｈを設計することで、制御器Ｈのパラメータを決定して出力する。 The control unit 10 inputs the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 after the bandpass filter processing from the bandpass filter 50. Then, the controller 10 controls the controller H in the processing of the controller design program based on the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 after the bandpass filter processing. By designing, the parameters of the controller H are determined and output.

これにより、制御部１０－１にて、所定周波数の帯域に対する制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0、遅延器Ｒｄの遅延時間Ｒd1～Ｒdm2及び遅延器Ｒｃの遅延時間Ｒc1～Ｒcm1が決定される。 As a result, in the control unit 10-1, the coefficient matrices A0, B0, C0, D0 of the controller H0, the delay times Rd1 to Rdm2 of the delay device Rd, and the delay times Rc1 to Rcm1 of the delay device Rc for the band of the predetermined frequency are changed. It is determined.

図１４は、図１３に示した第１構成例の制御部１０－１により設計された制御器Ｈの適用例１を示すブロック図である。このシステムは、音響信号の周波数帯域を２分割した例であり、帯域通過フィルタ５１－１，５１－２及び制御器Ｈ－１ａ，Ｈ－１ｂを備えている。 FIG. 14 is a block diagram showing application example 1 of the controller H designed by the control unit 10-1 of the first configuration example shown in FIG. This system is an example in which the frequency band of an acoustic signal is divided into two, and includes band pass filters 51-1 and 51-2 and controllers H-1a and H-1b.

尚、図１４は、音響信号を図示しないスピーカ１０１－１～１０１－ｍ１へ出力する際には、Ｄ／Ａコンバータ及びアンプを用いるが、これらの基本的な構成部は自明なものとして省略してある。後述する図１５、図１７及び図１８についても同様である。 In FIG. 14, a D / A converter and an amplifier are used when the acoustic signal is output to the speakers 101-1 to 101-m1 (not shown), but these basic components are omitted as they are obvious. There is. The same applies to FIGS. 15, 17, and 18 described later.

帯域通過フィルタ５１－１は、図１３に示した帯域通過フィルタ５０に対応し、音響信号u2-1～u2-m2を入力し、音響信号u2-1～u2-m2に対し、予め設定された第１のカットオフ周波数にて帯域通過フィルタの処理を行う。そして、帯域通過フィルタ５１－１は、第１のカットオフ周波数による帯域通過フィルタ処理後の音響信号u2-1～u2-m2を制御器Ｈ－１ａに出力する。 The bandpass filter 51-1 corresponds to the bandpass filter 50 shown in FIG. 13, and inputs acoustic signals u2-1 to u2-m2, and is preset for the acoustic signals u2-1 to u2-m2. The bandpass filter is processed at the first cutoff frequency. Then, the bandpass filter 51-1 outputs the acoustic signals u2-1 to u2-m2 after the bandpass filter processing by the first cutoff frequency to the controller H-1a.

制御器Ｈ－１ａは、図１３に示した制御部１０－１により設計され、決定されたパラメータにて構成された遅延器Ｒｄａ、制御器Ｈ0ａ及び遅延器Ｒｃａを備えている。遅延器Ｒｄａ、制御器Ｈ0ａ及び遅延器Ｒｃａは、図１３に示した帯域通過フィルタ５０の周波数として、帯域通過フィルタ５１－１と同じ第１のカットオフ周波数を用いた場合に得られたパラメータにて構成される。 The controller H-1a includes a delay device Rda, a controller H0a, and a delay device Rca designed by the control unit 10-1 shown in FIG. 13 and composed of the determined parameters. The delay device Rda, the controller H0a, and the delay device Rca have the parameters obtained when the same first cutoff frequency as the band pass filter 51-1 is used as the frequency of the band pass filter 50 shown in FIG. It is composed of.

制御器Ｈ－１ａは、帯域通過フィルタ５１－１から第１のカットオフ周波数による帯域通過フィルタ処理後の音響信号u2-1～u2-m2を入力し、図５に示した制御器Ｈと同じ処理を行い、音響信号u1a-1～u1a-m1を生成して出力する。これにより、第１のカットオフ周波数の帯域における音響信号u1a-1～u1a-m1が出力される。 The controller H-1a inputs the acoustic signals u2-1 to u2-m2 after the bandpass filter processing by the first cutoff frequency from the bandpass filter 51-1, and is the same as the controller H shown in FIG. Processing is performed to generate and output acoustic signals u1a-1 to u1a-m1. As a result, the acoustic signals u1a-1 to u1a-m1 in the band of the first cutoff frequency are output.

帯域通過フィルタ５１－２は、図１３に示した帯域通過フィルタ５０に対応し、音響信号u2-1～u2-m2を入力し、音響信号u2-1～u2-m2に対し、予め設定された第２のカットオフ周波数（第１のカットオフ周波数とは異なる周波数）にて帯域通過フィルタ処理を行う。そして、帯域通過フィルタ５１－２は、第２のカットオフ周波数による帯域通過フィルタ処理後の音響信号u2-1～u2-m2を制御器Ｈ－１ｂに出力する。 The bandpass filter 51-2 corresponds to the bandpass filter 50 shown in FIG. 13, and inputs acoustic signals u2-1 to u2-m2, and is preset for the acoustic signals u2-1 to u2-m2. The passband filter processing is performed at the second cutoff frequency (frequency different from the first cutoff frequency). Then, the bandpass filter 51-2 outputs the acoustic signals u2-1 to u2-m2 after the bandpass filter processing by the second cutoff frequency to the controller H-1b.

制御器Ｈ－１ｂは、図１３に示した制御部１０－１により設計され、決定されたパラメータにて構成された遅延器Ｒｄｂ、制御器Ｈ0ｂ及び遅延器Ｒｃｂを備えている。遅延器Ｒｄｂ、制御器Ｈ0ｂ及び遅延器Ｒｃｂは、図１３に示した帯域通過フィルタ５０の周波数として、帯域通過フィルタ５１－２と同じ第２のカットオフ周波数を用いた場合に得られたパラメータにて構成される。 The controller H-1b includes a delay device Rdb, a controller H0b, and a delay device Rcb designed by the control unit 10-1 shown in FIG. 13 and composed of the determined parameters. The delay device Rdb, the controller H0b, and the delay device Rcb have the parameters obtained when the same second cutoff frequency as the band pass filter 51-2 is used as the frequency of the band pass filter 50 shown in FIG. It is composed of.

制御器Ｈ－１ｂは、帯域通過フィルタ５１－２から第２のカットオフ周波数による帯域通過フィルタ処理後の音響信号u2-1～u2-m2を入力し、図５に示した制御器Ｈと同じ処理を行い、音響信号u1b-1～u1b-m1を生成して出力する。これにより、第２のカットオフ周波数の帯域における音響信号u1b-1～u1b-m1が出力される。 The controller H-1b inputs the acoustic signals u2-1 to u2-m2 after the bandpass filter processing by the second cutoff frequency from the bandpass filter 51-2, and is the same as the controller H shown in FIG. Processing is performed to generate and output acoustic signals u1b-1 to u1b-m1. As a result, the acoustic signals u1b-1 to u1b-m1 in the band of the second cutoff frequency are output.

そして、第１のカットオフ周波数の帯域における音響信号u1a-1～u1a-m1と第２のカットオフ周波数の帯域における音響信号u1b-1～u1b-m1とがそれぞれ加算され、加算された音響信号がスピーカ１０１－１～１０１－ｍ１へ出力される。この場合、音響信号u1a-1～u1a-m1のみがスピーカ１０１－１～１０１－ｍ１へ出力されるようにしてもよいし、音響信号u1b-1～u1b-m1のみがスピーカ１０１－１～１０１－ｍ１へ出力されるようにしてもよい。スピーカ１０１－１～１０１－ｍ１が再生周波数帯域の異なる複数のスピーカユニットを備える場合には、音響信号u1a-1～u1a-m1が、当該音響信号u1a-1～u1a-m1の通過帯域を再生周波数帯域に含むスピーカユニットへ出力されるようにしてもよい。また、音響信号u1b-1～u1b-m1が、当該音響信号u1b-1～u1b-m1の通過帯域を再生周波数帯域に含むスピーカユニットへ出力されるようにしてもよい。 Then, the acoustic signals u1a-1 to u1a-m1 in the band of the first cutoff frequency and the acoustic signals u1b-1 to u1b-m1 in the band of the second cutoff frequency are added, and the added acoustic signals are added. Is output to the speakers 101-1 to 101-m1. In this case, only the acoustic signals u1a-1 to u1a-m1 may be output to the speakers 101-1 to 101-m1, or only the acoustic signals u1b-1 to u1b-m1 may be output to the speakers 101-1 to 101. It may be output to −m1. When the speakers 101-1 to 101-m1 include a plurality of speaker units having different reproduction frequency bands, the acoustic signals u1a-1 to u1a-m1 reproduce the pass band of the acoustic signals u1a-1 to u1a-m1. It may be output to the speaker unit included in the frequency band. Further, the acoustic signals u1b-1 to u1b-m1 may be output to the speaker unit including the pass band of the acoustic signals u1b-1 to u1b-m1 in the reproduction frequency band.

図１５は、図１３に示した第１構成例の制御部１０－１により設計された制御器Ｈの適用例２を示すブロック図である。このシステムは、音響信号の周波数帯域を２分割した例であり、帯域通過フィルタ５１－１，５１－２、制御器Ｈ－１ａ及び他の制御器を備えている。このシステムには、本発明の実施形態にて設計された制御器Ｈ－１ａと従来技術にて設計された他の制御器とが混在している。 FIG. 15 is a block diagram showing an application example 2 of the controller H designed by the control unit 10-1 of the first configuration example shown in FIG. This system is an example in which the frequency band of an acoustic signal is divided into two, and includes band pass filters 51-1 and 51-2, a controller H-1a, and another controller. In this system, the controller H-1a designed according to the embodiment of the present invention and another controller designed by the prior art are mixed.

帯域通過フィルタ５１－１，５１－２及び制御器Ｈ－１ａは、図１４に示した帯域通過フィルタ５１－１，５１－２及び制御器Ｈ－１ａと同じであるから、ここでは説明を省略する。 Since the bandpass filters 51-1 and 51-2 and the controller H-1a are the same as the bandpass filters 51-1 and 51-2 and the controller H-1a shown in FIG. 14, the description thereof is omitted here. do.

他の制御器は、前述の非特許文献１等の従来の制御器設計装置により設計された制御器であり、帯域通過フィルタ５１－２から第２のカットオフ周波数による帯域通過フィルタ処理後の音響信号u2-1～u2-m2を入力し、音響信号u1’-1～u1’-m1を生成して出力する。 The other controller is a controller designed by the conventional controller design device such as Non-Patent Document 1 described above, and the sound after the bandpass filter processing by the bandpass filter 51-2 to the second cutoff frequency is performed. The signals u2-1 to u2-m2 are input, and the acoustic signals u1'-1 to u1'-m1 are generated and output.

これにより、制御器Ｈ－１ａから第１のカットオフ周波数の帯域における音響信号u1a-1～u1a-m1が出力され、他の制御器から第２のカットオフ周波数の帯域における音響信号u1’-1～u1’-m1が出力される。そして、音響信号u1a-1～u1a-m1と音響信号u1’-1～u1’-m1とが加算され、加算された音響信号がスピーカ１０１－１～１０１－ｍ１へ出力される。この場合、音響信号u1a-1～u1a-m1及び音響信号u1’-1～u1’-m1のうちのいずれかがスピーカ１０１－１～１０１－ｍ１へ出力されるようにしてもよい。スピーカ１０１－１～１０１－ｍ１が再生周波数帯域の異なる複数のスピーカユニットを備える場合には、音響信号u1a-1～u1a-m1が、当該音響信号u1a-1～u1a-m1の通過帯域を再生周波数帯域に含むスピーカユニットへ出力されるようにしてもよい。また、音響信号u1’-1～u1’-m1が、当該音響信号u1’-1～u1’-m1の通過帯域を再生周波数帯域に含むスピーカユニットへ出力されるようにしてもよい。 As a result, the acoustic signals u1a-1 to u1a-m1 in the band of the first cutoff frequency are output from the controller H-1a, and the acoustic signals u1'-in the band of the second cutoff frequency are output from the other controllers. 1 to u1'-m1 are output. Then, the acoustic signals u1a-1 to u1a-m1 and the acoustic signals u1'-1 to u1'-m1 are added, and the added acoustic signals are output to the speakers 101-1 to 101-m1. In this case, any one of the acoustic signals u1a-1 to u1a-m1 and the acoustic signals u1'-1 to u1'-m1 may be output to the speakers 101-1 to 101-m1. When the speakers 101-1 to 101-m1 include a plurality of speaker units having different reproduction frequency bands, the acoustic signals u1a-1 to u1a-m1 reproduce the pass band of the acoustic signals u1a-1 to u1a-m1. It may be output to the speaker unit included in the frequency band. Further, the acoustic signals u1'-1 to u1'-m1 may be output to the speaker unit including the pass band of the acoustic signals u1'-1 to u1'-m1 in the reproduction frequency band.

以上のように、本発明の実施形態による制御器設計装置１の制御部１０－１によれば、音響信号を複数の周波数帯域に分割した場合の帯域毎に、制御器設計プログラムの処理にて制御器Ｈを設計することで、制御器Ｈのパラメータを決定するようにした。 As described above, according to the control unit 10-1 of the controller design device 1 according to the embodiment of the present invention, in the processing of the controller design program for each band when the acoustic signal is divided into a plurality of frequency bands. By designing the controller H, the parameters of the controller H are determined.

これにより、前述の制御部１０と同様の効果を奏する。すなわち、帯域通過フィルタ５０にてフィルタ処理がなされた所定の周波数帯域の音響信号について、制御器Ｈの制御性能の低下がある程度許容され、設計に多少の誤差は生じるが、ゲインを小さくすることができる。 As a result, the same effect as that of the above-mentioned control unit 10 is obtained. That is, with respect to the acoustic signal in the predetermined frequency band filtered by the band pass filter 50, a decrease in the control performance of the controller H is allowed to some extent, and a slight error occurs in the design, but the gain can be reduced. can.

制御部１０－１により第１のカットオフ周波数の帯域にて設計された制御器Ｈ－１ａは、帯域通過フィルタ５１－１にてフィルタ処理がなされた第１のカットオフ周波数の音響信号u2-1～u2-m2に対し、ゲインが陽に抑制され、系の外乱及び制御器Ｈ－１ａ自身の摂動に対して高い頑健性を得ることが可能となる。 The controller H-1a designed by the control unit 10-1 in the band of the first cutoff frequency is an acoustic signal u2- of the first cutoff frequency filtered by the band pass filter 51-1. The gain is positively suppressed for 1 to u2-m2, and it becomes possible to obtain high robustness against the disturbance of the system and the perturbation of the controller H-1a itself.

尚、図１４及び図１５は、音響信号の周波数帯域を２分割した例であるが、音響信号の周波数帯域を３以上の周波数帯域に分割するようにしてもよい。この場合の制御器Ｈ－１ａ等は、周波数帯域毎に、当該周波数帯域の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に基づいて設計される。 Although FIGS. 14 and 15 show an example in which the frequency band of the acoustic signal is divided into two, the frequency band of the acoustic signal may be divided into three or more frequency bands. In this case, the controller H-1a and the like are designed for each frequency band based on the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 in the frequency band.

また、図１３に示した制御部１０－１は、測定した音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を入力し、帯域通過フィルタ５０にて、所定のカットオフ周波数による帯域通過フィルタの処理を行うようにした。これに対し、所定のカットオフ周波数により帯域制限された音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2が測定された場合、制御部１０－１は、当該所定のカットオフ周波数により帯域制限された音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を直接入力するようにしてもよい。この場合、制御部１０－１は、帯域通過フィルタ５０を備えていない。 Further, the control unit 10-1 shown in FIG. 13 inputs the measured acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2, and is predetermined by the bandpass filter 50. Changed to process the bandpass filter by the cutoff frequency. On the other hand, when the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 whose band is limited by the predetermined cutoff frequency are measured, the control unit 10-1 performs the predetermined. The acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 whose band is limited by the cutoff frequency of the above may be directly input. In this case, the control unit 10-1 does not include the bandpass filter 50.

（第２構成例：帯域分割及び間引き後の音響信号にて制御器Ｈを設計）
次に、第２構成例について説明する。前述のとおり、第２構成例は、音響信号を複数の周波数帯域に分割し、さらに音響信号のサンプルの間引きを行った場合の帯域毎に、制御器設計プログラムの処理にて制御器Ｈを設計することで、制御器Ｈのパラメータを決定する例である。 (Second configuration example: Controller H is designed based on the acoustic signal after band division and thinning)
Next, a second configuration example will be described. As described above, in the second configuration example, the controller H is designed by processing the controller design program for each band when the acoustic signal is divided into a plurality of frequency bands and the sample of the acoustic signal is thinned out. This is an example of determining the parameter of the controller H.

図１６は、制御部１０の他の例（第２構成例）を示すブロック図である。この制御部１０－２は、帯域通過フィルタ５０、間引き部５２及び図４に示した制御部１０を備えている。 FIG. 16 is a block diagram showing another example (second configuration example) of the control unit 10. The control unit 10-2 includes a bandpass filter 50, a thinning unit 52, and a control unit 10 shown in FIG.

帯域通過フィルタ５０は、図１３に示した帯域通過フィルタ５０と同じであるからここでは説明を省略する。 Since the band-passing filter 50 is the same as the band-passing filter 50 shown in FIG. 13, the description thereof is omitted here.

間引き部５２は、帯域通過フィルタ５０から帯域通過フィルタ処理後の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を入力する。そして、間引き部５２は、帯域通過フィルタ５０のカットオフ周波数に応じて、帯域通過フィルタ処理後の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に対しサンプルの間引き（ダウンサンプリング）を行う。間引き部５２は、間引き後の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を制御部１０に出力する。 The thinning unit 52 inputs the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 after the bandpass filter processing from the bandpass filter 50. Then, the thinning unit 52 is a sample for the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 after the bandpass filter processing according to the cutoff frequency of the bandpass filter 50. Perform thinning (downsampling). The thinning unit 52 outputs the thinned-out acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1 and y2 to the control unit 10.

これにより、間引き後の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2のサンプリング周波数が低い場合、制御部１０及び後述する制御器Ｈ－１ｃにおけるモデルの次数を低く抑えることができ、信号処理量を削減することができる。 As a result, when the sampling frequency of the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 after thinning is low, the order of the model in the control unit 10 and the controller H-1c described later is changed. It can be kept low and the amount of signal processing can be reduced.

制御部１０は、間引き部５２から間引き後の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を入力し、図１３に示した制御部１０と同じ処理を行う。 The control unit 10 inputs the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 after thinning from the thinning unit 52, and performs the same processing as the control unit 10 shown in FIG. ..

これにより、制御部１０－２にて、所定周波数の帯域に対する制御器Ｈ0の係数行列Ａ0，Ｂ0，Ｃ0，Ｄ0、遅延器Ｒｄの遅延時間Ｒd1～Ｒdm2及び遅延器Ｒｃの遅延時間Ｒc1～Ｒcm1が決定される。 As a result, in the control unit 10-2, the coefficient matrices A0, B0, C0, D0 of the controller H0, the delay times Rd1 to Rdm2 of the delay device Rd, and the delay times Rc1 to Rcm1 of the delay device Rc for the band of the predetermined frequency are changed. It is determined.

図１７は、図１６に示した第２構成例の制御部１０－２により設計された制御器Ｈの適用例１を示すブロック図である。このシステムは、音響信号の周波数帯域を２分割した例であり、帯域通過フィルタ５１－１，５１－２、間引き部５３－１，５３－２、制御器Ｈ－１ｃ、他の制御器及び増幅器５４－１，５４－２を備えている。このシステムには、本発明の実施形態にて設計された制御器Ｈ－１ｃと従来技術にて設計された他の制御器とが混在している。尚、本発明の実施形態にて設計された制御器Ｈ－１ｃと制御器Ｈ－１ｄとを並べて構成するようにしてもよい。 FIG. 17 is a block diagram showing application example 1 of the controller H designed by the control unit 10-2 of the second configuration example shown in FIG. This system is an example in which the frequency band of an acoustic signal is divided into two, a band pass filter 51-1, 51-2, a thinning unit 53-1, 53-2, a controller H-1c, another controller, and an amplifier. It is equipped with 54-1 and 54-2. In this system, the controller H-1c designed according to the embodiment of the present invention and another controller designed by the prior art are mixed. The controller H-1c and the controller H-1d designed in the embodiment of the present invention may be arranged side by side.

帯域通過フィルタ５１－１，５１－２は、図１４に示した帯域通過フィルタ５１－１，５１－２と同じであるから、ここでは説明を省略する。 Since the band-passing filters 51-1 and 51-2 are the same as the band-passing filters 51-1 and 51-2 shown in FIG. 14, the description thereof will be omitted here.

間引き部５３－１は、帯域通過フィルタ５１－１から第１のカットオフ周波数による帯域通過フィルタ処理後の音響信号u2-1～u2-m2を入力する。そして、間引き部５３－１は、帯域通過フィルタ５１－１の第１のカットオフ周波数に応じて、図１６に示した間引き部５２と同様に、帯域通過フィルタ処理後の音響信号u2-1～u2-m2に対しサンプルの間引きを行う。そして、間引き部５３－１は、間引き後の音響信号u2-1～u2-m2を制御器Ｈ－１ｃに出力する。 The thinning unit 53-1 inputs the acoustic signals u2-1 to u2-m2 after the bandpass filter processing by the first cutoff frequency from the bandpass filter 51-1. Then, the thinning section 53-1 responds to the first cutoff frequency of the bandpass filter 51-1, and similarly to the thinning section 52 shown in FIG. 16, the acoustic signals u2-1 to the bandpass filter processing are performed. Thin out the sample for u2-m2. Then, the thinning unit 53-1 outputs the acoustic signals u2-1 to u2-m2 after the thinning out to the controller H-1c.

制御器Ｈ－１ｃは、図１６に示した制御部１０－２により設計され、決定されたパラメータにて構成された遅延器Ｒｄ、制御器Ｈ0及び遅延器Ｒｃを備えている。遅延器Ｒｄ、制御器Ｈ0及び遅延器Ｒｃは、図１６に示した帯域通過フィルタ５０のカットオフ周波数及び間引き部５２の間引き間隔として、帯域通過フィルタ５１－１と同じ第１のカットオフ周波数、及び間引き部５３－１と同じ間引き間隔を用いた場合に得られたパラメータにて構成される。 The controller H-1c includes a delay device Rd, a controller H0, and a delay device Rc designed by the control unit 10-2 shown in FIG. 16 and composed of the determined parameters. The delayr Rd, the controller H0, and the delayer Rc have the same first cutoff frequency as the bandpass filter 51-1 as the cutoff frequency of the bandpass filter 50 and the thinning interval of the thinning section 52 shown in FIG. It is composed of the parameters obtained when the same thinning interval as that of the thinning portion 53-1 is used.

制御器Ｈ－１ｃは、間引き部５３－１から間引き後の音響信号u2-1～u2-m2を入力し、図５に示した制御器Ｈと同じ処理を行い、増幅前の音響信号u1-1～u1-m1を生成して増幅器５４－１に出力する。 The controller H-1c inputs the thinned-out acoustic signals u2-1 to u2-m2 from the thinned-out unit 53-1 and performs the same processing as the controller H shown in FIG. 5, and the acoustic signal u1- before amplification is performed. 1 to u1-m1 are generated and output to the amplifier 54-1.

増幅器５４－１は、制御器Ｈ－１ｃから増幅前の音響信号u1-1～u1-m1を入力し、間引き部５３－１の間引き量に応じて、増幅前の音響信号u1-1～u1-m1を増幅し、音響信号u1c-1～u1c-m1を生成して出力する。 The amplifier 54-1 inputs the acoustic signals u1-1 to u1-m1 before amplification from the controller H-1c, and the acoustic signals u1-1 to u1 before amplification are subjected to the thinning amount of the thinning unit 53-1. -Amplifies m1 to generate and output acoustic signals u1c-1 to u1c-m1.

間引き部５３－２は、帯域通過フィルタ５１－２から第２のカットオフ周波数による帯域通過フィルタ処理後の音響信号u2-1～u2-m2を入力する。そして、間引き部５３－２は、帯域通過フィルタ５１－２の第２のカットオフ周波数に応じて、図１６に示した間引き部５２と同様に、帯域通過フィルタ処理後の音響信号u2-1～u2-m2に対しサンプルの間引きを行う。そして、間引き部５３－２は、間引き後の音響信号u2-1～u2-m2を他の制御器に出力する。 The thinning unit 53-2 inputs the acoustic signals u2-1 to u2-m2 after the bandpass filter processing by the second cutoff frequency from the bandpass filter 51-2. Then, the thinning unit 53-2 responds to the second cutoff frequency of the bandpass filter 51-2, and similarly to the thinning unit 52 shown in FIG. 16, the acoustic signals u2-1 to after the bandpass filter processing Thin out the sample for u2-m2. Then, the thinning unit 53-2 outputs the acoustic signals u2-1 to u2-m2 after the thinning out to another controller.

他の制御器は、前述の非特許文献１等の従来の制御器設計装置により設計された制御器であり、間引き部５３－２から間引き後の音響信号u2-1～u2-m2を入力し、増幅前の音響信号u1-1～u1-m1を生成して出力する。 The other controller is a controller designed by the conventional controller design device such as the above-mentioned Non-Patent Document 1, and inputs the acoustic signals u2-1 to u2-m2 after thinning from the thinning unit 53-2. , Generates and outputs the acoustic signals u1-1 to u1-m1 before amplification.

増幅器５４－２は、他の制御器から増幅前の音響信号u1-1～u1-m1を入力し、間引き部５３－２の間引き量に応じて、増幅前の音響信号u1-1～u1-m1を増幅し、音響信号u1’-1～u1’-m1を生成して出力する。 The amplifier 54-2 inputs the acoustic signals u1-1 to u1-m1 before amplification from another controller, and the acoustic signals u1-1 to u1- before amplification are input according to the thinning amount of the thinning unit 53-2. Amplifies m1 and generates and outputs acoustic signals u1'-1 to u1'-m1.

これにより、増幅器５４－１から第１のカットオフ周波数の帯域における間引き及び増幅後の音響信号u1c-1～u1c-m1が出力され、増幅器５４－２から第２のカットオフ周波数の帯域における間引き及び増幅後の音響信号u1’-1～u1’c-m1が出力される。そして、音響信号u1c-1～u1c-m1と音響信号u1’-1～u1’-m1とが加算され、加算された音響信号がスピーカ１０１－１～１０１－ｍ１へ出力される。この場合、音響信号u1c-1～u1c-m1及び音響信号u1’-1～u1’-m1のいずれかがスピーカ１０１－１～１０１－ｍ１へ出力されるようにしてもよい。スピーカ１０１－１～１０１－ｍ１が再生周波数帯域の異なる複数のスピーカユニットを備える場合には、音響信号u1a-1～u1a-m1が、当該音響信号u1a-1～u1a-m1の通過帯域を再生周波数帯域に含むスピーカユニットへ出力されるようにしてもよい。また、音響信号u1’-1～u1’-m1が、当該音響信号u1’-1～u1’-m1の通過帯域を再生周波数帯域に含むスピーカユニットへ出力されるようにしてもよい。 As a result, the amplifiers 54-1 are thinned out in the band of the first cutoff frequency and the amplified acoustic signals u1c-1 to u1c-m1 are output, and the amplifiers 54-2 are thinned out in the band of the second cutoff frequency. And the amplified acoustic signals u1'-1 to u1'c-m1 are output. Then, the acoustic signals u1c-1 to u1c-m1 and the acoustic signals u1'-1 to u1'-m1 are added, and the added acoustic signals are output to the speakers 101-1 to 101-m1. In this case, any of the acoustic signals u1c-1 to u1c-m1 and the acoustic signals u1'-1 to u1'-m1 may be output to the speakers 101-1 to 101-m1. When the speakers 101-1 to 101-m1 include a plurality of speaker units having different reproduction frequency bands, the acoustic signals u1a-1 to u1a-m1 reproduce the pass band of the acoustic signals u1a-1 to u1a-m1. It may be output to the speaker unit included in the frequency band. Further, the acoustic signals u1'-1 to u1'-m1 may be output to the speaker unit including the pass band of the acoustic signals u1'-1 to u1'-m1 in the reproduction frequency band.

図１８は、図１６に示した第２構成例の制御部１０－２により設計された制御器Ｈの適用例２を示すブロック図である。このシステムは、音響信号の周波数帯域を２分割した例であり、帯域通過フィルタ５１－１，５１－２、間引き部５３－１，５３－２、制御器Ｈ－１ｃ、他の制御器、補間及びＬＰＦ部５５－１，５５－２、増幅器５４－１，５４－２及び加算器５６を備えている。このシステムには、本発明の実施形態にて設計された制御器Ｈ－１ｃと従来技術にて設計された他の制御器とが混在している。尚、本発明の実施形態にて設計された制御器Ｈ－１ｃと制御器Ｈ－１ｄとを並べて構成するようにしてもよい。 FIG. 18 is a block diagram showing an application example 2 of the controller H designed by the control unit 10-2 of the second configuration example shown in FIG. This system is an example in which the frequency band of an acoustic signal is divided into two, a band pass filter 51-1, 51-2, a thinning unit 53-1, 53-2, a controller H-1c, another controller, and an interpolation. It also includes an LPF unit 55-1, 55-2, amplifiers 54-1, 54-2, and an adder 56. In this system, the controller H-1c designed according to the embodiment of the present invention and another controller designed by the prior art are mixed. The controller H-1c and the controller H-1d designed in the embodiment of the present invention may be arranged side by side.

図１７に示したシステムと図１８に示すシステムとを比較すると、図１８のシステムは、図１７に示したシステムに加え、さらに、制御器Ｈ－１ｃと増幅器５４－１との間に補間及びＬＰＦ部５５－１を備え、他の制御器と増幅器５４－２との間に補間及びＬＰＦ部５５－２を備え、加算器５６を備えている点で相違する。図１８において、図１７と共通する部分には図１７と同一の符号を付し、その詳しい説明は省略する。 Comparing the system shown in FIG. 17 with the system shown in FIG. 18, the system of FIG. 18 further interpolates and interpolates between the controller H-1c and the amplifier 54-1 in addition to the system shown in FIG. The difference is that the LPF unit 55-1 is provided, the interpolation and LPF unit 55-2 are provided between the other controller and the amplifier 54-2, and the adder 56 is provided. In FIG. 18, the parts common to FIG. 17 are designated by the same reference numerals as those in FIG. 17, and detailed description thereof will be omitted.

補間及びＬＰＦ部５５－１は、制御器Ｈ－１ｃから音響信号u1-1～u1-m1を入力し、音響信号u1-1～u1-m1に対し、間引き部５３－１による間引きとは逆の処理である補間（アップサンプリング）を行う。これにより、元のサンプリング周波数に戻すことができる。そして、補間及びＬＰＦ部５５－１は、補間後の音響信号u1-1～u1-m1に対し、帯域通過フィルタ５１－１と同じ第１のカットオフ周波数のＬＰＦ処理を施し、ＬＰＦ処理後の音響信号u1-1～u1-m1を増幅器５４－１に出力する。 The interpolation and LPF unit 55-1 inputs the acoustic signals u1-1 to u1-m1 from the controller H-1c, and the acoustic signals u1-1 to u1-m1 are the opposite of the thinning by the thinning unit 53-1. Interpolation (upsampling), which is the processing of. This makes it possible to return to the original sampling frequency. Then, the interpolation and LPF unit 55-1 performs LPF processing of the first cutoff frequency same as that of the bandpass filter 51-1 on the acoustic signals u1-1 to u1-m1 after the interpolation, and after the LPF processing. The acoustic signals u1-1 to u1-m1 are output to the amplifier 54-1.

これにより、補間処理により生じる折り返し歪は、帯域通過フィルタ５１－１と同じカットオフ周波数を用いる帯域通過フィルタのＬＰＦ処理により削除することができる。 Thereby, the folding distortion caused by the interpolation processing can be removed by the LPF processing of the bandpass filter using the same cutoff frequency as the bandpass filter 51-1.

補間及びＬＰＦ部５５－２は、他の制御器から音響信号u1’-1～u1’-m1を入力し、音響信号u1’-1～u1’-m1に対し、間引き部５３－２による間引きとは逆の補間を行う。これにより、元のサンプリング周波数に戻すことができる。そして、補間及びＬＰＦ部５５－２は、補間後の音響信号u1’-1～u1’-m1に対し、帯域通過フィルタ５１－２と同じ第２のカットオフ周波数のＬＰＦを施し、ＬＰＦ後の音響信号u1’-1～u1’-m1を増幅器５４－２に出力する。 The interpolation and LPF section 55-2 inputs acoustic signals u1'-1 to u1'-m1 from other controllers, and thins out the acoustic signals u1'-1 to u1'-m1 by the thinning section 53-2. Interpolates the opposite of. This makes it possible to return to the original sampling frequency. Then, the interpolation and LPF unit 55-2 applies LPF of the same second cutoff frequency as the bandpass filter 51-2 to the acoustic signals u1'-1 to u1'-m1 after LPF, and after LPF. The acoustic signals u1'-1 to u1'-m1 are output to the amplifier 54-2.

これにより、補間処理により生じる折り返し歪は、間引き部５３－２と同じカットオフ周波数を用いる帯域通過フィルタのＬＰＦ処理により削除することができる。増幅器５４－１，５４－２の処理は、図１７にて説明したとおりである。 Thereby, the folding distortion caused by the interpolation processing can be deleted by the LPF processing of the bandpass filter using the same cutoff frequency as the thinning unit 53-2. The processing of the amplifiers 54-1 and 54-2 is as described with reference to FIG.

加算器５６は、増幅器５４－１から増幅後の音響信号u1c-1～u1c-m1を入力すると共に、増幅器５４－２から増幅後の音響信号u1’-1～u1’-m1を入力し、音響信号u1c-1～u1c-m1及び音響信号u1’-1～u1’-m1を加算し、加算後の音響信号u1d-1～u1d-m1を出力する。 The adder 56 inputs the amplified acoustic signals u1c-1 to u1c-m1 from the amplifier 54-1 and inputs the amplified acoustic signals u1'-1 to u1'-m1 from the amplifier 54-2. The acoustic signals u1c-1 to u1c-m1 and the acoustic signals u1'-1 to u1'-m1 are added, and the added acoustic signals u1d-1 to u1d-m1 are output.

以上のように、本発明の実施形態による制御器設計装置１の制御部１０－２によれば、音響信号を複数の周波数帯域に分割し、さらに音響信号のサンプルの間引きを行った場合の帯域毎に、制御器設計プログラムの処理にて制御器Ｈを設計することで、制御器Ｈのパラメータを決定するようにした。 As described above, according to the control unit 10-2 of the controller design device 1 according to the embodiment of the present invention, the band when the acoustic signal is divided into a plurality of frequency bands and the sample of the acoustic signal is thinned out. The parameters of the controller H are determined by designing the controller H by processing the controller design program each time.

これにより、前述の制御部１０と同様の効果を奏する。すなわち、帯域通過フィルタ５０によりフィルタ処理がなされた所定の周波数帯域の音響信号について、制御器Ｈの制御性能の低下がある程度許容され、設計に多少の誤差は生じるが、ゲインを小さくすることができる。 As a result, the same effect as that of the above-mentioned control unit 10 is obtained. That is, with respect to the acoustic signal in the predetermined frequency band filtered by the band pass filter 50, the control performance of the controller H is allowed to deteriorate to some extent, and the gain can be reduced although some errors occur in the design. ..

制御部１０－２により第１のカットオフ周波数の帯域にて設計された制御器Ｈ－１ｃは、帯域通過フィルタ５１－１及び間引き部５３－１にて生成された第１のカットオフ周波数及び間引き後の音響信号u2-1～u2-m2に対し、ゲインが陽に抑制され、系の外乱及び制御器Ｈ－１ｃ自身の摂動に対して高い頑健性を得ることが可能となる。 The controller H-1c designed by the control unit 10-2 in the band of the first cutoff frequency has the first cutoff frequency generated by the band pass filter 51-1 and the thinning unit 53-1. The gain is positively suppressed for the acoustic signals u2-1 to u2-m2 after thinning, and it becomes possible to obtain high robustness against the disturbance of the system and the perturbation of the controller H-1c itself.

尚、図１７及び図１８は、音響信号の周波数帯域を２分割した例であるが、図１４及び図１５と同様に、音響信号の周波数帯域を３以上の周波数帯域に分割するようにしてもよい。この場合の制御器Ｈ－１ｃ等は、周波数帯域毎に、当該周波数帯域の音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2に基づいて設計される。また、分割した全ての周波数帯域のそれぞれについて、制御部１０－２により設計された制御器Ｈを用いるようにしてもよい。 17 and 18 are examples in which the frequency band of the acoustic signal is divided into two, but similarly to FIGS. 14 and 15, the frequency band of the acoustic signal may be divided into three or more frequency bands. good. In this case, the controller H-1c and the like are designed for each frequency band based on the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 in the frequency band. Further, the controller H designed by the control unit 10-2 may be used for each of all the divided frequency bands.

また、図１６に示した制御部１０－２は、測定した音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を入力し、帯域通過フィルタ５０にて、所定のカットオフ周波数による帯域通過フィルタ処理を行い、間引き部５２にて間引き処理を行うようにした。これに対し、所定のカットオフ周波数かつ間引き処理がなされた音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2が測定された場合、制御部１０－２は、当該音響信号u1-1～u1-m1,u2-1～u2-m2,y1,y2を直接入力するようにしてもよい。この場合、制御部１０－２は、帯域通過フィルタ５０及び間引き部５２を備えていない。 Further, the control unit 10-2 shown in FIG. 16 inputs the measured acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1 and y2, and uses a bandpass filter 50 to determine a predetermined value. Bandpass filter processing based on the cutoff frequency was performed, and thinning processing was performed by the thinning unit 52. On the other hand, when the acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 having a predetermined cutoff frequency and being thinned out are measured, the control unit 10-2 is concerned. The acoustic signals u1-1 to u1-m1, u2-1 to u2-m2, y1, y2 may be directly input. In this case, the control unit 10-2 does not include the bandpass filter 50 and the thinning unit 52.

図５に示した制御器Ｈのハードウェア構成は、制御器設計装置１と同様に、通常のコンピュータを使用することができる。制御器Ｈは、ＣＰＵ、ＲＡＭ等の揮発性の記憶媒体、ＲＯＭ等の不揮発性の記憶媒体、及びインタフェース等を備えたコンピュータによって構成される。制御器Ｈに備えた遅延器Ｒｄ、制御器Ｈ0及び遅延器Ｒｃの各機能は、これらの機能を記述したプログラムをＣＰＵに実行させることによりそれぞれ実現される。これらのプログラムは、前記記憶媒体に格納されており、ＣＰＵに読み出されて実行される。また、これらのプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスク等）、光ディスク（ＣＤ－ＲＯＭ、ＤＶＤ等）、半導体メモリ等の記憶媒体に格納して頒布することもでき、ネットワークを介して送受信することもできる。図１４及び図１５に示した制御器Ｈ－１ａ、図１４に示した制御器Ｈ－１ｂ、並びに図１７及び図１８に示した制御器Ｈ－１ｃについても同様である。 As for the hardware configuration of the controller H shown in FIG. 5, a normal computer can be used as in the controller design device 1. The controller H is composed of a computer including a volatile storage medium such as a CPU and RAM, a non-volatile storage medium such as ROM, and an interface. Each function of the delay device Rd, the controller H0, and the delay device Rc provided in the controller H is realized by causing the CPU to execute a program describing these functions. These programs are stored in the storage medium, read by the CPU, and executed. In addition, these programs can be stored and distributed in storage media such as magnetic disks (floppy (registered trademark) disks, hard disks, etc.), optical disks (CD-ROM, DVD, etc.), semiconductor memories, etc., and can be distributed via a network. You can also send and receive. The same applies to the controller H-1a shown in FIGS. 14 and 15, the controller H-1b shown in FIG. 14, and the controller H-1c shown in FIGS. 17 and 18.

１ 制御器設計装置
１０ 制御部
１１ ＣＰＵ
１２ 記憶部
１３ 記憶装置
１４ 操作／入力部
１５ 表示出力インタフェース部
１６ 通信部
１７ システムバス
２０ 同定部
２１ 遅延時間分離部
２２ システム構成部
２３ パラメータ決定部
２４ 遅延時間補償部（遅延時間算出部）
３０，３１，４２ 遅延器
４０，４３，４４，４５ 乗算器
４１，４６ 加算器
４７ 減算器
５０，５１ 帯域通過フィルタ
５２，５３ 間引き部
５４ 増幅器
５５ 補間及びＬＰＦ部
５６ 加算器
１０１，１０３ スピーカ
１０２ マイクロホン
１０４ 聴取者
１０５ 耳
Ｇｐ 制御対象
Ｇｐ＃ 遅延時間が分離された制御対象（分離制御対象）
Ｈinv 逆システム
Ｇｔ 目標システム
Ｇｔ＃ 遅延時間が分離された目標システム（分離目標システム）
Ｇｗ 結合システム
Ｈall 全体システム
Ｈ 制御器
Ｈ0 制御器（仮制御器）
Ｒｄ 遅延器（第１遅延器）
Ｒｃ 遅延器（第２遅延器） 1 Controller design device 10 Control unit 11 CPU
12 Storage unit 13 Storage device 14 Operation / input unit 15 Display output interface unit 16 Communication unit 17 System bus 20 Identification unit 21 Delay time separation unit 22 System configuration unit 23 Parameter determination unit 24 Delay time compensation unit (delay time calculation unit)
30, 31, 42 Delayer 40, 43, 44, 45 Multiplier 41, 46 Adder 47 Subtractor 50, 51 Band pass filter 52, 53 Thinning unit 54 Amplifier 55 Interpolator and LPF unit 56 Adder 101, 103 Speaker 102 Microphone 104 Listener 105 Ear Gp Control target Gp # Control target with separated delay time (separate control target)
Hinv Inverse system Gt Target system Gt # Target system with separated delay time (separated target system)
Gw coupling system Hall whole system H controller H0 controller (temporary controller)
Rd delayer (first delayer)
Rc delayer (second delayer)

Claims (2)

信号が観察される所定数の制御点と前記所定数の制御点へ信号をそれぞれ提示する所定数の２次音源とを含む制御対象についての逆システムであって、前記所定数の２次音源へ信号を出力する前記逆システムの特性と、前記所定数の制御点へ信号をそれぞれ提示する所定数の仮想音源を含む目標システムの特性とを合わせた制御器を設計する音響信号用の制御器設計装置において、
前記制御器が、第１遅延器、仮制御器及び第２遅延器を備え、
前記制御対象を、前記所定数の２次音源から提示される第１信号と前記所定数の制御点にて観察される第２信号との間の関係について、直達項を持たない状態空間モデルで表し、
前記目標システムを、前記所定数の仮想音源から提示される第３信号と前記所定数の制御点にて観察される第４信号との間の関係について、直達項を持たない状態空間モデルで表した場合に、
前記第１信号、前記第２信号、前記第３信号及び前記第４信号に基づいて、前記制御対象及び前記目標システムを同定する同定部と、
前記第１信号及び前記第２信号に基づいて、前記制御対象から遅延時間を分離し、当該遅延時間を分離した分離制御対象を求めると共に、前記第３信号及び前記第４信号に基づいて、前記目標システムから遅延時間を分離し、当該遅延時間を分離した分離目標システムを求める遅延時間分離部と、
入力信号を前記第３信号とし、出力信号を、前記仮制御器のゲイン、及び前記分離制御対象と前記分離目標システムとの間の出力追従誤差とし、前記ゲインと前記出力追従誤差との間のバランスを定めるパラメータをバランス係数として、前記入力信号を入力する前記仮制御器、前記仮制御器により出力された信号を入力する前記分離制御対象、前記仮制御器により出力された信号に前記バランス係数を乗算し、乗算結果を前記ゲインとして出力する演算器、前記入力信号を入力する前記分離目標システム、及び、前記分離制御対象により出力された信号から前記分離目標システムにより出力された信号を減算し、減算結果を前記出力追従誤差として出力する減算器からなる全体システムを構成するシステム構成部と、
前記システム構成部により構成された前記全体システムのＨ∞ノルムが最小となるように、前記仮制御器のパラメータを決定するパラメータ決定部と、
前記遅延時間分離部により分離された前記遅延時間に基づいて、前記第３信号を入力する前記第１遅延器の遅延時間を算出すると共に、前記第１信号を出力する前記第２遅延器の遅延時間を算出する遅延時間算出部と、を備え、
前記全体システムを、直達項を持つ状態空間モデルで表し、前記全体システムの伝達関数をＨall（ｚ）、前記全体システムの係数行列をＡall、入力側の係数行列をＢall、出力側の係数行列をＣall及び直達項の係数行列をＤallとし、
前記仮制御器を、直達項を持つ状態空間モデルで表し、前記仮制御器の係数行列をＡ0、入力側の係数行列をＢ0、出力側の係数行列をＣ0及び直達項の係数行列をＤ0とし、
遅延時間を分離した前記制御対象である遅延分離制御対象を、直達項を持たない状態空間モデルで表し、前記遅延分離制御対象の係数行列をＡｗ＃、入力側の係数行列をＢｐ＃、出力側の係数行列をＣｗ＃とし、
遅延時間を分離した前記目標システムである遅延分離目標システムを、直達項を持たない状態空間モデルで表し、前記遅延分離目標システムの係数行列を前記Ａｗ＃、入力側の係数行列をＢｔ＃、出力側の係数行列を前記Ｃｗ＃とし、
前記バランス係数をαとして、以下の数式：

で表される場合に、
前記パラメータ決定部は、
前記全体システムＨallのＨ∞ノルムを最小化する条件の数式：

を線形行列不等式で表したときに、パラメータγが最小となるように、前記仮制御器の係数行列Ａ0、前記入力側の係数行列Ｂ0、前記出力側の係数行列Ｃ0及び前記直達項の係数行列Ｄ0を、前記パラメータとして決定する、ことを特徴とする制御器設計装置。 An inverse system for a controlled object including a predetermined number of control points where a signal is observed and a predetermined number of secondary sound sources that present signals to the predetermined number of control points, respectively, to the predetermined number of secondary sound sources. Controller design for acoustic signals to design a controller that combines the characteristics of the reverse system that outputs a signal and the characteristics of a target system that includes a predetermined number of virtual sound sources that present signals to the predetermined number of control points. In the device
The controller comprises a first delayer, a temporary controller and a second delayer.
The control target is a state space model having no direct term regarding the relationship between the first signal presented from the predetermined number of secondary sound sources and the second signal observed at the predetermined number of control points. Represent,
The target system is represented by a state-space model having no direct term regarding the relationship between the third signal presented by the predetermined number of virtual sound sources and the fourth signal observed at the predetermined number of control points. If you do
An identification unit that identifies the controlled object and the target system based on the first signal, the second signal, the third signal, and the fourth signal.
Based on the first signal and the second signal, the delay time is separated from the control target, the separation control target separated from the delay time is obtained, and the third signal and the fourth signal are used as the basis. A delay time separator that separates the delay time from the target system and obtains a separation target system that separates the delay time.
The input signal is the third signal, the output signal is the gain of the temporary controller, and the output tracking error between the separation control target and the separation target system, and the gain is between the output tracking error. With the parameter that determines the balance as the balance coefficient, the temporary controller that inputs the input signal, the separation control target that inputs the signal output by the temporary controller, and the balance coefficient for the signal output by the temporary controller. Is multiplied, and the arithmetic unit that outputs the multiplication result as the gain, the separation target system that inputs the input signal, and the signal output by the separation target system are subtracted from the signal output by the separation control target. , A system configuration unit that constitutes an entire system consisting of a subtractor that outputs the subtraction result as the output tracking error.
A parameter determination unit that determines the parameters of the temporary controller so that the H∞ norm of the entire system configured by the system configuration unit is minimized.
Based on the delay time separated by the delay time separation unit, the delay time of the first delay device for inputting the third signal is calculated, and the delay of the second delay device for outputting the first signal is calculated. It is equipped with a delay time calculation unit that calculates the time.
The whole system is represented by a state-space model having a direct term, the transfer function of the whole system is Hall (z), the coefficient matrix of the whole system is Aall, the coefficient matrix of the input side is Ball, and the coefficient matrix of the output side is. Let Dall be the coefficient matrix of the Call and the direct term, and let it be Dall.
The temporary controller is represented by a state-space model having a direct term, the coefficient matrix of the temporary controller is A0, the coefficient matrix of the input side is B0, the coefficient matrix of the output side is C0, and the coefficient matrix of the direct term is D0. ,
The delay separation control target, which is the control target with the delay time separated, is represented by a state-space model that does not have a direct term, the coefficient matrix of the delay separation control target is Aw #, the coefficient matrix on the input side is Bp #, and the output side. Let Cw # be the coefficient matrix of
The delay separation target system, which is the target system in which the delay time is separated, is represented by a state space model having no direct term, the coefficient matrix of the delay separation target system is Aw #, the coefficient matrix on the input side is Bt #, and the output. Let the coefficient matrix on the side be Cw #, and let it be.
With the balance coefficient as α, the following formula:

When represented by
The parameter determination unit is
The formula of the condition for minimizing the H∞ norm of the whole system Hall:

Is expressed by a linear matrix inequality, the coefficient matrix A0 of the temporary controller, the coefficient matrix B0 on the input side, the coefficient matrix C0 on the output side, and the coefficient matrix of the direct term so that the parameter γ is minimized. A controller design device characterized in that D0 is determined as the parameter. コンピュータを、請求項１に記載の制御器設計装置として機能させるプログラム。 A program that causes a computer to function as the controller design device according to claim 1.

Images