机械设计中病态问题的求解方法(在机构综合中的运用)

　　【摘要】
　　本文应用矩阵的“伪秩”概念，求解病态位置的机械设计问题。该方法计算稳定性、收敛性好。
　　
　　【关键词】矩阵 伪秩 求解 病态
　　
　　引言在机构综合、机器人运动学计算中，常常出现线性方程组的求解问题，由于机构参数、机器人逆运动学计算不同位置以及非线性方程迭代的初始值选择等会使要求解的方程组本身出现病态(奇异)或近病态(奇异)方程组。同时由于问题中的数据本身具有一定的误差，加上计算机进行计算的过程中， 由于计算机的字长有限，不可避免地要产生舍入误差，这两种误差都会使现行方程组出现病态(奇异)和近奇异问题。如对奇异位置和近病态位置问题求解不当，所得到的计算结果并不是原始问题的解。本论文应用矩阵“伪秩”概念求解病态(奇异)方程，该方法计算稳定性、收敛性好，计算速度快。
　　
　　一、计算问题的良态、病态
　　
　　线性方程组
　　式中： 是含n个未知量的n维向量。
　　b是一个m维向量。
　　A是mXn阶矩阵。
　　线性方程组<1>的求解问题有两种。一种是相容方程组的求解，它是指存在向量x∈R满足方程组<1>。另一种是不相容方程组的求解，方程组<1>，当m>n或m<n但矩阵A的秩小于m时，它是指不存在一般意义下的X∈R满足方程组<1>。
　　
　　这种方程又称为矛盾方程组。
　　
　　不论是相容方程组还是矛盾方程组，方程组<1>的解为；式中：A 是矩阵A的广义逆矩阵。当矩阵A是一个n×n阶矩阵，且矩阵A秩等于n时，A =A一，A是矩阵A的逆矩阵。
　　
　　在求解方程<1>时，由于实际问题的数据误差及计算机的舍入误差，使得在应用计算机求解方程组<l>时，实际上是求解一个摄动方程：
　　
　　(A+6 A)(i+6 )= +6 6 <3>
　　
　　式中6 A、6 和6 6方程式矩阵A，常数变量b和未知变量x的摄动量，根据Banach~l理的推论可知，从理论上讲只要I A I l 6 A I<1R当l 6 6 l、I 6 A l都趋于零时，就是说只要I 6 A I和l 6 6 I足够小，总可以使解的摄动量I 6 i I 4,N任意规定的程度，而逆矩阵(A+6 A) 也可以逼近到A到任意的精确程度。实际上它仅仅是一种纯理论的结果，实际并非如此。其原因一方面是摄动量是不可避免的。因实际问题中数据往往是由于实验测定或计算所得到的，数据误差是不可避免的，另外计算机的字长限制不可避免地有舍入误差存在；另一方面从计算机字长有限这一观点来看，摄动量的相对精度如l 6 A l／l A I或I 6 l／Ib I是有一定局限的，它们并不能任意小。
　　
　　由上分析可知，对于一个给定的计算问题所必须研究的一个重要问题是：问题中参数的微小变动，对问题的解会产生什么样的影响。这就是问题(理论)解对于参数变动的敏感性问题。如果问题中参数的微小摄动能所引起的相对摄动不大，称这种计算问题为良态问题。但是若问题中参数有微小的相对摄动时，解会引起“巨大”的相对摄动， 这种问题称之为病态问题或是奇异问题。病态问题是计算数学领域中普遍存在的问题。因此，在进行实际计算时，如何判别问题的良态、病态、病态程度以及如何解决病态问题的求解是十分重要的。
　　
　　二、病态程度的度量
　　
　　对广义逆矩阵的计算问题以及与其相关联的最小二乘问题的病态程度，可用相应的矩阵条件数k(A)来度量。
　　
　　k(A)=I A I I A l <5>
　　
　　k(A)越大， “病态”越严重，k(A)越小， “良态”越好。
　　
　　三、病态问题的处理方法
　　
　　在许多实际计算问题中， 所遇到的参数矩阵A常常是满秩的(理论证明是满秩的)，但如果按满秩进行计算，不仅得小到好的计算结果，甚至会使计算不能进行下去。为解决此问题常采用降秩方法，引用矩阵的“伪秩”概念解决病态问题。
　　
　　设矩阵A是一个IlXn实矩阵，，其中是矩阵A的列向量。
　　
　　定义设e为给定的一个正的小数，若称矩阵A的e秩为r，其中E是向量对向量系相对相关指标(1)从实际计算的角度来讲，病态方程的真秩往往是难以确定的。而实际在许多计算中，并不需要计算矩阵的真秩。重要的是，为使计算能得到较好的结果，应当设法确定它的某一个秩， 即矩阵的伪秩。
　　
　　为了控制矩阵A病态的危害程度，常根据实际精度适当选取一个正小数e来确定A的一个降秩近似， 在计算中采取这种措施，就可以避免舍入误差影响的恶化，从而提高算法的数值稳定性。
　　
　　四、矩阵伪秩的确定和线性方程的求解
　　
　　为计算矩阵A的伪秩，对矩阵A进行Lu分解。
　　
　　设A是一个Ill×n阶矩阵，并假定伪秩r=Rank(A)，NA分解成一个Ill×r的单位下三角阵L；~lr×n的上三角阵U， 即为了使计算不遇到不应有的中断，并且加强数值稳定性，应采取选主元的措施。本论文采取全面选主元法对矩阵A作满秩Lu分解。
　　
　　矩阵A的秩实际上事先并不知道，只能在计算过程中确定它的一个适当的秩。为此，根据实际计算要求，给定一个控制量￡(一般比计算机相对精度要大一些)，如果在第k+l步所选的主元的绝对值小于e时，就确定r=k，即为矩阵A的e一秩。
　　
　　五、例题及计算结果
　　
　　(一)机构刚体导引设计同理，在连杆平面上再另任选一点D ，固定铰链中心0。
　　
　　的坐标的计算同前(见图1)。
　　
　　由于C 、DI点是在连杆平面上任意选择的，所以给定连杆平面的三个位置设计四杆机构可有无穷多解。
　　
　　(二)机构函数再现设计图2(q)所示，连架杆l、3对机架0 0 的三组对应位置分别为I、II、III与I 、II 、III ，相对应位移是Q，和B 及n 。和B 两连架杆角位置的对应关系仅与各杆的相对长度有关，适当选取机架0 0 的长度。图中O 、0 为固定铰链中心， 再任选定一个连架杆的长度， 图2(b)所示， 在连架杆3上选取0BB 。连接0AB2、OAB3，将A OAB208刚化后绕0^反转一Q 2至A0^B2 0B ，同理将△0AB30 刚化后绕0 反转～a 3至△0 B。”0 ”。转化后原连架杆0 B 变为连杆，即将函数再现设计转化为刚体导引设计，转化后得连杆的三个位置分别为：O B 、0B B 、0 B。 由于固定铰链中一G'OA、0B已确定，现在的设计任务是确定铰链中心A和B的位置。为此通过B2$HB3坐标B2(X YB2)和B3(x Y 。)使用转动矩阵[RM]计算B2图2(C)，已知连杆三个位置0 B 、0 B。 、0 ”B。 确定铰链中心A、B的位置可用以上介绍的刚体导引设汁方法进行设计计算。由于设计时固定铰链中心O 、OB及连架杆3的长度0 B 是任意选定的，所以已知连架杆三组对应关系设计四杆机构可有无穷多解。
　　
　　六、结论
　　
　　本论文引用矩阵“伪秩”概念求解病态线性方程组，用全面选主元法进行矩阵Lu分解。该方法在计算过程中可自动确定方程系数矩阵的秩，可同时适用于良态和病态方程的计算，在机械设计计算中，稳定性好、收敛性好、计算速度快。

发表科技论文网专业提供机械论文代写服务，并提供大量论文资料，如有业务需求请咨询网站客服人员！