搜索
轮式移动机器人运动学建模方法*(二号黑体)常勇1,2马书根1,3王洪光1谈大龙1(四号仿宋)(1.中国科学院沈阳自动化研究所机器人学国家重点试验室沈阳110016；(五号宋体，文中数字和字母为TimesNewRoman)2.中国科学院研究生院北京100039；3.日本立命馆大学机器人系滋贺525-8577日本)摘要(小五黑体)：研究轮式移动机器人运动学建模方法问题。提出用于不规则地形下的轮式移动机器人运动学建模方法——轮心建模法(Wheel-centermodeling,WCM)在分析多刚体链式结构运动的速度特性以及不规则地形上轮式移动机构转动角速度特性的基础上，建立车轮轮心速度的矢量表达式，即确立WCM。以六轮摇臂转向架月球漫游车的运动学建模为例，分析各关节转动角速度矢量在车体坐标系下的投影，并根据转动矢量方向在车体坐标系中变化与否，将矢量叉乘的投影写成不同的形式，利用车轮轮心坐标系、轮地接触坐标系相对于车体中心坐标系的齐次变换矩阵中的相应元素，将车轮轮心的速度矢量表达式投影到车体中心坐标系下，建立车体运动学模型，从中阐述WCM的方法和过程，并分析WCM的特点。用试验和仿真验证该建模方法的正确性。(小五宋体，摘要长度300~500字)关键词(小五黑体)：双层辉光人工神经网络预测模型(小五宋体)中图分类号(小五黑体)：TG156(小五TimesNewRoman)ResearchonANN-basedPredictionModelUsedtoDoubleGlowPlasmaSurfaceAlloyingProcessing(小三)CHANG(姓大写)Yong1,2XXXXxxxxx1,3XXXXxxxxx1XXXXxxxxx1(小四)(1.StateKeyLaboratoryofRobotics,ShenyangInstituteofAutomation,ChineseAcademyofSciences,Shenyang110016;2.GraduateSchool,ChineseAcademyofSciences,Beijing100039;3.DepartmentofRobotics,RitsumeikanUniversity,Shiga-ken525-8577,Japan)(五号)Abstract(小五黑体)：Thekinematicmodelingmethodforawheeledmobilerobotisstudied.Amodelingmethodcalledwheel-centermodeling(WCM)forkinematicsofawheeledmobilerobotthatmovesonuneventerrainispresented.Thevectorexpressionofwheelcentervelocityisderivedbasedontheanalysisofthevelocitycharacteristicofamultirigidchain’slocomotionandthewheel-modelmovingmechanism’smovingonuneventerrain.ThemethodofWCMisintroducedasaresult.Theprojectionoftheangularvelocityoftherotaryjointsintherovercentreframeisanalysedandtheprojectionofvectorcrossproductisexpressedasdifferentformsaccordingtowhetherthedirectionoftheangularvelocityrelativelytorovercentreframechangesornot.Thewheelcentrevelocityisprojectedtotheroverbodycentrecoordinateframewiththeuseoftheelementsofthehomogeneoustransformationsfromthewheelcentreframeandthewheel-terraincontactframetotherovercenterframe.Thus,thekinematicsoftheroverismodelled.TheprocedureofkinematicmodellingusingWCMisintroducedandthecharactersofWCMisanalysed.Themethod’svalidityisverifiedbyexperimentandsimulation.(小五)Keywords(小五黑体)：DoubleglowArtificialneuralnetworkPredictionmodel0前言(四号宋体)(五号宋体)*离散频谱分析只能在截断的有限国家自然科学基金资助项目(7150080050)。xxxxxxxx收到初稿，xxxxxxxx收到修改稿（六号宋体，此处为角注，和正文分开）长度时域上进行，不可避免地存在能量泄漏，导致经过快速傅里叶变换(FastFouriertransform，FFT)得到的频谱，其频率、幅值和相位均可能产生较大的误差，大大限制了该技术的工程应用，因此需要研究离散频谱的校正技术以消除或大幅度减小这个误差，提高分析精度。对于单频率成分或间隔较远67的多频率成分的离散频谱，目前主要有比值法、能量重心法、FFT+FT连续细化傅里叶变换分析法和相位差法四种校正方法[1]，其中能量重心法是2001年丁康等[2]在三点卷积幅值校正法基础上提出来的一种通用的离散频谱校正方法，利用离散窗谱函数的能量重心是坐标原点的原理求出频率校正量，是一种校正精度很高的近似方法。在实际工程应用中所测得的信号通常伴随有噪声，噪声会影响频谱校正精度。因此从理论上系统的分析随机噪声背景下的各种离散频谱校正方法的精度，可为实际工程应用中提供理论依据，具有重要的指导意义。1992年SCHOUKENS等[3-4]研究了噪声对内插法的影响。2004年齐国清等[5]对插值FFT估计的频率精度进行了分析。2002年朱利民等[6]对一种特定时域平移相位差校正法(平移点数等于FFT变换点数)的校正精度的影响进行了理论分析和仿真计算。2007年杨志坚等[7-8]对高斯白噪声背景下时移相位差校正法和改变窗长相位差法校正精度进行了理论分析和仿真计算。目前还没有学者对噪声背景下能量重心法频率校正精度进行过系统分析。能量重心法是一种通用的校正方法，能对多段平均功率谱直接进行校正，特别适用于频率有微小波动的实际工况，在实际工程中得到很广泛的应用[9-10]，因此分析噪声背景下能量重心法校正精度有重要的意义。本文首先对高斯白噪声影背景下加对称窗的离散频谱利用能量重心校正法进行校正的频率校正精度进行了理论分析，推导其理论校正误差公式，进行了仿真验证，并且对能量重心法做了改进，使其具有更广的适用性和更高的抗噪能力。1试验方法和试验方案(四号宋体)1.1试验方法渗金属试验在自制双层辉光离子渗金属炉中进行，源极材料为HastelloyC—2000合金，尺寸为130mm×50mm×4mm，工件材料为20钢，尺寸为80mm×25mm×3mm。采用脉冲放电模式：源极采用直流电源，工件采用脉冲电源。源极材料HastelloyC—2000的质量分数：wNi=59%，wMo=16%，wCr=23%，wCu=1.6%，wC＜0.01%。1.2试验方案(五号黑体)为了选定正交试验各个工艺参数的取值范围，先结合以往试验研究的经验，然后又进行了20余炉的摸索性试验，确定了正交工艺参数。正交试验按照L16(45)正交表进行试验。指标项目为渗层表面合金元素成分及总质量分数、渗层厚度和吸收率。因素水平表如表1所示。表1因素水平表(小五黑体)因素(六号宋体)水平1234源极电压U／V10501000950900工件电压U／V275250350300气压p／Pa35304540极间距d／mm15202522.51.3频率校正误差理论分析没有噪声影响下的能量重心校正法的归一化频率校正公式[2]****1()nkiinnkkiinkiGfnG(1)式中*k——无噪声功率谱中最大值谱线*kiG——第*ki条功率谱线在白噪声背景下的归一化频率校正公式1()ˆrkrrnrkiinnkiinkiPfnP(2)式中rk——白噪声背景下功率谱中最大值谱线rkiP——第rki条功率谱线2数学模型在网络学习部分，采用三层BP神经网络来完成函数的映射。误差逆传播神经网络是一种具有三层或三层以上的阶层型神经网络。图1为一个三层前馈神经网络：它包括输入层、隐含层(中间层)、输出层；输入层有i个节点，隐含层有j个节点，输出层有t个节点。上、下层之间各神经元实现全连接，即下层的每一单元与上层的每一单元都实现权连接，而每层各神经元之间无连接。网络按有教师示教的方式进行学习，当一对学习模式提供给网络后，神经元激活值从输入层经各中间层向输出层传播，在输出层的各神经元获得网络的输入响应。这以后，按减小希望输出与实际输出之间误差的方向，从输出层经各中间层逐层修正各连接权值，最后回到输入层。算法步骤如下所述。(1)设置初始权系w(0)为较小的随机非零值。68(2)给定输入/输出样本对，计算网络的输出：设第p组样本输入、输出分别为up=(u1p,u2p,…,unp)dp=(d1p,d2p,…,dnp)p=1,2,…,L节点i在第p组样本输入时，输出为jjpijipipΙtwftxfy)()]([(3)式中，IjP为在第p组样本输入时，节点i的第j个输入。图1阻尼器屈服力矩—电流特性（图中文字为六号）f是激励函数，采用Sigmoid型，即)exp(11)(xxf(4)可由输入层经隐层至输出层，求得网络输出层节点的输入。(3)计算网络的目标函数J。设Ep为在第p组样本输入时网络的目标函数，取L2范数，则)(21)(21)(21)(2223tetydtydtEkkpkkpkpppp(5)式中，ykp(t)为在第p组样本输入时，经t次权值调整网络的输出，k是输出层第k个节点。网络的总目标函数为)()(tEtJpP(6)作为对网络学习状况的评价。判别：若J，则算法结束，否则，至步骤(4)。(4)反向传播计算。由输出层，依据J按“梯度下降法”反向计算，逐层调整权值(图3～5)。)()()()()1(twtwtJηtwtwijijijij)()()()(twtwtwtEηijijpijp(8)式中，为步长或称为学习率，n取1000000，取0.9。图2路径规划三维图(小五宋体)图3轮齿本体温度分布图4微压印法制备的微复制模板及仿生鲨鱼皮图5微生物细胞金属化约束成形工艺流程3计算结果与比较为检验程序的可靠性与实用性，同时又具有一定的推广能力，在正交试验16组数据中选择13组作为训练样本，余下的3组以及正交优化工艺作为69检测样本。经过1000000次训练的试验数据与计算结果的比较见表2、表3。表2预测渗层表面的成分编号质量分数w/%NiCrMoCu549.88114.69511.3651.605预测值51.73416.96211.2111.6611250.51416.69114.0710.912预测值51.17117.50913.0480.4151650.29019.1646.4801.308预测值50.82618.1047.81301.3331757.16819.65514.7321.353预测值59.89120.10712.1631.403对双层辉光离子多元共渗工艺参数：源极电压、工件电压、极间距、气压对渗层表面的合金元素总质量分数、渗层厚度、各合金元素质量分数、吸收率(工件增重/源网