CN101454190B - 车辆控制装置 - Google Patents

Abstract

一种车辆控制装置，以使由车辆模型决定的标准状态量与实车的实际状态量之间的差(状态量偏差)接近于0的形式，由FB分配法则决定实车执行装置操作用控制输入和模型操作用控制输入，并通过这些控制输入来分别操作实车的执行机构和车辆模型。设定车辆模型的参数的值对应于实际车辆的动作状态而可变，以使驾驶操作量(操舵角)产生阶梯状变化时的标准状态量的衰减性比实际状态量的衰减性高。由此，可以不受实车的动作状态的影响而恰当地控制实车的执行机构以使与实际车辆的动作有关的状态量(实际状态量)接近于动态特性模型(车辆模型(16))下的车辆的状态量(标准状态量)。

Description

车辆控制装置

技术领域

本发明涉及一种汽车(发动机汽车)、混合动力汽车、两轮摩托车等具有多个车轮的车辆的控制装置。 

背景技术

以往，如在日本特开平11-91608号公报(以下称专利文献1)、特开2000-159138号公报(以下称专利文献2)所见到的那样，下述技术已为人所知：即，使用表示车辆的动态特性的模型来设定车辆的横摆比率、横向速度等的动作的状态量的目标值，控制实际装备于车辆上的执行机构以使实际的车辆的状态量追随于该目标值。 

在专利文献1、2中所见的技术中，通过具有被作为状态量的车辆的横摆比率和横向速度的动态模型，根据车速的测定值和方向盘的操舵角的测定值来设定目标横摆比率以及目标横向速度。然后算出根据目标横摆比率而被要求的车辆后轮的舵角的目标值(第一目标后轮舵角)和根据目标横向速度而被要求的车辆后轮的舵角的目标值(第二目标后轮舵角)。进一步将第一目标后轮舵角和第二目标后轮舵角进行线性结合而得的值(加权平均值)决定为后轮舵角的最终的目标值。随后，以使实际车辆的后轮的舵角追随该目标值的形式来控制后轮的操舵用的执行装置。 

另外，例如PCT国际公布公报WO2006/013922A1(以下称为专利文献3)中所见到的，如本发明的申请人所提议的那样，不仅以使实际车辆的状态量接近于该车辆的动态特性模型下的状态量的形式来控制实际车辆的执行机构，还以使动态特定模型下的状态量接近于实际车辆的状态量的形式来操作动态特性模型(将附加性的控制输入赋予动态特性模型)。 

然而，在对使与实际车辆的动作相关的状态量追随于车辆的动态特性模型下的状态量的情况进行控制时，对于方向盘的操舵角等的输入的变化，车辆的动态特性模型下的状态量一旦产生振动性的响应，追随于该响 应的实际车辆的状态量也会表现出振动性的响应。因此，为了尽可能保持实际车辆良好的动作特性，认为理想的是以下一种情况：即，动态特性模型下的状态量对于方向盘的操舵角等的输入的变化的响应特性(过渡响应特性)能够不受实际车辆的动作状态的影响，而尽可能成为衰减性高的一种响应特性。 

这里，本说明书中的“衰减性高”是指：使对象***的输入呈阶梯状变化时的状态量(控制量)的振动分量的振幅值的衰减的时间常数较短，即，该振动分量的振幅值的衰减速度较快。另外，所谓临界制动(criticaldumping)或过度制动(over dumping)的响应特性虽说是一种不带有振动分量的非振动性的响应特性，这里则将其看作衰减性为最高的一种特性。 

然而，所述专利文献1、2中所见到的技术由于以下的原因而难于充分地满足如上的要求。即，一般来说，实际车辆的横摆比率和横向速度等状态量对于方向盘操舵角的变化的响应特性在行驶速度比较高的情况下具有成为振动性的响应特性(状态量一边振动一边收敛于恒定值的一种特性)的倾向。 

因此，利用专利文献1、2中所见的技术，假设一旦将车辆的动态特性模型的响应特性设定为不受实际车辆的行驶速度的影响(在任意的行驶速度中)的衰减性高的一种特性、例如是临界制动或是过度制动的特性(非振动性的特性)时，特别是在车辆的行驶速度比较高的状态下，该动态特性模型下的车辆的状态量的响应特性与实际车辆的状态量的响应特性会产生较大的乖离。其结果，专利文献1、2中所见技术中作为动态特性模型下的状态量的横摆比率和横向速度分别与实际车辆的横摆比率和横向速度的差变大。进而，用于使实际车辆的状态量追随动态特性模型下的状态量的执行机构的要求操作量则容易变得过大。其结果，极易产生事实上在该执行机构的能力范围内无法满足其要求操作量的情况。并且，在该情况下，事实上难以使实际车辆的状态量追随动态特性模型下的状态量。 

因此，通过专利文献1、2中所见的技术，事实上无法将车辆动态特性模型下的状态量的响应特性以不受实际车辆的行驶速度的影响的形式而设定为非振动性的特性等、衰减性高的特性。即，事实上，只有在诸如车辆的行驶速度较低的情况等的车辆受限制的动作状态(实际车辆的状态 量的响应特性成为衰减性较高的一种特性时的动作状态)下，才能将车辆的动态特性模型下的状态量的响应特性设定为衰减性较高的特性。所以，专利文献1、2中所见的技术无法充分满足所述的要求。 

另一方面，在所述专利文献3中所见的技术中，除了根据实际车辆的状态量与动态特性模型下的状态量之间的偏差来操作实际车辆的执行机构以外，还操作车辆的动态特性模型，这样，该动态特性模型下的状态量可以接近实际车辆的状态量。即，该偏差不仅被反馈给实际车辆，而且还被反馈于动态特性模型。因此，能够抑制实际车辆的状态量与动态特性模型下的状态量之间的偏差变得过大。所以，如专利文献3中所见到的那样，在将实际车辆的状态量与动态特性模型下的状态量之间的偏差反馈给实际车辆和动态特性模型这两者时，可以认为即使在实际车辆的状态量的响应特性与动态特性模型下的车辆状态量的响应特性产生比较大的乖离的时候，也能抑制执行机构的要求操作量变得过大。进而，可以认为能够使实际车辆的状态量恰当地追随于动态特性模型下的状态量。 

所以，如专利文献3中所见到的那样，在将实际车辆的状态量与动态特性模型下的状态量之间的偏差反馈给实际车辆和动态特性模型这两者时，认为是可以将车辆动态特性模型下的状态量的响应特性以不受实际车辆的动作状态的影响的形式而设定成衰减性较高的一种响应特性。 

发明内容

本发明是鉴于上述背景而完成的，目的在于提供以下一种车辆控制装置，该控制装置使用了衰减性较高的动态特性模型，能够在不受车辆的动作状态的影响、并以使与实际车辆的动作有关的状态量接近于动态特性模型下的车辆的状态量的形式来控制车辆的执行机构，进而，能够将车辆的动作特性保持为良好的状态。 

为达到上述目的，本发明的车辆控制装置具有：驾驶操作量测定机构，其测定驾驶操作量，该驾驶操作量表示由操纵者操纵具有多个车轮的车辆时的该车辆的驾驶操作状态；执行机构，其设置于该车辆上，可操作所述车辆的规定的动作；实际状态量掌握机构，其测定或推定实际状态量向量，而该实际状态量向量是与实际车辆的规定的动作有关的多种状态量的值 的组合；模型状态量决定机构，其决定模型状态量向量，而该模型状态量向量是作为表示所述车辆的动态特性的模型而被事先设定的车辆模型下的车辆的所述多种状态量的值的组合；以及状态量偏差计算机构，其算出状态量偏差，该状态量偏差是所述测定或推定而得的实际状态量向量的各种状态量的值与所述被决定的模型状态量向量的各种状态量的值之间的偏差，并且，该车辆控制装置至少控制所述执行机构的动作以使所述状态量偏差接近于0，该车辆控制装置的特征在于，还具有：状态量偏差应动控制机构，该状态量偏差应动控制机构至少根据所述算出的状态量偏差并以使该状态量偏差接近于0的形式来决定实车执行装置操作用控制输入和车辆模型操作用控制输入，其中，所述实车执行装置操作用控制输入用于操作实际车辆的所述执行机构，所述车辆模型操作用控制输入用于操作所述车辆模型下的车辆的动作；和执行机构控制机构，该执行机构控制机构至少根据所述被决定的实车执行装置操作用控制输入来控制所述执行机构的动作，其中，所述模型状态量决定机构是至少根据所述测定而得的驾驶操作量和所述被决定的车辆模型操作用控制输入来决定所述模型状态量向量的一种机构，同时，该车辆控制装置还具有车辆模型特性设定机构，该车辆模型特性设定机构将所述车辆模型的至少一个参数的值设定为对应于所述实际车辆的动作状态而为可变，以使所述模型状态量向量的各种状态量的值的衰减性成为比所述实际状态量向量的各种状态量的值的衰减性要高的一种高衰减特性，其中，所述模型状态量向量的各种状态量的值的衰减性对应于将所述车辆模型操作用控制输入维持在0的状态时的所述驾驶操作量的阶梯状变化，所述实际状态量向量的各种状态量的值的衰减性对应于将所述实车执行装置操作用控制输入维持在0的状态时的所述驾驶操作量的阶梯状变化(第一发明)。 

根据上述第一发明，通过将所述车辆模型的至少一个参数设定为对应实际车辆的动作状态而可变，在车辆任意的动作状态下，都能够使所述模型状态量向量的各种状态量的值的衰减性成为高于实际状态量向量的各种状态量的值的衰减性的所述高衰减特性。在该情况下，特别是诸如实际状态量向量的各种状态量的值的响应特性是一种振动性的特性时的车辆的动作状态下，实际状态量向量的各种状态量的值的响应特性与模型状态 量向量的各种状态量的值的响应特性产生乖离的、所述车辆模型操作用控制输入被输入于所述车辆模型中。因此，抑制了所述状态量偏差变得过大，进而，抑制所述实车执行装置操作用控制输入变得过大。所以，能够使实际状态量向量的各种状态量的值恰当地接近于(追随于)模型状态量向量的各种状态量的值。并且，由于模型状态量向量的各种状态量的值的衰减性是一种所述高衰减性特性，因而该状态量的值会迅速地收敛于恒定值。其结果，能够抑制所述驾驶操作量变化时的实际状态量向量的各种状态量的值的振动，并能够将车辆的动作特性保持为良好的状态。 

因而，根据第一发明，使用作为一种衰减性高的动态特性模型的所述车辆模型，能够不受车辆的动作状态的影响而恰当地控制车辆的执行机构，以使与实际车辆的动作有关的状态量接近于该车辆模型下的车辆的状态量。进而，将车辆的动作特性保持为良好的状态。 

另外，在本说明书中，如上所述，临界制动或过度制动的响应特性是衰减性最高的一种响应特性。因此，在本发明中，实际状态量向量的各种状态量的值的响应特性(对应了在将所述实车执行装置操作用控制输入维持在0的状态下的所述驾驶操作量的阶梯状变化的一种响应特性)成为临界制动或过度制动的一种车辆的动作状态下的所述高衰减特性即指临界制动或是过度制动的特性。 

在所述的第一发明中，对应了所述驾驶操作量的阶梯状变化的所述模型状态量向量的各种状态量的值的响应特性最好是在车辆任意的动作状态下，尽可能地不带有振动分量。所以，所述高衰减特性优选为是临界制动或是过度制动的一种响应特性(第二发明)。 

根据该第二发明，能够有效地抑制或消除所述驾驶操作量产生变化时的实际状态量向量的各种状态量的值的振动。其结果，能够将车辆的动作特性保持在较为良好的状态。 

此外，在所述第一或第二发明中，所述多种状态量例如是包含了与所述车辆的侧滑以及绕横摆轴的旋转有关的两种以上的状态量的情况时，对应所述驾驶操作量的阶梯状变化的实际状态量向量的各种状态量的值的衰减性则根据实际车辆的行驶速度而进行变化。所以，在该情况下，所述车辆模型特性设定机构将所述车辆模型的参数的值设定为对应于所述实 际车辆的行驶速度而可变(第三发明)。 

根据该第三发明，能够可靠地将模型状态量向量的各种状态量的值的衰减性设定为所述高衰减特性。另外，所述两种以上的状态量可以不用个别地包括与车辆的侧滑有关的状态量和与绕横摆轴的旋转有关的状态量。例如，所述两种以上的状态量也可以是包括与车辆的侧滑有关的状态量和与绕横摆轴的旋转有关的状态量的两种状态量的线性结合的值。 

在所述第三发明中，在为了提高所述模型状态量向量的各种状态量的值的衰减性而对车辆模型的参数的值设定为对应车辆的行驶速度而可变的情况下，特别是在该行驶速度为较高的速度时的车辆的动作状态时，该车辆模型的固有值的绝对值容易变得过大。另外，车辆模型的“固有值”是指该车辆模型是线性模型的时的通常意义上的“固有值”。而当该车辆模型为非线性的模型时，车辆模型的“固有值”是指：将所述车辆模型操作用控制输入设为0，在所述模型状态量向量的各种状态量的任意的值中将该车辆模型近似为线性的模型(在该状态量的值的附近具有线性的模型)的固有值。 

如上所述在行驶速度较高时的车辆的动作状态下，车辆模型的固有值的绝对值容易变得过大。因此，车辆的行驶速度在较高的速度范围内时，恐怕有对应了驾驶操作量变化的模型状态量向量的各种状态量的值的响应性(驾驶操作量刚变化后的各种状态量的值的变化速度)会比实际状态量向量的各种状态量的值的响应性要高得多的可能性。该情况下，依据所述模型操作用控制输入，可以将所述状态量偏差抑制为较小的状态量偏差。然而，特别是将模型操作用控制输入相对于该状态量偏差的增益设定为较小的增益时，驾驶操作量刚变化后的状态量偏差会变得过大，所述实车执行装置操作用控制输入有可能变得过大。 

因而，在第三发明中，优选为所述车辆模型特性设定机构设定所述车辆模型的参数的值并使其对应于所述行驶速度为可变，以使至少所述行驶速度比规定速度高的时候，所述车辆模型的固有值的绝对值成为规定值以下的数值(第四发明)。 

根据该第四发明，除了将所述模型状态量向量的各种状态量的值的衰减性设为所述高衰减特性以外，至少在所述行驶速度比规定速度要高的时 候，以使车辆模型的固有值的绝对值成为规定值以下的数值的形式，该车辆模型的参数的值对应该行驶速度而被设定。因此，能够提高对应了驾驶操作量的阶梯状变化的、模型状态量向量的各种状态量的值的衰减性的同时，能够防止该状态量的值的响应性相对实际状态量向量的各种状态量的值的响应性而变得过高。其结果，能够比较可靠地防止所述状态量偏差变得过大。进而，能够比较可靠地防止所述实车执行装置操作用控制输入变得过大。所以，能够较恰当地执行使实际状态量向量的各种状态量的值追随于模型状态量向量的各种状态量的值的控制。 

另外，在第四发明中，由于是对应行驶速度来操作模型状态量向量的各种状态量的值的衰减性和响应性，因而，对应该行驶速度而设定为可变的车辆模型的参数需要有2个以上。 

此外，在所述第一至第四发明中，优选为所述车辆模型特性设定机构设定所述车辆模型的参数值并使其可变，以使下述驾驶操作量的恒定值与所述模型状态量向量的各种状态量的恒定值的关系保持为固定的关系，其中，所述驾驶操作量的恒定值是将所述车辆模型操作用控制输入维持在0的状态下使所述驾驶操作量产生阶梯状变化时的该驾驶操作量的恒定值(第五发明)。 

根据该第五发明，使所述驾驶操作量产生阶梯状变化时的该驾驶操作量的恒定值与所述模型状态量向量的各种状态量的值的恒定值的关系即使在车辆模型的稳定特性设定所述车辆模型的参数的值为可变(不受该参数的值的变更的影响)的情况下也会被保持为固定。因此，可以使该车辆模型的稳定特性成为与实际车辆的稳定特性(具体是指，将所述实车执行装置操作用控制输入维持在0的状态下的驾驶操作量与所述实际状态量向量的各种状态量的恒定值的关系)大致相同的特性。其结果，可以将实际车辆的动作时的稳定状态下(在驾驶操作量为恒定、且路面等的环境状态为恒定或相同的情况下，等到过渡性的动作消失为止经过了足够长的时间的一种状态)的所述状态量偏差一直控制在极其小的范围内。进而，将实际车辆的稳定状态下的所述实车执行装置操作用控制输入控制在最小限度内，能够防止所述执行机构进行过剩的动作。

附图说明

图1是方框图，表示本发明实施方式中车辆的概略构成。 

图2是功能方框图，表示本发明第一实施方式中车辆所具有的控制装置概略的整体控制处理功能。 

图3是示意图，表示第一实施方式中的标准动态特性模型(车辆模型)下的车辆构造。 

图4(a)是例示图，例示分别在实际车辆和模型车辆中的横摆比率的响应特性；图4(b)是例示图，例示分别在实际车辆和车辆模型中的车辆重心点的侧滑角的响应特性。 

图5是示意图，表示第一实施方式中的模型特性调整用参数k2的设定例。 

图6是功能方框图，表示第一实施方式中的标准操作量决定部的处理功能。 

图7是说明图，用于说明第一实施方式中的标准操作量决定部所具有的离心力过大防止限幅器的处理。 

图8是说明图，用于说明第一实施方式中的离心力过大防止限幅器的处理的另一例。 

图9是说明图，用于说明第一实施方式中的离心力过大防止限幅器的处理的另一例。 

图10是功能方框图，表示第一实施方式中由标准操作量决定部决定第二限制完毕前轮舵角δf_ltd2的处理的另一例。 

图11是功能方框图，表示第一实施方式中FB分配法则部的处理功能。 

图12是功能方框图，表示第一实施方式中假想外力决定部的处理的另一例。 

图13是说明图，用于说明第一实施方式中γβ限幅器处理的另一例。 

图14是功能方框图，表示第一实施方式中执行装置动作FB目标值决定部的处理。 

图15是说明图，用于说明第一实施方式中执行装置动作FB目标值决定部的处理中所使用的变量。 

图16(a)、图16(b)是示意图，表示第一实施方式中执行装置动作FB目标值决定部的处理中所使用的分配增益的设定例。 

图17(a)～(e)是例示图，例示了第一实施方式的执行装置动作FB目标值决定部的处理的另一例中所使用的图表。 

图18(a)～(e)是例示图，例示了第一实施方式中执行装置动作FB目标值决定部的处理的另一例中所使用的图表。 

图19是功能方框图，表示第一实施方式中FF法则的处理。 

图20是功能方框图，表示第一实施方式中执行装置动作目标值合成部的处理。 

图21是流程图，表示第一实施方式中执行装置动作目标值合成部所具有的最佳目标第n轮驱动/制动力决定部的处理。 

图22是功能方框图，表示第一实施方式中执行装置动作目标值合成部所具有的最佳目标自动舵角决定部的处理。 

图23是表示第三实施方式中模型车辆的固有值与车辆行驶速度的关系的图表。 

图24(a)、(b)是示意图，表示第三实施方式中模型特性调整用参数k1、k2的设定例。 

具体实施方式

下面，说明本发明的车辆控制装置的实施方式。 

首先，参照图1，说明本说明书实施方式中的车辆的概略构成。图1是表示该车辆的概略构成的框图。另外，本说明书实施方式中所例举的车辆具有4个车轮(车辆前后各2个车轮)。由于汽车构造本身可以是公知技术中的构造，因此在本说明书中省略其详细图示以及说明。 

如图1所示，车辆1(汽车)与公知的普通汽车相同，具有：驱动/制动装置3A(驱动/制动***)，其将转动驱动力(作为车辆1推进力的转动力)施加给4个车轮W1、W2、W3、W4中的驱动轮，或者将制动力(作为车辆1制动力的转动力)施加给各车轮W1～W4；转向装置3B(转向***)，其操舵4个车轮W1～W4中的操舵轮；悬架装置3C(悬架***)，其将车身1B弹性支承在4个车轮W1～W4上。车轮W1、W2、W3、W4分别是车辆1的左前方、右前方、左后方、右后方的车轮。另外，驱动轮和操舵 轮在本说明书所说明的实施方式中是2个前轮W1、W2。因此，后轮W3、W4为从动轮，且不为操舵轮。 

其中，驱动轮可以是2个后轮W3、W4，或者也可以是前轮W1、W2及后轮W3、W4两者(4个车轮W1～W4)。另外，操舵轮也可以不仅是2个前轮W1、W2，还可以包括后轮W3、W4。 

这些装置3A、3B、3C具有操作车辆1动作的功能。例如，驱动/制动装置3A主要具有操作车辆1行进方向的动作(车辆1行进方向的位置、速度、加速度等)的功能。转向装置3B主要具有操作车辆1横摆方向的转动动作(车辆1横摆方向的姿势、角速度、角加速度等)的功能。悬架装置3C主要具有操作车辆1车身1B的前后方向及左右方向的动作(车身1B的前后方向及左右方向的姿势等)、或者车身1B上下方向上的动作(车体1B距路面的高度(车身1B相对于车轮W1～W4的上下方向上的位置)等)的功能。另外，在本说明书中，车辆1或车体1B的“姿势”是表示空间性的朝向。 

作为补充，一般情况下，车辆1在转弯等时，会产生车轮W1～W4侧滑。而且，该侧滑受到车辆1操舵轮的舵角、车辆1横摆比率(横摆方向的角速度)、各车轮W1～W4驱动/制动力等的影响。因此，驱动/制动装置3A和转向装置3B也具有操作车辆1横向(左右方向)的平移动作的功能。另外，车轮的“驱动/制动力”是表示：在从路面作用于该车轮的路面反力之中，该车轮前后方向(具体而言，是该车轮的转动面(通过车轮中心点而与该车轮的转轴正交的面)与路面或水平面间的交线的方向)的平移力分量。另外，将路面反力中车轮的宽度方向(与车轮的转轴平行的方向)的平移力分量称为“横向力”，将路面反力中与路面或水平面垂直的方向的平移力分量称为“着地负荷”。 

驱动/制动装置3A虽省略了详细图示，具体而言，其具有：驱动***和将制动力施加给各车轮W1～W4的刹车装置(制动***)，其中，所述驱动***由作为车辆1的动力产生源(车辆1传动力产生源)的发动机(内燃机)和将该发动机的输出(转动驱动力)传递给车轮W1～W4中的驱动轮的动力传递***构成。动力传递***中包括变速装置、差动齿轮装置等。 

另外，实施方式中说明的车辆1虽然是以发动机作为动力产生源的车 辆，但也可以是以发动机和电动机作为动力产生源的车辆(所谓并行型的混合动力车)，或者是以电动机作为动力产生源的车辆(所谓电动汽车、或者系列型的混合动力车)。 

并且，作为供驾驶者操纵车辆1(汽车)所操作用的操作器5(人为操作的操作器)，在车辆1的驾驶室内设置了驾驶盘(方向盘)、油门踏板、刹车踏板、变速杆等。另外，操作器5的各要素省略了图示。 

操作器5中的驾驶盘与所述转向装置3B的动作相关连。即，通过转动操作该驾驶盘，转向装置3B则与其对应进行动作，来操舵车轮W1～W4中的操舵轮W1、W2。 

操作器5中的油门踏板、刹车踏板以及变速杆与所述驱动/制动装置3A的动作相关连。即，根据油门踏板的操作量(踩入量)，发动机所具有的调节阀的开度发生变化，调整发动机的空气吸入量及燃料喷射量(进而调整发动机的输出)。另外，根据刹车踏板的操作量(踩入量)，刹车装置进行动作，与刹车踏板的操作量对应的制动力矩被施加给各车轮W1～W4。此外，通过操作变速杆，变速装置的变速比等该变速装置的动作状态发生变化，进行从发动机传递给驱动轮的驱动转矩的调整等。 

另外，通过适当的传感器测定由驾驶者(车辆1的操作者)操作的驾驶盘等各操作器5的驾驶操作状态，其中所述传感器省略了图示。以下，称该驾驶操作状态的测定值(传感器的测定输出)为驾驶操作输入。该驾驶操作输入包括：作为驾驶盘转角的转向角、作为油门踏板的操作量的油门踏板操作量、作为刹车踏板的操作量的刹车踏板操作量、以及作为变速杆的操作位置的换档位置的测定值。输出该驾驶操作输入的传感器相当于本发明中的驾驶操作量测定机构。 

在本说明书的实施方式中，所述驱动/制动装置3A以及转向装置3B为：不仅仅根据所述驾驶操作输入，还能根据该驾驶操作输入以外的主要因素(车辆1动作状态和环境状态等)可以能动地控制所述驱动/制动装置3A以及转向装置3B动作(进而控制车辆1的动作)的装置。此处，“可以能动地控制”表示可以将装置3A、3B的动作控制为：对与所述驾驶操作输入对应的基本动作(对应于驾驶操作输入而决定的基本目标动作)进行了修正后的动作。

具体而言，驱动/制动装置3A是具有下述功能的驱动/制动装置，即，该功能为：关于前轮W1、W2一组和后轮W3、W4一组中的至少任意一组，可以经该驱动/制动装置3A所具有的液压执行装置、电动机、电磁控制阀等执行装置，能动地控制左侧车轮W1、W3的驱动/制动力与右侧车轮W2、W4的驱动/制动力间的差或比率(以下，将该控制功能称为左右动力分配控制功能)。 

更具体而言，在本说明书的实施方式中，驱动/制动装置3A为下述驱动/制动装置：经所述刹车装置所具有的执行装置可以能动地控制由该刹车装置的动作而作用于各车轮W1～W4上的驱动/制动力(具体为车辆1的制动方向的驱动/制动力)的装置(即为在按照刹车踏板的操作量而被决定的基本驱动/制动力的基础之上可增加或可减少地控制由刹车装置作用于各车轮W1～W4上的驱动/制动力的装置)。因此，驱动/制动装置3A即为如下装置：经执行装置可以能动地控制在前轮W1、W2一组和后轮W3、W4一组的两组中由刹车装置决定的左侧车轮W1、W3的驱动/制动力与右侧车轮W2、W4的驱动/制动力间的差或比率的装置(在前轮W1、W2一组和后轮W3、W4一组的两组中具有左右动力分配控制功能的驱动/制动装置)。 

