基于mc9s12dg128微控制器智能汽车设计

第一届“飞思卡尔”杯全国大学生智能汽车邀请赛技术报告学校华北电力大学（保定）队伍名称睿智参赛队员刘廉隅郭丽君汪志佳带队教师林永君关于技术报告和研究论文使用授权的说明本人完全了解第一届“飞思卡尔”杯全国大学生智能汽车邀请赛关保留、使用技术报告和研究论文的规定，即参赛作品著作权归参赛者本人，比赛组委会和飞思卡尔半导体公司可以在相关主页上收录并公开参赛作品的设计方案、技术报告以及参赛模型车的视频、图像资料，并将相关内容编纂收录在组委会出版论文集中。参赛队员签名带队教师签名日期基于MC9S12DG128微控制器智能汽车设计林永君，刘廉隅，郭丽君，汪志佳（华北电力大学自动化系，保定071003）文摘本设计采用飞思卡尔16位微控制器MC9S12DG128作为小车模型的核心控制单元，路面黑线检测使用反射式红外光电传感器，车速检测使用断续式光电开关，通过接收管接收光的频率时时计算小车速度，利用PWM技术动态控制舵机转角以及电动机的转速。采用PID控制思想，用单片机C语言编程，智能化控制使小车沿黑线行走。通过仿真以及实验检测，PID控制思想可以对小车进行较好的控制，使其快速、平稳的沿黑线行驶。关键词微控制器；传感器；智能控制；汽车模型；PIDINTELLIGENTCARDESIGNBASEDONMC9S12DG128MICROCONTROLLERLINYONGJUNLIULIANYUGUOLIJUNWANGZHIJIANORTHCHINAELECTRICUNIVERSITYAUTOMATIONDEPARTMENT,BAODING071003ABSTRACTTHISDESIGNUSESFREESCALE16BITMICROCONTROLLERMC9S12DG128TOBETHECARMODELSCORECONTROLUNITTHEDETECTINGOFTHEBLACKLINEONTHEROADSSURFACEUSEREFLECTIVEINFRAREDPHOTOELECTRICSENSORSTHEDETECTINGOFVEHICLESSPEEDUSESOFFANDONLIKEPHOTOELECTRICSWITCH,ANDCALCULATESTHECARSSPEEDEVERANDAGAINBYRECEIVINGLIGHTFREQUENCYWITHTHERECEIVINGCELLS,ANDUSESTHEPWMTECHNOLOGYDYNAMICCONTROLSACTUATORSCORNERASWELLASTHEROTATIONALSPEEDOFTHEMOTORUSESTHEIDEAOFTHEPIDCONTROL,WITHCPROGRAMMINGINTHESINGLECHIPMICROCOMPUTER,THEINTELLECTUALIZEDCONTROLCAUSESTHECARTORUNALONGTHEBLACKLINETHROUGHTHESIMULATIONASWELLASTHEEXPERIMENTEXAMINATION,THEPIDCONTROLMAYCARRYONABETTERCONTROLTOTHECAR,ANDMAKEITRUNALONGTHEBLACKLINEFASTANDSTEADILYKEYWORDMICROCONTROLLERSENSORINTELLIGENTCONTROLCARMODELPID目录第一章引言1第二章智能汽车方案的选择与论证221电动机驱动调速模块222舵机驱动模块223路面黑线探测模块224车轮测速模块325控制模块326电源4小结4第三章系统的具体设计与实现631系统的硬件设计6311电动机PWM驱动模块的电路设计与实现6312路面黑线检测模块的电路设计与实现6313车速检测模块的设计与实现7314舵机驱动模块的设计与实现7315电源的设计与实现832系统的机械设计8321车模结构的改进与调整8322路况检测传感器的设计制作与安装9323测速传感器的设计制作与安装9324系统电路板的固定1033系统的软件设计10331路面黑线检测103311黑线检测原理103312黑线检测算法11332车速检测12333控制过程133331舵机控制原理133332电机控制原理14第四章智能汽车控制流程16第五章结论17参考文献附录A程序代码第一章引言自从开始制作智能汽车以来，我们已经在系统规划论证、方案选择实施以及实验验证等方面做了大量的工作，智能汽车制作基本完成。这个技术报告对智能汽车制作的方案选择、系统的软硬件设计、控制理论论证以及整体控制流程做了比较详细的说明。本设计的特色在于将光敏三极管检测到的电压转换成连续的模拟量，在传感器数量和布局受到限制的条件下对黑线位置进行精确定位；用简单的PID控制算法有效控制小车。第二章智能汽车方案的选择与论证路面检测单片机电动机驱动转速检测舵机驱动图21系统可以划分为几个基本模块，如图21。对各模块的实现，分别有以下一些不同的设计方案。21电动机驱动调速模块方案一采用专用的PWM电机驱动芯片,通过PWM精确调节车速。本方案优点在于电路简单、调速精确，但芯片功率较小、不易散热、稳定性没有保障。方案二采用达林顿管组成的PWM驱动电路,通达控制达林顿管使之工作在占空比可调的开关状态,调整电机转速。由于管子工作在饱和状态,效率较高并且管子本身功率较大,稳定性能也很好。基于上述理论分析，拟选择方案二。22舵机驱动模块方案一单片机PWM信号直接控制舵机。方案二单片机PWM信号通过PWM驱动电路驱动舵机。由于舵机自身内部带有驱动电路可采用PWM信号直接驱动舵机，但是为了防止单片机的过载、所以采用方案二。23路面黑线探测模块路面信息检测的任务就是在白色路面上检测出黑线相对于小车的位置，以确定小车的行使方向。这种检测主要是对黑白两色区分,常用的方式有两种一种是用摄像头采集路面信息并将图像送至单片机进行分析，一种是根据黑白对光反射系数的不同,对路面进行光的发射与接收,根据接收到光的强弱来定位黑线的位置。方案一安装CCD摄像头方式。在小车车头顶部安置摄像头，对车前方大范围的路况进行连续拍摄，将图象送至单片机进行图象识别。这种方式可以提前将路况了解得十分完整，达到预先判断的效果，如果能够实现，一定是非常有效的检测方案。但是这个方案的难点在于，现在市场上CCD原件的分辨率大多为320240，就算是只采集黑白两种颜色，每次采集的数据量大约为320240/89600字节，也就是96K字节，对于一个只有8KRAM的16位单片机来说处理量过大，不能及时处理数据就意味着不能对小车方向做出及时有效的控制，因此这种方案无法达到比赛要求。方案二红外光发光二级管与光敏三级管组成的发射接收器,将光敏三级管电流变化转变成电压变化信号,再将电压变化信号变成高低电平送至单片机。由于采用了红外光管代替可见光管有效防止了外界光源的干扰,输出信号为开关量有利于单片机对信号的处理；但是每对传感器只能输出两个信息量,在传感器数量和布局受到限制的条件下很难对黑线位置进行精确定位。方案三红外光发光二级管与光敏三级管组成的发射接收器,将光敏三级管电流变化转变成电压变化信号,再将电压信号放大送至单片机A/D转换口。这种方式不但有方案三的优点,而且传感器输出的是连续的模拟量,在一定的范围内可根据电压与黑线位置的关系曲线定位黑线的位置,所以通过这种方式可在传感器数量和布局受到限制的条件下对黑线位置进行精确定位。通过比较，方案三易实现且能够达到较高的检测精度，这更有利于采用PID控制策略对小车进行有效控制，因此采用方案三。24车轮测速模块方案一采用霍尔集成片，此器件内部由三片金属组成，当磁铁正对金属板时，由于霍尔效果使金属板发生横向导通，因此可以在车轮上安装磁片，将霍尔集成片固定在车轮旁边，通过对脉冲的计数进行车速测量。方案二反射式光电开关。在车轮一侧贴上黑白相间的贴片，将由红外光发光二级管与光敏三级管组成的发射接收器对准贴片。在车轮转动时光可间断地反射到三级管上。通过对间隔周期的测量，计算出车速。以上两种都是比较可行的测速方法,但是由于车模较小，方案一的磁片安装十分困难，也容易产生相互干扰，相反方案二适用于精度较高的场合，所以拟采用方案二。25控制模块方案一采用模糊控制。模糊控制不需要建立精确的数学模型，但是需要根据经验建立合适的隶属函数。它是解决不确定系统的有效途径，并且具有较强的鲁棒性。方案二采用神经网络。神经网络控制同样不需要建立精确数学模型，神经网络在控制系统中的应用主要有两个方面一是采用神经网络作为被控系统的辨识器；另一是神经网络在控制系统直接充当控制器。这些都需要大量的试验来让神经网络在不同环境下进行学习。方案三采用PID进行控制。这种方法建立在被控对象精确数学模型基础上的。考虑到被控对象小车的数学模型比较简单，可以用以下方程表示舵机/1SDKVEIST电机QPID算法可以对此控制对象进行简单有效的控制。由于模糊控制、神经网络控制比较复杂，需要大量的试验，且不易编程实现。基于上述理论分析，拟选择方案三。26电源电源可分为三个部分电机电源、单片机以及光电传感器和LED电源、舵机电源。电机电源为提高控制效率，采用直接将电池电压通过电机驱动模块加在电机上。单片机电源为了保证提供稳定的5V电压，采用稳压器件LM2940将电池电压恒降至5V。舵机电源由于舵机电源电压范围在48V至6V之间,不能直接将电池电压加在舵机上。方案一将电池电压用稳压器件稳定到舵机电源电压范围之内,但是由于舵机内部的动力机构是直流电机,此电机工作时处于频繁的启停状态,会对稳压器件产生很大冲击,容易将稳压器件击穿。方案二通过二级管降压将电池电压降至舵机电源电压范围之内，这种方法电路简单、功耗小、性能稳定。通过比较采用方案二。小结经过一番比较论证，我们决定系统各个主要模块最终方案如下电动机驱动与调速模块采用达林顿管组成的PWM驱动电路。舵机驱动模块单片机PWM信号通过PWM驱动电路驱动舵机。路面黑线检测模块采用红外发射接收器，经信号放大、A/D转换、编程精确定位黑线位置。车速检测模块采用光电断续开关构成的光电感应系统。控制模块采用PID控制。第三章系统的具体设计与实现31系统的硬件设计311电动机PWM驱动模块的电路设计与实现入入VCC510120K10K10K15K10K10F电机PWM驱动电路如图311所示。本电路是一个基于PWM驱动单向转动电机的一个低功耗、低压降、高效率、高驱动能力的驱动电路。电路主要采用大功率达林顿管TIP42C作为驱动管，以保证电机启动瞬间的大电流要求。电路中用双管驱动是因为我们在实际运行中发现用单管驱动时赛车全速前进速度不够理想，经测试发现单管运行时压降达到24V；加了双管驱动后，压降仅为08V，小车全速前进速度比单管驱动时大大提高。电路中二级管的作用是起到在输入信号由高到低时，通过二级管形成电机感应电势的回路，起到保护电机和驱动管的作用。电机加旁路电容是为了使电机运行更加平稳，减小电池电压的波动。图311电机PWM驱动电路312路面黑线检测模块的电路设计与实现为了检测到路面上黑线的具体位置，在车前端安装了一排光电传感器,由8个红外光发光二级管和7个光敏三级管排列组成。为了精确定位黑线位置将三级管的输出信号转换放大成电压的模拟量送至单片机A/D转换口转换成数字量后通过程序算法给出黑线的精确位置。发射、接收和信号处理的具体电路如图312。红外发光二级管限流电阻采用了510是为了提高发射光的强度，削弱外部光源的干扰，提高传感器的抗干扰能力。入入72V10K5101K100F901351010K入入1K20K20K5V10K10K20K图312黑线检测模块电路由于环境条件的不同和跑道路面材料的不同会引起光电传感器输出信号大小的不同，为了保证输出信号范围恒定，在电压信号放大电路中、反馈电阻采用了电位器。313车速检测模块的设计与实现在差速器齿轮的一侧贴上一张用电脑绘制打印出来的黑白相间的纸片，如图313所示，将配对好的光电开关固定在纸片的正前方，齿轮带动纸片转动时光电开关接收管便得到间断的红外光信号、从而输出间断的高低电平分信号，信号经整形送至单片，以实现车速的检测。图313将纸片贴在差速器齿轮上而没有贴在车轮上是由于两后轮的轴并不是一体的，两轴间是通过差速器连接的、在转弯时两车轮转速不同。入入入入51K5V51KLM2940470F470F001F001F图314光电开关及其信号处理电路光电开关及其信号处理电路如图314。在纸片转动时光电传感器的输出电压变化在08V42V之间，为了提高稳定性电压比较器的比较电压取两个极限电压的平均值25V。由于常用运放的输出电压幅度都很难升到电源电压、也很难得很低，所以在运放输出端又加了只三级管增加输出电压的幅度以防止单片机对输入信号的误判。314舵机驱动模块的设计与实现舵机驱动模块电路如图315所示，由于舵机自身内部带有驱动模块,而且驱动信号的输入阻抗相当大,所以采用了一个简单，低功耗的PWM信号驱动电路。图315舵机驱动电路315电源的设计与实现电机电源为提高控制效率，采用直接将电池电压通过电机驱动模块加在电机上。入入72VPWM图316单片机以及光电传感器电源电路单片机以及光电传感器电源电源具体电路如图316所示，采用低压差线性稳压器件LM2940。该稳压器件输入输出压差低输出电流大（最大输出1A）输出电压稳定精度高。可为单片机提供一个稳定的5V电源。二级管的作用是防止稳压器件的反向击穿。舵机电源电源电路图如图317所示降压管采用硅二级管1N4007，电池在充满电后空载电压能达到8V而在满载时电压又降至7V左右，为了保证给舵机的电压保持在48V至6V之间，以每只硅管压降07V计算、加三只管刚好能满足要求。图317舵机电源电路32系统的机械设计车模的机械部分分为车模结构的改进与调整；路况检测传感器的制作与安装；测速传感器的制作与安装；系统电路板的固定。321车模结构的改进与调整影响车模运动性能的机械结构主要有前轮的转向传动机构、后轮的转动传动机构，所以我们对与这几个机构有关的各部件进行了相应的调整。在实验中发现在舵机自锁的情况下前轮有很大的转向空程，经过仔细的调试发现舵机自身空程很小，主要原因来自于从舵机转轴到转向拉杆的传动机构上的转动空程，经过对这个传动机构的改进完全消除了这部分空程。改进后在舵机自锁的情况下前轮转向空程变得很小，大大提高了车模的转向精度。后轮转动传动机构的调整主要是对电机主轴齿轮与差速器齿轮配合的调整和对差速器预紧力的调整。在齿轮配合中若配合过紧则转动阻力拒太大影响车速、若配合过松则容易发生滑齿，在调整中力求使两齿轮的齿合线相重合以达到最佳配合效果。对于差速器的预紧力如果过大会减弱差速器的差速效果，在转弯时两后轮速度差不够会增加车模前进阻力，如果预紧力过小则在车模加速时差速器容易发生打滑。调整预紧力时，在保证最大加速度时差速器不发生打滑前提下尽量减小。322路况检测传感器的设计制作与安装图321传感器结构图传感器结构图如图321所示，由于路况信息的检测靠7路传感器输出的电压信号作为计算依据，所以各传感器的输也特性曲线必须相同，这就要求各管的间距和离地面的距离要相同。为了达到这一要求我们采用了传感器定们板，在定们板上开上15个直径48MM、间距762MM的孔，将发射管和接收管按顺序插入孔中使各管露出板的距离相等然后再将定位板用螺丝固定在电路板上，调整两边螺丝使定们板与电路板平行，最后再将各管焊接在电路板上。在车模前端底盘上反向固定有三颗340的螺丝，传感器靠螺母固定在三颗螺丝上，调节固定螺母在一平面内以使电路板与底盘平行，同时调节几个螺母的上下位置改变传感器距离地面的距离使传感器信号输出范围最大。323测速传感器的设计制作与安装图322测速传感器结构图测速传感器的结构图如图322所示，光电开关采用经配对的光电管焊接在电路板上，为了提向反射光强度两管中心线成30角。光电开关用铝合金支架固定在差速器齿轮正前方。差速器齿轮上贴有黑白相间的纸片，如图313所示，纸片分为24个区域12个黑色区域，12个白色区域，也就是说车轮每转动一周传感器输出12个周期的脉冲。324系统电路板的固定系统电路板靠螺丝直接固定在车模中后部的底盘上方，S12电路板直接插接在系统电路板上。经改造后，小车车身长30CM，宽16CM，高10CM，重量18KG。33系统的软件设计单片机MC9S12DG128是控制核心，实现对路面黑线的软件检测，小车方向控制，车速检测，电动机驱动。331路面黑线检测3311黑线检测原理黑线检测共用8个红外发光管，7个红外接收管，将7个红外接收管编号为（3，2，1，0，1，2，3）即中间管为0号管。采用MC9S12DG128自带的A/D转换功能，通过设置相应寄存器，可以将7个接收管接收的光强信号由模拟量转换成相应的数字量。当接收管位于黑线正中时，输出电压最高约45伏，转换成数字量大约为720，远离黑线的接收管输出电压最低约为11伏，转换成数字量大约为330。图331单个接收管的感光曲线经实验得到单个接收管的感光曲线如图331所示。因为采用PID控制策略，所以黑线位置反馈信号越精确，越有利于PID控制，控制效果也会越好。其中当接收管距黑线较远时，曲线在最小值处，输出的数字量大约在230300之间变化，且变化很缓慢。当接收管渐渐靠近黑线时，输出数字量变化很快，在300700之间变化，且曲线斜率几乎保持不变，因此可认为管子输出电压在300到700之间与移动的距离成线性关系，将其视为一段直线其斜率根据所测得数据计算可得，K700300/14。当接收管靠近黑线中心4MM以内，输出在700720之间变化，曲线又趋平缓，且有波动。黑线检测流程如图332。图332黑线检测流程图程序初始化后，将接收管程序开始初始化接收管数据A/D转换转换数据排序选出次值和次次值计算黑线位置数据A/D转换、排序，判断是否有大于600的数，600是经验数值，如果没有大于600的数，说明小车没在黑线上，那么继续A/D转换数据、排序，直到有大于600的数；当有大于600的数时，说明小车头传感器在黑线上，我们可以根据测得各个接收管A/D转换后的数值，通过下面的算法可以判断出黑线位置。图333拟合后感光曲线图3312黑线检测算法经过MATLAB处理，可以将曲线线性区域拟合为一条直线（如图333中直线），与曲线基本重合。A/D转换程序每转换一次，我们都可以得到一组7个240720的数，找到7个管子中电压值最大的管子，并设最大值为MAX，记录其接收管序号。再取其左右两个值，其一应为次大值MAX_1,另一个为次次大值MAX_2。如图334所示。图334最大的数值所对应的接收管最靠近黑线中心，且黑线应该在输出为最大值的接收管和输出为次大值的接收管之间。因为两个接收管输出值MAX_1、,MAX_2均在300道700之间，所以次大值和次次大值肯定都能落在斜率为K的这段直线上。我们把次大值和次次大值代入公式，即可得此接受管距离中670KYX心线的距离。用接收管编号乘上15便可得到此接收管距离0号接收管的距离，正号表示在0号管右边，负号表示在0号管左边。当输出为次小值的接收管在0号管右边时，小车中心偏离黑线的距离应该为此管距0号管距离加上它距黑线的距离。当它在0号管左边时，小车中心偏离黑线的距离应该为此管距0号管距离减去它距黑线的距离。按此法计算所得值为正说明黑线在小车左边，为负说明小车在黑线右边。即以小车中心位置，零号管的位置为坐标原点，当黑线在小车左边时输出为负值，在右边时为正值。由于只用次大值可能造成误差，我们这里对次次大值用同样的方法计算小车偏离中线距离，然后两者相加取均值。由于七个接收管性能相似，其最大值、最小值、线性区以及曲线走势极为相近，所以上面用中间接收管的算法也可以用于其它六个接收管。用MATLAB将三个接收管特性比较，其曲线如图335图335将一组实测数据（235、253、571、686、309、243、227），代入程序测得小车中心偏离黑线119MM，取整数后为12MM。实际黑线在小车左边12MM处。其误差范围在2MM。332车速检测（如图336）程序开始读取捕捉到的累计脉冲数脉冲数除以时间返回图336车速检测当车轮转动时，安装在车轴上的圆形折光板上的圆孔依次通过断续式光电开关，每通过一个圆孔产生一个脉冲。因此只需记录固定时间内捕捉到的脉冲数，然后除以时间，便能得到实际车速。其中，时间由定时器的时钟计算得到。333控制过程3331舵机控制原理采用PID控制，按照7个接收管建立一维坐标系，0号管为坐标原点，面向车前进方向，左边为负，右边为正，黑线位置为5050。为控制小车能在路中间行驶，转弯时不偏离跑道，将“0”作为给定值，实时测得的黑线位置作为反馈值。当小车偏离跑道时，调整舵机转角，使反馈值跟踪给定值，让小车一直沿黑线行驶。1被控对象小车舵机，根据经验对系统建模型，将其简化为积分，比例，一阶惯性加迟延环节。由于传动机构非常灵敏，所以迟延时间，和惯性时间都很短，根据经验设参数，迟延时间常数01S，惯性时间常数02S。T综上所述，其传递函数为SES10212调节器采用比例控制。经整定其参数可设置为，比例系数，3PK积分时间常数S。0T舵机控制原理方框图如图337PID012SES图337舵机控制原理方框图用MATLAB仿真，阶跃相应曲线如图338。图338控制舵机阶跃相应曲线由此图可以看出，小车舵机可以由比例环节P控制住。3332电机控制原理电机转速与小车速度成比例关系。因为小车在直线应当加速，在转弯处需要减速，小车舵机转角越大，应该减速越多。所以由舵机的12个转角度数（左右两侧个6个转角）给出电机6个档位（电机6个档位转速由实验与经验给定），根据舵机转角改变电机转速给定值，测量电机实际转速作为电机转速实际值，令实际值跟踪给定值，控制小车速度。采用反馈闭环控制1被控对象为小车电机，仍根据经验对系统建模型，将其简化为比例，一阶惯性加迟延环节。同舵机一样其传动机构非常灵敏，所以迟延时间，和惯性时间都很短，通过实验测得迟延时间常数01S，惯性时间常数02S。综T上所述其传递函数为。SES1022调节器采用比例，积分。经整定其参数可设置为，比例系数，积分时间常数S。1PK5T控制原理方框图如图339PID012SE图339电机控制原理方框图图3310电机阶跃相应曲线加阶跃相应，用MATLAB仿真其控制效果如图3310。由此曲线可以看出电机可以由PI控制得很好。由于小车舵机、电机模型是根据经验给出，而且实际小车舵机、电机有一定的耦合关系，控制器KP、T等参数都需要根据实验来整定、完善。第四章智能汽车控制流程智能汽车的控制流程如图41程序开始初始化接收管数据A/D转换转换数据排序判断黑线位置舵机PWM程序电机转速舵机转角舵机PID控制电机PWM程序根据舵机转角给出电机转速档位电机PID控制图41智能汽车控制流程图首先开始程序，初始化程序，将7个光敏接收管接收到的电压信号经过A/D转换为7个数字信号，对这7个数据进行由大到小排序，判断最大值是否有大于600（经验数值）的数，如果没有，那么通过PWM程序将电机转速置于0档，接收管继续采数；如果最大值大于600，则通过黑线位置判断程序计算黑线位置。将黑线位置给舵机PID控制模块后，继续判断是否有大于600的数，这样循环刷新黑线位置。然后，因为舵机PID给定值定为0，黑线位置（5050）为实际值，当实际值不为0时，经PWM程序控制舵机转角，直到实际值变为0，此时第0号管位于黑线中心。最后，根据舵机转角给出电机转速的档位，作为电机PID控制的给定值，电机PWM程序控制电机转速，使电机转速实际值跟踪给定值。整个控制小车的程序在这样的循环往复的进行中，控制小车沿黑线前进。第五章结论经过理论论证、仿真研究与实验验证，本控制方案可以使小车很好的沿黑线行驶。小车行驶最大速度约为15M/S。但是由于采用连续的数值（而非开关量）对小车舵机控制，黑线位置的检测、计算受外界环境尤其是光线强度影响较大，黑线位置的求取算法、程序还需要进一步改善。小车的模型参数、控制参数都还需要大量的实验，来继续完善对小车的控制，优化小车的速度。参考文献1刘会英，杨志强机械原理北京机械工业出版社，20032韩利竹，王华MATLAB电子仿真与应用2版北京国防工业出版社，200393于希宁，刘红军自动控制原理北京中国电力出版社，20014姚佩阳，曹锦，常永昌，张科英，李锋自动控制原理北京清华大学出版社，20055周斌，蒋荻南，黄开胜基于虚拟仪器技术的智能车仿真系统2006,02S1321346邵贝贝单片机嵌入式应用的再现开发方法北京清华大学出版社，2004107马忠梅，籍顺心，张凯，马岩单片机的C语言应用程序设计3版北京北京航空航天大学出版社,2003118徐叙瑢,苏勉曾发光学与发光材料北京化学工业出版社20049苏大图光学测量北京机械工业出版社1988610翟瑞彩，谢伟松数值分析M天津天津大学出版社，200011HOLLANDJHADAPTATIONINNATURALANDARTIFICIALSYSTEMSANINTRODUCTORYANALYSISWITHAPPLICATIONSTOBIOLOGY，CONTROLANDARTIFICIALINTELLIGENCEARMHARBORUNIVERSITYOFMICHIGANPRESS，1975附录A程序代码/MYDEFINEH/UNSIGNEDINTTIMER