JP4291863B2 - Method and apparatus for predicting intake pipe pressure of internal combustion engine - Google Patents

Method and apparatus for predicting intake pipe pressure of internal combustion engine Download PDF

Info

Publication number
JP4291863B2
JP4291863B2 JP2007322117A JP2007322117A JP4291863B2 JP 4291863 B2 JP4291863 B2 JP 4291863B2 JP 2007322117 A JP2007322117 A JP 2007322117A JP 2007322117 A JP2007322117 A JP 2007322117A JP 4291863 B2 JP4291863 B2 JP 4291863B2
Authority
JP
Japan
Prior art keywords
difference
intake pipe
pipe pressure
value
current
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
JP2007322117A
Other languages
Japanese (ja)
Other versions
JP2008075663A (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.)
Honda Motor Co Ltd
Original Assignee
Honda 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 Honda Motor Co Ltd filed Critical Honda Motor Co Ltd
Priority to JP2007322117A priority Critical patent/JP4291863B2/en
Publication of JP2008075663A publication Critical patent/JP2008075663A/en
Application granted granted Critical
Publication of JP4291863B2 publication Critical patent/JP4291863B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Combined Controls Of Internal Combustion Engines (AREA)

Description

本発明は、内燃機関の吸気管圧力予測方法および装置に関する。特にファジィアルゴリズムを使用した吸気管圧力予測方法およびファジィ推定器に関する。   The present invention relates to an intake pipe pressure prediction method and apparatus for an internal combustion engine. In particular, the present invention relates to an intake pipe pressure prediction method using a fuzzy algorithm and a fuzzy estimator.

内燃機関の効率的な燃焼を実現するように燃料噴射制御が行われる。図1は、本発明に関係する内燃機関の吸気部分の構成を示す。吸入空気は、スロットル1を介してシリンダに供給される。スロットル弁開度は、所望の値に制御されている。また、吸気管圧力は、センサ2により測定されている。ここで、燃料噴射制御を適切に行うためには、シリンダ内吸入空気量を推定する必要があり、そのために吸気管圧力を予測する必要がある。   Fuel injection control is performed so as to realize efficient combustion of the internal combustion engine. FIG. 1 shows a configuration of an intake portion of an internal combustion engine related to the present invention. The intake air is supplied to the cylinder via the throttle 1. The throttle valve opening is controlled to a desired value. The intake pipe pressure is measured by the sensor 2. Here, in order to appropriately perform the fuel injection control, it is necessary to estimate the intake air amount in the cylinder, and therefore it is necessary to predict the intake pipe pressure.

従来から、内燃機関の運転状態に基づいて現時点の吸気管圧力を予測すると共に、該予測された吸気管圧力と吸気管内圧力の検出値とに基づき所定期間先の予測時点における吸気管圧力を予測する吸気管圧力予測装置が既に提案されている(特開平2−42160号公報)。   Conventionally, the current intake pipe pressure is predicted based on the operating state of the internal combustion engine, and the intake pipe pressure at a prediction point ahead of a predetermined period is predicted based on the predicted intake pipe pressure and the detected value of the intake pipe pressure. An intake pipe pressure prediction device has been proposed (Japanese Patent Laid-Open No. 2-42160).

しかしながら、上記従来の予測装置においては、過渡時における吸気弁閉時の予測吸気管内圧力を内燃機関の運転状態とは無関係に吸気管内圧力の偏差ΔPに基づき一律に修正を施しているため、内燃機関の回転数の広い範囲にわたり十分な予測を行うことができなかった。   However, in the above-described conventional prediction device, the predicted intake pipe pressure when the intake valve is closed during the transition is uniformly corrected based on the deviation ΔP of the intake pipe pressure regardless of the operating state of the internal combustion engine. A sufficient prediction could not be made over a wide range of engine speeds.

特2886771号は、上記の問題点を解決するために、内燃機関の運転状態を考慮して予測を行い、低回転域や吸気管圧力が高い場合でも高精度な制御を行うことができる吸気管圧力予測装置を開示している。   In Japanese Patent No. 2886771, in order to solve the above-mentioned problems, an intake pipe that performs prediction in consideration of the operating state of the internal combustion engine and can perform high-precision control even in a low rotation range or when the intake pipe pressure is high. A pressure prediction device is disclosed.

上記の従来の方法においては、吸気管圧力値(以下、PBと呼称する)の差分(以下、ΔPBと呼称する)とスロットル弁開度値の差分(以下、ΔTHと呼称する)とに基づいてPBの予測値 (以下、HATPBと呼称する) を算出し、このHATPBに基づいて燃料噴射制御や燃料付着補正パラメータの検索などを行っていた。ここで、kを吸入行程(TDC)に同期した制御時刻としてΔTHおよびΔPBは、以下の式で表される。   In the above conventional method, based on the difference between the intake pipe pressure values (hereinafter referred to as PB) (hereinafter referred to as ΔPB) and the difference between the throttle valve opening values (hereinafter referred to as ΔTH). A predicted value of PB (hereinafter referred to as HATPB) is calculated, and fuel injection control and retrieval of a fuel adhesion correction parameter are performed based on the HATPB. Here, ΔTH and ΔPB are expressed by the following equations, where k is a control time synchronized with the suction stroke (TDC).

ΔTH(k)=TH(k)+TH(k−1) (1)
ΔPB(k)=PB(k)+PB(k−1) (2)
ΔTH (k) = TH (k) + TH (k−1) (1)
ΔPB (k) = PB (k) + PB (k−1) (2)

他方、最近では、PBセンサへの水の進入や進入した水の氷結による破壊防止などのため、PBセンサの検出部と吸気管との間の、吸気管側、またはPBセンサ内に、水の進入を防ぐためのラビリンス機構などが設けられている。このため、PBセンサ出力の実際の圧力値に対する遅れやむだ時間が大きくなってきた。   On the other hand, recently, in order to prevent water from entering the PB sensor and preventing damage due to freezing of the water that has entered the water, there is water on the intake pipe side or in the PB sensor between the detection part of the PB sensor and the intake pipe. A labyrinth mechanism is provided to prevent entry. For this reason, the delay and dead time of the PB sensor output with respect to the actual pressure value have increased.

そこで、このむだ時間を従来の予測アルゴリズムによって補償することが試みられたが、補償すべき遅れやむだ時間が大きくなったため、図2に示すように、HATPBが実際のPBに対してオーバーシュートを生じたり、不連続な挙動を示すようになった。オーバーシュートを生じる原因は、従来の予測アルゴリズムが、PBの予測値の精度不足を所定時間前のPBの予測値と現在のPBとの偏差をフィードバックすることによって補っていたことによる。また、不連続な挙動を示す原因は、従来の予測アルゴリズムが、ΔPBとΔTHのそれぞれに基づいて算出された予測値を所定の条件に応じて切り替えて使用していたことによる。このような挙動を示す従来の予測アルゴリズムによって、燃料制御が影響を受け、過渡運転時の変動が大きくなり、排ガスの有害成分の排出量が増加するという問題が生じてきた。
特開平2−42160号公報 特許第2886771号
Therefore, an attempt was made to compensate for this dead time by a conventional prediction algorithm. However, since delay and dead time to be compensated increased, HATPB overshoots actual PB as shown in FIG. It started to appear or showed discontinuous behavior. The cause of the overshoot is that the conventional prediction algorithm compensates for the insufficient accuracy of the PB prediction value by feeding back the deviation between the PB prediction value a predetermined time ago and the current PB. Moreover, the cause which shows discontinuous behavior is because the conventional prediction algorithm switched and used the predicted value calculated based on each of (DELTA) PB and (DELTA) TH according to predetermined conditions. The conventional prediction algorithm exhibiting such behavior has a problem in that fuel control is affected, fluctuation during transient operation becomes large, and emission of harmful components of exhaust gas increases.
JP-A-2-42160 Japanese Patent No. 2886771

したがって、HATPBのオーバーシュートや不連続挙動を生じずに、大きくなった遅れやむだ時間を補償することができる新たなPB予測アルゴリズムが求められていた。   Therefore, there has been a demand for a new PB prediction algorithm that can compensate for the increased delay and dead time without causing overshoot or discontinuous behavior of HATPB.

本発明の1態様による吸気管圧力予測方法は、吸気管圧力値の差分およびスロットル弁開度値の差分を求めるステップと、ファジィ推定アルゴリズムにより吸気管圧力の予測差分を求めるステップと、吸気管圧力値と吸気管圧力の予測差分の値を加算して吸気管圧力の予測値を求めるステップとを含む。ファジィ推定アルゴリズムは、吸気管圧力値の差分の大きさおよびスロットル弁開度値の差分の大きさに基づいて定めたファジィルールを含む。本方法においては、スロットル弁開度値とスロットル弁開度目標値を無駄時間要素と遅れ系を用いてモデル化し、当該モデルと当該目標値によって推定された値をスロットル弁開度値として用いる。   An intake pipe pressure prediction method according to an aspect of the present invention includes a step of obtaining a difference between intake pipe pressure values and a difference between throttle valve opening values, a step of obtaining a predicted difference in intake pipe pressure using a fuzzy estimation algorithm, and an intake pipe pressure Adding a value and a predicted difference value of the intake pipe pressure to obtain a predicted value of the intake pipe pressure. The fuzzy estimation algorithm includes a fuzzy rule determined based on the magnitude of the difference between the intake pipe pressure values and the magnitude of the difference between the throttle valve opening values. In this method, a throttle valve opening value and a throttle valve opening target value are modeled using a dead time element and a delay system, and a value estimated by the model and the target value is used as a throttle valve opening value.

