JPH11192959A - Kinetic condition quantity estimating device for vehicle - Google Patents

Kinetic condition quantity estimating device for vehicle

Info

Publication number
JPH11192959A
JPH11192959A JP10000164A JP16498A JPH11192959A JP H11192959 A JPH11192959 A JP H11192959A JP 10000164 A JP10000164 A JP 10000164A JP 16498 A JP16498 A JP 16498A JP H11192959 A JPH11192959 A JP H11192959A
Authority
JP
Japan
Prior art keywords
motion state
vehicle
yaw rate
vehicle motion
value
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.)
Pending
Application number
JP10000164A
Other languages
Japanese (ja)
Inventor
Takashi Nakajima
孝 中島
Kazutaka Adachi
和孝 安達
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.)
Nissan Motor Co Ltd
Original Assignee
Nissan Motor Co Ltd
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 Nissan Motor Co Ltd filed Critical Nissan Motor Co Ltd
Priority to JP10000164A priority Critical patent/JPH11192959A/en
Publication of JPH11192959A publication Critical patent/JPH11192959A/en
Pending legal-status Critical Current

Links

Abstract

PROBLEM TO BE SOLVED: To provide a kinetic condition quantity estimation device wherein an estimated value of a kinetic condition quantity can be coincided with an actual true value, in no relation to a delay of change of the kinetic condition quantity relating to steering. SOLUTION: To an observer with a yaw rate γ as a reference input, a detection value of car speed Vx , steering angle θ, yaw rate γ, and a road surface μ, is given, in the case of obtaining a yaw rate estimation value γH and side slip angle estimation value βH, a delay deviation γe between the yaw rate estimation value γB via a matrix C and a yaw rate detection value γ is returned by a gain K to a preceding step of an integrator 1/S, this γe , θ, Vx and a vehicle model A are used, the estimation value γH and βH are obtained. In the case of calculating the delay deviation γe , delaying by a waste time (lengthened in accordance with decreasing a road surface friction coefficient μ) exp(-ls) of a yaw rate relating to steering, the estimation value γH is fed, this delayed estimation value γH is further delayed by a time corresponding to a dynamic characteristic Gγ(s) of the yaw rate sensor, to be made to contribute to calculation of the deviation γe .

Description

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

【0001】[0001]

【発明の属する技術分野】本発明は、車両の横滑り角
や、横速度など、直接検出するのが困難な車両運動状態
量を推定するための装置に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an apparatus for estimating a vehicle motion state quantity which is difficult to directly detect, such as a side slip angle and a lateral speed of a vehicle.

【0002】[0002]

【従来の技術】車両の操舵角制御装置や、スタビリティ
ー制御装置は、検出が容易なヨーレート等の車両運動状
態量は直接これを検出し、また検出が容易でない車両の
横滑り角や、横速度などの車両運動状態量は、上記容易
に検出した車両運動状態量と、操舵角および車速などを
用いてこれを推定し、当該検出および推定した車両運動
状態量を制御情報として用いるため、運動状態量の推定
装置を必要とする。
2. Description of the Related Art A vehicle steering angle control device and a stability control device directly detect a vehicle motion state quantity such as a yaw rate, which can be easily detected. The vehicle motion state quantity, such as the vehicle motion state quantity, is estimated using the easily detected vehicle motion state quantity, the steering angle and the vehicle speed, and the detected and estimated vehicle motion state quantity is used as control information. Requires a quantity estimator.

【0003】車両の運動状態量を推定する装置としては
従来、例えば特開平2−45802号公報や、特開平8
−5650号公報に記載されたものがある。前者の公報
に記載の従来装置は、車両の車輪操舵角およびヨー角加
速度を検出し、これら検出した操舵角およびヨー角加速
度を用いた車両の運動方程式に基づくオブザーバ手法に
より車両のヨーレートおよび横加速度を推定するもので
ある。
[0003] As a device for estimating the motion state quantity of a vehicle, for example, Japanese Patent Application Laid-Open No. 2-45802 and Japanese Patent Application Laid-Open No.
There is one described in JP-A-5650. The conventional device described in the former publication detects a wheel steering angle and a yaw angular acceleration of a vehicle, and uses an observer method based on a vehicle equation of motion using the detected steering angle and yaw angular acceleration to obtain a yaw rate and a lateral acceleration of the vehicle. Is estimated.

【0004】また、後者の公報に記載された従来の車両
運動状態量推定装置は、車輪操舵角を入力変数とし、横
滑り角およびヨーレートを状態変数とし、ヨーレートを
参照入力とするオブザーバ手法を用いて車両の横滑り角
を推定するに際し、参照入力としてのヨーレートを検出
する際に検出信号中に含まれているセンサのノイズや、
ドリフト成分を除去するためのフィルタの特性を上記オ
ブザーバの状態方程式に加味し、これにより車両横滑り
角の推定制度を高めるというものである。
The conventional vehicle motion state quantity estimating device described in the latter publication uses an observer method in which the wheel steering angle is used as an input variable, the sideslip angle and the yaw rate are used as state variables, and the yaw rate is used as a reference input. When estimating the side slip angle of the vehicle, when detecting the yaw rate as a reference input, the sensor noise included in the detection signal,
The characteristic of the filter for removing the drift component is added to the state equation of the observer, thereby increasing the estimation accuracy of the vehicle side slip angle.

【0005】ここで上記のオブザーバは、ヨーレート
(推定値γH )および横滑り角(推定値βH )を推定す
る場合について示すと、一般的に図5のように構成さ
れ、その状態方程式は以下の(1)式により、また出力
方程式は以下の(2)式によって、それぞれ表されるこ
とが知られている。 (d/dt)XH =(A−KC)XH +Kγ+BU …(1) YH =CXH …(2) 従って、
[0005] The above-mentioned observer is generally constructed as shown in FIG. 5 when estimating the yaw rate (estimated value γ H ) and the sideslip angle (estimated value β H ). It is known that the output equation is expressed by the following equation (1), and the output equation is expressed by the following equation (2). (d / dt) X H = (A−KC) X H + Kγ + BU (1) Y H = CX H (2)

【数1】 となる。(Equation 1) Becomes