另外，驱动/制动装置3A除了具有能动地控制由刹车装置的动作决定的各车轮W1～W4的驱动/制动力的功能之外，还可以具有：经该驱动***所具备的执行装置，能动地控制通过驱动/制动装置3A的驱动***的动作而作用于驱动轮即前轮W1、W2上的驱动/制动力之差或比率的功能。 

作为具有这种左右动力分配控制功能的驱动/制动装置3A，可以使用公知的装置。 

作为补充，如上所述具有左右动力分配控制功能的驱动/制动装置3A根据其控制功能，还可以具有：能动地操作车辆1的横摆方向上的转动动作或横向上的平移动作的功能。 

另外，该驱动/制动装置3A除了包括有与左右动力分配控制功能有关的执行装置以外，还包括：用于产生刹车装置的制动转矩的执行装置、驱动发动机的节流阀的执行装置、驱动燃料喷射阀的执行装置以及执行变 速装置的变速驱动的执行装置等。 

另外，所述转向装置3B例如是下述这样的转向装置，其除了具有根据驾驶盘的转动操作、经齿轮齿条传动等操舵机构来机械性地操舵作为操舵轮的前轮W1、W2的功能之外，还可以是根据需要可以通过电动机等执行装置辅助性地操舵前轮W1、W2的转向装置(该转向装置可以在根据驾驶盘的转动角而机械性确定的舵角的基础上，可增加或可减少地控制前轮W1、W2的舵角的转向装置)。或者，转向装置3B是只使用执行装置的驱动力来操舵前轮W1、W2的转向装置(所谓线控转向的转向装置)。因此，转向装置3B是经执行装置可以能动地控制前轮W1、W2的舵角的转向装置(以下称主动式转向装置)。 

当转向装置3B是除了根据驾驶盘的转动操作而机械地操舵操舵轮W1、W2之外，还通过由执行装置辅助性地操舵操舵轮的主动式转向装置(以下，将这种主动式转向装置称为执行装置辅助型的转向装置)的情况下，通过驾驶盘的转动操作而机械地决定的操舵轮舵角和通过执行装置的动作决定的舵角(舵角的补正量)的合成角成为操舵轮的舵角。 

另外，当转向装置3B是只使用执行装置的驱动力而进行操舵轮W1、W2操舵的主动式转向装置(以下，称这种主动式转向装置为执行装置驱动型的转向装置)的情况下，至少根据转向角的测定值来决定操舵轮的舵角的目标值，并控制执行装置以使操舵轮的实际的舵角成为其目标值。 

作为这种经执行装置而可以能动地控制操舵轮W1、W2舵角的转向装置3B(主动式转向装置)，使用公知的技术即可。 

另外，本说明书的实施方式中的转向装置3B虽然是经执行装置可以能动地控制前轮W1、W2舵角的主动式转向装置，但也可以是根据驾驶盘的转动操作，只进行前轮W1、W2的机械性的操舵(以下称为机械式转向装置)。另外，在以所有车轮W1～W4作为操舵轮的车辆中，转向装置可以是经执行装置能动地控制前轮W1、W2及后轮W3、W4两者的舵角的转向装置。或者，该转向装置也可以是如下转向装置，即，只通过齿轮齿条传动等机械性的机构来进行与驾驶盘的转动操作对应的前轮W1、W2的操舵，并同时经执行装置只是能动地控制后轮W3、W4的舵角。 

所述悬架装置3C在本说明书的实施方式中，例如是根据车辆1的动 作而被动地进行动作的悬架装置。 

其中，悬架装置3C可以是下述的悬架装置，即，例如经电磁控制阀或电动机等执行装置，可变地控制设置在车身1B与车轮W1～W4之间的减振器的衰减力、硬度等。或者悬架装置3C也可以是下述的悬架装置，即，通过液压缸或气压缸可以直接控制悬架(悬架装置3C的弹簧等机构部分)的冲程(车身1B与各车轮W1～W4间的上下方向的位移量)或者在车身1B与车轮W1～W4之间产生的悬架上下方向的伸缩力(所谓电子控制悬架)。当悬架装置3C是如上所述的可以控制减振器的衰减力和硬度、悬架的冲程或伸缩力的悬架装置(以下称为主动式悬架装置)时，该悬架装置3C可以能动地控制其动作。 

在以下的说明中，在驱动/制动装置3A、转向装置3B以及悬架装置3C之中，有时会将如前所述那样可以能动地控制动作的装置统称为执行机构3。在本说明书的实施方式中，该执行机构3包括驱动/制动装置3A以及转向装置3B。另外，当悬架装置3C是主动式悬架装置时，执行机构3也包括该悬架装置3C。 

另外，车辆1具有控制装置10，该控制装置10根据所述驾驶操作输入等来决定所述各执行机构3所具有的执行装置的操作量(针对执行装置的控制输入。以下称执行装置操作量)，并通过该执行装置操作量来控制各执行机构3的动作。该控制装置10是由包括微机等在内的电路单元构成。该控制装置10从操作器5的传感器将所述驾驶操作输入进行输入，并同时将车辆1的行驶车速、横摆比率等车辆1的状态量测定值和车辆1的行驶环境的信息等从未图示的各种传感器进行输入。并且，该控制装置10，依据这些输入，以规定的控制处理周期逐步决定执行装置操作量，并逐步控制各执行机构3的动作。 

以上是本说明书的实施方式中的车辆1(汽车)的整体概略性构成。该概略性构成在以下说明的任一实施方式中均相同。 

作为补充，在本说明书的实施方式中，所述驱动/制动装置3A、转向装置3B以及悬架装置3C中的相当于本发明的执行装置(适用于本发明进行动作控制的执行装置)的是驱动/制动装置3A、或者该驱动/制动装置3A及转向装置3B。

另外，控制装置10根据其控制处理功能来实现本发明中的各种功能。 

(第一实施方式) 

下面，参照图2说明第一实施方式中的控制装置10的概略的控制处理。图2是表示控制装置10整体大致的控制处理功能的框图。此外，在以后的说明中，称实际车辆1为实车1。 

图2中除去了实车1后的部分(更准确地说，是指除去实车1和后述的传感器/推定器12所包含的传感器之后的部分)即为控制装置10的主要控制处理功能。图2中的实车1具有所述的驱动/制动装置3A、转向装置3B、以及悬架装置3C。 

如图所示，作为主要的处理功能部，控制装置10具有传感器/推定器12、标准操作量决定部14、标准动态特性模型16、减法器18、反馈分配则部(FB分配法则部)20、前馈则部(FF法则)22、执行装置动作目标值合成部24以及执行装置驱动控制装置26。另外，图2中实线的箭头符号表示对于各处理功能部的主要输入，虚线的箭头符号表示对于各处理功能部的辅助性输入。 

控制装置10以规定的控制处理周期执行这些处理功能部的处理，在每一该控制处理周期，逐步决定执行装置操作量。而且，根据其执行装置操作量，来逐步控制实车1的执行机构3的动作。 

以下，说明控制装置10的各处理功能部的概要和整体处理的概要。另外，以下关于由控制装置10的各控制处理周期所决定的变量值，将由现在(最新的)的控制处理周期的处理最终得到的数值称作此次值，将由前次控制处理周期的处理最终得到的值称作前次值。 

控制装置10在各控制处理周期，首先，通过传感器/推定器12测定或推定出实车1的状态量及实车1的行驶环境的状态量。在本实施方式中，传感器/推定器12的测定对象或推定对象包括例如：实车1横摆方向的角速度即横摆比率γact、实车1行驶速度Vact(对地速度)、实车1重心点的侧滑角即车辆重心点侧滑角βact、实车1前轮W1、W2的侧滑角的前轮侧滑角βf_act、实车1后轮W3、W4的侧滑角即后轮侧滑角βr_act、从路面作用于实车1各车轮W1～W4上的反向力即路面反力(驱动/制动力、 横向力、着地负荷)、实车1各车轮W1～W4的打滑比、实车1前轮W1、W2的舵角δf_act。 

在这些测定对象或推定对象中，车辆重心点侧滑角βact是指从上方观看实车1时(水平面上)的该实车1的行驶速度Vact的矢量相对于实车1前后方向所成的角度。另外，前轮侧滑角βf_act是指从上方观看实车1时(水平面上)的前轮W1、W2的行进速度矢量相对于前轮W1、W2前后方向所成的角度。另外，后轮侧滑角βr_act是指从上方观看实车1时(水平面上)的后轮W3、W4的行进速度矢量相对于后轮W3、W4前后方向所成的角度。另外，舵角δf_act，是从上方观看实车1时(水平面上)的前轮W1、W2的转动面相对于实车1前后方向所成的角度。 

另外，前轮侧滑角βf_act，虽然可以针对各个前轮W1、W2进行测定或推定，但也可以有代表性地测定或推定任意一方的前轮W1或W2的侧滑角作为βf_act，或者，还可以测定或推定两者的侧滑角的平均值作为βf_act。关于后轮侧滑角βr_act也同样。 

而且，作为传感器/推定器12的推定对象，包括有实车1车轮W1～W4和跟其接触的实际路面间的摩擦系数(以下，将该摩擦系数的推定值称为推定摩擦系数μestm)。另外，为了不使推定摩擦系数μestm产生频繁的变动，理想的是在摩擦系数的推定处理中具有低通特性的滤波处理等。 

传感器/推定器12具有被搭载在实车1上的各种传感器，用于测定或推定出上述的测定对象或推定对象。作为这些传感器，例如包括：测定实车1角速度的比率传感器、测定实车1前后方向及左右方向的加速度的加速度传感器、测定实车1行驶速度(对地速度)的速度传感器、测定实车1各车轮W1～W4的转速的转速传感器、测定从路面作用于实车1各车轮W1～W4上的路面反力的力传感器等。 

在该情况下，传感器/推定器12，关于在其测定对象或推定对象中通过搭载在实车1上的传感器无法直接测定的推定对象，则依据与该推定对象有关连的状态量的测定值、或者控制装置10决定的执行装置操作量的值或目标值，通过观测器等来进行推定。例如，依据搭载在实车1上的加速度传感器的测定值等来推定车辆重心点侧滑角βact。另外，例如依据加速度传感器的测定值等，并利用公知的手法来推定摩擦系数。

作为补充，传感器/推定器12具有作为本发明的实际状态量掌握机构的功能。在本实施方式中，作为与车辆的动作有关的状态量的种类，使用了车辆的横摆比率和车辆重心点侧滑角。此时，横摆比率具有作为与车辆的横摆方向的转动运动相关的状态量的意思，车辆重心点侧滑角具有作为与车辆侧滑(车辆横向的平移动作)有关的状态量的意思。所以，作为本发明中的实际状态量向量，所述横摆比率γact和车辆重心点侧滑角βact这组由传感器/推定器12进行测定或推定。 

以后，在由传感器/推定器12测定或推定出的实车1的状态量等名称中时常会标注“实”字。例如，将实车1的横摆比率γact、实车1的行驶速度Vact、实车1的车辆重心点侧滑角βact分别称为实横摆比率γact、实行驶速度Vact、实车辆重心点侧滑角βact。 

接着，控制装置10通过标准操作量决定部14，来决定作为对后述的标准动态特性模型16的输入的标准模型操作量。此时，由所述操作器5的传感器测定出的驾驶操作输入被输入给标准操作量决定部14，该标准操作量决定部14至少依据该驾驶操作输入来决定标准模型操作量。 

进一步具体而言，在本实施方式中，标准操作量决定部14决定的标准模型操作量即为后述的标准动态特性模型16下的车辆前轮的舵角(以下，称模型前轮舵角)。为了决定该模型前轮舵角，所述驾驶操作输入之中的转向角θh(此次值)被作为主要的输入量输入给标准操作量决定部14。并且，由传感器/推定器12测定或推定出的实行驶速度Vact(此次值)及推定摩擦系数μestm(此次值)、和标准动态特性模型16下的车辆的状态量(前次值)被输入给标准操作量决定部14。然后，标准操作量决定部14依据这些输入来决定模型前轮舵角。另外，模型前轮舵角基本上只要根据转向角θh来决定即可。但是，在本实施方式中，对输入给标准动态特性模型16的模型前轮舵角加了必要的限制。为了施加该限制，标准操作量决定部14中除了转向角θh以外还输入Vact、μestm等。 

作为补充，标准模型操作量的种类一般情况下是依靠于标准动态特性模型16的形态、或者通过该标准动态特性模型16想要决定的状态量的种类。另外，标准动态特性模型16也可以包含标准操作量决定部14。当在标准动态特性模型16是以驾驶操作输入作为必要输入来构成的情况下， 可以省略标准操作量决定部14。 

然后，控制装置10通过标准动态特性模型16来决定标准状态量并予以输出，其中该标准状态量是指作为实车1的标准的动作(以后称为标准动作)的状态量。标准动态特性模型16是表示车辆的动态特性且是被预先确定的模型，其依据包括所述标准模型操作量在内的必要的输入，并来逐步决定标准动作的状态量(标准状态量)。该标准动作表示：基本上是驾驶者所希望的实车1的理想动作或与之相近的动作。 

在该情况下，由标准操作量决定部14决定的标准模型操作量、和由后述的FB分配法则部20决定的标准动态特性模型16的操作用控制输入(反馈控制输入)Mvir、Fvir等被输入给标准动态特性模型16。随后，标准动态特性模型16根据这些输入量来决定标准动作(进而决定标准状态量的时序)。 

进一步具体而言，在本实施方式中，通过标准动态特性模型16决定并输出的标准状态量，是指与车辆的横摆方向的转动动作有关的标准状态量以及与车辆横向的平移动作(侧滑动作)有关的标准状态量的一组。与车辆的横摆方向的转动动作有关的标准状态量例如是横摆比率的标准值γd(以后，有时称作标准横摆比率γd)，与车辆的横向平移动作有关的标准状态量例如是车辆重心点侧滑角的标准值βd(以后，有时会称为标准车辆重心点侧滑角βd)。为了在每一控制处理周期逐步决定这些标准状态量γd、βd，输入作为标准模型操作量的所述模型前轮舵角(此次值)、和所述反馈控制输入Mvir、Fvir(前次值)。这种情况下，在本实施方式中，使标准动态特性模型16下的车辆的行驶速度与实行驶速度Vact一致。因此，由传感器/推定器12测定或推定出的实行驶速度Vact(此次值)也被输入给标准动态特性模型16。而且，标准动态特性模型16依据这些输入，来决定该标准动态特性模型16下的车辆的横摆比率及车辆重心点侧滑角，并将之作为标准状态量γd、βd而予以输出。 

另外，因为实车1行驶环境(路面状况等)的变化(在标准动态特性模型16没有被考虑到的变化)、标准动态特性模型16的模型化误差、或者传感器/推定器12的测定误差或推定误差等原因，会导致实车1动作与标准动作相背离(相乖离)，为了防止这一问题发生(使标准动作接近于实车1的 动作)，输入给标准动态特性模型16的反馈控制输入Mvir、Fvir是附加性地输入给标准动态特性模型16的反馈控制输入。该反馈控制输入Mvir、Fvir在本实施方式中，是假想性地作用于标准动态特性模型16下的车辆上的假想外力。该假想外力Mvir、Fvir中的Mvir，是围绕标准动态特性模型16下的车辆重心点而作用的横摆方向的假想性的力矩，Fvir是作用于该重心点的横向的假想性的平移力。 

作为补充，所述标准状态量γd、βd这一组相当于本发明中的模型状态量向量，标准动态特性模型16相当于本发明中的车辆模型。而且，根据标准操作量决定部14及标准动态特性模型16的处理，来构成本发明中的模型状态量决定机构。在本实施方式中，在标准动态特性模型16处理中包含有作为本发明的车辆模型特性设定机构的功能。 

此后，控制装置10通过减法器18计算出状态量偏差，该状态量偏差是指由传感器/推定器12测定或推定出的实际状态量(与标准状态量相同种类的实际状态量)、与由标准动态特性模型16决定的标准状态量间的差。 

进一步具体而言，在减法器18，求解出实横摆比率γact及实车辆重心点侧滑角βact的各自的值(此次值)、与由标准动态特性模型16决定的标准横摆比率γd及标准车辆重心点侧滑角βd的各自的值(此次值)之间的差γerr(＝γact—γd)、βerr(＝βact—γd)，并将其作为状态量偏差。 

作为补充，根据减法器18的处理来构成本发明中的状态量偏差计算机构。 

之后，控制装置10将如上所述求得的状态量偏差γerr、βerr输入给FB分配法则部20，通过该FB分配法则部20来决定：作为标准动态特性模型16的操作用反馈控制输入的所述假想外力Mvir、Fvir、和作为实车1的执行机构3的操作用反馈控制输入的执行装置动作反馈目标值(执行装置动作FB目标值)。 

另外，在本实施方式中，执行装置动作FB目标值包括：与驱动/制动装置3A的刹车装置的动作有关的反馈控制输入(进一步具体而言，是指操作通过该刹车装置的动作而作用于各车轮W1～W4上的驱动/制动力的反馈控制输入)。或者，执行装置动作FB目标值除了包括与驱动/制动装置3A的动作有关的反馈控制输入之外，还包括：与转向装置3B的动 作有关的反馈控制输入(进一步具体而言，是指操作通过转向装置3B的动作而产生的前轮W1、W2的横向力的反馈控制输入)。该执行装置动作FB目标值，换言之，是用于操作(修正)作为作用于实车1上的外力的路面反力的反馈控制输入。 

FB分配法则部20基本上是以能使输入的状态量偏差γerr、βerr接近于0的方式来决定假想外力Mvir、Fvir和执行装置动作FB目标值。但是，FB分配法则部20在决定假想外力Mvir、Fvir时，不仅仅只是使状态量偏差γerr、βerr接近于0，还以能抑制实车1或标准动态特性模型16下的车辆规定的限制对象量从规定的容许范围脱离的方式，来决定假想外力Mvir、Fvir。另外，为了绕实车1重心点产生能使状态量偏差γerr、βerr接近于0所需的横摆方向的力矩(更一般而言，以使能让状态量偏差γerr、βerr接近于0的所需外力(路面反力)作用于实车1)，FB分配法则部20决定与驱动/制动装置3A的刹车装置的动作有关的反馈控制输入，或者还决定该反馈控制输入和与转向装置3B的动作有关的反馈控制输入，来作为执行装置动作FB目标值。 

为了决定所述假想外力Mvir、Fvir和执行装置动作FB目标值，不仅将状态量偏差γerr、βerr、还将标准动态特性模型16的输出即标准状态量γd、βd和由传感器/推定器12测定或推定出的实际状态量γact、βact中至少任意一方输入给FB分配法则部20。此外，由传感器/推定器12测定或推定出的实行驶速度Vact、实前轮侧滑角βf_act、实后轮侧滑角βr_act等的实际状态量也被输入给FB分配法则部20。并且，FB分配法则部20依据这些输入，决定假想外力Mvir、Fvir和执行装置动作FB目标值。 

作为补充，假想外力Mvir、Fvir相当于本发明中的车辆模型操作用控制输入，执行装置动作FB目标值相当于本发明中的实车执行装置操作用控制输入。所以，FB分配法则部20具有作为状态量偏差应动控制机构的功能。 

另一方面，在以上说明的标准操作量决定部14、标准动态特性模型16、减法器18以及FB分配法则部20的控制处理的同时(或者根据时分处理)，控制装置10将所述驾驶操作输入输入至FF法则部22，通过该FF法则部22来决定作为执行机构3动作的前馈目标值(基本目标值)的执行装 置动作FF目标值。 

在本实施方式中，执行装置动作FF目标值包括：与由驱动/制动装置3A的刹车装置的动作决定的实车1各车轮W1～W4的驱动/制动力有关的前馈目标值、和与由驱动/制动装置3A的驱动***的动作决定的实车1驱动轮W1、W2的驱动/制动力有关的前馈目标值、和与驱动/制动装置3A的变速装置的减速比(变速比)有关的前馈目标值、和与由转向装置3B的实车1操舵轮W1、W2的舵角有关的前馈目标值。 

为了决定这些执行装置动作FF目标值，所述驾驶操作输入被输入给FF法则部22，并且同时由传感器/推定器12测定或推定出的实际状态量(实行驶速度Vact等)也被输入给FF法则部22。并且，FF法则部22依据这些输入，决定执行装置动作FF目标值。该执行装置动作FF目标值是不依靠所述状态量偏差γerr、βerr(第一状态量偏差)而被决定的、执行机构3的动作目标值。 

作为补充，当悬架装置3C为主动式悬架装置时，执行装置动作FF目标值在一般情况下，还包含与该悬架装置3C的动作有关的前馈目标值。 

接着，控制装置10将由FF法则部22决定的执行装置动作FF目标值(此次值)和由所述FB分配法则部20决定的执行装置动作FB目标值(此次值)输入给该执行装置动作目标值合成部24。然后，控制装置10通过该执行装置动作目标值合成部24合成执行装置动作FF目标值和执行装置动作FB目标值，来决定作为规定执行机构3动作的目标值的执行装置动作目标值。 

在本实施方式中，执行装置动作目标值包括：实车1各车轮W1～W4的驱动/制动力的目标值(由驱动/制动装置3A的驱动***及刹车装置的动作决定的总驱动/制动力的目标值)、实车1各车轮W1～W4的打滑比的目标值、由转向装置3B决定的实车1操舵轮W1、W2的舵角的目标值、由驱动/制动装置3A的驱动***的动作决定的实车1各驱动轮W1、W2的驱动/制动力的目标值、以及驱动/制动装置3A的变速装置的减速比的目标值。 

为了决定这些执行装置动作目标值，不仅将所述执行装置动作FF目标值及执行装置动作FB目标值、还将由传感器/推定器12测定或推定出 的实际状态量(前轮W1、W2的实侧滑角βf_act、推定摩擦系数μestm等)也输入给执行装置动作目标值合成部24。而且，执行装置动作目标值合成部24依据这些输入来决定执行装置动作目标值。 

作为补充，执行装置动作目标值，不仅仅限于上述种类的目标值。可以替换这些目标值，例如来决定对应于该目标值的各执行机构的执行装置操作量的目标值来作为执行装置动作目标值。执行装置动作目标值基本上是只要能规定执行装置的动作即可。例如，作为与刹车装置的动作有关的执行装置动作目标值，可以决定刹车压的目标值，或者决定与之对应的刹车装置的执行装置操作量的目标值。 

接着，控制装置10将由执行装置动作目标值合成部24决定的执行装置动作目标值输入给执行装置驱动控制装置26，并通过该执行装置驱动控制装置26来决定实车1各执行机构3的执行装置操作量。而且，通过所决定的执行装置操作量来控制实车1各执行机构3的执行装置。 

这种情况下，执行装置驱动控制装置26，以满足所输入的执行装置动作目标值、或按照该执行装置动作目标值的方式来决定执行装置操作量。并且，为了该决定，除了执行装置动作目标值以外，由传感器/推定器12测定或推定出的实车1的实际状态量也被输入给执行装置驱动控制装置26。另外，在执行装置驱动控制装置26的控制功能中，与驱动/制动装置3A的刹车装置有关的控制功能最好组装有所谓的防抱死制动***。 

以上是控制装置10的每一控制处理周期的控制处理的概要。 

另外，关于控制装置10的各控制处理功能部的处理，可以适当对这些顺序进行改变。例如，可以在各控制处理周期的最后执行传感器/推定器12的处理，在下次的控制处理周期的处理中使用由此得到的测定值或推定值。 

下面，说明本实施方式中控制装置10的控制处理功能部的更为详细的处理。 

[关于标准动态特性模型] 

首先，参照图3说明本实施方式中的所述标准动态特性模型16。图3是表示本实施方式中的标准动态特性模型16下的车辆构造图。该标准动 态特性模型16是通过下述特性来表示车辆动态特性的模型(所谓2轮模型)，即，该特性是指在前后具有1个前轮Wf和1个后轮Wr的车辆水平面上的动态特性(动力学特性)。以下，将标准动态特性模型16下的车辆(标准动态特性模型16下的对应于实车1的车辆)称为模型车辆。该模型车辆的前轮Wf相当于把实车1的2个前轮W1、W2进行一体化后的车轮，即模型车辆的操舵轮。后轮Wr相当于把实车1的后轮W3、W4进行一体化后的车轮，在本实施方式中属于非操舵轮。 

该模型车辆的重心点Gd在水平面上的速度矢量Vd(模型车辆行驶速度Vd的向量)相对于模型车辆的前后方向所成的角度βd(即，模型车辆的车辆重心点侧滑角βd)和模型车辆绕竖直轴的角速度γd(即，模型车辆的横摆比率γd)是标准状态量，并分别作为所述标准车辆重心点侧滑角、标准横摆比率，通过标准动态特性模型16而逐步决定。另外，模型车辆的前轮Wf的转动面和水平面间的交线相对于模型车辆的前后方向所成的角度δf_d是标准模型操作量，其作为所述模型前轮舵角被输入给标准动态特性模型16。另外，附加性地作用于模型车辆的重心点Gd的横向(模型车辆左右方向)的平移力Fvir、和附加性地绕该模型车辆的重心点Gd而作用的横摆方向的(绕竖直轴的)力矩Mvir是反馈控制输入，并作为所述假想外力被输入给标准动态特性模型16。 

另外，在图3中，Vf_d是模型车辆前轮Wf在水平面上的行进速度矢量，Vr_d是模型车辆后轮Wr在水平面上的行进速度矢量，βf_d是前轮Wf的侧滑角(前轮Wf的行进速度矢量Vf_d相对于前轮Wf的前后方向(前轮Wf的转动面和水平面的交线方向)所成的角度。以下，称为前轮侧滑角βf_d)，βr_d是后轮Wr的侧滑角(后轮Wr的行进速度矢量Vr_d相对于后轮Wr的前后方向(后轮Wr的转动面和水平面的交线方向)所成的角度。以下，称为后轮侧滑角βr_d)，βf0是模型车辆的前轮Wf的行进速度矢量Vf_d相对于模型车辆前后方向所成的角度(以下，称为车辆前轮位置侧滑角βf0)。 