本発明の別の態様による吸気管圧力予測装置は、吸気管圧力値の差分を求める差分器と、スロットル弁開度値の差分を求める差分器と、ファジィ推定器と、スロットル弁開度値とスロットル弁開度目標値を無駄時間要素と遅れ系を用いてモデル化したモデルとを含む。ファジィ推定器は、吸気管圧力値の差分およびスロットル弁開度値の差分を入力とし、ファジィ推定アルゴリズムにより吸気管圧力の予測差分を求め、当該予測差分を出力とする。ファジィ推定アルゴリズムは、吸気管圧力値の差分の大きさおよびスロットル弁開度値の差分の大きさに基づいて定めたファジィルールを含む。本装置において、スロットル弁開度値として、当該モデルと当該目標値によって推定された値を用いる。   An intake pipe pressure predicting apparatus according to another aspect of the present invention includes a differencer for obtaining a difference between intake pipe pressure values, a differencer for obtaining a difference between throttle valve opening values, a fuzzy estimator, a throttle valve opening value, The model which modeled the throttle valve opening target value using the dead time element and the delay system is included. The fuzzy estimator receives the difference between the intake pipe pressure values and the difference between the throttle valve opening values, obtains a predicted difference in the intake pipe pressure by a fuzzy estimation algorithm, and outputs the predicted difference as an output. The fuzzy estimation algorithm includes a fuzzy rule determined based on the magnitude of the difference between the intake pipe pressure values and the magnitude of the difference between the throttle valve opening values. In this device, the value estimated by the model and the target value is used as the throttle valve opening value.

このように、本発明の上記の2態様においては、PBセンサ出力とTH開度に基づいたファジィ推定アルゴリズムにより、PBセンサ出力の予測値を算出し、その予測値に基づいて内燃機関の燃料噴射料を決定する。特に、ΔPBの大きさとΔTHの大きさに基づいて定めたファジィルールを使用することにより、PBの変化に先行するTHの変化の情報を有効に取り込んだ制御が行えるようになる。これによって、PBセンサ出力が、実際の吸気管圧力(負圧)に対して大きなむだ時間や遅れを生じている場合にも、従来手法のように予測値が不連続的になることがなく、連続的な予測値を精度よく算出することができる。したがって、内燃機関の空燃比が不連続な挙動を示すことがなくなる。また、予測値の実際の吸気管圧力に対するオーバーシュートも従来手法に比べて飛躍的に小さくできる。   Thus, in the above two aspects of the present invention, the predicted value of the PB sensor output is calculated by the fuzzy estimation algorithm based on the PB sensor output and the TH opening, and the fuel injection of the internal combustion engine is performed based on the predicted value. Determine the fee. In particular, by using a fuzzy rule determined based on the magnitude of ΔPB and the magnitude of ΔTH, it is possible to perform control in which information on the change in TH preceding the change in PB is taken in effectively. As a result, even when the PB sensor output has a large dead time or delay with respect to the actual intake pipe pressure (negative pressure), the predicted value does not become discontinuous as in the conventional method. A continuous prediction value can be calculated with high accuracy. Therefore, the air-fuel ratio of the internal combustion engine does not show discontinuous behavior. Further, the overshoot of the predicted value with respect to the actual intake pipe pressure can be remarkably reduced as compared with the conventional method.

また、本発明の上記の2態様においては、スロットル弁開度値とスロットル弁開度目標値を無駄時間要素と遅れ系を用いてモデル化し、当該モデルと当該目標値によって推定された値をスロットル弁開度値として用いる。すなわち、時点kにおけるスロットル弁開度推定値およびスロットル弁開度目標値をTHHAT(k)およびTHCMD(k)、むだ時間相当値をddly、定数をKdlyとしたときに、式
THHAT(k)=Kdly× THHAT(k)+ (1−Kdly)× THCMD(k-ddly)
によって求めたスロットル弁開度推定値THHAT(k)を使用する。
In the two aspects of the present invention, the throttle valve opening value and the throttle valve opening target value are modeled using a dead time element and a delay system, and the value estimated by the model and the target value is throttled. Used as a valve opening value. That is, when the throttle valve opening estimated value and the throttle valve opening target value at time k are THHAT (k) and THCMD (k), the dead time equivalent value is ddly, and the constant is Kdly,
THHAT (k) = Kdly x THHAT (k) + (1-Kdly) x THCMD (k-ddly)
The throttle valve opening estimated value THHAT (k) obtained by the above is used.

したがって、THの読み込みまでに無効時間が生じる電子制御スロットルを使用する場合にも、精度の高い予測値を算出することができる。   Therefore, a highly accurate predicted value can be calculated even when using an electronically controlled throttle that has an invalid time before TH is read.

本発明の1実施形態によれば、吸気管圧力値の差分の大きさを正、零、負に区分し、スロットル弁開度値の差分の大きさを正、零、負に区分して、2種類の区分によって定まる領域に対してそれぞれファジィルールを定める。   According to one embodiment of the present invention, the magnitude of the difference in the intake pipe pressure value is divided into positive, zero, and negative, and the magnitude of the difference in the throttle valve opening value is divided into positive, zero, and negative, A fuzzy rule is defined for each area determined by two types of classification.

したがって、ファジィ推定演算を複雑にすることなく、簡単かつ効率的に行うことができる。   Therefore, the fuzzy estimation calculation can be performed easily and efficiently without making it complicated.

本発明の別の実施形態によれば、さらに吸気管圧力値の2回差分を求め、吸気管圧力値の差分の大きさ、スロットル弁開度値の差分の大きさおよび吸気管圧力値の2回差分の大きさに基づいてファジィルールを定める。   According to another embodiment of the present invention, the two-time difference between the intake pipe pressure values is further obtained, and the difference between the intake pipe pressure values, the difference between the throttle valve opening values, and the intake pipe pressure value of 2 are calculated. A fuzzy rule is determined based on the magnitude of the time difference.

したがって、ΔPBの差分(以下、ΔΔPBと呼称する)に含まれる、ΔPBの将来の挙動についての情報を有効に取り込んだ制御が行えるようになる。   Therefore, it becomes possible to perform control that effectively incorporates information about the future behavior of ΔPB included in the difference of ΔPB (hereinafter referred to as ΔΔPB).

本発明の別の実施形態によれば、吸気管圧力値の2回差分の大きさを正、零、負に区分し、吸気管圧力値の差分の大きさ、スロットル弁開度値の差分の大きさおよび吸気管圧力値の2回差分の大きさの合計3種類の区分によって定まる領域に対してそれぞれファジィルールを定める。   According to another embodiment of the present invention, the magnitude of the two-time difference in the intake pipe pressure value is divided into positive, zero, and negative, the magnitude of the difference in the intake pipe pressure value, and the difference in the throttle valve opening value. A fuzzy rule is defined for each region determined by a total of three types of size and the size of the difference between the two times of the intake pipe pressure value.

したがって、ファジィ推定演算を複雑にすることなく、簡単かつ効率的に行うことができる。   Therefore, the fuzzy estimation calculation can be performed easily and efficiently without making it complicated.

本発明の別の実施形態によれば、ファジィ推定アルゴリズムの後件部メンバーシップ関数は、シングルトンの棒状関数である。   According to another embodiment of the invention, the consequent membership function of the fuzzy estimation algorithm is a singleton bar function.

したがって、ファジィ推定演算を複雑にすることなく、簡単かつ効率的に行うことができる。特に、極低温、高温、多湿、振動などの車両の使用条件に耐えうるプロセッサの演算能力では、ファジィ推定アルゴリズムのミニマックス重心法の演算を実行することは困難であるが、後件部メンバーシップ関数をシングルトンの棒状関数とすることにより、その演算が可能となり、十分な精度を有する予測を算出することができる。   Therefore, the fuzzy estimation calculation can be performed easily and efficiently without making it complicated. In particular, it is difficult to execute the minimax centroid method of the fuzzy estimation algorithm with the computing power of the processor that can withstand the vehicle usage conditions such as cryogenic temperature, high temperature, high humidity, and vibration, but the consequent membership By making the function a singleton bar-like function, the calculation can be performed, and a prediction having sufficient accuracy can be calculated.

本発明の別の実施形態によれば、入力のフィルタリングが行われる。   According to another embodiment of the invention, input filtering is performed.

したがって、ファジィ推定アルゴリズムへの入力データをフィルタリングされたデータとすることによって、予測値にノイズが混入した場合に、そのノイズの影響によって予測値が振動的になるのを防止できる。   Therefore, by setting the input data to the fuzzy estimation algorithm as filtered data, when noise is mixed in the predicted value, it is possible to prevent the predicted value from oscillating due to the influence of the noise.

本発明の別の実施形態によれば、フィルタリングが適応フィルタによって行われる。   According to another embodiment of the invention, the filtering is performed by an adaptive filter.

したがって、予測演算のためのデータとして使用できるように位相遅れを最小限に保ちつつ、十分にノイズを除去することができる。   Therefore, it is possible to sufficiently remove noise while keeping the phase delay to a minimum so that it can be used as data for prediction calculation.

