JP3155154B2 - Wheel information estimation device - Google Patents

Wheel information estimation device

Info

Publication number
JP3155154B2
JP3155154B2 JP17125894A JP17125894A JP3155154B2 JP 3155154 B2 JP3155154 B2 JP 3155154B2 JP 17125894 A JP17125894 A JP 17125894A JP 17125894 A JP17125894 A JP 17125894A JP 3155154 B2 JP3155154 B2 JP 3155154B2
Authority
JP
Japan
Prior art keywords
value
correction
estimated
disturbance
wheel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP17125894A
Other languages
Japanese (ja)
Other versions
JPH0834215A (en
Inventor
秀樹 大橋
弘之 河井
弘義 小島
孝治 梅野
勝宏 浅野
俊治 内藤
伸好 小野木
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Denso Corp
Toyota Motor Corp
Toyota Central R&D Labs Inc
Original Assignee
Denso Corp
Toyota Motor Corp
Toyota Central R&D Labs Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Denso Corp, Toyota Motor Corp, Toyota Central R&D Labs Inc filed Critical Denso Corp
Priority to JP17125894A priority Critical patent/JP3155154B2/en
Publication of JPH0834215A publication Critical patent/JPH0834215A/en
Application granted granted Critical
Publication of JP3155154B2 publication Critical patent/JP3155154B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【産業上の利用分野】本発明は、車輪の回転速度に基づ
いて車輪に関する情報を推定する車輪情報推定装置に関
するものであり、特に、車輪速度の検出値に混入した周
期的外乱による推定精度の低下を抑制する技術に関する
ものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a wheel information estimating apparatus for estimating information on a wheel based on the rotational speed of the wheel, and more particularly, to a method for estimating the accuracy of estimation due to a periodic disturbance mixed in a detected value of a wheel speed. The present invention relates to a technique for suppressing the reduction.

【0002】[0002]

【従来の技術】上記車輪情報推定装置は一般に、(a) 車
輪の回転速度を検出する検出装置と、(b) その検出装置
による検出値に基づき、車輪に関する情報を推定する推
定手段とを含むように構成される。
2. Description of the Related Art The wheel information estimating device generally includes (a) a detecting device for detecting a rotational speed of a wheel, and (b) estimating means for estimating information on a wheel based on a value detected by the detecting device. It is configured as follows.

【0003】この車輪情報推定装置の一形式が特開平5
−133831号公報に記載されている。これは、検出
装置により検出された回転速度の複数の周波数成分のう
ち設定周波数範囲内において強度が実質的に最大となる
ものの周波数に基づいてタイヤ空気圧を車輪に関する情
報として推定するものであり、その周波数が低いほどタ
イヤ空気圧が低いと推定するものである。
One type of the wheel information estimating apparatus is disclosed in
-133831. This is to estimate the tire air pressure as information about the wheel based on the frequency of the one whose intensity is substantially maximum within a set frequency range among a plurality of frequency components of the rotational speed detected by the detection device, It is estimated that the lower the frequency, the lower the tire pressure.

【0004】[0004]

【発明が解決しようとする課題】本出願人はこの車輪情
報推定装置について研究を行い、その結果、次のような
事実を得た。検出装置による検出値には車輪のタイヤ不
均一性等に基づく周期的外乱が混入しており、この周期
的外乱をそのまま放置したのではその周期的外乱による
影響が推定手段による推定値に現れてしまい、推定精度
が低下するおそれがあるという事実を得たのである。そ
こで、本出願人は、検出値に混入した周期的外乱を除去
するため、前記検出装置および推定手段を含む車輪情報
推定装置に、検出値を補正対象として、検出値に混入し
た周期的外乱が除去されるように検出値を補正する検出
値補正手段を付加することを提案した。
The present applicant has conducted research on the wheel information estimating apparatus, and as a result, has obtained the following facts. The values detected by the detection device include periodic disturbances based on the tire non-uniformity of the wheels, etc. If the periodic disturbances are left as they are, the effects of the periodic disturbances appear in the estimation values by the estimation means. As a result, the fact that the estimation accuracy may be reduced has been obtained. In order to eliminate the periodic disturbance mixed in the detection value, the present applicant sets the wheel information estimation device including the detection device and the estimating means as a correction target for the detection value and the periodic disturbance mixed in the detection value. It has been proposed to add a detection value correction means for correcting the detection value so as to be removed.

【0005】しかし、この検出値補正手段は具体的に
は、周期的外乱の周期性を利用し、車輪が複数回回転す
る間に取得された検出値を用いることによって周期的外
乱を除去するものである。そのため、この検出値補正手
段による補正誤差が許容値以下となるまでに時間がかか
り、その間、検出値に混入した周期的外乱によって推定
精度が低下することを避け得ないという問題がある。
However, the detection value correcting means specifically removes the periodic disturbance by utilizing the periodicity of the periodic disturbance and using the detection values acquired during a plurality of rotations of the wheel. It is. Therefore, there is a problem that it takes time until the correction error by the detection value correction means becomes equal to or less than the allowable value, and during that time, it is inevitable that the estimation accuracy deteriorates due to the periodic disturbance mixed into the detection value.

【0006】このような事情に鑑み、請求項1の発明
は、検出値に混入した周期的外乱による推定誤差が車両
の走行状況に強く関連するという事実に着目し、検出値
ではなく推定値を補正する推定値補正手段を前記検出値
補正手段と併用することにより、検出値補正手段の作用
完了前でも周期的外乱による推定精度の低下を抑制可能
とすることを課題としてなされたものである。
In view of such circumstances, the invention of claim 1 pays attention to the fact that an estimation error due to a periodic disturbance mixed in a detection value is strongly related to a traveling state of a vehicle, and calculates an estimation value instead of a detection value. An object of the present invention is to make it possible to suppress a decrease in estimation accuracy due to a periodic disturbance even before the operation of the detection value correcting means is completed, by using the estimated value correcting means for correction together with the detection value correcting means.

【0007】請求項2の発明は、その請求項1の発明に
おいて、車両の走行状況のうち特に車体速度と路面凹凸
度という要素が周期的外乱による推定誤差に強く関連す
るという事実に着目し、車両の走行状況に基づく推定値
補正を効果的に実行可能とすることを課題としてなされ
たものである。
A second aspect of the present invention is based on the first aspect of the invention, which focuses on the fact that, among the driving conditions of the vehicle, the elements of the vehicle body speed and the degree of road surface irregularity are strongly related to the estimation error due to the periodic disturbance. An object of the present invention is to make it possible to effectively execute an estimated value correction based on a traveling state of a vehicle.

【0008】請求項3の発明は、請求項1または2の発
明における推定値補正手段を前記検出値補正手段と併用
する車輪情報推定装置において、両補正手段の作用が重
複して推定精度が悪化する事態を回避することを課題と
してなされたものである。
According to a third aspect of the present invention, there is provided a wheel information estimating apparatus in which the estimated value correcting means according to the first or second aspect of the present invention is used in combination with the detected value correcting means. The purpose of this is to avoid the situation of doing so.

【0009】[0009]

【課題を解決するための手段】それぞれの課題を解決す
るために、請求項1の発明は、図19に概念的に表すよ
うに、前記検出装置1および推定手段2を含む車輪情報
推定装置において、(a) 検出装置1による検出値を、そ
れに混入した周期的外乱が除去されるように補正して推
定手段2に供給する検出値補正手段3と、(b) 推定手段
2による推定値を、車両の走行状況に基づき、検出値補
正手段3により除去されないで検出装置1による検出値
に残る周期的外乱による影響が推定値に現れないように
補正する推定値補正手段4とを設けたことを特徴とす
る。
In order to solve each of the problems, the invention of claim 1 is directed to a wheel information estimating apparatus including the detecting device 1 and the estimating means 2 as conceptually shown in FIG. (A) The value detected by the detection device 1 is
And correct so that the periodic disturbance mixed into the
Detection value correction means 3 to be supplied to the determination means 2, and (b) estimation means
2 is added to the detected value based on the running condition of the vehicle.
The value detected by the detector 1 without being removed by the corrector 3
And an estimated value correcting means 4 for correcting the influence of the periodic disturbance remaining in the estimated value to not appear in the estimated value.

【0010】 なお、 「推定手段2」は例えば、前記公報
に記載されているように、検出装置1により検出された
回転速度の複数の周波数成分のうち設定周波数範囲内に
おいて強度が実質的に最大となるものの周波数に基づい
て車輪のタイヤ空気圧を車輪情報として推定する形式と
することができる。
[0010] Incidentally, "estimating means 2" for example, as described in the publication, the maximum intensity is substantially within the set frequency range out of a plurality of frequency components of the rotational speed detected by the detection device 1 The tire pressure of the wheel can be estimated as wheel information based on the frequency of the following.

【0011】 「推定手段2」はまた、 (a) 少なくとも、
車輪情報の基礎値である車輪情報基礎値と検出装置によ
り検出された回転速度とから、車輪に対する外乱を推定
する外乱オブザーバと、(b) 推定された外乱に基づき、
車輪情報の実際値である車輪情報実際値の、車輪情報基
礎値からの変化量を推定する車輪情報変化量推定部とを
含む形式とすることもできる。この形式については実施
例において詳述する。
[0011] "estimated means 2" is also, at least (a),
From a wheel information basic value, which is a basic value of wheel information, and a rotation speed detected by the detection device, a disturbance observer that estimates a disturbance to the wheel, and (b) based on the estimated disturbance,
It is also possible to adopt a form including a wheel information change amount estimating unit that estimates a change amount of the wheel information actual value, which is the actual value of the wheel information, from the wheel information base value. This format will be described in detail in an embodiment.

【0012】 なお、ここにおける「検出値補正手段
は前記のように、周期的外乱の周期性を利用し、車輪が
複数回回転する間に取得された検出値を用いることによ
って周期的外乱を除去する態様とすることが可能であ
り、例えば、次のような態様とすることができる。すな
わち、検出装置からの周期的信号に基づいて順に決定さ
れた過去の複数の回転速度のうち今回決定された回転速
度と予め定められた関係を有するものを今回決定された
回転速度に加算または減算することを含む演算により、
回転速度に混入した周期的外乱を除去する態様とするこ
とができる。
[0012] It should be noted that, here in the definitive "detection value correcting means 3"
As described above, utilizing the periodicity of the periodic disturbance, it is possible to remove the periodic disturbance by using a detection value obtained during multiple rotations of the wheel, for example, The following embodiment can be adopted. That is, among a plurality of past rotational speeds determined in order based on the periodic signal from the detection device, a rotational speed having a predetermined relationship with the rotational speed determined this time is added to or subtracted from the rotational speed determined this time. Operations that include
A mode in which periodic disturbance mixed in the rotation speed is removed can be adopted.

【0013】 さらに具体的に説明すれば、予め定めら
れたサンプリング周期が経過する毎に、その間に検出装
置からの周期的信号に基づいて演算された少なくとも1
個の元回転速度の平均値をサンプル回転速度として演算
するサンプル回転速度演算部と、今回演算されたサン
プル回転速度の変動成分をサンプル変動回転速度として
演算するサンプル変動回転速度演算部と、除去すべき
周期的外乱の周期的変化を表す誤差値(変動回転速度に
加算することによってその変動回転速度に混入した周期
的外乱を除去するもの)を少なくとも車輪1回転分、車
輪の回転位置に関連付けて記憶する1回転分メモリと、
サンプリング周期が経過する毎に1回転分メモリか
ら、今回のサンプル変動回転速度よりちょうど車輪1回
転前に取得された少なくとも1個の誤差値を旧誤差値群
として読み出し、その旧誤差値群の平均値である平均旧
誤差値と今回のサンプル変動回転速度との和を演算して
出力値とする出力値演算部と、その出力値とサンプル
変動回転速度の固定成分との和を最終的な回転速度とし
て演算する最終的回転速度演算部と、今回のサンプリ
ング変動回転速度よりちょうど車輪1回転前に取得され
た平均旧誤差値から今回の前記出力値を減算することに
より、今回の新誤差値を演算する誤差値演算部と、1
回転分メモリにおいて、今回の旧誤差値群の各値をすべ
て、演算された新誤差値に更新する1回転分メモリ更新
部とを含む態様とすることができるのである。
[0013] In more detail, each time a sampling cycle set in advance has elapsed, at least calculated based on the periodic signal from the detection device during 1
A sample rotation speed calculation unit that calculates the average value of the original rotation speeds as the sample rotation speed, a sample fluctuation rotation speed calculation unit that calculates the fluctuation component of the sample rotation speed calculated this time as the sample fluctuation rotation speed, and An error value representing a periodic change of the power periodic disturbance (which removes the periodic disturbance mixed in the variable rotational speed by adding to the variable rotational speed) is associated with at least one rotation of the wheel and the rotational position of the wheel. One rotation memory for storing,
Every time the sampling period elapses, at least one error value obtained just before one rotation of the wheel from the current sample variation rotation speed is read from the memory for one rotation as an old error value group, and the average of the old error value group is read. The output value calculation unit calculates the sum of the average old error value, which is the current value, and the current sample fluctuation rotation speed to obtain an output value, and outputs the sum of the output value and the fixed component of the sample fluctuation rotation speed to the final rotation. A final rotation speed calculation unit that calculates the speed, and a subtraction of the current output value from an average old error value obtained just before one rotation of the wheel from the current sampling fluctuation rotation speed, thereby obtaining the new error value of the current time. An error value calculator for calculating,
In this case, the rotation memory may include a one-rotation memory updating unit that updates all the values of the old error value group at this time with the calculated new error value.

【0014】 なお、ここにおける「車輪情報」には例え
ば、タイヤの空気圧,半径,慣性モーメント,ばね定
数,ダンパ係数,接地性,コーナリングパワー等を選ぶ
ことができる。
The "wheel information" here can be selected, for example, from tire pressure, radius, inertia moment, spring constant, damper coefficient, ground contact, cornering power, and the like.

【0015】 請求項2の発明は、請求項1の発明におい
て、推定値補正手段を、推定手段2による推定値の補
正量を走行状況としての車体速度が大きい場合において
小さい場合におけるより大きく決定することと、走行状
況としての路面凹凸度が低い場合において高い場合にお
けるより大きく決定することとの少なくとも一方を行う
ものとしたことを特徴とする。
According to a second aspect of the present invention, in the first aspect of the present invention, the estimated value correcting means 4 determines the correction amount of the estimated value by the estimating means 2 to be larger when the vehicle speed as the traveling condition is high and when the vehicle speed is small. And when the road surface unevenness degree as the running condition is low, the determination is made larger when the road surface unevenness degree is high.

【0016】 請求項3の発明は、請求項1または2の発
明において、図20に概念的に表すように、さらに、検
出値補正手段3により除去されないで検出装置1による
検出値に残る周期的外乱が許容値より大きい間は推定値
補正手段による補正を許可し、許容値以下になった後
には禁止する推定値補正許可・禁止手段5設けたこと
を特徴とする。
A third aspect of the present invention is the invention according to the first or second aspect, further comprising a detection step as conceptually shown in FIG.
By the detection device 1 without being removed by the output value correction means 3
Estimated value correction permitting / prohibiting means 5 for permitting the correction by the estimated value correcting means 4 while the periodic disturbance remaining in the detected value is larger than the allowable value, and prohibiting the correction when the periodic disturbance is less than the allowable value is provided. I do.

【0017】[0017]

【作用】請求項1の発明に係る車輪情報推定装置におい
ては、検出値補正手段3が検出値を、それに混入した周
期的外乱が除去されるように補正し、さらに、検出値に
混入した周期的外乱による推定誤差が車両の走行状況に
強く関連するという事実に着目し、推定値補正手段
が、推定手段2による推定値を車両の走行状況に基づい
て補正し検出装置1による検出値に混入した周期的外
乱の影響が推定値に現れないようにする。
In the wheel information estimating apparatus according to the first aspect of the present invention,
TheThe detection value correction means 3 converts the detection value into the
Correction to eliminate periodic disturbances,Detection value
Estimation errors due to mixed periodic disturbances may cause
Paying attention to the fact that they are strongly related,4
Calculates the estimated value by the estimating means 2 based on the traveling state of the vehicle.
Correct,Periodic outflow mixed with the value detected by the detection device 1
Estimate the effect of disturbanceNot to appear inI do.

【0018】本出願人の研究によれば、検出値に混入し
た周期的外乱が検出値から十分に除去されない状態で
は、車体速度が大きいほど推定誤差が大きくなり、ま
た、路面凹凸度が低いほど推定誤差が大きくなるという
傾向があることが判明した。車体速度が大きいほど、ま
た、路面凹凸度が低いほど、推定値に占める周期的外乱
による誤差値の割合が増加してしまうのである。 そこ
で、請求項2の発明に係る車輪情報推定装置において
は、推定値補正手段が、推定手段2による推定値の補
正量を走行状況としての車体速度が大きい場合において
小さい場合におけるより大きく決定することと、走行状
況としての路面凹凸度が低い場合において高い場合にお
けるより大きく決定することとの少なくとも一方を行
う。
According to the research of the present applicant, it has
In the state where the periodic disturbance is not sufficiently removed from the detected value
The higher the vehicle speed, the greater the estimation error.
Also, the lower the road surface roughness, the greater the estimation error
Turned out to be prone. The higher the vehicle speed, the better
In addition, the lower the degree of unevenness of the road surface, the more the periodic disturbance
Therefore, the ratio of the error value due to is increased. There
In the wheel information estimating device according to the second aspect of the present invention,
Is the estimated value correction means4Is the complement of the estimated value by the estimating means 2.
When the vehicle speed is high with the positive amount as the driving situation
Deciding bigger when small and running
When the road surface roughness is low or high
Make at least one of
U.

【0019】それら請求項1または2の発明における
「推定値補正手段」は一般に、車両の走行状況に基づ
いて推定値をオープンループ方式で補正する態様とされ
るため、補正速度(補正誤差を許容値まで低下させるの
にかかる時間に相当する)を向上させることが容易であ
るという利点がある反面、補正精度を十分に高めること
が困難であるという欠点がある。これに対し、前記検出
値補正手段を例えば、周期的外乱の周期性を利用し、
車輪が複数回回転する間に取得された検出値を用いるこ
とによって周期的外乱を除去する態様とする場合には、
補正精度を十分に高めることが容易であるという利点が
ある反面、補正速度を向上させることが困難であるとい
う欠点がある。要するに、推定値補正手段と検出値補
正手段とは一方の欠点を他方の利点とする関係にある
のであり、一方の欠点が他方の利点で補われるように両
者を併用すれば、補正速度の低下を抑制しつつ補正精度
を十分に高めることが可能となる。
Since the "estimated value correcting means 4 " in the first or second aspect of the present invention is generally adapted to correct the estimated value based on the running condition of the vehicle in an open loop manner, the correction speed (correction error is corrected). This is advantageous in that it is easy to improve the correction accuracy (equivalent to the time required to reduce the value to the allowable value), but has the disadvantage that it is difficult to sufficiently increase the correction accuracy. On the other hand, the detection value correction means 3 uses, for example, the periodicity of a periodic disturbance,
In the case of removing the periodic disturbance by using the detection value obtained during the rotation of the wheel a plurality of times,
There is an advantage that it is easy to sufficiently increase the correction accuracy, but there is a disadvantage that it is difficult to improve the correction speed. In short, the estimated value correcting means 4 and the detected value correcting means 3 have a relationship in which one of the disadvantages is used as the other advantage. , It is possible to sufficiently increase the correction accuracy while suppressing the decrease.

【0020】しかしながら、両補正手段3,4を併用す
る場合には次のような問題が生ずる。すなわち、検出値
補正手段の作用が完了した後にも推定値補正手段
作用することがあると、検出値補正手段によって周期
的外乱が十分に除去された検出値に基づいて推定手段2
により推定された推定値であってもはや周期的外乱の影
響がほとんど現れていないものに対して推定値補正手段
が作用してしまい、検出値補正手段の作用と推定値
補正手段の作用とが相互に干渉し、かえって推定精度
が低下してしまうという問題があるのである。
However, when both correction means 3 and 4 are used together, the following problem occurs. That is, if the estimated value correcting means 4 may operate even after the operation of the detected value correcting means 3 is completed, the estimating means 2 based on the detected value from which the periodic disturbance has been sufficiently removed by the detected value correcting means 3 .
Value correction means for the estimated value estimated by
4 acts, and the operation of the detection value correction means 3 and the operation of the estimation value correction means 4 interfere with each other, and the estimation accuracy is rather reduced.

【0021】この問題を解決するために請求項3の発明
検出値補正手段3と推定値補正手段4とを併用する形
式の車輪情報推定装置に対してなされたのであり、本発
明に係る車輪情報推定装置においては、推定値補正許可
・禁止手段5が、検出値補正手段3により除去されない
で検出装置1による検出値に残る周期的外乱、すなわ
ち、補正誤差が許容値より大きい間は推定値補正手段
による補正を許可し、許容値以下になった後には禁止す
る。
In order to solve this problem, the invention of claim 3 has been made for a wheel information estimating apparatus of a type using both the detected value correcting means 3 and the estimated value correcting means 4. In the information estimation device, the estimation value correction permission / prohibition unit 5 is not removed by the detection value correction unit 3.
The periodic disturbance remaining in the value detected by the detection device 1
That is, while the correction error is larger than the allowable value, the estimated value correcting means 4
Is permitted, and is prohibited after the value becomes equal to or less than the allowable value.

【0022】[0022]

【発明の効果】以上の説明から明らかなように、請求項
1の発明によれば、検出値補正手段と推定値補正手段と
が併用されるから、それらの一方の欠点が他方の利点で
補われ、それにより、補正速度の低下を抑制しつつ補正
精度を十分に高めることが可能となるという効果が得ら
れる。
As is apparent from the above description, according to the first aspect of the present invention, the detection value correcting means and the estimated value correcting means
Are used together, so the disadvantages of one of them are the advantages of the other.
Compensation, thereby compensating while suppressing a decrease in compensation speed.
The effect is obtained that the accuracy can be sufficiently increased .

【0023】また、請求項2の発明によれば、車両の走
行状況のうち周期的外乱による推定誤差に強く関連する
要素である車体速度と路面凹凸度との少なくとも一方に
基づいて推定値が補正されるため、車両の走行状況に基
づく推定値の補正を効果的に実行可能となるという効果
が得られる。
According to the second aspect of the present invention, the estimated value is corrected on the basis of at least one of the vehicle speed and the degree of road surface unevenness, which are factors strongly related to the estimation error due to the periodic disturbance in the running condition of the vehicle. Therefore, it is possible to effectively execute correction of the estimated value based on the traveling state of the vehicle.

【0024】また、請求項3の発明によれば、検出値補
正手段と推定値補正手段と併用されるにもかかわら
ず、両補正手段の相互干渉が回避され、補正の適正化が
確保されるという効果が得られる。
According to the third aspect of the present invention, the detection value correcting means and the estimated value correcting means are used in combination.
Therefore, the mutual interference between the two correction means is avoided, and the effect of ensuring the appropriate correction is obtained.

【0025】[0025]

【実施例】以下、発明を図示の実施例に基づいて具体
的に説明する。図2において10はロータ、12は電磁
ピックアップである。ロータ10は図3に示す車輪14
と共に回転するものであり、外周に多数の歯16を備え
ている。電磁ピックアップ12はそれらの歯16の通過
に応じて周期的に変化する電圧を発生する。この電圧は
波形整形器18によって矩形波に整形され、コンピュー
タ20のI/Oポート22に供給される。車輪14は4
個あり、それらに設けられている各電磁ピックアップ1
2が全て波形整形器18を経てコンピュータ20に接続
されるが、図2には代表的に1組のみが図示されてい
る。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, the present invention will be specifically described based on the illustrated embodiments. In FIG. 2, reference numeral 10 denotes a rotor, and 12 denotes an electromagnetic pickup. The rotor 10 is a wheel 14 shown in FIG.
, And has a number of teeth 16 on the outer circumference. The electromagnetic pickup 12 generates a voltage that changes periodically as the teeth 16 pass. This voltage is shaped into a rectangular wave by the waveform shaper 18 and supplied to the I / O port 22 of the computer 20. Wheel 14 is 4
And each electromagnetic pickup 1 provided on them.
2 are all connected to a computer 20 via a waveform shaper 18, but only one set is typically shown in FIG.

【0026】車輪14は図3に示すように、ホイール2
4の外周にタイヤ26が取り付けられたタイヤ付ホイー
ルであるが、図4に示すように、相対回転可能なリム側
部28とベルト側部30とがねじりばね32によって連
結されたものと考えることができる。上記ロータ10は
ホイール24と一体的に回転するように取り付けられる
ため、電磁ピックアップ12は厳密にはリム側部28の
角速度を検出することになる。
The wheel 14 is, as shown in FIG.
4 is a tire-equipped wheel having a tire 26 mounted on the outer periphery thereof. As shown in FIG. 4, it is considered that a rim side portion 28 and a belt side portion 30 that are relatively rotatable are connected by a torsion spring 32. Can be. Since the rotor 10 is mounted so as to rotate integrally with the wheel 24, the electromagnetic pickup 12 detects the angular velocity of the rim side portion 28 strictly.

【0027】コンピュータ20は図2に示すように処理
装置としてのCPU40,第一記憶装置としてのROM
42および第二記憶装置としてのRAM44を備えてお
り、ROM42に図5のフローチャートで表される回転
速度演算・補正ルーチンを始めとする種々の制御プログ
ラムが格納されることによって、図1に示すリム側部回
転速度演算・補正部45を構成している。このコンピュ
ータ20は別のコンピュータ47と接続されている。こ
のコンピュータ47は処理装置としてのCPU48,第
一記憶装置としてのROM49,第二記憶装置としての
RAM50および入出力装置としてのI/Oポート51
を備えており、ROM49に図18のフローチャートで
表される推定値補正ルーチンを始めとする種々の制御プ
ログラムが格納されることによって、図1に示す外乱オ
ブザーバ52,相関演算部54,正規化部56,ばね定
数変化率推定部58,推定値補正部60および異常判定
部62を構成している。
As shown in FIG. 2, the computer 20 has a CPU 40 as a processing device and a ROM as a first storage device.
1 and a RAM 44 serving as a second storage device. The ROM 42 stores various control programs including a rotational speed calculation / correction routine shown in the flowchart of FIG. The side rotation speed calculation / correction unit 45 is configured. This computer 20 is connected to another computer 47. The computer 47 includes a CPU 48 as a processing device, a ROM 49 as a first storage device, a RAM 50 as a second storage device, and an I / O port 51 as an input / output device.
The ROM 49 stores various control programs such as an estimated value correction routine shown in the flowchart of FIG. 18 so that the disturbance observer 52, the correlation calculator 54, and the normalizer shown in FIG. 56, a spring constant change rate estimating unit 58, an estimated value correcting unit 60, and an abnormality determining unit 62.

【0028】コンピュータ47のI/Oポート51には
異常判定部62の判定結果を運転者に知らせる表示装置
66が接続されている。表示装置66は本実施例におい
ては液晶ディスプレイであるが、点灯あるいは点滅する
ランプ等別の表示装置を用いることも可能であり、音声
で運転者に知らせる音声報知装置などを含めて種々の形
態の報知装置を採用することが可能である。
A display device 66 is connected to the I / O port 51 of the computer 47 to notify the driver of the result of the determination made by the abnormality determining section 62. The display device 66 is a liquid crystal display in the present embodiment, but it is also possible to use another display device such as a lighted or flashing lamp, and various forms including a voice notification device that notifies the driver by voice. It is possible to employ a notification device.

【0029】コンピュータ47のI/Oポート51には
さらに、ホイール24(リム側部28)に加えられる駆
動・制動トルクを、ホイール24の軸に取り付けられた
歪みゲージ等により検出する駆動・制動トルク検出装置
68が接続されている。
The I / O port 51 of the computer 47 is further provided with a driving / braking torque for detecting a driving / braking torque applied to the wheel 24 (rim side portion 28) by a strain gauge or the like attached to a shaft of the wheel 24. The detection device 68 is connected.

【0030】I/Oポート51にはさらにまた、車体速
度センサ70と路面凹凸センサ72とがそれぞれ接続さ
れている。車体速度センサ70は例えば、波のドップラ
効果を利用して対地車速を検出するドップラ方式とした
り、4個の回転速度vに基づいて車体速度を推定する回
転速度利用方式とすることができる。路面凹凸センサ7
2は例えば、路面凹凸度Rの影響を受けた信号を読み込
んで、それのP−P値を路面凹凸度Rとして検出した
り、回転速度信号の周波数特性を路面凹凸度Rとして検
出するものとすることができる。路面凹凸度Rの影響を
受けた信号には例えば、回転速度信号,車両ばね上部の
上下加速度信号,車両ばね下部の上下加速度信号等、路
面凹凸度Rを間接に表す信号を選んだり、路面から反射
した信号等、路面凹凸度Rを直接に表す信号を選ぶこと
ができる。
The I / O port 51 is further connected to a vehicle speed sensor 70 and a road surface unevenness sensor 72, respectively. The vehicle speed sensor 70 may be, for example, a Doppler system that detects the ground vehicle speed using the Doppler effect of waves or a rotation speed utilization system that estimates the vehicle speed based on four rotation speeds v. Road surface unevenness sensor 7
2, for example, reads a signal affected by the road surface unevenness R and detects the PP value thereof as the road surface unevenness R, or detects the frequency characteristic of the rotation speed signal as the road surface unevenness R. can do. As a signal affected by the road surface unevenness R, for example, a signal representing the road surface unevenness R indirectly, such as a rotation speed signal, a vertical acceleration signal of a vehicle spring upper portion, a vertical acceleration signal of a vehicle lower spring portion, or the like, is selected. It is possible to select a signal that directly indicates the road surface unevenness R, such as a reflected signal.

【0031】リム側部回転速度演算・補正部45は上記
4個の車輪14に対応する各電磁ピックアップ12およ
び波形整形器18から供給される矩形波信号(以下、パ
ルス信号ともいう)に基づいて各車輪14の回転速度を
算出するとともに、後述の信号遅延・重ね合わせ処理に
よって各車輪14の回転速度を補正する。各車輪14お
よびロータ10には製造,組立誤差が存在し、これら誤
差等に起因して周期的な回転速度誤差が発生するため、
この各車輪14に固有の固有回転速度誤差を除いた回転
速度を求めるのである。図6に、車体速度Vが120
〔km/h〕、タイヤ半径が約330〔mm〕、タイヤ空気圧
が0.2〔MPa〕である走行状況下において取得され
た回転速度信号の周波数特性(周波数とパワースペクト
ルとの関係)の一例をグラフで表す。グラフから明らか
なように、基本周波数を約16〔Hz〕とする複数の高調
波が発生しており、これが車輪14の固有回転速度誤差
すなわち周期的外乱であり、これを除去するために回転
速度補正が行われるのである。すなわち、本実施例にお
いては、車輪14の固有回転速度誤差が除去すべき周期
的外乱なのである。
The rim side rotational speed calculating / correcting unit 45 is based on rectangular wave signals (hereinafter also referred to as pulse signals) supplied from the electromagnetic pickups 12 and the waveform shaper 18 corresponding to the four wheels 14. The rotation speed of each wheel 14 is calculated, and the rotation speed of each wheel 14 is corrected by a signal delay / superposition process described later. Manufacturing and assembly errors are present in each wheel 14 and rotor 10, and periodic errors in rotational speed occur due to these errors and the like.
The rotation speed excluding the inherent rotation speed error unique to each wheel 14 is obtained. FIG. 6 shows that the vehicle speed V is 120
[Km / h], an example of frequency characteristics (relation between frequency and power spectrum) of a rotation speed signal acquired under a running condition in which a tire radius is about 330 [mm] and a tire air pressure is 0.2 [MPa]. Is represented by a graph. As is clear from the graph, a plurality of harmonics having a fundamental frequency of about 16 [Hz] are generated, which is a natural rotational speed error of the wheel 14, that is, a periodic disturbance. The correction is made. That is, in the present embodiment, the inherent rotational speed error of the wheel 14 is a periodic disturbance to be removed.

【0032】なお、車輪14の回転速度は周速度で演算
されるが、そのためにはタイヤ26の実質的な半径(タ
イヤが荷重で変形した状態における路面から車輪14の
中心までの距離)が必要であり、これはタイヤ26の空
気圧によって変わる。よって、当初は空気圧が正規であ
る場合の正規の半径が使用されるが、後に説明する処理
によってタイヤ26の空気圧変化が判明した場合は、予
めROM42に格納されているタイヤ径テーブルからそ
の空気圧変化に対応したタイヤ半径が読み出されて使用
される。なお、回転速度が角速度で演算される場合には
タイヤ径テーブルは不要である。
The rotational speed of the wheel 14 is calculated from the peripheral speed. For this, a substantial radius of the tire 26 (the distance from the road surface to the center of the wheel 14 when the tire is deformed by the load) is required. Which depends on the air pressure of the tire 26. Therefore, at first, the normal radius when the air pressure is normal is used. However, if the air pressure change of the tire 26 is found by the processing described later, the air pressure change is obtained from the tire diameter table stored in the ROM 42 in advance. Is read and used. When the rotational speed is calculated based on the angular velocity, the tire diameter table is not required.

【0033】リム側部回転速度演算・補正部45の機能
は図5の回転速度演算・補正ルーチンの実行により果た
されるが、まず、その原理を説明する。回転速度に侵入
する外乱は、車輪14と共に1回転する周期性を有す
る。したがって、この周期性を利用し、回転速度の時間
的変化のうち先行する部分をちょうど1周期分だけ遅延
させて後続する部分に重ね合わせることにより、周期的
外乱を除去することにする。信号遅延・重ね合わせ処理
の伝達関数は図7のブロック線図で表すことができ、こ
の場合の状態方程式は、z変数を用いれば、
The function of the rim side rotation speed calculation / correction section 45 is fulfilled by executing the rotation speed calculation / correction routine shown in FIG. 5. First, the principle will be described. The disturbance entering the rotation speed has a periodicity of making one rotation with the wheel 14. Therefore, by utilizing this periodicity, the periodic disturbance is removed by delaying the preceding portion of the temporal change in the rotational speed by exactly one cycle and superimposing it on the subsequent portion. The transfer function of the signal delay / superposition process can be represented by the block diagram of FIG. 7, and the state equation in this case is obtained by using the z variable.

【0034】[0034]

【数1】 (Equation 1)

【0035】となる。ただし、 V :入力としての車体速度 v :入力としての回転速度 y :出力としての回転速度 L :回転速度vのサンプリング周期Tを単位として記
述される遅延量
## EQU1 ## Here, V: body speed as input v: rotation speed as input y: rotation speed as output L: delay amount described in units of sampling period T of rotation speed v

【0036】なお、図において遅延演算部は、旧誤差値
2 を遅延量L(L回分のサンプリング周期T)だけ遅
延させて入力としての誤差値e(=V−y)に加算して
新誤差値e1 を得ることにより、入力vにおいてそれの
変動成分(すなわち除去すべき周期的外乱)を相対的に
強調して変動成分を抽出する機能を果たす。
In the figure, the delay calculating section delays the old error value e 2 by the delay amount L (L sampling cycles T), adds it to the input error value e (= V−y), and By obtaining the error value e 1 , the function of extracting the variable component is obtained by relatively emphasizing the variable component (that is, the periodic disturbance to be removed) at the input v.

【0037】入力vから出力yまでの伝達関数H
1 (z)は、
Transfer function H from input v to output y
1 (z) is

【0038】[0038]

【数2】 (Equation 2)

【0039】となる。ここで、伝達関数H1 (z)の周
波数応答H1 (ω)を求めるために、 z-1=e-jwn とおく。ただし、 j :虚数単位 ωn :正規化角周波数 なお、正規化角周波数ωn は、1秒間を単位として記述
される実角周波数ωAではなく、サンプリング周期Tを
単位として記述される角周波数である。そのため、正規
化角周波数ωn は、 ωn =ωA T で表される。
## EQU1 ## Here, in order to obtain the frequency response H 1 (ω) of the transfer function H 1 (z), z −1 = e −jwn is set. Here, j: imaginary unit ω n : normalized angular frequency Note that the normalized angular frequency ω n is not the actual angular frequency ω A described in units of 1 second, but the angular frequency described in units of the sampling period T. It is. Therefore, the normalized angular frequency omega n is expressed by ω n = ω A T.

【0040】その結果、周波数応答H1 (ω)が、As a result, the frequency response H 1 (ω) becomes

【0041】[0041]

【数3】 (Equation 3)

【0042】として求められる。この周波数応答H
1 (ω)の絶対値の二乗は、
Is obtained as follows. This frequency response H
The square of the absolute value of 1 (ω) is

【0043】[0043]

【数4】 (Equation 4)

【0044】となる。この式から明らかなように、cos
ωn L=1すなわちωn L=2nπ(n:整数)のと
き、周波数応答H1 (ω)が0、すなわち入力vが出力
yに全く現れない状態となる。また、周波数応答H
1 (ω)が0となる複数の角周波数ω n のうち互いに隣
接する2個のものの中間値(ωn L≠2nπである各区
間における代表点)において、入力vがそのまま出力y
に現れるようにするためには、例えば、ωn L=(2n
+1)πのときに周波数応答H1 (ω)の絶対値が1と
ならなければならないが、入力vをそのまま使用したの
では、
Is as follows. As is clear from this equation, cos
ωnL = 1, that is, ωnL = 2nπ (n: integer)
Frequency response H1(Ω) is 0, that is, input v is output
It is in a state that does not appear at all in y. Also, the frequency response H
1Multiple angular frequencies ω where (ω) is 0 nNext to each other
The median value of two touching objects (ωnEach ward where L ≠ 2nπ
Input point v at the output y
, For example, ωnL = (2n
+1) Frequency response H at π1The absolute value of (ω) is 1
Must use the input v as it is
Then

【0045】[0045]

【数5】 (Equation 5)

【0046】となってしまう。そこで、入力vをそのま
ま使用するのではなく、それに|1−g/2|なる調整
ゲインを掛け算した値v1 として使用することとする。
## EQU1 ## Therefore, instead of having to utilize the input v, it | 1-g / 2 | and the use of consisting adjustment gain as a value v 1 obtained by multiplication.

【0047】以上、入力vから出力yまでの伝達関数H
1 (z)について説明したが、出力yには入力vの影響
のみならず入力Vの影響も加えられる。そこで、出力y
における入力Vの影響を評価するため、入力Vから出力
yまでの伝達関数H2 (z)をも求めることとする。こ
の伝達関数H2 (z)は伝達関数H1 (z)に同様にし
て求められ、
As described above, the transfer function H from the input v to the output y is
Although 1 (z) has been described, the output y is affected not only by the input v but also by the input V. Then, output y
In order to evaluate the effect of the input V in the above, the transfer function H 2 (z) from the input V to the output y is also determined. This transfer function H 2 (z) is obtained in the same manner as the transfer function H 1 (z),

【0048】[0048]

【数6】 (Equation 6)

【0049】となる。これの周波数応答H2 (ω)は、Is as follows. The frequency response H 2 (ω) of this is

【0050】[0050]

【数7】 (Equation 7)

【0051】となり、この周波数応答H2 (ω)の絶対
値の二乗は、
The square of the absolute value of the frequency response H 2 (ω) is

【0052】[0052]

【数8】 (Equation 8)

【0053】となる。この式から明らかなように、ωn
L=2nπのとき、周波数応答H2 (ω)の絶対値は1
となり、その結果、 V=y となる。すなわち、出力yは、 y=H1 (z)v1 +H2 (z)V なる式で近似的に表すことができ、ωn L≠2nπのと
きには、 y=v1 となり、ωn L=2nπのときには、 y=V となるのである。ここで、入力v1 を、時間的変動の中
心である固定成分(オフセット値ともいう)v0 と変動
成分Δvとに分離して考えると、オフセットv0は本来
入力Vに一致すべきであるから、ωn L≠2nπのとき
には、 y=v0 +Δv1 となり、ωn L=2nπのときには、 y=v0 となる。
Is as follows. As is clear from this equation, ω n
When L = 2nπ, the absolute value of the frequency response H 2 (ω) is 1
And consequently V = y. That is, the output y can be approximately expressed by the following equation: y = H 1 (z) v 1 + H 2 (z) V. When ω n L ≠ 2nπ, y = v 1 and ω n L = When 2nπ, y = V. Here, when the input v 1 is considered as being separated into a fixed component (also referred to as an offset value) v 0 which is the center of the temporal variation and a variation component Δv, the offset v 0 should originally match the input V. Therefore, when ω n L ≠ 2nπ, y = v 0 + Δv 1 , and when ω n L = 2nπ, y = v 0 .

【0054】したがって、固定成分v0 を0とみなし、
入力v1 として変動成分Δv、入力Vとして0をそれぞ
れ使用すれば、出力yがΔvとなり、これが求めるべき
周期的外乱の時間的変化を表す物理量となる。また、
に〔0073〕欄において記載されているように、周期
的外乱を除去された出力yすなわち変動成分Δvに固定
成分v0 すなわちオフセット値を加算すれば、周期的外
乱が除去された回転速度vが取得されることになる。
Therefore, assuming that the fixed component v 0 is 0,
If the fluctuation component Δv is used as the input v 1 and 0 is used as the input V, the output y becomes Δv, which is a physical quantity representing the temporal change of the periodic disturbance to be obtained. Also after
As described in [0073], the period
If the fixed component v 0, ie, the offset value, is added to the output y, ie, the fluctuation component Δv, from which the target disturbance has been removed, the rotation speed v, from which the periodic disturbance has been removed, is obtained.

【0055】前記のように、遅延量Lについては、 ωn L=2nπ なる関係が成立するため、遅延量Lは、 L=(2nπ)/ωn で表される。ここで、n=1とすれば、 L=(2π)/ωn となり、また、前記のように、 ωn =ωA T なる関係が成立する。また、除去すべき周期的外乱の基
本波の周波数である基本周波数f1 は、 1/f1 =(2π)/ωA で表される。ただし、f1 =V/(2πR)である。し
たがって、遅延量Lは、 L=(1/f1 )/T で表される。
As described above, since the relationship of ω n L = 2nπ holds for the delay amount L, the delay amount L is represented by L = (2nπ) / ω n . Here, if n = 1, then L = (2π) / ω n , and, as described above, the relationship of ω n = ω A T is established. The fundamental frequency f 1, which is the frequency of the fundamental wave of the periodic disturbance to be removed, is represented by 1 / f 1 = (2π) / ω A. Here, f 1 = V / (2πR). Therefore, the delay amount L is represented by L = (1 / f 1 ) / T.

【0056】この式は物理的には、『遅延量Lが、除去
すべき周期的外乱が1周期進行する間に、すなわち、車
輪14が1回転する間に回転速度vのサンプリングが行
われる回数に等しい』ことを表している。
This equation is physically expressed as "the delay amount L is the number of times the rotation speed v is sampled while the periodic disturbance to be removed advances one cycle, that is, while the wheel 14 makes one rotation. Is equal to ”.

【0057】タイヤ半径Rが約330〔mm〕、車体速度
Vが120〔km/h〕である場合には、基本周波数f1
約16〔Hz〕となり、また、サンプリング周期Tを5
〔ms〕とすれば、この場合の遅延量Lは約12〔サンプ
リング周期〕となる。図8に、L=12,g=0.4と
した場合の伝達関数H1 (z)の周波数特性をグラフで
表す。このグラフは縦軸に周波数応答H1 (ω)がdB
単位でとられているため、周波数応答H1 (ω)が1で
あることがグラフでは0で表され、周波数応答H
1 (ω)が0であることがグラフでは−∞で表されるこ
とになる。したがって、グラフから明らかなように、1
6〔Hz〕を基本周波数f1 とする複数の高調波が発生す
る周波数の各位置において周波数応答H1 (ω)が十分
に0に近くなり、それ以外の各位置において周波数応答
1 (ω)が十分に1に近くなっている。
The tire radius R is about 330 [mm], and the vehicle speed is
When V is 120 [km / h], the fundamental frequency f1But
Approximately 16 [Hz], and the sampling period T is 5
[Ms], the delay amount L in this case is about 12 [sample
Ring cycle]. FIG. 8 shows that L = 12 and g = 0.4.
Transfer function H1Graph of frequency characteristics of (z)
Represent. In this graph, the vertical axis represents the frequency response H1(Ω) is dB
Frequency response H1(Ω) is 1
That is represented by 0 in the graph, and the frequency response H
1The fact that (ω) is 0 is represented by -∞ in the graph.
And Therefore, as is clear from the graph, 1
6 [Hz] is the fundamental frequency f1Multiple harmonics occur
Frequency response H at each frequency position1(Ω) is enough
And the frequency response at each other position
H1(Ω) is sufficiently close to 1.

【0058】次に、この原理に従って回転速度演算・補
正技術をコンピュータを用いて実施するための具体的技
術について説明する。
Next, a description will be given of a specific technique for implementing a rotational speed calculation / correction technique using a computer according to this principle.

【0059】まず、概略的に説明する。電磁ピックアッ
プ12から各パルス(一連の周期的信号である矩形波を
構成する各単位波)が供給される毎に、各パルスの時間
間隔に基づいて元回転速度vが順に演算される。各パル
スの時間間隔には例えば、各パルスの立ち上がり間隔時
間,立ち下がり間隔時間,パルス中点間隔時間等を選ぶ
ことができる。
First, a brief description will be given. Each time each pulse (each unit wave constituting a series of rectangular signals, which is a rectangular wave) is supplied from the electromagnetic pickup 12, the original rotational speed v is sequentially calculated based on the time interval between each pulse. As the time interval of each pulse, for example, a rising interval time, a falling interval time, and a pulse midpoint interval time of each pulse can be selected.

【0060】また、1サンプリング周期Tが経過する毎
に、その間に入力・演算された少なくとも1個の元回転
速度vの平均値がサンプル回転速度vとして演算され
る。さらに、過去に演算されたサンプル回転速度vと現
サンプル回転速度vとからその平均値が演算される。こ
の平均値は前記オフセット値として使用され、現サンプ
ル回転速度vからそのオフセット値が減算されることに
より、サンプル変動回転速度Δv(図7における「v」
に相当する)が求められる。
Every time one sampling period T elapses, an average value of at least one original rotation speed v input and calculated during that period is calculated as a sample rotation speed v. Further, the average value is calculated from the sample rotation speed v calculated in the past and the current sample rotation speed v. This average value is used as the offset value, and the offset value is subtracted from the current sample rotation speed v to obtain the sample variation rotation speed Δv (“v” in FIG. 7).
Is required).

【0061】また、RAM44には1回転分メモリが設
けられている。この1回転分メモリは誤差値E(図7に
おける「e(=V−y)」に相当する)をロータ10の
各歯に関連付けて記憶するものである。誤差値Eは、除
去すべき周期的外乱、すなわち、本実施例においては車
輪の不均一性に起因する外乱がロータ10の回転と共に
変化する状態を表す値であり、この取得方法については
後に詳述する。1回転分メモリはロータ10が1回転す
る間に電磁ピックアップ12から出力されるパルスの
数、すなわちロータ10の歯の数と同数の記憶アドレス
が設けられており、各記憶アドレスに各誤差値Eが順に
記憶されることとなる。
The RAM 44 is provided with a memory for one rotation. The memory for one rotation stores an error value E (corresponding to “e (= V−y)” in FIG. 7) in association with each tooth of the rotor 10. The error value E is a value representing a periodic disturbance to be removed, that is, a state in which the disturbance caused by the non-uniformity of the wheels changes with the rotation of the rotor 10 in this embodiment. Will be described. The memory for one rotation is provided with the same number of storage addresses as the number of pulses output from the electromagnetic pickup 12 during one rotation of the rotor 10, that is, the number of teeth of the rotor 10, and each error value E is assigned to each storage address. Are sequentially stored.

【0062】図9には、ロータ10の歯数が8である場
合を例にとり、1回転分メモリの構成が概念的に示され
ている。図において、「1−1」は第1回転時において
第1番目に演算された誤差値Eを示し、「1−2」は第
1回転時において第2番目に演算された誤差値Eを示
し、「2−1」は第2回転時において第1番目に演算さ
れた誤差値Eを示すというように、ハイフンで連結され
た2個の数字のうち左側の数字は車輪14が第何回転目
にあるかを示し、右側の数字は各パルスが各回転時にお
いて第何番目にあるかを示している。
FIG. 9 conceptually shows the structure of the memory for one rotation, taking as an example the case where the rotor 10 has eight teeth. In the figure, “1-1” indicates an error value E calculated first during the first rotation, and “1-2” indicates an error value E calculated second during the first rotation. , "2-1" indicates the error value E calculated first during the second rotation, and the left-hand number of the two numbers connected by a hyphen indicates the number of rotations of the wheel 14 , And the number on the right side indicates the number of each pulse in each rotation.

【0063】サンプリング周期Tが経過する毎に、この
1回転分メモリから誤差値Eが読み込まれる。1回転分
メモリに既に記憶されている少なくとも1個の誤差値E
のうち、今回のサンプル変動回転速度Δvに係る少なく
とも1個の変動回転速度Δvの各々とロータ10のパル
ス発生回転位置が共通する誤差値群Eが読み込まれる。
すなわち、今回のサンプル変動回転速度Δvよりちょう
どロータ1回転前に取得された誤差値群Eが旧誤差値群
Eとして読み出されるのである。読み出された旧誤差値
群Eについては平均値が求められ、平均旧誤差値EMEAN
(図7における遅延演算部の「e2 」に相当する)とさ
れる。
Every time the sampling period T elapses, the error value E is read from the memory for one rotation. At least one error value E already stored in the memory for one rotation
Among them, the error value group E in which each of the at least one variable rotational speed Δv relating to the current sample variable rotational speed Δv and the pulse generation rotational position of the rotor 10 are common is read.
That is, the error value group E obtained just before the rotation of the rotor by one rotation from the current sample fluctuation rotation speed Δv is read as the old error value group E. An average value is obtained for the read old error value group E, and the average old error value E MEAN is obtained.
(Corresponding to “e 2 ” of the delay calculation unit in FIG. 7).

【0064】例えば、図9の例を用いて説明すれば、各
回のサンプリングにおいて変動回転速度Δvが2個ずつ
得られると仮定すれば、図10に示すように、例えば、
図9における「1−1」と「1−2」とが「1−」、
すなわち、第1回転時における第1番目の旧誤差値群E
を構成し、「1−1」と「1−2」でそれぞれ示される
2個の旧誤差値Eの平均値が「1−」で示される平均
旧誤差値EMEANとされる。また、各回のサンプリングに
おいて変動回転速度Δvが4個ずつ得られると仮定すれ
ば、図11に示すように、例えば、図9における「1−
1」〜「1−4」が「1−」、すなわち、第1回転時
における第1番目の旧誤差値群Eを構成し、「1−1」
〜「1−4」でそれぞれ示される4個の旧誤差値Eの平
均値が「1−」で示される平均旧誤差値EMEANとされ
る。
For example, referring to the example of FIG. 9, assuming that two variable rotation speeds Δv can be obtained in each sampling, as shown in FIG.
“1-1” and “1-2” in FIG. 9 are “1-”,
That is, the first old error value group E during the first rotation
And the average value of the two old error values E respectively indicated by “1-1” and “1-2” is set as the average old error value E MEAN indicated by “1-”. Further, assuming that four variable rotation speeds Δv are obtained in each sampling, for example, as shown in FIG.
"1" to "1-4" constitute "1-", that is, a first old error value group E at the time of the first rotation, and "1-1".
The average value of the four old error values E respectively indicated by .about. "1-4" is defined as the average old error value E.sub.MEAN indicated by "1-".

【0065】このようにして求められた今回の平均旧誤
差値EMEANと前記ゲインgとの積が求められ、さらに、
今回のサンプル変動回転速度Δvと前記調整ゲイン|1
−g/2|との積(図7における「v1 」に相当する)
も求められる。さらに、それら2個の積の和が演算さ
れ、図7における出力yとされる。すなわち、出力y
が、 旧EMEAN・g+サンプルΔv・|1−g/2| として求められるのである。この出力yに前記オフセッ
ト値が加算されることにより、最終的な回転速度vが求
められることになる。
The product of the current average old error value E MEAN thus obtained and the gain g is obtained.
The current sample fluctuation rotation speed Δv and the adjustment gain | 1
Product with −g / 2 | (corresponding to “v 1 ” in FIG. 7)
Is also required. Further, the sum of the two products is calculated, and the result is taken as the output y in FIG. That is, the output y
Is obtained as old E MEAN · g + sample Δv · | 1-g / 2 |. By adding the offset value to the output y, the final rotational speed v is obtained.

【0066】今回の出力yと今回の平均旧誤差値EMEAN
とから新誤差値E(図7における遅延演算部の「e1
に相当する)が求められる。すなわち、新誤差値Eが、 e+e2 =(V−y)+e2 =(0−y)+e2 =−y
+旧EMEAN として求められるのである。
The current output y and the current average old error value E MEAN
From the new error value E (“e 1 ” of the delay calculation unit in FIG. 7)
Is required). That is, the new error value E is e + e 2 = (V−y) + e 2 = (0−y) + e 2 = −y
+ It is required as an old E MEAN .

