TOP > 外国特許検索 > System estimation method and program, recording medium, and system estimation device

System estimation method and program, recording medium, and system estimation device 新技術説明会 実績あり

外国特許コード F110005755
整理番号 Y0412WO
掲載日 2011年9月14日
出願国 アメリカ合衆国
出願番号 56751404
公報番号 20070185693
公報番号 7480595
出願日 平成16年8月5日(2004.8.5)
公報発行日 平成19年8月9日(2007.8.9)
公報発行日 平成21年1月20日(2009.1.20)
国際出願番号 JP2004011568
国際公開番号 WO2005015737
国際出願日 平成16年8月5日(2004.8.5)
国際公開日 平成17年2月17日(2005.2.17)
優先権データ
  • 特願2003-291614 (2003.8.11) JP
  • 2004WO-JP11568 (2004.8.5) WO
発明の名称 (英語) System estimation method and program, recording medium, and system estimation device 新技術説明会 実績あり
発明の概要(英語) (US7480595)
It is possible to establish an estimation method capable of logically and optimally deciding a forgetting coefficient and develop an estimation algorithm and a high-speed algorithm which are numerically stable.
Firstly, a Processing section reads out or receives an upper limit value gammaf from a storage section or an input section (S101).
The.processing section decides a forgetting coefficient rho by equation (15) (S103).
After this, according to the forgetting coefficient rho, the processing section executes a hyper H∞ filter of equations (10-13) (S105).
The processing section (101) calculates the existence condition of equation (17) (or equation (18) which will be given later) (S107).
When the existence condition is satisfied at all the times (S109), gammaf is decreased by Deltagamma and the same processing is repeated (S111).
On the other hand, when the existence condition is not satisfied by a certain gammaf (S109), the Deltagamma is added to the gammaf and the sum is output to an output section and/or stored in the storage section as an optimal value gammafOP of the gammaf (S113).
特許請求の範囲(英語) [claim1]
1. A system estimation method, for a communication system or a sound system or sound field reproduction or noise control, for making state estimation robust and optimizing a forgetting factor rho simultaneously in an estimation algorithm, in which for a state space model expressed by following expressions:
xk+1=Fkxk+Gkwk
yk=Hkxk+vk
zk=Hkxk
here, xk: a state vector or simply a state,wk: a system noise,vk: an observation noise,yk: an observation signal,zk: an output signal,Fk: dynamics of a system, andGk: a drive matrix,as an evaluation criterion, a maximum value of an energy gain which indicates a ratio of a filter error to a disturbance including the system noise wk and the observation noise vk and is weighted with the forgetting factor rho is suppressed to be smaller than a term corresponding to a previously given upper limit value gamma f, andthe system estimation method comprises:a step at which a processing section inputs the upper limit value gamma f, the observation signal yk as an input of a filter and a value including an observation matrix Hk from a storage section or an input section;
a step at which the processing section determines the forgetting factor rho ;
as a following function of gamma f;