作为补充，在本说明书的实施方式中，从车辆上方来看，关于车辆或者车轮的侧滑角、车轮的舵角、车辆的横摆比率、横摆方向的力矩，以逆时针方向为正方向。另外，假想外力Mvir、Fvir中的平移力Fvir，是以车 辆的左向朝向为正的方向。同时，车轮的驱动/制动力是以在车轮的转动面与路面或水平面的交线方向上、使车辆朝向行进方向加速的力(路面反力)的方向为正的方向。换言之，当相对于车辆的行进方向是驱动力时的朝向时，驱动/制动力为正值；当相对于车辆的行进方向是制动力时的朝向时，驱动/制动力为负值。 

该模型车辆的动态特性(连续性的动态特性)具体而言可由下式01表述。另外，将式子01右边的“K”看作单位矩阵，并且除去右边括号内第三项后(包含Fvir、Mvir的项)的式子与例如题名为“汽车的动作与控制”的公知文献(作者：安部正人，发行社：株式会社山海堂，2004年7月23日第二版第二次印刷发行。以后，称为非专利文件1)记载的公知的式子(3.12)、(3.13)相同。 

【数1】 

$\frac{d}{\mathrm{dt}}\left[\begin{array}{c}\mathrm{\βd}\\ \mathrm{\γd}\end{array}\right]=K\·(\left[\begin{array}{cc}a11& a12\\ a21& a22\end{array}\right]\·\left[\begin{array}{c}\mathrm{\βd}\\ \mathrm{\γd}\end{array}\right]+\left[\begin{array}{c}b1\\ b2\end{array}\right]\·\mathrm{\δf}\_d+\left[\begin{array}{cc}b11& 0\\ 0& b22\end{array}\right]\·\left[\begin{array}{c}\mathrm{Fvir}\\ \mathrm{Mvir}\end{array}\right])$ ……式子01 

其中，  $K=\left[\begin{array}{cc}k1& 0\\ 0& k2\end{array}\right]$

$a11=-\frac{2\·(\mathrm{Kf}+\mathrm{Kr})}{m\·\mathrm{Vd}}$              $a12=-\frac{m\·V{d}^2+2\·(\mathrm{Lf}\·\mathrm{Kf}-\mathrm{Lr}\·\mathrm{Kr})}{m\·V{d}^2}$

$a21=-\frac{2\·(\mathrm{Lf}\·\mathrm{Kf}-\mathrm{Lr}\·\mathrm{Kr})}{I}$          $a22=-\frac{2\·(L{f}^2\·\mathrm{Kf}+L{r}^2\·\mathrm{Kr})}{I\·\mathrm{Vd}}$

$b1=\frac{2\·\mathrm{Kf}}{m\·\mathrm{Vd}}$       $b2=\frac{2\·\mathrm{Lf}\·\mathrm{Kf}}{I}$       $b11=\frac{1}{m\·\mathrm{Vd}}$        $b22=\frac{1}{I}$

在该式子01的条件部分，m为模型车辆的总质量，Kf为将模型车辆前轮Wf视为2个左右前轮的连结体时的每1车轮的回转率(corneringpower)，Kr是将模型车辆后轮Wr视为2个左右后轮的连结体时的每1车轮的回转率，Lf是模型车辆的前轮Wf中心与车辆重心点Gd的前后方向上的距离(前轮Wf的舵角为0时的该前轮Wf的转轴与车辆重心点Gd在前后方向上的距离。参照图3)，Lr是模型车辆的后轮Wr的中心与车辆重心点Gd的前后方向上的距离(后轮Wr的转轴与车辆重心点Gd在前后 方向上的距离。参照图3)，I是模型车辆的重心点Gd的绕横摆轴的惯性(惯性力矩)。这些参数m、I、Lf、Lr、Kf、Kr的数值(或是式子01右边的参数a11、a12、a21、a22、b1、b2、b11、b22的数值)是被预先设定了的数值。此时，例如将m、I、Lf、Lr设定成与实车1中的这些值相同或大致相同。同时，分别考虑了实车1的前轮W1、W2、后轮W3、W4的轮胎特性(或者该轮胎所被要求的特性)来设定Kf、Kr。例如分别与实车1的前轮W1、W2、后轮W3、W4的轮胎回转率相同或大致相同的形式来设定Kf、Kr。 

更一般而言，以使稳定状态(转向角θh等的所述驾驶操作输入为恒定，并且路面状态等的环境状态为恒定或一致的情况下、在过渡性的动作消失为止经过了充足的时间的状态)中的模型车辆的动作特性(模型车辆的稳定特性)接近于实车1的稳定状态时的动作特性(实车1的稳定特性)的形式来设定m、I、Lf、Lr、Kf、Kr的数值(或是a11、a12、a21、a22、b1、b2、b11、b22的数值)。在该情况下，本实施方式中的模型车辆的稳定特性由转向角θh与模型车辆的车辆重心点侧滑角βd以及横摆比率γd之间的关系(其中设Vd＝Vact)来表示，其中，上述关系是在该模型车辆的动作状态量βd、γd的时间性变化(式01的左边)为0、且假想外力Fvir、Mvir为0时的状态。所以，以使该关系与稳定状态中的实车1的转向角θh与实际车辆的重心点侧滑角βact以及实际横摆比率γact之间的关系(执行装置动作FB目标值设为0时的关系)大致相一致的形式来设定模型车辆的m、I、Lf、Lr、Kf、Kr的数值(或是a11、a12、a21、a22、b1、b2、b11、b22的数值)。 

另外，式子01中的矩阵式K是用于调整特性的矩阵式(对角矩阵式)，其为了调整模型车辆的动态特性(更详细而言是指模型车辆的状态量βd、γd的过渡响应特性)，而分别在式子01右边括号内的各项的第一行、第二行的值上乘上了K1倍和K2倍(其中，K1≠0、K2≠0)。在以后的说明中，将特性调整用矩阵式K的各对角分量K1、K2称为模型特性调整用参数。 

在本实施方式中，为了调整模型车辆的状态量βd、γd的过渡响应特性，在标准动态特性模型16的处理中，将特性调整用矩阵式K的模型特 性调整用参数k1、k2中的k2的值设定为可变的数值。并且，k1的值被固定设为1。因此，在本实施方式中，通过将模型特性调整用参数k2的值设为可变的数值，来调整模型车辆的状态量βd、γd的过渡响应特性，其中，模型特性调整用参数k2是将式子01右边括号内的各项的第二行的值乘以k2倍之后得到的值。另外，模型特性调整用参数k2的值的设定方法在后面进行叙述。 

作为补充，由于特性调整用矩阵式K是使式子01右边括号内的同行的所有项成为等倍数的矩阵式，只要特性调整用矩阵式K的各对角分量k1、k2的值不为0，模型车辆的稳定特性将不受k1、k2值的影响而维持恒定。所以，特性调整用矩阵式K(模型特性调整用参数k1、k2)是将模型车辆的稳定特性维持在恒定的同时并调整过渡响应特性的矩阵式。 

在本实施方式的标准动态特性模型16的处理中，将前式01的δf_d、Fvir、Mvir作为输入值，通过以控制装置10的控制处理周期逐次执行该式子01的运算处理(详细而言，将式子01用离散时间***来表示后的式子的运算处理)，由此时序地逐次计算出βd、γd。此时，在各控制处理周期中，由所述传感器/推定器12检测或推定的实行驶速度Vact的最新值(此次值)被用来作为模型车辆的行驶速度Vd的值。即，模型车辆的行驶速度Vd被设为总是与实行驶速度Vact一致。并且，根据如上述方式被设定为与实行驶速度Vact相一致的模型车辆的行驶速度Vd(换言之，根据实行驶速度Vact)，来设定所述模型特性调整用参数k1、k2中的k2值，使其成为可变动的值。另外，作为Fvir、Mvir的值，使用由FB分配法则部20按照后述的方式决定的假想外力的最新值(前次值)。另外，作为δf_d的值，使用由标准操作量决定部14按照后述的方式决定的模型前轮舵角的最新值(此次值)。此外，为了计算新的βd、γd的值(此次值)，还使用βd、γd的前次值。 

作为补充，由下式02a、02b、02c来表示模型车辆的βf0、βd、βf_d、βr_d、γd、δf_d之间的关系。

βf_d＝βd+Lf·γd/Vd-δf_d        ......式02a 

βr_d＝βd-Lr·γd/Vd              ......式02b 

βf0＝βf_d+δf_d＝βd+Lf·γd/Vd  ......式02c 

另外，如图3所示，将模型车辆的前轮Wf的回转力(≈前轮Wf的横向力)设为Ffy_d，将模型车辆的后轮Wr的回转力(≈后轮Wr的横向力)设为Fry_d，并由下式03a、03b来表示Ffy_d与βf_d的关系以及Fry_d与βr_d的关系。 

Ffy_d＝-2·Kf·βf_d  ......式03a 

Fry_d＝-2·Kr·βr_d  ......式03b 

以下对本实施方式中的模型特性调整用参数k2的值的设定方法进行说明。 

在实车1的动态特性中(将执行装置动作FB目标值维持在0的状态时的动态特性)中，一般，当行驶速度Vact在某个数值Vd_critical以下时，使驾驶操作输入中的转向角θh成阶梯状变化时(进而使实前轮舵角δf_act呈阶梯状变化时)的实际状态量βact、γact的响应特性(过渡响应特性)为非振动性的响应特性(过度制动或临界制动的响应特性)。即，实际状态量βact、γact相对于转向角θh的阶梯型变化(相对于实前轮舵角δf_act的阶梯型变化)，不含带以最终的恒定值为中心的振动分量而趋向于收敛至该恒定值。另外，当Vact＝Vd_critical时，实际状态量βact、γact的响应特性成为临界制动的特性。在以后的说明中，称Vd_critical为临界制动速度。 

另外，当实行驶速度Vact一旦大于临界制动速度Vd_critical时，使转向角θh呈阶梯状变化时的实际状态量βact、γact的响应特性为含带振动分量的响应特性(振动性的特性)。即，如图4(a)、(b)中分别用虚线a、b例示的那样，实际状态量βact、γact相对于转向角θh的阶梯型变化，在产生以最终的恒定值为中心的振动的同时，逐渐收敛于该恒定值。并且， 该振动分量的衰减速度随着Vact的上升而变慢(振动分量的振幅衰减的时间常数变长)。所以，实行驶速度Vact一旦超过临界制动速度Vd_critical，实际状态量βact、γact的衰减性随着Vact的上升而降低。 

另一方面，如前所述，在本实施方式中，经执行机构3控制实车1的动作以使所述状态量偏差βerr、γerr接近于0。因此，模型车辆的状态量βd、γd一旦产生振动，实车1的状态量βact、γact也容易产生振动。所以，当使操舵角θ产生变化时，实车1的状态量βact、γact会有产生振动的情况，但为了使该情况尽可能不发生，而使实车1的动作特性良好，最好是使模型车辆的状态量βd、γd的响应特性成为非振动的特性、或者使状态量βd、γd的振动分量的衰减速度尽可能地快速。更一般而言，最好是采用下述情况，即：相对转向角θh的阶梯状变化，模型车辆的状态量βd、γd的响应特性中的衰减性至少要比实车1的状态量βact、γact的响应特性中的衰减性要高。 

因此，在本实施方式中，通过根据行驶速度Vact(＝Vd)而将所述模型特性调整用参数k2的值设定为可变，使下述情况时的模型车辆的状态量βd、γd的响应特性(详细而言，将式子01的假想外力Fvir、Mvir恒定设为0时的响应特性)成为非振动(成为临界制动或过度制动的特性)的形式，所述情况是指：针对任意的行驶速度Vact，使转向角θh成阶梯状变化时的情况。并且，由此使相对于转向角θh变化的模型车辆的状态量βd、γd的衰减性比实车1的状态量βact、γact的衰减性要高。 

在本实施方式中，具体为：对于Vd≤Vd_critical(Vact≤Vd_critical)的任意的行驶速度Vd(＝Vact)，模型特性调整用参数k1的值被设定为1。同时，对于Vd>Vd_critical(Vact>Vd_critical)的任意的行驶速度Vd(＝Vact)，以使模型车辆(由包含特性调整用矩阵式K的式子01表述的***)的状态量βd、γd的响应特性成为非振动性的特性、例如临界制动的特性的形式，设定模型特性调整用参数k2的值。即，对于Vd>Vd_critical的Vd(＝Vact)的任意值，由包含了特性调整用矩阵式K(其中，在本实施方式中k1＝1)的式子01所表述的***(模型车辆)的特性方程式det(λ·I—K·A)＝0(其中，λ：标量变量；det()：矩阵式；I：单位矩阵)的解、即模型车辆的固有值成为重根的形式根据Vd(＝Vact)来设定k2的值。另外，“A”是式子01右边的括号内第一项的矩阵式(以a11、a12、a21、a22为组成部分的2行2列的矩阵式)。 

针对Vd＞Vd_critical的Vd(＝Vact)的任意值，如上被设定的k2的值是满足了下式101的k2的值。 

(a11+k2·a22)²-4·k2·(a11·a22+a12·a21)＝0···式子101 

作为补充，临界制动速度Vd_critical是接近于下述Vd值的一种数值，该Vd值是诸如使式子01右边括号内第一项的矩阵式A的固有值会成为重根的数值。即，临界制动速度Vd_critical成为接近于下述Vd值的一种数值，该Vd值是指：在将式子01的特性调整用矩阵式K看作单位矩阵式，并且由除去了包含假想外力Fvir、Mvir的项的式子(该式子表示使实车1的动态特性近似为线性的一种式子)所表述的***的特性方程式det(λ·I-A)＝0(其中，λ：标量变量；det()：矩阵式；I：单位矩阵)的解、即该***的固有值会成为重根的Vd的值。所以，Vd_critical被设定为满足了下式102的Vd的值(更具体而言，是满足了下式102的Vd的值中正的值)。 

(a11+a22)²-4·(a11·a22+a12·a21)＝0···式子102 

图5是表示根据如上所述的行驶速度Vact(＝Vd)而设定的k2的值变化的图表。如图所示，k2的值在Vact＞Vd_critical的情况下，被设定为大于1的数值。此时，伴随Vact(＝Vd)的上升，k2的值呈单调递增。 

如上所示，通过将模型特性调整用参数k的值设定为根据行驶速度Vact(＝Vd)而可变，能够不受行驶速度Vact的影响而将模型车辆的状态量γd、βd的响应特性设定为非振动性的特性，其中，所述模型车辆的状态量γd、βd对应了转向角θh的阶梯状的变化。例如，如图4(a)、(b)中实线c、d所例示的那样，比临界制动速度Vd_critical高的行驶速度Vact时的模型车辆的状态量γd、βd的响应特性会成为临界制动的特性(非制动性的特性)。另外，在图4(a)、(b)中，实车1的车轮与路面之间的摩擦系数 μ与模型车辆的相一致的情况下，以使实车1的稳定特性与模型车辆的相一致的形式来设定模型车辆的动作特性。此外，临界制动速度Vd_critical以下的行驶速度Vact时的模型车辆的状态量γd、βd的响应特性成为临界制动、或过度制动的一种特性。 

所以，在大于等于Vd_critical的速度范围，能够使模型车辆的状态量γd、βd各自的衰减性比实车1的状态量γact、βact各自的衰减性要高，其中，模型车辆的状态量γd、βd各自的衰减性对应了转向角θh的阶梯状的变化。 

而且，如上所述的那样，由于模型车辆的稳定特性不依存于k2的值，因此能够将模型车辆的稳定特性维持在与实车1的稳定特性相近的一种特性上。所以，能够在将模型车辆的稳定特性维持在与实车1的稳定特性相近的特性上的同时，能够防止对应了转向角θh的阶梯状的变化的模型车辆状态量γd、βd的响应特性成为振动性的一种特性。 

如上述说明，在本实施方式中，在控制装置10的各控制处理周期中的标准动态特性模型16的处理中，模型特性调整用参数k2的值被设定为对应行驶速度Vact(＝Vd)可变。该k2的值的设定处理相当于本发明中的车辆模型特性设定机构。另外，在标准动态特性模型16的处理中，实际设定k2的值时，只要事先确定表示k2的值与实车1的行驶速度Vact的关系(图5所示的关系)的图表或运算式，使用各控制处理周期中的Vact的值(此次值)，基于该图表或运算式来决定k2的值即可。 

作为补充，在Vact>Vd_critical的情况下的k2值不一定要设成使模型车辆的状态量γd、βd的响应特性成为临界制动的特性的值，也可以将k2值设定成该响应特性成为过度制动的特性的值。即，在Vact>Vd_critical的情况下的k2值可以设定成下述值：使所述特性方程式det(λ·I—K·A)＝0的解(模型车辆的固有值)成为两个不同数值的实数的值。这种情况下，k2的值被设定成大于比如使模型车辆的状态量γd、βd的响应特性成为临界制动的特性时的K2的值。但是，在Vact>Vd_critical时的实车1的动作状态中，在避免模型车辆的状态量γd、βd的响应特性(衰减性)较大地乖离于实际状态量γact、βact的响应特性这一方面考虑，k2的值最好是一种接近于使模型车辆的状态量γd、βd的响应特性成为临界制动的 特性的k2的值。 

另外，可以将Vact>Vd_critical的实车1的动作状态下的k2的值设定为：在下述范围内，比使模型车辆的状态量γd、βd的响应特性成为临界制动的特性的k2的值要稍微接近于1的值，其中所述范围是指使模型车辆的状态量γd、βd的衰减性比实际状态量γact、βact的衰减性要高的范围。通过这种做法，能够防止模型车辆的状态量γd、βd的响应特性(衰减性)较大地乖离于实际状态量γact、βact的响应特性。但是，这种情况下的模型车辆的状态量γd、βd的响应特性中，状态量γd、βd的衰减性比k2＝1的情况时的高，但比临界制动的特性要低。 

此外，在由式子01表述的模型车辆中，当行驶速度Vact(＝Vd)是比临界制动速度Vd_critical要低得多的极低速的情况时，相对于转向角θh的变化(前轮舵角δd的变化)的模型车辆的状态量γd、βd的响应性(快速响应性)与实车1相比有变高的倾向。因此，在行驶速度Vact为极低速(在比临界制动速度Vd_critical还要小的规定值以下的行驶速度)的情况下，可以将k2的值设为比1小的值。但是，在该情况下，将极低速时的k2的值设定为大于等于使模型车辆的状态量γd、βd的响应特性成为临界制动的特性时的k2值的值，从而使模型车辆的状态量γd、βd的响应特性不会成为一种振动性的特性。以该种方式设定极低速时的k2的值，在极低速时，可以防止模型车辆的状态量γd、βd的响应特性成为振动性的特性，并同时能够使相对于转向角θh的变化的模型车辆的状态量γd、βd的响应性(γd、βd相对于转向角θh的变化的启动速度)接近于实际状态量γact、βact的响应性。 

但是，式子01的括号内的各项的第二行乘上k2倍的意思相当于在模型车辆的惯量I上乘以(1/k2)倍。所以，作为使用特性调整用矩阵式K的替代，还可以使用由I’＝I/k2定义的惯量I’并通过式子103来表述模型车辆的动态特性。 

【数2】

$\frac{d}{\mathrm{dt}}\left[\begin{array}{c}\mathrm{\βd}\\ \mathrm{\γd}\end{array}\right]=K\·(\left[\begin{array}{cc}a11& a12\\ a21& a22\end{array}\right]\·\left[\begin{array}{c}\mathrm{\βd}\\ \mathrm{\γd}\end{array}\right]+\left[\begin{array}{c}b1\\ b2\end{array}\right]\·\mathrm{\δf}\_d+\left[\begin{array}{cc}b11& 0\\ 0& b22\end{array}\right]\·\left[\begin{array}{c}\mathrm{Fvir}\\ \mathrm{Mvir}\end{array}\right])$  ......式子103 

其中， 

$a11=-\frac{2\·(\mathrm{Kf}+\mathrm{Kr})}{m\·\mathrm{Vd}}$               $a12=-\frac{m\·V{d}^2+2\·(\mathrm{Lf}\·\mathrm{Kf}-\mathrm{Lr}\·\mathrm{Kr})}{m\·V{d}^2}$

$a21=-\frac{2\·(\mathrm{Lf}\·\mathrm{Kf}-\mathrm{Lr}\·\mathrm{Kr})}{I}$           $a22=-\frac{2\·(L{f}^2\·\mathrm{Kf}+L{r}^2\·\mathrm{Kr})}{I\·\mathrm{Vd}}$

$b1=\frac{2\·\mathrm{Kf}}{m\·\mathrm{Vd}}$        $b2=\frac{2\·\mathrm{Lf}\·\mathrm{Kf}}{I}$       $b11=\frac{1}{m\·\mathrm{Vd}}$         $b22=\frac{1}{I}$

$I\′=\frac{I}{k2}$

在该情况时，如上所述只需根据行驶速度Vact(＝Vd)来设定k2的值即可。这时，在Vact>Vd_critical时的情况下，惯量I’伴随Vact的上升而呈单调递减。 

另外，本实施方式中的实车1的动态特性表示的是下述一种中间性的特性，即：不适用本发明的实车1的开放特性(将所述执行装置FB动作目标值恒定维持在0时的实车1的动态特性)和将假想外力Mvir、Fvir恒定维持在0时的标准动态特性模型16的动态特性之间的一种特性。因此，一般最好将标准动态特性模型16设定成一种表示了与实车1的开放特性相比驾驶者认为更理想的动态特性的模型。具体而言，最好将标准动态特性模型16设定为与实车1相比为高线性的模型。例如，以使模型车辆的车轮的侧滑角或打滑比与从路面作用于该车轮的路面反力(横向力或是驱动/制动力)之间的关系成为线性关系或与线性关系相近的关系的形式来设定标准动态特性模型16为佳。通过所述式子01表述了动态特性的标准动态特性模型16即为满足这些要求的模型的其中一个例子。 

[关于标准操作量决定部] 

接着，参照图6及图7说明所述标准操作量决定部14的详细处理。 图6是表示所述标准操作量决定部14的具体的处理功能的功能方框图，图7是用于说明标准操作量决定部14所具有的防止离心力过大的限幅器14f的处理的图表。 

参照图6，首先，在处理部14a，标准操作量决定部14通过所输入的驾驶操作输入中的转向角θh(此次值)除以总转向比is来决定无限制时前轮舵角δf_unltd。该无限制时前轮舵角δf_unltd具有：作为与转向角θh对应的模型前轮舵角δf_d的基本要求值的意思。 

在此，总转向比is是转向角θh与模型车辆前轮Wf的操舵角的比率，例如，是对照实车1的转向角θh和与该转向角θh对应的实车1前轮W1、W2的操舵角的前馈值之间的关系而被设定的。 

另外，也可以不使总转向比is为一定值(固定值)，而根据由传感器/推定器12测定或推定出的实车1的行驶速度Vact，可变性地进行设定。这种情况下，最好是总转向比is随着实车1的行驶速度Vact的变大而变大的方式来设定is。 

接着，通过βf0计算部14b求解出标准动态特性模型16下的模型车辆的车辆前轮位置侧滑角βf0。由标准动态特性模型16决定的标准横摆比率γd和标准车辆重心点侧滑角βd的前次值被输入给该βf0计算部14b，并根据这些值，通过上述式子02c的运算(式子02c的第二个等号右边的运算)来求解出βf0的前次值。另外，此时，使用实行驶速度Vact的前次值作为式子02c运算中所需的Vd值。因此，由βf0计算部14b计算出的βf0即为前次的控制处理周期中的模型车辆的车辆前轮位置侧滑角βf0的值。 

另外，也可以根据γd、βd的前次值、由标准操作量决定部14决定的模型前轮舵角δf_d的前次值以及实行驶速度Vact的前次值，通过上述式子02a的运算来求解出模型车辆的前轮侧滑角βf_d的前次值，并将由标准操作量决定部14决定的模型前轮舵角δf_d的前次值加上所求得的βf_d(进行式子02c的第一个等号右边的运算)，从而来求解出βf0。另外，在各控制处理周期，也可以通过标准动态特性模型16的处理来执行βf0的计算，将所计算出的βf0的前次值输入给标准操作量决定部14。这种情况下，不需要标准操作量决定部14的βf0计算部14b的运算处理。 

此后，通过减法器14c，从上述所求得的车辆前轮位置侧滑角βf0减 去无限制时前轮舵角δf_unltd，由此来求得无限制时前轮侧滑角。该无限制时前轮侧滑角，是表示将模型车辆的模型前轮舵角δf_d从前次值瞬间地控制在无限制时前轮舵角δf_unltd(此次值)时所产生的模型车辆的前轮侧滑角βf_d的瞬间预测值。 

之后，标准操作量决定部14使该无限制时前轮侧滑角通过前轮侧滑角限幅器14d，来决定限制完毕前轮侧滑角。在此，图中所示的前轮侧滑角限幅器14d的坐标是例示无限制时前轮侧滑角与限制完毕前轮侧滑角之间关系的坐标，与该坐标相关的横轴方向的值为无限制时前轮侧滑角的值，纵轴方向的值为限制完毕前轮侧滑角的值。 

该前轮侧滑角限幅器14d是用于抑制模型车辆的前轮侧滑角βf_d变得过大(进而使得对实车1所要求的前轮W1、W2的横向力不会变得过大)的限幅器。 

在本实施方式中，前轮侧滑角限幅器14d根据从传感器/推定器12输入给标准操作量决定部14的推定摩擦系数μestm(此次值)及实行驶速度Vact(此次值)来设定前轮侧滑角βf_d的容许范围(具体而言，是指该容许范围的上限值βf_max(>0)及下限值βf_min(<0))。这种情况下，一般是推定摩擦系数μestm越小、或者实行驶速度Vact越高，使得容许范围[βf_min、βf_max]越窄(使βf_max、βf_min接近于0)的方式来设定该容许范围。此时，该容许范围[βf_min、βf_max]例如设定在：实车1的前轮W1、W2的侧滑角与横向力或回转力之间的关系大致被维持在线性关系(正比关系)下的侧滑角的值的范围内。 

另外，该容许范围[βf_min、βf_max]可以根据μestm和Vact中任意一方来设定，或者也可以不根据μestm和Vact而是设定在预先固定的容许范围内。 