【0006】ここでオブザーバの収束を考慮して、その
設定極p1 ,p2 を任意に定め、オブザーバのゲイン行
列Kを求めると、k1,k2 はそれぞれ次式のごとくにな
り、 k1={ p1 ・p2- a11( a11 -p1-p2) / a21 }+a12 …(5) k2=a11+a22−p1-p2 …(6) これにより図5のオブザーバを構成することができる。
Here, considering the convergence of the observer, the set poles p 1 and p 2 are arbitrarily determined, and the gain matrix K of the observer is obtained. Then, k 1 and k 2 are as follows. 1 = {p 1 · p 2 -a 11 (a 11 -p 1 -p 2 ) / a 21 } + a 12 (5) k 2 = a 11 + a 22 -p 1 -p 2 ... (6) The observer of FIG. 5 can be configured.

【0007】[0007]

【発明が解決しようとする課題】しかし従来の車両運動
状態量の推定手法では、車輪操舵があったら直ちにこれ
に対応する車両運動状態量の変化があるとの仮定のもと
で当該推定を行っていたから、以下の問題を生じてい
た。つまり車両は通常、操舵周波数に対するヨーレート
ゲイン特性および位相特性を示す図6により説明する
と、車輪タイヤの変形などに起因してタイヤの破線で示
すように実際は操舵があってから無駄時間Δtの後にヨ
ーレート等の運動状態量の変化を惹起する。
However, in the conventional method of estimating the vehicle motion state quantity, the estimation is performed on the assumption that there is a change in the vehicle motion state quantity corresponding to the wheel steering immediately after the wheel steering. Had the following problems. That is, the vehicle is usually described with reference to FIG. 6 showing the yaw rate gain characteristic and the phase characteristic with respect to the steering frequency. Cause a change in the amount of exercise state.