rho =1-chi (gamma f)
where chi (gamma f) denotes a monotonicaly damping function of gamma f to satisfy chi (1)=1 and chi (infin )=0,;
a step of executing a hyper Hinfin filter at which the processing section reads out an initial value or a value including the observation matrix Hk at a time from the storage section and obtains a filter gain Ks,k by using the forgetting factor rho and a gain matrix Kk and by following expressions (20) to (22):
(Equation image 34 not included in text) THETA (k) denotes a J-unitary matrix, that is, satisfiesTHETA (k) JTHETA H(k)T=J, J=(J1+I), I denotes a unit matrix, Kk(:, 1) denotes a column vector of a first column of the matrix Kk
(Equation image 35 not included in text)
here, x^k|k: the estimated value of the state xk at the time k using the observation signals y0 to yk,yk: the observation signal,Fk: the dynamics of the system, Fk=I for simplification,Ks,k: the filter gain,Hk: the observation matrix,SIGMA ^k|k: corresponding to a covariance matrix of an error of x^k|k,THETA (k): the J-unitary matrix, andRe,k: an auxiliary variable,a step at which the processing section stores an estimated value of the state xk by the hyper Hinfin , filter into the storage section;a step at which the processing section calculates an existence condition based on the upper limit value gamma f and the forgetting factor rho by the obtained observation matrix Hi or the observation matrix Hi and the filter gain Ks,i, anda step at which the processing section decreases the upper limit value gamma f by a factor of DELTA gamma and stores the resultant value into the storage section while the existence condition is satisfied in the step of executing the hyper Hinfin , filter,wherein the Hinfin filter equation is applied to obtain the state estimated value x^k|K=[h^1[k],. . . , h^N[k]]T, where h^[k] is the estimated value of impulse response,a pseudo-echo is estimated by a following expression:
(Equation image 36 not included in text)
and an actual echo is cancelled by the obtained pseudo-echo.
[claim2]
2. The system estimation method according to claim 1, wherein the processing section calculates the existence condition in accordance with a following expression:
(Equation image 37 not included in text)
[claim3]
3.
The system estimation method according to claim 1, wherein the processing section calculates the existence condition in accordance with a following expression:
(Equation image 38 not included in text)
where the forgetting factor rho and the upper limit value gamma f have a following relation: 0<rho =1-chi (gamma f) <= 1, where chi (gamma f) denotes a monotonically damping function of gamma f to satisfy chi (1)=1 and chi (infin )=0.
[claim4]
4. The system estimation method according to claim 1, wherein the step of executing the hyper Hinfin filter includes: a step at which the processing section calculates SIGMA ^k+1|k1/2 by using the expression (22);a step at which the processing section calculates the filter gain Ks,k based on an initial condition of SIGMA ^k|k-1 and an initial condition of Ck, by using the expression (21);a step at which the processing section updates a filter equation of the Hinfin filter of the expression (20);
and
a step at which the processing section repeatedly executes the step of calculating by using the expression (20), the step of calculating by using the expression (21) and, the step of updating while advancing the time k.
[claim5]
5. A system estimation method, for a communication system or a sound system or sound field reproduction or noise control, for making state estimation robust and optimizing a forgetting factor rho simultaneously in an estimation algorithm, in which for a state space model expressed by following expressions:
xk+1=Fkxk+Gkwk
yk=Hkxk+vk
zk=Hkxk
here, xk: a state vector or simply a state,wk: a system noise,vk: an observation noise,yk: an observation signal,zk: an output signal,Fk: dynamics of a system, andGk: a drive matrix,as an evaluation criterion, a maximum value of an energy gain which indicates a ratio of a filter error to a disturbance including the system noise wk and the observation noise vk and is weighted with the forgetting factor rho is suppressed to be smaller than a term corresponding to a previously given upper limit value gamma f, andthe system estimation method comprises:a step at which a processing section inputs the upper limit value gamma f, the observation signal yk as an input of a filter and a value including an observation matrix Hk from a storage section or an input section;
a step at which the processing section determines the forgetting factor rho relevant to the state space model in accordance with the upper limit value gamma f;
as a following function of gamma f,
rho =1-chi (gamma f)
where chi (gamma f) denotes a monotonicaly damping function of gamma f to satisfy chi (1)=1 and chi (infin )=0, a step of executing a hyper Hinfin filter at which the processing section reads out an initial value or a value including the observation matrix Hk at a time from the storage section and obtains a filter gain Ks,k by using the forgetting factor rho and a gain matrix Kk and by following expressions:
(Equation image 39 not included in text)
here, THETA (k) denotes an arbitrary J-unitary matrix, and {hacek over (C)}k={hacek over (C)}k+1PSI is established, where
(Equation image 40 not included in text)
(Equation image 41 not included in text)
here, x^k|k: the estimated value of the state xk at the time k using the observation signals y0 to yk,yk: the observation signal,Ks,k: the filter gain,Hk: the observation matrix,THETA (k): the J-unitary matrix, andRe,k: an auxiliary variable,a step at which the processing section stores an estimated value of the state xk by the hyper Hinfin filter into the storage section;a step at which the processing section calculates an existence condition based on the upper limit value gamma f and the forgetting factor rho by the obtained observation matrix Hi or the observation matrix Hi and the filter gain Ks,i, anda step at which the processing section decreases the upper limit value gamma f by a factor of DELTA gamma and stores the resultant value into the storage section while the existence condition is satisfied in the step of executing the hyper Hinfin filter,wherein the Hinfin filter equation is applied to obtain the state estimated value x^k|K=[h^1[k],. . . , h^N[k]]T, where h^[k] is the estimated value of impulse response,a pseudo-echo is estimated by a following expression:
(Equation image 42 not included in text)
and an actual echo is cancelled by the obtained pseudo-echo.
[claim6]
6. The system estimation method according to claim 5, wherein the step of executing the hyper Hinfin filter includes: a step at which the processing section calculates K-k based on an initial condition of Re,k+1, Rr,k+1 and L{tilde over ( )}k+1 by using the expression (63);a step at which the processing section calculates the filter gain Ks,k based on the initial condition and by using the expression (62);a step at which the processing section updates a filter equation of the Hinfin filter of the expression (61);
and
a step at which the processing section repeatedly executes the step of calculating by using the expression (63), the step of calculating by using the expression (62), and, the step of updating while advancing the time k.
[claim7]
7. A system estimation method, for a communication system or a sound system or sound field reproduction or noise control, for making state estimation robust and optimizing a forgetting factor rho simultaneously in an estimation algorithm, in which for a state space model expressed by following expressions:
xk+1=Fkxk+Gkwk
yk=Hkxk+vk
zk=Hkxk
here, xk: a state vector or simply a state,wk: a system noise,vk: an observation noise,yk: an observation signal,zk: an output signal,Fk: dynamics of a system, andGk: a drive matrix,as an evaluation criterion, a maximum value of an energy gain which indicates a ratio of a filter error to a disturbance including the system noise wk and the observation noise vk and is weighted with the forgetting factor rho is suppressed to be smaller than a term corresponding to a previously given upper limit value gamma f, andthe system estimation method comprises:a step at which a processing section inputs the upper limit value gamma f, the observation signal yk as an input of a filter and a value including an observation matrix Hk from a storage section or an input section;a step at which the processing section determines the forgetting factor rho as a following function gamma f,
rho =1-chi (gamma f)
where chi (gamma f) denotes a monotonicaly damping function of gamma f to satisfy chi (1)=1 and chi (infin )=0, a step of executing a hyper Hinfin , filter at which the processing section reads out an initial value or a value including the observation matrix Hk at a time from the storage section and obtains a filter gain Ks,k by using the forgetting factor rho and a gain matrix K-k and by following expressions:
(Equation image 43 not included in text)
here, yk: the observation signal,Fk: the dynamics of the system, Fk=I for simplification,Hk: the observation matrix,x^k|k: the estimated value of the state xk at the time k using the observation signals y0 to yk, Ks,k: the filter gain, obtained from the gain matrix K-k, andRe,k,L{tilde over ( )}k: an auxiliary variable,a step at which the processing section stores an estimated value of the state xk by the hyper Hinfin filter into the storage section;a step at which the processing section calculates an existence condition based on the upper limit value gamma f and the forgetting factor rho by the obtained observation matrix Hi or the observation matrix Hi and the filter gain Ks,i, anda step at which the processing section decreases the upper limit value gamma f by a factor of DELTA gamma and stores the resultant value into the storage section while the existence condition is satisfied in the step of executing the hyper Hinfin filter,wherein the Hinfin filter equation is applied to obtain the state estimated value x^k|K=[h^1[k],. . . , h^N[k]]T, where h^[k] is the estimated value of impulse response,a pseudo-echo is estimated by a following expression:
(Equation image 44 not included in text)
and an actual echo is cancelled by the obtained pseudo-echo.
[claim8]
8. The system estimation method according to claim 1, wherein an estimated value zvk|k of the output signal is obtained from the state estimated value x^k|k at the time k by a following expression:
zvk|k=Hkx^k|k.
[claim9]
9. A system estimation program product, for a communication system or a sound system or sound field reproduction or noise control, embodied on a computer-readable medium and comprising code that, when executed, causes a computer to make state estimation robust and to optimize a forgetting factor rho simultaneously in an estimation algorithm, in which for a state space model expressed by following expressions:
xk+1=Fkxk+Gkwk
yk=Hkxk+vk
zk=Hkxk
here, xk: a state vector or simply a state,wk: a system noise,vk: an observation noise,yk: an observation signal,zk: an output signal,Fk: dynamics of a system, andGk: a drive matrix,as an evaluation criterion, a maximum value of an energy gain which indicates a ratio of a filter error to a disturbance including the system noise wk and the observation noise vk and is weighted with the forgetting factor rho is suppressed to be smaller than a term corresponding to a previously given upper limit value gamma f, andthe system estimation program causes the computer to execute:a step at which a processing section inputs the upper limit value gamma f, the observation signal yk as an input of a filter and a value including an observation matrix Hk from a storage section or an input section;a step at which the processing section determines the forgetting factor rho as a following function gamma f,
rho =1-chi (gamma f)
where rho (gamma f) denotes a monotonicaly damping function of gamma f to satisfy chi (1)=1 and chi (infin )=0, a step of executing a hyper Hinfin filter at which the processing section reads out an initial value or a value including the observation matrix Hk at a time from the storage section and obtains a filter gain Ks,k by using the forgetting factor rho and a gain matrix K-k and by following expressions:
(Equation image 45 not included in text)
here, yk: the observation signal,Fk: the dynamics of the system, Fk=I for simplification,Hk: the observation matrix,x^k|k: the estimated value of the state xk at the time k using the observation signals y0 to yk, Ks,k: the filter gain, obtained from the gain matrix K-k, andRe,k, L{tilde over ( )}k: an auxiliary variable,a step at which the processing section stores an estimated value of the state xk by the hyper Hinfin filter into the storage section;a step at which the processing section calculates an existence condition based on the upper limit value gamma f and the forgetting factor rho by the obtained observation matrix Hi or the observation matrix Hi and the filter gain Ks,i, anda step at which the processing section decreases the upper limit value gamma f by a factor of DELTA gamma and stores the resultant value into the storage section while the existence condition is satisfied in the step of executing the hyper Hinfin filter;wherein the Hinfin filter equation is applied to obtain the state estimated value x^k|K=[h^1[k],. . . , h^N[k]]T, where h^[k] is the estimated value of impulse response,a pseudo-echo is estimated by a following expression:
(Equation image 46 not included in text)
and an actual echo is cancelled by the obtained pseudo-echo.
[claim10]
10. A computer readable recording medium recording a system estimation program product, for a communication system or a sound system or sound field reproduction or noise control, embodied on a computer-readable medium and comprising code that, when executed, causes a computer to make state estimation robust and to optimize a forgetting factor rho simultaneously in an estimation algorithm, in which for a state space model expressed by following expressions:
xk+1=Fkxk+Gkwk
yk=Hkxk+vk
zk=Hkxk
here, xk: a state vector or simply a state,wk: a system noise,vk: an observation noise,yk: an observation signal,zk: an output signal,Fk: dynamics of a system, andGk: a drive matrix,as an evaluation criterion, a maximum value of an energy gain which indicates a ratio of a filter error to a disturbance including the system noise wk and the observation noise vk and is weighted with the forgetting factor rho is suppressed to be smaller than a term corresponding to a previously given upper limit value gamma f, andthe computer readable recording medium recording the system estimation program causes the computer to execute:a step at which a processing section inputs the upper limit value gamma f, the observation signal yk as an input of a filter and a value including an observation matrix Hk from a storage section or an input section;a step at which the processing section determines the forgetting factor rho relevant to the state space model in accordance as a following function gamma f,
rho =1-chi (gamma f)
where chi (gamma f) denotes a monotonicaly damping function of gamma f to satisfy chi (1)=1 and chi (infin )=0, a step of executing a hyper Hinfin filter at which the processing section reads out an initial value or a value including the observation matrix Hk at a time from the storage section and obtains a filter gain Ks,k by using the forgetting factor rho and a gain matrix K-k and by following expressions:
(Equation image 47 not included in text)
here, yk: the observation signal,Fk: the dynamics of the system, Fk=I for simplification,Hk: the observation matrix,x^k|k: the estimated value of the state xk at the time k using the observation signals y0 to yk,Ks,k: the filter gain, obtained from the gain matrix K-k, andRe,k, L{tilde over ( )}k: an auxiliary variable,a step at which the processing section stores an estimated value of the state xk by the hyper Hinfin filter into the storage section;a step at which the processing section calculates an existence condition based on the upper limit value gamma f and the forgetting factor rho by the obtained observation matrix Hi or the observation matrix Hi and the filter gain Ks,i, anda step at which the processing section decreases the upper limit value gamma f by a factor of DELTA gamma and stores the resultant value into the storage section while the existence condition is satisfied in the step of executing the hyper Hinfin filter,wherein the Hinfin filter equation is applied to obtain the state estimated value x^k|K=[h^1[k],. . . , h^N[k]]T, where h^[k] is the estimated value of impulse response,a pseudo-echo is estimated by a following expression:
(Equation image 48 not included in text)
and an actual echo is cancelled by the obtained pseudo-echo.
[claim11]
11. A system estimation device, for communication system or a sound system or sound field reproduction or noise control, for making state estimation robust and optimizing a forgetting factor rho simultaneously in an estimation algorithm, in which for a state space model expressed by following expressions:
xk+1=Fkxk+Gkwk
yk=Hkxk+vk
zk=Hkxk
here, xk: a state vector or simply a state,wk: a system noise,vk: an observation noise,yk: an observation signal,zk: an output signal,Fk: dynamics of a system, andGk: a drive matrix,as an evaluation criterion, a maximum value of an energy gain which indicates a ratio of a filter error to a disturbance including the system noise wk and the observation noise vk and is weighted with the forgetting factor rho is suppressed to be smaller than a term corresponding to a previously given upper limit value gamma f, andthe system estimation device comprises:a processing section to execute the estimation algorithm;
and
a storage section to which reading and/or writing is performed by the processing section and which stores respective observed values, set values, and estimated values relevant to the state space model,further comprising:a means at which the processing section inputs the upper limit value gamma f, the observation signal yk as an input of a filter and a value including an observation matrix Hk from the storage section or an input section;a means at which the processing section determines the forgetting factor rho as a following function gamma f,
rho =1-chi (gamma f)
where chi (gamma f) denotes a monotonicaly damping function of gamma f to satisfy chi (1)=1 and chi (infin )=0, a means of executing a hyper Hinfin filter at which the processing section reads out an initial value or a value including the observation matrix Hk at a time from the storage section and obtains a filter gain Ks,k by using the forgetting factor rho and a gain matrix K-k and by following expressions:
(Equation image 49 not included in text)
here, yk: the observation signal,Fk: the dynamics of the system, Fk=I for simplification,Hk: the observation matrix,x^k|k: the estimated value of the state xk at the time k using the observation signals y0 to yk, Ks,k: the filter gain, obtained from the gain matrix K-k, andRe,k, L{tilde over ( )}k: an auxiliary variable,a means at which the processing section stores an estimated value of the state xk by the hyper Hinfin filter into the storage section;a means at which the processing section calculates an existence condition based on the upper limit value gamma f and the forgetting factor rho by the obtained observation matrix Hi or the observation matrix Hi and the filter gain Ks,i anda means at which the processing section decreases the upper limit value gamma f by a factor of DELTA gamma and stores the resultant value into the storage section while the existence condition is satisfied in the means of executing the hyper Hinfin filter,wherein the Hinfin filter equation is applied to obtain the state estimated value x^k|K=[h^1[k],. . . , h^N[k]]T, where h^[k] is the estimated value of impulse response,a pseudo-echo is estimated by a following expression:
(Equation image 50 not included in text)
and an actual echo is cancelled by the obtained pseudo-echo.
[claim12]
12. The system estimation method according to claim 5, wherein the processing section calculates the existence condition in accordance with a following expression:
(Equation image 51 not included in text)
where the forgetting factor rho and the upper limit value gamma f have a following relation: 0<rho =1-chi (gamma f) <= 1, where chi (gamma f) denotes a monotonically damping function of gamma f to satisfy chi (1)=1 and chi (infin )=0.
[claim13]
13. The system estimation method according to claim 5, wherein an estimated value zvk|k of the output signal is obtained from the state estimated value x^k|k at the time k by a following expression:
zvk|k=Hkx^k|K.
[claim14]
14. The system estimation method according to claim 7, wherein the processing section calculates the existence condition in accordance with a following expression:
(Equation image 52 not included in text)
where the forgetting factor rho and the upper limit value gamma f have a following relation: 0<rho =1-chi (gamma f) <= 1, where chi (gamma f) denotes a monotonically damping function of gamma f to satisfy chi (1)=1 and chi (infin )=0.
[claim15]
15. The system estimation method according to claim 7, wherein an estimated value zvk|k of the output signal is obtained from the state estimated value x^k|k at the time k by a following expression:
zvk|k=Hkx^k|k.
  • 発明者/出願人(英語)
  • NISHIYAMA KIYOSHI
  • JAPAN SCIENCE AND TECHNOLOGY AGENCY
国際特許分類(IPC)
ライセンスをご希望の方、特許の内容に興味を持たれた方は、問合せボタンを押してください。

PAGE TOP

close
close
close
close
close
close