而且，在所输入的无限制时前轮侧滑角是在如上所述那样设定的容许范围[βf_min、βf_max]内的值时(βf_min≤无限制时前轮侧滑角≤βf_max之时)，前轮侧滑角限幅器14d直接将无限制时前轮侧滑角的值作为限制完毕前轮侧滑角来输出。另外，在所输入的无限制时前轮侧滑角的值脱离出容许范围时，该前轮侧滑角限幅器14d则将容许范围[βf_min、βf_max]的下限值βf_min或上限值βf_max作为限制完毕前轮侧滑角来输出。具体而言， 当无限制时前轮侧滑角>βf_max时，βf_max作为限制完毕前轮侧滑角被输出，当无限制时前轮侧滑角<βf_min时，βf_min作为限制完毕前轮侧滑角被输出。据此，限制完毕前轮侧滑角在容许范围[βf_min、βf_max]内被决定为与无限制时前轮侧滑角一致、或者成为最接近于该无限制时前轮侧滑角的值。 

而后，通过利用减法器14e从由所述βf0计算部14b求得的车辆前轮位置侧滑角βf0减去如上所述求得的限制完毕前轮侧滑角，求解出第一限制完毕前轮舵角δf_ltd1。如此求得的第一限制完毕前轮舵角δf_ltd1具有作为下述模型前轮舵角δf_d的意思，即该模型前轮舵角δf_d是以使模型车辆的前轮侧滑角βf_d不会脱离容许范围[βf_min、βf_max]而对无限制时前轮舵角δf_unltd加以限制而成的角度。 

然后，标准操作量决定部14通过使该第一限制完毕前轮舵角δf_1td1通过离心力过大防止限幅器14f来决定第二限制完毕前轮舵角δf_ltd2。该δf_ltd2是作为输入给标准动态特性模型16的模型前轮舵角δf_d的值被使用。在此，图中所示的离心力过大防止限幅器14f的坐标是例示第一限制完毕前轮舵角δf_ltd1和第二限制完毕前轮舵角δf_ltd2之间关系的坐标，与该坐标相关的横轴方向的值为δf_ltd1值，纵轴方向的值为δf_ltd2值。 

该离心力过大防止限幅器14f是用于使得模型车辆产生的离心力不至于变得过大(进而使得对实车1所要求的离心力不会变得过大)的限幅器。 

在本实施方式中，离心力过大防止限幅器14f根据输入给标准操作量决定部14的推定摩擦系数μestm(此次值)与实行驶速度Vact(此次值)来设定模型前轮舵角δf_d的容许范围(具体而言，是指该容许范围的上限值δf_max(>0)及下限值δf_min(<0))。该容许范围[δf_min、δf_max]，是在假想外力Mvir、Fvir一直被保持在0时，模型车辆不会超越与路面间的摩擦界限即可进行正常转圈的模型前轮舵角δf_d的容许范围。 

具体而言，首先，依据输入给标准操作量决定部14的Vact、μestm的值(此次值)求解出满足下式05的横摆比率即正常转圈时的最大横摆比率γmax(>0)。

m·γmax·Vact＝C1·μestm·m·g    ……式05 

在此，式子05中的m如前所述，是模型车辆的总质量。同时，g是重力加速度，C1是1以下的正系数。该式子05左边表示离心力，该离心力是在将模型车辆的横摆比率γd和行驶速度Vd分别保持在γmax、Vact进行该模型车辆正常转圈时而在该模型车辆上产生的离心力(进一步具体而言，是该离心力的收敛预想值)。另外，式子05右边演算结果的值是根据μestm而确定的路面反力(具体而言，是可从路面经车轮Wf、Wr而作用于模型车辆的总摩擦力(路面反力的平移力水平分量的总和))的大小的界限值乘上系数C1之后的值(≤界限值)。因此，正常转圈时最大横摆比率γmax被决定为：在下述情况下在模型车辆上产生的离心力不会超过对应于推定摩擦系数μestm而可作用于模型车辆的总摩擦力(路面反力的平移力水平分量的总和)的界限值，即该情况为，将作用于模型车辆上的假想外力Mvir、Fvir保持在0的同时，将模型车辆的横摆比率γd及行驶速度Vd分别保持在γmax、Vact，并进行该模型车辆的正常转圈。 

另外，对应于μestm、Vact中的至少任意一方的值，式子05的系数C1值可以进行变动性设定。此时，最好是μestm越小、或Vact越高，将C1的值设得越小。 

接着，模型车辆的正常转圈时与γmax对应的模型前轮舵角δf_d的值是作为正常转圈时限界舵角δf_max_c(>0)而被求得的。在此，在由上述式子01表示的标准动态特性模型16方面，在正常转圈时的模型车辆的横摆比率γd和模型前轮舵角δf_d之间，下式06的关系成立。 

【数3】 

$\mathrm{\&gamma;d}=\frac{1}{1-\frac{m}{2\&CenterDot;L^2}\&CenterDot;\frac{\mathrm{Lf}\&CenterDot;\mathrm{Kf}-\mathrm{Lr}\&CenterDot;\mathrm{Kr}}{\mathrm{Kf}\&CenterDot;\mathrm{Kr}}\&CenterDot;V{d}^2}\&CenterDot;\frac{\mathrm{Vd}}{L}\&CenterDot;\mathrm{\&delta;}\&CenterDot;f\_d$      ......式06 

其中、L＝Lf+Lr 

另外，当Vd足够小时，式子06可以近似地转换为下式07。

γd＝(Vd/L)·δf_d     ......式07 

因此，在本实施方式中，将式子06或式子07中的γd、Vd各值作为γmax、Vact，对δf_d进行求解，由此求出与γmax对应的正常转圈时限界舵角δf_max_c。 

用于使得产生于模型车辆的离心力不至于变得过大的模型前轮舵角δf_d的容许范围[δf_min、δf_max]，一般只要设定在容许范围[—δf_min_c，δf_max_c]即可。但是，此时，在实车1的逆向操舵状态(朝着与实车1的横摆比率的极性相反的极性方向来操舵前轮W1、W2的状态)，有时模型前轮舵角δf_d受到不必要的限制。 

因此，在本实施方式中，根据模型车辆的横摆比率γd和γmax，通过下式08a、08b对δf_max_c、—δf_max_c进行修正，由此设定模型前轮舵角δf_d的容许范围的上限值δf_max及下限值δf_min。 

δf_max＝δf_max_c+fe(γd、γmax)  ……式08a 

δf_min＝—δf_max_c—fe(—γd、—γmax)  ……式08b 

式子08a、08b中的fe(γd、γmax)、fe(—γd、—γmax)是γd、γmax的函数，其函数值例如如图7(a)、(b)的坐标所示，是根据γd、γmax的值而变化的函数。在该例子中，函数fe(γd、γmax)的值如图7(a)的坐标所示，当γd为比0稍微大些的规定值γ1以下的值时(包含γd<0的情况)，则为正的一定值fex。同时，fe(γd、γmax)的值在γd>γ1时，随着γd变大而呈单调减少，在γd达到γmax以下的规定值即γ2(>γ1)时，则变为0。而且，fe(γd、γmax)的值在γd>γ2时(包含γd≥γmax的情况在内)，则被维持在0。 

另外，因为函数fe(—γd、—γmax)是使函数fe(γd、γmax)的变量γd、γmax的极性反转了的函数，因此该函数fe(—γd、—γmax)的值相对于图7(b)的坐标所示的γd而变化。即，在γd为比0略小的规定负值—γ1以上的值时(包含γd>0的情况在内)，为正的一定值fex。而且，fe(—γd、—γmax)的值在γd<—γ1时，随着γd变小而单调减少，在γd达到—γmax以上的规定值即-γ2时，则变为0。此外，fe(-γd、-γmax)的值在γd＜-γ2时(包含γd≤-γmax时)，被维持在0。 

另外，作为决定函数fe(γd、γmax)、fe(-γd、-γmax)的值所需的γd值，使用由标准动态特性模型16决定的标准横摆比率γd的前次值即可。 

另外，函数fe(γd、γmax)在坐标中转折点处的γd的值γ1、γ2、或者所述正的一定值fex，也可以根据推定摩擦系数μestm、实行驶速度Vact而进行可变性更改。 

如上所述，通过函数fe的值，对δf_max_c进行补正，设定模型前轮舵角δf_d的容许范围[δf_min、δf_max]，由此，与γd朝向相反方向的模型前轮舵角δf_d的界限值δf_max或δf_min的大小(绝对值)被设定为比正常转圈时限界舵角δf_max_c较大，其中该正常转圈时限界舵角δf_max_c是指与产生在模型车辆上的离心力的界限对应的操舵角。由此，在实车1的逆向操舵状态中，可以防止模型前轮舵角δf_d受到不必要的限制。另外，实行驶速度Vact越高、或者推定摩擦系数μestm越小，则该容许范围[-δf_min、δf_max]为越窄。 

如上所述，在设定模型前轮舵角δf_d的容许范围后，当被输入的第一限制完毕前轮舵角δf_ltd1为容许范围[δf_min、δf_max]内的值时(δf_min≤δf_ltd1≤δf_max时)，离心力过大防止限幅器14f将δf_ltd1的值直接作为第二限制完毕前轮舵角δf_ltd2(＝输入给标准动态特性模型16的模型前轮舵角δf_d)来输出。同时，当被输入的δf_ltd1值脱离出容许范围[δf_min、δf_max]时，该离心力过大防止限幅器14f则将强制限制了其输入值之后的值作为第二限制完毕前轮舵角δf_ltd2来进行输出。具体而言，当δf_ltd1＞δf_max时，δf_max作为第二限制完毕前轮舵角δf_ltd2被输出，当δf_1td1＜δf_min时，δf_min作为第二限制完毕前轮舵角δf_ltd2被输出。由此，δf_ltd2在容许范围[δf_min、δf_max]内，被决定为与第一限制完毕前轮舵角δf_ltd1一致、或是成为最接近于第一限制完毕前轮舵角δf_ltd1的值。 

另外，在上述式子01所表示的标准动态特性模型16中，在模型车辆的正常转圈时，βd与γd之间下式09的关系成立。 

【数4】 

$\mathrm{\&beta;d}=(1-\frac{m}{2\&CenterDot;L}\&CenterDot;\frac{\mathrm{Lf}}{\mathrm{Lr}\&CenterDot;\mathrm{Kr}}\&CenterDot;V{d}^2)\&CenterDot;\frac{\mathrm{Lr}}{\mathrm{Vd}}\&CenterDot;\mathrm{\&gamma;d}$      ......式09 

另外，当Vd足够小时，式子09可以近似地转换为下式10。 

βd＝(Lr/Vd)·γd     ……式10 

因此，模型车辆的正常转圈时的γd或者γmax的值，可以通过式子09或式子10转换为βd的值(其中，Vd＝Vact)。因此，可以替代如上所述根据横摆比率γd、γmax的值来设定模型前轮舵角δf_d的容许范围的方式，而采用根据与横摆比率γd、γmax对应的车辆重心点侧滑角βd的值来设定模型前轮舵角δf_d的容许范围。 

以上是标准操作量决定部14的详细处理。 

通过以上说明的标准操作量决定部14的处理，标准动态特性模型16下的模型车辆的前轮侧滑角βf_d的瞬间值不会变得过大，并且使得产生在模型车辆上的离心力不会变得过大，同时根据驾驶操作输入中的转向角θh，在每一控制处理周期，将第二限制完毕前轮舵角δf_ltd2决定为输入给标准动态特性模型16的模型前轮舵角δf_d。 

作为补充，在离心力过大防止限幅器14f中，如上所述对输入给标准动态特性模型16的模型前轮舵角δf_d进行限制，使产生在模型车辆上的离心力不会变得过大，这与使得模型车辆的车辆重心点侧滑角βd(或后轮侧滑角βr_d)不会变得过大的方式来限制模型前轮舵角δf_d的说法是一样的。另外，一般情况下，由于车辆的离心力、车辆重心点侧滑角(或者后轮侧滑角)相对于转向操作会产生延迟，因此，由离心力过大防止限幅器14f进行的模型前轮舵角δf_d的限制处理，可以说是依据车辆的离心力、车辆重心点侧滑角(或者后轮侧滑角)的收敛预想值来限制模型前轮舵角δf_d的一种处理。对此，前轮侧滑角限幅器14d的限制处理，可以说是为了使得模型车辆的前轮侧滑角βf_d的瞬间值不至于变得过大而被用来限制模型前轮舵角δf_d的一种处理。

另外，在本实施方式中，在由离心力过大防止限幅器14f来设定容许范围[δf_min、δf_max]时所使用的函数fe，虽然是如所述图7(a)、(b)所示那样来进行设定，但也不仅限于此。 

例如，也可以如图8中的实线所示那样来设定函数fe(γd、γmax)。在该例子中，fe(γd、γmax)的值随着γd值的增加(从负侧的值向正侧的值增加)而单调减少，而且当γd＝γmax时则变成0。同时，此时，函数fe(—γd、—γmax)则为图8中的虚线所示。这种情况下，一旦γd超过γmax，随着γd的增加，由上述式子08a决定的模型前轮舵角δf_d的容许范围的上限值δf_max则比正常转圈时限界舵角δf_max_c更接近于0。同样，一旦γd在负向一侧超过—γmax时，随着γd的减少(数值大小上的增加)，由上述式子08b决定的模型前轮舵角δf_d的容许范围的下限值δf_min比—δf_max更接近于0。 

另外，也可以替换上述式子08a、08b，而通过下式11a、11b，来设定δf_d的容许范围的上限值δf_max及下限值δf_min，并同时，例如如图9实线、虚线所示那样来分别设定函数fe(γd、γmax)、fe(—γd、—γmax)。 

δf_max＝δf_max_c·fe(γd、γmax)  ……式11a 

δf_min＝—δf_max_c·fe(—γd、—γmax)  ……式11b 

在该例子中，fe(γd、γmax)、fe(—γd、—γmax)的值一直处于1以上，并以与图7(a)、(b)相同的形态，根据γd而变化。而且，通过分别将这些fe(γd、γmax)、fe(—γd、—γmax)的值乘上δf_max_c、δf_min_c来设定上限值δf_max和下限值δf_min。 

另外，也可以取代通过函数fe的值补正δf_max_c来设定模型前轮舵角δf_d的容许范围[δf_min、δf_max]，例如通过下述的处理，来决定第二限制完毕前轮舵角δf_ltd2。图10是用于说明其处理功能的功能方框图。 

即，在处理部14g，根据模型车辆的横摆比率γd(前次值)来决定前轮舵角补正量Δδf，该前轮舵角补正量Δδf用于对由所述减法器14e(参照图6)算出的第一限制完毕前轮舵角δf_ltd1进行补正。此时，如处理部14g内的坐标所表示的那样，Δδf一般是以下述方式来决定的：随着γd在正向一 侧的增加，Δδf的值在正向一侧呈单调增加，同时，随着γd在负向一侧的减少，Δδf的值在负向一侧呈单调减少。另外，在处理部14g的坐标中，Δδf的值设置了上限值(>0)及下限值(<0)。在这种情况下，上限值及下限值被设定为例如其绝对值是与所述图7(a)、(b)所示的一定值fex相同的值。 

接着，通过加法器14h将如上所述决定的前轮舵角补正量Δδf加上由所述减法器14e(参照图6)计算出的第一限制完毕前轮舵角δf_ltd1，来决定附带有输入补正的第一限制完毕前轮舵角。在这种情况下，在δf_ltd1的方向与γd的方向相互反向时，附带有输入补正的第一限制完毕前轮舵角的大小小于δf_ltd1的大小。但是，在δf_ltd1方向与γd的方向相同时，附带有输入补正的第一限制完毕前轮舵角的大小大于δf_ltd1的大小。 

此后，使该附带有输入补正的第一限制完毕前轮舵角通过离心力过大防止限幅器14f，由此来决定附带有输入补正的第二限制完毕前轮舵角，而该附带有输入补正的第二限制完毕前轮舵角是将附带有输入补正的第一限制完毕前轮舵角限制在模型前轮舵角δf_d的容许范围[δf_min、δf_max]内的值而成的操舵角。即，当附带有输入补正的第一限制完毕前轮舵角是容许范围内的值时，该附带有输入补正的第一限制完毕前轮舵角就原样被决定为附带有输入补正的第二限制完毕前轮舵角。另外，当附带有输入补正的第一限制完毕前轮舵角脱离出容许范围时，δf_max及δf_min之中的接近于附带有输入补正的第一限制完毕前轮舵角的值被决定为附带有输入补正的第二限制完毕前轮舵角。 

此时，离心力过大防止限幅器14f的模型前轮舵角δf_d的容许范围的上限值δf_max(>0)被设定为下述的值，即，预测δf_ltd1方向与γd的方向相同情况下的δf_ltd1的补正量，且比所述正常转圈时舵角界限值δf_max_c大的值(例如δf_max_c+fex)。同样，模型前轮舵角δf_d的容许范围的下限值δf_min(<0)被设定为其绝对值比δf_max_c大的值(例如—δf_max_c—fex)。 

此后，通过减法器14i，从所述决定的附带有输入补正的第二限制完毕前轮舵角之中减去所述前轮舵角补正量Δδf，由此来决定第二限制完毕前轮舵角δf_ltd2。 

即使如上所述方法来决定第二限制完毕前轮舵角δf_ltd2，也可以一面 使产生在模型车辆上的离心力不至于变得过大、且防止在实车1逆向操舵状态下受到不必要的限制，一面决定输入给标准动态特性模型16的模型前轮舵角δf_d(＝δf_ltd2)。 

另外，在本实施方式中，虽然为了决定输入给标准动态特性模型16的模型前轮舵角δf_d而进行了所述前轮侧滑角限幅器14d和离心力过大防止限幅器14f的处理，但是，也可以省略任意一方或者两者的处理。即，可以将由处理部14a决定的无限制时前轮舵角δf_unltd、或者将该δf_unltd输入给离心力过大防止限幅器14f之后得到的值、或者由所述减法器14e决定的第一限制完毕前轮舵角δf_ltd1决定为输入给标准动态特性模型16的模型前轮舵角δf_d。 

如上所述由标准操作量决定部14决定的模型前轮舵角δf_d的此次值(＝δf_ltd2的此次值)被输入给标准动态特性模型16，根据该输入值与如后述的由FB分配法则部20决定的假想外力Fvir、Mvir(前次值)，通过该标准动态特性模型16(按照上述式子01或式子103)，来重新决定标准横摆比率γd和标准车辆重心点侧滑角βd的此次值。另外，因为该处理实际上是通过用离散时间系列表述式子01后的式子来进行的，因此为了决定γd、βd的此次值，还需要使用γd、βd的前次值。另外，式子01中的特性调整用矩阵式K的模型特性调整用参数k2的值、或者式子103中的惯性量I’的值如前所述根据行驶速度Vd(＝Vact)来设定。 

在该情况下，因为输入给标准动态特性模型16的模型前轮舵角δf_d，如上所述，受到标准操作量决定部14限制，因此可防止模型车辆的旋转和极端性侧滑的产生。 

[关于FB分配法则部] 

下面，参照图11～图18，说明FB分配法则部20的详细处理。 

图11是表示FB分配法则部20的处理功能的功能方框图。如图所示的FB分配法则部20若大致区分其处理功能，则由执行决定假想外力Mvir、Fvir的处理的假想外力决定部20a和执行决定执行装置动作FB目标值的处理的执行装置动作FB目标值决定部20b构成。 

首先，参照图11，说明假想外力决定部20a，该假想外力决定部20a 的处理功能被大致分成假想外力暂定值决定部201和γβ限幅器202。 

在假想外力决定部20a的处理中，首先，对应于从所述减法器18输入的状态量偏差γerr(＝γact—γd)、βerr(＝βact—βd)，通过假想外力暂定值决定部201来决定假想外力的暂定值Mvirtmp、Fvirtmp。暂定值Mvirtmp、Fvirtmp中的Mvirtmp是表示为使状态量偏差γerr、βerr接近于0而应绕标准动态特性模型16的模型车辆重心点Gd附加性地产生的力矩(横摆方向的力矩)，Fvirtmp是表示为使状态量偏差γerr、βerr接近于0而应附加性地作用于标准动态特性模型16的模型车辆重心点Gd的平移力(模型车辆的横向平移力)。 

具体而言，如下式15所示，通过将规定的增益矩阵Kfvir乘上由被输入的状态量偏差γerr、βerr组成的矢量(γerr、βerr)^T(脚标T表示倒置)，来决定假想外力的暂定值Mvirtmp、Fvirtmp(以下，称作假想外力暂定值Mvirtmp、Fvirtmp)。 

【数5】 

          $\left[\begin{array}{c}\mathrm{Fvirtmp}\\ \mathrm{Mvirmp}\end{array}\right]=\mathrm{Kfvir}\&CenterDot;\left[\begin{array}{cc}\mathrm{\&beta;}& \mathrm{err}\\ \mathrm{\&gamma;}& \mathrm{err}\end{array}\right]$     ......式15 

其中 

          $\mathrm{Kfvir}\&equiv;\left[\begin{array}{cc}\mathrm{Kfvir}11& \mathrm{Kfvir}12\\ \mathrm{Kfvir}21& \mathrm{Kfvir}22\end{array}\right]$

根据该式15，为使状态量偏差γerr、βerr接近于0，作为向标准动态特性模型16反馈的控制输入暂定值，决定假想外力暂定值Mvirtmp、Fvirtmp。 

另外，如果想要只在模型车辆的车辆重心点侧滑角βd或实车1的实车辆重心点侧滑角βact快要超过规定的容许范围、以及超过规定的容许范围时，通过以下进行详细说明的γβ限幅器202强力产生使βd或βact返回到容许范围的作用，则最好用与时间常数较小的1次延迟特性相接近的特性来使βerr收敛于0。因此，例如只需将增益矩阵Kfvir组成部分中的Kfvir12设定为0、并以使Kfvir11的绝对值变大的形式进行设定即可。 

接着，通过γβ限幅器202来执行修正假想外力的暂定值Mvirtmp、 Fvirtmp的处理，以抑制标准动态特性模型16下的模型车辆的横摆比率γd和车辆重心点侧滑角βd分别从规定的容许范围脱离。 

具体而言，γβ限幅器202首先执行预测演算部203的处理。该预测演算部203预测规定时间后(1个以上的规定数的控制处理周期的时间后)的模型车辆的横摆比率γd和车辆重心点侧滑角βd，并将这些预测值分别作为预测横摆比率γda、预测车辆重心点侧滑角βda进行输出。 

此时，由标准动态特性模型16决定的标准横摆比率γd(此次值)及标准车辆重心点侧滑角βd(此次值)、由传感器/推定器12测定或推定出的实行驶速度Vact(此次值)、由标准操作量决定部14决定的第二限制完毕前轮舵角δf_ltd2(此次值)、以及由假想外力暂定值决定部201如上决定的假想外力暂定值Mvirtmp、Fvirtmp(此次值)被输入给预测演算部203。而且，该预测演算部203假定将模型前轮舵角δf_d保持在所输入的δf_ltd2，并将作用于模型车辆的假想外力Mvir、Fvir保持在所输入的Mvirtmp、Fvirtmp，且将模型车辆的行驶速度Vd保持在所输入的Vact，然后依据上述式子01或者式子103来计算出预测横摆比率γda及预测车辆重心点侧滑角βda。 

接着，γβ限幅器202使如上所述由预测演算部203计算出的γda、βda分别通过γ不灵敏区处理部204、β不灵敏区处理部205，由此来求解出γda、βda分别从规定的容许范围脱离出来的脱离量γover、βover。图中所示的γ不灵敏区处理部204的坐标是例举γda和γover关系的坐标，该坐标的横轴方向的值为γda值，纵轴方向的值为γover值。同样，图中所示的β不灵敏区处理部205的坐标是例举βda和βover关系的坐标，该坐标的横轴方向的值为βda值，纵轴方向的值为βover值。 

在此，γ不灵敏区处理部204的容许范围是以γdamin(<0)、γdamax(>0)分别为其下限值、上限值的容许范围(横摆比率γd的容许范围)；β不灵敏区处理部205的容许范围是以βdamin(<0)、βdamax(>0)分别为其下限值、上限值的容许范围(车辆重心点侧滑角βd的容许范围)。 

在本实施方式中，关于横摆比率γd的容许范围[γdamin、γdamax]如下所述进行设定，例如在将模型车辆的行驶速度Vd保持在Vact(此次值)、且同时将模型车辆的横摆比率γd保持在γdamin或γdamax而进行正常转 圈时，产生在模型车辆上的离心力不超过对应于推定摩擦系数μestm(此次值)的摩擦力的界限值。即，根据Vact(此次值)和μestm(此次值)，以满足下式16a、16b的形式来设定γdamax、γdamin。 

m·Vact·γdamax<μestm·m·g  ……式16a 

m·Vact·γdamin>—μestm·m·g    ……式16b 

例如只要使得各自的绝对值是与由上述式子05决定的正常转圈时最大横摆比率γmax相同值的方式来设定γdamax、γdamin即可(设γdamax＝γmax、γdamin＝—γmax)。但是，也可以将γdamax、γdamin设定为其绝对值与γmax不同的值(例如小于γmax的值)。 

另外，关于车辆重心点侧滑角βd的容许范围[βdamin、βdamax]例如被设定在：将实车1的车辆重心点侧滑角和作用于实车1重心点的横向的平移力之间的关系大致维持在线性关系(正比关系)的车辆重心点侧滑角的范围内。此时，优选是根据Vact(此次值)和μestm(此次值)中至少任意一方来设定βdamin、βdamax。 

而且，在γ不灵敏区处理部204的处理中，具体而言，当所输入的γda为规定的容许范围[γdamin、γdamax]内的值时(γdamin≤γda≤γdamax时)，则γover＝0，当γda<γdamin时，则γover＝γda—γdamin，当γda>γdamax时，则γover＝γda—γdamax。据此，可以求得预测横摆比率γda从容许范围[γdamin、γdamax]脱离出来的脱离量γover。 

同样，β不灵敏区处理部205的处理，当所输入的βda的值为规定的容许范围[βdamin、βdamax]内的值时(βdamin≤βda≤βdamax时)，则βover＝0，当βda<βdamin时，则βover＝βda—βdamin，当βda>βdamax时，则βover＝βda—βdamax。据此，可以求得预测车辆重心点侧滑角βda从容许范围[βdamin、βdamax]脱离出来的脱离量βover。 