【0008】しかして従来は当該無駄時間を考慮してい
なかったことから、参照入力(図5では実ヨーレート
γ)に対して推定値(図5ではヨーレート推定値γH
が誤差を持ったものとなり、特に高い操舵周波数領域に
おいてオブザーバが不安定になり易いという問題を生ず
る。
Conventionally, the dead time has not been taken into account, so that the estimated value (the estimated yaw rate value γ H in FIG. 5) is used for the reference input (actual yaw rate γ in FIG. 5).
Has an error, which causes a problem that the observer tends to be unstable particularly in a high steering frequency region.

【0009】これがため、当該不安定になるのを避ける
意味合いにおいてオブザーバのゲインを高くすることが
できず、推定値を速やかに実際の真値に持ち来し得なか
った。従って、車速Vx =40km/hでの走行中に操舵角
θを±120度の範囲で図7に示すようにさせた場合に
つき説明すると、同図に実線で示す実際のヨーレート真
値γ、横速度真値Vy 、横滑り角真値βに対して、破線
で示すヨーレート推定値γH 、横速度推定値VyH、横滑
り角推定値βH が過渡期において1点鎖線で示すような
誤差を持ち、推定値の真値への収束性が悪くなるのを免
れない。よって、従来の手法により推定した車両運動状
態量を車両のスタビリティー制御などに用いると、制御
が不安定、若しくは不正確になり、狙い通りの制御効果
を達成し得ないのが実情であった。
As a result, the gain of the observer cannot be increased in the sense of avoiding the instability, and the estimated value cannot be immediately brought to the actual true value. Therefore, the case where the steering angle θ is made as shown in FIG. 7 in the range of ± 120 degrees during traveling at the vehicle speed V x = 40 km / h will be described. The actual yaw rate true value γ shown by the solid line in FIG. The yaw rate estimated value γ H , the lateral speed estimated value V yH , and the estimated side slip angle β H indicated by the broken line are different from the true lateral velocity V y and the true side slip angle β in the transient period as indicated by the one-dot chain line. It is inevitable that the convergence of the estimated value to the true value deteriorates. Therefore, when the vehicle motion state quantity estimated by the conventional method is used for the stability control of the vehicle, the control becomes unstable or inaccurate, and the actual control effect cannot be achieved. .

【0010】請求項1に記載の第1発明は、過渡期も含
めて車両運動状態量の推定値を速やかに実際の運動状態
量に一致させ得るようにして、フィードバックゲインを
高くすることができるようにし、もって上述の問題を解
消することを目的とする。
According to the first aspect of the present invention, the feedback gain can be increased by making it possible to quickly match the estimated value of the vehicle motion state quantity including the transition period with the actual motion state quantity. Thus, an object of the present invention is to solve the above-mentioned problem.

【0011】請求項2に記載の第2発明は、第1発明の
目的を更に確実に達成し得るようにすることを目的とす
る。
A second aspect of the present invention is to achieve the object of the first aspect more reliably.

【0012】請求項3に記載の第3発明は、路面摩擦係
数の如何にかかわらず第1発明の目的を確実に達成し得
るようにすることを目的とする。
A third object of the present invention is to ensure that the object of the first invention can be achieved irrespective of the coefficient of road friction.

【0013】請求項4に記載の第4発明は、車両運動状
態量の推定値と実際の真値との間における遅れ偏差を求
めるに先立って車両運動状態量を検出する手段の動特性
の如何にかかわらず第1発明の目的を確実に達成し得る
ようにすることを目的とする。
According to a fourth aspect of the present invention, the dynamic characteristic of the means for detecting the vehicle motion state quantity prior to obtaining the delay deviation between the estimated value of the vehicle motion state quantity and the actual true value is determined. Regardless, the object of the first invention is to be achieved reliably.

【0014】請求項5に記載の第5発明は、オブザーバ
を用いて車両の運動状態量としての横滑り角を推定する
のに好適な車両運動状態量推定装置を提案することを目
的とする。
A fifth aspect of the present invention provides a vehicle motion state quantity estimating apparatus suitable for estimating a side slip angle as a vehicle motion state quantity using an observer.

【0015】[0015]

【課題を解決するための手段】これらの目的の第1発明
による車両の運動状態量推定装置は、操舵角と、車速
と、或る車両運動状態量の検出値とを用いて、該或る車
両運動状態量および他の車両運動状態量を推定するよう
にした装置において、操舵角と、車速と、前記或る車両
運動状態量の検出値および推定値間における遅れ偏差
と、車両モデルとを用いて、前記或る車両運動状態量お
よび他の車両運動状態量を推定し、該或る車両運動状態
量の推定値を所定時間遅らせて前記遅れ偏差の演算に資
するよう構成したことを特徴とするものである。
According to a first aspect of the present invention, there is provided an apparatus for estimating a motion state of a vehicle using a steering angle, a vehicle speed, and a detected value of a certain vehicle motion state. In a device for estimating a vehicle motion state quantity and another vehicle motion state quantity, a steering angle, a vehicle speed, a delay deviation between a detected value and an estimated value of the certain vehicle motion state quantity, and a vehicle model are calculated. The vehicle motion state quantity and the other vehicle motion state quantity are estimated, and the estimated value of the certain vehicle motion state quantity is delayed for a predetermined time to contribute to the calculation of the delay deviation. Is what you do.

【0016】第2発明による車両の運動状態量推定装置
は、第1発明において、車輪操舵に対する前記或る車両
運動状態量の遅れ時間を無駄時間として予め測定してお
き、前記或る車両運動状態量の推定値をこの無駄時間に
相当する所定時間だけ遅らせて前記遅れ偏差の演算に資
するよう構成したことを特徴とするものである。
According to a second aspect of the present invention, in the first aspect, the delay time of the certain vehicle motion state quantity with respect to wheel steering is measured in advance as a dead time, and the certain vehicle motion state quantity is estimated. The method is characterized in that the estimated value of the amount is delayed by a predetermined time corresponding to the dead time to contribute to the calculation of the delay deviation.

【0017】第3発明による車両の運動状態量推定装置
は、第2発明において、前記無駄時間を路面摩擦係数が
低いほど長くしたことを特徴とするものである。
A third aspect of the present invention is directed to the vehicle motion state estimating apparatus according to the second aspect, wherein the dead time is increased as the road surface friction coefficient decreases.

【0018】第4発明による車両の運動状態量推定装置
は、第2発明または第3発明において、前記無駄時間だ
け遅らせた或る車両運動状態量の推定値を、該或る車両
運動状態量の検出に用いた車両運動状態量検出手段の動
特性に対応する時間だけ更に遅らせて前記遅れ偏差の演
算に資するように構成したことを特徴とするものであ
る。
According to a fourth aspect of the present invention, there is provided the vehicle motion state quantity estimating apparatus according to the second or third invention, wherein the estimated value of the certain vehicle motion state quantity delayed by the dead time is replaced with the certain vehicle motion state quantity. The present invention is characterized in that it is configured to further delay by a time corresponding to the dynamic characteristic of the vehicle motion state amount detection means used for the detection and to contribute to the calculation of the delay deviation.

【0019】第5発明による車両の運動状態量推定装置
は、第1発明乃至第4発明のいずれかにおいて、前記或
る車両運動状態量として車両のヨーレートを、また、前
記他の車両運動状態量として車両の横滑り角を状態変数
とし、ヨーレートを参照入力とするオブザーバを構成
し、このオブザーバを用いて車両の横滑り角を推定する
よう構成したことを特徴とするものである。
According to a fifth aspect of the present invention, there is provided the vehicle motion state quantity estimation device according to any one of the first to fourth inventions, wherein the certain vehicle motion state quantity is a yaw rate of the vehicle, and the other vehicle motion state quantity. , An observer using the vehicle's side slip angle as a state variable and a yaw rate as a reference input is configured, and the side slip angle of the vehicle is estimated using the observer.

【0020】[0020]

【発明の効果】第1発明においては、操舵角と、車速
と、或る車両運動状態量の検出値とを用いて、当該或る
車両運動状態量および他の車両運動状態量を推定するに
際し、操舵角と、車速と、上記或る車両運動状態量の検
出値および推定値間における遅れ偏差と、車両モデルと
を用いて、当該或る車両運動状態量および他の車両運動
状態量を推定し、合わせて、上記或る車両運動状態量の
推定値を所定時間遅らせて上記遅れ偏差の演算に資する
ことから、車両が実際は操舵から無駄時間の遅れをもっ
て対応した挙動変化を生ずるのに符合して、車両運動状
態量の推定値を実際値に一層良く一致させ得ることとな
る。
According to the first aspect of the invention, when a certain vehicle motion state quantity and another vehicle motion state quantity are estimated using a steering angle, a vehicle speed, and a certain vehicle motion state quantity detection value. Estimating the certain vehicle motion state quantity and the other vehicle motion state quantity using the steering angle, the vehicle speed, the delay deviation between the detected value and the estimated value of the certain vehicle motion state quantity, and the vehicle model. In addition, since the estimation value of the certain vehicle motion state amount is delayed for a predetermined time to contribute to the calculation of the delay deviation, it is concluded that the vehicle actually produces a corresponding behavior change with a delay of dead time from steering. As a result, the estimated value of the vehicle motion state quantity can be made to better match the actual value.

【0021】これがため第1発明においては、高い操舵
周波数領域においても上記の推定が不安定になることが
なく、これを懸念することなしに当該推定に際してその
ゲインを高くすることができ、推定値の真値への収束性
を向上させ得る。
Therefore, in the first invention, the above estimation does not become unstable even in a high steering frequency range, and the gain can be increased at the time of the estimation without worrying about this. Can be improved to the true value.

【0022】第2発明においては、車輪操舵に対する前
記或る車両運動状態量の遅れ時間を無駄時間として予め
測定しておき、前記或る車両運動状態量の推定値をこの
無駄時間に相当する所定時間だけ遅らせて前記遅れ偏差
の演算に資することから、上記の所定時間が車両運動状
態量の遅れ時間に対応して、一層実情にマッチしたもの
なり、第1発明の目的を更に確実に達成することができ
る。
In the second invention, the delay time of the certain vehicle motion state amount with respect to wheel steering is measured in advance as a dead time, and the estimated value of the certain vehicle motion state amount is determined as a predetermined value corresponding to the dead time. Since the time is delayed by a time to contribute to the calculation of the delay deviation, the above-mentioned predetermined time is more matched to the actual situation in accordance with the delay time of the vehicle motion state quantity, and the object of the first invention is more reliably achieved. be able to.

【0023】第3発明においては、上記無駄時間を路面
摩擦係数が低いほど長くしたことから、低μ路面では車
輪操舵に対する車両運動状態量の発生遅れが大きくなる
のに符合することとなり、路面摩擦係数の如何にかかわ
らず第1発明の目的を確実に達成することができる。
In the third aspect of the present invention, since the above-mentioned dead time is increased as the road surface friction coefficient becomes lower, the delay of the vehicle motion state amount with respect to wheel steering becomes larger on a low μ road surface. The object of the first invention can be reliably achieved regardless of the coefficient.

【0024】第4発明においては、上記無駄時間だけ遅
らせた或る車両運動状態量の推定値を、該或る車両運動
状態量の検出に用いた車両運動状態量検出手段の動特性
に対応する時間だけ更に遅らせて前記遅れ偏差の演算に
資することから、当該車両運動状態量検出手段の動特性
の如何にかかわらず第1発明の目的を確実に達成するこ
とができる。
In the fourth invention, the estimated value of the certain vehicle motion state amount delayed by the dead time corresponds to the dynamic characteristic of the vehicle motion state amount detection means used for detecting the certain vehicle motion state amount. Since the delay is further delayed by a time to contribute to the calculation of the delay deviation, the object of the first invention can be reliably achieved regardless of the dynamic characteristics of the vehicle motion state amount detecting means.

【0025】第5発明においては、前記或る車両運動状
態量として車両のヨーレートを、また、前記他の車両運
動状態量として車両の横滑り角を状態変数とし、ヨーレ
ートを参照入力とするオブザーバを構成し、このオブザ
ーバを用いて車両の横滑り角を推定することから、オブ
ザーバを用いて横滑り角を推定する、実際的で、且つ、
好適な車両運動状態量推定装置を構築することができ
る。
According to a fifth aspect of the present invention, there is provided an observer which uses the yaw rate of the vehicle as the certain vehicle motion state quantity, the vehicle side slip angle as the other vehicle motion state quantity, and the yaw rate as a reference input. Then, since the sideslip angle of the vehicle is estimated using the observer, the sideslip angle is estimated using the observer, and it is practical and
A suitable vehicle motion state amount estimation device can be constructed.

【0026】[0026]

【発明の実施の形態】以下、本発明の実施の形態を図面
に基づき詳細に説明する。図1は、本発明による車両の
運動状態量推定装置をマイクロコンピュータ10で構成
した一実施の形態を示し、この場合、マイクロコンピュ
ータ10には、ステアリングホイールの操舵角θを検出
する操舵角センサ11からの信号、車速V x を検出する
車速センサ12からの信号、ヨーレートγを検出するヨ
ーレートセンサ13からの信号、および路面摩擦係数μ
を検出する路面μセンサ14(雨天時にONするワイパ
スイッチからの信号でもよい)からの信号を、それぞれ
の信号処理回路15〜18を経て入力する。なお各信号
処理回路15〜18は、入力されてくる検出信号をフィ
ルタリング処理およびディジタル変換処理してマイクロ
コンピュータ10に入力するものとする。
Embodiments of the present invention will be described below with reference to the drawings.
This will be described in detail based on FIG. FIG. 1 shows a vehicle according to the invention.
Exercise state amount estimation device is composed of microcomputer 10
In this case, a micro computer
The steering angle θ of the steering wheel
Signal from the steering angle sensor 11 xDetect
A signal from the vehicle speed sensor 12, a yaw rate γ
-Signal from the rate sensor 13 and coefficient of road friction μ
Road surface sensor 14 (wiper that turns ON in rainy weather)
Signal from the switch).
Through the signal processing circuits 15-18. Each signal
The processing circuits 15 to 18 filter the inputted detection signals.
Filtering and digital conversion
It is assumed to be input to the computer 10.