【0067】新誤差値Eが求められたならば、1回転分
メモリにおいて誤差値Eの更新が行われる。すなわち、
今回の旧誤差値群Eを構成する各旧誤差値Eの各々がい
ずれも同じ新誤差値Eに書き換えられるのである。
When the new error value E is obtained, the error value E is updated in the memory for one rotation. That is,
Each of the old error values E constituting the current old error value group E is rewritten to the same new error value E.

【0068】図9の例を用いて説明すれば、例えば、
「2−」に対応する新誤差値Eが演算された場合に
は、「1−」で示される旧誤差値群Eを構成する「1
−1」と「1−2」でそれぞれ示される旧誤差値Eがい
ずれも同じ新誤差値Eに書き換えられることになる。
With reference to the example of FIG. 9, for example,
When the new error value E corresponding to “2-” is calculated, “1” constituting the old error value group E indicated by “1-”
The old error values E indicated by “−1” and “1-2” are both rewritten to the same new error value E.

【0069】ただし、上記の処理が何回も繰り返される
ことによって1回転分メモリに記憶されている誤差値E
が除去すべき周期的外乱を正確に反映するに至れば、出
力yが0となり、ひいては新誤差値Eがすわなち旧誤差
値Eとなるから、この段階においては誤差値Eの書換え
が行われても誤差値Eは実質的に固定されることにな
る。
However, the error value E stored in the memory for one rotation is obtained by repeating the above processing many times.
, Accurately reflects the periodic disturbance to be removed, the output y becomes 0, and the new error value E becomes the old error value E. Therefore, the error value E is rewritten at this stage. In this case, the error value E is substantially fixed.