接着，γβ限幅器202在处理部206计算出作为假想外力暂定值Mvirtmp、Fvirtmp的补正量的暂定值操作量Mvir_over、Fvir_over以使这些脱离量γover、βover接近于0。 

具体而言，如下式17所示，将规定的增益矩阵Kfov乘上由γover、βover 组成的矢量(γover、βover)^T，由此来决定Mvir_over、Fvir_over。 

【数6】 

          $\left[\begin{array}{c}\mathrm{Fvir}\_\mathrm{over}\\ \mathrm{Mvir}\_\mathrm{over}\end{array}\right]=\mathrm{Kfov}\&CenterDot;\left[\begin{array}{cc}\mathrm{\&beta;}& \mathrm{over}\\ \mathrm{\&gamma;}& \mathrm{over}\end{array}\right]$     ......式17 

其中 

             $\mathrm{Kfov}\&equiv;\left[\begin{array}{cc}\mathrm{Kfov}11& \mathrm{Kfov}12\\ \mathrm{Kfov}21& \mathrm{Kfov}22\end{array}\right]$

此后，γβ限幅器202通过减法器207分别从假想外力暂定值Mvirtmp、Fvirtmp减去该暂定值操作量Mvir_over、Fvir_over，由此来决定假想外力Mvir、Fvir的此次值。即，通过下式18a、18b来决定假想外力Mvir、Fvir。 

Mvir＝Mvirtmp—Mvir_over   ……式18a 

Fvir＝Fvirtmp—Fvir_over    ……式18b 

通过执行如上所述的假想外力决定部20a的处理，可以抑制预测横摆比率γda及预测车辆重心点侧滑角βda分别从容许范围[γdamin、γdamax]、[βdamin、βdamax]脱离出来，同时以使状态量偏差γerr、βerr接近于0的方式来决定假想外力Mvir、Fvir。 

另外，以上说明的假想外力决定部20a的γβ限幅器202，虽然是通过暂定值操作量Mvir_over、Fvir_over对假想外力暂定值Mvirtmp、Fvirtmp进行补正，由此来决定假想外力Mvir、Fvir(更一般而言，是通过Mvir_over与Mvirtmp的线性组合、和Fvir_over与Fvirtmp的线性组合来分别决定Mvir、Fvir)，但也可以以如下方式来决定假想外力Mvir、Fvir。图12是用于说明其处理的功能方框图。 

参照该图，在该例子中，假想外力暂定值决定部201、预测演算部203、γ不灵敏区处理部204、β不灵敏区处理部205、处理部206的处理与图11中的处理相同。另一方面，在本例中，由处理部206求得的暂定值操作量Fvir_over、Mvir_over被分别输入给处理部208、209，在该处理部208、209中，决定用于分别补正假想外力暂定值Mvirtmp、Fvirtmp的补正系数 Katt1(≥0)、Katt2(≥0)。这些补正系数Katt1、Katt2是分别乘以假想外力暂定值Mvirtmp、Fvirtmp的补正系数。另外，与图中所示的处理部208有关的坐标是例举Mvir_over与Katt1关系的坐标，该坐标的横轴方向的值为Mvir_over的值，纵轴方向的值为Katt1的值。同样，与图中所示的处理部209有关的坐标是例举Fvir_over与Katt2关系的坐标，该坐标的横轴方向的值为Fvir_over值，纵轴方向的值为Katt2值。 

在处理部208的处理中，如图中的坐标所示进行Katt1值的设定，即，当Mvir_over为0时，Katt1＝1，随着Mvir_over的绝对值从0开始增加，Katt1的值从1到0单调减少。而且，一旦Mvir_over的绝对值超过规定值(Katt1达到0的值)，Katt1值则被维持在0。 

同样，在处理部209的处理中，如图中坐标所示进行Katt2值的设定，即，当Fvir_over为0时，Katt2＝1，随着Fvir_over的绝对值从0开始增加，Katt2的值从1到0单调减少。而且，一旦Fvir_over的绝对值超过规定值(Katt2达到0的值)，Katt2的值被维持在0。 

接着，如上所述决定的补正系数Katt1、Katt2分别在乘法器210、211，与假想外力暂定值Mvirtmp、Fvirtmp进行乘法运算，据此，决定假想外力Mvir、Fvir的此次值。 

这样，在图12的例子中，决定假想外力Mvir，使得随着脱离量Mvir_over的绝对值变大，使假想外力Mvir的大小相对于假想外力暂定值Mvirtmp而减小(使之接近于0)。同样，决定假想外力Fvir，使得随着脱离量Fvir_over的绝对值变大，使假想外力Fvir的大小相对于假想外力暂定值Fvirtmp而减小(使之接近于0)。如此决定假想外力Mvir、Fvir意味着：认为γda、βda从容许范围脱离出是起因于假想外力Mvir、Fvir，在抑制γda、βda从容许范围[γdamin、γdamax]、[βdamin、βdamax]脱离出的同时，决定假想外力Mvir、Fvir以使状态量偏差γerr、βerr接近于0。此时，在标准操作量决定部14，最好是如前所述，对输入给标准动态特性模型16的模型前轮舵角δf_d进行限制。 

另外，在以上说明的γβ限幅器202中，将由预测演算部203并使用如上所述的式子01或者式子103求得的预测横摆比率γda及预测车辆重心点侧滑角βda分别输入向γ不灵敏区处理部204、β不灵敏区处理部205， 求出脱离量γover、βover，但这时，也可以取代γda、βda而使用标准横摆比率γd、标准车辆重心点侧滑角βd的此次值、或者实横摆比率γact、实车辆重心点侧滑角βact的此次值、或者对这些值施以滤波处理后的值。 

例如，也可以在各控制处理周期，取代γda将γd的此次值输入给γ不灵敏区处理部204，同时取代βda将下述值输入给β不灵敏区处理部205，即该值是指对通过标准动态特性模型16逐步计算出的βd施以滤波处理(T1、T2为某时间常数，s为拉普拉斯运算符)之后而得到的值，其中该滤波处理是以(1+T1·s)/(1+T2·s)形式表示传递函数的滤波处理。此时，例如，若以T1>T2的方式设定时间常数T1、T2，则该滤波处理起到作为相位推移补偿要素的作用。此时，推移一定程度的高频领域中的βd的频率分量的相位，提高相对于该频率分量的增益，由此可以在各控制处理周期决定的βd值自身从容许范围[βdamin、βdamax]脱离出之前，根据βover来限制假想外力Mvir、Fvir。 

另外，在预测演算部203，可以如下式19a、19b所示，使用适当的系数cij，对γd、βd的此次值进行线性组合，并将由此得到的值作为γda、βda求出。 

γda＝c11·γd+c12·βd   ……式19a 

βda＝c21·γd+c22·βd   ……式19b 

或者，也可以如下式20a、20b所示，使用适当的系数cij，对γd、βd、Mvirtmp、Fvirtmp、以及δf_ltd2的此次值进行线性组合并求出，并将由此得到的值作为γda、βda求出。 

γda＝c11·γd+c12·βd+c13·Mvirtmp 

+c14·Fvirtmp+c15·δf_ltd2   ……式20a 

βda＝c21·γd+c22·βd+c23·Mvirtmp 

+c24·Fvirtmp+c25·δf_ltd2   ……式20b 

另外，这些式子20a、20b是对上述的预测演算部203的处理进行了 更一般化表现的式子。 

或者，也可以如下式21a、21b所示那样，使用适当的系数cij，对γact、βact的此次值进行线性组合，并将由此得到的值作为γda、βda求出。 

γda＝c11·γact+c12·βact   ……式21a 

βda＝c21·γact+c22·βact   ……式21b 

或者，也可以如下式22a、22b所示那样，使用适当的系数cij，对γd、βd、βd的时间微分值dβd/dt、γact、βact、βact的时间微分值dβact/dt、Mvirtmp、Fvirtmp、以及δf_ltd2的此次值进行线性组合，将由此得到的值作为γda、βda求出。 

γda＝c11·γd+c12·βd+c13·dβd/dt+c14·γact 

+c15·βact+c16·dβact/dt+c17·Mvirtmp 

+c18·Fvirtmp+c19·δf_ltd2    ……式22a 

γda＝c21·γd+c22·βd+c23·dβd/dt+c24·γact 

+c25·βact+c26·dβact/dt+c27·Mvirtmp 

+c28·Fvirtmp+c29·δf_ltd2    ……式22b 

或者，也可以将式子20a右边演算结果的值与式子21a右边演算结果的值的加权平均值、以及式子20b右边演算结果的值与式子21b右边演算结果的值的加权平均值分别作为γda、βda来求解。另外，这是根据式子22a、式子22b来求解γda、βda的一个例子。另外，可以省略式子20a、式子20b、或者式子22a、式子22b中的Mvirtmp、Fvirtmp的项。 

或者，可以依据上述式子01或式子103，求解出到规定时间的各控制处理周期的γd、βd的预测值，并将所求得的γd、βd中的峰值决定为γda、βda。 

此外，即使是使用组合式20a、20b、或者组合式21a、21b、或者组 合式22a、22b中的任意一个组合式，来求解γda、βda时，也可以使这些式子的系数cij具有频率特性(换言之，对乘以cij的变量值施以低通滤波器等的滤波处理)。或者，也可以对乘以系数cij的变量值加以该变量的时间性变化率的限制。 

作为补充，在根据上述式子21a、式子21b、或者式子22a、式子22b来决定γda、βda时，优选设定各系数cij，使得该γda、βda具有作为某规定时间后的实车1的实横摆比率γact、实车辆重心点侧滑角βact的预测值的意思。 

另外，在标准动态特性模型16如上述式子01或是式子103所示为线形模型时，即使使用式子20a、式子20b、或者式子21a、式子21b、或者式子22a、式子22b中的任意一个，也可以适当地求解出作为实车1或模型车辆的某规定时间后的横摆比率及车辆重心点侧滑角的预测值的γda、βda。 

另外，代替γda、βda，而使用γact、βact的此次值、或者对γact、βact施以滤波处理而得到的值时，或者根据上述式子21a、式子21b、或式子22a、式子22b来决定γda、βda时，可以在抑制实车1的实横摆比率γact及实车辆重心点侧滑角βact的此次值、或者滤波值或预测值分别从容许范围[γdamin、γdamax]、[βdamin、βdamax]脱离的同时，以使状态量偏差γerr、βerr接近于0的方式来决定假想外力Mvir、Fvir。 

作为补充，在假想外力决定部20a的处理中，更一般而言，也可以根据下式200来决定假想外力Mvir、Fvir。 

【数7】 

$\left[\begin{array}{c}\mathrm{Fvir}\\ \mathrm{Mvir}\end{array}\right]=\left[\begin{array}{cccccc}\mathrm{Kfb}11& \mathrm{Kfb}12& \mathrm{Kfb}13& \mathrm{Kfb}14& \mathrm{Kfb}15& \mathrm{Kfb}16\\ \mathrm{Kfb}21& \mathrm{Kfb}22& \mathrm{Kfb}23& \mathrm{Kfb}24& \mathrm{Kfb}25& \mathrm{Kfb}26\end{array}\right]\&CenterDot;\left[\begin{array}{cc}\mathrm{\&beta;}& d\\ \mathrm{\&gamma;}& d\\ \mathrm{\&beta;}& \mathrm{act}\\ \mathrm{\&gamma;}& \mathrm{act}\\ \mathrm{\&beta;}& \mathrm{over}\\ \mathrm{\&gamma;}& \mathrm{over}\end{array}\right]$

                      $+\left[\begin{array}{cc}\mathrm{Kfb}\_\mathrm{\&delta;}& 1\\ \mathrm{Kfb}\_\mathrm{\&delta;}& 2\end{array}\right]\&CenterDot;\mathrm{\&delta;f}\_\mathrm{ltd}2$   ......式200

另外，虽然在所述γβ限幅器202的γ不灵敏区处理部204及β不灵敏区处理部205分别各自设定γda、βda的容许范围[γdamin、γdamax]、[βdamin、βdamax]，来决定脱离量γover、βover，但是，考虑到γda和βda之间的相关性，也可以针对γda、βda一组来设定容许范围(容许领域)，决定脱离量γover、βover。 

例如如图13所示，在以γda为横轴、以βda为纵轴的坐标平面上，由直线1～4围起来的区域A(平行四边形的区域)设定为γda、βda一组的容许区域A。此时，直线1、3是分别规定γda的下限值、上限值的直线，其下限值、上限值被设定为例如与所述γ不灵敏区处理部204的容许范围[γdamin、γdamax]的下限值γdamin、上限值γdamax相同。另外，直线2、4是分别规定βda的下限值、上限值的直线，在该例子中，设定该下限值及上限值分别对应γda呈线性变化。而且，例如，如下这样来决定脱离量γover、βover。即，γda、βda一组，如图13中的点P1所示处于容许区域A内时，γover＝βover＝0。另一方面，例如图13中的点P2所示，当γda、βda一组从容许区域A脱离出来时，决定点P3，该点P3是指：在穿过点P2并具有规定倾斜度的直线5上的点之中距离点P2最近的容许区域A的边界上的点(在直线5上且处于容许区域A内的点中与P2最近的点P3)。而且，将点P2处的γda的值和点P3处的γda的值之间的差决定为脱离量γover，同时将点P2处的βda的值和点P3处的βda的值之间的差决定为脱离量βover。另外，当与γda、βda一组对应的点例如是图13所示的点P4时，即，当穿过与γda、βda一组对应的点P4且具有规定的倾斜度(与直线5相同的倾斜度)的直线6与容许区域A不相交时(在直线6上不存在容许范围A内的点时)，则将P4点的γda的值与容许区域A内的点中距离该直线6最近的点P5的γda的值之间的差值决定为脱离量γover，而将点P4的βda的值和点P5的βda的值之间的差值决定为脱离量βover即可。 

作为补充，γda、βda一组的容许区域不一定为平行四边形的区域，例如，如图13中的虚线所示，也可以是光滑地形成临界部(形成为不具有锐角的角部)的区域A′。 

另外，在所述γβ限幅器202，虽然关于γda、βda两者，求解出从[γdamin、γdamax]、[βdamin、βdamax]脱离出来的脱离量γover、βover，并对应该脱 离量对暂定值Mvirtmp、Fvirtmp进行了补正，但也可以只对应γover、βover中的任意一方而对暂定值Mvirtmp、Fvirtmp进行补正。这种情况下，在所述处理部206的处理中，只需将γover、βover的任意一方的值固定为0来求解暂定值操作量Mvir_over、Fvir_over即可。 

下面，参照图14～图16说明执行装置动作FB目标值决定部20b的处理。另外，在以下的说明中，有时会将各车轮W1～W4称为第n轮Wn(n＝1、2、3、4)。 

图14是表示该执行装置动作FB目标值决定部20b的处理的功能方框图。参照该图，执行装置动作FB目标值决定部20b首先在处理部220，根据所输入的状态量偏差γerr、βerr，将反馈横摆力矩基本要求值Mfbdmd决定为针对实车1的执行机构3而言的反馈控制输入的基本要求值，而该反馈横摆力矩基本要求值Mfbdmd是指：为了使该状态量偏差γerr、βerr接近于0而绕实车1重心点G产生的横摆方向的力矩的基本要求值。 

具体而言，如下式23所示，将规定的增益矩阵Kfbdmd乘上由βerr、γerr组成的矢量(βerr、γerr)^T(对βerr、γerr进行线性组合)，由此决定Mfbdmd。 

【数8】 

         $\mathrm{Mfbdmd}=\mathrm{Kfbdmd}\&CenterDot;\left[\begin{array}{cc}\mathrm{\&beta;}& \mathrm{err}\\ \mathrm{\&gamma;}& \mathrm{err}\end{array}\right]$    ......式23 

其中 

     Kfbdmd≡[Kfbdmd1  Kfbdmd2] 

另外，也可以根据βerr、γerr和βerr的1阶微分值dβerr/dt来决定Mfbdmd。例如，可以将适当的增益矩阵乘上由βerr、γerr、dβerr/dt组成的矢量(通过适当的系数对βerr、γerr、dβerr/dt进行线性组合)来决定Mfbdmd。 

另外，可以将以(1+Tc1·s)/(1+Tc2·s)形式表示传递函数的相位补偿要素乘上增益矩阵Kfbdmd的要素Kfbdmd1及Kfbdmd2之中的至少任意一方。例如，设定时间常数Tc1、Tc2值，使得乘以βerr的Kfbdmd1又被乘上所述相位补偿要素，并且使Tc1>Tc2。在这种情况下，Kfbdmd1乘上 βerr后的项与下述的值等值，该值是指对βerr以及其微分值进行线性组合后的项又被通过高截止滤波器之后的值。 

接着，执行装置动作FB目标值决定部20b使该Mfbdmd通过不灵敏区处理部221，由此来决定不灵敏区超过反馈横摆力矩要求值Mfbdmd_a。另外，图中的不灵敏区处理部221的坐标是例举Mfbdmd和Mfbdmd_a之间关系的坐标，该坐标中的横轴方向的值为Mfbdmd值，纵轴方向的值为Mfbdmd_a值。 

在本实施方式中，在实车1的执行机构3的反馈控制中，为了使状态量偏差γerr、βerr接近于0，主要操作执行机构3之中的驱动/制动装置3A的刹车装置。此时，若根据如上所述决定的Mfbdmd来操作刹车装置，则该刹车装置恐怕会被频繁操作。在本实施方式中，为了防止该频繁操作，使Mfbdmd通过不灵敏区处理部221，并根据由此得到的不灵敏区超过反馈横摆力矩要求值Mfbdmd_a来操作刹车装置。 

该不灵敏区处理部221的处理具体如下所述执行。即，当Mfbdmd的值是处于被决定在0附近的规定的不灵敏区时，该不灵敏区处理部221的Mfbdmd_a＝0，当Mfbdmd大于该不灵敏区的上限值(>0)时，Mfbdmd_a＝Mfbdmd—上限值；当Mfbdmd小于该不灵敏区的下限值(<0)时，Mfbdmd_a＝Mfbdmd—下限值。换言之，将从Mfbdmd的不灵敏区超出的超过部分决定为Mfbdmd_a。根据由此决定的Mfbdmd_a来操作驱动/制动装置3A的刹车装置，由此可以抑制对应了状态量偏差γerr、βerr的刹车装置的频繁操作。 

然后，根据该不灵敏区超过反馈横摆力矩要求值Mfbdmd_a，决定所述执行装置动作FB目标值(针对执行机构3的反馈控制输入)的处理通过执行装置动作FB目标值分配处理部222执行。 

概略性地对该执行装置动作FB目标值分配处理部222的处理进行如下说明：即，决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n(n＝1、2、3、4)，以使绕实车1重心点产生Mfbdmd_a(进而使γerr、βerr接近于0)，其中，该FB目标第n轮刹车驱动/制动力Fxfbdmd_n是指由驱动/制动装置3A的刹车装置的动作而产生的各车轮W1～W4的驱动/制动力的反馈目标值(用于使γerr、βerr接近于0的刹车装置的反馈控制输入)。或者，除了决定Fxfbdmd_n(n＝1、2、3、4)之外，还决定自动操舵用FB目标横向力Fyfbdmd_f，而该自动操舵用FB目标横向力Fyfbdmd_f是指由转向装置3B的动作而产生的前轮W1、W2的横向力的反馈目标值。 

此时，在本实施方式中，当不灵敏区超过反馈横摆力矩要求值Mfbdmd_a为正向的力矩(从实车1上方观看时，逆时针转动方向的力矩)时，一般是在制动方向上增加实车1左侧车轮W1、W3的驱动/制动力，由此来决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n(n＝1、2、3、4)，以使绕实车1重心点G产生Mfbdmd_a。另外，当Mfbdmd_a为负向的力矩(从实车1上方观看时为顺时针方向的力矩)时，一般是在制动方向上增加实车1右侧车轮W1、W3的驱动/制动力，由此决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n(n＝1、2、3、4)，以使绕实车1重心点G产生Mfbdmd_a。 

在以后的说明中，如图15所示，实车1前轮W1、W2的间隔(即前轮W1、W2的轮距)设为df，后轮W3、W4的间隔(即后轮W3、W4的轮距)设为dr，前轮W1、W2的实操舵角(实前轮舵角)设为δf_act。另外，从上方观看实车1时，在与第n轮Wn的前后方向正交的方向(在水平面上正交的方向)上的该第n轮Wn与实车1重心点G之间的距离设为Ln(n＝1、2、3、4)。另外，在本实施方式中，后轮W3、W4是非操舵轮，因而省略了图示，但将后轮W3、W4的实操舵角(实后轮舵角)设为δr_act。在本实施方式中，δr_act＝0，L3＝L4＝dr/2。 

另外，图15中的Lf是实车1重心点G与前轮W1、W2的车轴之间的前后方向距离，Lr是实车1重心点G与后轮W3、W4的车轴之间的前后方向距离。这些Lf、Lr的值与所述图3所示的模型车辆所相关的Lf、Lr的值相同。 

下面具体说明执行装置动作FB目标值分配处理部222的处理。首先，设想实车1的前进行驶状态(δf_act＝0时的行驶状态)，在该前进行驶状态下，分别通过处理部222a_n(n＝1、2、3、4)来决定第n轮驱动/制动力最大要求值Fxfullfbdmd_n，而该第n轮驱动/制动力是指为了绕实车1重心点G产生与Mfbdmd_a相等的横摆方向的力矩而所需的第n轮Wn(n＝1、2、3、4)的驱动/制动力。 

具体而言，Fxfullfbdmd_n(n＝1、2、3、4)在各处理部222a_n，通过下式24a～24d的乘法演算决定。 

Fxfullfbdmd_1＝—(2/df)·Mfbdmd_a    ……式24a 

Fxfullfbdmd_2＝(2/df)·Mfbdmd_a      ……式24b 

Fxfullfbdmd_3＝—(2/dr)·Mfbdmd_a    ……式24c 

Fxfullfbdmd_4＝(2/dr)·Mfbdmd_a      ……式24d 

接着，执行装置动作FB目标值分配处理部222根据实前轮舵角δf_act，分别在处理部222b_1、222b2决定第一轮分配比率补正值K1_str及第二轮分配比率补正值K2_str，而且同时根据实后轮舵角δr_act，分别在处理部222b3、222b4决定第三轮分配比率补正值K3_str及第四轮分配比率补正值K4_str。这些第n轮分配比率补正值Kn_str(n＝1、2、3、4)是分别乘在Fxfullfbdmd_n上的补正系数。 

在此，实前轮舵角δf_act一旦从0开始变化，绕实车1的重心点G产生与Mfbdmd_a相等的横摆方向的力矩的第一轮W1及第二轮W2的驱动/制动力则分别从由上述式子24a、24b决定的Fxfullfbdmd_1、Fxfullfbdmd_2开始进行变化。同样，当后轮W3、W4为操舵轮时，实后轮舵角δr_act一旦从0开始变化，绕实车1重心点G产生与Mfbdmd_a相等的横摆方向的力矩的第三轮W3及第四轮W4的驱动/制动力则分别从由上述式子24c、24d决定的Fxfullfbdmd_3、Fxfullfbdmd_4开始进行变化。第n轮分配比率补正值Kn_str，一般是考虑到这种操舵角的影响而对Fxfullfbdmd_n(n＝1、2、3、4)进行补正，是用于决定下述第n轮Wn的驱动/制动力的补正系数，该第n轮Wn的驱动/制动力是指绕实车1重心点G产生与Mfbdmd_a相等或与之接近的横摆方向的力矩的驱动/制动力。 

但是，在本实施方式中，因为后轮W3、W4是非操舵轮，一直是δr_act＝0。因此，K3_str及K4_str实际一直被设定在“1”。所以，可以省略处理部222b3、222b4。 

另一方面，分别在处理部222b_1、222b_2，以下述方式决定与前轮 W1、W2有关的K1_str、K2_str。即，首先，根据预先设定的df、Lf的值和δf_act的值，并通过下式25a、25b的几何学演算来算出图15所示的L1、L2值。另外，作为该演算中的δf_act值，虽然可以使用由传感器/推定器12来测定或推定出的值(此次值)，但是，也可以使用实车1前轮W1、W2的操舵角的目标值(在各控制处理周期中最终决定的目标值)的前次值。另外，当转向装置3B是机器式转向装置时，也可以根据该机器式转向装置的总转向比和所述驾驶操作输入中的转向角θh来决定。或者，也可以使用由所述标准操作量决定部14的处理部14a决定的无限制时前轮舵角δf_unltd的此次值。 

L1＝(df/2)·cosδf_act-Lf·sinδf_act  ……式25a 

L2＝(df/2)·cosδf_act+Lf·sinδf_act  ……式25b 

此处，由于在前轮W1、W2的各自驱动/制动力上乘上L1、L2的值是绕实车1重心点G产生的横摆方向的力矩的值，因而，一般是使得K1_str＝(df/2)/L1，K2_str＝(df/2)/L2，并通过这两值分别乘上Fxfullfbdmd_1、Fxfullfbdmd_2来决定用于绕重心点G产生与Mfbdmd_a相等的横摆方向的力矩的前轮W1、W2的驱动/制动力。 

但是，这样的话，当L1或L2较小时，K1_str或K2_str则变得过大，与状态量偏差γerr、βerr对应的实车1整体的反馈环增益变得过大，容易产生控制***的振荡等。 

因此，在本实施方式中，通过下式26a、26b来决定K1_str、k2_str。 

K1_str＝(df/2)/max(L1，Lmin)……式26a 

K2_str＝(df/2)/max(L2，Lmin)……式26b 

此处，在式子26a、26b中，max(a、b)(a、b为一般变量)是对变量a、b中较大的值进行输出的函数，Lmin是比df/2小的正的常数。据此，防止K1_str、k2_str变得过大。换言之，在本实施方式中，将(df/2)/Lmin(＞1)作为K1_str、k2_str的上限值，根据实前轮舵角δf_act，将K1_str、k2_str设定在该上限值以下。 

另外，在本实施方式中，因为后轮W3、W4是非操舵轮，如前所述，K3_str＝K4_str＝1。但是，当后轮W3、W4为操舵轮时，最好是与对应于实前轮舵角δf_act而设定上述的K1_str、k2_str时同样，对应于实后轮舵角δr_act来设定K3_str、k4_str。 