【0027】マイクロコンピュータ10はこれら入力情
報をもとに、図2に示す機能別ブロック線図および図3
のプログラムに沿う、詳しくは後述する演算により横滑
り角およびヨーレートを推定し、これら横滑り角推定値
βH およびヨーレート推定値γH に関する信号をそれぞ
れ、D/A変換器19,20によりアナログ信号に変換
し、その後、増幅器21,22により増幅して出力する
ものとする。
Based on the input information, the microcomputer 10 uses the functional block diagram shown in FIG.
In accordance with the program, the slip angle and the yaw rate are estimated by the calculation described later in detail, and signals relating to the estimated slip angle β H and the estimated yaw rate γ H are converted into analog signals by D / A converters 19 and 20, respectively. Thereafter, the signals are amplified by the amplifiers 21 and 22 and output.

【0028】ここでマイクロコンピュータ10は、図2
に示すごとく同一次元オブザーバを機能的に構成するよ
うになし、車速Vx および操舵角θ(ステアリングギヤ
比をNとすると、前輪舵角δF はδF =θ/Nで表され
る)により決定される(但し、後輪操舵を行わず、後輪
舵角δR は常に0とする)車両の運動状態量としてヨー
レートγおよび横滑り角βを推定することとする。先ず
推定理論を説明するに、図2のごとく入力変数をUと
し、状態変数をXHとすると、これらの間には車両の運
動方程式から、車両モデルを表すシステム行列A、観測
行列Bを用いて、 (d/dt)XH =AXH +BU …(7) の関係があることが知られている。
Here, the microcomputer 10 corresponds to FIG.
As shown in FIG. 6, the same-dimensional observer is functionally constructed, and the vehicle speed V x and the steering angle θ (where the steering gear ratio is N, the front wheel steering angle δ F is represented by δ F = θ / N) The yaw rate γ and the side slip angle β are estimated as the motion state quantities of the vehicle to be determined (however, the rear wheel steering is not performed and the rear wheel steering angle δ R is always 0). First, explaining the estimation theory, the as input variables and U 2, when the state variable and X H, using the motion equation of the vehicle between them, the system matrix A which represents the vehicle model, the observation matrix B Therefore, it is known that there is a relationship of (d / dt) X H = AX H + BU (7).

【0029】ところで入力変数Uが前輪舵角δF であ
り、状態変数XH が横滑り角βおよびヨーレートγであ
ることから、車両の運動状態方程式は次式により表され
る。
Since the input variable U is the front wheel steering angle δ F and the state variables X H are the sideslip angle β and the yaw rate γ, the vehicle motion state equation is expressed by the following equation.

【数2】 ただし、システム行列の各要素及び観測行列の各要素は
下記のとおりである。
(Equation 2) Here, each element of the system matrix and each element of the observation matrix are as follows.

【数3】 (Equation 3)

【数4】 M:車両質量 Iz:ヨー慣性モーメント KF :前輪等価コーナリングパワー KR :後輪等価コーナリングパワー LF :車両重心から前輪車軸までの距離 LR :車両重心から後輪車軸までの距離(Equation 4) M: vehicle mass Iz: yaw moment of inertia K F : front-wheel equivalent cornering power K R : rear-wheel equivalent cornering power L F : distance from vehicle center of gravity to front wheel axle L R : distance from vehicle center of gravity to rear wheel axle