【0070】ここで、1回転分メモリと遅延量Lとの関
係について説明する。本実施例においては、各パルスが
発生する毎に元回転速度vの演算が行われるが、サンプ
ル変動回転速度Δvの演算,平均誤差値EMEANの演算,
出力yの演算および最終的回転速度vの演算は各パルス
毎にではなく、サンプリング周期Tが経過する毎に行わ
れる。したがって、1回転分メモリは実質的には、車輪
14が1回転する間にサンプリングされるパルスの数と
同数の誤差値E(図10の例では4個、図11の例では
2個)が記憶されるメモリであるということができる。
Here, the relationship between the memory for one rotation and the delay amount L will be described. In the present embodiment, the calculation of the original rotation speed v is performed every time each pulse is generated. However, the calculation of the sample fluctuation rotation speed Δv, the calculation of the average error value E MEAN ,
The calculation of the output y and the calculation of the final rotational speed v are performed not every pulse but every time the sampling period T elapses. Therefore, the memory for one rotation has substantially the same number of error values E (four in the example of FIG. 10 and two in the example of FIG. 11) as the number of pulses sampled during one rotation of the wheel 14. It can be said that it is a memory to be stored.

【0071】また、本実施例においては、先に演算され
た元回転速度vをロータ1回転分遅延させて後に演算さ
れる元回転速度vに重ね合わせる際のその遅延量が、1
回転分メモリに記憶されている誤差値群Eの数(すなわ
ち、サンプル変動回転速度Δvの数)と等しくされる。
In this embodiment, when the previously calculated original rotational speed v is delayed by one rotation of the rotor and then superimposed on the original calculated rotational speed v, the delay amount is 1
The number is equal to the number of error value groups E stored in the rotation memory (that is, the number of sample fluctuation rotation speeds Δv).