本発明の別の態様による吸気管圧力予測方法は、現時点の吸気管圧力値の差分、現時点の吸気管圧力値の2回差分および現時点のスロットル弁開度値の差分を求めるステップと、ファジィ推定アルゴリズムにより吸気管圧力の予測差分を求めるステップと、現時点の吸気管圧力値と前記吸気管圧力の予測差分の値を加算して吸気管圧力の予測値を求めるステップとを含む。ファジィ推定アルゴリズムは、前記現時点の吸気管圧力値の差分の大きさ、前記現時点のスロットル弁開度値の差分の大きさおよび前記現時点の吸気管圧力値の2回差分の大きさに基づいて定めたファジィルールを含む。   An intake pipe pressure prediction method according to another aspect of the present invention includes a step of obtaining a difference between a current intake pipe pressure value, a two-time difference between current intake pipe pressure values, and a difference between current throttle valve opening values, and fuzzy estimation. A step of obtaining a prediction difference of the intake pipe pressure by an algorithm, and a step of obtaining a predicted value of the intake pipe pressure by adding a current intake pipe pressure value and a value of the prediction difference of the intake pipe pressure. The fuzzy estimation algorithm is determined based on the magnitude of the difference between the current intake pipe pressure values, the magnitude of the difference between the current throttle valve opening values, and the magnitude of the two-time difference between the current intake pipe pressure values. Includes fuzzy rules.

本発明の別の態様による吸気管圧力予測装置は、現時点の吸気管圧力値の差分を求める差分器と、現時点の吸気管圧力値の2回差分を求める差分器と、現時点のスロットル弁開度値の差分を求める差分器と、前記現時点の吸気管圧力値の差分、前記現時点のスロットル弁開度値の差分および前記現時点の吸気管圧力値の2回差分を入力として、ファジィ推定アルゴリズムにより吸気管圧力の予測差分を求め、当該予測差分を出力とするファジィ推定器と、を含む。ファジィ推定アルゴリズムは、前記現時点の吸気管圧力値の差分の大きさ、前記現時点のスロットル弁開度値の差分の大きさおよび前記現時点の吸気管圧力値の2回差分の大きさに基づいて定めたファジィルールを含む。   An intake pipe pressure predicting apparatus according to another aspect of the present invention includes a differencer for obtaining a difference between current intake pipe pressure values, a differencer for obtaining a difference between current intake pipe pressure values twice, and a current throttle valve opening degree. A differencer for obtaining a difference between the values, a difference between the current intake pipe pressure value, a difference between the current throttle valve opening value and a second difference between the current intake pipe pressure values and an intake air by a fuzzy estimation algorithm A fuzzy estimator that calculates a prediction difference of the tube pressure and outputs the prediction difference. The fuzzy estimation algorithm is determined based on the magnitude of the difference between the current intake pipe pressure values, the magnitude of the difference between the current throttle valve opening values, and the magnitude of the two-time difference between the current intake pipe pressure values. Includes fuzzy rules.

上記の2態様によれば、吸気管圧力値の差分の大きさ、スロットル弁開度値の差分の大きさおよび吸気管圧力値の2回差分の大きさに基づいてファジィルールを定めるので、PBの変化(ΔPB)に先行するTHの変化(ΔTH)の情報を有効に取り込むとともに、ΔΔPBに含まれる、ΔPBの将来の挙動についての情報を有効に取り込んだ制御が行えるようになる。   According to the above two aspects, the fuzzy rule is determined based on the magnitude of the difference between the intake pipe pressure values, the magnitude of the difference between the throttle valve opening values, and the magnitude of the twice difference between the intake pipe pressure values. It is possible to effectively take in the information on the change (ΔTH) of the TH preceding the change (ΔPB) in the control, and effectively take in the information on the future behavior of ΔPB included in the ΔΔPB.

本発明の別の態様による吸気管圧力予測方法は、現時点の吸気管圧力値の差分および現時点のスロットル弁開度値の差分を求めるステップと、前記吸気管圧力値の差分および前記スロットル弁開度値の差分をフィルタリングするステップと、ファジィ推定アルゴリズムにより吸気管圧力の予測差分を求めるステップと、現時点の吸気管圧力値と前記吸気管圧力の予測差分の値を加算して吸気管圧力の予測値を求めるステップとを含む。ファジィ推定アルゴリズムは、フィルタリングされた前記現時点の吸気管圧力値の差分の大きさおよびフィルタリングされた前記現時点のスロットル弁開度値の差分の大きさに基づいて定めたファジィルールを含む。   An intake pipe pressure prediction method according to another aspect of the present invention includes a step of obtaining a difference between a current intake pipe pressure value and a current throttle valve opening value, and a difference between the intake pipe pressure value and the throttle valve opening degree. A step of filtering the difference between the values, a step of obtaining a prediction difference of the intake pipe pressure by a fuzzy estimation algorithm, and a predicted value of the intake pipe pressure by adding a value of the current intake pipe pressure value and the prediction difference of the intake pipe pressure Determining. The fuzzy estimation algorithm includes a fuzzy rule determined based on the magnitude of the difference between the filtered current intake pipe pressure values and the magnitude of the filtered difference between the current throttle valve opening values.

本発明の別の態様による吸気管圧力予測装置は、現時点の吸気管圧力値の差分を求める第1の差分器と、現時点のスロットル弁開度値の差分を求める第2の差分器と、前記現時点の吸気管圧力値の差分および前記現時点のスロットル弁開度値の差分を入力として、ファジィ推定アルゴリズムにより吸気管圧力の予測差分を求め、当該予測差分を出力とするファジィ推定器と、を含む。ファジィ推定アルゴリズムは、前記現時点の吸気管圧力値の差分の大きさおよび前記現時点のスロットル弁開度値の差分の大きさに基づいて定めたファジィルールを含む。さらに、本装置は、前記第1の差分器と前記ファジィ推定器との間に配置された第1のフィルタと、前記第2の差分器と前記ファジィ推定器との間に配置された第2のフィルタと、を含む。   An intake pipe pressure predicting device according to another aspect of the present invention includes a first differentiator that obtains a difference between current intake pipe pressure values, a second differentiator that obtains a difference between current throttle valve opening values, A fuzzy estimator that calculates a difference in the intake pipe pressure by a fuzzy estimation algorithm using the difference in the current intake pipe pressure value and the difference in the current throttle valve opening value as an input, and outputs the prediction difference . The fuzzy estimation algorithm includes a fuzzy rule determined based on the magnitude of the difference in the current intake pipe pressure value and the magnitude of the difference in the current throttle valve opening value. Furthermore, the apparatus includes a first filter disposed between the first differencer and the fuzzy estimator, and a second filter disposed between the second differencer and the fuzzy estimator. And a filter.

上記の2態様によれば、PBの変化に先行するTHの変化の情報を有効に取り込んだ制御が行えるようになる。また、ファジィ推定アルゴリズムへの入力データをフィルタリングされたデータとすることによって、予測値にノイズが混入した場合に、そのノイズの影響によって予測値が振動的になるのを防止できる。   According to the above two aspects, it is possible to perform control that effectively incorporates information on the change in TH preceding the change in PB. Further, by making the input data to the fuzzy estimation algorithm into filtered data, when noise is mixed in the predicted value, it is possible to prevent the predicted value from becoming oscillating due to the influence of the noise.

本発明は、ΔPBの大きさおよびΔTHの大きさに基づいて定めたファジィルールを含むファジィ推定アルゴリズムを使用してPBのファジィ推定による予測変化量(以下、ΔFZPBと呼称する)を求める。つぎに以下の式により、ファジィ推定による予測値(以下、FZPBと呼称する)を算出する。   The present invention obtains a predicted change amount (hereinafter referred to as ΔFZPB) by fuzzy estimation of PB using a fuzzy estimation algorithm including a fuzzy rule determined based on the magnitude of ΔPB and the magnitude of ΔTH. Next, a predicted value by fuzzy estimation (hereinafter referred to as FZPB) is calculated by the following equation.

FZPB(k)=PB(k)+ΔFZPB(k) (3)
すなわち、現在のPBサンプル値PB(k)にΔFZPB(k)を加算することによりFZPB(k)を求める。ここで、kは、吸入行程(TDC)に同期した制御時刻を表す。図3は、PB、TH、ΔPB、ΔTHおよびΔFZPBの間の関係を示す。
FZPB (k) = PB (k) + ΔFZPB (k) (3)
That is, FZPB (k) is obtained by adding ΔFZPB (k) to the current PB sample value PB (k). Here, k represents a control time synchronized with the suction stroke (TDC). FIG. 3 shows the relationship between PB, TH, ΔPB, ΔTH and ΔFZPB.

本実施形態で使用するファジィルールを図6の表に示す。ΔPBの大きさを正、零、負に区分し、ΔTHの大きさを正、零、負に区分して、2種類の区分によって定まる9個の領域に対してそれぞれファジィルールを定めている。ここで留意すべき点は、THの変化は、PBの変化に先行しており、将来のPBの挙動の情報を含む点である。したがって、ΔPBの大きさとΔTHの大きさに基づいて定めたファジィルールを使用することにより、PBの変化に先行するTHの変化の情報を有効に取り込んだ制御が行えるようになる。   The fuzzy rules used in this embodiment are shown in the table of FIG. The magnitude of ΔPB is divided into positive, zero, and negative, and the magnitude of ΔTH is divided into positive, zero, and negative, and fuzzy rules are defined for each of nine regions determined by two types of classification. It should be noted here that the change in TH precedes the change in PB and includes information on the future behavior of PB. Therefore, by using a fuzzy rule determined based on the magnitude of ΔPB and the magnitude of ΔTH, it is possible to perform control that effectively incorporates information on the change in TH preceding the change in PB.

本発明のファジィ推定アルゴリズムで使用する前件部メンバーシップ関数と後件部メンバーシップ関数とをそれぞれ図4および5に示す。ΔPBおよびΔTHに対する前件部メンバーシップ関数は、正(P)および負(N)に対しては台形状とし、零(Z)に対しては三角形状とした。後件部メンバーシップ関数としては、ファジィ推定演算の簡素化を実現するようにシングルトンの棒状関数を用いた。   The antecedent part membership function and the consequent part membership function used in the fuzzy estimation algorithm of the present invention are shown in FIGS. 4 and 5, respectively. The antecedent part membership functions for ΔPB and ΔTH were trapezoidal for positive (P) and negative (N) and triangular for zero (Z). As the consequent part membership function, a singleton bar function was used to simplify the fuzzy estimation calculation.