然后，执行装置动作FB目标值分配处理部222在处理部222c_n(n＝1、2、3、4)，根据实前轮侧滑角βf_act(此次值)或者实后轮侧滑角βr_act(此次值)来决定第n轮分配增益Kn。该Kn是通过将它乘上第n轮驱动/制动力最大要求值Fxfullfbdmd_n来对Fxfullfbdmd_n进行补正的补正系数(比1小的正值)。 

这种情况下，第n轮分配增益Kn在各处理部222c_n被决定如下。 

关于实车1左侧前后配置的第一轮W1及第三轮W3的第一轮分配增益K1和第三轮分配增益K3分别如图16(a)、(b)中实线图所示，根据βf_act、βr_act而被决定为使其实质上呈连续变化。另外，关于实车1右侧前后配置的第二轮W2及第四轮W4的第二轮分配增益K2和第四轮分配增益K4分别如图16(a)、(b)中虚线图所示，根据βf_act、βr_act而被决定为使其实质上呈连续变化。此外，Kn的任意一个均是比1小的正值。并且，所谓“实质上连续”是表示：在用离散***表示模拟量时所必然产生的值的跳跃(量子化)不会损坏模拟量的连续性。 

这种情况下，进一步具体而言，关于第一轮分配增益K1及第三轮分配增益K3，K1如图16(a)中实线图所示，是根据βf_act的值来决定的，即，随着βf_act从负的值向正的值增加，K1从规定的下限值单调地增加到规定的上限值。因此，K1的值被决定为其值在βf_act为正值时比βf_act为负值时大。 

另一方面，K3如图16(b)中坐标实线所示，根据βr_act的值来决定，即，随着βr_act从负的值向正的值增加，K3从规定的上限值单调地减少到规定的下限值。因此，K3的值被决定为其值在βr_act为负的值时比βr_act为正的值时大。 

另外，图16(a)、(b)中的坐标实线以下述方式来设定，即，当βf_act、βr_act相互一致或者大体上一致时，与这些βf_act、βr_act对应的K1、k3 值之和大致为1。 

另外，关于第二轮分配增益K2及第四轮分配增益K4，K2如图16(a)中的坐标虚线所示，根据βf_act值来决定，即，随着βf_act从负的值向正的值增加，K2从规定的上限值单调地减少到规定的下限值。此时，表示K2和βf_act关系的虚线的坐标与将表示K1和βf_act关系的坐标实线以纵轴(βf_act＝0线)为中心进行左右翻转而成的坐标相同。因此，与βf_act各值对应的K2值被决定为：与βf_act正负值翻转后的值对应的K1值相等。 

另外，K4如图16(b)中虚线图所示，是根据βr_act的值来被决定的，即，随着βr_act从负的值增加到正的值，K4从规定的下限值单调地增加到规定的上限值。此时，表示K4与βr_act关系的虚线坐标，是与将表示K3与βr_act关系的坐标实线以纵轴(βr_act＝0线)为中心进行左右翻转而成的坐标相同。因此，与βr_act的各值对应的K4值被决定为：与翻转βr_act的正负值之后的值对应的K3值相等。 

通过如上所述地决定第n轮分配增益Kn(n＝1、2、3、4)，在实车1正常行驶等时βf_act与βr_act为大体上相同值的状况下，与前轮W1对应的第一轮分配增益K1和与该前轮W1正后方的后轮W3对应的第三轮分配增益K2之间的比率，一边将K1和K3的和几乎保持在一定，一边相对于βf_act及βr_act的变化而单调变化。同样，对应于前轮W2的第二轮分配增益K2与对应于该前轮W2正后方的后轮W4的第四轮分配增益K4之间的比率，一边将K2和K4的和几乎保持在一定，一边相对于βf_act及βr_act的变化而单调变化。 

对于根据βf_act、βr_act如上所述来决定第n轮分配增益Kn(n＝1、2、3、4)的理由将在后面进行说明。 

如上所述在决定Kn_str、Kn(n＝1、2、3、4)之后，执行装置动作FB目标值分配处理部222通过在处理部222b_n、222c_n将Kn_str、Kn分别乘上各第n轮驱动/制动力最大要求值Fxfullfbdmd_n(n＝1、2、3、4)，来决定第n轮分配驱动/制动力基本值Fxfb_n。即，通过下式27a～27d来决定第n轮分配驱动/制动力基本值Fxfb_n(n＝1、2、3、4)。

Fxfb_1＝Fxfullfbdmd_1·K1_str·K1  ……式27a 

Fxfb_2＝Fxfullfbdmd_2·K2_str·K2  ……式27b 

Fxfb_3＝Fxfullfbdmd_3·K3_str·K3  ……式27c 

Fxfb_4＝Fxfullfbdmd_4·K4_str·K4  ……式27d 

另外，如此决定Fxfb_n(n＝1、2、3、4)后，当Mfbdmd_a>0时，左侧车轮W1、W3的Fxfb_1、Fxfb_3为制动方向的驱动/制动力(负的驱动/制动力)；右侧车轮W2、W4的Fxfb_2、Fxfb_4为驱动方向的驱动/制动力(正的驱动/制动力)。同时，当Mfbdmd_a<0时，左侧车轮W1、W3的Fxfb_1、Fxfb_3为驱动方向的驱动/制动力(正的驱动/制动力)；右侧车轮W2、W4的Fxfb_2、Fxfb_4为制动方向的驱动/制动力(负的驱动/制动力)。 

此后，执行装置动作FB目标值分配处理部222，使如上决定的第n轮分配驱动/制动力基本值Fxfb_n(n＝1、2、3、4)分别通过与第n轮Wn对应的限幅器222d_n，来分别决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n，而该FB目标第n轮刹车驱动/制动力Fxfbdmd_n是指由驱动/制动装置3A的刹车装置的动作而产生的第n轮Wn的驱动/制动力的反馈目标值。 

此处，图14中的各限幅器222d_n(n＝1、2、3、4)的坐标是表示Fxfb_n和Fxfbdmd_n关系的坐标，该坐标的横轴方向的值为Fxfb_n值，纵轴方向的值为Fxfbdmd_n的值。 

该限幅器222d_n只有在输入给它的Fxfb_n的值为0或负值时，将Fxfb_n原样作为Fxfbdmd_n进行输出，当Fxfb_n为正值时，使不依据Fxfb_n值而进行输出的Fxfbdmd_n值为0。换而言之，以0为上限值对Fxfb_n加以限制，由此来决定Fxfbdmd_n。 

如上所述，通过分别决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n，如前述当Mfbdmd_a>0时，在制动方向上增加实车1左侧车轮W1、W3的驱动/制动力(设Fxfbdmd_1<0，Fxfbdmd_3<0)，由此来决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n(n＝1、2、3、4)以使得绕 实车1重心点G产生Mfbdmd_a。另外，在该情况下，在本实施方式中，关于右侧车轮W2、W4，设定Fxfbdmd_2＝Fxfbdmd_4＝0。 

另外，当Mfbdmd_a<0时，在制动方向上增加实车1右侧车轮W2、W4的驱动/制动力(Fxfbdmd_2<0，Fxfbdmd_4<0)，据此，以使得绕实车1重心点G产生Mfbdmd_a的方式来决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n(n＝1、2、3、4)。另外，在该情况下，在本实施方式中，关于左侧车轮W1、W3，使得Fxfbdmd_1＝Fxfbdmd_3＝0。 

并且，在任一情况下，因为所述第n轮分配增益Kn(n＝1、2、3、4)被决定为对应于βf_act或βr_act而实质上连续变化，因而，防止了Fxfbdmd_n会进行不连续变化的问题。 

此处，对第n轮分配增益Kn(n＝1、2、3、4)以根据βf_act、βr_act并以如前所述的倾向来进行决定的理由说明如下。 

首先，当Mfbdmd_a>0时，如前所述，决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n，以使得在制动方向上增加实车1的左侧车轮即第一轮W1及第三轮W3的驱动/制动力。 

而且，在这种情况下，设想βf_act<0、βr_act<0时的状况。在该状况下，暂时较大地设定K1值(进而使得Fxfbdmd_1在制动方向上变大)，同时较小地设定K3值(进而抑制Fxfbdmd_3在制动方向上变大)，这样，第一轮W1的横向力(其作用在于使得绕该实车1重心点产生与Mfbdmd_a相同方向的力矩)变小，另外，第三轮W3的横向力(其作用在于使得绕实车1重心点产生与Mfbdmd_a相反方向的力矩)变大。因此，有可能难以通过Mfbdmd_a充分产生绕实车1重心点G的所需的正向力矩(绕横摆轴的力矩)。因此，在βf_act<0、βr_act<0的状况下，将第一轮分配增益K1决定为较小的值，同时将第三轮分配增益K3决定为较大的值。 

此外，当Mfbdmd_a>0时，设想βf_act>0、βr_act>0时的状况。在这种状况下，暂时较小地设定K1值(进而抑制Fxfbdmd_1在制动方向上变大)，同时较大地设定K3值(进而使得Fxfbdmd_3在制动方向上变大)，这样第一轮W1的横向力(其作用在于使得绕实车1重心点产生与Mfbdmd_a反向的力矩)变大，另外，第三轮W3的横向力(其作用在于使得绕实车1重心点产生与Mfbdmd_a同向的力矩)变小。因此，有可能难以通过 Mfbdmd_a充分产生绕实车1重心点G的所需的负向力矩(绕横摆轴的力矩)。所以，在βf_act＞0、βr_act＞0的状况下，将第一轮分配增益K1决定为较大的值，同时将第三轮分配增益K3决定为较小的值。 

另外，当Mfbdmd_a＜0时，如前所述，决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n，以使得在制动方向上增加实车1的右侧车轮即第二轮W2及第四轮W4的驱动/制动力。 

而且，在这种情况下，设想βf_act＜0、βr_act＜0时的状况。在这种状况下，暂时较小地设定K2值(进而抑制Fxfbdmd_2在制动方向上变大)，同时较大地设定K4值(进而使得Fxfbdmd_4在制动方向上变大)，这样，第二轮W2的横向力(其作用在于使得绕实车1重心点产生与Mfbdmd_a反向的力矩)变大，另外，第四轮W4的横向力(其作用在于使得绕该实车1重心点产生与Mfbdmd_a同向的力矩)变小。因此，有可能难以通过Mfbdmd_a充分产生绕实车1重心点G的所需的负向力矩(绕横摆轴的力矩)。所以，在βf_act＜0、βr_act＜0的状况下，将第二轮分配增益K2决定为较大的值，同时将第四轮分配增益K4决定为较小的值。 

此外，当Mfbdmd_a＜0时，设想βf_act＞0、βr_act＞0时的状况。在这种状况下，暂时较大地设定K2值(进而使得Fxfbdmd_2在制动方向上变大)，同时较小地设定K4值(进而抑制Fxfbdmd_4在制动方向上变大)，这样第二轮W2的横向力(其作用在于使得绕实车1重心点产生与Mfbdmd_a同向的力矩)变小，另外，第四轮W4的横向力(其作用在于使得绕实车1重心点产生与Mfbdmd_a反向的力矩)变大。因此，有可能难以通过Mfbdmd_a充分产生绕实车1重心点G的所需的负向力矩(绕横摆轴的力矩)。所以，在βf_act＞0、βr_act＞0的状况下，将第二轮分配增益K2决定为较小的值，同时将第四轮分配增益K4决定为较大的值。 

如上所述，通过前述那样决定第n轮分配增益Kn(n＝1、2、3、4)，能够使对绕实车1重心点G产生Mfbdmd_a的横摆方向力矩有用的横向力不至于变得过小，并同时能使对绕实车1重心点G产生Mfbdmd_a的横摆方向力矩形成干扰的横向力不至于变得过大。 

另外，通过如前述那样决定第n轮分配增益Kn(n＝1、2、3、4)，正如实车1正常转圈时或正常前进时那样，在βf_act与βr_act一致或大体上呈一 致的状况下，K1值与K3值之和、以及K2值与K4值之和分别几乎为1。这就意味着，只要驱动/制动装置3A的刹车装置依照FB目标第n轮刹车驱动/制动力Fxfbdmd_n忠实地进行动作，则从Mfbdmd_a到实际绕实车1重心点G产生的力矩(横摆方向的力矩)的传递函数的增益大致为1(实际产生的横摆方向的力矩几乎等于Mfbdmd_a)。 

作为补充，在实车1的过渡性的动作状况等时，βf_act与βr_act间的差有时会变大。而且，在这种情况下，K1值与K3值之和、以及K2值与K4值之和会分别大幅偏离于1。为了解决这个问题，在如前所述那样决定K1、K3值之后，一面将这些值的比保持为一定，一面对K1、k3值进行修正，最好是使其修正后的K1、K3值之和大致为1。同样，在如前所述那样决定K2、K4值之后，一面将这些值的比保持为一定，一面对K2、K4值进行修正，最好是使修正后的K2、K4值之和大致为1。具体而言，只要在依照所述图16(a)、(b)的坐标决定第n分配增益Kn(n＝1、2、3、4)之后，通过K1′＝K1/(K1+K3)、K3′＝K3/(K1+K3)、k2′＝K2/(K2+K4)、K4′＝K4/(K2+K4)，求解出K1′、K2′、K3′、K4′，并将它们分别重新决定为K1、K2、K3、K4的值即可。 

另外，本实施方式中的执行装置动作FB目标值分配处理部222，除了如上述方式来决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n之外，还将所述反馈横摆力矩要求值Mfbdmd输入至处理部222e，通过该处理部222e来决定自动操舵用FB目标横向力Fyfbdmd_f，其中，该自动操舵用FB目标横向力Fyfbdmd_f是指由转向装置3B的动作而产生的前轮W1、W2的横向力的反馈目标值。在此，图中的处理部222e的坐标图是表示Mfbdmd与Fyfbdmd_f间关系的坐标图，该坐标图中的横轴方向的值为Mfbdmd的值，纵轴方向的值为Fyfbdmd_f的值。由该坐标图可知，在处理部222e处，基本上是随着Mfbdmd的增加，Fyfbdmd_f呈单调增加的方式来决定Fyfbdmd_f。这种情况下，Fyfbdmd_f是根据输入给处理部222e的Mfbdmd的值，例如使用图表而被决定的。 

另外，Fyfbdmd_f也可以通过将规定的增益乘以Mfbdmd来进行决定。另外，Fyfbdmd_f还可以在规定的上限值(>0)与下限值(<0)间的范围内根据Mfbdmd来决定。

作为补充，无论转向装置3B是主动式转向装置还是机械式转向装置，均可以省略处理部222e的处理。在通过处理部222e的处理决定自动操舵用FB目标横向力Fyfbdmd_f、并与之对应地操作转向装置3B的动作的情况下，更好的方案是，以使得下述两力矩之和与所述反馈横摆力矩基本要求值Mfbdmd大致相等的方式来决定Fxfbdmd_n(n＝1、2、3、4)及Fyfbdmd_f，其中，力矩之一是指：可通过FB目标第n轮刹车驱动/制动力Fxfbdmd_n(n＝1、2、3、4)而产生绕实车1重心点G的横摆方向的力矩，力矩之二是指：通过自动操舵用FB目标横向力Fyfbdmd_f而产生绕实车1重心点G的横摆方向的力矩。例如，也可以根据Mfbdmd与Mfbdmd_a之差来决定自动操舵用FB目标横向力Fyfbdmd_f。这种情况下最好是，当Mfbdmd_a＝0时，通过Fyfbdmd_f，以绕实车1重心点G产生与Mfbdmd大致相等的横摆方向的力矩的方式来决定Fyfbdmd_f。 

以上是本实施方式中的执行装置动作FB目标值决定部20b的详细处理。如上所述，通过该处理以使Mfbdmd接近于0(进而使状态量偏差γerr、βerr接近于0)的方式，将FB目标第n轮刹车驱动/制动力Fxfbdmd_n(n＝1、2、3、4)、或者将Fxfbdmd_n(n＝1、2、3、4)和自动操舵用FB目标横向力Fyfbdmd_f决定为执行装置动作FB目标值。 

另外，所述限幅器222d_n(n＝1、2、3、4)也可以将下述值作为Fxfbdmd_n进行输出，该值是指将输入给所述限幅器222d_n的Fxfb_n限制在比0稍微大的规定的正上限值以下而得到的值。例如，当Fxfb_n是该上限值以下的值时，将Fxfb n原样作为Fxfbdmd_n进行输出，当Fxfb_n是比上限值大的正的值时，将该上限值作为Fxfbdmd_n进行输出。在进行如此设定的情况下，正的值的Fxfbdmd_n是起到下述作用的反馈控制输入，该作用为使由刹车装置产生的第n轮Wn的制动方向的驱动/制动力的大小减少。 

另外，针对各车轮Wn(n＝1、2、3、4)，也可以对以下这些处理进行变更，使得利用所述这些处理所需要的输入值并使用图表或函数式来决定输出，所述这些处理包括：从处理部222a_n到限幅器222d_n的处理(是依据Mfbdmd_a、δf_act或δr_act、βf_act或βr_act来决定Fxfbdmd_n的处理)；或者从处理部222b_n到限幅器222d_n的处理(是依据Fxfullfbdmd_n、 δf_act或δr_act、βf_act或βr_act来决定Fxfbdmd_n的处理)；或者从处理部222c_n到限幅器222d_n的处理(是依据处理部222b_n的输出、βf_act或βr_act来决定Fxfbdmd_n的处理)；或者从处理部222a_n到限幅器222d_n的处理中结合了两个以上的部分的处理(例如从处理部222b_n到处理部222c_n的处理)。 

例如，在使用图表进行从处理部222c_n到限幅器222d_n的处理时，例如可以如图17(a)～(e)所示那样设定第一轮用的图表、如图18(a)～(e)所示那样设定第三轮用的图表。此时，图17(a)～(e)各自的图表是分别对应于βf_act的有代表性的各种值，并以各自的值作为坐标的横轴方向值、纵轴方向的值，来表示处理部222b_1的输出(＝Fxfullfbdmd_1·K1_str)与Fxfbdmd_1的关系。另外，图18(a)～(e)各自的图表是分别对应于βr_act的有代表性的各种值，并以各自的值作为坐标的横轴方向的值、纵轴方向的值，来表示处理部222b3的输出(＝Fxfullfbdmd_3·K3_str)与Fxfbdmd_3的关系。另外，在图17中，关于βf_act的值，“βf--”表示绝对值较大的负值，“βf-”表示绝对值较小的负值，“βf+”是表示绝对值较小的正值，“βf++”是表示绝对值较大的正值。同样，在图18中，关于βr_act的值，“βr--”是表示绝对值较大的负值，“βr-”是表示绝对值较小的负值，“βr+”是表示绝对值较小的正值，“βr++”是表示绝对值较大的正值。 

另外，第二轮用的图表虽省略了图示，但是，处理部222b2的输出(＝Fxfullfbdmd_2·K2_str)与Fxfbdmd_2的关系只要设定为与下述第一轮用图表相同即可，即该第一轮用图表是指与βf_act的各值符号翻转后的值相对应的图表(例如βf_act＝βf—时的处理部222b_2的输出(＝Fxfullfbdmd_2·K2_str)与Fxfbdmd_2的关系，跟βf_act＝βf+时的处理部222b_1的输出与Fxfbdmd_1的关系(图17(c)的图表所示的关系)相同)。同样，第四轮用的图表虽省略了图示，但处理部222b4的输出(＝Fxfullfbdmd_4·K4_str)与Fxfbdmd_4的关系只要设定为与下述第三轮用图表相同即可，即该第三轮用图表是指与βr_act的各值符号翻转后的值相对应的第三轮用图表(例如βr_act＝βr—时的处理部222b_4的输出(＝Fxfullfbdmd_4·K4_str)与Fxfbdmd_4的关系，跟βr_act＝βr+时的处理部 222b_3的输出与Fxfbdmd_3的关系(图18(c)的图表所示的关系)相同)。 

另外，在该例子中，当处理部222b_n(n＝1、2、3、4)的输出为0以下的值时，与所述图14所示的同样地决定Fxfbdmd_n。另一方面，当处理部222b_n(n＝1、2、3、4)的输出为正的值时，与将如上所述的限幅器222d_n的上限值设定为正的值的情况同样，Fxfbdmd_n在较小值的范围内为正的值。 

作为补充，关于第三轮W3和第四轮W4的所述处理部222b_3、222b_4，由于它们的输入值和输出值均相等，因而关于第三轮W3和第四轮W4，使用如上所述的图表进行从处理部222c_3到限幅器222d_3的处理、以及从处理部222c_4到限幅器222d_4的处理，这跟使用图表进行从处理部222b_3到限幅器222d_3的处理、以及从处理部222b_4到限幅器222d_4的处理是相同的。 

而且，在决定所述第n轮分配增益Kn(n＝1、2、3、4)时，还可以取代βf_act、βr_act而是根据实际车辆重心点侧滑角βact来决定所述第n轮分配增益Kn(n＝1、2、3、4)。在这种情况下，只要设定βact与第n轮分配增益Kn之间的关系，使其成为和前述的βf_act或βr_act与第n轮分配增益Kn之间的关系呈相同的倾向的关系即可。例如，只需按照如下所述的图表来决定第n轮分配增益Kn，该图表是将图16(a)、(b)的图表中坐标横轴的数值βf_act、βr_act分别置换成βact后的图表。 

或者，也可以根据实际车辆重心点侧滑角βact、实际横摆比率γact以及实行驶速度Vact、或是βact、γact、Vact以及实前轮舵角δf_act并使用图表或函数式来决定所述第n轮分配增益Kn(n＝1、2、3、4)。例如，将关于所述模型车辆的所述式子02a的βf_d、βd、γd、Vd、δf_d分别置换成βf_act、βact、γact、Vact、δf_act，并基于置换后得到式子，将所述βf_act与第一轮分配增益K1和第二轮分配增益K2的关系(用所述图16(a)的图表表示的关系)变换成βact、γact、Vact、δf_act与K1和K2的关系，并基于该变换后的关系，根据βact、γact、Vact、δf_act来决定K1以及K2。同样，将关于所述模型车辆的所述式子02b的βr_d、βd、γd、Vd分别置换成βr_act、βact、γact、Vact，基于置换后得到式子，将所述βr_act与第三轮分配增益K3和第四轮分配增益K4的关系(用所述图16(b)的图表 表示的关系)变换成βact、γact、Vact与K3和K4的关系，并基于该变换后的关系，根据βact、γact、Vact来决定K3以及K4。 

而且，如上所述，关于在根据βf_act、βr_act、或βact来决定第n轮分配增益Kn的情况时，可以代替βf_act、βr_act、βact而使用标准动态特性模型16的模型车辆的βf_d、βr_d、βd的值来决定第n轮分配增益Kn。在本实施方式中，由于是以使所述状态量偏差γerr、βerr接近于0的形式来操作实车1和模型车辆这两者的动作的，所以这些动作的状态量不会有大的乖离。所以，可以代替βf_act、βr_act、βact而使用模型车辆的βf_d、βrd、βd。 

或者，还可以比如根据实车1的βf_act、βr_act各自的加权平均值以及模型车辆的βf_d、βr_d各自的加权平均值来决定第n轮分配增益Kn，以及比如根据实车1的βact与模型车辆的βd的加权平均值来决定第n轮分配增益Kn。这种情况下，可以使权重带有频率特性(例如起到相位补偿要素作用的频率特性)。 

或者，根据βf_act或βr_act或βact决定第n轮分配增益Kn(n＝1、2、3、4)的第一暂定值，与此同时根据βf_d或βr_d或βd决定第n轮分配增益Kn的第二暂定值，然后可以将这些暂定值的加权平均值或加权平均值等的合成值决定为第n轮分配增益Kn。例如，根据βf_act并如所述图16(a)所示的图表那样来决定关于第一轮W1的K1的第一暂定值，同时根据βf_d并以与第一暂定值同样的方式决定K1的第二暂定值。在这种情况下，第二暂定值相对βf_d的变化的倾向可以与第一暂定值相对βf_act的变化的倾向相同。然后，将所述的第一暂定值和第二暂定值的加权平均值决定为第一轮分配增益K1。对于其他的第n轮分配增益K2、K3、K4也以同样的方法进行决定。 

此外，更好的方案是，不仅根据βf_act、βr_act或βact等使第n轮分配增益Kn(n＝1、2、3、4)的值发生变化、还根据推定摩擦系数μestm使第n轮分配增益Kn(n＝1、2、3、4)发生变化的方式来决定第n轮分配增益Kn(n＝1、2、3、4)的值。例如，在本实施方式中，如上所述，在根据βf_act、βr_ct来决定第n轮分配增益Kn时，最好是决定K1，使得μestm越小，使在βf_act为其绝对值较大的负值时的第一轮分配增益K1越小。另外， 最好是决定K3，使得μestm越小，使在βr_act为绝对值较大的正值时的第三轮分配增益K3越小。同样，最好是决定K2，使得μestm越小，使在βf_act为绝对值大的正值时的第二轮分配增益K2越小。另外，最好是决定K4，使得μestm变得越小，使在βr_act为绝对值大的负值时的第四轮分配增益K4越小。这是因为μestm越小，使第n轮Wn(n＝1、2、3、4)的制动方向的驱动/制动力增加时的该第n轮Wn的横向力的降低会变得越显著。 

另外，也可以根据第n轮的实着地负荷(作用于第n轮的路面反力中的竖直方向或与路面垂直方向的平移力的测定值或推定值)，来调整第n轮分配增益Kn(n＝1、2、3、4)的值(根据βf_act、βr_act(或根据βact、βf_d、βr_d、βd中任意一值)来进行设定的值)。这种情况下，最好是第n轮Wn的实着地负荷越小，将第n轮分配增益Kn的值决定为越小。 