【0030】従って、上記(7)式の状態方程式に基づ
き、ヨーレートγを観測可能な状態変数(参照入力)と
するオブザーバを構成し、車両諸元を予め測定した既知
の定数として設定し、車速Vx 、操舵角θ(前輪舵角δ
F )およびヨーレートγを検出して与えれば、上記のオ
ブザーバを用いてヨーレート推定値γH および横滑り角
推定値βH を求めることができる。
Therefore, based on the state equation of the above equation (7), an observer is set up using the yaw rate γ as a state variable (reference input) that can be observed, and the vehicle parameters are set as known constants measured in advance, and the vehicle speed is set. V x , steering angle θ (front wheel steering angle δ
If F ) and the yaw rate γ are detected and given, the estimated yaw rate γ H and the estimated side slip angle β H can be obtained using the above-described observer.

【0031】ところで本実施の形態においては特に、観
測行列Cを経由したヨーレート推定値γH と、実際のヨ
ーレート検出値γとの間における遅れ偏差γe をフィー
ドバックゲインKで、積分器1/sの前段に戻すように
し、これにより操舵角δF と、車速Vx と、上記のヨー
レート遅れ偏差γe と、車両モデルAとを用いて、ヨー
レート推定値γH および横滑り角推定値βH を求めるこ
ととする。
In this embodiment, the delay deviation γ e between the estimated yaw rate value γ H via the observation matrix C and the actual detected yaw rate value γ is represented by the feedback gain K and the integrator 1 / s of the back to front, and thereby the steering angle [delta] F, and the vehicle speed V x and the yaw rate delay deviation gamma e above, by using the vehicle model a, the yaw rate estimated value gamma H and sideslip angle estimates beta H I will ask for it.

【0032】ところで本実施の形態においては更に、ヨ
ーレート遅れ偏差γe を算出するに際してヨーレート推
定値γH をそのまま当該演算に資することとせず、図6
につき前述した、車輪操舵に対するヨーレートの発生遅
れ時間Δtに相当する無駄時間 exp(-ls) を予め測定し
ておき、この無駄時間exp(-ls)だけ遅らせてヨーレート
推定値γH を送出する。そしてこの遅らせたヨーレート
推定値γH を、ヨーレートセンサ13(図1参照)の動
特性Gγ(s)に相当するフィルタに通して、当該動特
性Gγ(s)に対応する時間だけ更に遅らせて上記ヨー
レート遅れ偏差γe の演算に資する。
Incidentally, in the present embodiment, when calculating the yaw rate delay deviation γ e , the estimated yaw rate γ H is not directly used for the calculation.
As described above, the dead time exp (-ls) corresponding to the yaw rate generation delay time Δt for the wheel steering is measured in advance, and the yaw rate estimated value γ H is transmitted with a delay by the dead time exp (-ls). Then, the delayed yaw rate estimated value γ H is passed through a filter corresponding to the dynamic characteristic Gγ (s) of the yaw rate sensor 13 (see FIG. 1), and further delayed by the time corresponding to the dynamic characteristic Gγ (s). This contributes to the calculation of the yaw rate delay deviation γ e .

【0033】なお本実施の形態においては、無駄時間修
正部により無駄時間 exp(-ls) を路面摩擦係数μが小さ
くなるにつれて長くするように修正する。その理由は、
路面摩擦係数μが小さくなるにつれて車輪タイヤの変形
が起こりにくく、結果として路面摩擦係数μが小さくな
るにつれ、車輪操舵に対するヨーレートの発生遅れ時間
が大きくなるからである。
In this embodiment, the dead time correcting section corrects the dead time exp (-ls) so as to be longer as the road friction coefficient μ becomes smaller. The reason is,
This is because, as the road surface friction coefficient μ decreases, deformation of the wheel tires hardly occurs. As a result, as the road surface friction coefficient μ decreases, the delay time of the yaw rate generation with respect to wheel steering increases.

【0034】以上の構成になる本実施の形態において運
動状態方程式は、ヨーレート遅れ偏差γe をゲインKで
フィードバックする系が前記(7)式に付加されて、次
式により表される。 (d/dt)XH =AXH −K〔γH ・ exp(-ls) ・Gγ(s)−γ〕 +BU …(11)
In the present embodiment having the above-described configuration, the motion state equation is represented by the following equation with a system for feeding back the yaw rate delay deviation γ e with the gain K added to the equation (7). (d / dt) X H = AX H -K [γ H · exp (-ls) · Gγ (s) -γ] + BU (11)

【0035】つまり本実施の形態においてマイクロコン
ピュータ10(図1参照)は、ヨーレート推定値γ
H を、路面摩擦係数μの低下につれて長くなるよう修正
した無駄時間 exp(-ls) だけ遅らせ、この遅らせたヨー
レート推定値γH を、ヨーレートセンサ13(図1参
照)の動特性Gγ(s)に対応する時間だけ更に遅らせ
て、ヨーレート検出値γと比較する。そして、これらヨ
ーレート推定値γH とヨーレート検出値γとの間におけ
るヨーレート遅れ偏差γe を演算し、これをゲインKで
フィードバックしてヨーレート推定値γH および横滑り
角推定値βH を求める。なお、行列A,B,Kは車速V
x の関数としてメモリに記憶しておき、必要な時に読み
出して使用すること勿論である。
That is, in the present embodiment, the microcomputer 10 (see FIG. 1)
H is delayed by a dead time exp (-ls) corrected so as to become longer as the road surface friction coefficient μ decreases, and the delayed yaw rate estimation value γ H is calculated as the dynamic characteristic Gγ (s) Is further delayed by a time corresponding to the above, and is compared with the yaw rate detection value γ. Then, a yaw rate delay deviation γ e between the estimated yaw rate value γ H and the detected yaw rate value γ is calculated, and this is fed back with a gain K to obtain an estimated yaw rate value γ H and an estimated side slip angle value β H. Note that the matrices A, B, and K represent the vehicle speed V
Needless to say, it is stored in a memory as a function of x and read out and used when necessary.