つぎに各領域におけるルール(図6におけるルール1から9)のそれぞれについて図を使用して説明する。   Next, rules in each area (rules 1 to 9 in FIG. 6) will be described with reference to the drawings.

図7は、ルール1が適用される、ΔPBおよびΔTHがともに負の状況を示す。ΔPBおよびΔTHがともに負であるので、後件部メンバーシップ関数のΔFZPBも負とする。   FIG. 7 shows a situation in which ΔPB and ΔTH are both negative, where rule 1 is applied. Since ΔPB and ΔTH are both negative, ΔFZPB of the consequent part membership function is also negative.

図8は、ルール2が適用される、ΔPBが負およびΔTHが零の状況を示す。先行するΔTHが零であるがΔPBが負であるので、後件部メンバーシップ関数のΔFZPBは負とする。   FIG. 8 illustrates a situation where rule 2 is applied and ΔPB is negative and ΔTH is zero. Since the preceding ΔTH is zero but ΔPB is negative, ΔFZPB of the consequent part membership function is negative.

図9は、ルール3が適用される、ΔPBが負およびΔTHが正の状況を示す。この状況は、エンジンブレーキ時、TH増加によるスロットル通過空気量の増加よりもエンジン回転数上昇が速い場合に対応する。後件部メンバーシップ関数のΔFZPBは零とする。   FIG. 9 shows a situation where rule 3 is applied and ΔPB is negative and ΔTH is positive. This situation corresponds to a case where the engine speed rises faster than the increase in the amount of air passing through the throttle due to the increase in TH during engine braking. The consequent part membership function ΔFZPB is assumed to be zero.

図10は、ルール4が適用される、ΔPBが零およびΔTHが負の状況を示す。先行するΔTHが負であるので、後件部メンバーシップ関数のΔFZPBは負とする。   FIG. 10 shows a situation where rule 4 is applied and ΔPB is zero and ΔTH is negative. Since the preceding ΔTH is negative, ΔFZPB of the consequent part membership function is negative.

図11は、ルール5が適用される、ΔPBおよびΔTHがともに零の状況を示す。ΔPBおよびΔTHがともに零であるので、後件部メンバーシップ関数のΔFZPBも零とする。   FIG. 11 shows a situation in which ΔPB and ΔTH are both zero, where rule 5 is applied. Since both ΔPB and ΔTH are zero, ΔFZPB of the consequent part membership function is also zero.

図12は、ルール6が適用される、ΔPBが零およびΔTHが正の状況を示す。先行するΔTHが正であるので、後件部メンバーシップ関数のΔFZPBは正(P)とする。 FIG. 12 shows a situation where rule 6 is applied and ΔPB is zero and ΔTH is positive. Since the preceding ΔTH is positive, ΔFZPB of the consequent part membership function is positive (P 6 ).

図13は、ルール7が適用される、ΔPBが正およびΔTHが負の状況を示す。この状況は、エンジン回転数が低下する際、TH減少による低下よりも、外力によって大幅にエンジン回転数が低下した場合に対応する。後件部メンバーシップ関数のΔFZPBは零とする。   FIG. 13 shows a situation where rule 7 is applied and ΔPB is positive and ΔTH is negative. This situation corresponds to a case where the engine speed is greatly reduced by an external force when the engine speed is reduced, rather than by a decrease in TH. The consequent part membership function ΔFZPB is assumed to be zero.

図14は、ルール8が適用される、ΔPBが正およびΔTHが零の状況を示す。後件部メンバーシップ関数のΔFZPBは正(P)とするが、先行するΔTHが零であるのでルール6のPよりも小さな値とする。 FIG. 14 shows a situation where rule 8 is applied and ΔPB is positive and ΔTH is zero. The consequent part membership function ΔFZPB is positive (P 8 ), but since the preceding ΔTH is zero, the value is smaller than P 6 of rule 6.

図15は、ルール9が適用される、ΔPBおよびΔTHがともに正の状況を示す。後件部メンバーシップ関数のΔFZPBは正(P10)とするが、ΔPBが既に正であるので、ΔPBが零であるルール6のPよりも小さな値とする。 FIG. 15 shows a situation in which ΔPB and ΔTH are both positive, to which rule 9 is applied. Although ΔFZPB of the consequent part membership function is positive (P 10 ), since ΔPB is already positive, it is set to a value smaller than P 6 of rule 6 in which ΔPB is zero.

最終的なΔFZPBは、上述のメンバーシップ関数とファジィルールを用いてミニマックス重心法により算出される。以下にルール6を例として具体的に説明する。   The final ΔFZPB is calculated by the minimax centroid method using the above-described membership function and fuzzy rules. The rule 6 will be specifically described below as an example.

まず、ミニマックス選択について説明する。現在のΔPB、ΔTHのサンプル値をΔPB(k)、ΔTH(k)とする。これらの値の、ルール6に対する適合度を求める。図6および図12に示すように、ルール6は、ΔPBが零でΔTHが正である領域に適用される。そこで、ΔPBの、前件部メンバーシップ関数(零)に対する適合度は、図16(a)に示すようにmΔPB(6)となる。また、ΔTHの、前件部メンバーシップ関数(正)に対する適合度は、図16(b)に示すようにmΔTH(6)となる。ここで、ミニマックス選択では、前件部の適合度は、同一ルールにおける適合度のうち最も小さいものをそのルールの適合度m(i)とするため、ルール6に対する適合度は以下の通り、mΔPB(6)となる。   First, minimax selection will be described. Let the current sample values of ΔPB and ΔTH be ΔPB (k) and ΔTH (k). The degree of conformity of these values to rule 6 is obtained. As shown in FIGS. 6 and 12, rule 6 is applied to a region where ΔPB is zero and ΔTH is positive. Therefore, the degree of fitness of ΔPB with respect to the antecedent part membership function (zero) is mΔPB (6) as shown in FIG. Further, the degree of fitness of ΔTH with respect to the antecedent part membership function (positive) is mΔTH (6) as shown in FIG. Here, in the minimax selection, since the degree of conformity of the antecedent part is the smallest degree of conformance in the same rule as the conformity m (i) of the rule, mΔPB (6).

mΔPB(6)<mΔTH(6) (4)
したがって
m(6)= mΔPB(6) (5)
となる。
mΔPB (6) <mΔTH (6) (4)
Therefore, m (6) = mΔPB (6) (5)
It becomes.

さらに、ルール6の後件部のメンバーシップ関数は、図17に示すように、位置pP6、基準重み(棒状関数の長さ)wP6であるので、ΔFZPBの推定におけるルール6の重み、すなわち、位置pPにおける重みw(6)は、
w(6)=m(6)×wP6 (6)
となる。
Further, as shown in FIG. 17, the membership function of the consequent part of rule 6 is a position pP6 and a reference weight (the length of a bar function) wP6. Therefore, the weight of rule 6 in estimating ΔFZPB, The weight w (6) at pP is
w (6) = m (6) × wP6 (6)
It becomes.

つぎに、各ルールについて、同様に重みw(i)を求める。このように求めた各ルールの重みが全て推定に使用されるため、これはマックス選択となる。   Next, the weight w (i) is similarly obtained for each rule. Since all the weights of the respective rules thus obtained are used for estimation, this is a maximum selection.

つぎに、求められた各ルールに対する重みに基づいてΔFZPBの推定、すなわちファジィルールの非ファジィ化を以下の式に示す重心法に基づいて行う。

Figure 0004291863
Next, ΔFZPB is estimated based on the calculated weights for each rule, that is, fuzzy rules are defuzzified based on the centroid method shown in the following equation.
Figure 0004291863

したがって、式3は具体的に以下の式によって表され、ΔFZPB(k)が算出される。

Figure 0004291863
Therefore, Formula 3 is specifically expressed by the following formula, and ΔFZPB (k) is calculated.
Figure 0004291863

図18は、上述のファジィ推定アルゴリズムを使用したFZPBの推定結果である。図からわかるようにFZPBは、PBに対して多少のオーバーシュートを伴うものの、おおよそ推定すべき予測値を再現しており、図2に示す、従来の制御アルゴリズムによる推定に比較して大幅に推定精度が高くなっている。   FIG. 18 shows FZPB estimation results using the above-described fuzzy estimation algorithm. As can be seen from the figure, FZPB reproduces the estimated value that should be estimated roughly, with some overshooting compared to PB, and is greatly estimated compared to the estimation by the conventional control algorithm shown in FIG. The accuracy is high.

つぎに、本発明の別の実施形態について説明する。上述の実施形態によるファジィ予測では、図18に示すように多少のオーバーシュートの発生が見られた。このようなオーバーシュートは、特に急加速時やスナップ時に生じやすい。そこで、このようなオーバーシュートを消去するために、ΔPBの変化量、すなわちPBの2回差分(以下、ΔΔPBと呼称する)を、kを吸入行程(TDC)に同期した制御時刻として以下の式で定義した。   Next, another embodiment of the present invention will be described. In the fuzzy prediction according to the above-described embodiment, some overshoot was observed as shown in FIG. Such an overshoot is particularly likely to occur during sudden acceleration or snapping. Therefore, in order to eliminate such overshoot, the amount of change of ΔPB, that is, the two-time difference of PB (hereinafter referred to as ΔΔPB) is set as the control time in which k is synchronized with the suction stroke (TDC). Defined in

ΔΔPB(k)=ΔPB(k)−ΔPB(k−1) (9)   ΔΔPB (k) = ΔPB (k) −ΔPB (k−1) (9)

そして、ΔPBの大きさおよびΔTHの大きさに加えてΔΔPBの大きさに基づいてファジィルールを定めるようにした。すなわち、ΔPBの大きさ、ΔTHの大きさおよびΔΔPBの大きさの合計3種類の区分によって定まる領域に対してそれぞれファジィルールを定めた。本実施形態によるファジィルールを図20に示す。図6のファジィルールが2次元であるのに対し、図20のファジィルールは3次元である。ここで、留意すべき点は、ΔΔPBは、ΔPBの将来の挙動についての情報を含む点である。したがって、ΔΔPBに基づいたファジィルールを使用することによりこの情報を有効に取り込んだ制御が行えるようになる。また、本実施形態による前件部メンバーシップ関数を図19に示す。   A fuzzy rule is determined based on the magnitude of ΔΔPB in addition to the magnitude of ΔPB and the magnitude of ΔTH. In other words, fuzzy rules were respectively defined for areas determined by three types of divisions, that is, the size of ΔPB, the size of ΔTH, and the size of ΔΔPB. A fuzzy rule according to this embodiment is shown in FIG. While the fuzzy rule in FIG. 6 is two-dimensional, the fuzzy rule in FIG. 20 is three-dimensional. Here, it should be noted that ΔΔPB includes information about the future behavior of ΔPB. Therefore, by using a fuzzy rule based on ΔΔPB, it is possible to perform control that effectively incorporates this information. FIG. 19 shows an antecedent part membership function according to this embodiment.

図20のファジィルールは、図6のファジィルールにおけるルール8とルール9のみが、ΔΔPBが正と負の場合に変化している。他のルールは、ΔΔPBに関係なく変化していない。   In the fuzzy rule of FIG. 20, only the rules 8 and 9 in the fuzzy rule of FIG. 6 are changed when ΔΔPB is positive and negative. The other rules are unchanged regardless of ΔΔPB.

つぎに新たなルール(図12におけるルール10および11)のそれぞれについて図を使用して説明する。   Next, each of the new rules (rules 10 and 11 in FIG. 12) will be described with reference to the drawings.

図22は、ルール10が適用される、ΔPB、ΔTHおよびΔΔPBがともに正の状況を示す。ΔPB、ΔTHおよびΔΔPBがともに正であるので、後件部メンバーシップ関数のΔFZPBは最も大きな正の値(P10)とする。 FIG. 22 shows a situation where ΔPB, ΔTH, and ΔΔPB are all positive where rule 10 is applied. Since ΔPB, ΔTH, and ΔΔPB are all positive, ΔFZPB of the consequent part membership function is set to the largest positive value (P 10 ).

図23は、ルール11が適用される、ΔPBが正、ΔTHが零であり、ΔΔPBが負の状況を示す。後件部メンバーシップ関数のΔFZPBは零とする。   FIG. 23 shows a situation where rule 11 is applied, ΔPB is positive, ΔTH is zero, and ΔΔPB is negative. The consequent part membership function ΔFZPB is assumed to be zero.

図24は、本実施形態においてルール6が適用される、ΔPBが零でΔTHが正の状況を示す。上述のように、ΔPB、ΔTHおよびΔΔPBがともに正の状況には、ルール10が適用されるので、ルール6が適用される範囲は、ΔΔPB未使用の場合と比較して縮小される。   FIG. 24 shows a situation where rule 6 is applied in the present embodiment and ΔPB is zero and ΔTH is positive. As described above, when ΔPB, ΔTH, and ΔΔPB are all positive, rule 10 is applied. Therefore, the range to which rule 6 is applied is reduced as compared with the case where ΔΔPB is not used.

図25は、本実施形態においてルール9が適用される、ΔPBおよびΔTHがともに正の状況を示す。上述のように、ΔPB、ΔTHおよびΔΔPBがともに正の状況には、ルール10が適用されるので、ルール9が適用される範囲は、ΔΔPB未使用の場合と比較して縮小される。   FIG. 25 shows a situation in which ΔPB and ΔTH are both positive, to which rule 9 is applied in the present embodiment. As described above, when ΔPB, ΔTH, and ΔΔPB are all positive, rule 10 is applied. Therefore, the range to which rule 9 is applied is reduced as compared with the case where ΔΔPB is not used.

図21に本実施形態の後件部メンバーシップ関数を示す。上述のように、P10は、最も大きな値とされている。図5の後件部メンバーシップ関数と比較すると、ルール10が新たに設けられ、ルール6およびルール9の範囲が縮小したので、pP6(Pの位置)およびpP9(Pの位置)が零に近づいている。 FIG. 21 shows the consequent part membership function of this embodiment. As described above, P 10 is the largest value. Compared with the consequent membership functions of FIG. 5, the rule 10 is newly provided, the range of the rule 6 and rule 9 is reduced, PP6 (position of P 6) and PP9 (position of P 9) is zero Is approaching.

図26は、本実施形態のファジィ推定アルゴリズムを使用したFZPBの推定結果を示す。図からわかるように、図18に見られるオーバーシュートが解消している。   FIG. 26 shows the FZPB estimation result using the fuzzy estimation algorithm of this embodiment. As can be seen, the overshoot seen in FIG. 18 is eliminated.

なお、本実施形態においては、補償すべきむだ時間が少ない場合は、零に近くなったPおよびPは、データ設定の簡素化のために零とすることもできる。また、図20のΔΔPBが正の場合のルール8は、現実には現れにくい状況であるので、このルールを消去したり、対応する後件部を零へ変更して予測への寄与度の無効化を行うこともある。 In the present embodiment, when the dead time to be compensated is small, P 6 and P 9 that are close to zero can be set to zero in order to simplify data setting. Further, rule 8 in the case where ΔΔPB is positive in FIG. 20 is a situation that is unlikely to appear in reality, so this rule is deleted or the corresponding consequent part is changed to zero to invalidate the contribution to the prediction. May also be performed.

つぎに、入力のフィルタリングを行う実施形態について説明する。本発明は、ΔPB、ΔTHおよびΔΔPBを入力とするファジィ推定アルゴリズムを使用する。したがって、THやPBにノイズが混入すると、これらの差分値や2回差分値が振動的またはスパイク的になり、ファジィ推定によって推定されたFZPBも振動的またはスパイク的になることがあった。   Next, an embodiment for performing input filtering will be described. The present invention uses a fuzzy estimation algorithm with ΔPB, ΔTH and ΔΔPB as inputs. Therefore, when noise is mixed into TH or PB, the difference value or the two-time difference value becomes vibrational or spike-like, and the FZPB estimated by fuzzy estimation sometimes becomes vibrational or spike-like.

したがって、本実施形態では、図27に示すように、ファジィ推定器への入力をΔPB、ΔTH、ΔΔPBとせずに、図28に示すようにこれらの値をフィルタリングした後の値を入力するようにした。図27において、2711乃至2713は、差分器を、2701はファジィ推定器を示す。図28において、2811乃至2813は差分器を、2821乃至2823はフィルタを、2801はファジィ推定器を示す   Therefore, in this embodiment, as shown in FIG. 27, the input to the fuzzy estimator is not set to ΔPB, ΔTH, and ΔΔPB, and values after filtering these values as shown in FIG. 28 are input. did. In FIG. 27, reference numerals 2711 to 2713 denote differentiators, and 2701 denotes a fuzzy estimator. In FIG. 28, 2811 to 2813 indicate differentiators, 2821 to 2823 indicate filters, and 2801 indicates a fuzzy estimator.

ΔPB、ΔTH、ΔΔPBのフィルタとしては、以下の式で表す適応フィルタを使用した。これによって図29に示すようにFZPBがノイズの影響により振動的スパイク的になることが防止される。   As the filters for ΔPB, ΔTH, and ΔΔPB, adaptive filters represented by the following equations were used. As a result, as shown in FIG. 29, FZPB is prevented from becoming a vibrational spike due to the influence of noise.

X_f (k) = X_f (k-1) + KP(k)・ide(k) (10)
KP(k) = P(k-1)/(1+P(k-1)) (11)
ide(k) = X_f (k-1) - X(k) (12)
P(k+1) = (1/λ1)[1-λ2・P(k)/( λ1 + λ2 ・P(k))] (13)
ここで、X_fは、ΔPB、ΔTH、ΔΔPBの適応フィルタ値を表し、Xは、ΔPB、ΔTH、ΔΔPBのサンプル値を表す。λ1,λ2は重みパラメータを表す。
X_f (k) = X_f (k-1) + KP (k) ・ ide (k) (10)
KP (k) = P (k-1) / (1 + P (k-1)) (11)
ide (k) = X_f (k-1)-X (k) (12)
P (k + 1) = (1 / λ 1 ) [1-λ 2 · P (k) / (λ 1 + λ 2 · P (k))] (13)
Here, X_f represents adaptive filter values of ΔPB, ΔTH, and ΔΔPB, and X represents sample values of ΔPB, ΔTH, and ΔΔPB. λ 1 and λ 2 represent weight parameters.

つぎに、THとして、スロットル弁開度目標値からむだ時間を考慮して求めた値を使用する実施形態について説明する。最近、燃費向上のためトランスミッションとの協調制御や操縦安定制御の要求から電子制御スロットルの採用が多くなっている。この電子制御スロットルは、電子制御ユニットとは、別のドライバによって駆動されることが多く、さらに電子制御ユニットとは車載ネットワーク(CAN等)によって接続されている。   Next, an embodiment in which a value obtained in consideration of the dead time from the throttle valve opening target value is used as TH will be described. In recent years, electronic throttles have been increasingly used to improve fuel efficiency because of demands for cooperative control with the transmission and steering stability control. The electronic control throttle is often driven by a driver different from the electronic control unit, and is further connected to the electronic control unit by an in-vehicle network (CAN or the like).

したがって、電子制御ユニットによって算出されるスロットル弁開度指令値THCMDと電子制御スロットルの実弁開度THACTとの間には、電子制御ユニットとドライバの間の通信周期(10msec)に起因するむだ時間が生じる。すなわち、THCMDがTHに反映されるまでに無効時間が生じる。また、THACTが、CAN等を介して電子制御ユニットによって観測されるTHに反映されるまでに無効時間が生じる。電子制御スロットル使用時の状況を図30に示す。図30において、上述のように、THCMDは、スロットル弁開度指令値を、THACTは、実弁開度を、THは弁開度値(観測値)を示す。   Therefore, the time delay caused by the communication cycle (10 msec) between the electronic control unit and the driver is between the throttle valve opening command value THCMD calculated by the electronic control unit and the actual valve opening THACT of the electronic control throttle. Occurs. That is, an invalid time occurs before THCMD is reflected in TH. In addition, an invalid time occurs until THACT is reflected in TH observed by the electronic control unit via CAN or the like. FIG. 30 shows the situation when the electronic control throttle is used. In FIG. 30, as described above, THCMD indicates the throttle valve opening command value, THACT indicates the actual valve opening, and TH indicates the valve opening value (observed value).

このため、PBの変化よりも、電子制御ユニットによって観測されるTHの変化の方が遅れてしまうため、上述のファジィ推定アルゴリズムのように、ΔTHを使用してPBを予測することはできない。   For this reason, since the change in TH observed by the electronic control unit is delayed from the change in PB, it is impossible to predict PB using ΔTH as in the above-described fuzzy estimation algorithm.

したがって、本実施形態では、以下の式に示すように通信によるむだ時間と電子制御スロットルの応答遅れを考慮して、THCMDからTHHATを推定し、このTHHATの差分ΔTHHATを、ΔTHの代わりに使用した。   Therefore, in the present embodiment, THHAT is estimated from THCMD in consideration of the dead time due to communication and the response delay of the electronic throttle as shown in the following equation, and this THHAT difference ΔTHHAT is used instead of ΔTH. .

THHAT(k)=Kdly× THHAT(k)+ (1−Kdly)× THCMD(k-ddly) (14)
ΔTHHAT(k)= THHAT(k)− THHAT(k-1) (15)
ここで、ddlyは、むだ時間相当値、Kdlyは、定数を示す。
THHAT (k) = Kdly x THHAT (k) + (1-Kdly) x THCMD (k-ddly) (14)
ΔTHHAT (k) = THHAT (k) − THHAT (k-1) (15)
Here, ddly is a dead time equivalent value, and Kdly is a constant.

図31に本実施形態による装置の構成を示し、図32に予測結果を示す。図31において、3111乃至3113は差分器を、3121乃至3123は適応フィルタを、3131は式14の演算を行うモジュールを、3101はファジィ推定器を示す。   FIG. 31 shows the configuration of the apparatus according to the present embodiment, and FIG. 32 shows the prediction result. In FIG. 31, 3111 to 3113 denote differentiators, 3121 to 3123 denote adaptive filters, 3131 denotes a module for performing the calculation of Expression 14, and 3101 denotes a fuzzy estimator.

つぎに、THにむだ時間を挿入する実施形態について説明する。最近の低速トルクの増大要求により、内燃機関のインマニボリュームが過大となっている場合が多い。このとき、THの変化に対する実際の吸気管圧力の変化にむだ時間dthが生じる。この場合、現在のTHを使用して実際の吸気管圧力の推定を行うと、実際の吸気管圧力の変化に対してTHの変化時刻が速すぎるため、予測演算にTHの情報を使用できなくなる。したがって、図33に示すように、予測演算に使用するTHにむだ時間dthを挿入することによってこの問題を解決する。図33において、3311乃至3313は差分器を、3341はむだ時間モジュールを、3301はファジィ推定器を示す。   Next, an embodiment in which a dead time is inserted into TH will be described. Due to recent demands for increasing low-speed torque, the intake manifold volume of the internal combustion engine is often excessive. At this time, a dead time dth occurs in the actual change in the intake pipe pressure with respect to the change in TH. In this case, if the actual intake pipe pressure is estimated using the current TH, the TH change time is too fast with respect to the actual intake pipe pressure change, and therefore the TH information cannot be used for the prediction calculation. . Therefore, as shown in FIG. 33, this problem is solved by inserting the dead time dth into TH used for the prediction calculation. 33, reference numerals 3311 to 3313 denote differentiators, 3341 denotes a dead time module, and 3301 denotes a fuzzy estimator.

本発明の適応フィルタを使用する実施形態の手順を図34の流れ図によって説明する。ステップS10でPBセンサが活性であるかチェックする。PBセンサが活性でなければステップS80に進み、PBのサンプル値にPBの代替値を代入する。さらに、ステップS90でPBの代替値であるPBのサンプル値をFZPBとして終了する。PBセンサが活性であれば、ステップS20に進み、PBセンサが正常であるかチェックする。正常でなければ、ステップS80に進む。正常であればステップS30に進み、THの推定値THHATを演算する(式14)。なお、THの代わりにTHHATを使用するのは電子制御スロットルの場合である。つぎに、ステップS40に進み、ΔPB、ΔΔPBおよびΔTHHATを演算する(式2、7、14)。さらに,ステップS50に進み、適応フィルタの演算を行う(式10乃至13)。つぎに、ステップS60においてΔFZPBの推定を行う。ステップS70において、PBのサンプル値にΔFZPBを加算してFZPBを求めて終了する。   The procedure of the embodiment using the adaptive filter of the present invention will be described with reference to the flowchart of FIG. In step S10, it is checked whether the PB sensor is active. If the PB sensor is not active, the process proceeds to step S80, and the substitute value of PB is substituted for the sample value of PB. Further, in step S90, the sample value of PB, which is an alternative value of PB, is ended as FZPB. If the PB sensor is active, the process proceeds to step S20 to check whether the PB sensor is normal. If not normal, the process proceeds to step S80. If normal, the process proceeds to step S30 to calculate an estimated value THHAT of TH (Formula 14). Note that THHAT is used instead of TH in the case of an electronically controlled throttle. Next, the process proceeds to step S40, and ΔPB, ΔΔPB and ΔTHHAT are calculated (equations 2, 7, and 14). In step S50, the adaptive filter is calculated (Equations 10 to 13). Next, in step S60, ΔFZPB is estimated. In step S70, ΔFZPB is added to the sample value of PB to obtain FZPB, and the process ends.

内燃機関の吸気部分の構成を示す。1 shows a configuration of an intake portion of an internal combustion engine. 従来の吸気管圧力予測アルゴリズムの挙動を示す。The behavior of a conventional intake pipe pressure prediction algorithm is shown. PB、TH、ΔPB、ΔTHおよびΔFZPBの間の関係を示す。The relationship between PB, TH, ΔPB, ΔTH and ΔFZPB is shown. 本発明の1実施形態による前件部メンバーシップ関数を示す。Fig. 5 shows an antecedent part membership function according to an embodiment of the present invention. 本発明の1実施形態による後件部メンバーシップ関数を示す。6 illustrates a consequent part membership function according to one embodiment of the present invention. 本発明の1実施形態によるファジィルールを示す。2 illustrates a fuzzy rule according to an embodiment of the present invention. ルール1が適用される状況を示す。The situation where rule 1 is applied is shown. ルール2が適用される状況を示す。The situation where rule 2 is applied is shown. ルール3が適用される状況を示す。The situation where rule 3 is applied is shown. ルール4が適用される状況を示す。The situation where rule 4 is applied is shown. ルール5が適用される状況を示す。The situation where rule 5 is applied is shown. ルール6が適用される状況を示す。The situation where rule 6 is applied is shown. ルール7が適用される状況を示す。The situation where rule 7 is applied is shown. ルール8が適用される状況を示す。The situation where rule 8 is applied is shown. ルール9が適用される状況を示す。The situation where rule 9 is applied is shown. 前件部メンバーシップ関数からルール6の適合度を求める方法を示す。A method for obtaining the conformity of rule 6 from the antecedent part membership function will be described. 後件部メンバーシップ関数からルール6の重みを求める方法を示す。A method of obtaining the weight of rule 6 from the consequent part membership function is shown. 本発明の1実施形態によるファジイ推定アルゴリズムを使用したFZPBの推定結果を示す。FIG. 6 shows an estimation result of FZPB using a fuzzy estimation algorithm according to an embodiment of the present invention. FIG. 本発明の別の実施形態による前件部メンバーシップ関数を示す。Fig. 5 illustrates an antecedent membership function according to another embodiment of the present invention. 本発明の別の実施形態によるファジィルールを示す。6 illustrates a fuzzy rule according to another embodiment of the present invention. 本発明の別の実施形態による後件部メンバーシップ関数を示す。6 illustrates a consequent membership function according to another embodiment of the present invention. ルール10が適用される状況を示す。The situation where rule 10 is applied is shown. ルール11が適用される状況を示す。The situation where rule 11 is applied is shown. 本発明の別の実施形態においてルール6が適用される状況を示す。FIG. 6 shows a situation where rule 6 is applied in another embodiment of the present invention. FIG. 本発明の別の実施形態においてルール9が適用される状況を示す。The situation where rule 9 is applied in another embodiment of the present invention is shown. 本発明の別の実施形態によるファジイ推定アルゴリズムを使用したFZPBの推定結果を示す。Fig. 6 shows the estimation results of FZPB using a fuzzy estimation algorithm according to another embodiment of the present invention. 本発明の実施形態によるファジイ推定器の構成を示す。2 shows a configuration of a fuzzy estimator according to an embodiment of the present invention. 本発明の実施形態による、適応フィルタを備えたファジイ推定器の構成を示す。2 shows a configuration of a fuzzy estimator with an adaptive filter according to an embodiment of the present invention. 適応フィルタを使用しない場合と使用した場合のFZPBの挙動を示す。The behavior of FZPB when the adaptive filter is not used and when it is used is shown. 電子制御スロットル使用時の状況を示す。Shows the situation when using electronically controlled throttle. 電子制御スロットルに対応したファジィ推定器の構成を示す。A configuration of a fuzzy estimator corresponding to an electronically controlled throttle is shown. 電子制御スロットルに対応したファジィ推定器の予測結果を示す。The prediction result of the fuzzy estimator corresponding to the electronic control throttle is shown. 過大インマニボリュームに対応したファジィ推定器の予測結果を示す。The prediction result of the fuzzy estimator corresponding to the excessive intake manifold volume is shown. 本発明の適応フィルタを使用する実施形態の手順を示す。Fig. 4 shows a procedure of an embodiment using the adaptive filter of the present invention.