【0072】この遅延量Lは前記の理論の説明において
は、サンプリング周期Tを単位として、 L=2πR/V/T で表され、また、サンプリング周期Tは不変である。し
たがって、周期的外乱をその周期性を利用して除去する
ためには、遅延量Lを車体速度Vすなわち回転速度vに
応じて変化させることが必要となる。これに対し、本実
施例においては、前記のように、1回転分メモリの数が
ロータ10の歯数と同数とされ、常に最新かつロータ1
回転分の誤差値Eが記憶されるようにされる。しかも、
サンプリング周期Tが経過する毎に、その間に演算され
た少なくとも1個の誤差値Eについて平均値が求められ
て唯一の代表値として用いられることにより、1回転分
メモリの記憶アドレスが実質的にはサンプリング周期T
毎に得られた元回転速度vの数に応じて分割され、結
局、その分割数が遅延量Lと等しくなっている。したが
って、本実施例においては、回転速度v毎に遅延量Lの
異なる遅延演算部を設けることなく、1個の遅延演算部
で遅延量Lが回転速度vに応じて自動的に変化させられ
ることになる。
In the description of the above theory, the amount of delay L is represented by L = 2πR / V / T using the sampling period T as a unit, and the sampling period T is unchanged. Therefore, in order to remove the periodic disturbance using its periodicity, it is necessary to change the delay amount L according to the vehicle speed V, that is, the rotation speed v. On the other hand, in the present embodiment, as described above, the number of memories for one rotation is the same as the number of teeth of the rotor 10, and the memory is always the latest
The error value E for the rotation is stored. Moreover,
Each time the sampling period T elapses, the average value of at least one error value E calculated during the sampling period T is obtained and used as the only representative value. Sampling period T
The division is performed in accordance with the number of the original rotational speeds v obtained each time. As a result, the division number is equal to the delay amount L. Therefore, in the present embodiment, the delay amount L can be automatically changed according to the rotation speed v by one delay calculation unit without providing a delay operation unit having a different delay amount L for each rotation speed v. become.