【0036】これがため、車輪操舵から車両運動状態量
が変化するまでに無駄時間 exp(-ls) があっても、ま
た、この無駄時間が路面摩擦係数μに応じて変化して
も、更に、ヨーレート検出値γがヨーレートセンサ13
の動特性Gγ(s)により遅れることがあっても、図7
と同じ条件で行った図4のシミュレーション結果から明
らかなように、ヨーレート推定値γH を絶えず実際の検
出値γにほぼ一致させることができると共に、横速度推
定値VyHおよび車速Vx から(VyH/Vx )の演算によ
り求める横滑り角推定値βH も実際の横滑り角推定値β
にほぼ一致させることができる。
Therefore, even if there is a dead time exp (-ls) from the wheel steering to the change of the vehicle motion state quantity, and even if this dead time changes according to the road surface friction coefficient μ, The yaw rate detection value γ is the yaw rate sensor 13
7 may be delayed due to the dynamic characteristic Gγ (s) of FIG.
As is clear from the simulation result of FIG. 4 performed under the same conditions as in FIG. 4, the estimated yaw rate value γ H can be constantly made to substantially match the actual detected value γ, and the estimated yaw rate value V yH and the vehicle speed V x ( V yH / V x ) is also calculated as the actual slip angle estimated value β H
Can be almost matched.

【0037】マイクロコンピュータ10が行う機能を図
2では便宜上ブロック線図として示したが、実際は一定
時間ΔTごとの定時割り込みにより実行する図3のよう
な制御プログラムである。以下このフローチャートを説
明するに、ステップ31においては、操舵角θ、車速V
x 、ヨーレートγ、および路面摩擦係数μを読み込み、
これらを今回の操舵角検出値θ(i)、車速検出値Vx
(i)、ヨーレート検出値γ(i)、および路面摩擦係
数検出値μ(i)としてメモリする。
Although the functions performed by the microcomputer 10 are shown in FIG. 2 as a block diagram for convenience, it is actually a control program as shown in FIG. 3 which is executed by a periodic interruption every predetermined time ΔT. The flowchart will be described below. In step 31, the steering angle θ and the vehicle speed V
x , yaw rate γ, and road friction coefficient μ,
These are used as the current detected steering angle θ (i) and the detected vehicle speed V x
(I), the yaw rate detection value γ (i), and the road surface friction coefficient detection value μ (i) are stored in the memory.

【0038】次いでステップ32において、入力変数U
(i)に操舵角検出値θ(i)を代入すると共に、ヨー
レート推定値γH を、路面摩擦係数検出値μ(i)に応
じて前記のごとく修正した無駄時間 exp(-ls) だけ遅ら
せて得られるヨーレート推定値の遅れ演算値γH (i−
n)を、ヨーレートセンサ13(図1参照)の動特性G
γ(s)に対応する時間だけ更に遅らせて求めたγ
H (i−n)・Gγ(s)に対するヨーレート検出値γ
(i)の偏差、つまり前記のヨーレート遅れ偏差γ
e (i)を演算する。
Next, at step 32, the input variable U
(I) is substituted for the detected steering angle θ (i),
Rate estimate γHAccording to the road friction coefficient detection value μ (i).
Delay by exp (-ls)
Calculated value γ of the estimated yaw rate valueH(I-
n) is the dynamic characteristic G of the yaw rate sensor 13 (see FIG. 1).
γ obtained by further delaying by the time corresponding to γ (s)
H(In) · Gγ (s), the detected yaw rate value γ
(I), that is, the yaw rate delay deviation γ
e(I) is calculated.

【0039】次にステップ33において、車速検出値V
x (i)に応じたA行列をテーブル1の参照により作成
し、ステップ34において、車速検出値Vx (i)に応
じたB行列をテーブル2の参照により作成し、ステップ
35において、車速検出値Vx (i)に応じたK行列を
テーブル3の参照により作成する。
Next, at step 33, the detected vehicle speed V
x The A matrix in accordance with (i) prepared by reference to Table 1, in step 34, the B matrix corresponding to the vehicle speed detection value V x (i) prepared by reference to Table 2, in step 35, the vehicle speed detection A K matrix corresponding to the value V x (i) is created by referring to Table 3.

【0040】次のステップ36では、前記(11)式に
基づくオブザーバ行列 (d/dt)XH (i)=AXH (i−1)−Kγe (i)+BU(i) を作成する。このとき状態変数ベクトルXH は、前回の
演算周期における、つまりΔT時間前における(i−
1)回目の演算値XH (i−1)を用いる。ここでXH
は、横滑り角推定値βH およびヨーレート推定値γH
らなる状態変数ベクトルである。更にステップ37で、
前記(11)式の演算を行って状態変数ベクトルX
H (i)の微分値(d/dt)XH (i)を求める。
In the next step 36, an observer matrix (d / dt) X H (i) = AX H (i−1) −Kγ e (i) + BU (i) based on the above equation (11) is created. At this time, the state variable vector XH is calculated in the previous calculation cycle, that is, (i−
1) Use the calculated value X H (i-1) for the first time. Where X H
Is a state variable vector including the estimated sideslip angle value β H and the estimated yaw rate value γ H. Further, in step 37,
By performing the operation of the above equation (11), the state variable vector X
Differential value of H (i) the (d / dt) X H ( i) determined.

【0041】ステップ38では、ステップ37で演算し
た状態変数ベクトルの微分値(d/dt)XH (i)を積分す
る、 XH (i)=∫(d/dt)XH (i)・dt …(12) の演算を行って、今回の横滑り角推定値βH (i)およ
びヨーレート推定値γH(i)からなる状態変数の推定
値XH (i)を算出する。具体的には、1演算周期ΔT
前の(i−1)回目における推定値XH (i−1)を用
い、今回のi回目の演算の推定値を XH (i)=XH (i−1)+ΔT・(d/dt)XH (i) …(13) の式から求める。
In step 38, the differential value (d / dt) X H (i) of the state variable vector calculated in step 37 is integrated. X H (i) = ∫ (d / dt) X H (i) · The calculation of dt (12) is performed to calculate the estimated value X H (i) of the state variable composed of the current estimated slip angle β H (i) and the estimated yaw rate γ H (i). Specifically, one operation cycle ΔT
Using previous the (i-1) times the estimated value at the first X H (i-1), the estimated value X H of this i-th operation (i) = X H (i -1) + ΔT · (d / dt ) X H (i) Calculated from the equation (13).

