基于模糊逻辑的车速控制系统设计

基于模糊逻辑的车速控制设计摘要：汽车巡航控制系统具有强非线性、时变不确定性，并受到外界扰动、复杂的运行工况等影响，采用传统PID控制很难取得满意的效果。本文设计了一种基于模糊PI和模糊PD控制算法的汽车巡航控制系统，并给出了系统的仿真结果。研究结果显示，所设计的巡航控制系统算法可行有效，性能指标均达到预定工程实用要求。关键词：巡航控制；模糊控制；PI控制；PD控制一、概述汽车巡航控制系统，简称COS(CruiseControlSystem)。它实际上是一种辅助驾驶系统。在装备巡航控制系统的汽车上，当汽车行驶速度超过一定值(一般为40km／h)时，如果驾驶员利用巡航控制开关设定一个车速，那么在巡航控制期间，随着道路坡度的变化以及汽车行驶中所可能遇到的阻力，车辆自动变换节气门开度或自动进行档位转换，以按存储在微机内的最佳燃料经济性规律或动力性规律稳定行驶。汽车巡航控制系统可以减轻驾驶员的负担，减少了不必要的车速变化，最大限度地节省燃料，降低排气污染，提高发动机的使用效率，并可在一定程度上提高汽车的动力性能和乘坐的舒适性。目前，很多车辆，特别是高级轿车已经把巡航控制系统作为配备设备或备选设备。由于国内对汽车巡航的研究起步较晚，并且技术相对落后，所以国内对汽车巡航控制系统的研究应用主要是以单车定速控制为主。虽然国内的一些规构已经开始了对电子式巡航控制装置的研究，但从总体上来说，目前国内对汽车巡航控制系统的研究还不是很成熟。目前，发动机和车辆系统的电子控制系统在不断增加，这有助于改善燃油经济性，减少排放，提高驾驶安全性以及降低生产成本。但是，车辆系统环境恶劣：特别是在高温、潮湿、振动强烈、充满电子干扰的发动机室，会导致电子系统各种形式的失效，进而甚至导致交通事故和人员伤亡。本文应用模糊控制逻辑设计了车速控制器，调节车辆速度到设定值并保持恒定。二、车辆模型构建根据研究角度选择汽车动力学模型如式(1)所示：式(1)式中：为控制输入（>0代表油门输入，<0代表刹车输入）；车辆质量m为1300Kg；空气阻力系数=0.3；恒定摩擦力=100N；为驱动/制动力；=0.2s。应用模糊控制方法的目标为快速调节车辆速度达到驾驶员设定值并保持恒定。三、基于模糊逻辑的速度控制器设计及仿真试验模糊控制与神经网络是控制工程师为了应对十分复杂的系统而提出的控制方法论，模糊控制在巡航控制领域有着优良的应用。3.1模糊PI控制器设计假如我们想精确的跟踪驾驶员设定的阶跃或倾斜车速变化，将应用图1所示的PI模糊车速控制器。图中，g0、g1、g2为信号增益，b(t)为积分器输入。图1模糊PI车速控制系统框图由于车辆系统为一阶常微分方程，因此可以用龙格库塔方法进行闭环仿真试验。实施以下不同参考输入下的仿真分析：1)当0<t<10时，初始车速设定为18m/s；当10<t<30时，目标车速设定为22m/s2)当0<t<10时，初始车速设定为18m/s；当t>10，目标车速在25s内由18m/s至22m/s线性倾斜变化；当25<t<30时，目标车速保持22m/s。图7测试输入二响应曲线图8测试输入三响应曲线可以看出，所设计的模糊PD控制器各项指标达到预定设计要求。在测试输入一响应曲线中，上升时间为9秒，在7至10秒的范围内；超调量远小于2%；稳定时间成功控制在10秒以内。对于测试输入二，稳态误差控制在1mps内。在测试输入三仿真试验中，上升时间符合7至10秒范围，稳态误差小，性能优良。总结与展望在车辆巡航控制系统中，为了保持恒定的车速，设计了一个基于模糊逻辑的控制器。通过仿真结果可以看出，当采用模糊逻辑控制方法联合经典PID控制算法来设计汽车巡航控制系统时，系统的超调小，反应速度快，具有良好的鲁棒性，可以达到预期控制效果。控制器设计过程可总结为：控制对象建模；模糊控制器设计；计算机数值仿真试验。参考文献[1]PassinoKM,YurkovichS.FuzzyControl,1stedn,AddisionWesleyLongman,Colifornia,1997[2]WardD.Berlitzcompleteguidetocruisingandcruiseships2000.Princeton,NewJersey:BerlitzPublishingCompany,1999[3]IoannouPA,ChienCC.”AutonomousIntelligentCruiseControl,”IEEETrans.onVehicularTechnology,1993,42(4):657~672[4]MayrR.”Intelligentcruisecontrolforvehiclebasedonfeedbacklinearization”.Proc.ofAmericanControlConference,1994:16~20[5]MayrR,BauerO.”Safetyissuesinintelligentcruisecontrol”.Proc.Of1999IEEEIntelligentTransportationSystems,1999:970~975.附录一文中仿真所用程序clear%VehicleParameterm=1300;Ar=0.3;tau=0.2;d=100;%Initialiazeparameterforthefuzzycontrollernume=11;numie=11;g1=1;g2=.01;g0=1000;we=0.2*(1/g1);wie=0.2*(1/g2);base=0.4*g0;ce=[-1-0.8-0.6-0.4-0.200.20.40.60.81]*(1/g1);cie=[-1-0.8-0.6-0.4-0.200.20.40.60.81]*(1/g2);rules=[-1-1-1-1-1-1-0.8-0.6-0.4-0.20;-1-1-1-1-1-0.8-0.6-0.4-0.200.2;-1-1-1-1-0.8-0.6-0.4-0.200.20.4;-1-1-1-0.8-0.6-0.4-0.200.20.40.6;-1-1-0.8-0.6-0.4-0.200.20.40.60.8;-1-0.8-0.6-0.4-0.200.20.40.60.81;-0.8-0.6-0.4-0.200.20.40.60.811;-0.6-0.4-0.200.20.40.60.8111;-0.4-0.200.20.40.60.81111;-0.200.20.40.60.811111;00.20.40.60.8111111]*g0;t=0;index=1;tstop=60;step=0.01;x=[18;197.2;20];whilet<=tstopv(index)=x(1);%Test1ift<=10,vd(index)=18;endift>10,vd(index)=22;end%Test2%ift<=10vd(index)=18;end%ift>10,vd(index)=vd(index-1)+(2/750);end%ift>25,vd(index)=22;end%Test3%vd(index)=22;ie_count=0;e_count=0;e(index)=vd(index)-v(index);b(index)=x(3);%%ife(index)<=ce(1)mfe=[10000000000];e_count=e_count+1;e_int=1;elseife(index)>=ce(nume)mfe=[00000000001];e_count=e_count+1;e_int=nume;elsefori=1:numeife(index)<=ce(i)mfe(i)=max([0,1+(e(index)-ce(i))/we]);ifmfe(i)~=0e_count=e_count+1;e_int=i;endelsemfe(i)=max([0,1+(ce(i)-e(index))/we]);ifmfe(i)~=0e_count=e_count+1;e_int=i;endendendend%%ifb(index)<=cie(1)mfie=[10000000000];ie_count=ie_count+1;ie_int=1;elseifb(index)>=cie(numie)mfie=[00000000001];ie_count=ie_count+1;ie_int=numie;elsefori=1:numieifb(index)<=cie(i)mfie(i)=max([0,1+(b(index)-cie(i))/wie]);ifmfie(i)~=0ie_count=ie_count+1;ie_int=i;endelsemfie(i)=max([0,1+(cie(i)-b(index))/wie]);ifmfie~=0ie_count=ie_count+1;ie_int=i;endendendend%%num=0;den=0;fork=(e_int-e_count+1):e_intforl=(ie_int-ie_count+1):ie_intprem=min([mfe(k)mfie(l)]);num=num+rules(k,l)*base*(prem-(prem)^2/2);den=den+base*(prem-(prem)^2/2);endendu(index)=num/den;time(index)=t;F=[(1/m)*(-Ar*x(1)^2-d+x(2));(1/tau)*(-x(2)+u(index));vd(index)-x(1)];k1=step*F;xnew=x+k1/2;F=[(1/m)*(-Ar*xnew(1)^2-d+xnew(2));(1/tau)*(-xnew(2)+u(