1.一种基于自适应摩擦补偿的机械臂内模控制方法，其特征在于，包括以下步骤：步骤1：建立n自由度机械臂系统的动力学模型；步骤2：设计内模，将机械臂系统的位置跟踪和干扰抑制问题转化为由该系统和内模组成的增广系统的镇定问题；步骤3：设计最终控制器，对摩擦力进行补偿并解决增广系统的镇定问题；步骤1中，建立n自由度机械臂系统的动力学模型，其过程如下：1.1，根据欧拉—拉格朗日方法建立n自由度机械臂系统的动力学模型如下：其中 分别为机械臂的关节位置和速度向量， 为正定对称的惯性矩阵， 为向心力和科里奥利力矩， 为重力矩， 为输入力矩， 为摩擦力矩， 为外界扰动；上述机械臂系统具有以下性质：性质1：矩阵 和 满足以下条件：性质2：机械臂模型可进行下列参数线性化变换：其中 均为已知矩阵，满足 为机械臂中不确定参数；1.2，假设机械臂关节位置期望值为q d ，且q d , 均为连续有界信号，外部力矩干扰τ d (t)＝[τ d1 (t),…,τ dn (t)] T 为有限个阶跃信号和正弦信号的组合，形式如下：其中a i0 为常数，k i 为正弦信号的个数，a ij ,ω ij ,φ ij 分别为正弦信号的幅值，频率和相位；1.3，定义跟踪误差e＝q-q d 和 其中α为一个正常数，进一步可得 取 则有：1.4，系统(1)可以写为如下形式：结合性质2，(6)式可以进一步写为：1.5，引入Stribeck摩擦力模型对机械臂关节摩擦力矩进行描述，第i关节的摩擦力矩形式如下：其中对于i＝1,…,n，f ci 为库伦摩擦系数，f si 为静摩擦系数，f vi 为黏性摩擦系数， 为Stribeck速度，δ为经验参数；定义则机械臂摩擦力可表示为 其中 为已知矩阵，β为未知参数向量。
2.根据权利要求1所述的一种基于自适应摩擦补偿的机械臂内模控制方法，其特征在于，步骤2中，设计内模，将机械臂系统的位置跟踪和干扰抑制问题转化为由该系统和内模组成的增广系统的镇定问题，其过程如下：2.1，对τ di (t)，i＝1,…,n进行微分，可得：其中 均为实数，s i 为正常数；2.2，令 则有：其中 2.3，选择一对可控矩阵(M i ,N i )，其中 为Hurwitz矩阵， 则存在一个非奇异矩阵T i 满足以下Sylvester方程：2.4，定义如下记号：令ρ＝-Tυ,Ψ＝ΓT -1 ，根据公式(11)，可得：τ d ＝-Ψρ (15)2.5，设计内模为如下形式：2.6，进行如下坐标变换：根据性质2，得到如下误差方程：其中 2.7，设计如下动态补偿器：2.8，令 可得：则(7)式可以写为如下形式：其中 此时，机械臂系统的位置跟踪和干扰抑制问题被转化为由该系统和内模组成的增广系统的镇定问题。
3.根据权利要求2所述的一种基于自适应摩擦补偿的机械臂内模控制方法，其特征在于，步骤3中，设计最终控制器，对摩擦力进行补偿并解决增广系统的镇定问题；3.1，设计系统的状态反馈控制律为如下形式：其中K为正定对角常数矩阵， 为θ的估计值， 为β的估计值；3.2，取李雅普诺夫方程：其中Δ 1 为一个正定矩阵满足M T Δ 1 +Δ 1 M＝-I，I为单位矩阵，κ,λ为正常数， Λ 1 ,Λ 2 为任意的正定对角常数矩阵，则有：3.3，由于θ,β为未知常量， 选择自适应律为：此时，公式(24)可以化简为如下形式：设计Ψ为有界对角矩阵，存在正常数a 1 ,b 1 满足 选择合适的κ和K，使得κ>b 1 ，(K-a 1 I n )为正定矩阵，则式(26)可写为：其中a 0 为矩阵(K-a 1 I n )的最小特征值，b 0 ＝κ-b 1 ；经过稳定性分析可知，在所提出控制器的作用下，当t→∞时，闭环系统状态 s, 有界；对公式(27)两边对时间积分后可知： 根据Barbalat引理可知 则位置跟踪误差e和速度跟踪误差 最终可以收敛到零；3.4，则基于李雅普诺夫方法，可设计最终控制器如下：