或者，在设各第n轮Wn的实着地负荷为Fzact_n(n＝1、2、3、4)、设它们的总和为∑Fzact(＝Fzact_1+Fzact_2+Fzact_3+Fzact_4)时，可以根据各前轮W1、W2的实着地负荷之和(＝Fzact_1+Fzact_2)来调整与前轮W1、W2有关的第n轮分配增益K1、K2的值，或者可以根据它们之和相对于∑Fzact的比例(＝(Fzact_1+Fzact_2)/∑Fzact)来进行调整。同样，根据各后轮W3、W4的实着地负荷之和(＝Fzact_3+Fzact_4)来调整与后轮W3、W4有关的第n轮分配增益K3、K4，或者可以根据它们之和相对于∑Fzact的比例(＝(Fzact_3+Fzact_4)/∑Fzact)来调整。或者可以分别根据第n轮Wn相对实着地负荷的∑Fzact的比例(＝Fzact_n/∑Fzact)来调整各第n轮分配增益Kn(n＝1、2、3、4)的值。 

另外，在本实施方式中，作为针对驱动/制动装置3A的刹车装置而言的反馈控制输入(作为执行装置动作FB目标值)，来决定FB目标第n轮刹车驱动/制动力Fxfbdmd_n(n＝1、2、3、4)，但是也可以取代Fxfbdmd_n而决定由刹车装置产生的各车轮Wn(n＝1、2、3、4)的目标打滑比，或者决定该目标打滑比与Fxfbdmd_n这两者。 

另外，为了决定Fxfbdmd_n等的执行装置动作FB目标值，可以不决定中间变量Mfbdmd和Mfbdmd_a、而是根据状态量偏差γerr、βerr并使用图表等来直接决定执行装置动作FB目标值。例如，可以使用将γerr、 γd(或γact)、βd(或βact)、Vact、μestm等变量作为输入的多次元的图表，来决定执行装置动作FB目标值。 

另外，不仅使状态量偏差γerr、βerr接近于0的方式来决定反馈横摆力矩基本要求值Mfbdmd，还可以使由所述假想外力决定部20a的γβ限幅器202求得的脱离量γover、βover接近于0地(进而抑制所述γda、βda从各自的容许范围[γdamin、γdamax]、[βdamin、βdamax]脱离出来的方式)来决定Mfbdmd。例如，可以使用适当的系数Kfbdmd1～Kfbdmd4，通过下式28a来决定Mfbdmd。 

Mfbdmd＝Kfbdmd1·γerr+Kfbdmd2·βerr 

—Kfbdmd3·γover—Kfbdmd4·βover  ……式28 

作为补充，在本实施方式中，如上所述，通过γβ限幅器202，使γover、βover接近于0的方式操控假想外力暂定值Mvirtmp、Fvirtmp，来决定假想外力Mvir、Fvir。即使仅此而已，模型车辆的γd、βd也分别会以不脱离出容许范围[γdamin、γdamax]、[βdamin、βdamax]的方式进行变化。而且，伴随与此，执行装置动作FB目标值发生变化以使实车1的γact、βact分别接近于γd、βd，因此，即使是在决定执行装置动作FB目标值，以使γerr、βerr接近于0的情况下，其结果也可以抑制γact、βact从容许范围[γdamin、γdamax]、[βdamin、βdamax]中脱离出来。其中，如上所述，除了γerr、βerr之外，通过以使γover、βover也接近于0的方式来决定Mfbdmd(进而决定执行装置动作FB目标值)，因而可更有效地抑制γact、βact分别从容许范围[γdamin、γdamax]、[βdamin、βdamax]中脱离。 

[关于FF法则] 

下面，参照图19，进一步详细说明所述FF法则部22的处理。图19是表示FF法则部22处理的功能方框图。 