【0073】以上のようにして周期的外乱が除去された
出力yが演算されたならば、これに前記オフセット値が
加算されることにより、周期的外乱が除去された最終的
な回転速度vが演算される。
After the output y from which the periodic disturbance has been removed is calculated as described above, the final rotational speed v from which the periodic disturbance has been removed can be obtained by adding the offset value to the output y. Is calculated.

【0074】なお、図12には変動回転速度Δvの時間
的変化の二例がそれぞれグラフで表されている。上側の
グラフは回転速度vが遅いときの時間的変化の一例を示
し、下側のグラフは回転速度vが速いときの時間的変化
の一例を示している。それらグラフを対比すれば、回転
速度vが速いほど1回のサンプリング中に1回転分メモ
リに記憶される変動回転速度Δvの数(図において〜
で表す)が増加することが分かり、さらに、回転速度
vが速いほど1回の遅延量Lに対応する変動回転速度群
Δvの数が減少することも分かる。
FIG. 12 is a graph showing two examples of the temporal change of the variable rotation speed Δv. The upper graph shows an example of a temporal change when the rotation speed v is slow, and the lower graph shows an example of a temporal change when the rotation speed v is fast. Comparing these graphs, as the rotation speed v increases, the number of variable rotation speeds Δv stored in the memory for one rotation during one sampling (in FIG.
It is also found that the number of the variable rotation speed groups Δv corresponding to one delay amount L decreases as the rotation speed v increases.

【0075】以上、回転速度演算・補正技術の内容を概
略的に説明したが、次に、図5の回転速度演算・補正ル
ーチンに基づいて具体的に説明する。
The contents of the rotation speed calculation / correction technique have been described briefly above. Next, the rotation speed calculation / correction routine of FIG. 5 will be specifically described.

【0076】まず、ステップS51(以下、単にS51
で表す。他のステップについても同じ)において、整数
iの値が1に初期化され、次に、S52において、今回
のサンプリング周期が開始されるのを待つ状態となる。
今回のサンプリング周期が開始されたならば、S53に
おいて、今回のサンプリング周期が終了したか否かが判
定される。今回は開始されたばかりであるから、判定が
NOとなり、S54において、電磁ピックアップ12か
ら第i番目のパルスが出されるのを待つ状態となる。出
されたならば判定がYESとなり、S55において、そ
のパルスの継続時間が計測され、S56において、その
継続時間に基づいて第i番目の元回転速度v(i) が演算
される。すなわち、コンピュータ20のこのS56を実
行する部分によって「元回転速度演算部」が構成されて
いるのである。演算値はRAM44の1サンプリング分
メモリに記憶される。
First, in step S51 (hereinafter simply referred to as S51
Expressed by In the other steps as well, the value of the integer i is initialized to 1, and then, in S52, a state of waiting for the start of the current sampling cycle is set.
If the current sampling cycle has started, it is determined in S53 whether or not the current sampling cycle has ended. Since this time has just been started, the determination is NO, and the flow waits for the electromagnetic pickup 12 to output the i-th pulse in S54. If so, the determination is YES, the duration of the pulse is measured in S55, and the i-th original rotational speed v (i) is calculated based on the duration in S56. That is, the part that executes S56 of the computer 20 constitutes the “original rotational speed calculation unit”. The calculated value is stored in the memory for one sampling of the RAM 44.

【0077】その後、S57において、整数iの値が1
増加させられ、S58において、整数iの現在値が最大
値iMAX を超えたか否かが判定される。最大値iMAX
ロータ10の歯数と同じ値に設定されているから、この
判定は結局、1回転分の元回転速度v(i) が取得される
ごとにYESとなることになる。今回は整数iの現在値
が最大値iMAX を超えてはいないと仮定すれば判定がN
Oとなり、直ちにS53に戻るが、超えていると仮定す
れば判定がYESとなり、S66において整数iの値が
1にリセットされた後にS53に戻ることになる。
Thereafter, in S57, the value of the integer i is 1
It is increased, in S58, the current value of the integer i whether exceeds the maximum value i MAX is determined. Since the maximum value iMAX is set to the same value as the number of teeth of the rotor 10, this determination is eventually YES every time the original rotation speed v (i) for one rotation is obtained. This time, if it is assumed that the current value of the integer i does not exceed the maximum value i MAX , the determination is N
The result is O, and the process immediately returns to S53. However, if the value is exceeded, the determination becomes YES, and the process returns to S53 after the value of the integer i is reset to 1 in S66.

【0078】その後S53〜S58およびS66の実行
が繰り返されるうちに今回のサンプリング周期が終了す
れば、S53の判定がYESとなり、S59以下のステ
ップに移行する。
Thereafter, if the current sampling period ends while the execution of S53 to S58 and S66 is repeated, the determination in S53 becomes YES, and the process proceeds to S59 and subsequent steps.

【0079】S59においては、今回のサンプリング周
期の間に取得されて1サンプリング分メモリに記憶され
ている元回転速度vの平均値が演算され、サンプル回転
速度vとされる。すなわち、コンピュータ20のこのS
59を実行する部分によって「サンプル回転速度演算
部」が構成されているのである。
In S59, the average value of the original rotational speed v obtained during the current sampling period and stored in the memory for one sampling is calculated, and is set as the sample rotational speed v. That is, this S of the computer 20
That is, the “sample rotation speed calculation unit” is constituted by the part that executes step 59.

【0080】S60において、過去に得られたサンプル
回転速度vと現サンプル回転速度vとからその平均値が
演算され、これが前記オフセット値とされる。さらに、
現サンプル回転速度vからこの平均値が減算されること
により、今回のサンプル変動回転速度がΔvとして演算
される。すなわち、コンピュータ20のこのS60を実
行する部分によって「サンプル変動回転速度演算部」が
構成されているのである。
In S60, the average value is calculated from the sample rotation speed v obtained in the past and the current sample rotation speed v, and this is set as the offset value. further,
By subtracting this average value from the current sample rotation speed v, the current sample rotation speed is calculated as Δv. That is, the portion that executes S60 of the computer 20 constitutes a “sample variation rotation speed calculation unit”.

【0081】なお、このサンプル変動回転速度演算部を
例えば、遮断周波数が1〔Hz〕程度のハイパスフィルタ
手段として構成し、このハイパスフィルタ手段にサンプ
ル回転速度vを供給することによってもサンプル変動回
転速度Δvを取得可能である。また、遮断周波数が1
〔Hz〕程度のローパスフィルタ手段にサンプル回転速度
vを供給することによっても前記オフセット値を取得可
能である。
It is to be noted that the sample fluctuation rotation speed calculating section may be constituted as high-pass filter means having a cutoff frequency of about 1 [Hz], for example, and the sample rotation speed v may be supplied to the high-pass filter means. Δv can be obtained. When the cutoff frequency is 1
The offset value can also be obtained by supplying the sample rotation speed v to a low-pass filter means of about [Hz].

【0082】S61においては、1回転分メモリに既に
記憶されている旧誤差値Eのうち、今回の変動回転速度
群Δvと記憶アドレスが一致するものが読み出される。
すなわち、ちょうどロータ1回転分前に取得された旧誤
差値群Eが読み出されるのである。さらに、今回の旧誤
差値群Eの平均値が平均旧誤差値EMEANとして演算され
る。さらにまた、今回のサンプル変動回転速度Δvと調
整ゲイン|1−g/2|との積と、平均旧誤差値EMEAN
とゲインgとの積との和が出力yとして演算される。す
なわち、コンピュータ20のうちこのS61を実行する
部分によって「出力値演算部」が構成されているのであ
る。
At S61, the old error value E which has been stored in the memory for one rotation and whose storage address matches the current variable rotation speed group Δv is read out.
That is, the old error value group E obtained just one rotation before the rotor is read. Further, the average value of the old error value group E this time is calculated as the average old error value E MEAN . Furthermore, the current sample variation speed Δv adjustment gain | 1-g / 2 | and the product of the average old error value E MEAN
The sum of the product of the gain and the gain g is calculated as the output y. That is, the “output value calculation unit” is configured by the part of the computer 20 that executes S61.

【0083】S62においては、その出力yに前記オフ
セット値が加算されることにより、今回の最終的な回転
速度v、すなわち、周期的外乱が除去された回転速度v
が求められる。すなわち、コンピュータ20のうちこの
S62を実行する部分によって「最終的回転速度演算
部」が構成されているのである。
In step S62, the offset value is added to the output y, whereby the current final rotational speed v, that is, the rotational speed v from which the periodic disturbance has been removed, is obtained.
Is required. That is, the “final rotation speed calculation unit” is configured by the part of the computer 20 that executes S62.

【0084】また、このS62においては、信号遅延・
重ね合わせ処理による周期的外乱除去状況が判定される
ようになっている。具体的には、最終的な回転速度vが
時間的にほぼ安定したか否かが判定され、安定した場合
には周期的外乱の除去が実質的に完了したと判定される
ようになっている。さらに、周期的外乱の除去が実質的
に完了したと判定した場合には、ON状態で周期的外乱
の除去が実質的に完了したことを表し、OFF状態で実
質的に完了していないことを表す周期的除去完了フラグ
(RAM44に設けられている)がON状態とされ、一
方、周期的外乱の除去が実質的に完了してはいないと判
定した場合には、周期的外乱除去完了フラグがOFF状
態とされるようにもなっている。なお、ここに「最終的
な回転速度vが時間的にほぼ安定したか否かの判定」は
例えば、ロータ10が1回転するごとに、ロータ10の
各歯16に対応する最終的な回転速度vの前回値から今
回値を差し引いた値の絶対値の和を回転速度vの変動量
δとして取得し、その変動量δが基準値δ0 以下である
場合には最終的な回転速度vが時間的にほぼ安定したと
判定するものとすることができる。
In S62, the signal delay and
The periodic disturbance removal state by the superimposition process is determined. Specifically, it is determined whether or not the final rotation speed v is substantially stable in time, and when it is stabilized, it is determined that the removal of the periodic disturbance is substantially completed. . Further, when it is determined that the removal of the periodic disturbance is substantially completed, it indicates that the removal of the periodic disturbance is substantially completed in the ON state, and it is determined that the removal of the periodic disturbance is not substantially completed in the OFF state. When the periodic disturbance removal completion flag (provided in the RAM 44) is turned on, and it is determined that the removal of the periodic disturbance has not been substantially completed, the periodic disturbance removal completion flag is set. It is also turned off. Here, “determination of whether or not the final rotation speed v is substantially stable in time” is performed, for example, every time the rotor 10 makes one rotation, the final rotation speed corresponding to each tooth 16 of the rotor 10 The sum of the absolute values of the values obtained by subtracting the current value from the previous value of v is obtained as the variation δ of the rotation speed v. If the variation δ is equal to or less than the reference value δ 0 , the final rotation speed v is It can be determined that the time is almost stable.

【0085】S63においては、今回の出力yが平均旧
誤差値EMEANから減算されることにより今回の新誤差値
Eが求められる。すなわち、コンピュータ20のうちこ
のS63を実行する部分によって「誤差値演算部」が構
成されているのである。
In S63, the current output y is subtracted from the average old error value E MEAN to obtain the current new error value E. That is, the “error value calculation unit” is configured by the part of the computer 20 that executes S63.

【0086】S64においては、1回転分メモリにおい
て今回の旧誤差値群Eの各値がすべて新誤差値Eに書き
換えられる。すなわち、コンピュータ20のうちこのS
64を実行する部分によって「1回転分メモリ更新部」
が構成されているのである。その後、S52に戻る。
In S64, all the values of the old error value group E at this time are rewritten to the new error value E in the memory for one rotation. That is, this S in the computer 20
64 performs "one rotation memory update unit"
Is configured. Then, the process returns to S52.

【0087】以上、回転速度演算・補正技術を詳細に説
明したが、この技術を実施した場合の効果の一例を図1
3にグラフで表す。このグラフには、この技術を実施し
た後の回転速度vの周波数特性が実線で示され、これと
比較するために、この技術を実施する前の回転速度vの
周波数特性、すなわち、図6のグラフと同じものが破線
で示されている。それらグラフから明らかなように、こ
の技術を実施すれば、車輪14に固有の周期的外乱を良
好に除去されることになる。
The rotation speed calculation / correction technique has been described in detail above. One example of the effect of implementing this technique is shown in FIG.
3 shows a graph. In this graph, the frequency characteristic of the rotational speed v after the implementation of this technique is shown by a solid line, and for comparison with this, the frequency characteristic of the rotational speed v before the implementation of this technique, that is, FIG. The same thing as a graph is shown by the broken line. As is apparent from these graphs, the implementation of this technique effectively removes the periodic disturbance inherent in the wheels 14.

【0088】次に、外乱オブザーバ52について説明す
る。外乱オブザーバ52は、車輪14の図4に示すモデ
ルに基づいて構成されている。以下、この外乱オブザー
バ52の構成について説明する。車輪14を、慣性モー
メントJR のリム側部28と慣性モーメントJB のベル
ト側部30とがばね定数Kのねじりばね32により接続
されたものとしてモデル化すれば、(2) 〜(4) の状態方
程式が成立し、これによって線形システムが構成され
る。 JR ωR ′=−KθRB+T1 ・・・(2) JB ωB ′= KθRB−Td ・・・(3) θRB′=ωR −ωB ・・・(4) ただし、 ωR :リム側部28の角速度 ωR ′:リム側部28の角加速度 ωB :ベルト側部30の角速度 ωB ′:ベルト側部30の角加速度 θRB :リム側部28とベルト側部30とのねじり角 T1 :駆動・制動トルク検出装置68により検出され
る駆動・制動トルク Td :路面からの外乱トルク
Next, the disturbance observer 52 will be described. The disturbance observer 52 is configured based on the model of the wheel 14 shown in FIG. Hereinafter, the configuration of the disturbance observer 52 will be described. The wheels 14, when modeled as being connected by a torsion spring 32 of the belt side 30 Togabane constant K of inertia J R of the rim 28 and moment of inertia J B, (2) ~ ( 4) Is established, which constitutes a linear system. J R ω R '= -Kθ RB + T 1 ··· (2) J B ω B' = Kθ RB -T d ··· (3) θ RB '= ω R -ω B ··· (4) provided that , Ω R : angular velocity of rim side 28 ω R ': angular acceleration of rim side 28 ω B : angular velocity of belt side 30 ω B ': angular acceleration of belt side 30 θ RB : rim side 28 and belt Torsion angle with side 30 T 1 : Driving / braking torque detected by driving / braking torque detecting device 68 T d : Disturbance torque from road surface

【0089】なお、実際にはリム側部28とベルト側部
30との間にはダンパが存在するが、その影響は比較的
小さいため、本実施例においてはその存在が無視されて
いる。
Although a damper actually exists between the rim side portion 28 and the belt side portion 30, its influence is relatively small, and therefore its existence is neglected in this embodiment.

【0090】上記状態方程式をベクトルおよび行列を用
いて表せば(5) 式となる。
If the above state equation is expressed using a vector and a matrix, the following equation (5) is obtained.

【0091】[0091]

【数9】 (Equation 9)

【0092】ここで、タイヤ26の空気圧が変化し、ね
じりばね32のばね定数が正規値であるKからK+ΔK
に変化したときの車輪14の運動は(6) 式で表される。
Here, the air pressure of the tire 26 changes, and the spring constant of the torsion spring 32 changes from K, which is a normal value, to K + ΔK.
The movement of the wheel 14 when the state changes to is expressed by equation (6).

【0093】[0093]

【数10】 (Equation 10)

【0094】すなわち、ばね定数KがΔKだけ変化する
ことは正常なタイヤ26に(6) 式の右辺の最終項で表さ
れる外乱が加えられるのと等価である。この外乱にはば
ね定数Kの変化量ΔKの情報が含まれており、かつ、ば
ね定数Kはタイヤ26の空気圧に応じて変化するので、
この外乱を推定することによってタイヤの空気圧の変化
量を推定することができる。この外乱の推定に外乱オブ
ザーバの手法を用いるのであり、いま路面からのトルク
d をも外乱として扱うことにすれば、推定すべき外乱
wは(7) 式で表される。
That is, the change of the spring constant K by ΔK is equivalent to the addition of a disturbance represented by the last term on the right side of the equation (6) to the normal tire 26. Since this disturbance includes information on the change amount ΔK of the spring constant K, and the spring constant K changes according to the air pressure of the tire 26,
By estimating the disturbance, the amount of change in the tire air pressure can be estimated. The disturbance observer is used for the estimation of the disturbance. If the torque Td from the road surface is also treated as a disturbance, the disturbance w to be estimated is expressed by the following equation (7).

【0095】[0095]

【数11】 [Equation 11]

【0096】しかし、理論上、外乱[w]の中の一つの
要素しか推定することができないため、第2要素である
2 を推定することとする。外乱w2 を(8) 式で定義す
れば、車輪14の状態方程式は(9) 式のようになるた
め、この(9) 式に基づいて外乱オブザーバを構成する。 w2 =(−1/JB )Td +(ΔK/JB )θRB・・・(8)
However, theoretically, only one element of the disturbance [w] can be estimated, so that the second element w 2 is estimated. By defining the disturbance w 2 in (8), for the state equation of the wheel 14 is made to equation (9), constitutes a disturbance observer based on the equation (9). w 2 = (− 1 / J B ) T d + (ΔK / J B ) θ RB (8)

【0097】[0097]

【数12】 (Equation 12)

【0098】外乱オブザーバは外乱をシステムの状態変
数の一つとして推定するものである。そこで、(8) 式の
外乱w2 をシステムの状態に含めるために、推定すべき
外乱のダイナミクスを(10)式で近似する。 w2 ′=0・・・(10) これは図14に示すように連続して変化する外乱を階段
状に近似(零次近似)することを意味し、外乱オブザー
バ52の外乱推定速度を推定すべき外乱の変化に比べて
十分速くすれば、この近似は十分に許容される。(10)式
より、外乱w2をシステムの状態に含めると(11)式の拡
張系が構成される。
The disturbance observer estimates disturbance as one of the state variables of the system. Therefore, (8) for inclusion in the state of the disturbance w 2 system equation to approximate the disturbance dynamics to be estimated by equation (10). w 2 ′ = 0 (10) This means that the continuously changing disturbance is approximated stepwise (zero order approximation) as shown in FIG. 14, and the disturbance estimation speed of the disturbance observer 52 is estimated. This approximation is well tolerated if it is fast enough compared to the change in disturbance to be made. (10) from the equation, the inclusion of the disturbance w 2 to the state of the system (11) extension system is constructed.

【0099】[0099]

【数13】 (Equation 13)

【0100】(11)式において[ωB θRB2T
検出することができない状態となる。したがって、この
システムに基づいて外乱オブザーバ52を構成すれば、
外乱w2 と元々測定できない状態変数ωB ,θRBとを推
定することができる。記述を簡単にするために、(11)式
のベクトルおよび行列を分解して次のように表すことと
する。
In equation (11), [ω B θ RB w 2 ] T cannot be detected. Therefore, if the disturbance observer 52 is configured based on this system,
The disturbance w 2 and the state variables ω B and θ RB that cannot be measured originally can be estimated. In order to simplify the description, the vectors and matrices in equation (11) are decomposed and expressed as follows.

【0101】[0101]

【数14】 [Equation 14]

【0102】このとき、状態[z]=[ωB θRB
2T を推定する最小次元オブザーバの構成は(12)式で
表される。 [zp ′]=[A21][xa ]+[A22][zp ]+[B2 ][u]+[G]{ [xa ′]−([A11][xa ]+[A12][zp ]+[B1 ][u])}=( [A21]−[G][A11])[xa ]+([A22]−[G][A12])[zp ] +[G][xa ′]+([B2 ]−[G][B1 ])[u]・・・(12) ただし、 [zp ] :[z]の推定値 [zp ′]:推定値[zp ]の変化率 [G] :外乱オブザーバ52の推定速度を決めるゲ
イン この方程式をブロック図で表わすと図15のようにな
る。なお、図において[I]は単位行列、sはラプラス
演算子である。また、真値[z]と推定値[zp ]との
誤差[e]を[e]=[z]−[zp]とおき、誤差
[e]の変化率を[e′]とすると、(13)式の関係を得
る。 [e′]=([A22]−[G][A12])[e]・・・(13) これは外乱オブザーバ52の推定特性を表しており、行
列([A22]−[G][A12])の固有値がすなわち外
乱オブザーバ52の極となる。したがって、この固有値
がs平面の左半面において原点から離れるほど外乱オブ
ザーバ52の推定速度が速くなる。オブザーバゲイン
[G]は希望の推定速度になるように決定すればよい。
At this time, the state [z] = [ω B θ RB w
2 ] The configuration of the minimum dimension observer for estimating T is expressed by equation (12). [Z p ′] = [A 21 ] [x a ] + [A 22 ] [z p ] + [B 2 ] [u] + [G] {[x a ′] − ([A 11 ] [x a ] + [A 12 ] [z p ] + [B 1 ] [u])} = ([A 21 ] − [G] [A 11 ]) [x a ] + ([A 22 ] − [G] [ A 12]) [z p] + [G] [x a '] + ([B 2] - [G] [B 1]) [u] ··· (12) However, [z p]: [z ] [Z p ']: change rate of the estimated value [z p ] [G]: gain that determines the estimated speed of the disturbance observer 52 This equation is represented in a block diagram as shown in FIG. In the figure, [I] is a unit matrix, and s is a Laplace operator. If the error [e] between the true value [z] and the estimated value [z p ] is [e] = [z] − [z p ], and the rate of change of the error [e] is [e ′]. , (13) are obtained. [E '] = ([A 22] - [G] [A 12]) [e] ··· (13) which represents the estimated property of the disturbance observer 52, the matrix ([A 22] - [G [A 12 ]) is the pole of the disturbance observer 52. Therefore, the estimated speed of the disturbance observer 52 increases as the eigenvalue moves away from the origin on the left half surface of the s-plane. The observer gain [G] may be determined so as to achieve a desired estimated speed.

【0103】上記外乱推定値w2pおよびねじり角推定値
θRBp を用いて相関演算部54において相関演算が行わ
れ、正規化部56で正規化が行われて、ねじりばね32
のばね定数Kの変化が求められる。
A correlation operation is performed in the correlation operation section 54 using the disturbance estimation value w 2p and the torsion angle estimation value θ RBp , and normalization is performed in the normalization section 56.
Of the spring constant K is determined.

【0104】まず、相関演算部54において、図16の
フローチャートで表されるばね定数変化取得用相関演算
ルーチンが実行される。S21の初期設定において、整
数iが1にリセットされ、前記(8) 式で表される外乱w
2 の推定値w2pとねじり角推定値θRBp との相互相関C
(w2p,θRBp)とねじり角推定値θRBp の自己相関C
(θRBp ,θRBp )とが0にリセットされる。RAM5
0の相互相関メモリおよび自己相関メモリの内容が0に
されるのである。
First, in the correlation calculating section 54, a correlation calculation routine for acquiring a change in spring constant shown in the flowchart of FIG. 16 is executed. In the initial setting of S21, the integer i is reset to 1, and the disturbance w represented by the above equation (8) is obtained.
2 between the estimated value w 2p and the estimated torsion angle θ RBp C
Autocorrelation C between (w 2p , θ RBp ) and estimated torsion angle θ RBp
RBp , θ RBp ) are reset to 0. RAM5
The contents of the zero cross-correlation memory and the auto-correlation memory are set to zero.

【0105】続いて、S22で現時点の外乱推定値w
2p(i) およびねじり角推定値θRBp(i)が読み込まれ、S
23で外乱推定値w2p(i) とねじり角推定値θRBp(i)
の積が演算され、相互相関C(w2p,θRBp )に加算さ
れる。ただし、最初にS23が実行される際には相互相
関C(w2p,θRBp )が0であるため、相互相関メモリ
に外乱推定値w2p(i) とねじり角推定値θRBp(i)との積
が格納されるのみである。同様にS24でねじり角推定
値θRBp(i)の二乗が演算され、自己相関メモリの自己相
関C(θRBp ,θRBp )に加算される。
Subsequently, in S22, the current estimated disturbance value w
2p (i) and the estimated torsion angle θ RBp (i) are read and S
At 23, the product of the disturbance estimation value w 2p (i) and the torsion angle estimation value θ RBp (i) is calculated and added to the cross-correlation C (w 2p , θ RBp ). However, when S23 is first executed, the cross-correlation C (w 2p , θ RBp ) is 0, and therefore, the estimated disturbance w 2p (i) and the estimated torsion angle θ RBp (i) are stored in the cross-correlation memory. Is merely stored. Similarly, the square of the estimated torsion angle θ RBp (i) is calculated in S24 and added to the autocorrelation C (θ RBp , θ RBp ) in the autocorrelation memory.

【0106】S25において整数iが予め定められた整
数M以上になったか否かが判定されるが、当初は判定が
NOであるため、S26で整数iが1増加させられ、再
びS22〜S24が実行される。この実行がM回繰り返
されたときS25の判定がYESとなり、ばね定数変化
取得用相関演算ルーチンの1回の実行が終了する。
At S25, it is determined whether or not the integer i has become equal to or greater than the predetermined integer M. Since the determination is initially NO, the integer i is incremented by 1 at S26, and S22 to S24 are repeated. Be executed. When this execution is repeated M times, the determination in S25 becomes YES, and one execution of the correlation calculation routine for acquiring a spring constant change ends.

【0107】相関演算部54において以上のようにして
相互相関C(w2p,θRBp )と自己相関C(θRBp ,θ
RBp )とが求められた後、正規化部56において(21)式
によりLK 値が求められ、RAM50のLK 値メモリに
格納される。 Lk =C(w2p,θRBp )/C(θRBp ,θRBp )・・・(21) このLK 値は前記(8) 式に基づき、(22)式で表される。 Lk =(−1/JB )C0 +ΔK/JB ・・・(22) ただし、C0 はC(Tdp,θRBp )/C(θRBp ,θ
RBp )で表される値であり、ばね定数Kの変化とは無関
係であるので、タイヤ空気圧が正常の状態で予め求めて
おくことによって補償することができる。また、C(T
dp,θRBp )は外乱トルクTd の推定値とねじり角θRB
の推定値との相互相関を表している。
In the correlation calculation section 54, the cross-correlation C (w 2p , θ RBp ) and the auto-correlation C (θ RBp , θ
RBp ), the L K value is calculated by the normalization unit 56 according to the equation (21) and stored in the L K value memory of the RAM 50. L k = C (w 2p , θ RBp ) / C (θ RBp , θ RBp ) (21) This L K value is expressed by equation (22) based on equation (8). L k = (− 1 / J B ) C 0 + ΔK / J B (22) where C 0 is C (T dp , θ RBp ) / C (θ RBp , θ
RBp ), which is independent of the change in the spring constant K, and can be compensated for by obtaining the tire pressure in a normal state in advance. Also, C (T
dp , θ RBp ) is the estimated value of the disturbance torque T d and the torsion angle θ RB
Represents the cross-correlation with the estimated value.

【0108】ばね定数変化率推定部58においては、以
上のようにして取得され、各L値メモリに格納されてい
るLK =C(w2p,θRBp )/C(θRBp ,θRBp )に
基づいてねじりばね32のばね定数変化量ΔKが推定さ
れる。LK は前述のように、 LK =(−1/JB )C0 +ΔK/JB で表されるため、予めLK とΔKの関係がばね定数変化
テーブルとしてROM49に格納されており、このテー
ブルに基づいてばね定数変化量ΔKが推定されるのであ
る。さらに、ばね定数変化量ΔKを正規値であるKで割
り算することによってばね定数変化率γが求められる。
In the spring constant change rate estimating section 58, L K = C (w 2p , θ RBp ) / C (θ RBp , θ RBp ) obtained as described above and stored in each L value memory. Is used to estimate the spring constant change amount ΔK of the torsion spring 32. As described above, since L K is represented by L K = (− 1 / J B ) C 0 + ΔK / J B , the relationship between L K and ΔK is stored in advance in the ROM 49 as a spring constant change table. The spring constant change amount ΔK is estimated based on this table. Further, the spring constant change rate γ is obtained by dividing the spring constant change amount ΔK by the normal value K.

【0109】推定値補正部60において、回転速度vの
検出値に混入した周期的外乱による影響がばね定数変化
率γの推定値から除去されるように推定値γが補正され
る。回転速度vの検出値自体に混入した周期的外乱は前
記のように、リム側部回転速度演算・補正部45による
信号重ね合わせ処理によって除去可能であるが、周期的
外乱が実質的に完全に除去されるまでに時間がかかり、
その間、回転速度vの検出値に残存する周期的外乱によ
る影響が推定値γに現れることとなる。一方、タイヤ2
6の不均一性に起因する周期的外乱が回転速度vに混入
し、それが回転速度vから除去されない状態において
は、路面凹凸度Rが同じでも車体速度Vが大きい高速走
行状態においてばね定数変化率γの推定値が真の値より
小さくなる傾向があり、また、車体速度Vが同じでも路
面凹凸度Rが低い平坦路走行状態においてばね定数変化
率γの推定値が真の値より小さくなる傾向がある。車体
速度Vが大きい状態または路面凹凸度Rが低い状態にお
いて推定値γに占めるタイヤ26の不均一性による周期
的外乱による誤差値の割合が増加し、周期的外乱の影響
が強く推定値γに現れることになるからである。そこ
で、推定値補正部60においては、信号遅延・重ね合わ
せ処理によって周期的外乱が十分に除去されるまでの
間、すなわち、前記周期的外乱除去完了フラグがON状
態になるまでの間、車体速度Vおよび路面凹凸度Rに基
づいて推定値γが補正されるのである。
In the estimated value correcting section 60, the estimated value γ is corrected so that the influence of the periodic disturbance mixed in the detected value of the rotational speed v is removed from the estimated value of the spring constant change rate γ. As described above, the periodic disturbance mixed in the detected value of the rotational speed v itself can be removed by the signal superimposition process by the rim side rotational speed calculation / correction unit 45, but the periodic disturbance is substantially completely eliminated. It takes time to be removed,
In the meantime, the influence of the periodic disturbance remaining on the detected value of the rotational speed v appears in the estimated value γ. On the other hand, tire 2
In a state in which the periodic disturbance caused by the non-uniformity of No. 6 is mixed with the rotational speed v and is not removed from the rotational speed v, the spring constant changes in a high-speed running state in which the vehicle body speed V is large even if the road surface roughness R is the same. The estimated value of the rate γ tends to be smaller than the true value, and the estimated value of the spring constant change rate γ becomes smaller than the true value in a flat road running state where the road surface unevenness R is low even when the vehicle speed V is the same. Tend. When the vehicle speed V is high or the road surface unevenness R is low, the ratio of the error value due to the periodic disturbance due to the unevenness of the tire 26 to the estimated value γ increases, and the influence of the periodic disturbance is strong to the estimated value γ. Because it will appear. Therefore, the estimated value correction unit 60 controls the vehicle speed until the periodic disturbance is sufficiently removed by the signal delay / superposition process, that is, until the periodic disturbance removal completion flag is turned on. The estimated value γ is corrected based on V and the road surface unevenness R.

【0110】具体的には、推定値γに補正値Δが加算さ
れることによって推定値γの補正が行われるとともに、
その補正値Δが図17(a)にグラフで概念的に表され
ているように、車体速度Vが基準値VS (例えば、10
0〔km/h〕)より大きい高速走行状態においてα(例え
ば、0.1)、基準値VS 以下である低速走行状態にお
いて0となるように決定され、そのαが、図17(b)
に表されているように、路面凹凸度Rが基準値RS 以下
である平坦路走行状態においてα1 、基準値R S より高
い凹凸路走行状態においてα2 (<α1 )となるように
決定される。
Specifically, the correction value Δ is added to the estimated value γ.
The correction of the estimated value γ is performed by
The correction value Δ is conceptually represented by a graph in FIG.
As shown in FIG.S(For example, 10
0 (km / h)) at high speeds greater than α (eg
0.1), reference value VSThe following low-speed driving conditions
Is determined to be 0, and the α is determined as shown in FIG.
As shown in the figure, the road surface roughness R is equal to the reference value R.SLess than
Α on a flat road1, Reference value R SHigher
Α when driving on uneven roadsTwo(<Α1)
It is determined.

【0111】このような推定値補正部60の機能はコン
ピュータ47が図18の推定値補正ルーチンを実行する
ことによって果たされる。まず、S111において、コ
ンピュータ20から周期的外乱除去完了フラグが読み込
まれる。次に、S112において、その周期的外乱除去
完了フラグがON状態にあるか否か、すなわち、周期的
外乱の除去が完了したか否かが判定される。今回は未だ
完了してはいないと仮定すれば判定がNOとなり、S1
13以下のステップに移行する。
The function of the estimated value correcting section 60 is fulfilled by the computer 47 executing the estimated value correcting routine of FIG. First, in step S111, a periodic disturbance removal completion flag is read from the computer 20. Next, in S112, it is determined whether or not the periodic disturbance removal completion flag is in the ON state, that is, whether or not the removal of the periodic disturbance is completed. If it is not completed this time, the determination is NO and S1
Go to step 13 and below.

【0112】S113においては、車体速度センサ70
から車体速度Vが読み込まれ、次に、S114におい
て、その車体速度Vが基準値VS より大きいか否か、す
なわち、車両が高速走行状態にあるか否かが判定され
る。今回は高速走行状態にあると仮定すれば判定がYE
Sとなり、S115において、路面凹凸センサ72から
路面凹凸度Rが読み込まれる。続いて、S116におい
て、補正値Δの決定が行われる。路面凹凸度Rが基準値
S 以下である場合には今回の補正値Δがα1 に、基準
値RS より高い場合にはα2 に決定されるのである。そ
の後、S117において、RAM50からばね定数変化
率γの最新の推定値が読み込まれ、S118において、
その推定値γに決定された補正値Δを加算することによ
って推定値γの補正が行われる。補正後の推定値γはR
AM50に戻される。その後、S111に戻る。
In S113, vehicle speed sensor 70
From the vehicle speed V is read, then, in S114, whether the vehicle speed V is greater than the reference value V S, that is, whether the vehicle is in a high-speed running state is determined. This time, if it is assumed that the vehicle is traveling at high speed, the determination is YE
In S115, the road surface unevenness degree R is read from the road surface unevenness sensor 72 in S115. Subsequently, in S116, the correction value Δ is determined. The correction value Δ is alpha 1 of this time when the road asperity R is less than the reference value R S, is higher than the reference value R S is being determined in alpha 2. Then, in S117, the latest estimated value of the spring constant change rate γ is read from the RAM 50, and in S118,
The estimated value γ is corrected by adding the determined correction value Δ to the estimated value γ. The corrected estimated value γ is R
Returned to AM50. Then, the process returns to S111.

【0113】これに対し、車体速度Vが基準値VS 以下
であり、車両が低速走行状態にある場合には、S114
の判定がNOとなり、S115以下のステップがスキッ
プされ、推定値γの補正が省略される。車両が低速走行
状態にある場合には、推定値γの精度がそれほどには低
下しておらず、補正を行う実益に乏しいからである。
On the other hand, if the vehicle speed V is equal to or lower than the reference value V S and the vehicle is running at a low speed, S114
Is NO, the steps after S115 are skipped, and the correction of the estimated value γ is omitted. This is because, when the vehicle is in the low-speed running state, the accuracy of the estimated value γ does not decrease so much, and the benefit of performing the correction is poor.

【0114】また、周期的外乱除去フラグがON状態に
あって、周期的外乱の除去が完了した場合には、S11
2の判定がYESとなり、S113以下のステップがス
キップされることにより、推定値γの補正が禁止され
る。周期的外乱の除去が完了した後に推定値γの補正を
行うことは、もはや周期的外乱による影響が現れないた
めに補正を行う必要がない推定値γに対して不要な補正
を行うことを意味し、かえって推定値γの精度が低下し
てしまうからである。
If the periodic disturbance elimination flag is in the ON state and the elimination of the periodic disturbance is completed, the process goes to S11.
The determination of 2 is YES, and the steps of S113 and subsequent steps are skipped, so that the correction of the estimated value γ is prohibited. Correcting the estimated value γ after the removal of the periodic disturbance means performing unnecessary correction on the estimated value γ that does not need to be corrected because the influence of the periodic disturbance no longer appears. On the contrary, the accuracy of the estimated value γ decreases.

【0115】異常判定部62においては、推定値補正部
60から供給された推定値γの絶対値が基準値γ0 以上
であるか否かが判定され、基準値γ0 以上である場合に
はタイヤ26の空気圧が異常であると判定され、基準値
γ0 より小さい場合には空気圧が正常であると判定され
る。空気圧が異常であると判定された場合には、その旨
が表示装置66により運転者に知らされる。
[0115] In the abnormality determining unit 62, when the absolute value of the estimated values supplied from the estimated value correction section 60 gamma is determined whether or not the reference value gamma 0 or is the reference value gamma 0 or more tire pressure 26 is determined to be abnormal, when the reference value gamma 0 is smaller than it is determined that the air pressure is normal. If it is determined that the air pressure is abnormal, the driver is notified by the display device 66 to that effect.

【0116】以上の説明から明らかなように、本実施例
においては、電磁ピックアップ12,波形整形器18お
よびコンピュータ20のうち図5のS55およびS56
を実行する部分が、各請求項の発明における「検出装
置」を構成し、また、コンピュータ47のうち外乱オブ
ザーバ52,相関演算部54,正規化部56およびばね
定数変化率推定部58の機能を果たす部分が、各請求項
の発明における「推定手段」を構成し、また、コンピュ
ータ47のうち推定値補正部60の機能を果たす部分が
車体速度センサ70および路面凹凸センサ72と共同し
て、各請求項の発明における「推定値補正手段」を構成
している。さらに、コンピュータ20のうち図5のS5
9〜S64を実行する部分が、各請求項の発明における
「検出値補正手段」を構成しているのである。さらに、
コンピュータ47のうち図18のS112を実行する部
分が、請求項3の発明における「推定値補正許可・禁止
手段」を構成しているのである。
As is apparent from the above description, in the present embodiment, the electromagnetic pickup 12, the waveform shaper 18, and the computer 20 of S55 and S56 of FIG.
Performs the function of the disturbance observer 52, the correlation calculation unit 54, the normalization unit 56, and the spring constant change rate estimation unit 58 of the computer 47. The part that fulfills constitutes the “estimating means” in the invention of each claim, and the part of the computer 47 that fulfills the function of the estimated value correction unit 60 cooperates with the vehicle body speed sensor 70 and the road surface unevenness sensor 72 to It constitutes "estimated value correcting means" in the claimed invention. Further, in the computer 20, S5 in FIG.
The part executing steps 9 to S64 constitutes "detection value correcting means" in the invention of each claim . further,
A part of the computer 47 that executes S112 of FIG.
In the invention of claim 3, "Estimated value correction permission / prohibition
It constitutes "means".

【0117】本実施例においては、ばね定数変化率γの
推定値が、車体速度Vが基準値VSより大きい場合のみ
補正され、その補正値が、路面凹凸度Rが基準値RS
下であるか否かによってα1 ,α2 の2段階に変えられ
るようになっているが、図21(a),(b)に示すよ
うに、車体速度Vが大きいほど、また路面凹凸度Rが小
さいほど補正値Δが大きい値に決定されるようにするこ
とも可能である。例えば、車体速度Vに応じて決定され
た暫定補正値に、路面凹凸度Rに応じて変わる補正係数
を掛けて最終補正値を得たり、ROM49に、車体速度
Vと路面凹凸度Rとをそれぞれ複数領域に分け、それら
領域の組合わせ毎に1個の補正値を割り当てたテーブル
を設けておき、このテーブルに基づいて各時点における
補正値Δが決定されるようにするのである。前者の場合
には補正値Δが連続的に変わり、後者の場合には段階的
に変わることとなる。これらいずれの場合にも、車体速
度Vと路面凹凸度Rとの少なくとも一方において、補正
が行われる領域と行われない領域とを設けることが可能
であり、また、複数の領域(あるいはそれに相当する広
さの領域)について補正値Δが同じ値に決定されるよう
にすることも可能である。さらに、車体速度Vと路面凹
凸度Rとの一方については段階的に、他方については連
続的に補正値が変えられるようにすることも可能であ
る。その他にも、いちいち例示はしないが、各請求項の
発明は種々の改良,変形を加えた態様で実施することが
できる。
In this embodiment, the estimated value of the spring constant change rate γ is corrected only when the vehicle speed V is greater than the reference value V S , and the correction value is obtained when the road surface roughness R is equal to or less than the reference value R S. Depending on whether or not there is, there are two stages, α 1 and α 2. However, as shown in FIGS. 21 (a) and 21 (b), as the vehicle speed V increases, the road surface roughness R increases. It is also possible that the smaller the value, the larger the correction value Δ is determined. For example, the provisional correction value determined according to the vehicle speed V is multiplied by a correction coefficient that varies according to the road surface roughness R to obtain a final correction value, or the ROM 49 stores the vehicle speed V and the road surface roughness R respectively. A table is provided which is divided into a plurality of areas and one correction value is assigned for each combination of these areas, and the correction value Δ at each time is determined based on this table. In the former case, the correction value Δ changes continuously, and in the latter case, it changes stepwise. In any of these cases, it is possible to provide a region where correction is performed and a region where correction is not performed in at least one of the vehicle body speed V and the road surface roughness R, and a plurality of regions (or equivalent regions). It is also possible that the correction value Δ is determined to be the same value for a wide area). Further, the correction value can be changed stepwise for one of the vehicle speed V and the road surface unevenness R, and continuously for the other. Although not specifically exemplified, the invention of each claim can be embodied with various improvements and modifications.

【図面の簡単な説明】[Brief description of the drawings]

【図1】発明の一実施例である車輪情報推定装置の機
能ブロック図である。
FIG. 1 is a functional block diagram of a wheel information estimation device according to an embodiment of the present invention.

【図2】上記車輪情報推定装置の構成ブロック図であ
る。
FIG. 2 is a configuration block diagram of the wheel information estimation device.

【図3】上記車輪情報推定装置により外乱を検出される
車輪の一部を示す断面図である。
FIG. 3 is a sectional view showing a part of a wheel whose disturbance is detected by the wheel information estimating device.

【図4】上記車輪の力学モデルを示す図である。FIG. 4 is a diagram showing a dynamic model of the wheel.

【図5】上記車輪情報推定装置の一構成要素であるコン
ピュータのROMに格納されている制御プログラムを示
すフローチャートである。
FIG. 5 is a flowchart showing a control program stored in a ROM of a computer which is a component of the wheel information estimation device.

【図6】車輪の回転速度信号の周波数特性の一例を示す
グラフである。
FIG. 6 is a graph showing an example of a frequency characteristic of a wheel rotation speed signal.

【図7】上記車輪情報推定装置における周期的外乱除去
部(信号遅延・重ね合わせ処理部)の機能ブロック図で
ある。
FIG. 7 is a functional block diagram of a periodic disturbance remover (signal delay / superposition processor) in the wheel information estimating device.

【図8】上記周期的外乱除去部の特性の一例を示すグラ
フである。
FIG. 8 is a graph showing an example of a characteristic of the periodic disturbance removing unit.

【図9】上記車輪情報推定装置の一構成要素である1回
転分メモリの構成とその使用方法とを概念的に示す図で
ある。
FIG. 9 is a diagram conceptually showing a configuration of a memory for one rotation, which is a component of the wheel information estimation device, and a method of using the same.

【図10】上記1回転分メモリにおける記憶内容とサン
プリング周期との関係を説明するための図である。
FIG. 10 is a diagram for explaining the relationship between the contents stored in the memory for one rotation and the sampling period.

【図11】上記1回転分メモリにおける記憶内容とサン
プリング周期との関係を説明するための別の図である。
FIG. 11 is another diagram for explaining the relationship between the contents stored in the one-rotation memory and the sampling period.

【図12】上記周期的外乱除去部の原理を概念的に示す
グラフである。
FIG. 12 is a graph conceptually showing the principle of the periodic disturbance removing unit.

【図13】上記周期的外乱除去部の効果の一例を示すグ
ラフである。
FIG. 13 is a graph showing an example of the effect of the periodic disturbance removing unit.

【図14】上記車輪情報推定装置における外乱のダイナ
ミクスの近似を説明するためのグラフ図である。
FIG. 14 is a graph for explaining approximation of disturbance dynamics in the wheel information estimation device.

【図15】上記車輪情報推定装置における外乱オブザー
バを示すブロック線図である。
FIG. 15 is a block diagram showing a disturbance observer in the wheel information estimation device.

【図16】上記車輪情報推定装置の一構成要素であるコ
ンピュータのROMに格納されている制御プログラムを
示すフローチャートである。
FIG. 16 is a flowchart showing a control program stored in a ROM of a computer which is a component of the wheel information estimating apparatus.

【図17】上記車輪情報推定装置における推定値の補正
の一例を概念的に示すグラフである。
FIG. 17 is a graph conceptually showing an example of correction of an estimated value in the wheel information estimating device.

【図18】上記車輪情報推定装置の一構成要素であるコ
ンピュータのROMに格納されている制御プログラムを
示すフローチャートである。
FIG. 18 is a flowchart showing a control program stored in a ROM of a computer which is a component of the wheel information estimation device.

【図19】請求項1および2の各発明の構成を概念的に
示すブロック図である。
FIG. 19 is a block diagram conceptually showing the configuration of each of the first and second aspects of the present invention.

【図20】請求項3の発明の構成を概念的に示すブロッ
ク図である。
FIG. 20 is a block diagram conceptually showing a configuration of the invention of claim 3;

【図21】前記車輪情報推定装置における推定値の補正
の別の例を概念的に示すグラフである。
FIG. 21 is a graph conceptually showing another example of the correction of the estimated value in the wheel information estimating device.

【符号の説明】[Explanation of symbols]

10 ロータ 12 電磁ピックアップ 14 車輪(タイヤ付ホイール) 20 コンピュータ 26 タイヤ 47 コンピュータ 52 外乱オブザーバ 58 ばね定数変化率推定部 60 推定値補正部 70 車体速度センサ 72 路面凹凸センサ Reference Signs List 10 rotor 12 electromagnetic pickup 14 wheels (wheels with tires) 20 computer 26 tires 47 computer 52 disturbance observer 58 spring constant change rate estimating unit 60 estimated value correcting unit 70 vehicle speed sensor 72 road surface unevenness sensor

───────────────────────────────────────────────────── フロントページの続き (72)発明者 河井 弘之 愛知県豊田市トヨタ町1番地 トヨタ自 動車株式会社内 (72)発明者 小島 弘義 愛知県豊田市トヨタ町1番地 トヨタ自 動車株式会社内 (72)発明者 梅野 孝治 愛知県愛知郡長久手町大字長湫字横道41 番地の1株式会社豊田中央研究所内 (72)発明者 浅野 勝宏 愛知県愛知郡長久手町大字長湫字横道41 番地の1株式会社豊田中央研究所内 (72)発明者 内藤 俊治 愛知県刈谷市昭和町1丁目1番地 日本 電装株式会社内 (72)発明者 小野木 伸好 愛知県刈谷市昭和町1丁目1番地 日本 電装株式会社内 (56)参考文献 特開 平5−294119(JP,A) 特開 平5−213018(JP,A) 特開 平5−133831(JP,A) 特開 平7−172121(JP,A) 特開 平6−320923(JP,A) (58)調査した分野(Int.Cl.7,DB名) B60C 23/00 - 23/08 G01L 17/00 B60T 8/00 - 8/96 ──────────────────────────────────────────────────続 き Continuing on the front page (72) Inventor Hiroyuki Kawai 1 Toyota Town, Toyota City, Aichi Prefecture Inside Toyota Motor Corporation (72) Inventor Hiroyoshi Kojima 1 Toyota Town, Toyota City, Aichi Prefecture Inside Toyota Motor Corporation ( 72) Inventor Koji Umeno 41 Toyota Chuo R & D Laboratories Co., Ltd. at 41, Chuchu-ji, Yokomichi, Nagakute-cho, Aichi-gun, Aichi Prefecture (72) Inventor Katsuhiro Asano 41 Co., Ltd. Inside the Central Research Laboratory (72) Inventor Shunji Naito 1-1-1, Showa-cho, Kariya-shi, Aichi, Japan Inside Denso Co., Ltd. (72) Inventor Nobuyoshi Onoki 1-1-1, Showa-cho, Kariya-shi, Aichi Japan Inside Denso Corporation (56 ) References JP-A-5-294119 (JP, A) JP-A-5-213018 (JP, A) JP-A-5-133831 (JP, A) JP-A-7-172121 (JP, A) JP-A-6-320923 (JP, A) (58) Fields investigated (Int. Cl. 7 , DB name) B60C 23/00-23/08 G01L 17 / 00 B60T 8/00-8/96

Claims (3)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】車輪の回転速度を検出する検出装置と、 その検出装置による検出値に基づき、前記車輪に関する
情報を推定する推定手段とを含む車輪情報推定装置にお
いて、前記検出装置による検出値を、それに混入した周期的外
乱が除去されるように補正して前記推定手段に供給する
検出値補正手段と、 前記推定手段による推定値を、車両の走行状況に基づ
き、前記検出値補正手段により除去されないで前記検出
装置による検出値に残る周期的外乱による影響が推定値
現れないように補正する推定値補正手段を設けたこ
とを特徴とする車輪情報推定装置。
1. A wheel information estimating device comprising: a detecting device for detecting a rotational speed of a wheel; and estimating means for estimating information on the wheel based on a detected value by the detecting device. , Periodic out mixed with it
Correction is made so that disturbance is removed and supplied to the estimating means.
A detection value correction unit configured to detect the value estimated by the estimation unit without being removed by the detection value correction unit based on a traveling state of the vehicle;
Estimated value of the effect of periodic disturbance remaining in the value detected by the device
Wheel information estimating apparatus is characterized by providing an estimated value correction means for correcting so as not to appear in.
【請求項2】請求項1に記載の車輪情報推定装置であっ
て、前記推定値補正手段が、前記推定手段による推定値
の補正量を前記走行状況としての車体速度が大きい場合
において小さい場合におけるより大きく決定すること
と、前記走行状況としての路面凹凸度が低い場合におい
て高い場合におけるより大きく決定することとの少なく
とも一方を行うものである車輪情報推定装置。
2. The wheel information estimating device according to claim 1, wherein the estimated value correcting means determines whether or not the amount of correction of the estimated value by the estimating means is small when the vehicle speed as the traveling condition is large. A wheel information estimating device that performs at least one of determining a larger value and determining a larger value when the road surface unevenness degree as the traveling condition is low and high.
【請求項3】請求項1または2に記載の車輪情報推定装
置であって、さらに、前記検出値補正手段により除去さ
れないで前記検出装置による検出値に残る周期的外乱
許容値より大きい間は前記推定値補正手段による補正を
許可し、許容値以下になった後には禁止する推定値補正
許可・禁止手段を含む車輪情報推定装置。
3. The wheel information estimating device according to claim 1, further comprising:
The estimated value correction permitting / prohibiting means for permitting the correction by the estimated value correcting means while the periodic disturbance remaining in the detected value by the detecting device is larger than the permissible value, and prohibiting the correction after the value becomes equal to or less than the permissible value. Including wheel information estimation device.
JP17125894A 1994-07-22 1994-07-22 Wheel information estimation device Expired - Fee Related JP3155154B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP17125894A JP3155154B2 (en) 1994-07-22 1994-07-22 Wheel information estimation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP17125894A JP3155154B2 (en) 1994-07-22 1994-07-22 Wheel information estimation device

Publications (2)

Publication Number Publication Date
JPH0834215A JPH0834215A (en) 1996-02-06
JP3155154B2 true JP3155154B2 (en) 2001-04-09

Family

ID=15919992

Family Applications (1)

Application Number Title Priority Date Filing Date
JP17125894A Expired - Fee Related JP3155154B2 (en) 1994-07-22 1994-07-22 Wheel information estimation device

Country Status (1)

Country Link
JP (1) JP3155154B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5068094B2 (en) * 2007-03-27 2012-11-07 カヤバ工業株式会社 Pressure control device
KR101558975B1 (en) 2013-12-19 2015-10-12 현대오트론 주식회사 Pretreatment apparatus and method of wheel speed
CN104369629B (en) * 2014-10-21 2017-01-18 中山市六源通电子科技有限公司 Road condition judgment prompt method based on tire pressure sensing

Also Published As

Publication number Publication date
JPH0834215A (en) 1996-02-06

Similar Documents

Publication Publication Date Title
JP3300601B2 (en) Tire pressure detector
JP2952151B2 (en) Wheel disturbance detection device and its use
CN101596916B (en) Tire attitude control device and method
JP4474475B2 (en) Tire pressure drop detection device and method, and tire pressure drop detection program
US20110107828A1 (en) Apparatus and method for detecting tire having decreased pressure, and program for detecting tire having decreased pressure
JP3155154B2 (en) Wheel information estimation device
JP2022100865A (en) Tyre wear state estimation device
JP2984546B2 (en) Wheel information estimation device
JP3135472B2 (en) Wheel speed detector
JP3329171B2 (en) Rotation speed detector
JP2002082123A (en) Wheel rotation information detector
JP3001372B2 (en) Wheel speed detector
JP2003326927A (en) Detecting method and device for tire pneumatic pressure drop and program to determine tire pressure reduction
JP3300572B2 (en) Tire pressure estimation device
JP2004203214A (en) Method and device for detecting drop of pneumatic pressure of tire and program for determining pressure drop of tire
JP3305878B2 (en) Wheel information acquisition device
JPH09295566A (en) Vehicle control method
JP3391486B2 (en) Tire pressure detector
JP2005153544A (en) Calculation method and device for loading sensitivity of dynamic load radius of tire, and program for calculating loading sensitivity of tire
JPH10193933A (en) Tire pneumatic pressure estimation device
JP3300532B2 (en) Tire pressure detector
JP3131201B2 (en) Tire pressure drop detection method and device
JPH0783949A (en) Apparatus for detecting yaw rate of vehicle
JP7434909B2 (en) Road surface condition determination device, determination method, and determination program
JP2003154828A (en) Tire air pressure detector

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees