"人間の上体に取り付けられる第1装具と、前記人間の左右の大腿のそれぞれに対して取り付けられる一対の第2装具と、一対のアクチュエータと、前記人間の左右の股関節角度のそれぞれに応じた信号を出力するように構成されている左股関節角度センサ及び右股関節角度センサと、少なくとも前記左股関節角度センサ及び前記右股関節角度センサの出力信号に基づいて前記一対のアクチュエータのそれぞれの動作を制御するように構成されている制御装置とを備え、前記一対のアクチュエータのそれぞれの動作によって前記第1装具に対して前記一対の第2装具のそれぞれを動かすことにより、前記上体に対する前記左右の大腿の相対的な周期的運動を伴う前記人間の歩行運動を補助する装置であって、前記制御装置は、左屈曲係数、左伸展係数、右屈曲係数及び右伸展係数のそれぞれの値の大小に応じて、前記人間の左大腿の屈曲運動及び伸展運動並びに右大腿の屈曲運動及び伸展運動のそれぞれを補助するための前記一対のアクチュエータのそれぞれの動作振幅の大小を制御するように構成され、前記左股関節角度センサの出力の時間変化態様を表わす波形信号である左運動振動子と、前記右股関節角度センサの出力の時間変化態様を表わす波形信号である右運動振動子との非対称度を評価するように構成されている非対称度評価要素と、前記非対称度評価要素により評価された前記非対称度を低下させるように、前記左屈曲係数、前記左伸展係数、前記右屈曲係数及び前記右伸展係数のうち少なくとも1つの値を調節するように構成されている調節要素とを備えていることを特徴とする歩行補助装置。請求項1記載の歩行補助装置において、前記調節要素は、前記左屈曲係数及び前記右伸展係数が同値になり、かつ、前記右屈曲係数及び前記左伸展係数が同値になるように、前記左屈曲係数、前記左伸展係数、前記右屈曲係数及び前記右伸展係数のうち少なくとも1つの値を調節するように構成されていることを特徴とする歩行補助装置。請求項1又は2記載の歩行補助装置において、前記調節要素は、前記左屈曲係数、前記左伸展係数、前記右屈曲係数及び前記右伸展係数のうち少なくとも1つの値を、前記非対称度を変数とする増加関数によって定義される補正量に応じて調節するように構成されていることを特徴とする歩行補助装置。請求項1~3のうちいずれか1つに記載の歩行補助装置において、前記制御装置は、左時定数及び右時定数のそれぞれの値の大小に応じて、前記人間の左大腿及び右大腿のそれぞれの運動を補助するための前記アクチュエータの動作周波数の高低を制御するように構成され、前記調節要素は、前記左股関節角度センサ及び前記右股関節角度センサのそれぞれの出力から得られる、前記人間の左右の股関節角度の時間変化態様を表わす波形信号に基づき、前記左時定数及び前記右時定数のそれぞれの値を調節するように構成されていることを特徴とする歩行補助装置。請求項4記載の歩行補助装置において、前記左股関節角度センサ及び前記右股関節角度センサのそれぞれの出力信号に基づき、前記人間の左右の股関節角度の差分の指定期間にわたるサンプリングにより得られる波形信号である差分振動子に窓を掛ける窓処理を実行するように構成されている窓処理要素と、前記窓が掛けられた前記差分振動子を周波数解析することによりパワースペクトルを取得するように構成されている周波数解析処理要素と、前記パワースペクトルにおいて閾値以上の高さを有するとともに、最低の周波数帯域に位置するピークを示す基本周波数を決定するように構成されているスペクトル解析処理要素とを備え、前記調節要素は、前記基本周波数の逆数に比例するように、前記左時定数及び前記右時定数のそれぞれの値を調節するように構成されていることを特徴とする歩行補助装置。請求項1~5記載の歩行補助装置において、前記左運動振動子及び前記右運動振動子のそれぞれにおける、前記人間の上体に対する大腿の屈曲運動振幅を表わす屈曲振幅と、前記人間の大腿の伸展運動振幅を表わす伸展振幅と、前記屈曲振幅及び前記伸展振幅の和である総合振幅とのうちいずれか1つの偏差又は当該偏差の複数周期にわたる平均値を前記非対称度として評価するように構成されていることを特徴とする歩行補助装置。請求項1~6記載の歩行補助装置において、前記制御装置は、人間の各大腿の屈曲運動状態及び伸展運動状態のそれぞれを表わす複数の状態変数の連立微分方程式により定義され、入力波形信号に基づいて定まる前記連立微分方程式の解の時間変化態様により表現され、出力波形信号を生成する状態振動子モデルに対して、前記左運動振動子及び前記右運動振動子を当該入力波形信号として入力することにより、前記左屈曲係数、前記左伸展係数、前記右屈曲係数及び前記右伸展係数のそれぞれの値に応じた振幅にしたがって変化する左屈曲振動子、左伸展振動子、右屈曲振動子及び右伸展振動子を当該出力波形信号として生成するように構成されている状態振動子生成要素と、前記左屈曲振動子及び前記左伸展振動子を組み合わせることにより左側の前記アクチュエータに対する制御指令信号としての左制御振動子を生成するとともに、前記右屈曲振動子及び前記右伸展振動子を組み合わせることにより右側の前記アクチュエータに対する制御指令信号としての右制御振動子を生成するように構成されている制御振動子生成要素とを備えていることを特徴とする歩行補助装置。しかし、人間の左右片側の身体機能の低下等のために、当該人間の左右の運動形態が不均等になる状況において、歩幅等をその目標値に一律に一致させようとすると、その不均等が是正されないだけではなく、かえって当該不均等を助長する可能性がある。そこで、本発明は、人間の左右の運動形態の均等を図りながら、当該人間の歩行運動を補助することができる装置を提供することを解決課題とする。本発明は、人間の上体に取り付けられる第1装具と、前記人間の左右の大腿のそれぞれに対して取り付けられる一対の第2装具と、一対のアクチュエータと、前記人間の左右の股関節角度のそれぞれに応じた信号を出力するように構成されている左股関節角度センサ及び右股関節角度センサと、少なくとも前記左股関節角度センサ及び前記右股関節角度センサの出力信号に基づいて前記一対のアクチュエータのそれぞれの動作を制御するように構成されている制御装置とを備え、前記一対のアクチュエータのそれぞれの動作によって前記第1装具に対して前記一対の第2装具のそれぞれを動かすことにより、前記上体に対する前記左右の大腿の相対的な周期的運動を伴う前記人間の歩行運動を補助する装置に関する。当該歩行補助装置において、前記制御装置は、左屈曲係数、左伸展係数、右屈曲係数及び右伸展係数のそれぞれの値の大小に応じて、前記人間の左大腿の屈曲運動及び伸展運動並びに右大腿の屈曲運動及び伸展運動のそれぞれを補助するための前記アクチュエータの動作振幅の大小を制御するように構成されている。前記課題を解決するための本発明の本発明の歩行補助装置は、前記左股関節角度センサの出力の時間変化態様を表わす波形信号である左運動振動子と、前記右股関節角度センサの出力の時間変化態様を表わす波形信号である右運動振動子との非対称度を評価するように構成されている非対称度評価要素と、前記非対称度評価要素により評価された前記非対称度を低下させるように、前記左屈曲係数、前記左伸展係数、前記右屈曲係数及び前記右伸展係数のうち少なくとも1つの値を調節するように構成されている調節要素とを備えていることを特徴とする。本発明の歩行補助装置によれば、人間の左右の運動形態の不均等の度合が、左右の運動振動子の非対称度に現われることに鑑みて、当該非対称度の低減が図られる。その結果、人間の左右の運動形態の均等を図りながら、当該人間の歩行運動が補助されうる。前記調節要素は、前記左屈曲係数及び前記右伸展係数が同値になり、かつ、前記右屈曲係数及び前記左伸展係数が同値になるように、前記左屈曲係数、前記左伸展係数、前記右屈曲係数及び前記右伸展係数のうち少なくとも1つの値を調節するように構成されていることが好ましい。前記調節要素は、前記左屈曲係数、前記左伸展係数、前記右屈曲係数及び前記右伸展係数のうち少なくとも1つの値を、前記非対称度を変数とする増加関数によって定義される補正量に応じて調節するように構成されていることが好ましい。前記制御装置は、左時定数及び右時定数のそれぞれの値の大小に応じて、前記人間の左大腿及び右大腿のそれぞれの運動を補助するための前記アクチュエータの動作周波数の高低を制御するように構成され、前記調節要素は、前記左股関節角度センサ及び前記右股関節角度センサのそれぞれの出力から得られる、前記人間の左右の股関節角度の時間変化態様を表わす波形信号に基づき、前記左時定数及び前記右時定数のそれぞれの値を調節するように構成されていることが好ましい。前記左股関節角度センサ及び前記右股関節角度センサのそれぞれの出力信号に基づき、前記人間の左右の股関節角度の差分の指定期間にわたるサンプリングにより得られる波形信号である差分振動子に窓を掛ける窓処理を実行するように構成されている窓処理要素と、前記窓が掛けられた前記差分振動子を周波数解析することによりパワースペクトルを取得するように構成されている周波数解析処理要素と、前記パワースペクトルにおいて閾値以上の高さを有するとともに、最低の周波数帯域に位置するピークを示す基本周波数を決定するように構成されているスペクトル解析処理要素とを備え、前記調節要素は、前記基本周波数の逆数に比例するように、前記左時定数及び前記右時定数のそれぞれの値を調節するように構成されていることが好ましい。前記左運動振動子及び前記右運動振動子のそれぞれにおける、前記人間の上体に対する大腿の屈曲運動振幅を表わす屈曲振幅と、前記人間の大腿の伸展運動振幅を表わす伸展振幅と、前記屈曲振幅及び前記伸展振幅の和である総合振幅とのうちいずれか1つの偏差又は当該偏差の複数周期にわたる平均値を前記非対称度として評価するように構成されていることが好ましい。前記制御装置は、人間の各大腿の屈曲運動状態及び伸展運動状態のそれぞれを表わす複数の状態変数の連立微分方程式により定義され、入力波形信号に基づいて定まる前記連立微分方程式の解の時間変化態様により表現され、出力波形信号を生成する状態振動子モデルに対して、前記左運動振動子及び前記右運動振動子を当該入力波形信号として入力することにより、前記左屈曲係数、前記左伸展係数、前記右屈曲係数及び前記右伸展係数のそれぞれの値に応じた振幅にしたがって変化する左屈曲振動子、左伸展振動子、右屈曲振動子及び右伸展振動子を当該出力波形信号として生成するように構成されている状態振動子生成要素と、前記左屈曲振動子及び前記左伸展振動子を組み合わせることにより左側の前記アクチュエータに対する制御指令信号としての左制御振動子を生成するとともに、前記右屈曲振動子及び前記右伸展振動子を組み合わせることにより右側の前記アクチュエータに対する制御指令信号としての右制御振動子を生成するように構成されている制御振動子生成要素とを備えていることが好ましい。歩行補助装置の構成以下、脚体等の左右を区別するために符号L及びRを用いるが、左右を区別する必要がない場合や左右成分を有するベクトルを表現する場合には当該符号を省略する。また、上体に対する各大腿の屈曲運動前方運動及び伸展運動後方運動を区別するために符号及び-を用いる。図1に示されている歩行補助装置10は、第1装具11と、左右一対の第2装具12と、左右一対のアクチュエータ14と、バッテリ16と、制御装置20と、股関節角度センサ202とを備えている。第1装具11は、人間又は利用者の上体又は腰部第1身体部分に巻き付けられるように装着される。第1装具11のうち少なくとも人間の背中に当接する背後部は、軽量合金、硬質樹脂又はカーボンファイバ等の剛性素材により構成され、その他の部分は繊維等の柔軟性素材により構成されている。第2装具12は、繊維等の柔軟性素材により構成され、人間の大腿第2身体部分に巻き付けられるように装着される。第2装具12は左右両側ではなく左右片側のみに設けられてもよい。アクチュエータ14は電動モータにより構成され、必要に応じてモータに加え減速機及びコンプライアンス機構のうち一方又は両方により構成される。アクチュエータ14は第1装具11が上体に取り付けられたとき、上体の左右両側に配置されるように第1装具11に連結されている。アクチュエータ14は、軽量の軽量合金、硬質樹脂又はカーボンファイバ等の剛性素材により形成されている連結部材15を介して大腿に装着される第2装具12に連結される。これにより、アクチュエータ14が動作することで、上体及び各大腿の相対運動が補助されるように上体及び各大腿に力が作用する。上体及び各大腿の相対運動には、離床している脚体の大腿の上体に対する前後運動が含まれ、着床している脚体に対する上体の前後運動が含まれる。バッテリ16は制御装置20とともに第1装具11の背後部に取り付けられるケース13に収納されており、アクチュエータ14及び制御装置20等に対して電力を供給する。なお、バッテリ16及び制御装置20のそれぞれ又はこれらを収納するケース13の配置箇所は適宜変更されてもよい。股関節角度センサ202は人間の腰部左右両側に配置されるロータリーエンコーダにより構成され、股関節角度に応じた信号を出力する。股関節角度は、大腿が基本前額面の前方にある場合には正値になり、大腿が基本前額面の後方にある場合には負値になるように定義される。制御装置20はコンピュータCPU、ROM、RAM、信号入力回路、信号出力回路等により構成される。と、このコンピュータのメモリ又は記憶装置に格納されているソフトウェアとにより構成されている。制御装置20は、バッテリ16からアクチュエータ14に対する供給電力を調節するほか、アクチュエータ14の動作を制御する。図2に示されているように、制御装置20は、後述の演算処理を実行する又は機能を発揮する運動振動子出力要素21と、非対称度評価要素22と、基本周波数推定要素23と、状態振動子生成要素24と、モデル調節要素242本発明の調節要素に該当する。と、制御振動子生成要素25とを備えている。基本周波数推定要素23は、差分振動子出力要素231と、窓処理要素232と、周波数解析処理要素233と、スペクトル解析処理要素234とを備えている。制御装置20の各構成要素の各構成要素は、記憶装置に格納されているプログラムを読み取り、当該プログラムにしたがって担当する演算処理を実行する演算処理装置により構成されている。各構成要素は、共通の演算処理装置により構成されていてもよく、物理的に別個の演算処理装置により構成されていてもよい。たとえば、基本周波数推定要素23が、その他の構成要素とは別個の演算処理装置により構成されていてもよい。運転スイッチ図示略が操作され、バッテリ16から制御装置20に対して電力が供給されることにより、制御装置20は諸機能を発揮しうる。歩行補助装置及び歩行状態推定装置の機能運動振動子出力要素21が、股関節角度センサ202の出力に基づき、人間の左右の股関節角度の時間変化を表わす波形信号を運動振動子φ=φL,φRとして生成する図3/STEP1。φLは左運動振動子を表わし、φRは右運動振動子を表わしている。続いて、非対称度評価要素22が、左運動振動子φL及び右運動振動子φRの非対称度sを評価する図3/STEP2。たとえば、図4に一点鎖線及び二点鎖線のそれぞれにより示されているように、左運動振動子φL及び右運動振動子φRのそれぞれが変化している場合について考察する。左運動振動子φLの極性が正であることは、上体に対して左大腿が屈曲していることを表わし、左運動振動子φLの極性が負であることは、上体に対して左大腿が伸展していることを表わしている。各周期における左運動振動子φLの最大値正値である屈曲振幅|φL|は、上体に対して左大腿が最も屈曲している時点での左股関節角度に相当する。左運動振動子φLの各周期における最小値負値の絶対値である伸展振幅|φL|は、上体に対して左大腿が最も伸展している状態での左股関節角度に相当する。同様に、右運動振動子φRの極性が正であることは、上体に対して右大腿が屈曲していることを表わし、右運動振動子φRの極性が負であることは、上体に対して右大腿が伸展していることを表わしている。右運動振動子φRの各周期における最大値正値である屈曲振幅|φR|は、上体に対して右大腿が最も屈曲している時点での右股関節角度に相当する。右運動振動子φRの各周期における最小値負値の絶対値である伸展振幅|φL|は、上体に対して右大腿が最も伸展している状態での右股関節角度に相当する。たとえば、左屈曲振幅|φL|と左伸展振幅|φL|との和である左総合振幅|φL|と、右屈曲振幅|φR|と右伸展振幅|φR|との和である右総合振幅|φR|との偏差|φL|-|φR|又は当該偏差の複数周期にわたる平均値が非対称度sとして評価される。図4に示されている状況では、右総合振幅|φR|が左総合振幅|φL|より大きいので、非対称度sは負値になる。また、基本周波数推定要素23が後述する演算処理を実行することにより、人間の歩行周期の逆数に相当する基本周波数f0を推定する図3/STEP3。さらに、モデル調節要素242が、非対称度s及び基本周波数f0正確には後述する補正後基本周波数f0dに基づき、状態振動子モデルを定義する係数の値を調節する図3/STEP4。状態振動子モデルは、各大腿の屈曲運動状態及び伸展運動状態のそれぞれを表わす複数の状態変数uii=L,L-,R,R-及び各大腿の屈曲運動状態及び伸展運動状態のそれぞれの順応性を表現するための自己抑制因子viの連立微分方程式010により定義されている。τ1L=cL-uLwLLξLwLRξR-λLvLf1f2KφLτ1L=cL-uLwLLξLwLRξR-λLvLf1f2KφLτ1R=cR-uRwRLξLwRRξR-λRvRf1f2KφRτ1R=cR-uRwRLξLwRRξR-λRvRf1f2KφRτ2i=-viξiξi=H=0。cLは、左大腿の屈曲運転状態を表わす状態変数uLの振幅の大小を定める左屈曲係数である。cLは、左大腿の伸展運転状態を表わす状態変数uLの振幅の大小を定める左伸展係数である。cRは、右大腿の屈曲運転状態を表わす状態変数uRの振幅の大小を定める右屈曲係数である。cRは、右大腿の伸展運転状態を表わす状態変数uRの振幅の大小を定める右伸展係数である。τ1Lは状態変数uL及びuLの周波数の高低を定める第1左時定数である。τ1Rは状態変数uR及びuRの周波数の高低を定める第1右時定数である。第1左時定数τ1L及び第1右時定数τ1Rは同値であってもことなる値であってもよい。τ2iは抑制因子viの周波数の高低を定める第2時定数である。第2時定数τ2iは同値であってもよく、異なる値であってもよい。たとえば、第2左屈曲時定数τ2Lと第2左伸展時定数τ2Lとが同値であり、第2右屈曲時定数τ2Rと第2右伸展時定数τ2Rとが同値である一方、第2左屈曲時定数τ2Lと第2右屈曲時定数τ2Rとが異なる値であってもよい。wijは状態変数ui及びujの相関関係を表わす負値の相関係数である。λLおよびλRは慣れ係数である。Kは運動振動子φに応じたフィードバック係数である。ωは人間の歩行周期の長短に応じて定まる固有角速度である。固有角速度ωは、たとえば基本周波数f0又は補正後基本周波数f0dに2πを乗じた値に設定される。固有角速度ωには左成分ωL及び右成分ωRが含まれるが、基本周波数f0には左右の別がないので、ωL及びωRは同値である。f1は正係数cを用いて式011により定義される固有角速度ωの1次関数である。f1=cω。f2は係数c0,c1およびc2を用いて式012により定義される固有角速度ωの2次関数である。f2=c0c1ωc2ω2。具体的には、モデル調節要素242は、非対称度sを変数とし、少なくとも非対称度sの標準的な定義域において増加関数fsにより定義される補正量に基づき、関係式021にしたがって各状態係数ciの値を設定する。cL=cR=χ1χ2fcL=cR=χ1χ2f。係数χ1>0は、歩行補助装置10の動作による人間の歩幅又は歩行率を定めるための係数である。係数χ1の値が大きく設定されることにより人間の歩幅が大きくなるように歩行運動が補助される一方、係数χ1の値が小さく設定されることにより人間の歩幅が小さくなるように歩行運動が補助される。係数χ2>0は、非対称度sを低下させるために左右の歩幅のそれぞれを別個に補正の刻み量又は補正速度を表わすゲイン係数である。係数χ2の値が大きく設定されることにより人間の歩幅が早く変化するように歩行運動が補助される一方、係数χ2の値が小さく設定されることにより人間の歩幅が遅く変化するように歩行運動が補助される。増加関数fsはたとえば関係式022により定義される。f=s1exp。関係式022により定義される関数fsは、図5aに示されているような変化特性を示す。このため、非対称度sが正値であってその絶対値が大きくなるにつれて、左屈曲係数cL及び右伸展係数cRの値が減少するように調節される一方、左伸展係数cL及び右屈曲係数cRの値が増加するように調節される関係式021参照。また、非対称度sが負値であってその絶対値が大きくなるにつれて、左屈曲係数cL及び右伸展係数cRの値が増加するように調節される一方、左伸展係数cL及び右屈曲係数cRの値が減少するように調節される。図4に示されている状況では前記のように非対称度sは負値なので、図5bに示されているように、左屈曲係数cL及び右伸展係数cRの値が、左伸展係数cL及び右屈曲係数cRの値よりも相対的に大きくなるように各係数の値が調節される。モデル調節要素242は、補正後基本周波数f0dと、係数α0>0及びβ0>0とに基づき、関係式023にしたがって第1時定数τ1の値を調節又は設定する。τ1=α0f0dβ0。モデル調節要素242は、第1時定数τ1と、正係数γ0たとえば2とに基づき、関係式024にしたがって第2時定数τ2の値を調節又は設定する。τ2=γ0τ1。関係式023により定義される関数は、図6に示されているような変化特性を示す。このため、補正後基本周波数f0dが高くなる歩行周期が短くなるにつれて、第1時定数τ1及び第2時定数τ2の値が減少するように調節される。また、補正後基本周波数f0dが低くなる歩行周期が長くなるにつれて、第1時定数τ1及び第2時定数τ2の値が増加するように調節される。続いて、状態振動子生成要素24が、連立微分方程式010により定義される状態振動子モデルに対して、入力波形信号として運動振動子出力要素21から出力される運動振動子を入力することにより、出力波形信号として状態振動子ξを生成する図3/STEP5。状態振動子ξには、左屈曲振動子ξL、左伸展振動子ξL、右屈曲振動子ξR及び右伸展振動子ξRが含まれる。歩行補助装置10の動作開始から、非対称度sの評価に必要な運動振動子が得られるまでの期間において、各係数ciは予め定められている初期値に設定される。同様に、歩行補助装置10の動作開始から、基本周波数f0の推定に必要な差分振動子後述が得られるまでの期間において、固有角速度ωは予め定められている初期値に設定される。状態振動子モデルに入力される運動振動子φは、非対称度sの評価に用いられた運動振動子φと同一であっても異なっていてもよい。たとえば、状態振動子モデルに入力される運動振動子φとして、人間の左右の股関節角度の時間変化態様を表わす波形信号に代えて、左右の股関節角速度、左右の肩関節角度又は左右の肩関節角速度等、人間の歩行周期と緊密に関連した周期で変化する左右一対の変数の時間変化態様を表わす波形信号が採用されてもよい。肩関節角度は、股関節角度と同様に、ロータリーエンコーダにより構成されている一対の肩関節角度センサが人間の左右の肩の外側に配置され、当該肩関節角度センサの出力に基づいて測定されうる。状態振動子ξiは、状態変数uiの値が閾値uth未満である場合は0、状態変数uiの値が閾値uth以上である場合はこのuiの値をとる連立微分方程式010参照。これにより、左大腿の屈曲運動状態を表わす状態変数uLが大きくなると左屈曲振動子ξLの振幅が左伸展振動子ξLよりも大きくなる。また、右大腿の屈曲運動状態を表わす状態変数uRが大きくなると右屈曲振動子ξRの振幅が右伸展振動子ξRの振幅よりも大きくなる。さらに、左大腿の伸展運動状態を表わす状態変数uLが大きくなると左伸展振動子ξLの振幅が左屈曲振動子ξLよりも大きくなる。また、右大腿の伸展運動状態を表わす状態変数uRが大きくなると右伸展振動子ξRの振幅が右屈曲振動子ξRの振幅よりも大きくなる。結果として、連立微分方程式010に含まれる左屈曲係数cL及び左伸展係数cLのそれぞれの値に応じた振幅と、第1左時定数τ1Lの値に応じた周波数とにしたがって変化する左屈曲振動子ξL及び左伸展振動子ξLが生成される。また、連立微分方程式010に含まれる右屈曲係数cR及び右伸展係数cRのそれぞれの値に応じた振幅と、第1右時定数τ1Rの値に応じた周波数とにしたがって変化する右屈曲振動子ξR及び右伸展振動子ξRが生成される。そして、制御振動子生成要素25が状態振動子ξに基づき、関係式040にしたがって制御振動子η=ηL,ηRを設定する図3/STEP6。ηL=χLξL-χLξL,ηR=χRξR-χRξR。左制御振動子ηLは、左屈曲振動子ξL及び係数χLの積と、左伸展振動子ξL及び係数χLの積との差として算定される。右制御振動子ηRは、右屈曲振動子ξR及び係数χRの積と、右伸展振動子ξR及び係数χRの積との差として算定される。当該4つの係数χiは同じ値に設定されてもよい。図4に示されているように左脚の歩幅が右脚の歩幅よりも小さくなっている状況では、図5bに示されているように左屈曲係数cLの値が右屈曲係数cRの値よりも大きく、かつ、右伸展係数cRの値が左伸展係数cLの値よりも大きく設定される。これにより、図5cに示されているように左制御振動子ηLの屈曲振幅|χLξL|が右制御振動子ηRの屈曲振幅|χRξR|よりも大きくなる。また、左制御振動子ηLの伸展振幅|χLξL|が右制御振動子ηRの伸展振幅|χRξR|よりも小さくなる。そして、制御装置20により制御振動子ηに基づいてバッテリ16から左右のアクチュエータ14L,14Rにそれぞれ供給される電流I=IL,IRが調節される。これにより、第1装具11及び第2装具12を介して上体に対する左右の大腿の屈曲運動及び伸展運動を補助する力又は股関節回りの回転力F=FL,FRが調節される。補助力Fは電流Iに基づき、たとえばFt=GIG:比例係数と表現される。エージェントの歩行運動はトレッドミルの上で実施されてもよい。その後、運転スイッチがONからOFFに切り替えられたこと、又は、動作異常が検知されたこと等の動作終了条件が満たされたか否かが判定される図3/STEP7。そして、当該判定結果が否定的である場合図3/STEP7..NO、前記一連の処理が繰り返される一方、当該判定結果が肯定的である場合図3/STEP7..YES、前記一連の処理が終了する。基本周波数の推定基本周波数推定要素23による基本周波数f0の推定処理図3/STEP3参照について説明する。まず、差分振動子出力要素231が、左右の股関節角度センサ202R,202Lの出力に基づき、股関節角度の差分の時間変化を表わす波形信号を差分振動子として出力する図7/STEP31。これにより、図8aに示されているような波形信号が差分振動子として得られる。なお、股関節角度の差分に代えて肩関節角度の差分の時間変化を表わす波形信号が差分振動子として出力されてもよい。続いて、窓処理要素232が、指定期間にわたる差分振動子に窓を掛ける図7/STEP32。窓関数としてはハン窓が用いられ、これにより図8bに破線で示されている指定期間にわたる差分振動子が、実線で示されているように変形される。窓関数としては、ハン窓のほか、矩形窓、ガウス窓、ハミング窓、ブラックマン窓、カイザー窓、バートレット窓又は指数窓などのさまざまな窓関数が用いられうる。窓処理要素232は、スペクトル解析処理要素234により決定された前回基本周波数f0に基づき、当該基本周波数f0を変数とする減少関数にしたがって、今回の窓の幅wを設定する。この減少関数はたとえば関係式050により定義される。係数w01及びw02は歩行周期が少なくとも2周期だけ含まれるように設定されている。なお、窓処理要素232は、初回の窓幅wとして、歩行周期が少なくとも2周期含まれるような定数を採用する。w=w01f0w02。窓処理要素232は、差分振動子波形信号をダウンサンプリングすることによって、差分振動子から第1指定周波数を超える高周波数成分を除去した上で、窓処理を実行してもよい。窓処理要素232は、差分振動子波形信号をハイパスフィルタに通すことにより、差分振動子から第2指定周波数以下の低周波成分を除去した上で、窓処理を実行してもよい。周波数解析処理要素233は、窓処理された今回指定期間における差分振動子をFFT等、周波数解析処理することにより、パワースペクトルを作成する図7/STEP33。これにより、図8cに示されているように周波数fが横軸で表わされ、周波数解析結果としてのパワーが縦軸で表わされているパワースペクトルが生成される。スペクトル解析処理要素234は、閾値以上の高さを有するとともに、最低の周波数帯域にあるピークの位置を表わす周波数を基本周波数f0として決定又は推定する図7/STEP34。閾値は、ノイズを除去する観点から、パワースペクトルにおけるピーク最高値の0.05~0.20倍の範囲の値に設定される。これにより、図8cに示されているように、閾値pth以上の高さを有するピークP1及びP2のうち、最低の周波数帯域にあるピークP1の位置を表わす周波数が基本周波数f0として決定される。さらに、図8dに示されている曲線を表わす関係式060にしたがって基本周波数f0が補正されることにより、補正後基本周波数f0dが決定される。f0d=f0fth。歩行補助装置の作用効果前記機能を発揮する歩行補助装置10によれば、人間の左右の運動形態の不均等の度合が、左運動振動子φL及びφRの非対称度sに現われることに鑑みて、当該非対称度sの低減が図られる。その結果、人間の左右の運動形態の均等を図りながら、当該人間の歩行運動が補助されうる。たとえば、図4に示されているように左脚の歩幅が右脚の歩幅よりも小さくなっている状況では、図5cに示されているように変化する左制御振動子ηL及び右制御振動子ηRが生成される。歩幅は、利用者である人間の各大腿の長さ及び必要に応じて各下腿の長さメモリに格納される。と、股関節角度センサ202の出力により表わされる左右の股関節角度の最大値及び最小値とに基づき、幾何学的関係にしたがって算出されうる。これにより、右脚が立脚状態足が着床している状態にあるとき、遊脚状態足が離床している状態にある左脚の前方への踏み出し量の増加、ひいては左脚の歩幅の増加が図られる。一方、左脚が立脚状態にあるとき、遊脚状態にある右脚の前方への踏み出し量の維持又は低減、ひいては右脚の歩幅の維持又は低下が図られる。その結果、左脚については現状よりも大股での歩行が誘導される一方、右脚については現状と同じ又は現状よりも小股での歩行が誘導されるように、歩行補助装置10がその歩行運動を補助するので、人間の左右の脚の歩幅及び歩行率の均等化が図られる。また、左屈曲係数cL及び右伸展係数cRが同値になり、かつ、右屈曲係数cR及び左伸展係数cLが同値になるように、各係数ciの値が調節される関係式021参照。これにより、左制御振動子ηLの屈曲振幅|χLξL|及び右制御振動子ηRの伸展振幅|χRξR|を一致させ又は近似させ、かつ、右制御振動子ηRの屈曲振幅|χRξR|及び左制御振動子ηLの伸展振幅|χLξL|を一致させ又は近似させることができる図5c参照。このため、立脚状態にある一方の脚の伸展量に鑑みて、遊脚状態にある他方の脚の屈曲量が著しく乖離したものになり、歩行補助装置10の利用者である人間に違和感を覚えさせる事態が回避されうる。本発明の他の実施形態左総合振幅|φL|と右総合振幅|φR|との偏差に代えて、左屈曲振幅|φL|と右屈曲振幅|φR|との偏差、左伸展振幅|φL|と右伸展振幅|φR|との偏差、若しくは、左屈曲振幅|φL|及び右伸展振幅|φR|の和と、右屈曲振幅|φR|及び左伸展振幅|φL|の和との偏差又は当該偏差の複数周期にわたる平均値が非対称度sとして評価されてもよい。左屈曲係数cL、左伸展係数cL、右屈曲係数cR及び右伸展係数cRのうち一部のみが補正量f-s又はfsに応じて補正されてもよい関係式021022参照。たとえば、左屈曲係数cL及び右伸展係数cRの組と、左伸展係数cL及び右屈曲係数cRの組とのうち、一方の組に属する係数のみが補正量fsに応じて増減調節されることにより、他方の組に属する係数との相対的な偏差が調節されてもよい。身体機能の回復という観点から、一方の脚の歩幅の増加のみにより、両脚の歩幅の均等化が図られてもよい。その一方、非対称度sが顕著に高い場合、前記のように一方の脚の歩幅の増加に加えて、他方の脚の歩幅の低減によって、両脚の歩幅の均等化又は歩み寄りが図られることが好ましい。この観点から、非対称度sが閾値を超えているか否かが判定され、非対称度sが閾値未満である場合には一方の脚の歩幅の増加のみにより、両脚の歩幅の均等化が図られ、非対称度sが閾値以上である場合には両脚の歩幅の歩み寄りが図られてもよい。係数D関係式022参照の値が変更されることにより、非対称度sに応じた左屈曲係数cL及び右伸展係数cRの補正量f-sと、左伸展係数cL及び右屈曲係数cRの補正量fsとの比率が変更されてもよい。係数Dがsを変数とする関数であってもよい。制御装置20は、左屈曲係数cL、左伸展係数cL、右屈曲係数cR及び右伸展係数cRのそれぞれの値の大小に応じて、人間の左大腿の屈曲運動及び伸展運動並びに右大腿の屈曲運動及び伸展運動のそれぞれを補助するためのアクチュエータ14の動作振幅の大小を制御するように構成されていれば、その構成は前記実施形態の構成に限定されない。アクチュエータ14L,14Rの動作周期は、固有角速度ωL、ωRに応じて2π/ωL、2π/ωRに一致するように制御されてもよい。たとえば、制御装置20は、左屈曲係数cL、左伸展係数cL、右屈曲係数cR及び右伸展係数cRのそれぞれの値を、左股関節角度の屈曲側目標値、左股関節角度の伸展側目標値、右股関節角度の屈曲側目標値及び右股関節角度の伸展側目標値のそれぞれとして、PID制御等のフィードバック制御則にしたがって、左右の股関節角度を制御するように構成されていてもよい。この場合、運動振動子出力要素21、状態振動子生成要素24及び制御振動子生成要素25は省略されてもよい。基本周波数推定要素23が省略されてもよい。特許第3930399号、特許第3950149号、特許第4008464号又は特許第4271711号等の公報に記載されているように、ファンデルポル方程式等により定義される第1モデルにしたがって第1振動子が生成され、かつ、当該第1振動子に基づいて固有角速度ω連立微分方程式010参照が設定された上で、当該固有角速度ωの逆数に比例するように時定数τ1が設定されてもよい。当該公報における第2振動子が本願発明における状態振動子に相当する。制御振動子ηは特許第4271711号等の公報に記載されているように、仮想的な弾性要素による弾性力及び仮想的な減衰要素による減衰力のうち一方又は両方を表わすように生成されてもよい。基本周波数推定要素23は、左運動振動子φLの周波数、右運動振動子φRの周波数、左運動振動子φLの周波数と右運動振動子φRの周波数との平均周波数、又は、当該周波数の複数周期にわたる平均値を基本周波数f0として推定するように構成されていてもよい。"