基于软计算的疲劳驾驶融合检测方法

1.一种基于软计算的疲劳驾驶融合检测方法，其特征在于，包括如下步骤：
步骤1被测驾驶员疲劳特征参数xi的计算：对被测驾驶员的疲劳特征参数xi进行采集，i＝1、2、3或4，所采集的疲劳特征参数包括：驾驶员面部图像特征信号、车辆前方车道标志线信号、方向盘转动信号，其中，驾驶员面部图像特征信号包括驾驶员眨眼信号及打哈欠信号，并根据所采集的数据，计算眨眼频率x1、打哈欠频率x2、车辆异常偏离车道的频率x3及方向盘转动异常频率x4；
步骤2将眨眼频率x1、打哈欠频率x2、车辆异常偏离车道频率x3及方向盘转动异常频率x4代入高斯隶属函数 i＝1，2，3或4，j为正整数，j＝1，
2，3，…，(cR-1)或cR，cR为模糊规则个数，式中，ωij、σij和 分别表示疲劳特征参数xi的第j个模糊规则的隶属函数的中心、宽度和隶属度；
步骤3利用步骤2得到的隶属度 计算激励强度αj，即： i＝1，2，3或4，
再计算归一化激励强度 即： 然后计算第j个模糊规则的后件加权值yj，
即： 其中 为权值系数，d＝0，1，2，3或4，再计算归一化后件加
权值 即： 最后计算疲劳度检测值 即：
得到所述cR、ωij、σij和 的方法是：首先对6名男性和4名女性共10名驾驶员进行疲劳度量化实验，选择200组包含疲劳特征参数x1，x2，x3，x4及相应的疲劳度量化值Y的数据作为实验样本，为保证样本的时效性和完备性，疲劳度量化值Y的量化实验应在疲劳特征参数采集之后1分钟内开始，同时样本应尽量全面地包含驾驶员的各种疲劳程度，所述疲劳程度为不疲劳、轻微疲劳、中等疲劳或重度疲劳，其中，x1＝n1/N1，N1为Δt1时间内采集的图片帧数，n1为其中眼睛闭合的图片帧数，x2＝n2/N2，N2为Δt2时间内采集的图片帧数，n2为其中打呵欠的图片帧数，x3＝n3/N3，N3为Δt3时间内采集的图片帧数，n3为其中异常偏离车道的图片帧数，x4＝t/10，t为每10秒内方向盘持续不动的时间，单位为秒，疲劳度量化值 δb为对第b次弹出的交通标志图片识别错误的当量系数，
tb为疲劳度量化时对第b次弹出的交通标识图片进行识别和响应所用的时间，单位为秒，上述疲劳度量化值Y的量化采用编制的测试程序进行计算，实验前先让驾驶员对4张配有文字说明的交通标志图片进行学习，编程时每张图片被用一个图片标识符进行定义，与其相对应的文字说明按钮也被定义一个跟图片标识符相同的文字说明按钮标识符；实验时，程序界面每隔2秒随机弹出4张交通标志图片中的任一张，此时驾驶员须在限定的1秒时间内对图片所代表的含义进行识别和响应，并用鼠标迅速点击相应的文字说明按钮，根据每次出现的图片的图片标识符跟鼠标点击的文字说明按钮标识符是否相同来判断每次识别是否正确，当出现的图片的图片标识符跟鼠标点击的文字说明按钮标识相同时，此次识别正确，否则识别错误，同时计算从图片出现到鼠标击中文字说明按钮之间的响应时间，用tb表示，根据识别结果是否正确以及响应时间的长短确定每次识别错误的当量系数δb；
每次实验图片共随机出现100次，实验结束时程序会自动计算本次实验图片识别的错误率并将其作为疲劳程度的量化值在程序界面上进行显示；
然后利用减法聚类确定模糊规则个数cR，以及ωij和σij的初始值，减法聚类的步骤为：
(1)计算第k个样本Xk中xi和Y的归一化值 和 (j′＝1，
2，3或4)， 其中，对每一个样本值Xk，Xk＝{xk，1，xk，2，xk，3，xk，4，Yk}，k＝1，
2，…，200，前4项为第k个样本对应的四个疲劳特征参数值x1，x2，x3和x4，第5项为相应的疲劳度量化值Y，归一化样本值 经比较确定：
(2)计算每一个 对应的密度值Pk， Ra＝0.4，然后经统计比
较得到具有最大密度值 的数据点 定义 并初始化聚类个数cN＝0，
(3)令
(4)如果Rp＞εA，执行(6)，其中，接收率εA＝0.8，
(5)如果Rp＜εR，输出聚类中心Ci″及聚类个数cN，并将聚类个数作为模糊规则个数cR的值，即cR＝cN，程序结束，其中，拒绝率εR＝0.15，否则，统计计算 跟已求得的聚类中心Ci″之间的最小距离dmin， i″为正整数，1≤i″≤cN，C表示由聚
类中心组成的聚类中心集合，C＝{Ci″}，Ci″表示第i″个聚类中心，如果 拒绝成为聚类中心，并令其密度值Pk＝0，重新进行统计比较，选择具有最大密度值 的数据点作为新的 转至(3)，否则转至(6)，
(6)令 将 作为新的聚类中心依次添加到集合C中，并令cN＝cN+1，接着根据
对其余未聚类的 的密度值Pk进行更新，即 其中，Rb＝
1.5Ra，以避免出现距离相近的聚类中心，
(7)统计并找到新的具有最高密度值 的 返回步骤(3)继续运行程序，
减法聚类完成后，将聚类中心集合C中的cR个聚类中心Cj进行转换求得归一化前的聚类中心C′j，其中，
再将 作为隶属函数 中相应ωij的初始值 即
在集合C中统计计算跟第一个聚类中心C′1的欧氏距离最近的另一个聚类中
心C′s，s＝2，3，…，(cR-1)或cR，σij的初始值 被确定为： 其中，
为归一化前的聚类中心
C′j中的前4个元素， 为聚类中心Cj中的前4个元素， 为距第一个聚类中心C′1的欧氏距离最近的另一个聚类中心C′s中的前4个元素；
根据遗传算法，将TS模糊神经网络需要优化的ωij、σij和 按实数编码格式连接起来构成一条染色体，共随机产生100条染色体，组成群体Z，Z＝{z1，z2，…，z100}，其中，每条染色体的基因值由ω′ij、σ′ij和 组成，第h条染色体zh可表示为：
h为正整数，且h＝1，2，…，或100，其中，ω′ij、σ′ij和
为组成每条染色体基因值的元素，
计算网络在训练集上产生的均方误差ferr(Z)， 式中，
为代入第k个样本的疲劳特征参数xi和染色体中的ω′ij，σ′ij和 后TS模糊神经网络计算出的网络实际输出值，并根据ferr(Z)构造适应度函数，计算每条染色体的适应度ffit(Z)， 通过遗传算法的选择，交叉，变异操作对染色体中的ω′ij，
*
σ′ij和 进行优化，将在最大进化代数范围内满足最小均方误差 的染色体z 作为最佳染色体，
*
根据最佳染色体z，令 并将
ωij，σij和 作为最佳网络参数，其中，遗传算法中设定的最小均方误差 最大进化代数为Ne＝1000次，
选择操作采用轮盘赌法选择优质个体，交叉运算时利用公式(1)以线性交叉的方式产生子代
为待交叉的两个染色体，η为(0，1)区间产生的随机数，e为交叉运算的代数，变异操作采用非均匀变异，根据公式(2)和(3)，对原有的基因值做一随机扰动，以扰动后的结果做为变异后的新基因值d′(Zm)，
式中κ＝gc/gm，γ为随机数，γ＝rand(0，1)，g为形状参数，g＝3，al和ar分别为染色体中各基因位取值区域的左右边界、gc为当前进化代数、gm为最大进化代数，sign为随机数，取值为0或1，Zm为第m个染色体，
在交叉和变异操作过程中，交叉率pc和变异率pm采用一种自适应的方法获得，其表达式为：
式中fmax、favg分别表示群体中的最大适应度和平均适应度，f′表示两个待交叉个体中较大的适应度，pc1、pc2、pm1、pm2为小于1的常数，其中，交叉率pc1＝0.9，pc2＝0.6，变异率pm1＝0.1，pm2＝0.001；
步骤4将最佳网络参数ωij、σij和 以及实时计算的疲劳特征参数xi代入到TS模糊神经网络中计算疲劳度检测值 依据疲劳度检测值 的大小，将驾驶员的驾驶行为划分为四种状态： 最后根据 的大小检测驾驶员是否疲劳并对疲
劳的严重程度进行有效分类。

