发明名称 一种基于多阶相关分析的弧面凸轮运动曲线辨识方法
摘要 本发明公开了一种基于多阶相关分析的弧面凸轮运动曲线辨识方法:(1)弧面凸轮运动曲线辨识;(2)对获取的弧面凸轮运动采集数据进行相关分析并拟合出实际运动曲线,进而确定弧面凸轮实际运动曲线与理想的运动规律的相似性;(3)多阶相关分析,用于对获取的弧面凸轮运动采集数据进行多阶相关分析,进行弧面凸轮运动规律的一阶导数、二阶导数和多阶导数与理想的运动规律的相关分析,得到弧面凸轮实际运动曲线在弧面凸轮的位置、速度和加速度的相关系数及其相似形态类别,最终辨识出弧面凸轮实际运动曲线。本发明不仅有利于提高机械行业中机械运动机构的位置精度,而且有利于提高机械运动机构的速度和加速度精度,进而提高运动平稳性。
申请公布号 CN103268378B 申请公布日期 2016.03.30
申请号 CN201310182841.4 申请日期 2013.05.16
申请人 西安交通大学 发明人 要义勇;赵丽萍;王亚菲;赵虎;王龙轩
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 汪人和
主权项 一种基于多阶相关分析的弧面凸轮运动曲线辨识方法,其特征在于,该辨识方法包括以下步骤:(1)弧面凸轮运动曲线的最小二乘拟合:对采集到的弧面凸轮运动数据,采用最小二乘法进行曲线拟合,得到弧面凸轮的实际运动曲线;(2)相关分析:根据弧面凸轮的实际运动曲线与理想的弧面凸轮运动曲线,进行弧面凸轮运动规律的相关分析得到相关系数;弧面凸轮的实际运动曲线与理想的弧面凸轮运动曲线进行相关分析,包括如下具体流程:由实际运动曲线采集点和理想运动曲线点计算相关系数γ<sub>xy</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msub><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><msub><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msub><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub></mrow></mrow><mrow><msqrt><mrow><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msup><msub><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><mn>2</mn></msup></mrow><mo>-</mo><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>&CenterDot;</mo><msqrt><mrow><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msup><msub><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><mn>2</mn></msup></mrow><mo>-</mo><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000822117670000011.GIF" wi="1398" he="255" /></maths>其中,y<sub>ci</sub>表示实际运动曲线上对应x<sub>i</sub>的值,y<sub>ki</sub>表示理想的弧面凸轮运动曲线k上对应x<sub>i</sub>的值,n表示采集点个数;当K种理想的弧面凸轮运动曲线和弧面凸轮的实际运动曲线的相关系数γ<sub>xy</sub>≈1,说明弧面凸轮的实际运动曲线和K种理想的弧面凸轮运动曲线在位置曲线的相似程度非常高,需要进行一阶导数或二阶导数相关分析;由实际运动曲线采集点和理想运动曲线点计算相关系数,如果实际运动速度曲线和理想运动速度曲线的相关系数γ′<sub>xy</sub>还不能区分弧面凸轮的实际运动曲线的速度属性,则进行二阶求导和计算加速度的相关系数γ″<sub>xy</sub>,直到辨识出来弧面凸轮的实际运动规律:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>&gamma;</mi><mrow><mi>x</mi><mi>y</mi></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><mrow><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msubsup><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow><mo>&prime;</mo></msubsup><msubsup><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow><mo>&prime;</mo></msubsup></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msubsup><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow><mo>&prime;</mo></msubsup></mrow></mrow><mrow><msqrt><mrow><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msup><msubsup><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mn>2</mn></msup></mrow><mo>-</mo><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>&CenterDot;</mo><msqrt><mrow><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msup><msubsup><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mn>2</mn></msup></mrow><mo>-</mo><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000822117670000021.GIF" wi="1406" he="247" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>&gamma;</mi><mrow><mi>x</mi><mi>y</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>=</mo><mfrac><mrow><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msubsup><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><msubsup><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup></mrow><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msubsup><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup></mrow></mrow><mrow><msqrt><mrow><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msup><msubsup><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mn>2</mn></msup></mrow><mo>-</mo><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>y</mi><mrow><mi>c</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>&CenterDot;</mo><msqrt><mrow><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><msup><msubsup><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mn>2</mn></msup></mrow><mo>-</mo><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>y</mi><mrow><mi>k</mi><mi>i</mi></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000822117670000022.GIF" wi="1413" he="287" /></maths>其中,y′<sub>ki</sub>表示理想的弧面凸轮运动曲线k一阶导数上对应x<sub>i</sub>的值,y′<sub>ci</sub>表示实际运动曲线一阶导数上对应x<sub>i</sub>的值,y″<sub>ci</sub>表示实际运动曲线二阶导数上对应x<sub>i</sub>的值,y″<sub>ki</sub>表示理想的弧面凸轮运动曲线k二阶导数上对应x<sub>i</sub>的值,n表示采集点个数;(3)根据相关系数的差异完成弧面凸轮的实际运动曲线的识别过程。
地址 710049 陕西省西安市咸宁西路28号