符号の説明Explanation of symbols

1 スロットル弁
2 吸気管圧力センサ
2701,2801,3101、3301 ファジィ推定器
1 Throttle Valve 2 Intake Pipe Pressure Sensor 2701, 2801, 3101, 3301 Fuzzy Estimator

Claims (18)

現時点の吸気管圧力値の差分および現時点のスロットル弁開度値の差分を求めるステップと、
前記現時点の吸気管圧力値の差分の大きさおよび前記現時点のスロットル弁開度値の差分の大きさに基づいて定めたファジィルールを含むファジィ推定アルゴリズムにより吸気管圧力の予測差分を求めるステップと、
現時点の吸気管圧力値と前記吸気管圧力の予測差分の値を加算して吸気管圧力の予測値を求めるステップとを含み、
スロットル弁開度値とスロットル弁開度目標値を無駄時間要素と遅れ系を用いてモデル化し、当該モデルと当該目標値によって推定された値をスロットル弁開度値として用いる吸気管圧力予測方法。
Obtaining a difference between a current intake pipe pressure value and a current throttle valve opening value;
Obtaining a predicted difference of the intake pipe pressure by a fuzzy estimation algorithm including a fuzzy rule determined based on the magnitude of the difference in the current intake pipe pressure value and the difference in the current throttle valve opening value; and
Adding a current intake pipe pressure value and a predicted difference value of the intake pipe pressure to obtain a predicted value of the intake pipe pressure,
An intake pipe pressure prediction method in which a throttle valve opening value and a throttle valve opening target value are modeled using a dead time element and a delay system, and a value estimated from the model and the target value is used as a throttle valve opening value.
前記現時点の吸気管圧力値の差分の大きさを正、零、負に区分し、前記現時点のスロットル弁開度値の差分の大きさを正、零、負に区分して、それぞれが正、零、負のいずれかであるかによって定まる9個の領域に対してそれぞれファジィルールを定める請求項1に記載の吸気管圧力予測方法。   The magnitude of the difference in the current intake pipe pressure value is divided into positive, zero, and negative, and the magnitude of the difference in the current throttle valve opening value is divided into positive, zero, and negative, each being positive, The intake pipe pressure prediction method according to claim 1, wherein fuzzy rules are respectively defined for nine regions determined depending on whether they are zero or negative. 差分を求めるステップにおいて、さらに現時点の吸気管圧力値の2回差分を求め、吸気管圧力の予測差分を求めるステップにおいて、前記現時点の吸気管圧力値の差分の大きさ、前記現時点のスロットル弁開度値の差分の大きさおよび前記現時点の吸気管圧力値の2回差分の大きさに基づいてファジィルールを定める請求項1に記載の吸気管圧力予測方法。   In the step of obtaining the difference, a difference between the current intake pipe pressure values is further obtained, and in the step of obtaining a predicted difference in intake pipe pressure, the magnitude of the difference in the current intake pipe pressure value, the current throttle valve opening is determined. The intake pipe pressure prediction method according to claim 1, wherein a fuzzy rule is determined based on a magnitude of a difference in degree value and a magnitude of a two-time difference in the current intake pipe pressure value. 差分を求めるステップにおいて、さらに現時点の吸気管圧力値の2回差分を求め、吸気管圧力の予測差分を求めるステップにおいて、前記現時点の吸気管圧力値の2回差分の大きさを正、零、負に区分し、前記現時点の吸気管圧力値の差分の大きさ、前記現時点のスロットル弁開度値の差分の大きさおよび前記現時点の吸気管圧力値の2回差分の大きさのそれぞれが正、零、負のいずれかであるかによって定まる27個の領域に対してそれぞれファジィルールを定める請求項2に記載の吸気管圧力予測方法。   In the step of obtaining the difference, a difference between the current intake pipe pressure values is further obtained twice, and in the step of obtaining a predicted difference in the intake pipe pressure, the magnitude of the second difference in the current intake pipe pressure value is positive, zero, The difference between the current intake pipe pressure value, the current throttle valve opening value difference, and the current double difference of the intake pipe pressure value is positive. The intake pipe pressure prediction method according to claim 2, wherein fuzzy rules are respectively defined for 27 regions determined depending on whether they are zero, zero, or negative. ファジィ推定アルゴリズムの後件部メンバーシップ関数は、シングルトンの棒状関数である請求項1から4のいずれか一項に記載の吸気管圧力予測方法。   The intake pipe pressure prediction method according to any one of claims 1 to 4, wherein the consequent part membership function of the fuzzy estimation algorithm is a singleton bar function. 差分を求めるステップと予測差分を求めるステップとの間に、フィルタリングするステップをさらに含む請求項1から5のいずれか一項に記載の吸気管圧力予測方法。   The intake pipe pressure prediction method according to any one of claims 1 to 5, further comprising a filtering step between the step of obtaining the difference and the step of obtaining the prediction difference. フィルタリングが適応ファイルタによって行われる請求項6に記載の吸気管圧力予測方法。   The intake pipe pressure prediction method according to claim 6, wherein the filtering is performed by an adaptive filter. 現時点の吸気管圧力値の差分を求める差分器と、現時点のスロットル弁開度値の差分を求める差分器と、前記現時点の吸気管圧力値の差分および前記現時点のスロットル弁開度値の差分を入力として、前記現時点の吸気管圧力値の差分の大きさおよび前記現時点のスロットル弁開度値の差分の大きさに基づいて定めたファジィルールを含むファジィ推定アルゴリズムにより吸気管圧力の予測差分を求め、当該予測差分を出力とするファジィ推定器と、スロットル弁開度値とスロットル弁開度目標値を無駄時間要素と遅れ系を用いてモデル化したモデルとを含み、当該モデルと当該目標値によって推定された値をスロットル弁開度値として用いる吸気管圧力予測装置。   A differencer for obtaining a difference in the current intake pipe pressure value, a differencer for obtaining a difference in the current throttle valve opening value, a difference in the current intake pipe pressure value, and a difference in the current throttle valve opening value As an input, a predicted difference in the intake pipe pressure is obtained by a fuzzy estimation algorithm including a fuzzy rule determined based on the difference in the current intake pipe pressure value and the difference in the current throttle valve opening value. A fuzzy estimator that outputs the prediction difference, and a model in which the throttle valve opening value and the throttle valve opening target value are modeled using a dead time element and a delay system, depending on the model and the target value. An intake pipe pressure predicting device that uses an estimated value as a throttle valve opening value. ファジィ推定器において、吸気管圧力値の差分の大きさを正、零、負に区分し、スロットル弁開度値の差分の大きさを正、零、負に区分して、それぞれが正、零、負のいずれかであるかによって定まる9個の領域に対してそれぞれファジィルールを定める請求項8に記載の吸気管圧力予測装置。   In the fuzzy estimator, the difference in intake pipe pressure value is divided into positive, zero, and negative, and the difference in throttle valve opening value is divided into positive, zero, and negative. 9. The intake pipe pressure prediction apparatus according to claim 8, wherein fuzzy rules are respectively defined for nine regions determined depending on whether they are negative or negative. ファジィ推定器において、さらに現時点の吸気管圧力値の2回差分を入力とし、吸気管圧力値の差分の大きさ、スロットル弁開度値の差分の大きさおよび吸気管圧力値の2回差分の大きさに基づいてファジィルールを定める請求項8に記載の吸気管圧力予測装置。   In the fuzzy estimator, the difference between the current intake pipe pressure values and the difference between the intake pipe pressure values, the throttle valve opening value difference, and the intake pipe pressure value are calculated. The intake pipe pressure prediction apparatus according to claim 8, wherein a fuzzy rule is determined based on the size. ファジィ推定器において、さらに現時点の吸気管圧力値の2回差分を入力とし、前記現時点の吸気管圧力値の2回差分の大きさを正、零、負に区分し、前記現時点の吸気管圧力値の差分の大きさ、前記現時点のスロットル弁開度値の差分の大きさおよび前記現時点の吸気管圧力値の2回差分の大きさのそれぞれが正、零、負のいずれかであるかによって定まる27個の領域に対してそれぞれファジィルールを定める請求項9に記載の吸気管圧力予測装置。   In the fuzzy estimator, the difference between the current intake pipe pressure values is further inputted, the magnitude of the difference between the current intake pipe pressure values is divided into positive, zero, and negative, and the current intake pipe pressure is divided. Depending on whether the magnitude of the difference in the value, the magnitude of the difference in the current throttle valve opening value, and the magnitude of the two-time difference in the current intake pipe pressure value are positive, zero, or negative The intake pipe pressure prediction apparatus according to claim 9, wherein fuzzy rules are respectively defined for 27 fixed regions. ファジィ推定アルゴリズムの後件部メンバーシップ関数は、シングルトンの棒状関数である請求項8から11のいずれか一項に記載の吸気管圧力予測装置。   The intake pipe pressure prediction apparatus according to any one of claims 8 to 11, wherein the consequent part membership function of the fuzzy estimation algorithm is a singleton bar function. 入力をフィルタリングするためのフィルタをさらに備える請求項8から12のいずれか一項に記載の吸気管圧力予測装置。   The intake pipe pressure prediction apparatus according to any one of claims 8 to 12, further comprising a filter for filtering an input. フィルタが適応ファイルタである請求項13に記載の吸気管圧力予測装置。   The intake pipe pressure prediction apparatus according to claim 13, wherein the filter is an adaptive filter. 現時点の吸気管圧力値の差分、現時点の吸気管圧力値の2回差分および現時点のスロットル弁開度値の差分を求めるステップと、
前記現時点の吸気管圧力値の差分の大きさ、前記現時点のスロットル弁開度値の差分の大きさおよび前記現時点の吸気管圧力値の2回差分の大きさに基づいて定めたファジィルールを含むファジィ推定アルゴリズムにより吸気管圧力の予測差分を求めるステップと、
現時点の吸気管圧力値と前記吸気管圧力の予測差分の値を加算して吸気管圧力の予測値を求めるステップとを含む吸気管圧力予測方法。
Obtaining a difference in the current intake pipe pressure value, a difference in the current intake pipe pressure value twice and a difference in the current throttle valve opening value;
A fuzzy rule determined based on the magnitude of the difference between the current intake pipe pressure values, the magnitude of the difference between the current throttle valve opening values, and the magnitude of the second difference between the current intake pipe pressure values Obtaining a prediction difference of the intake pipe pressure by a fuzzy estimation algorithm;
An intake pipe pressure prediction method comprising: adding a current intake pipe pressure value and a predicted difference value of the intake pipe pressure to obtain a predicted value of the intake pipe pressure.
現時点の吸気管圧力値の差分および現時点のスロットル弁開度値の差分を求めるステップと、
前記吸気管圧力値の差分および前記スロットル弁開度値の差分をフィルタリングするステップと、
フィルタリングされた前記現時点の吸気管圧力値の差分の大きさおよびフィルタリングされた前記現時点のスロットル弁開度値の差分の大きさに基づいて定めたファジィルールを含むファジィ推定アルゴリズムにより吸気管圧力の予測差分を求めるステップと、
現時点の吸気管圧力値と前記吸気管圧力の予測差分の値を加算して吸気管圧力の予測値を求めるステップとを含む吸気管圧力予測方法。
Obtaining a difference between a current intake pipe pressure value and a current throttle valve opening value;
Filtering the difference between the intake pipe pressure values and the difference between the throttle valve opening values;
Prediction of intake pipe pressure by a fuzzy estimation algorithm including a fuzzy rule determined based on the magnitude of the difference between the filtered current intake pipe pressure values and the filtered difference between the current throttle valve opening values Obtaining a difference;
An intake pipe pressure prediction method comprising: adding a current intake pipe pressure value and a predicted difference value of the intake pipe pressure to obtain a predicted value of the intake pipe pressure.
現時点の吸気管圧力値の差分を求める差分器と、現時点の吸気管圧力値の2回差分を求める差分器と、現時点のスロットル弁開度値の差分を求める差分器と、前記現時点の吸気管圧力値の差分、前記現時点のスロットル弁開度値の差分および前記現時点の吸気管圧力値の2回差分を入力として、前記現時点の吸気管圧力値の差分の大きさ、前記現時点のスロットル弁開度値の差分の大きさおよび前記現時点の吸気管圧力値の2回差分の大きさに基づいて定めたファジィルールを含むファジィ推定アルゴリズムにより吸気管圧力の予測差分を求め、当該予測差分を出力とするファジィ推定器と、を含む吸気管圧力予測装置。   A differencer for obtaining a difference between the current intake pipe pressure values, a differencer for obtaining a difference between the current intake pipe pressure values twice, a differencer for obtaining a difference between current throttle valve opening values, and the current intake pipe The difference in the pressure value, the difference in the current throttle valve opening value, and the two-time difference in the current intake pipe pressure value are input, and the magnitude of the difference in the current intake pipe pressure value is determined. A fuzzy estimation algorithm including a fuzzy rule determined based on the magnitude of the difference in the degree value and the magnitude of the two-time difference in the current intake pipe pressure value is used to obtain a predicted difference in the intake pipe pressure, and the predicted difference is output as And a fuzzy estimator. 現時点の吸気管圧力値の差分を求める第1の差分器と、現時点のスロットル弁開度値の差分を求める第2の差分器と、前記現時点の吸気管圧力値の差分および前記現時点のスロットル弁開度値の差分を入力として、前記現時点の吸気管圧力値の差分の大きさおよび前記現時点のスロットル弁開度値の差分の大きさに基づいて定めたファジィルールを含むファジィ推定アルゴリズムにより吸気管圧力の予測差分を求め、当該予測差分を出力とするファジィ推定器と、前記第1の差分器と前記ファジィ推定器との間に配置された第1のフィルタと、前記第2の差分器と前記ファジィ推定器との間に配置された第2のフィルタと、を含む吸気管圧力予測装置。   A first subtractor for obtaining a difference between the current intake pipe pressure values; a second subtractor for obtaining a difference between the current throttle valve opening values; the difference between the present intake pipe pressure values and the current throttle valve; An intake pipe is obtained by a fuzzy estimation algorithm including a fuzzy rule determined based on the difference in the current intake pipe pressure value and the difference in the current throttle valve opening value. A fuzzy estimator that obtains a prediction difference of pressure, and outputs the prediction difference; a first filter disposed between the first differencer and the fuzzy estimator; and the second differencer; And a second filter disposed between the fuzzy estimator and an intake pipe pressure prediction device.
JP2007322117A 2007-12-13 2007-12-13 Method and apparatus for predicting intake pipe pressure of internal combustion engine Expired - Fee Related JP4291863B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007322117A JP4291863B2 (en) 2007-12-13 2007-12-13 Method and apparatus for predicting intake pipe pressure of internal combustion engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007322117A JP4291863B2 (en) 2007-12-13 2007-12-13 Method and apparatus for predicting intake pipe pressure of internal combustion engine

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2002315473A Division JP4104425B2 (en) 2002-10-30 2002-10-30 Method and apparatus for predicting intake pipe pressure of internal combustion engine

