JP3228168B2 - Rotation speed detection device and tire generation force detection device - Google Patents

Rotation speed detection device and tire generation force detection device

Info

Publication number
JP3228168B2
JP3228168B2 JP04610097A JP4610097A JP3228168B2 JP 3228168 B2 JP3228168 B2 JP 3228168B2 JP 04610097 A JP04610097 A JP 04610097A JP 4610097 A JP4610097 A JP 4610097A JP 3228168 B2 JP3228168 B2 JP 3228168B2
Authority
JP
Japan
Prior art keywords
error
speed
value
torque
estimated
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
JP04610097A
Other languages
Japanese (ja)
Other versions
JPH10239187A (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.)
Toyota Central R&D Labs Inc
Original Assignee
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 Toyota Central R&D Labs Inc filed Critical Toyota Central R&D Labs Inc
Priority to JP04610097A priority Critical patent/JP3228168B2/en
Publication of JPH10239187A publication Critical patent/JPH10239187A/en
Application granted granted Critical
Publication of JP3228168B2 publication Critical patent/JP3228168B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Force Measurement Appropriate To Specific Purposes (AREA)
  • Control Of Electric Motors In General (AREA)

Description

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

【0001】[0001]

【発明の属する技術分野】本発明は回転速度検出装置及
びタイヤ発生力検出装置に係り、特に、回転体の力学的
モデルを用いて構成されたオブザーバにより回転体の瞬
時速度を推定する回転速度検出装置及び該装置を用いて
車輪に作用する制動力を推定するタイヤ発生力検出装置
に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a rotational speed detecting device and a tire generating force detecting device, and more particularly to a rotational speed detecting device for estimating an instantaneous speed of a rotating body by an observer constructed using a dynamic model of the rotating body. The present invention relates to a device and a tire generation force detecting device for estimating a braking force acting on wheels using the device.

【0002】[0002]

【従来の技術】従来より、電動機や車両等の回転体の回
転速度を検出する方法として、いわゆるパルスエンコー
ダ、パルスジェネレータといった位置検出器を用いて回
転体の回転位置(回転角度)を検出し、検出された回転
位置(以下、単に「位置」ともいう)の値に基づいて回
転速度を検出する方法がある。一般的には、ある時間周
期で位置検出器の値が読み込まれ、その読み込みタイミ
ングにおける検出値と1タイミング前における検出値と
の差分によって、回転速度(以下、単に「速度」ともい
う)が検出される。
2. Description of the Related Art Conventionally, as a method of detecting the rotation speed of a rotating body such as an electric motor or a vehicle, the position (rotation angle) of the rotating body is detected by using a position detector such as a so-called pulse encoder or pulse generator. There is a method of detecting a rotation speed based on a value of a detected rotation position (hereinafter, also simply referred to as “position”). Generally, the value of the position detector is read at a certain time period, and the rotation speed (hereinafter, also simply referred to as “speed”) is detected by the difference between the detected value at the read timing and the detected value one timing before. Is done.

【0003】しかし、上記方法の場合、本来、連続量で
ある位置が、位置検出器の分解能に規定される離散値で
検出されるため、検出速度の時間遅れや量子化ノイズと
いった問題が生じる。さらに、速度が低速になると、位
置読み込み周期間で十分なパルス数が得られないため、
上記問題はより深刻なものとなる。
However, in the case of the above method, since a position which is originally a continuous quantity is detected by a discrete value defined by the resolution of the position detector, problems such as a time delay of the detection speed and quantization noise occur. Furthermore, when the speed is low, a sufficient number of pulses cannot be obtained between the position reading cycles.
The above problem is exacerbated.

【0004】そこで、上記問題を回避するため、回転体
のモデルに基づいて、いわゆるオブザーバ(状態観測
器)を構成し、該オブザーバによりパルスが来ない間の
瞬時速度を推定する「瞬時速度オブザーバ」の手法が提
案されている(今野,堀:「高次外乱補償機能を有する
瞬時速度オブザーバ」,電気学会論文誌D,112,No6,平
成4 年)。以下、本手法について説明する。なお、上記
文献記載の本手法は、電動機制御における速度検出につ
いて述べられているので、以下の説明では回転体を電動
機として説明する。
Therefore, in order to avoid the above problem, a so-called observer (state observer) is constructed based on a model of the rotating body, and an "instantaneous speed observer" for estimating an instantaneous speed during which no pulse comes by the observer. (Konno, Hori: "Instantaneous speed observer with higher-order disturbance compensation function", IEEJ Transactions on Electronics, D, 112, No6, 1994). Hereinafter, this method will be described. Since the method described in the above-mentioned document describes speed detection in electric motor control, the following description will be made assuming that the rotating body is an electric motor.

【0005】図5には、瞬時速度オブザーバ10の構
成、電動機12及び位置検出器14が示されている。同
図に示すように、入力された電流iにより回転駆動した
電動子の位置を求めるとき、電動機12は、電流iから
電動子の駆動トルクへの変換係数を示すトルク係数K
と、該駆動トルク及び電動子に外部から作用する外乱ト
ルクTd の加算要素と、電動子の慣性J及び積分要素1
/sからなる変換要素(1/Js)と、積分要素1/s
と、により表される。なお、sはラプラス変換の演算子
に相当している。
FIG. 5 shows the structure of the instantaneous speed observer 10, the electric motor 12, and the position detector 14. As shown in the figure, when obtaining the position of the motor driven to rotate by the input current i, the motor 12 generates a torque coefficient K indicating a conversion coefficient from the current i to the driving torque of the motor.
An adding element of the driving torque and a disturbance torque Td acting externally on the motor, an inertia J of the motor and an integral element 1
/ S conversion element (1 / Js) and integral element 1 / s
And is represented by Note that s corresponds to the Laplace transform operator.

【0006】この電動機12の作用を説明すると、まず
電流iを入力すると、電流iのエネルギーは、電動子の
駆動トルク(Ki)に変換される。この駆動トルクに対
し、電動機外部から作用する外乱トルクTd が加わる
と、実際に電動子に作用するトルクは全トルクTm とな
る。次に、演算要素(1/Js)を慣性Jの逆数(1/
J)と積分要素(1/s)とに分解して考えると、全ト
ルクTm は、(1/J)により電動子の回転角速度に変
換され、さらに、この回転角速度は積分要素(1/s)
により電動子の回転速度ωに変換される。そして、回転
速度ωは、後段の積分要素(1/s)により連続量とし
ての実際の回転位置に変換される。
The operation of the motor 12 will be described. First, when a current i is input, the energy of the current i is converted into a driving torque (Ki) of the motor. For this driving torque, the disturbance torque T d applied from the electric motor externally applied, the torque acting on the actual electric element becomes total torque T m. Next, the operation element (1 / Js) is calculated as the reciprocal of inertia J (1 / Js).
J) and the integral element (1 / s), the total torque Tm is converted into the rotational angular velocity of the armature by (1 / J). s)
Is converted into the rotational speed ω of the armature. Then, the rotation speed ω is converted into an actual rotation position as a continuous amount by an integral element (1 / s) at the subsequent stage.

【0007】電動子の実際の回転位置は、位置検出器1
4により離散値の回転位置θ(パルス計数値)として位
置読み取り周期T1 毎に検出される。
The actual rotational position of the armature is determined by the position detector 1
4 is detected as a discrete rotational position θ (pulse count value) at each position reading cycle T 1 .

【0008】瞬時速度オブザーバ10は、上記した電動
機12の等価モデル16を有しており、この等価モデル
16は、トルク係数Kn (トルク係数Kのノミナル値を
意味する)の定数器20と、駆動トルクの演算値及び外
乱トルクの推定値を加算する加算器22と、電動子に作
用する全トルクの推定値に対して(1/Jn s)を演算
する演算器24(Jn は慣性Jのノミナル値)と、演算
器24の演算結果に対して(1/s)を演算する積分器
26と、から構成される。なお、図5の瞬時速度オブザ
ーバ10では、回転速度ωの推定値、回転位置θの推定
値、外乱トルクTd の推定値及び全トルクTm の推定値
を、各変数に∧を付与した変数で表している。
The instantaneous speed observer 10 has an equivalent model 16 of the electric motor 12 described above. The equivalent model 16 includes a constant unit 20 for a torque coefficient K n (meaning a nominal value of the torque coefficient K), an adder 22 for adding the estimated value of the calculated value and the disturbance torque of the drive torque, calculator 24 (J n inertia for calculating a (1 / J n s) for the estimated value of the total torque acting on the armature J) (nominal value of J) and an integrator 26 that calculates (1 / s) for the calculation result of the calculator 24. In the instantaneous speed observer 10 of FIG. 5, the estimated value of the rotation speed omega, the estimated value of the rotational position theta, an estimate of the disturbance torque T d and the estimated value of the total torque T m, was applied a ∧ to each Variable It is represented by

【0009】等価モデル16では、電流i及び外乱トル
クの推定値の入力により、制御系の制御周期T2 毎に演
算器24が電動子の速度の推定値を出力し、さらに積分
器26が位置の推定値を出力する。すなわち、等価モデ
ル16によれば、制御周期T 2 毎に上記推定値を出力す
るので、位置検出器14の読み取り周期T1 の間の位置
情報が得られない間でも位置及び速度の推定値を得るこ
とができる。
In the equivalent model 16, the current i and the disturbance torque
The control period T of the control systemTwoEvery performance
The calculator 24 outputs the estimated value of the speed of the armature and further integrates it.
The detector 26 outputs the position estimate. That is, the equivalent model
According to rule 16, the control cycle T TwoOutput the above estimated value every time
Therefore, the reading cycle T of the position detector 141Position between
Obtain position and velocity estimates while no information is available.
Can be.

【0010】さらに、瞬時速度オブザーバ10は、等価
モデル16により演算された位置の推定値と位置検出器
14により検出された位置θとの偏差Δθを演算する偏
差器28と、演算されたΔθに基づいて、読み取り周期
1 毎に、等価モデル16で用いられる外乱トルクの推
定値及び該等価モデル16で推定される速度の推定値を
修正するためのゲインを各々演算する修正ゲイン演算部
18と、を備えている。
Further, the instantaneous speed observer 10 includes a deviation unit 28 for calculating a deviation Δθ between the position estimated value calculated by the equivalent model 16 and the position θ detected by the position detector 14, A correction gain calculator 18 for calculating the estimated value of the disturbance torque used in the equivalent model 16 and the gain for correcting the estimated value of the speed estimated by the equivalent model 16 for each reading cycle T 1 , Is provided.

【0011】瞬時速度オブザーバ10では、修正ゲイン
演算部18により演算されたゲインにより外乱トルクの
推定値及び速度の推定値が、Δθが零に一致するように
修正されるので、等価モデル16の各推定値は、現実の
電動機12の各出力値と略一致するようになり、よっ
て、回転速度及び外乱トルクの正確な推定が可能とな
る。
In the instantaneous speed observer 10, the estimated value of the disturbance torque and the estimated value of the speed are corrected by the gain calculated by the correction gain calculator 18 so that Δθ becomes equal to zero. The estimated value substantially coincides with the actual output value of the electric motor 12, so that the rotation speed and the disturbance torque can be accurately estimated.

【0012】図5の瞬時速度オブザーバ10として、上
記文献には、2つの態様が示されている。第1の態様は
「位置読み込み形」であり、第2の態様は「平均速度読
み込み形」である。以下、各態様についてさらに詳細に
説明する。
As the instantaneous speed observer 10 shown in FIG. 5, the above document discloses two modes. The first mode is a “position reading type”, and the second mode is an “average speed reading type”. Hereinafter, each embodiment will be described in more detail.

【0013】まず、第1の態様「位置読み込み形」につ
いて説明する。図6には、位置θに対応するパルス計数
値、位置検出器14による位置の読み込みタイミング、
及び制御周期が示されている。ここで、読み取り周期T
1 と制御周期T2 との間には、 T1 = mT2 ,mは整数 (1) という関係があるものとする。読み取り周期T1 は、速
度分解能が十分得られるように数[ms]〜数十[m
s]に設定されている。また、本手法では、DSP(Di
gital Signal Processor)などの高速の演算デバイスの
使用を前提としているので、制御周期T2 は数十[μ
s]程度となる。なお、第1の態様では、位置検出器1
4から発生するパルスの間隔は、制御周期T2 よりも短
いと仮定している。
First, the first mode "position reading type" will be described. FIG. 6 shows the pulse count value corresponding to the position θ, the timing of reading the position by the position detector 14,
And the control cycle. Here, the reading cycle T
It is assumed that T 1 = mT 2 , where m is an integer (1), between 1 and the control cycle T 2 . Reading cycle T 1, the number as velocity resolution can be obtained sufficiently [ms] ~ several tens [m
s]. In this method, the DSP (Di-
Since assumes the use of high-speed computing device gital Signal Processor), etc., the control period T 2 are tens [mu
s]. In the first embodiment, the position detector 1
Pulse interval generated from 4 assumes shorter than the control period T 2.

【0014】位置検出器14の読み込み周期内では、時
間tは、 t = jT1 +kT2 ,0≦k<m (2) となる。但し、jは読み取り周期T1 の各区間に順番に
付与されたインデックス(番号)、kは読み取り周期T
1 の1区間内で生じたm個の制御周期T2 の各区間に順
番に付与されたインデックスである。j、kを指定する
と、(2) 式より時間tが定まるので、以下では、時間t
における一般の変数f(t)をf[j,k]と表すこと
にする。なお、f[j,m]=f[j+1,0]となる
ので、(2)式では、kの定義域に0を含ませている。
In the reading cycle of the position detector 14, the time t is as follows: t = jT 1 + kT 2 , 0 ≦ k <m (2) Here, j is an index (number) sequentially assigned to each section of the reading cycle T 1 , and k is a reading cycle T
The m each section of the control period T 2 of the occurring within 1 of the first interval is an index assigned to the order. When j and k are specified, the time t is determined from the equation (2).
Is represented as f [j, k]. Since f [j, m] = f [j + 1,0], 0 is included in the domain of k in equation (2).

【0015】また、回転速度の推定のためには、制御周
期T2 毎に外乱の推定値が必要となるが、簡単のため、
Further, for the estimation of the rotation speed is the estimated value of the disturbance in each control period T 2 is required, for simplicity,

【0016】[0016]

【数1】 (Equation 1)

【0017】と仮定する。(3) 式は、間隔T2 において
外乱が一定値をとることを意味している。この仮定はT
2 が十分短い場合には、有効な近似となる。
Suppose that (3) means that the disturbance takes a constant value in the interval T 2. This assumption is T
If 2 is short enough, it is a valid approximation.

【0018】ここで、時点[j,k]で電動子に作用す
る全トルクTm [j,k]の推定値は、定数器20及び
加算器22により、
Here, the estimated value of the total torque T m [j, k] acting on the armature at the time point [j, k] is calculated by the constant unit 20 and the adder 22.

【0019】[0019]

【数2】 (Equation 2)

【0020】と演算される。但し、i[j,k]は、時
点[j,k]での電流値である。そこで、演算器24
が、全トルクの推定値を台形積分することにより、次式
のような速度の推定値を得る。
Is calculated. Here, i [j, k] is the current value at the time point [j, k]. Therefore, the arithmetic unit 24
However, by performing trapezoidal integration of the estimated value of the total torque, an estimated value of the speed as in the following equation is obtained.

【0021】[0021]

【数3】 (Equation 3)

【0022】次に、積分器26が(5) 式の速度推定値を
台形積分することにより、次式のような位置の推定値を
得る。
Next, the integrator 26 trapezoidally integrates the velocity estimation value of the equation (5) to obtain a position estimation value as shown in the following equation.

【0023】[0023]

【数4】 (Equation 4)

【0024】なお、(5) 、(6) 式の台形積分は、上述し
たように制御周期T2 毎に行われる。
The trapezoidal integration of the equations (5) and (6) is performed at every control cycle T 2 as described above.

【0025】次に、(5) 式による速度ω[j,m]の推
定値(∧付)及び(6) 式による位置θ[j,m]の推定
値(∧付)には、推定誤差が含まれているので以下のよ
うに修正を行う。
Next, the estimated value of the velocity ω [j, m] according to the equation (5) (appended) and the estimated value of the position θ [j, m] according to the equation (6) (appended) are estimated errors. Is included, so modify as follows.