【0042】ステップ39においては、次回の演算のた
めに今回の状態変数推定値XH (i)をXH (i−1)
に更新するが、本実施の形態においては状態変数推定値
H(i)のうちヨーレート推定値γH (i)を所定時
間遅らせるために、γH (i)をメモリに格納する。そ
してステップ40で、所定時間(ΔT・n)だけ遅れた
γH (i−n)を算出する処理を行い、当該算出値をメ
モリに格納する。最後にステップ41において、状態変
数推定値XH (i)、すなわち横滑り角推定値β
H (i)およびヨーレート推定値γH (i)に応じた信
号を、図1のD/A変換器19,20および増幅器2
1,22を順次経て出力する。
In step 39, the current state variable estimated value X H (i) is converted to X H (i-1) for the next calculation.
To update, but to retard yaw rate estimated value gamma H a (i) a predetermined time of the state variable estimates X H (i) in the present embodiment, stores gamma H (i) to the memory. Then, in step 40, a process of calculating γ H (in) delayed by a predetermined time (ΔT · n) is performed, and the calculated value is stored in the memory. Finally, in step 41, the estimated state variable X H (i), ie, the estimated slip angle β
H (i) and a signal corresponding to the yaw rate estimated value γ H (i) are converted into the D / A converters 19 and 20 and the amplifier 2 shown in FIG.
The data is sequentially output through 1 and 22.

【0043】以上の処理は前記した通り一定時間ΔTご
との定時割り込みにより行うことから、この演算周期Δ
Tを適切に設定することによりヨーレート推定値γ
H (i)を遅らせる所定時間を自由に設定することがで
き、これによりヨーレート推定値γH (i)を、車輪操
舵から車両運動状態量が変化するまでに無駄時間 exp(-
ls) があっても、また、この無駄時間が路面摩擦係数μ
に応じて変化しても、更に、ヨーレート検出値γがヨー
レートセンサ13の動特性Gγ(s)により遅れること
があっても、図4のシミュレーションの結果から明らか
なように実際のヨーレートγ(i)の発生にほぼ調時さ
せることができる。
Since the above processing is performed by the periodic interruption every predetermined time ΔT as described above, this calculation cycle Δ
By setting T appropriately, the estimated yaw rate γ
The predetermined time for delaying H (i) can be set freely, whereby the estimated yaw rate value γ H (i) is changed to a dead time exp (−) from the wheel steering to the change in the vehicle motion state quantity.
ls), this dead time is also the coefficient of road friction μ
Even if the yaw rate detection value γ is further delayed by the dynamic characteristic Gγ (s) of the yaw rate sensor 13 as shown in the simulation result of FIG. ) Can be almost timed.

【0044】これがため横滑り角の推定値βH (i)
も、同じく図4に示すように実際の検出値γ(i)にほ
ぼ一致させることができ、これがため、高い操舵周波数
領域においても上記の推定が不安定になることがなく、
これを懸念することなしに当該推定に際してそのゲイン
を高くすることができ、推定値の真値への収束性を向上
させることができる。
As a result, the estimated value of the sideslip angle β H (i)
Also, as shown in FIG. 4, it is possible to substantially match the actual detection value γ (i), so that the above estimation does not become unstable even in a high steering frequency region.
The gain can be increased during the estimation without worrying about this, and the convergence of the estimated value to the true value can be improved.

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

【図1】本発明一実施の形態になる車両の運動状態量推
定装置を示すブロック線図である。
FIG. 1 is a block diagram showing a vehicle motion state quantity estimation device according to an embodiment of the present invention.

【図2】同一実施の形態においてマイクロコンピュータ
が実行する運動状態量推定機能の機能別ブロック線図で
ある。
FIG. 2 is a functional block diagram of a motion state amount estimating function executed by a microcomputer in the same embodiment.

【図3】同一実施の形態においてマイクロコンピュータ
が実行する運動状態量推定制御のフローチャートであ
る。
FIG. 3 is a flowchart of exercise state amount estimation control executed by the microcomputer in the same embodiment.

【図4】同一実施の形態において運動状態量を推定した
場合における各種検出信号および推定信号の変化タイム
チャートである。
FIG. 4 is a change time chart of various detection signals and estimation signals when a motion state amount is estimated in the same embodiment.

【図5】従来の一般的な車両の運動状態量推定装置を示
すブロック線図である。
FIG. 5 is a block diagram showing a conventional general vehicle motion state amount estimation device.

【図6】車両のヨーレートゲイン特性および位相特性を
示す線図である。
FIG. 6 is a diagram showing a yaw rate gain characteristic and a phase characteristic of a vehicle.

【図7】従来の運動状態量推定装置により車両の運動状
態量を推定した場合における各種検出信号および推定信
号の変化タイムチャートである。
FIG. 7 is a change time chart of various detection signals and estimated signals when a motion state amount of a vehicle is estimated by a conventional motion state amount estimation device.

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

10 マイクロコンピュータ 11 操舵角センサ 12 車速センサ 13 ヨーレートセンサ 14 路面μセンサ 15 信号処理回路 16 信号処理回路 17 信号処理回路 18 信号処理回路 19 D/A 変換器 20 D/A 変換器 21 増幅器 22 増幅器 10 Microcomputer 11 Steering angle sensor 12 Vehicle speed sensor 13 Yaw rate sensor 14 Road surface μ sensor 15 Signal processing circuit 16 Signal processing circuit 17 Signal processing circuit 18 Signal processing circuit 19 D / A converter 20 D / A converter 21 Amplifier 22 Amplifier

Claims (5)