Publications (2)

Publication Number Publication Date
JP2008075663A JP2008075663A (en) 2008-04-03
JP4291863B2 true JP4291863B2 (en) 2009-07-08

Family

ID=39347966

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007322117A Expired - Fee Related JP4291863B2 (en) 2007-12-13 2007-12-13 Method and apparatus for predicting intake pipe pressure of internal combustion engine

Country Status (1)

Country Link
JP (1) JP4291863B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102066946B (en) * 2008-06-20 2016-08-31 味之素株式会社 The evaluation methodology of prostatosis
WO2009154296A1 (en) * 2008-06-20 2009-12-23 味の素株式会社 Female genital cancer evaluation method
JPWO2018066620A1 (en) 2016-10-04 2019-07-18 味の素株式会社 Evaluation method for pancreatic cancer, evaluation device, evaluation program, evaluation system, and terminal device

Also Published As

Publication number Publication date
JP2008075663A (en) 2008-04-03

Similar Documents

Publication Publication Date Title
EP1775451B1 (en) A method for estimating engine friction torque
JP5754514B2 (en) Supercharged engine control device
JP2008138630A (en) Control device of internal combustion engine
US20070100535A1 (en) Controller for controlling a plant
JPH10122017A (en) Engine control system
RU2719355C1 (en) Control device and control method
JP2004150424A (en) Controlling device of internal combustion engine
JP4432722B2 (en) Control device
JPH0431641A (en) Fuel injection control method of engine
JP4291863B2 (en) Method and apparatus for predicting intake pipe pressure of internal combustion engine
JP2007040266A (en) Suction air amount estimating device for internal combustion engine
WO2009107372A1 (en) Apparatus for controlling fuel injection amount for internal combustion engine
JP4104425B2 (en) Method and apparatus for predicting intake pipe pressure of internal combustion engine
JPH10184431A (en) Engine control system
CN108350820B (en) Engine control device for construction machine
JP3898114B2 (en) Intake air amount estimation method, estimation device, intake air amount control method and control device for internal combustion engine
JP4605049B2 (en) Control device for internal combustion engine
JP4536104B2 (en) Control device for internal combustion engine
JP2005105822A (en) Combustion state estimating device of internal combustion engine
JP2005069020A (en) Control device of internal combustion engine
JP2005330837A (en) Control device for internal combustion engine
US6176218B1 (en) Stabilizing function for torque based idle control
JP4702085B2 (en) Combustion state estimation device for internal combustion engine
US7546760B2 (en) Device for pressure-based load detection
CN107084058B (en) Intake air temperature estimation system for turbocharged engine

Legal Events

Date Code Title Description
TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20090331

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090403

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120410

Year of fee payment: 3

LAPS Cancellation because of no payment of annual fees