【0026】[j,m]の時点において、位置検出器1
4の読み込みにより正確な位置θ[j,m]=θ[j+
1,0]がわかるので、偏差器28により位置の推定誤
At the time point [j, m], the position detector 1
4 to obtain an accurate position θ [j, m] = θ [j +
[0], the position estimation error is obtained by the deviation device 28.

【0027】[0027]

【数5】 (Equation 5)

【0028】が得られる。そこで、この推定誤差が生じ
る原因を、推定外乱ベクトルの初期値Td [j,0]
(∧付)及び速度推定値の初期値ω[j,0](∧付)
の誤差によるものと考える。Δθを各初期値誤差に分け
て書くと次式のように表せる。
Is obtained. The cause of the estimation error is determined by the initial value T d [j, 0] of the estimated disturbance vector.
(∧) and initial value ω [j, 0] of estimated speed (∧)
It is thought to be due to the error of If Δθ is written separately for each initial value error, it can be expressed as the following equation.

【0029】[0029]

【数6】 (Equation 6)

【0030】である。ここに、λは誤差Δθを各初期値
誤差に分配するための重み係数ベクトルである。(8) 式
に基づいてxerr を計算し、これを以下のように次の速
度推定期間[j+1]における初期値の修正に用いる。
## EQU1 ## Here, λ is a weight coefficient vector for distributing the error Δθ to each initial value error. X err is calculated based on the equation (8), and is used for correcting the initial value in the next speed estimation period [j + 1] as follows.

【0031】[0031]

【数7】 (Equation 7)

【0032】である。従って、修正演算式は次式のよう
になる。
Is as follows. Therefore, the correction operation expression is as follows.

【0033】[0033]

【数8】 (Equation 8)

【0034】(10)式の演算は、位置検出器14の読み込
み周期毎に行う。また、同式右辺第2項の係数ベクトル
が、修正ゲイン演算部18が演算する修正ゲインに相当
する。
The calculation of equation (10) is performed for each reading cycle of the position detector 14. In addition, the coefficient vector of the second term on the right side of the equation corresponds to the correction gain calculated by the correction gain calculator 18.

【0035】なお、上記演算に必要なλは以下のように
して決定する。瞬時速度オブザーバ10は、制御理論的
観点からみると一種の最小次元オブザーバとみなすこと
ができる。すなわち、電動機12の状態量は、位置、速
度、外乱の3つであるが、位置は、位置検出器14の読
み込み時点でわかるので、速度と外乱のみを推定する最
小の構成となっている。
Note that λ required for the above calculation is determined as follows. The instantaneous speed observer 10 can be regarded as a kind of minimum-dimensional observer from the viewpoint of control theory. That is, the state quantity of the motor 12 is three, that is, position, speed, and disturbance. However, since the position is known at the time of reading the position detector 14, it has a minimum configuration that estimates only the speed and disturbance.

【0036】この瞬時速度オブザーバ10の推定誤差の
収束特性は、
The convergence characteristics of the estimation error of the instantaneous speed observer 10 are as follows.

【0037】[0037]

【数9】 (Equation 9)

【0038】で与えられる。ここに、A、Cは電動機1
2の物理パラメータで決まる定数、Gはオブザーバゲイ
ン行列、x[j]はシステムの状態変数(∧付きは該状
態変数の推定値)である。
Is given by Here, A and C are electric motors 1
G is an observer gain matrix, and x [j] is a state variable of the system (with an ∧ indicates an estimated value of the state variable).

【0039】オブザーバの極は、次の特性方程式の根と
して求められる。なお、Iは単位行列である。
The pole of the observer is obtained as the root of the following characteristic equation. Here, I is a unit matrix.

【0040】 |zI−(A−GCA)| = 0 (12) ここで、上記特性方程式の根が、z平面における原点を
中心とした半径1の円内に入るようにゲイン行列Gを設
計すれば、推定誤差は時間と共に0に収束し、瞬時速度
及び外乱を推定できるようになる。そして、ゲイン行列
Gと誤差配分行列λとの間には、G=LK-1λという関
係があるので、この関係からλを決定することができ
る。
| ZI− (A−GCA) | = 0 (12) Here, the gain matrix G is designed so that the root of the above characteristic equation falls within a circle having a radius of 1 centered on the origin on the z plane. For example, the estimation error converges to 0 with time, and the instantaneous speed and the disturbance can be estimated. Since there is a relationship G = LK -1 λ between the gain matrix G and the error distribution matrix λ, λ can be determined from this relationship.

【0041】次に、第2の態様「平均速度読み込み形オ
ブザーバ」について説明する。平均速度読み込み形の位
置検出器14は、以下で述べるように、エンコーダのパ
ルス間隔(隣接するパルス間の時間間隔)を測定する機
能を有しており、測定されたパルス間隔より回転位置θ
を求める。このように平均速度読み込み形オブザーバで
は、パルス間隔毎に位置θを求めるので、パルス間隔を
1 としている。この場合、T1 は可変となる。また、
平均速度読み込み形オブザーバでは、エンコーダのパル
ス間隔が制御周期T2 より長くなる場合を想定してお
り、エンコーダパルスが発生した直後の制御時点で推定
値の修正を行う。
Next, a second mode "average speed reading type observer" will be described. The average speed reading type position detector 14 has a function of measuring the pulse interval of the encoder (the time interval between adjacent pulses) as described below, and the rotational position θ is determined based on the measured pulse interval.
Ask for. As described above, in the average speed reading type observer, since the position θ is obtained for each pulse interval, the pulse interval is set to T 1 . In this case, T 1 is made variable. Also,
The average speed reading type observer assumes the case where the pulse interval of the encoder is longer than the control period T 2, to correct the estimated value at the control point immediately after the encoder pulse occurs.

【0042】図7に、平均速度読み込み形オブザーバの
タイミングチャートを示す。同図において、エンコーダ
のパルス間隔を測定するための高速カウンタはパルスの
到着と同時にリセットされるものとする。高速カウンタ
によりパルス間隔の測定ができ、また、1パルス当たり
の回転角度も予めわかっているので、両者よりその間の
平均速度<ω[j+1]>が以下のように計算される。
FIG. 7 is a timing chart of the average speed reading type observer. In the figure, it is assumed that the high-speed counter for measuring the pulse interval of the encoder is reset at the same time as the arrival of the pulse. Since the pulse interval can be measured by the high-speed counter and the rotation angle per pulse is also known in advance, the average speed <ω [j + 1]> between them is calculated as follows.

【0043】[0043]

【数10】 (Equation 10)

【0044】そして、パルスの到着直後の制御時点から
次のパルスの到着直後の制御時点までの間隔T1 と、そ
の間の平均速度<ω[j+1]>を用いて、時点[j+
1,0]における位置が、 θ[j+1,0]=θ[j,0]+<ω[j+1]>T1 (14) で近似できる。この発明では、パルス発生時点と制御時
点の時刻のずれは、大きくてT2 程度なので無視できる
と仮定している。
Then, using the interval T 1 from the control time immediately after the arrival of the pulse to the control time immediately after the arrival of the next pulse and the average speed <ω [j + 1]> therebetween, the time [j +
[1,0] can be approximated by θ [j + 1,0] = θ [j, 0] + <ω [j + 1]> T 1 (14). In the present invention, the deviation of the time of the control point and the pulse generation timing is assumed to be negligible because large T 2 degree.

【0045】(13)、(14)式より、パルスが発生した直後
の制御時点での真の回転位置θがわかるので、以後、誤
差修正等のアルゴリズムはT1 が変動する以外は前記第
1の態様の位置読み込み形の瞬時速度オブザーバ10と
全く同じである。
[0045] (13), (14) from the equation, the pulse is found true rotation position θ at the control point immediately after the occurrence, thereafter, the first except algorithm such as error correction is that T 1 is varied This is exactly the same as the position reading type instantaneous speed observer 10 of the embodiment.

【0046】ところで、車両の車輪速度を検出するに
は、現状、以下のような手法が用いられている(例え
ば、特開平6−92116)。すなわち、車軸に取り付
けられた、車輪速に比例した周波数成分を持つ交流信号
を発生する車輪速センサからの出力波形をパルス波形に
整形し、図8に示すように、このパルス信号を用いて、
決められた制御周期Ts に基づくある読み込み区間kで
発生したパルス個数N(k)及びサンプルタイミング直
前の立ち上がりパルス間隔Tw (k)を求める。そし
て、次式により区間kでの車輪速度ω(k)を求める。
Incidentally, the following method is currently used to detect the wheel speed of a vehicle (for example, Japanese Patent Laid-Open No. Hei 6-92116). That is, the output waveform from the wheel speed sensor that generates an AC signal having a frequency component proportional to the wheel speed attached to the axle is shaped into a pulse waveform, and as shown in FIG.
Request-determined control period T s pulse number occurred in certain loading period k based on N (k) and sample timing immediately before the rising pulse interval T w (k). Then, the wheel speed ω (k) in the section k is obtained by the following equation.

【0047】ω(k)=Kw ・N(k)/Tw (k) なお、Kw は変換係数である。Ω (k) = K w · N (k) / T w (k) where K w is a conversion coefficient.

【0048】[0048]

【発明が解決しようとする課題】しかしながら、上記従
来技術には、以下のような問題がある。
However, the above prior art has the following problems.

【0049】まず、上記発明のうち、瞬時速度オブザー
バの第1の態様では、位置検出器から発生するパルスの
間隔は制御周期T2 より短いものと仮定しているため、
この仮定が成立しないような低速度域では適用すること
ができない。
First, in the first aspect of the instantaneous speed observer in the above invention, since the interval between pulses generated from the position detector is assumed to be shorter than the control period T 2 ,
It cannot be applied in a low speed range where this assumption does not hold.

【0050】また、瞬時速度オブザーバの第2の態様で
は、低速度域に適用可能であるが、パルス発生時点と制
御時点との時刻のずれは、大きくてT2 程度なので無視
できると考えている。これは、DSPのような高速演算
デバイスを用いた場合に成り立つ条件であり、本発明で
は、T2 は数十[μs]程度という極めて速い周期を仮
定している。しかし、このような高速演算デバイスは一
般に高価であるため、装置全体が高額になるという問題
がある。さらに、安価な低速演算デバイスを用いると、
制御周期T2 が長くなり、パルス発生時点と制御時点の
時刻のずれは無視できなくなるため、真の検出位置に基
づく誤差の修正が難しくなる。
[0050] In the second aspect of the instantaneous speed observer, but is applicable to a low speed range, displacement of the time between the control point and the pulse generation time is considered negligible since large T 2 degree . This is a condition that is satisfied when a high-speed operation device such as a DSP is used. In the present invention, it is assumed that T 2 has an extremely fast cycle of about several tens [μs]. However, since such a high-speed arithmetic device is generally expensive, there is a problem that the entire apparatus becomes expensive. Furthermore, if an inexpensive low-speed arithmetic device is used,
Control period T 2 becomes longer, since the displacement time of the control point and the pulse generation time can not be ignored, the error correction is difficult based on the true position detection.

【0051】つまり、図9に示すように、パルス発生時
点と制御時点(k=m)との時刻のずれΔtが大きくな
り、その制御時点において真の位置θ* と(13)、(14)式
に基づいて計算された近似位置θ[j+1,0]との間
の差、すなわち量子化誤差が無視できなくなるほど大き
くなる。これに伴って、量子化ノイズの増大、パラメー
タ変動に対するロバスト性の劣化、オブザーバの不安定
化といった不具合が生じる。さらに、第2の態様では、
パルス発生毎に誤差修正を行うため、パルスが多く発生
する高速度域では、誤差修正の演算が増加する。
That is, as shown in FIG. 9, the time difference Δt between the pulse generation time and the control time (k = m) increases, and the true position θ * and (13), (14) at the control time. The difference from the approximate position θ [j + 1,0] calculated based on the equation, that is, the quantization error becomes so large that it cannot be ignored. Along with this, problems such as an increase in quantization noise, deterioration of robustness against parameter fluctuation, and instability of the observer occur. Further, in the second aspect,
Since the error correction is performed every time a pulse is generated, the calculation of the error correction increases in a high speed region where many pulses are generated.

【0052】一方、従来の車輪速度検出方法では、制御
周期間の平均速度を検出するものであって、その読み取
り時点での瞬時速度ではない。また、パルス発生時点と
読み取り時点との間にずれがあるため、それに起因する
検出遅れ及び量子化ノイズの増大を招いている。さらに
は、低速度域では、制御周期間でパルスが1個も入らな
い状態が存在するため、こうした問題はより顕著とな
る。車輪速度は、アンチスキッド装置のような車両制御
装置に用いられるが、本来こうした制御は、瞬時制御を
用いるべきであり、従来法による速度検出では、検出速
度の信頼性からより高度なきめの細かい車両制御の実現
を阻害しているという問題がある。
On the other hand, the conventional wheel speed detecting method detects the average speed during the control cycle, and does not indicate the instantaneous speed at the time of reading. In addition, since there is a difference between the pulse generation time and the reading time, the detection delay and quantization noise increase due to the difference. Further, in the low speed region, such a problem becomes more remarkable because there is a state in which no pulse enters between control periods. Wheel speed is used in a vehicle control device such as an anti-skid device, but such control should be performed by using instantaneous control, and the speed detection by the conventional method requires a higher level of fineness from the reliability of the detected speed. There is a problem that the realization of vehicle control is hindered.

【0053】より高度な車両制御を実現するためには、
タイヤと路面間に発生する力を知る必要がある。制動
時、この発生力は制動力と呼ばれており、次式で表せ
る。
In order to realize more advanced vehicle control,
It is necessary to know the force generated between the tire and the road. At the time of braking, this generated force is called a braking force and can be expressed by the following equation.

【0054】[0054]

【数11】 [Equation 11]

【0055】・ここに、Jは車輪の慣性モーメント、ω
は車輪加速度、Tb は制動トルク、Rはタイヤの動荷重
半径である。(15)式からも明らかなように、制動力を求
めるためには、車輪加速度が必要である。従来技術で
は、車輪速度の時間差分によって車輪加速度を算出する
が、差分という操作は、ノイズに対して敏感であるた
め、検出遅れや量子化ノイズを含んでいる車輪速度から
車輪加速度を正確に求めるのは事実上不可能である。従
って、制動力も正しい値を得ることができない。
Where J is the moment of inertia of the wheel, ω
Is the wheel acceleration, Tb is the braking torque, and R is the dynamic load radius of the tire. As is clear from equation (15), wheel acceleration is required to determine the braking force. In the related art, the wheel acceleration is calculated based on the time difference between the wheel speeds. Since the operation of the difference is sensitive to noise, the wheel acceleration is accurately obtained from the wheel speed including the detection delay and the quantization noise. Is virtually impossible. Therefore, a correct value cannot be obtained for the braking force.

【0056】車輪の瞬時速度、加速度を検出する従来手
段として、車輪速センサが出力した車輪速に比例した周
波数成分をもつ交流信号をアナログ回路を用いて車輪速
度に変換し、アンチスキッド装置に用いるものがある
(特開昭56−34551号公報)。しかし、センサ信
号の処理を含む制御装置全体がデジタルシステムとして
設計される現在においては、アナログ回路の追加はコス
トアップにつながる。さらに、車輪速センサからの交流
信号は正確な正弦波になっておらず、また、多くのノイ
ズ成分を含んでいるため、正しい瞬時速度、加速度が検
出されないという問題がある。加えて、この従来技術で
は、車輪速に比例した周波数成分をもつ交流信号を発生
する車輪速センサにしか適用できず、回転角に応じてパ
ルス信号などの非正弦波信号が発生するセンサには対応
できないという問題がある。
As a conventional means for detecting the instantaneous speed and acceleration of a wheel, an AC signal output from a wheel speed sensor and having a frequency component proportional to the wheel speed is converted into a wheel speed using an analog circuit and used for an anti-skid device. (JP-A-56-34551). However, at present, when the entire control device including the processing of the sensor signal is designed as a digital system, the addition of an analog circuit leads to an increase in cost. Furthermore, since the AC signal from the wheel speed sensor is not an accurate sine wave and contains many noise components, there is a problem that a correct instantaneous speed and acceleration cannot be detected. In addition, this conventional technique is applicable only to a wheel speed sensor that generates an AC signal having a frequency component proportional to the wheel speed, and is not applicable to a sensor that generates a non-sine wave signal such as a pulse signal according to the rotation angle. There is a problem that cannot be dealt with.

【0057】本発明は、上記事実に鑑みて成されたもの
で、低速度域においても、また、制御系の制御周期が長
く、パルス発生時点と制御時点の時刻のずれが無視でき
ないような状況においても、回転体の瞬時速度を精度良
く推定できる回転速度検出装置及び該装置を用いること
により精度良く制動力を推定できるタイヤ発生力検出装
置を提供することを目的とする。
The present invention has been made in view of the above-mentioned circumstances. Even in a low speed range, the control cycle of the control system is long, and the time lag between the pulse generation time and the control time cannot be ignored. It is another object of the present invention to provide a rotation speed detection device capable of accurately estimating the instantaneous speed of a rotating body and a tire generation force detection device capable of accurately estimating a braking force by using the rotation speed detection device.

【0058】[0058]

【課題を解決するための手段】上記目的を達成するため
に請求項1の発明は、回転体に取り付けられたパルス生
成器が出力したパルスの計数値に基づいて回転体の回転
位置を検出する位置検出手段と、回転体を制御するため
の制御トルクを検出する制御トルク検出手段と、前記制
御トルク検出手段により検出された値の制御トルク及び
推定された値の外乱トルクが作用する回転体のモデルに
基づいて、与えられた回転速度及び回転位置の各初期値
とから回転体の回転速度及び回転位置を推定する状態推
定手段と、少なくとも前記パルス生成器が出力したパル
スの数に基づいて、誤差修正の時点を判断する誤差修正
判断手段と、前記誤差修正判断手段により判断された誤
差修正の時点において、前記位置検出手段により検出さ
れた回転位置と前記状態推定手段により推定された回転
位置との間に生じた位置の誤差を、外乱トルクの推定値
の誤差、回転速度の初期値の誤差、及び検出された回転
位置に含まれる量子化誤差によるものとみなし、前記位
置の誤差に基づいて、前記状態推定手段が用いる外乱ト
ルクの推定値と回転速度及び回転位置の各初期値とを修
正する推定誤差修正手段と、を含んで構成したものであ
る。
According to a first aspect of the present invention, a rotational position of a rotating body is detected based on a count value of pulses output by a pulse generator attached to the rotating body. Position detecting means, control torque detecting means for detecting control torque for controlling the rotating body, and a rotating body on which the control torque of the value detected by the control torque detecting means and the disturbance torque of the estimated value act. Based on the model, state estimation means for estimating the rotational speed and rotational position of the rotating body from the given initial value of the rotational speed and rotational position, based on at least the number of pulses output by the pulse generator, Error correction determining means for determining a time point of error correction; and a rotational position detected by the position detecting means at a time point of error correction determined by the error correction determining means. The position error generated between the rotational position estimated by the state estimating means is determined by the error of the disturbance torque estimated value, the error of the initial value of the rotational speed, and the quantization error included in the detected rotational position. And an estimation error correction unit that corrects the estimated value of the disturbance torque used by the state estimation unit and the initial values of the rotation speed and the rotation position based on the position error. .

【0059】請求項1の発明では、回転体が回転する
と、この回転体に取り付けられたパルス生成器がパルス
を出力し、位置検出手段は、出力されたパルスの計数値
に基づいて回転体の回転位置を検出する。例えば、この
パルス生成器を、回転速度に比例した周波数成分を有す
る交流信号をパルス波形に変換するパルスジェネレータ
として実現することができる。また、位置検出手段を、
パルス発生回数を積算することによりパルスの計数値を
求め、このパルス計数値を回転位置に変換する手段とし
て実現することができる。なお、パルス生成器は、上記
のような回転速度に比例した周波数成分をもつ交流信号
を発生するものに限らず、回転角に応じてパルス等の何
らかの信号を発生するものでも良い。
According to the first aspect of the present invention, when the rotating body rotates, the pulse generator attached to the rotating body outputs a pulse, and the position detecting means detects the rotation of the rotating body based on the output pulse count value. Detect the rotational position. For example, this pulse generator can be realized as a pulse generator that converts an AC signal having a frequency component proportional to the rotation speed into a pulse waveform. Also, the position detecting means
This can be realized as a means for calculating the pulse count value by integrating the number of times of pulse generation, and converting the pulse count value into a rotational position. The pulse generator is not limited to generating an AC signal having a frequency component proportional to the rotation speed as described above, but may generate a signal such as a pulse according to the rotation angle.

【0060】そして、制御トルク検出手段が、回転体を
制御するための制御トルクを検出する。この制御トルク
は、回転体が車両の車輪の場合はホイールシリンダへの
ブレーキ圧による車輪への制動トルクであり、電動機の
場合はモータを駆動するための駆動トルクに相当する。
Then, the control torque detecting means detects a control torque for controlling the rotating body. When the rotating body is a vehicle wheel, the control torque is a braking torque applied to the wheel by a brake pressure applied to a wheel cylinder, and corresponds to a driving torque for driving a motor when the rotating body is an electric motor.

【0061】また、状態推定手段が、制御トルク検出手
段により検出された値の制御トルク及び推定された値の
外乱トルクが作用する回転体のモデルに基づいて、与え
られた回転速度及び回転位置の各初期値とから回転体の
回転速度及び回転位置を推定する。例えば、制御トルク
と外乱トルクとを加算して回転体に作用する全トルクを
求め、この全トルクが作用した場合の回転体の回転角加
速度を、回転体のモデルに含まれる物理量(慣性モーメ
ント)から求めることができる。そして、回転角加速度
と回転速度の初期値とから積分演算によりこの時点の回
転体の回転速度を求め、さらに、この回転速度と回転位
置の初期値とから積分演算によりこの時点の回転位置を
求めることができる。なお、外乱トルクの推定値と回転
速度及び回転位置の各初期値は、最初は、所定の測定値
や推定値により得られた値を用いても良いが、後述する
推定誤差修正手段により修正された場合、状態推定手段
は、修正された外乱トルク、回転速度及び回転位置の各
初期値に基づいて回転体の回転速度及び回転位置を推定
する。
Further, the state estimating means detects the given rotational speed and rotational position based on the model of the rotating body on which the control torque of the value detected by the control torque detecting means and the disturbance torque of the estimated value act. The rotation speed and the rotation position of the rotating body are estimated from each of the initial values. For example, the total torque acting on the rotating body is obtained by adding the control torque and the disturbance torque, and the rotational angular acceleration of the rotating body when the full torque acts is calculated as a physical quantity (moment of inertia) included in the model of the rotating body. Can be obtained from Then, the rotational speed of the rotating body at this time is obtained from the rotational angular acceleration and the initial value of the rotational speed by an integral operation, and the rotational position at this time is obtained from the rotational speed and the initial value of the rotational position by an integral operation. be able to. As the estimated value of the disturbance torque and each initial value of the rotational speed and the rotational position, initially, a value obtained from a predetermined measured value or an estimated value may be used, but the initial value is corrected by an estimation error correction unit described later. In this case, the state estimating means estimates the rotation speed and the rotation position of the rotating body based on the corrected initial values of the disturbance torque, the rotation speed, and the rotation position.

【0062】このように状態推定手段が回転速度等を推
定している間、誤差修正手段が、少なくともパルス生成
器が出力したパルスの数に基づいて、誤差修正の時点を
判断する。例えば、高速カウンタによりパルスが所定回
数カウントされた時点(この時点で高速カウンタはリセ
ットされる)を誤差修正の時点と判断しても良い(上記
所定回数が1回の場合は上記従来の第2の態様に相当す
る)。また、予め定められた時間間隔におけるパルスの
計数値の変化量が所定値以上となった時点を、誤差修正
時点と判断するようにしても良い。
As described above, while the state estimating means is estimating the rotational speed or the like, the error correcting means determines the time of error correction based on at least the number of pulses output from the pulse generator. For example, the point in time when the pulse is counted a predetermined number of times by the high-speed counter (the high-speed counter is reset at this point) may be determined to be the point in time of error correction. Corresponding to the embodiment). Further, a point in time at which the amount of change in the pulse count value at a predetermined time interval becomes equal to or greater than a predetermined value may be determined as an error correction point.

【0063】従来では、上記位置の誤差の原因として、
外乱トルクの推定値の誤差及び回転速度の初期値の誤差
のみを考え、位置の誤差に基づいて外乱トルクの推定値
及び回転速度の初期値を修正し、さらに回転位置の初期
値を、誤差修正の判断時点(第2の態様では、パルス発
生時点)で検出された回転位置に修正する。しかし、実
際に修正された初期値により状態の指定が実行される時
点は、誤差修正の判断時点の直後の制御時点であり、両
時点間には時間差があるので、回転位置の検出値と制御
時点での真の位置との間に差、いわゆる量子化誤差が発
生する。
Conventionally, as a cause of the position error,
Considering only the error of the estimated value of the disturbance torque and the error of the initial value of the rotational speed, correct the estimated value of the disturbance torque and the initial value of the rotational speed based on the position error, and further correct the error of the initial value of the rotational position. (In the second mode, at the time of pulse generation). However, the time point at which the state is specified by the actually corrected initial value is the control time point immediately after the error correction determination time point, and there is a time difference between the two time points. A difference between the true position at the time and a so-called quantization error occurs.

【0064】そこで、本発明では、誤差修正手段により
誤差修正の時点が判断されたとき、推定誤差修正手段
が、この時点において、位置検出手段により検出された
回転位置と状態推定手段により推定された回転位置との
間に生じた位置の誤差を、外乱トルクの推定値の誤差、
回転速度の初期値の誤差、及び回転位置の検出値に含ま
れる量子化誤差によるものとみなし、位置の誤差に基づ
いて、外乱トルクの推定値と回転速度及び回転位置の各
初期値とを修正する。修正された値は、状態推定手段で
用いられるので、各指定値は真の値に接近していく。な
お、位置の誤差を各物理量の誤差に分配する方法は、例
えば、推定誤差が時間と共に零に収束するように制御系
を設計する際の拘束条件などから得られる。
Therefore, according to the present invention, when the time of error correction is determined by the error correcting means, the estimated error correcting means estimates the rotational position detected by the position detecting means and the state estimating means at this time. The error of the position generated between the rotational position and the error of the estimated value of the disturbance torque,
Assuming that it is due to the error of the initial value of the rotational speed and the quantization error included in the detected value of the rotational position, the estimated value of the disturbance torque and the initial values of the rotational speed and the rotational position are corrected based on the position error I do. Since the corrected values are used by the state estimating means, each designated value approaches a true value. The method of distributing the position error to the error of each physical quantity can be obtained, for example, from a constraint condition when designing a control system such that the estimation error converges to zero with time.

【0065】このように本発明では、外乱トルクの推定
値と回転速度の初期値の誤差と共に回転位置の初期値に
含まれる量子化誤差を考慮し、位置の誤差に基づいて、
外乱トルクの推定値、回転速度の初期値及び回転位置の
初期値を修正するようにしたので、従来技術に比べて推
定速度を速く設定でき、かつ推定誤差も小さく抑えるこ
とができる。特に、演算デバイスの能力によって制御系
の制御周期が十分に短くできず、よって、回転位置の量
子化誤差が増大する場合においても、正確に回転速度を
推定することができると共に、量子化ノイズの増大、パ
ラメータ変動に対するロバスト性の劣化、及びオブザー
バの不安定化などを防止することができる。
As described above, in the present invention, the quantization error included in the initial value of the rotational position together with the error between the estimated value of the disturbance torque and the initial value of the rotational speed is considered, and based on the position error,
Since the estimated value of the disturbance torque, the initial value of the rotational speed, and the initial value of the rotational position are corrected, the estimated speed can be set faster than in the related art, and the estimation error can be reduced. In particular, even when the control cycle of the control system cannot be sufficiently shortened due to the capability of the arithmetic device, and thus the quantization error of the rotation position increases, the rotation speed can be accurately estimated and the quantization noise can be reduced. It is possible to prevent an increase, a deterioration in robustness against parameter fluctuation, and an unstable observer.

【0066】さらに、誤差修正の時点をパルス数に基づ
いて決定することにより、制御周期でパルスが発生しな
いような低速度域においても、正確に回転速度を推定す
ることができる。
Further, by determining the time of error correction based on the number of pulses, it is possible to accurately estimate the rotational speed even in a low speed range where no pulse is generated in the control cycle.

【0067】また、本発明は、請求項1の前記誤差修正
判断手段を、制御系の制御周期が一定回数以上繰り返さ
れる間に前記パルス生成器が出力したパルスの数が基準
値以上となった時点を誤差修正の時点と判断することも
できる。
Further, according to the present invention, the error correction judging means according to claim 1, wherein the number of pulses output by the pulse generator during the control cycle of the control system is repeated a predetermined number of times or more becomes equal to or more than a reference value. The time point may be determined to be the time point of error correction.

【0068】これによって、パルスの数の上記基準値を
適当な値に設定すれば、パルス発生毎に推定誤差の修正
を行う従来技術に見られる高速度域での演算負荷の増大
を回避することができる。
Thus, by setting the reference value of the number of pulses to an appropriate value, it is possible to avoid an increase in calculation load in a high-speed region, which is seen in the prior art in which an estimation error is corrected every time a pulse is generated. Can be.

【0069】請求項2の発明は、車輪に取り付けられた
パルス生成器が出力したパルスの計数値に基づいて車輪
の回転位置を検出する位置検出手段と、車輪の制動トル
クを検出する制動トルク検出手段と、前記制動トルク検
出手段により検出された値の制動トルク及び推定された
値の外乱トルクが作用する車輪のモデルに基づいて、与
えられた回転速度及び回転位置の各初期値とから車輪の
回転速度及び回転位置を推定する状態推定手段と、少な
くとも前記パルス生成器が出力したパルスの数に基づい
て、誤差修正の時点を判断する誤差修正判断手段と、前
記誤差修正判断手段により判断された誤差修正の時点に
おいて、前記位置検出手段により検出された回転位置と
前記状態推定手段により推定された回転位置との間に生
じた位置の誤差を、外乱トルクの推定値の誤差、回転速
度の初期値の誤差、及び検出された回転位置に含まれる
量子化誤差によるものとみなし、前記位置の誤差に基づ
いて、前記状態推定手段が用いる外乱トルクの推定値と
回転速度及び回転位置の各初期値とを修正する推定誤差
修正手段と、前記推定誤差修正手段により修正された外
乱トルクと車輪の動半径とに基づいて、路面から反力と
して車輪に作用する制動力を推定するタイヤ発生力推定
手段と、を含んで構成したものである。
According to a second aspect of the present invention, there is provided a position detecting means for detecting a rotational position of a wheel based on a count value of pulses output from a pulse generator attached to a wheel, and a braking torque detecting means for detecting a braking torque of the wheel. Means, based on a model of the wheel on which the braking torque detected by the braking torque detection means and the estimated disturbance torque act, and from the respective initial values of the given rotation speed and rotation position, State estimation means for estimating a rotation speed and a rotation position, error correction judgment means for judging a time point of error correction based on at least the number of pulses output by the pulse generator, and judgment by the error correction judgment means At the time of error correction, a position error generated between the rotation position detected by the position detection means and the rotation position estimated by the state estimation means is calculated. The error of the estimated value of the disturbance torque, the error of the initial value of the rotational speed, and the quantization error included in the detected rotational position are considered to be due to the error of the position. Estimation error correction means for correcting the estimated value and each initial value of the rotation speed and the rotation position, and the disturbance torque and the moving radius of the wheel corrected by the estimation error correction means, as a reaction force from the road surface to the wheel. Tire generating force estimating means for estimating the acting braking force.

【0070】請求項2の発明では、タイヤ発生力推定手
段が、推定誤差修正手段により修正された外乱トルクと
車輪の動半径とに基づいて、路面から反力として車輪に
作用する制動力を推定する。本発明では、請求項1の回
転速度検出装置と同じ原理で、制御周期にパルスが発生
しないような低速度域及び高速度域のいずれにおいて
も、正確かつ安定に外乱トルクを推定できるので、常に
良好に制動力を推定することができる。また、従来技術
にあるような車輪速度の差分を用いることがないので、
種々のノイズに強く、安定して制動力を推定することが
できる。
According to the second aspect of the present invention, the tire generating force estimating means estimates the braking force acting on the wheels as a reaction force from the road surface based on the disturbance torque corrected by the estimation error correcting means and the moving radius of the wheels. I do. According to the present invention, the disturbance torque can be accurately and stably estimated in both the low speed range and the high speed range in which no pulse is generated in the control cycle, based on the same principle as the rotation speed detection device of claim 1. The braking force can be satisfactorily estimated. Also, since there is no need to use the difference in wheel speed as in the prior art,
It is resistant to various noises and can stably estimate the braking force.

【0071】[0071]

【発明の実施の形態】以下、本発明の実施の形態に係る
回転速度検出装置及びタイヤ発生力検出装置を図面に基
づいて各々詳細に説明する。 (回転速度検出装置)本発明の実施の形態に係る回転速
度検出装置の構成ブロック図を図1に示す。同図に示す
ように、本回転速度検出装置は、計測対象である回転体
を制御するために加えられる制御トルクを検出する制御
トルク検出手段30と、回転体の力学的モデルに基づい
て、検出された制御トルク及び推定された外乱トルクに
より回転体の瞬時位置及び瞬時速度を推定する状態推定
手段32と、回転体の所定の回転角の回転毎にパルスを
生成するパルスジェネレータ37を備え、該パルスジェ
ネレータが出力したパルスの計数値に基づいて回転体の
位置を検出する位置検出手段38と、を含んで構成され
ている。なお、制御トルク検出手段30及び状態推定手
段32は、従来の瞬時速度オブザーバ10の等価モデル
16(図5)により実現できるが、この例に限定される
ものではない。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, a rotational speed detecting device and a tire generating force detecting device according to an embodiment of the present invention will be described in detail with reference to the drawings. (Rotation Speed Detecting Apparatus) FIG. 1 is a block diagram showing the configuration of the rotation speed detecting apparatus according to the embodiment of the present invention. As shown in the figure, the present rotational speed detecting device detects a control torque based on a dynamic model of a rotating body, and a control torque detecting means 30 for detecting a control torque applied to control a rotating body to be measured. A state estimating means 32 for estimating the instantaneous position and instantaneous speed of the rotating body based on the obtained control torque and the estimated disturbance torque; and a pulse generator 37 for generating a pulse for each rotation of the rotating body at a predetermined rotation angle. And position detecting means 38 for detecting the position of the rotating body based on the count value of the pulse output from the pulse generator. The control torque detecting means 30 and the state estimating means 32 can be realized by the equivalent model 16 (FIG. 5) of the conventional instantaneous speed observer 10, but the present invention is not limited to this example.

【0072】さらに、本回転速度検出装置は、パルスジ
ェネレータ37からのパルス数と制御周期T2 毎にカウ
ントされるカウント数とに基づいて誤差修正演算のタイ
ミングを判断する誤差修正判断手段34と、誤差修正演
算のタイミングと判断されたときに、状態推定手段32
により推定された位置と位置検出手段38により検出さ
れた位置との差分に基づいて、速度、外乱トルク及び位
置の推定値の修正を行う推定誤差修正手段36と、から
構成される。
Further, the present rotational speed detecting device includes an error correction judging means 34 for judging an error correction calculation timing based on the number of pulses from the pulse generator 37 and the count number counted in each control cycle T 2 . When it is determined that the timing of the error correction operation has been reached, the state estimating means 32
And an estimation error correction unit 36 that corrects the estimated values of the speed, the disturbance torque, and the position based on the difference between the position estimated by the above and the position detected by the position detection unit 38.

【0073】次に、図1の回転速度検出装置の作用を図
3のフローチャートに基づいて説明する。なお、図3の
フローチャートにおいて、繰り返しループは、制御周期
2のタイミングで繰り返されるものであり、定数M、
Nは、予め定められた零以上の定数である。また、jは
誤差修正タイミングのインデックス、kは制御周期T 2
の制御タイミングのインデックスであり、本実施の形態
でも、上述のように時刻tを[j,k]で表すこととす
る。
Next, the operation of the rotation speed detecting device of FIG. 1 will be described.
3 will be described. In addition, FIG.
In the flowchart, the repetition loop is a control cycle.
TTwoIs repeated at the timing of
N is a predetermined constant of zero or more. And j is
Index of error correction timing, k is control cycle T Two
Index of the control timing of the present embodiment.
However, the time t is represented by [j, k] as described above.
You.

【0074】図3のフローチャートに示すように、ま
ず、ステップ99において各変数の初期化を行い、次
に、制御周期T2 毎に繰り返される繰り返しループに移
行する。このループの始めにステップ100においてk
を1つカウントアップする。
[0074] As shown in the flowchart of FIG. 3, first, initializes the respective variables in step 99, then proceeds to repeat the loop which is repeated for each control period T 2. At the beginning of this loop, in step 100, k
Is counted up by one.

【0075】次のステップ101では、時点[j,k]
における回転体の制御トルクTq [j,k]を検出す
る。制御トルクは、回転体を制御するために加えられる
トルクであり、回転体が電動機の電動子の場合は、時点
[j,k]の電動子電流をセンサ等で検出し検出値に制
御トルクへの変換係数を乗ずることによって得られる。
また、回転体が車両の車輪の場合は、ホイールシリンダ
圧を検出し、検出値にブレーキトルクへの変換係数を乗
ずることによって得られる。そして、次のステップ10
2では、誤差修正判断手段34が、時点[j,k−1]
と[j,k]との間でパルスジェネレータ37から発生
したパルス数np を読み込む。
At the next step 101, the time [j, k]
, The control torque T q [j, k] of the rotating body is detected. The control torque is a torque applied to control the rotating body. When the rotating body is an electric motor, the motor current at the point [j, k] is detected by a sensor or the like, and the detected value is converted to the control torque. Is obtained by multiplying by the conversion coefficient of
When the rotating body is a wheel of a vehicle, it is obtained by detecting a wheel cylinder pressure and multiplying the detected value by a conversion coefficient into a brake torque. And the next step 10
In 2, the error correction determining means 34 determines whether the time [j, k-1]
The number of pulses n p generated from the pulse generator 37 between the time and [j, k] is read.

【0076】ところで、既に述べたように、回転体の角
加速度は、制御トルクと回転体に加わる外乱Td との総
和トルクを回転体の慣性モーメントで割ることによって
求まり、回転体の角速度は、この角速度を時間積分する
ことによって求まり、回転体の回転位置は、さらに角速
度を時間積分することによって求まる。
As described above, the angular acceleration of the rotating body is obtained by dividing the total torque of the control torque and the disturbance Td applied to the rotating body by the moment of inertia of the rotating body. The angular velocity is obtained by time integration, and the rotational position of the rotating body is obtained by further integrating the angular velocity over time.

【0077】そこで、状態推定手段32が、ステップ1
03において、ステップ101で得られた制御トルクと
外乱トルクの推定値とから回転体に加わる総和トルクを
(4)式により演算する。なお、外乱トルクの推定値は、
後述するように、誤差修正演算により修正されて真の外
乱トルクに接近するが、誤差修正演算を初めて実行する
前の外乱トルクの初期値は、検出された値を用いても良
い。例えば、外乱トルクとして車輪に作用する制動力を
推定する場合、上記(15)式に基づいて検出された制動力
を初期値として用いる。
Then, the state estimating means 32 executes step 1
In 03, the total torque applied to the rotating body is calculated from the control torque obtained in step 101 and the estimated value of the disturbance torque.
Calculate using equation (4). The estimated value of the disturbance torque is
As will be described later, the correction value is corrected by the error correction calculation to approach the true disturbance torque, but the detected value may be used as the initial value of the disturbance torque before the error correction calculation is performed for the first time. For example, when estimating the braking force acting on the wheels as the disturbance torque, the braking force detected based on the above equation (15) is used as the initial value.

【0078】次のステップ105では、状態推定手段
が、上記の(5) 、(6) 式により台形積分を用いて回転体
の速度の推定値及び回転位置の推定値を演算する。
In the next step 105, the state estimating means calculates the estimated value of the speed of the rotating body and the estimated value of the rotational position using the trapezoidal integration by the above equations (5) and (6).

【0079】次のステップ106では、誤差修正判断手
段34が、時点[j,0]から時点[j,k]までの間
にパルスジェネレータ37が出力したパルス数の総和n
を演算する。この演算では、ステップ102で読み込ま
れたパルス数npを、パルス数nに累積的に加算する
(n←n+np)。そして、ステップ107及びステッ
プ108において誤差修正判断を行う。
In the next step 106, the error correction judging means 34 calculates the sum n of the number of pulses output by the pulse generator 37 from the time point [j, 0] to the time point [j, k].
Is calculated. In this calculation, the pulse number np read in step 102 is cumulatively added to the pulse number n (n ← n + np). Then, in steps 107 and 108, an error correction determination is made.

【0080】すなわち、ステップ107では、インデッ
クスkが規定値M以上となったか否かを判定し、否定判
定であれば、再びステップ100から同様の演算を繰り
返す。演算の繰り返しによりkが増加し、やがてステッ
プ107で肯定判定となったとき、ステップ108にお
いて、ステップ106で加算されるパルス数の総和nが
規定値N以上であるか否かを判定する。ステップ108
が否定判定であれば、再びステップ100から同様の演
算を繰り返す。このループを繰り返し演算中にパルス数
が増加し、やがてステップ108で肯定判定となったと
き、誤差修正演算を行う。
That is, in step 107, it is determined whether or not the index k has become equal to or greater than the specified value M. If the determination is negative, the same calculation is repeated from step 100 again. When k is increased by the repetition of the calculation, and a positive determination is made in step 107, it is determined in step 108 whether or not the total number n of the pulses added in step 106 is equal to or greater than a specified value N. Step 108
Is negative, the same calculation is repeated from step 100 again. When the number of pulses increases during the repetitive calculation of this loop and a positive determination is made in step 108, an error correction calculation is performed.

【0081】すなわち、最初1にセットされたインデッ
クスkがM以上となるまでの間にカウントされたパルス
数の総和nがN以上となったときに誤差修正演算を行
う。なお、上記の規定値M、Nの値は、制御周期T2
長さやパルスジェネレータ37の分解能などに応じて任
意好適に定められる。
That is, the error correction operation is performed when the total number n of pulses counted until the index k initially set to 1 becomes M or more becomes N or more. The above predetermined value M, the value of N is determined arbitrarily suitably depending on the resolution of the length and pulse generator 37 of the control period T 2.

【0082】誤差修正演算では、以下のような演算を行
う。上記処理までで加算されたkの値の時点で誤差修正
を行うため、その定義より誤差修正周期T1 は制御周期
2のk倍となる。そこで、ステップ109において、 T1 = kT2 (16) より、誤差修正期間T1 を求める。なお、T1 は出力さ
れるパルス数の発生頻度に応じて変動する。そして、こ
のタイミングでの位置情報を、ステップ110におい
て、位置検出手段38が、 θ[j+1,0]=θ[j,0]+nθs (17) より演算する。ここに、θs は位置検出器から発生する
1 パルス当たりの回転角である。
In the error correction calculation, the following calculation is performed. Since the error correction is performed at the time of the value of k added up to the above processing, the error correction cycle T 1 is k times the control cycle T 2 according to the definition. Therefore, in step 109, from T 1 = kT 2 (16) , determining the error correction period T 1. Note that T 1 varies according to the frequency of occurrence of the number of pulses to be output. Then, in step 110, the position detecting means 38 calculates the position information at this timing from θ [j + 1,0] = θ [j, 0] + nθ s (17). Where θ s is generated from the position detector
The rotation angle per pulse.

【0083】次に、ステップ111において、(17)式に
より得られた位置θ[j+1,0]がθ1 以上となった
か否かを判定する。ここに、θ1 は回転体の1回転分の
回転角であり、ラジアンの単位で表現すれば、θ1 =2
π[rad]となる。すなわち、ステップ111は、回
転体が1回転したか否かの判定となる。
[0083] Next, in step 111, whether it is the position θ [j + 1,0] is theta 1 or obtained by the equation (17). Here, θ 1 is a rotation angle for one rotation of the rotating body, and expressed in units of radians, θ 1 = 2
π [rad]. That is, step 111 is to determine whether the rotating body has made one rotation.

【0084】ステップ111で、位置θ[j+1,0]
がθ1 以上であると肯定判定した場合には、ステップ1
12で、θ[j+1,0]及びθ[j,k]の推定値か
らそれぞれθ1 を減じ、次のステップ113に移行す
る。この処理を施すことによって、位置は常に
At step 111, the position θ [j + 1,0]
If a positive determination is made that θ 1 or more,
At 12, θ 1 is subtracted from the estimated values of θ [j + 1,0] and θ [j, k], and the process proceeds to the next step 113. By performing this process, the position is always

【0085】[0085]

【数12】 (Equation 12)

【0086】の範囲に入るため、θ[j+1,0],θ
[j,k]の推定値(∧付)のオーバーフローを防ぐこ
とができる。なお、ステップ111で、否定判定の場合
には、ステップ112の処理を実行することなく直ちに
ステップ113に移行する。
[0108] θ [j + 1, 0], θ
It is possible to prevent an overflow of the estimated value (appended) of [j, k]. If a negative determination is made in step 111, the process immediately proceeds to step 113 without executing the process of step 112.

【0087】次のステップ113では、ステップ105
で求められた位置の推定値と、ステップ110で求めら
れた位置検出手段からの位置情報との差分を次式のよう
に求める。
In the next step 113, step 105
The difference between the estimated position value obtained in step (1) and the position information from the position detecting means obtained in step 110 is obtained as in the following equation.

【0088】[0088]

【数13】 (Equation 13)

【0089】次のステップ114では、推定誤差修正手
段36が、ステップ113で求められた差分Δθに基づ
いて、次の時点[j+1,0]からの位置推定に行われ
る位置、速度、外乱のそれぞれの推定値の初期値を、次
式のように修正する。
In the next step 114, the estimation error correcting means 36 determines the position, velocity, and disturbance to be used for position estimation from the next time point [j + 1,0] based on the difference Δθ obtained in step 113. Is modified as in the following equation.

【0090】[0090]

【数14】 [Equation 14]

【0091】ここに、λは本実施の形態の特徴となる誤
差の分配を行うための重み係数ベクトルである。このλ
は以下のような役割を持つ。
Here, λ is a weight coefficient vector for distributing an error which is a feature of the present embodiment. This λ
Has the following roles:

【0092】すなわち、従来技術では、制御周期T2
十分に短くできることを仮定しているため、パルスジェ
ネレータ37のパルス発生時点と、パルス読み込み時点
とのずれは、高々T2 で無視できるものとしている。こ
の前提によれば位置検出手段38が検出した位置を正確
なものとみなすことができ、よって、差分Δθの生じる
原因を推定外乱の初期値Td [j,0](∧付)及び速
度推定値の初期値ω[j,0](∧付)の誤差によるも
のとすることができる。
That is, in the prior art, since it is assumed that the control period T 2 can be made sufficiently short, the difference between the pulse generation time of the pulse generator 37 and the pulse reading time can be ignored at most at T 2. I have. According to this premise, the position detected by the position detecting means 38 can be regarded as accurate. Therefore, the cause of the difference Δθ is determined by the initial value T d [j, 0] of the estimated disturbance and the speed estimation. It can be attributed to the error of the initial value ω [j, 0] (appended).

【0093】しかし、既に述べたように、制御周期T2
が十分に短くできない場合は、位置検出手段のパルス発
生時点と位置読み込み時点との間のずれは大きくなり、
従って、位置読み込み時点における真の位置と位置検出
手段からの位置の間の量子化誤差が増大する。このよう
な要因が無視できない場合は、Δθの生じる原因を従来
のように推定外乱の初期値と速度推定値の初期値に帰す
る前提の下では、量子化ノイズの増大、パラメータ変動
に対するロバスト性の劣化、オブザーバの不安定化とい
った不具合を生じる。
However, as described above, the control cycle T 2
Is not sufficiently short, the difference between the pulse generation time of the position detection means and the position reading time becomes large,
Therefore, the quantization error between the true position at the time of reading the position and the position from the position detecting means increases. If such factors cannot be ignored, under the premise that the cause of Δθ is attributable to the initial value of the estimated disturbance and the initial value of the speed estimation value as in the conventional case, the robustness against the increase in quantization noise and parameter fluctuation And the observer becomes unstable.

【0094】そこで、本実施の形態では、位置検出手段
38からの位置には、量子化誤差などの位置読み取り誤
差が含まれているとみなし、差分Δθの生じる原因に
は、推定外乱の初期値、速度推定値の初期値と共に、位
置読み取り誤差も含まれていると考える。この場合、Δ
θを各誤差に分けて書くと次式のようになる。
Therefore, in the present embodiment, it is considered that the position from the position detecting means 38 includes a position reading error such as a quantization error, and the cause of the difference Δθ is the initial value of the estimated disturbance. It is considered that the position reading error is included together with the initial value of the speed estimation value. In this case, Δ
When θ is written for each error, the following equation is obtained.

【0095】[0095]

【数15】 (Equation 15)

【0096】である。但し、Δθq が位置読み取り誤差
を表しており、推定誤差Δθのλ0 の分がΔθq の影響
であることを意味している。(20)式に基づいてxerr
計算し、これを次の推定期間[j+1]における初期値
の修正に用いるには、以下のようにする。
Is as follows. However, Δθ q represents the position reading error, which means that the estimated error Δθ of λ 0 is the effect of Δθ q . In order to calculate x err based on equation (20) and use it for correcting the initial value in the next estimation period [j + 1], the following is performed.

【0097】[0097]

【数16】 (Equation 16)

【0098】である。(21)式を整理すると、Is as follows. Rearranging equation (21),

【0099】[0099]

【数17】 [Equation 17]

【0100】となるので、これをステップ114の誤差
修正演算に用いるのである。一方、誤差修正時点におけ
る推定誤差の収束特性は、
This is used for the error correction calculation in step 114. On the other hand, the convergence characteristic of the estimation error at the time of error correction is

【0101】[0101]

【数18】 (Equation 18)

【0102】となる。従って、上記推定誤差が時間と共
に零に収束するには、方程式 |zI−(A−GCA)| = 0 (24) の根が、z平面における原点を中心とした半径1の円内
に入るようにゲイン行列Gを設計すればよい。
Is obtained. Therefore, in order for the above estimation error to converge to zero with time, the root of the equation | zI- (A-GCA) | May be designed as the gain matrix G.

【0103】Gとλとの間には、Between G and λ,

【0104】[0104]

【数19】 [Equation 19]

【0105】で示される関係が導き出されるため、(25)
式を用いて(24)式を整理すると、
Since the relationship shown in the following is derived, (25)
By rearranging equation (24) using the equation,

【0106】[0106]

【数20】 (Equation 20)

【0107】を得る。よって、オブザーバのz平面上の
極を(z1 ,z2 ,z3 )に設定するには、次の連立方
程式を解けばよいことになる。
Is obtained. Therefore, in order to set the pole on the z-plane of the observer to (z 1 , z 2 , z 3 ), the following simultaneous equations may be solved.

【0108】[0108]

【数21】 (Equation 21)

【0109】従って、λは、希望の特性(z1 ,z2
3 )を設定し、(27)式を解くことによって得ることが
できる。
Therefore, λ is the desired characteristic (z 1 , z 2 ,
z 3 ) and solving equation (27).

【0110】そして、ステップ115では、誤差修正タ
イミングのインデックスjを1つカウントアップすると
同時に、ステップ116、ステップ117において、k
及びnを零にリセットし、ステップ100に戻り再び繰
り返しループの演算を実行する。
In step 115, the index j of the error correction timing is counted up by one, and at the same time, in steps 116 and 117, k
, And n are reset to zero, and the flow returns to step 100 to execute the operation of the iterative loop again.

【0111】このように本実施の形態では、位置検出手
段38からの位置には、量子化誤差などの位置読み取り
誤差が含まれているとみなし、推定値誤差の生じる原因
を推定外乱の初期値と速度推定値の初期値の誤差、及び
位置読み取り誤差によるものと考え、推定値の修正を外
乱、速度に加えて、位置に対しても実施している。
As described above, in the present embodiment, it is assumed that the position from the position detecting means 38 includes a position reading error such as a quantization error, and the cause of the estimated value error is determined by the initial value of the estimated disturbance. It is considered that the error is caused by the error of the initial value of the speed estimation value and the position reading error, and the estimation value is corrected for the position in addition to the disturbance and the speed.

【0112】従って、演算デバイスの能力によって制御
周期T2 が十分に短くできない場合に生じる位置読み取
り誤差の増大に対しても、安定かつ正確に外乱トルク及
び位置を推定することができる。また、推定誤差修正タ
イミングを位置検出手段からのパルス発生回数に応じて
変化させているので、パルス発生毎に推定誤差の修正を
行う従来技術(第2の態様)にみられる高速度域での演
算負荷の増大を回避することができる。
Therefore, the disturbance torque and the position can be stably and accurately estimated even if the position reading error increases when the control period T 2 cannot be sufficiently shortened due to the capability of the arithmetic device. Further, since the estimation error correction timing is changed in accordance with the number of pulse generations from the position detection means, the estimation error correction is performed for each pulse generation in the high-speed range seen in the prior art (second mode). An increase in calculation load can be avoided.

【0113】さらに、パルス計数値から求められる位置
(回転角)と真の位置との差を推定誤差要因の一つとし
て考慮し、位置に対しても誤差修正を行うため、制御周
期内にパルスが発生しないような低速時においても良好
に速度を推定できる。 (タイヤ発生力検出装置)本発明の実施の形態に係るタ
イヤ発生力検出装置の構成ブロック図を図2に示す。同
図に示すように、本タイヤ発生力検出装置は、図1の回
転速度検出装置と同様の構成要件(同一符号)と共に、
さらに、推定誤差修正手段36により修正された外乱ト
ルクの推定値を、タイヤと路面との間に発生するタイヤ
発生力(制動力)として検出するタイヤ発生力推定手段
40と、を備えている。
Further, the difference between the position (rotation angle) obtained from the pulse count value and the true position is considered as one of the estimation error factors, and the position is corrected for error. The speed can be satisfactorily estimated even at a low speed at which no occurrence occurs. (Tire Generating Force Detecting Apparatus) FIG. 2 is a block diagram showing the configuration of the tire generating force detecting apparatus according to the embodiment of the present invention. As shown in the figure, the tire generating force detecting device has the same components (same reference numerals) as those of the rotational speed detecting device of FIG.
Further, there is provided a tire generating force estimating means 40 for detecting the estimated value of the disturbance torque corrected by the estimation error correcting means 36 as a tire generating force (braking force) generated between the tire and the road surface.

【0114】次に、図2のタイヤ発生力検出装置の作用
を図4のフローチャートに基づいて説明する。なお、図
4のフローチャートにおいて、図3のフローチャートと
同様の作用部分(ステップ)については同一の符号を付
して説明を省略し、異なる作用部分についてのみ説明す
る。
Next, the operation of the tire generating force detection device of FIG. 2 will be described with reference to the flowchart of FIG. In the flowchart of FIG. 4, the same operation parts (steps) as those in the flowchart of FIG. 3 are denoted by the same reference numerals, and description thereof will be omitted. Only different operation parts will be described.

【0115】図4のフローチャートに示すように、ステ
ップ104において、車輪に作用する外乱トルクT
d [j,k]の推定値に基づいて以下のように制動力を
推定する。
As shown in the flowchart of FIG. 4, in step 104, the disturbance torque T acting on the wheels
d Based on the estimated value of [j, k], the braking force is estimated as follows.

【0116】車輪の制動時の運動方程式は、The equation of motion at the time of braking the wheel is:

【0117】[0117]

【数22】 (Equation 22)

【0118】で表せる。ここに、Wは輪荷重、Pはホイ
ールシリンダ圧、Kはホイールシリンダ圧から制動トル
クTb への変換係数である。今、ステップ101の制動
トルク検出において、 Tq = −KP (29) と制動トルクが検出されるので、
Can be expressed by Here, W is the wheel load, P is a conversion coefficient of the wheel cylinder pressure, K is the wheel cylinder pressure to the braking torque T b. Now, the braking torque detecting step 101, the braking torque is detected to T q = -KP (29),

【0119】[0119]

【数23】 (Equation 23)

【0120】という関係が成り立ち、外乱トルクの推定
値が制動力μWにタイヤの動荷重半径Rを乗じたものに
対応していることがわかる。そこで、ステップ104に
おいて、外乱トルクの推定値からタイヤ動荷重半径Rを
割れば制動力μWを推定することができる。なお、外乱
トルクの初期値が、例えば(15)式等に基づいて与えられ
たとしても、外乱トルクの推定値は、誤差修正演算(ス
テップ109〜114)により順次修正されるので正確
な推定値となり、よって制動力の正確な推定が可能とな
る。
It can be seen that the estimated value of the disturbance torque corresponds to a value obtained by multiplying the braking force μW by the dynamic load radius R of the tire. Then, in step 104, the braking force μW can be estimated by dividing the tire dynamic load radius R from the estimated value of the disturbance torque. Even if the initial value of the disturbance torque is given based on, for example, equation (15) or the like, the estimated value of the disturbance torque is sequentially corrected by the error correction calculation (steps 109 to 114). Therefore, it is possible to accurately estimate the braking force.

【0121】さらに、輪荷重Wをセンサや演算等により
検出すれば、タイヤと路面間の摩擦係数μも推定するこ
とができる。推定した制動力、摩擦係数等は、アンチス
キッド装置や車両安定化制御装置などに用いることがで
きる。
Further, if the wheel load W is detected by a sensor or calculation, the friction coefficient μ between the tire and the road surface can be estimated. The estimated braking force, friction coefficient, and the like can be used for an anti-skid device, a vehicle stabilization control device, and the like.

【0122】[0122]

【実施例】以下、図面を参照して、本発明の実施例に係
る回転速度検出装置及びタイヤ発生力検出装置の動作実
験結果を説明する。 (回転速度検出装置)本実施例に係る回転速度検出装置
は、電動機、車両の車輪といった回転体の回転速度検出
装置(図1)に関するものであるが、本実施例では、車
両の車輪を回転体として扱う。なお、電動機の電動子及
び車両の車輪の間には、表1に示したような対応関係が
あり、本実施例に係る車輪の回転速度検出装置から、電
動機の回転速度検出装置への適用は、図10の変換表に
基づいて諸量を変換することにより簡単に得られる。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Referring to the drawings, description will be given below of experimental results of operation of a rotation speed detecting device and a tire generating force detecting device according to an embodiment of the present invention. (Rotation Speed Detection Device) The rotation speed detection device according to the present embodiment relates to a rotation speed detection device for a rotating body such as an electric motor or a vehicle wheel (FIG. 1). Treat as a body. Note that there is a correspondence relationship between the motor of the motor and the wheels of the vehicle as shown in Table 1, and the application from the wheel rotation speed detecting device according to the present embodiment to the motor rotation speed detecting device is not applicable. Can be easily obtained by converting various quantities based on the conversion table of FIG.

【0123】また、本実施例に係る回転速度検出装置
は、図3のフローチャートに基づいて回転速度を検出す
るが、本実施例の動作実験では、図3のフローチャート
において、M=1、N=1としている。車両制御におい
ては、一般に制御周期T2 は数msと長いため、Mを高
く設定すると、推定誤差周期T1 が長くなり過ぎるた
め、Mは1〜2程度が望ましいからである。また、現状
の車両に装着されている車輪速センサの分解能は、一回
転当たり48パルスと極めて低いため、誤差修正周期T
1 の期間でパルスが1回でも発生すれば誤差修正を行っ
た方が良好な推定が実現できるので、Nも1〜2回程度
が望ましいからである。
Further, the rotation speed detecting device according to the present embodiment detects the rotation speed based on the flowchart of FIG. 3, but in the operation experiment of the present embodiment, in the flowchart of FIG. It is set to 1. In a vehicle control, generally the control period T 2 for a long and a few ms, the set high M, since the estimation error period T 1 is too long, because M is desirably about 1-2. Further, the resolution of the wheel speed sensor currently mounted on the vehicle is extremely low at 48 pulses per rotation, so that the error correction period T
This is because, if a pulse is generated even once in one period, better estimation can be realized by performing error correction. Therefore, it is desirable that N is about 1 to 2 times.

【0124】本実施例の動作実験では、50km/hで
定速走行している状態から、ステップ状に制動トルクを
印加し、この条件下での車両の前輪の瞬時速度を推定し
た。なお、制御周期T2 は5msである。以下、この動
作実験の結果について図11〜図14を用いて説明す
る。
In the operation experiment of this embodiment, the braking torque was applied in a stepwise manner while the vehicle was running at a constant speed of 50 km / h, and the instantaneous speed of the front wheels of the vehicle under these conditions was estimated. The control period T 2 are a 5 ms. Hereinafter, results of the operation experiment will be described with reference to FIGS.

【0125】図11には、上記条件下での車体速度、車
輪速度及び制御周期内で発生するパルス数の時間的変化
が示されている。制動トルクを印加しているため、時間
と共に、車体速度及び車輪速度が減少し、制御周期内の
パルス数も減少していくことがわかる。特に、車輪速度
(実線)が25(km/h)付近まで低下すると、パル
ス発生期間が制御周期より長くなり、制御周期の間でパ
ルスが入らなくなる期間も現れる。
FIG. 11 shows changes over time in the vehicle speed, the wheel speed, and the number of pulses generated in the control cycle under the above conditions. Since the braking torque is applied, the vehicle speed and the wheel speed decrease with time, and the number of pulses in the control cycle also decreases. In particular, when the wheel speed (solid line) decreases to around 25 (km / h), the pulse generation period becomes longer than the control cycle, and a period during which the pulse does not enter during the control cycle also appears.

【0126】図12(a)には、本実施例に係る回転速
度検出装置の速度推定結果が示されている。同図では、
破線が真の瞬時速度、実線が速度推定値であるが、動作
実験の速度範囲内できわめて良好に速度を推定している
ため、破線及び実線が一本に重なっている。
FIG. 12A shows the speed estimation result of the rotation speed detecting device according to the present embodiment. In the figure,
The dashed line is the true instantaneous speed, and the solid line is the estimated speed value. Since the speed is estimated very well within the speed range of the operation experiment, the dashed line and the solid line overlap.

【0127】一方、図12(b)には、従来技術の瞬時
速度オブザーバによる同条件下での速度推定結果が示さ
れている。図12(a)の本実施例の結果では、車輪速
度が低速度域まで良好に推定されているのに対し、図1
2(b)の従来技術の推定結果では、低速度域において
推定値が振動し、不安定な結果となっていることがわか
る。これは、低速になるほどパルス計数値から求められ
る位置(回転角)と真の位置との差が大きくなるが、従
来技術では、パルス計数値から求められる位置を真の位
置として扱い、誤差修正しているためである。一方、本
発明の実施例では、パルス計数値から求められる位置
(回転角)と真の位置との差を推定誤差要因の一つとし
て考慮し、位置に対しても誤差修正を行うため、低速時
にも良好に速度を推定できるのである。
On the other hand, FIG. 12B shows the result of speed estimation under the same conditions by the instantaneous speed observer of the prior art. According to the results of the present embodiment in FIG. 12A, the wheel speed is well estimated up to the low speed range, whereas FIG.
It can be seen from the estimation result of the prior art of FIG. 2 (b) that the estimation value fluctuates in the low speed range, resulting in an unstable result. This is because the difference between the position (rotation angle) determined from the pulse count value and the true position increases as the speed decreases, but in the related art, the position determined from the pulse count value is treated as the true position and error correction is performed. Because it is. On the other hand, in the embodiment of the present invention, the difference between the position (rotation angle) obtained from the pulse count value and the true position is considered as one of the estimation error factors, and the position is corrected for error. Sometimes the speed can be estimated well.

【0128】また、図13、図14には、それぞれ、本
発明の実施例及び従来技術における瞬時速度推定値の推
定誤差と推定速度との関係が示されている。これらの図
において、横軸は連続時間系に変換したオブザーバの極
(rad/s)であり、左に行くほど推定速度の設定が
速くなる。また、縦軸は推定誤差を表しており、上に行
くほど誤差が大きいことを意味している。また、本発明
においては規定値Mを、従来技術においては、(1) 式に
おける整数mをパラメータとして変化させており、M、
mの値が小さいほど誤差修正周期T1 が短くなる。
FIGS. 13 and 14 show the relationship between the estimation error of the instantaneous speed estimation value and the estimated speed in the embodiment of the present invention and the prior art, respectively. In these figures, the horizontal axis is the pole (rad / s) of the observer converted into a continuous time system, and the setting of the estimated speed becomes faster toward the left. The vertical axis represents the estimation error, which means that the error increases as going upward. Further, in the present invention, the specified value M is changed using the integer m in the equation (1) as a parameter in the prior art.
more error correction cycle T 1 value of m is small is shortened.

【0129】従来技術の結果では、図14に示すよう
に、m=1、極=−8(rad/s)で誤差が最小値を迎える
が、本実施例の結果では、図13に示すように、M=
1、極=−23(rad/s)で誤差が最小となっている。し
かも、誤差の最小値は、本実施例の方が従来技術よりも
小さな値となっている。従って、本発明によれば、従来
技術と比べて、推定速度を速く設定でき、かつ推定誤差
も小さく抑えられることがわかる。 (タイヤ発生力検出装置)本実施例に係るタイヤ発生力
検出装置の動作実験では、本実施例の回転速度検出装置
と同様に、車両が50km/hで定速走行している状態
から、ステップ状に制動トルクを印加し、この条件下で
車輪に発生した制動力を推定した。なお、本動作実験で
は、制御周期T2 を5ms、図4のフローチャートにお
いてM=N=1とした。以下、この動作実験の結果につ
いて図15を用いて説明する。
In the result of the prior art, as shown in FIG. 14, the error reaches the minimum value at m = 1 and pole = −8 (rad / s), but in the result of the present embodiment, as shown in FIG. And M =
1, the error is minimum at pole = −23 (rad / s). Moreover, the minimum value of the error is smaller in the present embodiment than in the prior art. Therefore, according to the present invention, it can be seen that the estimation speed can be set faster and the estimation error can be reduced as compared with the related art. (Tire generating force detecting device) In the operation experiment of the tire generating force detecting device according to the present embodiment, as in the case of the rotational speed detecting device of the present embodiment, the vehicle is traveling at a constant speed of 50 km / h and the step is started. The braking torque was applied in the same manner, and the braking force generated on the wheels under these conditions was estimated. In the present operation experiment, the control period T 2 to 5 ms, and M = N = 1 in the flowchart of FIG. Hereinafter, the result of the operation experiment will be described with reference to FIG.

【0130】図15には、上記条件下において本実施例
のタイヤ発生力検出装置が推定した制動力の推定値の時
間的変化が示されている。なお、同図において、破線が
実測された制動力の真値、実線が制動力の推定値を示し
ている。また、参考までに、制御周期間に発生するパル
ス数も示してある。
FIG. 15 shows the change over time of the estimated value of the braking force estimated by the tire generating force detection device of this embodiment under the above conditions. In the figure, the broken line indicates the actual value of the actually measured braking force, and the solid line indicates the estimated value of the braking force. Also, for reference, the number of pulses generated during the control cycle is also shown.

【0131】図15に示すように、制動トルクを印加し
ているため、時間と共に車輪速が低下し、やがてパルス
発生期間が制御周期より長くなり、制御周期の間でパル
スが入らなくなる期間も現れる。このような低速域にお
いても図15では、真の制動力の値と推定された制動力
の値との一致度がきわめて良く、制動力が良好に推定さ
れていることがわかる。
As shown in FIG. 15, since the braking torque is applied, the wheel speed decreases with time, the pulse generation period becomes longer than the control cycle, and a period during which no pulse is input appears during the control cycle. . Even in such a low-speed range, in FIG. 15, the degree of coincidence between the true braking force value and the estimated braking force value is extremely good, and it can be seen that the braking force is well estimated.

【0132】[0132]

【発明の効果】以上説明したように請求項1の発明によ
れば、位置の誤差に基づいて、外乱トルクの推定値、回
転速度の初期値及び回転位置の初期値を修正するように
したので、従来技術に比べて推定速度を速く設定でき、
かつ推定誤差も小さく抑えることができる。特に、演算
デバイスの能力によって制御系の制御周期が十分に短く
できず、よって、回転位置の量子化誤差が増大する場合
においても、正確に回転速度を推定することができると
共に、量子化ノイズの増大、パラメータ変動に対するロ
バスト性の劣化、及びオブザーバの不安定化などを防止
することができる、という効果が得られる。
As described above, according to the first aspect of the present invention, the estimated value of the disturbance torque, the initial value of the rotational speed, and the initial value of the rotational position are corrected based on the position error. , The estimation speed can be set faster than the prior art,
In addition, the estimation error can be kept small. In particular, even when the control cycle of the control system cannot be sufficiently shortened due to the capability of the arithmetic device, and thus the quantization error of the rotation position increases, the rotation speed can be accurately estimated and the quantization noise can be reduced. The effect is obtained that it is possible to prevent the increase, the deterioration of the robustness to the parameter fluctuation, and the instability of the observer.

【0133】さらに、誤差修正の時点をパルス数に基づ
いて決定することにより、制御周期でパルスが発生しな
いような低速度域においても、正確に回転速度を推定す
ることができる。
Further, by determining the time of error correction based on the number of pulses, it is possible to accurately estimate the rotational speed even in a low speed range where no pulse is generated in the control cycle.

【0134】請求項2の発明では、請求項1の回転速度
検出装置と同じ原理で、制御周期にパルスが発生しない
ような低速度域及び高速度域のいずれにおいても、正確
かつ安定に外乱トルクを推定できるので、常に良好に制
動力を推定することができる。また、従来技術にあるよ
うな車輪速度の差分を用いることがないので、種々のノ
イズに強く、安定して制動力を推定することができる、
という効果が得られる。
According to the second aspect of the present invention, the disturbance torque can be accurately and stably obtained in both the low speed range and the high speed range in which no pulse is generated in the control cycle, based on the same principle as the rotation speed detecting device of the first aspect. Therefore, the braking force can always be satisfactorily estimated. Further, since the difference in wheel speed as in the related art is not used, it is resistant to various noises, and the braking force can be stably estimated.
The effect is obtained.

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

【図1】本発明の実施の形態に係る回転速度検出装置の
構成を示すブロック図である。
FIG. 1 is a block diagram illustrating a configuration of a rotation speed detection device according to an embodiment of the present invention.

【図2】本発明の実施の形態に係るタイヤ発生力検出装
置の構成を示すブロック図である。
FIG. 2 is a block diagram showing a configuration of a tire generation force detection device according to an embodiment of the present invention.

【図3】本実施の形態に係る回転速度検出装置の処理の
流れを示すフローチャートである。
FIG. 3 is a flowchart showing a processing flow of the rotation speed detecting device according to the present embodiment.

【図4】本実施の形態に係るタイヤ発生力検出装置の処
理の流れを示すフローチャートである。
FIG. 4 is a flowchart showing a processing flow of the tire generating force detection device according to the present embodiment.

【図5】従来の瞬時速度オブザーバの構成を示すブロッ
ク図である。
FIG. 5 is a block diagram showing a configuration of a conventional instantaneous speed observer.

【図6】従来の瞬時速度オブザーバの第1態様「位置読
み込み形」における位置読み込みタイミングを示す図で
ある。
FIG. 6 is a diagram showing a position reading timing in the first mode “position reading type” of the conventional instantaneous speed observer.

【図7】従来の瞬時速度オブザーバの第2態様「平均速
度読み込み形」における平均速度読み込みタイミングを
示す図である。
FIG. 7 is a diagram showing an average speed reading timing in a second mode “average speed reading type” of the conventional instantaneous speed observer.

【図8】従来の車輪速度演算法を説明するための図であ
る。
FIG. 8 is a diagram for explaining a conventional wheel speed calculation method.

【図9】従来の瞬時速度オブザーバの第2態様「平均速
度読み込み形」における平均速度読み込みタイミングと
誤差修正タイミングのずれを説明するための図である。
FIG. 9 is a diagram for explaining a difference between an average speed reading timing and an error correction timing in the second mode “average speed reading type” of the conventional instantaneous speed observer.

【図10】電動機のパラメータと車輪のパラメータとの
対応関係を示す変換表である。
FIG. 10 is a conversion table showing the correspondence between motor parameters and wheel parameters.

【図11】本発明の実施例に係る回転速度検出装置の動
作実験で用いられる車両の車体/車輪速度及び制御周期
内のパルス数の時間的変化を示す図である。
FIG. 11 is a diagram showing a temporal change of the vehicle body / wheel speed and the number of pulses in a control cycle used in an operation experiment of the rotation speed detecting device according to the embodiment of the present invention.

【図12】本実施例での動作実験において推定された回
転速度の時間的変化を示す図であって、(a)は本実施
例に係る回転速度検出装置の推定結果、(b)は従来の
瞬時速度オブザーバでの推定結果を示す図である。
12A and 12B are diagrams showing a temporal change of a rotation speed estimated in an operation experiment in the present embodiment, wherein FIG. 12A is a result of estimation by the rotation speed detection device according to the embodiment, and FIG. FIG. 7 is a diagram showing an estimation result of the instantaneous speed observer of FIG.

【図13】本実施例に係る回転速度検出装置での瞬時速
度推定値の推定誤差と極(rad/s)との関係を示す
図である。
FIG. 13 is a diagram showing a relationship between an estimation error of an instantaneous speed estimation value and a pole (rad / s) in the rotation speed detection device according to the embodiment.

【図14】従来の瞬時速度オブザーバでの瞬時速度推定
値の推定誤差と極(rad/s)との関係を示す図であ
る。
FIG. 14 is a diagram showing a relationship between an estimation error of an estimated instantaneous speed value in a conventional instantaneous speed observer and a pole (rad / s).

【図15】本実施例に係るタイヤ発生力検出装置の制動
力の推定結果と車輪速パルス数とを示す図である。
FIG. 15 is a diagram illustrating a result of estimating a braking force of the tire generating force detection device according to the present embodiment and a wheel speed pulse number.

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

30 制御トルク検出手段 32 状態推定手段 34 誤差修正判断手段 36 推定誤差修正手段 37 パルスジェネレータ 38 位置検出手段 40 タイヤ発生力推定手段 Reference Signs List 30 control torque detecting means 32 state estimating means 34 error correction determining means 36 estimation error correcting means 37 pulse generator 38 position detecting means 40 tire generating force estimating means

───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 昭60−171464(JP,A) 特開 平8−233840(JP,A) 特開 平7−89304(JP,A) (58)調査した分野(Int.Cl.7,DB名) G01L 5/00 G01P 3/489 H02P 5/00 ────────────────────────────────────────────────── ─── Continuation of the front page (56) References JP-A-60-17164 (JP, A) JP-A-8-233840 (JP, A) JP-A-7-89304 (JP, A) (58) Field (Int.Cl. 7 , DB name) G01L 5/00 G01P 3/489 H02P 5/00

Claims (2)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 回転体に取り付けられたパルス生成器が
出力したパルスの計数値に基づいて回転体の回転位置を
検出する位置検出手段と、 回転体を制御するための制御トルクを検出する制御トル
ク検出手段と、 前記制御トルク検出手段により検出された値の制御トル
ク及び推定された値の外乱トルクが作用する回転体のモ
デルに基づいて、与えられた回転速度及び回転位置の各
初期値とから回転体の回転速度及び回転位置を推定する
状態推定手段と、 少なくとも前記パルス生成器が出力したパルスの数に基
づいて、誤差修正の時点を判断する誤差修正判断手段
と、 前記誤差修正判断手段により判断された誤差修正の時点
において、前記位置検出手段により検出された回転位置
と前記状態推定手段により推定された回転位置との間に
生じた位置の誤差を、外乱トルクの推定値の誤差、回転
速度の初期値の誤差、及び検出された回転位置に含まれ
る量子化誤差によるものとみなし、前記位置の誤差に基
づいて、前記状態推定手段が用いる外乱トルクの推定値
と回転速度及び回転位置の各初期値とを修正する推定誤
差修正手段と、 を含む回転速度検出装置。
1. A position detecting means for detecting a rotational position of a rotating body based on a count value of pulses output from a pulse generator attached to the rotating body, and a control for detecting a control torque for controlling the rotating body. Based on a model of the rotating body on which the control torque of the value detected by the control torque detection means and the disturbance torque of the estimated value act, based on the initial value of the given rotation speed and rotation position, State estimating means for estimating the rotational speed and rotational position of the rotating body from an error correction determining means for determining the time of error correction based on at least the number of pulses output by the pulse generator; and the error correction determining means At the time of the error correction determined by the above, the position generated between the rotational position detected by the position detecting means and the rotational position estimated by the state estimating means. The error of the estimated value of the disturbance torque, the error of the initial value of the rotational speed, and the quantization error included in the detected rotational position, based on the position error, the state estimating means An estimation error correction unit that corrects the estimated value of the disturbance torque to be used and each initial value of the rotation speed and the rotation position.
【請求項2】 車輪に取り付けられたパルス生成器が出
力したパルスの計数値に基づいて車輪の回転位置を検出
する位置検出手段と、 車輪の制動トルクを検出する制動トルク検出手段と、 前記制動トルク検出手段により検出された値の制動トル
ク及び推定された値の外乱トルクが作用する車輪のモデ
ルに基づいて、与えられた回転速度及び回転位置の各初
期値とから車輪の回転速度及び回転位置を推定する状態
推定手段と、 少なくとも前記パルス生成器が出力したパルスの数に基
づいて、誤差修正の時点を判断する誤差修正判断手段
と、 前記誤差修正判断手段により判断された誤差修正の時点
において、前記位置検出手段により検出された回転位置
と前記状態推定手段により推定された回転位置との間に
生じた位置の誤差を、外乱トルクの推定値の誤差、回転
速度の初期値の誤差、及び検出された回転位置に含まれ
る量子化誤差によるものとみなし、前記位置の誤差に基
づいて、前記状態推定手段が用いる外乱トルクの推定値
と回転速度及び回転位置の各初期値とを修正する推定誤
差修正手段と、 前記推定誤差修正手段により修正された外乱トルクと車
輪の動半径とに基づいて、路面から反力として車輪に作
用する制動力を推定するタイヤ発生力推定手段と、 を含むタイヤ発生力検出装置。
2. A position detecting means for detecting a rotational position of a wheel based on a count value of a pulse output from a pulse generator attached to the wheel; a braking torque detecting means for detecting a braking torque of the wheel; Based on the model of the wheel on which the braking torque of the value detected by the torque detecting means and the disturbance torque of the estimated value act, the rotational speed and the rotational position of the wheel are obtained from the given initial values of the rotational speed and the rotational position. State estimation means for estimating, at least based on the number of pulses output by the pulse generator, error correction determination means for determining the time of error correction, and at the time of error correction determined by the error correction determination means The error of the position generated between the rotational position detected by the position detecting means and the rotational position estimated by the state estimating means is calculated as the disturbance torque. Based on the error of the constant value, the error of the initial value of the rotational speed, and the quantization error included in the detected rotational position, the estimated value of the disturbance torque and the rotational Estimated error correcting means for correcting the initial values of the speed and the rotational position; and a braking force acting on the wheels as a reaction force from the road surface based on the disturbance torque and the moving radius of the wheels corrected by the estimated error correcting means. A tire generation force estimating means for estimating the tire generation force.
JP04610097A 1997-02-28 1997-02-28 Rotation speed detection device and tire generation force detection device Expired - Fee Related JP3228168B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP04610097A JP3228168B2 (en) 1997-02-28 1997-02-28 Rotation speed detection device and tire generation force detection device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP04610097A JP3228168B2 (en) 1997-02-28 1997-02-28 Rotation speed detection device and tire generation force detection device

Publications (2)

Publication Number Publication Date
JPH10239187A JPH10239187A (en) 1998-09-11
JP3228168B2 true JP3228168B2 (en) 2001-11-12

Family

ID=12737585

Family Applications (1)

Application Number Title Priority Date Filing Date
JP04610097A Expired - Fee Related JP3228168B2 (en) 1997-02-28 1997-02-28 Rotation speed detection device and tire generation force detection device

Country Status (1)

Country Link
JP (1) JP3228168B2 (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070264564A1 (en) 2006-03-16 2007-11-15 Infinite Power Solutions, Inc. Thin film battery on an integrated circuit or circuit board and method thereof
US7993773B2 (en) 2002-08-09 2011-08-09 Infinite Power Solutions, Inc. Electrochemical apparatus with barrier layer protected substrate
US7102271B2 (en) * 2004-01-14 2006-09-05 Infinite Power Solutions, Inc. Method and apparatus for a high output sensor system
TWI268347B (en) 2005-07-25 2006-12-11 Sin Etke Technology Co Ltd Devices of speed of a motor vehicle with concurrent selector accumulative value being substituted into the convertible function to compute a practical value
JP4677953B2 (en) * 2006-05-26 2011-04-27 トヨタ自動車株式会社 Traveling apparatus and control method thereof
KR20100102180A (en) 2007-12-21 2010-09-20 인피니트 파워 솔루션스, 인크. Method for sputter targets for electrolyte films
KR101606183B1 (en) 2008-01-11 2016-03-25 사푸라스트 리써치 엘엘씨 Thin film encapsulation for thin film batteries and other devices
US8906523B2 (en) 2008-08-11 2014-12-09 Infinite Power Solutions, Inc. Energy device with integral collector surface for electromagnetic energy harvesting and method thereof
CN102576828B (en) 2009-09-01 2016-04-20 萨普拉斯特研究有限责任公司 There is the printed circuit board (PCB) of integrated thin film battery
JP2013528912A (en) 2010-06-07 2013-07-11 インフィニット パワー ソリューションズ, インコーポレイテッド Rechargeable high density electrochemical device
CN114465533A (en) * 2022-01-26 2022-05-10 山东新一代信息产业技术研究院有限公司 Method and tool for measuring speed of ultralow rotating speed of brushless direct current motor

Also Published As

Publication number Publication date
JPH10239187A (en) 1998-09-11

Similar Documents

Publication Publication Date Title
US10144445B2 (en) Modified static tire model for providing assist without a torque sensor for zero to low vehicle speeds
JP3228168B2 (en) Rotation speed detection device and tire generation force detection device
JP2000270600A (en) Drive of induction motor, and estimating method of its parameters
EP0514847B1 (en) System and method for controlling speed of electric motor in extremely low speed range using rotary pulse encoder
JPS61164481A (en) Speed control system
US20070043528A1 (en) Speed measurement system for speed control of high-speed motors
JPS6294249A (en) Speed controlling system
US4977525A (en) Highly efficient method and apparatus for calculating rotational speed
JP3277454B2 (en) Detection position correction device and rotation speed detection device
JP3617309B2 (en) Road friction coefficient estimation device
US6543278B1 (en) Apparatus for estimating road friction coefficient
JP2009095154A (en) Motor controller and its speed detection method
JP3412492B2 (en) Wheel speed detection device
JP3297182B2 (en) Electric vehicle mileage calculation device
JP2001033472A (en) Rotary speed-detecting device
JP3826204B2 (en) Electric vehicle tangential force coefficient estimation device
JP3441562B2 (en) Yaw rate estimator and yaw rate sensor failure diagnosis device
JP3271964B1 (en) Road surface friction coefficient estimation device for vehicles
JP2000055790A (en) Apparatus for estimating friction coefficient ?? of road- surface
JP3402017B2 (en) Speed detection control device
JP3339208B2 (en) Motor control device
JP3198690B2 (en) Speed detector for rotating body
JP3329171B2 (en) Rotation speed detector
JPH1038906A (en) Detector for pulse width and detector for rotational speed
KR100244245B1 (en) Method for detecting velocity of an electronic motor

Legal Events

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