【特許請求の範囲】[Claims] 【請求項1】 操舵角と、車速と、或る車両運動状態量
の検出値とを用いて、該或る車両運動状態量および他の
車両運動状態量を推定するようにした装置において、 操舵角と、車速と、前記或る車両運動状態量の検出値お
よび推定値間における遅れ偏差と、車両モデルとを用い
て、前記或る車両運動状態量および他の車両運動状態量
を推定し、 該或る車両運動状態量の推定値を所定時間遅らせて前記
遅れ偏差の演算に資するよう構成したことを特徴とする
車両の運動状態量推定装置。
1. An apparatus for estimating a certain vehicle motion state quantity and another vehicle motion state quantity using a steering angle, a vehicle speed, and a detected value of a certain vehicle motion state quantity. Angle, vehicle speed, a delay deviation between the detected value and the estimated value of the certain vehicle motion state amount, and using a vehicle model, estimating the certain vehicle motion state amount and other vehicle motion state amounts, An apparatus for estimating a motion state amount of a vehicle, wherein the estimation value of the certain vehicle motion state amount is delayed for a predetermined time to contribute to the calculation of the delay deviation.
【請求項2】 請求項1において、車輪操舵に対する前
記或る車両運動状態量の遅れ時間を無駄時間として予め
測定しておき、前記或る車両運動状態量の推定値をこの
無駄時間に相当する所定時間だけ遅らせて前記遅れ偏差
の演算に資するよう構成したことを特徴とする車両の運
動状態量推定装置。
2. The method according to claim 1, wherein a delay time of the certain vehicle motion state quantity with respect to wheel steering is measured in advance as a dead time, and the estimated value of the certain vehicle motion state quantity corresponds to the dead time. An apparatus for estimating a motion state amount of a vehicle, wherein the apparatus is configured to contribute to the calculation of the delay deviation by delaying by a predetermined time.
【請求項3】 請求項2において、前記無駄時間を路面
摩擦係数が低いほど長くしたことを特徴とする車両の運
動状態量推定装置。
3. The vehicle motion state estimating device according to claim 2, wherein the dead time is made longer as the road surface friction coefficient becomes lower.
【請求項4】 請求項2または3において、前記無駄時
間だけ遅らせた或る車両運動状態量の推定値を、該或る
車両運動状態量の検出に用いた車両運動状態量検出手段
の動特性に対応する時間だけ更に遅らせて前記遅れ偏差
の演算に資するよう構成したことを特徴とする車両の運
動状態量推定装置。
4. The dynamic characteristic of a vehicle motion state quantity detecting means according to claim 2, wherein an estimated value of the certain vehicle motion state quantity delayed by the dead time is used for detecting the certain vehicle motion state quantity. The vehicle motion state quantity estimating device is further configured to contribute to the calculation of the delay deviation by further delaying the time by a time corresponding to the following.
【請求項5】 請求項1乃至4のいずれか1項におい
て、前記或る車両運動状態量として車両のヨーレート
を、また、前記他の車両運動状態量として車両の横滑り
角を状態変数とし、ヨーレートを参照入力とするオブザ
ーバを構成し、このオブザーバを用いて車両の横滑り角
を推定するよう構成したことを特徴とする車両の運動状
態量推定装置。
5. The yaw rate according to claim 1, wherein a yaw rate of the vehicle is used as the certain vehicle motion state quantity, and a side slip angle of the vehicle is used as the other vehicle motion state quantity. A motion state estimating device for a vehicle, characterized in that the observer is configured to use the observer as a reference input, and the side slip angle of the vehicle is estimated using the observer.
JP10000164A 1998-01-05 1998-01-05 Kinetic condition quantity estimating device for vehicle Pending JPH11192959A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP10000164A JPH11192959A (en) 1998-01-05 1998-01-05 Kinetic condition quantity estimating device for vehicle

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP10000164A JPH11192959A (en) 1998-01-05 1998-01-05 Kinetic condition quantity estimating device for vehicle

Publications (1)

Publication Number Publication Date
JPH11192959A true JPH11192959A (en) 1999-07-21

Family

ID=11466403

Family Applications (1)

Application Number Title Priority Date Filing Date
JP10000164A Pending JPH11192959A (en) 1998-01-05 1998-01-05 Kinetic condition quantity estimating device for vehicle

Country Status (1)

Country Link
JP (1) JPH11192959A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006151071A (en) * 2004-11-26 2006-06-15 Toyota Motor Corp Device obtaining vehicle turn-traveling assist-yawing moment
JP2010167817A (en) * 2009-01-20 2010-08-05 Toyota Central R&D Labs Inc Steering control device and program
US8041480B2 (en) 2007-02-01 2011-10-18 Toyota Jidosha Kabushiki Kaisha Apparatus for estimating yaw rate
JP2016199080A (en) * 2015-04-07 2016-12-01 株式会社日本自動車部品総合研究所 Parking support device

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006151071A (en) * 2004-11-26 2006-06-15 Toyota Motor Corp Device obtaining vehicle turn-traveling assist-yawing moment
JP4604685B2 (en) * 2004-11-26 2011-01-05 トヨタ自動車株式会社 Device for determining vehicle turning travel assist yaw moment
US8041480B2 (en) 2007-02-01 2011-10-18 Toyota Jidosha Kabushiki Kaisha Apparatus for estimating yaw rate
JP2010167817A (en) * 2009-01-20 2010-08-05 Toyota Central R&D Labs Inc Steering control device and program
JP2016199080A (en) * 2015-04-07 2016-12-01 株式会社日本自動車部品総合研究所 Parking support device

Similar Documents

Publication Publication Date Title
US8457832B2 (en) Vehicle body speed estimating device
JP3060923B2 (en) Vehicle state estimation device
KR960033976A (en) Method for estimating the amount of motion state that changes according to turning of a vehicle
JPH11101732A (en) Road surface frictional coefficient detector
JP2007182209A (en) Vehicle physical quantity estimation apparatus and program
JP2007216942A (en) Sideslip angle estimation apparatus, automobile, and sideslip angle estimation method
JP3236003B1 (en) Road surface friction coefficient estimation device for vehicles
JPH10267685A (en) Method for estimating side slide angle of vehicle
JP2008094375A (en) Vehicle physical quantity estimating apparatus and program
JPH08202403A (en) Vehicle state quantity estimation device
JPH11192959A (en) Kinetic condition quantity estimating device for vehicle
JP3271945B2 (en) Road surface friction coefficient estimation device for vehicles
JP4887721B2 (en) Vehicle running state estimation device
JP3022167B2 (en) Vehicle motion state estimation method
US20040153216A1 (en) Method for estimating a vehicle's velocity
JP3217700B2 (en) Vehicle side slip angle detection device
JPH08119130A (en) Rear wheel steering angle control device for four-wheel drive vehicle
JPH03195973A (en) Car body speed estimating apparatus
JPH0245802A (en) Estimating device for vehicle state quantity
JP3271955B2 (en) Road surface friction coefficient estimation device for vehicles
US20230182523A1 (en) System and method for estimating lateral acceleration
JP2962026B2 (en) Body slip angle estimation device
JP3417099B2 (en) Road surface friction coefficient state detector
JP2008024233A (en) Vehicle side slip angle calculation unit
JPH07304460A (en) Vehicle motion control device

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050913

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20051026

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20060221