基于软计算的疲劳驾驶融合检测方法\n技术领域\n[0001] 本发明涉及一种驾驶员疲劳驾驶融合检测方法，尤其涉及一种基于软计算的疲劳驾驶融合检测方法，它是利用TS模糊神经网络融合两个面部疲劳特征和两个车辆行为特征，并运用减法聚类和遗传算法对网络进行优化和训练，进而对驾驶员疲劳驾驶进行融合检测的方法，属于驾驶员疲劳驾驶检测技术领域。\n背景技术\n[0002] 疲劳驾驶检测和预警已成为汽车主动安全领域中的研究热点。其中，基于物理传感器的非接触式疲劳驾驶检测方法在理论研究和应用领域近年来引起广泛重视。但现有的疲劳驾驶检测方法大都只针对驾驶员某一方面的单一疲劳特征，如频繁的眨眼和打呵欠，头部转动异常，或车辆行驶状态异常等。随着信息融合技术的发展，不少检测方法虽然已开始考虑融合多个疲劳特征，但它们大多只融合驾驶员面部的几个疲劳特征，而对于间接反应驾驶员是否疲劳的车辆行为特征，如因疲劳驾驶而产生的车辆异常偏离车道，方向盘转动和车速变化异常等却被忽略，导致目前这些方法的检测效果不够理想，容易造成误检，漏检等。同时，疲劳驾驶是一个非常复杂的生理现象，存在诱发原因多，症状复杂，难以检测和识别等问题，这些都给传统的疲劳驾驶检测方法带来了极大地挑战。\n发明内容\n[0003] 本发明的目的是克服现有基于单一特征疲劳驾驶检测方法的不足，提供一种可靠性高的基于软计算的疲劳驾驶融合检测方法。\n[0004] 本发明采用如下技术方案：\n[0005] 一种基于软计算的疲劳驾驶融合检测方法，该方法的主要步骤如下：\n[0006] 步骤1被测驾驶员疲劳特征参数xi的计算：对被测驾驶员的疲劳特征参数xi进行采集，i＝1、2、3或4，所采集的疲劳特征参数包括：驾驶员面部图像特征信号、车辆前方车道标志线信号、方向盘转动信号，其中，驾驶员面部图像特征信号包括驾驶员眨眼信号及打哈欠信号，并根据所采集的数据，计算眨眼频率x1、打哈欠频率x2、车辆异常偏离车道的频率x3及方向盘转动异常频率x4，其中，x1＝n1/N1，N1为Δt1时间内采集的图片帧数，n1为其中眼睛闭合的图片帧数，x2＝n2/N2，N2为Δt2时间内采集的图片帧数，n2为其中打呵欠的图片帧数，x3＝n3/N3，N3为Δt3时间内采集的图片帧数，n3为其中异常偏离车道的图片帧数，x4＝t/10，t为每10秒内方向盘持续不动的时间，单位为秒；\n[0007] 步骤2将眨眼频率x1、打哈欠频率x2、车辆异常偏离车道频率x3及方向盘转动异常频率x4代入TS模糊神经网络的高斯隶属函数μij中， i＝1，\n2，3或4，j为正整数，j＝1，2，3，…，(cR-1)或cR，cR为模糊规则个数，式中，ωij、σij和j\nμi 分别表示疲劳特征参数xi的第j个模糊规则的隶属函数的中心、宽度和隶属度；\n[0008] 步骤3利用步骤2得到的隶属度μij计算激励强度αj，即： i＝1，2，\n3或4，再计算归一化激励强度 即： 然后计算第j个模糊规则的后件加权\n值yi，即： 其中ρdj为权值系数，d＝0，1，2，3或4，再计算归一\n化后件加权值 即： 最后计算疲劳度检测值 即：\n[0009] 得到所述cR、ωij、σij和ρdj的方法是：首先对6名男性和4名女性共10名驾驶员进行疲劳度量化实验，选择200组包含疲劳特征参数x1，x2，x3，x4及相应的疲劳度量化值Y的数据作为实验样本，为保证样本的时效性和完备性，疲劳度量化值Y的量化实验应在疲劳特征参数采集之后1分钟内开始，同时样本应尽量全面地包含驾驶员的各种疲劳程度，如不疲劳、轻微疲劳、中等疲劳，重度疲劳，其中，x1＝n1/N1，N1为Δt1时间内采集的图片帧数，n1为其中眼睛闭合的图片帧数，x2＝n2/N2，N2为Δt2时间内采集的图片帧数，n2为其中打呵欠的图片帧数，x3＝n3/N3，N3为Δt3时间内采集的图片帧数，n3为其中异常偏离车道的图片帧数，x4＝t/10，t为每10秒内方向盘持续不动的时间，单位为秒，疲劳度量化值δb为对第b次弹出的交通标志图片识别错误的当量系数，\n[0010] \n[0011] tb为疲劳度量化时对第b次弹出的交通标识图片进行识别和响应所用的时间，单位为秒，上述疲劳度量化值Y的量化采用编制的测试程序进行计算，实验前先让驾驶员对\n4张配有文字说明的交通标志图片进行学习，编程时每张图片被用一个图片标识符进行定义，与其相对应的文字说明按钮也被定义一个跟图片标识符相同的文字说明按钮标识符；\n实验时，程序界面每隔2秒随机弹出4张交通标志图片中的任一张，此时驾驶员须在限定的1秒时间内对图片所代表的含义进行识别和响应，并用鼠标迅速点击相应的文字说明按钮，根据每次出现的图片的图片标识符跟鼠标点击的文字说明按钮标识符是否相同来判断每次识别是否正确，当出现的图片的图片标识符跟鼠标点击的文字说明按钮标识相同时，此次识别正确，否则识别错误，同时计算从图片出现到鼠标击中文字说明按钮之间的响应时间，用tb表示，根据识别结果是否正确以及响应时间的长短确定每次识别错误的当量系数δb；每次实验图片共随机出现100次，实验结束时程序会自动计算本次实验图片识别的错误率并将其作为疲劳程度的量化值在程序界面上进行显示；\n[0012] 然后利用减法聚类确定模糊规则个数cR，以及ωij和σij的初始值，减法聚类的步骤为：\n[0013] (1)计算第k个样本Xk中xi和Y的归一化值 和 (j′＝\n1，2，3或4)， 其中，对每一个样本值Xk，Xk＝{xk，1，xk，2，xk，3，xk，4，Yk}，k＝1，\n2，…，200，前4项为第k个样本对应的四个疲劳特征参数值x1，x2，x3和x4，第5项为相应的疲劳度量化值Y，归一化样本值 经比较确定：\n[0014] (2)计算每一个 对应的密度值Pk， Pa＝0.4，然后经\nmax\n统计比较得到具有最大密度值Pk 的数据点 定义 并初始化聚类个数cN＝0；\n[0015] (3)令\n[0016] (4)如果Rp＞εA，执行(6)，其中，接收率εA＝0.8；\n[0017] (5)如果Rp＜εR，输出聚类中心Ci″及聚类个数cN，并将聚类个数作为模糊规则个数cR的值，即cR＝cN，程序结束，其中，拒绝率εR＝0.15；否则，统计计算 跟已求得的聚类中心Ci″之间的最小距离dmin， i″为正整数，1≤i″≤cN，C表示由\n聚类中心组成的聚类中心集合，C＝{Ci″}，Ci″表示第i″个聚类中心，如果 拒绝 成为聚类中心，并令其密度值Pk＝0，重新进行统计比较，选择具有最大密度值Pkmax的数据点作为新的 转至(3)，否则转至(6)；\n[0018] (6)令 将 作为新的聚类中心依次添加到集合C中，并令cN＝cN+1，接着根据 对其余未聚类的 的密度值Pk进行更新，即 其中，\nRb＝1.5Ra，以避免出现距离相近的聚类中心；\n[0019] (7)统计并找到新的具有最高密度值Pkmax的 返回步骤(3)继续运行程序；\n[0020] 减法聚类完成后，将聚类中心集合C中的cR个聚类中心Cj进行转换求得归一化前的聚类中心Cj′，其中，\n* j\n再将xi，j 作为隶属函数μi 中相应ωij的初始值 即\n在集合C中统计计算跟第一个聚类中心C1′的欧氏距离最近的另一个聚类中\n心Cs′，s＝2，3，…，(cR-1)或cR，σij的初始值 被确定为： 其中，\n[0021] 根据遗传算法，将TS模糊神经网络需要优化的ωij、σij和ρdj按实数编码格式连接起来构成一条染色体，共随机产生100条染色体，组成群体Z，Z＝{z1，z2，…，z100}，其j\n中，每条染色体的基因值由ωij′、σij′和ρd′ 组成，第h条染色体zh可表示为：\n[0022] \n[0023] \n[0024] \nh为正 整数，且h＝1，2，…，或100，其中，\n计算网络在训练集上产\n生的均方误差ferr(Z)，\n[0025] 式中， 为代入第k个样本的疲劳特征参数xi和染色体中的ωij′，σij′和j\nρd′ 后TS模糊神经网络计算出的网络实际输出值，并根据ferr(Z)构造适应度函数，计算每条染色体的适应度ffit(Z)， 通过遗传算法的选择，交叉，变异操作对\nj\n染色体中的ωij′，σij′和ρd′ 进行优化，将在最大进化代数范围内满足最小均方误差*\n的染色体z 作为最佳染色体，\n[0026] \n[0027] \n[0028] \n*\n根据最佳染色体z，令 并将\nj\nωij，σij和ρd 作为最佳网络参数，其中，遗传算法中设定的最小均方误差 最大进化代数为Ne＝1000次；\n[0029] 选择操作采用轮盘赌法选择优质个体，交叉运算时利用公式(1)以线性交叉的方e+1 e+1\n式产生子代ZA ，ZB ，\n[0030] \n[0031] ZAe、ZBe为待交叉的两个染色体，η为(0，1)区间产生的随机数，e为交叉运算的代数，\n[0032] 变异操作采用非均匀变异，根据公式(2)和(3)，对原有的基因值做一随机扰动，以扰动后的结果做为变异后的新基因值d′(Zm)，\n[0033] \n[0034] \n[0035] 式中κ＝gc/gm，γ为随机数，γ＝rand(0，1)，g为形状参数，g＝3，al和ar分别为染色体中各基因位取值区域的左右边界、gc为当前进化代数、gm为最大进化代数，sign为随机数，取值为0或1，Zm为第m个染色体，\n[0036] 在交叉和变异操作过程中，交叉率pc和变异率pm采用一种自适应的方法获得，其表达式为：\n[0037] \n[0038] \n[0039] 式中fmax、favg分别表示群体中的最大适应度和平均适应度，f′表示两个待交叉个体中较大的适应度，pc1、pc2、pm1、pm2为小于1的常数，其中，交叉率pc1＝0.9，pc2＝0.6，变异率pm1＝0.1，pm2＝0.001；\n[0040] 步骤4将最佳网络参数ωij、σij和ρdj以及实时计算的疲劳特征参数xi代入到TS模糊神经网络中计算疲劳度检测值 依据疲劳度检测值 的大小，将驾驶员的驾驶行为划分为四种状态： 最后根据 的大小检测驾驶员是否疲劳\n并对疲劳的严重程度进行有效分类。\n[0041] 与现有技术相比，本发明具有如下优点：\n[0042] (1)本发明综合考虑反映疲劳驾驶的两个直接特征和两个间接特征，其中，两个直接特征包括频繁眨眼和打呵欠，两个间接特征包括车辆异常偏离车道和方向盘转动异常。\n上述四个疲劳特征能够全面、客观、准确的衡量驾驶员的疲劳驾驶行为。\n[0043] (2)运用TS模糊神经网络融合上述四个疲劳疲劳特征参数对驾驶员疲劳驾驶进行实时检测，克服基于单一特征疲劳驾驶检测方法的局限性，降低其误检率和漏检率，提高疲劳驾驶检测的可靠性和准确率，适合于实时性强的驾驶疲劳精确检测场合。\n[0044] (3)利用减法聚类确定TS模糊神经网络的结构和初始网络参数，包括TS模糊神经网络的模糊规则个数以及隶属函数的中心和宽度的初始值，从而优化了TS模糊神经网络的结构，简化了计算过程。\n[0045] (4)采用遗传算法训练TS模糊神经网络的初始网络参数，进而确定网络参数的最佳值。避免了传统方法训练TS模糊神经网络参数时容易陷入局部极小值的弊端，使得TS模糊神经网络的收敛性和检测精度均得到显著提高。\n附图说明\n[0046] 图1是本发明方法的检测流程图；\n[0047] 图2是疲劳度的量化流程图；\n[0048] 图3是基于TS模糊神经网络的结构图；\n[0049] 图4是减法聚类的流程图；\n[0050] 图5是遗传算法的流程图；\n[0051] 图6是基于软计算的模糊神经网络训练效果图；\n[0052] 图7是疲劳特征参数计算及驾驶行为分类输出框图。\n具体实施方式\n[0053] 下面参照说明书附图，对本发明的具体实施方式做出更为详细地说明：\n[0054] 本发明基于软计算的疲劳驾驶融合检测方法可通过两个阶段来实现，一是离线训练阶段，二是在线检测阶段。该方法的检测流程图如图1所示。\n[0055] 1、通过对数据样本进行离线训练确定最佳网络结构及网络参数\n[0056] 最佳网络结构及网络参数可通过六步完成：\n[0057] (1)数据采集\n[0058] 在实验车上安装两个CCD摄像头C1和C2，C1负责驾驶员面部图像特征信号的采集，C2负责车辆前方车道标志线信号的采集，另外在方向盘上安装光电式转角传感器负责对方向盘转动信号的采集。\n[0059] (2)特征参数计算\n[0060] 对被测驾驶员的疲劳特征参数xi进行采集，i＝1、2、3或4，所采集的疲劳特征参数包括：驾驶员面部图像特征信号、车辆前方车道标志线信号、方向盘转动信号，其中，驾驶员面部图像特征信号包括驾驶员眨眼信号及打哈欠信号，并根据所采集的数据，计算眨眼频率x1、打哈欠频率x2、车辆异常偏离车道的频率x3及方向盘转动异常频率x4，其中，x1＝n1/N1，N1为Δt1时间内采集的图片帧数，n1为其中眼睛闭合的图片帧数，x2＝n2/N2，N2为Δt2时间内采集的图片帧数，n2为其中打呵欠的图片帧数，x3＝n3/N3，N3为Δt3时间内采集的图片帧数，n3为其中异常偏离车道的图片帧数，x4＝t/10，t为每10秒内方向盘持续不动的时间，单位为秒；\n[0061] (3)疲劳度的量化\n[0062] 首先对6名男性和4名女性共10名驾驶员进行疲劳度量化实验，选择200组包含疲劳特征参数x1，x2，x3，x4及相应的疲劳度量化值Y的数据作为实验样本，为保证样本的时效性和完备性，疲劳度量化值Y的量化实验应在疲劳特征参数采集之后1分钟内开始，同时样本应尽量全面地包含驾驶员的各种疲劳程度，如不疲劳、轻微疲劳、中等疲劳，重度疲劳，其中，x1＝n1/N1，N1为Δt1时间内采集的图片帧数，n1为其中眼睛闭合的图片帧数，x2＝n2/N2，N2为Δt2时间内采集的图片帧数，n2为其中打呵欠的图片帧数，x3＝n3/N3，N3为Δt3时间内采集的图片帧数，n3为其中异常偏离车道的图片帧数，x4＝t/10，t为每10秒内方向盘持续不动的时间，单位为秒，疲劳度量化值 δb为对第b次弹出的交通标志图片识别错误的当量系数，\n[0063] \n[0064] tb为疲劳度量化时对第b次弹出的交通标识图片进行识别和响应所用的时间，单位为秒，上述疲劳度量化值Y的量化采用编制的测试程序进行计算，实验前先让驾驶员对\n4张配有文字说明的交通标志图片进行学习，编程时每张图片被用一个图片标识符进行定义，与其相对应的文字说明按钮也被定义一个跟图片标识符相同的文字说明按钮标识符；\n实验时，程序界面每隔2秒随机弹出4张交通标志图片中的任一张，此时驾驶员须在限定的1秒时间内对图片所代表的含义进行识别和响应，并用鼠标迅速点击相应的文字说明按钮，根据每次出现的图片的图片标识符跟鼠标点击的文字说明按钮标识符是否相同来判断每次识别是否正确，当出现的图片的图片标识符跟鼠标点击的文字说明按钮标识相同时，此次识别正确，否则识别错误，同时计算从图片出现到鼠标击中文字说明按钮之间的响应时间，用tb表示，根据识别结果是否正确以及响应时间的长短确定每次识别错误的当量系数δb；每次实验图片共随机出现100次，实验结束时程序会自动计算本次实验图片识别的错误率并将其作为疲劳程度的量化值在程序界面上进行显示，疲劳度的量化流程如图2所示。\n[0065] (4)构建TS模糊神经网络模型\n[0066] 本文采用TS模糊神经网络模型对疲劳特征进行融合检测，该模型由前件网络和后件网络组成，具有收敛速度快，样本需要量少的优点，其网络结构图如图3所示。\n[0067] 前件网络由四层组成，各层的功能及运算如下：\n[0068] 1)输入层：节点数为4，对应4个输入疲劳特征参数x1、x2、x3和x4；\n[0069] 2)模糊化层：对每一个节点利用高斯隶属函数计算其隶属度μij，其运算为：\ni＝1，2，3或4，j为正整数，j＝1，2，3，…，(cR-1)\n[0070] 或cR，cR为模糊规则个数，式中，ωij、σij分别表示疲劳特征参数xi的第j个模糊规则的隶属函数的中心、宽度；\n[0071] 3)规则层：该层用来计算每条规则的激励强度αj，\n[0072] 4)归一化层：该层用来计算规则的归一化激励强度\n[0073] 以下是后件网络各层的功能与计算：\n[0074] 1)输入层：比前件网络的输入层多一个节点输入x0＝1，用来补偿模糊规则后件中的常数项；\n[0075] 2)函数层：计算每一条规则的后件加权值yj，节点间的权值系数为后件参数。该层的输出为输入的线性组合，即：\n[0076] \nj\n[0077] 其中ρd 为权值系数，d＝0，1，2，3或4；\n[0078] 3)结合层：用来归一化后件加权值\n[0079] 4)输出层：该层有唯一节点，它将输入信号求和，计算疲劳度检测值 即：\n[0080] (5)减法聚类确定网络结构\n[0081] 利用减法聚类确定模糊规则个数cR，以及ωij和σij的初始值，减法聚类的步骤为：\n[0082] (1)计算第k个样本Xk中xi和Y的归一化值 和 (j′＝\n1，2，3或4)， 其中，对每一个样本值Xk，Xk＝{xk，1，xk，2，xk，3，xk，4，Yk}，k＝1，\n2，…，200，前4项为第k个样本对应的四个疲劳特征参数值x1，x2，x3和x4，第5项为相应的疲劳度量化值Y，归一化样本值 经比较确定：\n[0083] 2)计算每一个 对应的密度值Pk， Ra＝0.4，然后经统\n计比较得到具有最大密度值Pkmax的数据点 定义 并初始化聚类个数cN＝0；\n[0084] 3)令\n[0085] 4)如果Rp＞εA，执行6)，其中，接收率εA＝0.8；\n[0086] 5)如果Rp＜εR，输出聚类中心Ci″及聚类个数cN，并将聚类个数作为模糊规则个数cR的值，即cR＝cN，程序结束，其中，拒绝率εR＝0.15；否则，统计计算 跟已求得的聚类中心Ci″之间的最小距离dmin， i″为正整数，1≤i″≤cN，C表示由\n聚类中心组成的聚类中心集合，C＝{Ci″}，Ci″表示第i″个聚类中心，如果 拒绝 成为聚类中心，并令其密度值Pk＝0，\n[0087] 重新进行统计比较，选择具有最大密度值Pkmax的数据点作为新的 转至3)，否则转至6)；\n[0088] 6)令 将 作为新的聚类中心依次添加到集合C中，并令cN＝cN+1，接着根据 对其余未聚类的 的密度值Pk进行更新，即 其中，\nRb＝1.5Ra，以避免出现距离相近的聚类中心；\n[0089] 7)统计并找到新的具有最高密度值Pkmax的 返回步骤3)继续运行程序；\n减法聚类的流程如图4所示，减法聚类完成后，将聚类中心集合C中的cR个聚类中心Cj进行转换求得归一化前的聚类中心Cj′，其中，\n再将xi，j*作为隶属函数\nj\nμi 中相应ωij的初始值 即 在集合C中统计计算跟第一个聚类中心C1′的欧氏距离最近的另一个聚类中心Cs′，s＝2，3，…，(cR-1)或cR，σij的初始值 被确定为：\n其中，\n[0090] (6)遗传算法优化网络参数\n[0091] 根据遗传算法，将TS模糊神经网络需要优化的ωij、σij和ρdj按实数编码格式连接起来构成一条染色体，共随机产生100条染色体，组成群体Z，Z＝{z1，z2，…，z100}，其j\n中，每条染色体的基因值由ωij′、σij′和ρd′ 组成，第h条染色体zh可表示为：\n[0092] \n[0093] \n[0094] \nh为 正整数，且h＝1，2， …，或100，其中，\n[0095] 按如图5所示的遗传算法计算流程计算网络在训练集上产生的均方误差ferr(Z)， 式中， 为代入第k个样本的疲劳特征\nj\n参数xi和染色体中的ωij′，σij′和ρd′ 后TS模糊神经网络计算出的网络实际输出值，并根据ferr(Z)构造适应度函数，计算每条染色体的适应度ffit(Z)，通过遗传算法的选择，交叉，变异操作对染色体中的ωij′，σij′\nj\n和ρd′ 进行优化，将在最大进化代数范围内满足最小均方误差 的染色体\nz* 作 为 最 佳 染 色 体，\n*\n根据最佳染色体z，令\nj\n从而确定ωij，σij和ρd 并将其作为最佳网络参数；\n[0096] 选择操作采用轮盘赌法选择优质个体，交叉运算时利用公式(1)以线性交叉的方e+1 e+1\n式产生子代ZA ，ZB ，\n[0097] \n[0098] ZAe、ZBe为待交叉的两个染色体，η为(0，1)区间产生的随机数，e为交叉运算的代数，\n[0099] 变异操作采用非均匀变异，根据公式(2)和(3)，对原有的基因值做一随机扰动，以扰动后的结果做为变异后的新基因值d′(Zm)，\n[0100] \n[0101] \n[0102] 式中κ＝gc/gm，γ为随机数，γ＝rand(0，1)，g为形状参数，g＝3，al和ar分别为染色体中各基因位取值区域的左右边界、gc为当前进化代数、gm为最大进化代数，sign为随机数，取值为0或1，Zm为第m个染色体，\n[0103] 在交叉和变异操作过程中，交叉率pc和变异率pm采用一种自适应的方法获得，其表达式为：\n[0104] \n[0105] \n[0106] 式中fmax、favg分别表示群体中的最大适应度和平均适应度，f′表示两个待交叉个体中较大的适应度，pc1、pc2、pm1、pm2为小于1的常数，其中，交叉率pc1＝0.9，pc2＝0.6，变异率pm1＝0.1，pm2＝0.001；\n[0107] 选择200个实验样本，利用减法聚类算法对其进行聚类，最终确定的模糊规则个数cR＝5，遗传算法中设定的最小均方误差 最大进化代数Ne＝1000次，网络\n训练过程中获得的均方误差曲线如图6所示；\n[0108] 2、将离线训练阶段获得的最佳网络参数ωij，σij和ρdj代入TS模糊神经网络中，根据实时采集和计算的疲劳特征参数对驾驶员是否疲劳驾驶进行在线检测。\n[0109] (1)将离线训练阶段获得的最佳网络参数ωij，σij和ρdj代入TS模糊神经网络中；\n[0110] (2)数据采集\n[0111] 在线检测阶段的数据采集跟离线训练阶段的数据采集方式相同；\n[0112] (3)特征参数计算\n[0113] 在线检测阶段的特征参数计算跟离线训练阶段的特征参数计算方法相同；\n[0114] (4)疲劳度分类输出\n[0115] 将最佳网络参数ωij、σij和ρdj以及实时计算的疲劳特征参数xi代入到TS模糊神经网络中计算疲劳度检测值 依据疲劳度检测值 的大小，将驾驶员的驾驶行为划分为四种状态： 最后根据 的大小检测驾驶员是否疲劳并对疲\n劳的严重程度进行有效分类，疲劳特征参数计算及驾驶行为分类输出框图如图7所示。