如上所述，在本实施方式中，FF法则部22决定的前馈目标值(对应了驾驶操作输入的执行机构3的基本目标值)包括：由驱动/制动装置3A的刹车装置产生的实车1的各车轮W1～W4的驱动/制动力的前馈目标值 (以下，称为FF目标第n轮刹车驱动/制动力(n＝1、2，3、4))、由驱动/制动装置3A的驱动***产生的实车1的驱动轮W1、W2的驱动/制动力的前馈目标值(以下，称为FF目标第n轮驱动***驱动/制动力(n＝1、2)、驱动/制动装置3A的变速装置的减速比(变速比)的前馈目标值(以下，称为FF目标变速减速比)、以及由转向装置3B产生的实车1的操舵轮W1、W2的舵角的前馈目标值(以下，称为FF目标前轮舵角δf_ff)。 

如图19所示，根据驾驶操作输入之中的转向角θh(或根据θh和Vact)，并通过处理部230决定FF目标前轮舵角δf_ff。在图19中，设想转向装置3B为所述执行装置驱动型的转向装置的情形。在这种情况下，处理部230通过与所述标准操作量决定部14的处理部14a的处理相同的处理，来决定FF目标前轮舵角δf_ff。即，用转向角θh除以规定的总转向比is或者除以根据Vact设定的总转向比is来决定δf_ff。这样所决定的δf_ff的值，与通过所述标准操作量决定部14的处理部14a决定的无限制时前轮舵角δf_unltd的值相同。 

另外，在转向装置3B为所述执行装置辅助型的转向装置时，或者为机械式转向装置时，不需要决定δf_ff。或者也可以将δf_ff一直设定为0。但是，在转向装置3B为执行装置辅助型的转向装置，并具有：根据Vact，而对根据转向角θh机械性确定的前轮W1、W2的操舵角进行补正的功能时，可以根据Vact来决定其补正量，并将该补正量决定为δf_ff。 

作为补充，在转向装置3B为执行装置辅助型的转向装置时，因为前轮W1、W2的基本的舵角(δf_act的基本值)是根据转向角θh而机械性地确定的，故而δf_ff具有作为由执行装置确定的前轮W1、W2的舵角的补正量的前馈目标值的意思。 

另外，根据驾驶操作输入中的刹车踏板操作量，并分别通过处理部231a_n(n＝1、2、3、4)来决定FF目标第n轮刹车驱动/制动力(n＝1、2、3、4)。图中的各处理部231a_n所示的图表是分别例举刹车踏板操作量与FF目标第n轮刹车驱动/制动力(n＝1、2、3、4)之间关系的图表，该图表中的横轴方向的值为刹车踏板操作量的值，纵轴方向的值为FF目标第n轮刹车驱动/制动力。如图中的图表所示，FF目标第n轮刹车驱动/制动力(<0)基本上是随着刹车踏板操作量的增加而其大小(绝对值)呈单调增加的方式被决定的。另外，在图示的例子中，为使FF目标第n轮刹车驱动/制动力的大小不至于变得过大，一旦刹车踏板操作量超过规定量(相对于刹车踏板操作量的增加，FF目标第n轮刹车驱动/制动力的绝对值的增加率接近于0或成为0)时，FF目标第n轮刹车驱动/制动力会呈现饱和。 

根据驾驶操作输入中的油门踏板操作量、Vact和变速杆位置，并通过驱动系执行装置动作FF目标值决定部232来决定FF目标第n轮驱动系驱动/制动力(n＝1、2)和FF目标变速减速比。由于该驱动系执行装置动作FF目标值决定部232的处理可以与下述手法相同，故而本说明书中省略其详细说明，该手法是指，在公知的一般汽车中，根据油门踏板操作量和Vact和变速装置的变速杆位置，来决定从发动机传递给驱动轮的驱动力和变速装置的减速比。 

以上是本实施方式中的FF法则部22的具体处理的内容。 

[关于执行装置动作目标值合成部] 

下面，详细说明所述执行装置动作目标值合成部24的处理。图20是表示该执行装置动作目标值合成部24的处理的功能方框图。 

参照该图，执行装置动作目标值合成部24，关于第一轮W1，用加法器240来求解出所述执行装置动作FF目标值之中的FF目标第一轮刹车驱动/制动力与FF目标第一轮驱动***驱动/制动力之和，并将该和作为FF综合目标第一轮驱动/制动力FFtotal_1输入给最佳目标第一驱动/制动力决定部241a_1。此外，用加法器242来求解出该FFtotal_1与所述执行装置动作FB目标值之中的FB目标第一轮刹车驱动/制动力Fxfbdmd_1之和，并将该和作为无限制目标第一轮驱动/制动力Fxdmd_1输入给最佳目标第一驱动/制动力决定部241a_1。 

另外，执行装置动作目标值合成部24关于第二轮W2用加法器243求解出所述执行装置动作FF目标值之中的FF目标第二轮刹车驱动/制动力与FF目标第二轮驱动***驱动/制动力之和，并将上述和作为FF综合目标第二轮驱动/制动力FFtotal_2输入给最佳目标第二驱动/制动力决定部241a_2。此外，用加法器244来求解出该FFtotal_2与所述执行装置动作FB目标值之中的FB目标第二轮刹车驱动/制动力Fxfbdmd_2之和， 并将该和作为无限制目标第二轮驱动/制动力Fxdmd_2输入给最佳目标第二驱动/制动力决定部241a_2。 

另外，执行装置动作目标值合成部24关于第三轮W3将所述执行装置动作FF目标值之中的FF目标第三轮刹车驱动/制动力原样地作为FF综合目标第三轮驱动/制动力FFtotal_3输入给最佳目标第三驱动/制动力决定部241a_3。此外，用加法器245来求解出该FFtotal_3与所述执行装置动作FB目标值之中的FB目标第三轮刹车驱动/制动力Fxfbdmd_3之和，并将该和作为无限制目标第三轮驱动/制动力Fxdmd_3输入给最佳目标第三驱动/制动力决定部241a_3。 

另外，执行装置动作目标值合成部24关于第四轮W4将所述执行装置动作FF目标值中的FF目标第四轮刹车驱动/制动力原样地作为FF综合目标第四轮驱动/制动力FFtotal_4输入给最佳目标第四驱动/制动力决定部241a_4。此外，用加法器246来求解出该FFtotal_4与所述执行装置动作FB目标值之中的FB目标第四轮刹车驱动/制动力Fxfbdmd_4之和，并将该和作为无限制目标第四轮驱动/制动力Fxdmd_4输入给最佳目标第四驱动/制动力决定部241a_4。 

在此，所述FF综合目标第n轮驱动/制动力FFtotal_n(n＝1、2、3、4)，将其一般化而言，是表示由驱动/制动装置3A的驱动***动作产生的第n轮Wn的驱动/制动力的前馈目标值(FF目标第n轮驱动***驱动/制动力)与由刹车装置动作产生的第n轮Wn的驱动/制动力的前馈目标值(FF目标第n轮刹车驱动/制动力)之总和。此时，在本说明书的实施方式中，因为设前轮W1、W2为实车1的驱动轮，后轮W3、W4为从动轮，故而关于后轮W3、W4，FF目标第n轮刹车驱动/制动力(n＝3、4)原样地被决定为FF综合目标第n轮驱动/制动力FFtotal_n。 

另外，所述无限制目标第n轮驱动/制动力Fxdmd_n(n＝1、2、3、4)由于是所述FF综合目标第n轮驱动/制动力FFtotal_n与所述FB第n轮刹车驱动/制动力之和，故而是表示被驱动/制动装置3A的前馈控制动作(至少与驾驶操作输入对应的前馈控制动作)和反馈控制动作(至少与状态量偏差γerr、βerr对应的反馈控制动作)所要求的第n轮的总驱动/制动力。

而且，执行装置动作目标值合成部24通过最佳目标第n驱动/制动力决定部241a_n(n＝1、2、3、4)，分别决定第n轮Wn的驱动/制动力的最终目标值即目标第n轮驱动/制动力Fxcmd_n，而且同时决定第n轮的打滑比的最终目标值即目标第n轮打滑比。 

这种情况下，除了FFtotal_n和Fxdmd_n之外，还有第n轮Wn的实侧滑角(具体而言，当n＝1、2时，为实前轮侧滑角βf_act，当n＝3、4时，为实后轮侧滑角βr_act)的最新值(此次值)和推定摩擦系数μestm的最新值(此次值)也被输入给最佳目标第n驱动/制动力决定部241a_n(n＝1、2、3、4)。另外，虽省略图示，但是，实前轮舵角δf_act的最新值(此次值)也被输入给与前轮W1、W2有关的最佳目标第n驱动/制动力决定部241a_n(n＝1、2)。而且，最佳目标第n驱动/制动力决定部241a_n(n＝1、2、3、4)依据提供给各决定部的输入，以如后所述的方式来决定目标第n轮驱动/制动力Fxcmd_n和目标第n轮打滑比。 

另外，执行装置动作目标值合成部24将所述执行装置动作FB目标值之中的自动操舵用FB目标横向力Fyfbdmd_f和所述执行装置动作FF目标值之中的FF目标前轮舵角δf_ff输入给最佳目标自动舵角决定部247，通过该最佳目标自动舵角决定部247来决定前轮W1、W2的最终舵角的目标值即目标前轮舵角δ_fcmd。另外，当转向装置3B为所述执行装置驱动型的转向装置时，该δ_fcmd是表示由执行装置动作产生的前轮W1、W2的舵角本身(以实车1的前后方向为基准的舵角)的最终目标值，而当转向装置3B为所述执行装置辅助型的转向装置时，该δfcmd是表示由执行装置的动作产生的前轮W1、W2的操舵角的补正量的最终目标值。 

另外，执行装置动作目标值合成部24将所述执行装置动作FF目标值之中的FF目标第n轮驱动***驱动/制动力(n＝1、2)原样地作为目标第n轮驱动***驱动/制动力予以输出，该目标第n轮驱动***驱动/制动力是指由驱动/制动装置3A的驱动***动作产生的第n轮Wn的驱动/制动力的最终目标值。此外，执行装置动作目标值合成部24将所述执行装置动作FF目标值之中的FF目标变速减速比原样地作为目标变速减速比予以输出，该目标变速减速比是指驱动/制动装置3A的变速装置的减速比(变速比)的最终目标值。 

以下详细说明所述最佳目标第n驱动/制动力决定部241a_n(n＝1、2、3、4)的处理。图21是表示各最佳目标第n驱动/制动力决定部241a_n的处理的流程图。 

参照该图，首先，在S100中，当第n轮Wn(n＝1、2、3、4)的侧滑角为实侧滑角(具体而言，当n＝1、2时，为实前轮侧滑角βf_act，当n＝3、4时，为实后轮侧滑角βr_act)，路面摩擦系数(第n轮Wn和路面之间的摩擦系数)为推定摩擦系数μestm，并以此为前提条件，并在该前提条件的基础上，求解出：作为与无限制目标第n轮驱动/制动力Fxdmd_n最接近(包括相一致的情况)的第n轮Wn的驱动/制动力的值、即第n轮驱动/制动力候补Fxcand_n和与该第n轮驱动/制动力候补Fxcand_n对应的第n轮Wn的打滑比的值、即第n轮打滑比候补Scand_n。 

此处，在一般情况下，在各车轮的侧滑角、路面反力(驱动/制动力、横向力以及着地负荷)、打滑比以及路面摩擦系数之间，具有与该车轮轮胎的特性或悬架装置的特性对应的一定的相关关系。例如，在各车轮的侧滑角、路面反力(驱动/制动力、横向力以及着地负荷)、打滑比以及路面摩擦系数之间，具有由所述非专利文件1中的式子(2.57)、(2.58)、(2.72)、(2.73)表示的相关关系。另外，例如，当使着地负荷及路面摩擦系数为一定的情况时，在各车轮的侧滑角、驱动/制动力、横向力以及打滑比之间，具有所述非专利文件1中的图2.36所示的相关关系。因此，在侧滑角及路面摩擦系数分别为某值时的各车轮的路面反力和打滑比不是分别采用独立的值，而是各自的值按照所述相关关系(以下，称作车轮特性关系)而变化。另外，打滑比在驱动/制动力是驱动方向的驱动/制动力(>0)时为负的值，而在驱动/制动力是制动方向的驱动/制动力(<0)时为正的值。 

所以，在本实施方式的S100的处理中，根据表示第n轮Wn的侧滑角、路面摩擦系数、驱动/制动力以及打滑比之间关系的预先制作的图表，并根据第n轮Wn的实侧滑角βf_act或βr_act(最新值)和推定摩擦系数μestm(最新值)来求解出：与无限制目标第n轮驱动/制动力Fxdmd_n最接近或与之一致的驱动/制动力(与Fxdmd_n之差的绝对值为最小时的驱动/制动力)、以及与该驱动/制动力对应的打滑比。而且，将如此求得的驱动/制动力和打滑比分别决定为第n轮驱动/制动力候补Fxcand_n、和 第n轮打滑比候补Scand_n。 

另外，在该处理中使用的图表，例如只要通过各种实验等，或者依据车轮W1～W4的轮胎特性或悬架装置3C的特性，预先确定或者推定所述车轮特性关系，并依据该确定或推定的车轮特性关系来制作即可。另外，也可以将第n轮Wn的着地负荷作为变量参数加入到该图表中。在这种情况下，只要将第n轮Wn的实着地负荷Fzact_n输入给最佳目标第n驱动/制动力决定部241a_n，根据第n轮Wn的实侧滑角βf_act或βr_act、推定摩擦系数μestm以及实着地负荷Fzact_n，来决定Fxcand_n、Scand_n即可。其中，由于实着地负荷Fzact_n的变动一般情况下较小，故而可以将该实着地负荷Fzact_n视为一定值。 

作为补充，对应于第n轮Wn的实侧滑角βf_act或βr_act与推定摩擦系数μestm一组，或者对应于它们与实着地负荷Fzact_n一组，当Fxdmd_n处于可能在第n轮Wn产生(可从路面作用来的)的驱动/制动力(按照所述车轮特性关系可能产生的驱动/制动力)的值的范围内时，可以将该Fxdmd_n原样地决定为Fxcand_n。并且，当Fxdmd_n脱离出该范围时，可以将该范围中的上限值(>0)及下限值(<0)之中与Fxdmd_n接近的一方决定为Fxcand_n。 

另外，对应于第n轮Wn的实侧滑角βf_act或βr_act、与推定摩擦系数μestm一组，或者对应于它们与着地负荷Fzact_n一组，可能在第n轮Wn产生的打滑比与驱动/制动力的关系(按照所述车轮特性关系可能产生的打滑比与驱动/制动力的关系)一般情况下，是驱动/制动力相对于该打滑比的变化而具有峰值(极限值)的关系(是指在以打滑比为横轴的值、以驱动/制动力的大小为纵轴的值时，向上凸起的图表)。因此，有时会存在2种与绝对值小于其峰值的驱动/制动力值对应的打滑比的值。当与Fxcand_n对应的打滑比的值有2种存在时，只要将这2种打滑比的值之中更接近于0的一方的打滑比的值决定为第n轮打滑比候补Scand_n即可。换言之，在第n轮Wn的打滑比和驱动/制动力的关系(依照所述车轮特性关系的关系)中，只要是在驱动/制动力为峰值时的打滑比的值和0之间的范围内决定第n轮打滑比候补Scand_n即可。 

作为补充，在驱动/制动力为峰值时的打滑比的值和0之间的范围内，随着打滑比的绝对值从0开始增加，驱动/制动力的绝对值呈单调增加。 

接着，进入S102，在与S100相同的前提条件的基础之上，决定最大力矩产生时第n轮驱动/制动力Fxmmax_n和作为与之对应的打滑比的最大力矩产生时第n轮打滑比Smmax_n。在此，最大力矩产生时第n轮驱动/制动力Fxmmax_n是指：当第n轮Wn的侧滑角为实侧滑角βf_act或βr_act而路面摩擦系数为推定摩擦系数μestm时，可能在第n轮Wn产生的路面反力(具体而言，按照所述车轮特性关系，从路面可能作用于第n轮Wn上的驱动/制动力和横向力之合力)之中，通过该路面反力而绕实车1重心点G产生的横摆方向的力矩朝着与所述反馈横摆力矩基本要求值Mfbdmd的极性相同的极性(朝向)成为最大时的路面反力的驱动/制动力分量的值。另外，在这种情况下，Fxmmax_n、Smmax_n是在第n轮Wn的驱动/制动力和打滑比的关系(按照所述车轮特性关系的关系)中，在驱动/制动力的绝对值随着打滑比的绝对值从0开始增加而呈单调增加的区域内被决定。因此，Smmax_n是在驱动/制动力为峰值时的打滑比的值与0之间被决定。 

在S102中，关于前轮W1、W2(n＝1或2时)，例如根据实前轮侧滑角βf_act、推定摩擦系数μestm以及实前轮舵角δf_act，依据预先制作的图表(表示前轮侧滑角、路面摩擦系数、前轮舵角、最大力矩产生时驱动/制动力和最大力矩产生时打滑比之间的关系(按照所述车轮特性关系的关系)的图表)来决定最大力矩产生时第n轮驱动/制动力Fxmmax_n和与之对应的最大力矩产生时第n轮打滑比Smmax_n。或者，也可以用下述方式来决定Fxmmax_n和Smmax_n；即，依据表示前轮侧滑角、路面摩擦系数、打滑比、驱动/制动力和横向力之间关系的图表以及实前轮舵角δf_act，从对应于βf_act和μestm一组而可能产生的第n轮Wn(n＝1或2)的驱动/制动力和横向力的一组中，探索性地来决定这些合力绕实车1重心点G产生的横摆方向的力矩为最大时的驱动/制动力和横向力的一组，并将与该组对应的驱动/制动力和打滑比分别决定为Fxmmax_n、Smmax_n。 

另外，关于后轮W3、W4(n＝3或4时)，例如，根据实后轮侧滑角βr_act和推定摩擦系数μestm，并依据预先制作的图表(表示后轮侧滑角、路面摩擦系数、最大力矩产生时驱动/制动力和最大力矩产生时打滑比之间关系 (依照所述车轮特性关系的关系)的图表)，来决定最大力矩产生时第n轮驱动/制动力Fxmmax_n和与之对应的最大力矩产生时第n轮打滑比Smmax_n。或者，也可以用下述方式来决定Fxmmax_n和Smmax_n：即，依据表示后轮侧滑角、路面摩擦系数、打滑比、驱动/制动力和横向力之间关系的图表，从对应于βr_act和μestm一组而可能产生的第n轮Wn(n＝3或4)的驱动/制动力和横向力的组中，探索性地来决定这些合力绕实车1重心点G产生的横摆方向的力矩为最大时的驱动/制动力和横向力的一组，并将与该组对应的驱动/制动力和打滑比分别决定为Fxmmax_n、Smmax_n。 

另外，在S102的处理中，与就所述S100的处理进行说明的情况同样，也可以包括作为变量参数的第n轮Wn的实着地负荷Fzact_n。 

接着，如后面所述执行S104～S112的处理，决定目标第n轮驱动/制动力Fxcmd_n。在这种情况下，以满足下述条件(1)～(3)的方式来决定目标第n轮驱动/制动力Fxcmd_n。其中，条件(1)～(3)是以条件(1)、(2)、(3)的顺序作为优先顺序高低的条件，并且，当不能决定出满足全部条件(1)～(3)的目标第n轮驱动/制动力Fxcmd_n时，则以先满足优先顺序高的条件的方式来决定目标第n轮驱动/制动力Fxcmd_n。 

条件(1)：当FF综合目标第n轮驱动/制动力FFtotal_n和目标第n轮驱动/制动力Fxcmd_n为制动方向的驱动/制动力时，使得目标第n轮驱动/制动力Fxcmd_n的大小(绝对值)不小于FF综合目标第n轮驱动/制动力FFtotal_n的大小(绝对值)。换言之，不会出现0＞Fxcmd_n＞FFtotal_n的情况。 

条件(2)：当目标第n轮驱动/制动力Fxcmd_n变为与最大力矩产生时第n轮驱动/制动力Fxmmax_n相同极性时，使得Fxcmd_n的大小(绝对值)不会超过Fxmmax_n的大小(绝对值)。换言之，不会出现Fxcmd_n＞Fxmmax_n＞0、或者Fxcmd_n＜Fxmmax_n＜0的情况。 

条件(3)：使得目标第n轮驱动/制动力Fxcmd_n尽可能与第n轮驱动/制动力候补Fxcand_n一致(更加准确而言，是使得Fxcmd_n与Fxcand_n的差的绝对值为最小)。 

此处，条件(1)是为了使得目标第n轮驱动/制动力Fxcmd_n不会小 于实车1驾驶者操作刹车踏板所要求的实车1的第n轮Wn的制动方向的驱动/制动力(其相当于FFtotal_n)的条件。作为补充，在本说明书的实施方式中，因为后轮W3、W4是从动轮，故而与后轮W3、W4有关的FF综合目标第n轮驱动/制动力FFtotal_n(n＝3、4)及目标第n轮驱动/制动力Fxcmd_n(n＝3、4)一直为0以下的值。因此，关于后轮W3、W4，条件(1)与“使得目标第n轮驱动/制动力Fxcmd_n的大小(绝对值)不小于FF综合目标第n轮驱动/制动力FFtotal_n的大小(绝对值)。”的条件相同。 

另外，条件(2)是为了使得对应于目标第n轮驱动/制动力Fxcmd_n而在第n轮Wn产生的横向力不至于变得过小的条件。 

另外，条件(3)是为了尽可能满足通过所述执行装置动作FB目标值决定部20b和FF法则部22来决定的、执行机构3的动作的控制要求(目标)的条件。另外，Fxcand_n如上所述是：按照所述车轮特性关系(第n轮Wn的侧滑角为实侧滑角βf_act或βr_act而路面摩擦系数为推定摩擦系数μestm，并以此为前提条件时的车轮特性关系)，可能在第n轮Wn产生的驱动/制动力的值的范围内，与所述无限制目标第n轮驱动/制动力Fxdmd_n最接近(包括相一致的情况)的驱动/制动力的值。因此，条件(3)若换种说法，与下述的条件等同，即该条件为：目标第n轮驱动/制动力Fxcmd_n是按照所述车轮特性关系(以下述情况为前提条件时的车轮特性关系，即指第n轮Wn的侧滑角为实侧滑角βf_act或βr_act而路面摩擦系数为推定摩擦系数μestm)，可能在第n轮Wn产生的驱动/制动力的值的范围内的值；而且，尽可能与无限制目标第n轮驱动/制动力Fxdmd_n(依照控制要求的驱动/制动力)一致或接近(与Fxdmd_n之差的绝对值为最小)的条件。 

所述S104～S112的处理具体而言，按如下方式执行。首先，进入S104，判断在S100决定的Fxcand_n与在S102决定的Fxmmax_n的大小关系是否为0>Fxmmax_n>Fxcand_n或0<Fxmmax_n<Fxcand_n。当该判断结果为否时，进入S106，将Fxcand_n的值代入目标第n轮驱动/制动力Fxcmd_n。即，当Fxcand_n与Fxmmax_n为相互不同的极性时，或者当Fxcand_n与Fxmmax_n为相同极性，且Fxcand_n的大小(绝对值)在Fxmmax_n的大小(绝对值)以下时，将Fxcand_n的值原样地代入Fxcmd_n。另外，当 Fxcand_n＝0时(此时，Fxdmd_n也为0)，也将Fxcand_n的值代入Fxcmd_n(设Fxcmd_n＝0)。 

另一方面，当S104的判断结果为是时，进入S108，将Fxmmax_n的值(在S102决定的值)代入目标第n轮驱动/制动力Fxcmd_n。 

通过至此为止的处理，以满足所述条件(2)、(3)的方式(其中，优先条件(2))来决定Fxcmd_n。 

在S106或S108的处理后，进入S110，判断所述FF综合目标第n轮驱动/制动力FFtotal_n与当前的目标第n轮驱动/制动力Fxcmd_n(在S106或S108决定的值)间的大小关系是否为0>Fxcmd_n>FFtotal_n。当该判断结果为是时，进入S112，将FFtotal_n重新代入目标第n轮驱动/制动力Fxcmd_n。即，FF综合目标第n轮驱动/制动力FFtotal_n与在S106或S108决定的第n轮驱动/制动力候补Fxcmd_n为制动方向上的驱动/制动力，而且，在Fxcmd_n的大小(绝对值)小于FFtotal_n的大小(绝对值)时，将FFtotal_n的值代入Fxcmd_n。另外，当S110的判断结果为否时，原样维持此时的Fxcmd_n的值。 

通过以上的S104～S112的处理，如上所述，基本上是以满足所述条件(1)～(3)的方式来决定目标第n轮驱动/制动力Fxcmd_n。而且，当无法决定出满足全部条件(1)～(3)的目标第n轮驱动/制动力Fxcmd_n时，则以先满足优先顺序高的条件的方式来决定目标第n轮驱动/制动力Fxcmd_n。 

当S110的判断结果为否时，或者在S112的处理后，执行S114中的处理。在该S114，将与如上所述在S106～S112的处理中决定的Fxcmd_n相对应的打滑比决定为目标第n轮打滑比Scmd_n。在这种情况下，根据所述S104～S112的处理，Fxcmd_n是Fxcand_n、Fxmmax_n、FFtotal_n中的任意一值。而且，当Fxcmd_n＝Fxcand_n时，在S100求得的第n轮打滑比候补Scand_n被决定为Scmd_n；当Fxcmd_n＝Fxmmax_n时，在S102决定的最大力矩产生时第n轮打滑比Smmax_n被决定为Scmd_n。另外，当Fxcmd_n＝FFtotal_n时，例如只要依据在所述S100的处理中使用的图表，求解出与FFtotal_n对应的打滑比，并将所求得的打滑比决定为Scmd_n即可。这种情况下，当存在有2种与FFtotal_n对应的打滑比的值时，只要 将接近于0的一方的打滑比的值(第n轮Wn的驱动/制动力为峰值时的打滑比的值与0之间的范围内的值)决定为Scmd_n即可。另外，FFtotal_n在该图表中，当脱离出可能在第n轮Wn产生的驱动/制动力的值的范围时，在其范围内，将与最接近于FFtotal_n的驱动/制动力的值对应的打滑比决定为Scmd_n即可。 

以上是最佳目标第n驱动/制动力决定部241a_n(n＝1、2、3、4)的详细处理。 

另外，在本实施方式中，虽然是在决定目标第n轮驱动/制动力Fxcmd_n之后来决定与之对应的目标第n轮打滑比Scmd_n的，但也可以与之相反，在决定目标第n轮打滑比Scmd_n之后，再来决定与之对应的目标第n轮驱动/制动力Fxcmd_n。在这种情况下，依据与所述条件(1)～(3)对应的目标第n轮打滑比Scmd_n的相关条件，通过与所述S104～S112相同的处理来决定目标第n轮打滑比Scmd_n，随后，只要决定与该Scmd_n对应的Fxcmd_n即可。另外，在这种情况下，在依照第n轮Wn的所述车轮特性关系的打滑比与驱动/制动力之间关系中，在驱动/制动力为峰值时的打滑比的值与0之间的范围内决定Scmd_n。 

下面，说明所述最佳目标自动舵角决定部247的处理。图22是表示该最佳目标自动舵角决定部247的处理的功能方框图。 

参照该图，最佳目标自动舵角决定部247首先依据Fyfbdmd_f，通过处理部247a来决定FB自动舵角δf_fb，而该FB自动舵角δf_fb是指：为了使得实车1在前轮W1、W2产生由所述执行装置动作FB目标值决定部20b决定的自动操舵用FB目标横向力Fyfbdmd_f(具体而言，仅使前轮W1的横向力和前轮W2的横向力之合力变化Fyfbdmd_f)而被要求的前轮W1、W2的舵角的变化量。在这种情况下，在处理部247a，例如根据第一轮W1的实着地负荷Fzact_1，并通过规定的函数式或图表来求解出第一轮W_1的回转率Kf_1，而且同时根据第二轮W2的实着地负荷Fzact_2，并通过规定的函数式或图表，来求解出第二轮W2的回转率Kf_2。所述函数式或图表是依据实车1的前轮W1、W2的轮胎特性而预先设定的。并且，使用该回转率Kf_1、Kf_2，通过下式30，来决定FB自动舵角δf_fb。

δf_fb＝(1/(Kf_1+Kf_2))·Fyfbdmd_f  ......式30 

如此求得的FB自动舵角δf_fb相当于：为了使得前轮W1、W2的横向力之合力变化Fyfbdmd_f而所需的前轮侧滑角的修正量。 

另外，在通常情况下，因为实着地负荷Fzact_1、Fzact_2的变化较小，故而可以使在式子30中乘在Fyfbdmd_f上的系数(1/(Kf_1+Kf_2))为一定值。另外，也可以根据实着地负荷Fzact_1、Fzact_2和推定摩擦系数μestm来决定回转率Kf_1、Kf_2。 

接着，最佳目标自动舵角决定部247通过用加法器247b将如上所述决定的δf_fb加在FF目标前轮舵角δf_ff上，来决定目标前轮舵角δfcmd。 

另外，在不决定与所述状态量偏差γerr、βerr对应的自动操舵用FB目标横向力Fyfbdmd_f、或者一直使Fyfbdmd_f＝0时，只需将δf_ff原样地决定为目标前轮舵角δf_cmd即可。 

以上是所述执行装置动作目标值合成部24的详细处理。 

[关于执行装置驱动控制装置] 

所述执行装置驱动控制装置26为满足由所述执行装置动作目标值合成部24决定的目标值而使实车1的执行机构3动作。例如，以使得由驱动/制动装置3A的驱动***动作产生的第一轮W1的驱动/制动力(驱动方向的驱动/制动力)达到所述目标第一轮驱动***驱动/制动力的方式来决定该驱动***的执行装置操作量，并与之对应使该驱动***动作。此外，以使得第一轮W1的实路面反力之中的驱动/制动力(由驱动***动作产生的第一轮W1的驱动/制动力和由刹车装置动作产生的第一轮W1的驱动/制动力(制动方向的驱动/制动力)之和)达到所述目标第一轮驱动/制动力Fxcmd_1的方式来决定刹车装置的执行装置操作量，并与之对应使该刹车装置动作。而且，在这种情况下，以使第一轮W1的实打滑比和所述目标第一轮打滑比Scmd_1之差接近于0的方式来调整驱动***或刹车装置的动作。关于其他车轮W2～W4也与之相同。 

另外，当转向装置3B为执行装置驱动型的转向装置时，以使实前轮 舵角δf_act与所述目标前轮舵角δfcmd相一致的方式来决定转向装置3B的执行装置操作量，并与之对应来控制转向装置3B的动作。另外，当转向装置3B为执行装置辅助型的转向装置时，以使实前轮舵角δf_act与下述之和相一致的方式来控制转向装置3B的动作，该和为：所述目标前轮舵角δf_cmd与对应于转向角θh的机械性舵角量之和。 

另外，按照所述目标变速减速比，控制驱动/制动装置3A的驱动***的变速装置的减速比。 

此外，各车轮W1～W4的驱动/制动力或横向力等控制量容易引起驱动/制动装置3A、转向装置3B、悬架装置3C的动作的互相干涉。在这种情况下，为了将该控制量控制为目标值，最好是通过非干涉化处理综合性控制驱动/制动装置3A、转向装置3B、悬架装置3C的动作。 

[第二实施方式] 

接着，以下说明本发明的第二实施方式。而且，本实施方式与所述第一实施方式相比，不同的仅是用于调整模型车辆的过渡响应特性的特性调整用矩阵式K。因此以该不同部分为中心进行说明，关于与第一实施方式相同的构成和相同的处理，则省略其说明。 

在所述第一实施方式中，提出了以下的例子：即，为了提高模型车辆的状态量γd、βd的衰减性，将所述式子01的特性调整用矩阵K的对角成分即模型特性调整用参数k1、k2中的k2的值设定为对应于行驶速度Vact(＝Vd)可变，并将k1固定设为“1”，其中，所述模型车辆的状态量γd、βd对应于转向角θh的阶梯状变化。 

相对第一实施方式，在本实施方式中，模型特性调整用参数k1、k2中的k2的值被固定设为“1”。并且，为了提高模型车辆的状态量βd、γd的衰减性，在标准动态特性模型16的处理中，将k1的值设定为对应于行驶速度Vact可变，其中，所述模型车辆的状态量βd、γd对应于转向角θh的阶梯状变化。 

具体而言，例如与所述第一实施方式中的k2的值的设定的方式相同，将模型特性调整用参数k1的值设定为对应于行驶速度Vact可变。即，在Vact≤Vd_critical的情况下，设k1＝1。而在Vact>Vd_critical的情况下， 根据Vact设定k1的值，以使对应于转向角θh的阶梯状变化的模型车辆状态量βd、γd的响应特性成为临界制动的特性。在该情况下，Vact>Vd_critical时的k1的值是满足下式105的数值。 

(k1·a11+a22)²—4·k1·(a11·a22+a12·a21)＝0   ……式105 

如此设定的k1的值在Vact>Vd_critical的情况下，会被设定为比1大的值。此时，随着Vact的上升，k1的值呈单调增加。 

这样，通过将模型特性调整用参数k1的值设定为对应于行驶速度Vact可变，与第一实施方式同样，能够将模型车辆的稳定特性维持在接近于实车1的稳定特性的一种特性，同时使模型车辆的状态量βd、γd的过渡性的响应特性成为一种非振动性的特性，其中所述模型车辆的状态量βd、γd对应于转向角θh的阶梯状变化。 

除了以上所说明的构成以及处理，此外的构成和处理与第一实施方式中的相同。 

作为补充，也可以将Vact>Vd_critical情况下的k1的值设定为比下述情况下的k1值稍大的数值，所述情况是指：由式子01表达的模型车辆的状态量βd、γd的响应特性成为临界制动的一种特性的情况。即，也可以以使模型车辆的状态量βd、γd的响应特性成为过度制动的一种特性的方式来设定k1的值。另外，在行驶速度Vact为极低速(比临界制动速度Vd_critical小的规定值以下的行驶速度)的情况下，在模型车辆的响应性不成为振动特性的范围内，可将k1的值设为小于“1”。 

[第三实施方式] 

接着，说明本发明的第三实施方式。另外，本实施方式与所述第一实施方式或第二实施方式相比，不同的仅是用于调整模型车辆的过渡响应特性的特性调整用矩阵K。因此以该不同部分为中心进行说明，关于与第一实施方式或第二实施方式相同的构成和相同的处理，则省略其说明。 

在所述第一实施方式以及第二实施方式中，为了提高与转向角θh的阶梯状变化对应的模型车辆的状态量βd、γd的衰减性，仅将所述式子01的特性调整用矩阵K的模型特性调整用参数k1、k2中的一方设定为可变。在该情况下，当行驶速度Vact(＝Vd)比临界制动速度Vd_critical要高时，模型车辆的特性方程式det(λ·I-K·A)＝0的解(固有值)的绝对值有变得过大的倾向。进而，模型车辆的状态量γd、βd相对于转向角θh的变化的响应性(快速响应性)与实车1的相比会有变得过高的倾向。例如，如图23的虚线所例示的那样，上述第一实施方式中的模型车辆的特性方程式的解(固有值)的绝对值会对应于行驶速度Vact变化。如图所示，在行驶速度Vact比规定速度Vx(＞Vd_critical)还要高的高速区域，模型车辆的特性方程式的解(固有值)的绝对值伴随着行驶速度Vact的上升而增大。 

因此，特别是在增益值设得较小的情况下，在实车1的高速行驶时的过度期中状态量偏差γerr、βerr有时会变得较大，其中，所述增益值是假想外力Fvir、Mvir(模型操作用控制输入)相对于所述状态量偏差γerr、βerr的增益值。并且，在这种情况下，对应于状态量偏差γerr、βerr的执行机构3(驱动/制动装置3A、转向装置3B)的反馈控制所引起的执行装置的操作量恐怕会变得过大。 

因而，在第三实施方式中，将式子01中的特性调整用矩阵K的模型特性调整用参数k1、k2两者设定为可变。由此，在提高模型车辆的状态量βd、γd的衰减性的同时，防止了模型车辆的特性方程式det(λ·I-K·A)＝0的解、即模型车辆的固有值的绝对值变得过大。 

具体而言，在本实施方式中，模型特性调整用参数k1、k2两者的值根据行驶速度Vact而被设为可变。在这种情况下，k1、k2的值按下述方式被设定，该方式为：由式子01表述的模型车辆的状态量γd、βd的响应特性不成为振动性的特性(该响应特性成为临界制动或过度制动的特性)，且，该模型车辆的特性方程式的解的绝对值收敛在规定值以下。能满足这种条件的k1、k2的值的组，例如根据行驶速度Vact而预先被制成图表。并且，在各控制处理周期中的标准特性模型16的处理中，根据行驶速度Vact的值(此次值)，并基于所述图表来设定k1、k2的值。 

图24(a)、(b)是分别表示本实施方式中k1、k2的值对应于实行驶速度Vact的设定例的图表。 

在该例中，Vact≤Vd_critical时的k1、k2的值都被设定为1。而当Vact>Vd_critical时，如图24(b)所示，k2的值以与所述第一实施方式大致相同的方式，被设定为随着Vact的上升呈单调增加。另一方面，如图24(a)所示，k1的值在Vact上升到所述规定速度Vx之前被一直维持在1。并且，一旦Vact>Vx，则k1的值被设定为随着Vact的上升呈单调递减。这时，在图24的例子中，在Vact>Vx时，模型车辆的特性方程式det(λ·I—K·A)＝0的解(固有值)的绝对值与Vact＝Vx时该特性方程式的解的绝对值相等，且设定k1的值使得模型车辆的状态量γd、βd的响应特性成为临界制动的特性。 

通过以这样的方式来设定k1、k2的值，在行驶速度Vact的、Vact>Vx时的任意取值中，模型车辆的特性方程式的解(固有值)的绝对值成为Vact＝Vx时的该特性方程式的解的绝对值以下的数值。在本实施方式中，如图23中实线所示，在行驶速度Vact的、Vact>Vx时的任意取值中，模型车辆的特性方程式的解(固有值)的绝对值被维持在与Vact＝Vx时的该特性方程式的解的绝对值相等的数值。 

作为补充，例如本实施方式，在将特性调整用矩阵K所有的对角成分k1、k2的值设为可变的情况下，就可以任意设定模型车辆的特性方程式的解的值。 

在本实施方式中，通过将模型特性调整用参数k1、k2两者的值对应于行驶速度Vact而设定为可变，能够防止该模型车辆的状态量βd、γd的响应特性成为振动性的一种特性，并同时能够防止模型车辆的特性方程式的解的绝对值变得过大，进而防止模型车辆的状态量βd、γd相对于转向角θh变化的响应性(快速响应性)变得过大。另外，特性调整用矩阵K不会对模型车辆的稳定特性产生影响。因此，能够不受模型特性调整用参数k1、k2的值的变化的影响而将模型车辆的稳定特性维持在与实车1的稳定特性相接近的特性上，该情况与第一实施方式或第二实施方式相同。 

此外，在以上说明的实施方式中，以使用二阶***的动态特性模型的情况为例进行了说明，但是本发明也适用于使用三阶以上的动态特性模型的情况。 

例如，也可以使用相对于车轮侧滑角的变化考虑了横向力变化的响应滞迟的一种模型来作为标准动态特性模型。

这种情况下，标准动态特性模型的动态特性例如可以由下式110来表示。 

【数9】 

$\frac{d}{\mathrm{dt}}\left[\begin{array}{c}\mathrm{\&beta;d}\\ \mathrm{\&gamma;d}\\ \mathrm{Ffy}\_d\\ \mathrm{Fry}\_d\end{array}\right]=K_{4a}\&CenterDot;(A_{4a}\&CenterDot;\left[\begin{array}{c}\mathrm{\&beta;d}\\ \mathrm{\&gamma;d}\\ \mathrm{Ffy}\_d\\ \mathrm{Fry}\_d\end{array}\right]+\left[\begin{array}{c}0\\ 0\\ -\mathrm{kyf}\&CenterDot;\mathrm{Vd}\\ 0\end{array}\right]\&CenterDot;\mathrm{\&delta;f}\_d+\left[\begin{array}{c}0\\ 0\\ \mathrm{Fvir}1\\ \mathrm{Fvir}2\end{array}\right])$

                                      ……式子110 

其中， 

          $K_{4a}=\left[\begin{array}{cccc}k1& 0& 0& 0\\ 0& k2& 0& 0\\ 0& 0& k3& 0\\ 0& 0& 0& k4\end{array}\right]$

         $A_{4a}=\left[\begin{array}{cccc}0& -1& 2/(m\&CenterDot;\mathrm{Vd})& 2/(m\&CenterDot;\mathrm{Vd})\\ 0& 0& 2\&CenterDot;\mathrm{Lf}/I& -2\&CenterDot;\mathrm{Lf}/I\\ \mathrm{kyf}\&CenterDot;\mathrm{Vd}& \mathrm{kyf}\&CenterDot;\mathrm{Lf}& -\mathrm{kyf}\&CenterDot;\mathrm{Vd}/\mathrm{Kf}& 0\\ \mathrm{kyr}\&CenterDot;\mathrm{Vd}& \mathrm{kyr}\&CenterDot;\mathrm{Lr}& 0& -\mathrm{kyr}\&CenterDot;\mathrm{Vd}/\mathrm{Kr}\end{array}\right]$

另外，在式子110的条件部分，kyf、kyr分别为前轮的横向刚性、后轮的横向刚性。其中，m、Lf、Lr、I的含义与前述式子01中的相同。并且，Ffy_d、Fry_d分别是前轮的每一个轮的横向力、后轮的每一个轮的横向力。另外，Fvir1、Fvir2是一种假想外力(假想性的平移力)，作为模型操作量控制输入被附加性地作用于模型车辆上，用于使实车1与式子110的模型车辆之间的状态量偏差(横摆比率的偏差和车辆重心点侧滑角的偏差)接近于0。Fvir1、Fvir2是指假想性的横向平移力，分别作用在模型车辆的前轮位置和后轮位置。 

而且，与所述式子01的特性调整用矩阵K相同，式子110的K_4a是用于调整标准动态特性模型的动态特性的对角矩阵。通过将式子110的特性调整用矩阵K_4a的对角分量k1、k2、k3、k4中的一个以上的对角分量的 值设定为除了0以外的可变的数值，可以将式子110的模型的稳定特性维持为恒定，并同时可以使状态量βd、γd相对于转向角θh的阶梯状变化(前轮舵角δf_d的阶梯状变化)的衰减性比实车1的要高。 

另外，也可以使用诸如在绕车辆前后轴的动作与绕上下轴的动作之间具有干涉性的模型作为标准特性模型。这样的模型例如可以用下式111的形式来表示，其中，式子111是在所述非专利文献1中的式子(6.29)’、(6.30)’、(6.31)’上附加了模型操作用的假想外力(模型操作用控制输入)和特征调整用的参数之后的式子。 

d/dt(βd，γd，φd，φd′)^T＝K_4b·(A_4b·(βd，γd，φd，φd′)^T

             +B·δf_d+Bvir·(Fyvir，Mzvir，Mxvir)^T) 

                                     ……式子111 

另外，式子111中的φd是车身1B的侧翻角(绕前后方向轴的倾斜角)，φd′是φd的微分值(侧翻角的角速度)，A_4b是4行4列的正方矩阵，B是4行1列的矩阵，Bvir是4行3列的矩阵。并且，Fyvir、Mzvir、Mxvir是作为模型操作用控制输入的假想外力，该模型操作用控制输入用于使状态量偏差(例如车辆重心点侧滑角的偏差、横摆比率的偏差、侧翻角或其角速度的偏差)接近于0。在该情况下，Fyvir是横向假想性的平移力，Mzvir是绕横摆轴的假想性的力矩，Mxvir是绕前后方向轴的假想性的力矩。这些假想外力Fyvir、Mzvir、Mxvir根据该状态量偏差而决定。 

而且，式子111中的K_4b与所述式子01的特性调整用矩阵K相同，是用于调整标准动态特性模型的动态特性的对角矩阵(4行4列的对角矩阵)。将式子111的特性调整用矩阵K_4b的4个对角成分中一个以上的对角成分作为用于调整式子111的模型的响应特性的参数，通过将该参数的值设定为0以外的值、且可变，由此可以将式子111的模型的稳定特性维持为恒定，并同时可以使状态量βd、γd相对于转向角θh的阶梯状变化(前轮舵角δf_d的阶梯状变化)的衰减性比实车1的要高。 

此外，在所述第一～第三实施方式中，虽然例示了使用与车辆的侧滑相关的状态量以及与绕横摆轴的转动相关的状态量来作为控制对象的状 态量，也可以使用与车辆其他动作相关的状态量。例如，也可以使用与绕前后方向轴转动的动作相关的状态量(例如侧翻角和其角速度)。此时，在车辆的行驶速度较高的情况下，也可以提高侧翻角以及其角速度的衰减性。 

另外，在所述第一～第三实施方式中，虽然以4轮的车辆为例进行了说明，但本发明还能够适用于摩托车等车辆方面。 

工业实用性 

综上所述可知，本发明在能够提供以下汽车等车辆的方面起到作用。所述汽车等车辆能够不受车辆动作状态的影响，将对应了操舵角等的驾驶操作量变化的车辆响应特性控制在衰减性较高的一种特性，并显示出良好的动作特性。

Claims (5)

1.一种车辆控制装置，该车辆控制装置具有：

驾驶操作量测定机构，其测定驾驶操作量，该驾驶操作量表示由操纵者操纵具有多个车轮的车辆时的该车辆的驾驶操作状态；

执行机构，其设置于该车辆上，可操作所述车辆的规定的动作；

实际状态量掌握机构，其测定或推定实际状态量向量，该实际状态量向量是与实际车辆的规定的动作有关的多种状态量的值的组合；

模型状态量决定机构，其决定模型状态量向量，该模型状态量向量是作为表示所述车辆的动态特性的模型而被事先设定的车辆模型下的车辆的所述多种状态量的值的组合；以及

状态量偏差计算机构，其算出状态量偏差，该状态量偏差是所述测定或推定而得的实际状态量向量的各种状态量的值与所述被决定的模型状态量向量的各种状态量的值之间的偏差，

并且，该车辆控制装置至少控制所述执行机构的动作以使所述状态量偏差接近于0，该车辆控制装置的特征在于，其还具有：

状态量偏差应动控制机构，该状态量偏差应动控制机构至少根据所述算出的状态量偏差并以使该状态量偏差接近于0的形式来决定实车执行装置操作用控制输入和车辆模型操作用控制输入，其中，所述实车执行装置操作用控制输入用于操作实际车辆的所述执行机构，所述车辆模型操作用控制输入用于操作所述车辆模型下的车辆的动作；和

执行机构控制机构，该执行机构控制机构至少根据所述被决定的实车执行装置操作用控制输入来控制所述执行机构的动作，

其中，所述模型状态量决定机构是至少根据所述测定而得的驾驶操作量和所述被决定的车辆模型操作用控制输入来决定所述模型状态量向量的一种机构，

同时，该车辆控制装置还具有车辆模型特性设定机构，该车辆模型特性设定机构将所述车辆模型的至少一个参数的值设定为对应于所述实际车辆的动作状态而可变，以使所述模型状态量向量的各种状态量的值的衰减性成为比所述实际状态量向量的各种状态量的值的衰减性要高的一种高衰减特性，其中，所述模型状态量向量的各种状态量的值的衰减性对应于将所述车辆模型操作用控制输入维持在0的状态时的所述驾驶操作量的阶梯状变化，所述实际状态量向量的各种状态量的值的衰减性对应于将所述实车执行装置操作用控制输入维持在0的状态时的该驾驶操作量的阶梯状变化。

2.根据权利要求1所述的车辆控制装置，其特征在于，所述高衰减特性是临界制动或过度制动的响应特性。

3.根据权利要求1所述的车辆控制装置，其特征在于，所述多种状态量包含有与所述车辆的侧滑以及绕横摆轴的旋转有关的两种以上的状态量；所述车辆模型特性设定机构将所述车辆模型的参数的值设定为对应于所述实际车辆的行驶速度而可变。

4.根据权利要求3所述的车辆控制装置，其特征在于，所述车辆模型特性设定机构设定所述车辆模型的参数的值并使其对应于所述行驶速度为可变，以使至少在所述行驶速度比规定速度高的时候，所述车辆模型的固有值即模型车辆的特性方程式的解的绝对值成为规定值以下的数值。

5.根据权利要求1所述的车辆控制装置，其特征在于，所述车辆模型特性设定机构设定所述车辆模型的参数值并使其可变，以使下述驾驶操作量的恒定值与所述模型状态量向量的各种状态量的恒定值的关系保持为固定的关系，其中，所述驾驶操作量的恒定值是将所述车辆模型操作用控制输入维持在0的状态下使所述驾驶操作量产生阶梯状变化时的该驾驶操作量的恒定值。

Images