矩阵论研究报告.doc

矩阵论在方程解耦及最小二乘法中的应用摘要：模态（也称为固有振动模态，或主模态）是多自由度线性系统的一种固有属性，可由系统的特征值（也称为固有值）与系统的特征矢量（也称为固有矢量，或者主振型）二者共同来表示的；它们分别从时空两个方面来刻画系统的振动特性。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型，其可以使得耦合方程组解耦。作用于一个n维自由度系统，可以转换到模态坐标下来解耦，确定在模态坐标下响应，然后通过线性变换得到物理坐标下的响应。惯常使用中，将线性定常系统振动微分方程组中的物理坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，以便求出系统的模态参数1。在科学实验和工程计算中，我们希望从给定的数据出发,构造一个近似函数，使数据点均在离曲线的上方或下方不远处，所求的曲线称为拟合曲线，它既能反映数据的总体分布，又不至于出现局部较大的波动，更能反映被逼近函数的特性，使求得的逼近函数与已知函数从总体上来说其偏差按某种方法度量达到最小，这就是最小二乘法。最小二乘法（又称最小平方法）是一种数学优化技术，它通过最小化误差的平方和寻找数据的最佳函数匹配，使这些求得的数据与实际数据之间误差的平方和为最小2，则需要范数的知识。关键字：模态，方程解耦，最小二乘一、引言数学中解耦是指使含有多个变量的数学方程变成能够用单个变量表示的方程组，即变量不再同时共同直接影响一个方程的结果，从而简化分析计算。通过适当的控制量的选取，坐标变换等手段将一个多变量系统化为多个独立的单变量系统的数学模型，即解除各个变量之间的耦合。对离散型函数(即数表形式的函数)考虑数据较多的情况.若将每个点都当作插值节点，则插值函数是一个次数很高的多项式，比较复杂，而且由于龙格振荡现象，这个高次的插值多项式可能并不接近原函数。最小二乘法在实际工程数据处理中应用广泛，在工程问题中，使用最小二乘法根据两个变量的几组实验数据可以轻松的找出这两个变量的函数关系的近似表达式,拟合成一条曲线来反映所给数据点总趋势5。二、预备知识2.1 坐标变换设是线性空间上的线性变换，（，）和（，）是的两组基，在这两个基下的表示矩阵分别为A，B则：（，）=（，）A ；（，）=（，）B原文内容:选择参考点或坐标的原则是为了描述方便，对于振动系统一般采用自然空间坐标。设基变换公式为（，）=（，）C ，C为变换矩阵则 B=2.2 范数如果是数域上的线性空间，且对于的任以向量，对应于一个实数函数，它满足如下三个条件。i. 非负性 当时；当时，；ii. 齐次性 ；iii. 三角不等式 ；则称为上的范数。可以证明对于向量的长度是一种范数，我们称为2-范数，记为。三、坐标变换和2-范数在工程实践中运用3.1 坐标变换在多自由度振动系统解耦中的运用3.1.1 多自由度系统运动方程描述多自由度系统一般运动方程 ,M为质量矩阵，K为刚度矩阵，M、K且都为正定矩阵。振动响应：代入运动方程：有非零解的充分必要条件： （3-1）得到特征多项式 得出每一阶固有圆频率 i=1,2.n其中：为特征值（固有频率）， 为特征向量（模态）。推出 （3-2）描述了系统做第i阶主振动时具有的振动形态，称为第i阶主振型，或第i阶模态。3.1.2 模态关于质量矩阵和刚度矩阵的正交性由（3-2）式，得：转置后右乘 =左乘 = 两式相减：（-）=0 恒成立（3-3）当时， 模态关于质量矩阵正交， 模态关于刚度矩阵正交当ij时 第i 阶模态主质量， 第i 阶模态主刚度，其中为第i 阶主模态。3.1.3 运动方程的解耦将组成矩阵, 称为模态矩阵。分别对M、K处理，得以下等式 （3-4） （3-5）即 为主质量矩阵，为主刚度矩阵坐标变换： 则运动方程变为 （3-6） （3-7）物理空间 耦合主模态空间 解耦在主坐标下成为n个独立的单自由度运动方程，可见实现了解耦，展开写即：， 第1阶固有频率 第2阶固有频率 第n阶固有频率 3.2 2-范数在最小二乘问题中的运用3.2.1 实际问题描述一颗导弹从敌国发射，通过雷达我们观测到了它的飞行轨迹，具体有如下数据：表1水平距离/m02505007501000高度/m08151920我国军情处分析得出该导弹沿抛物线轨道飞行。问题：预测该导弹在什么水平距离着地？3.2.2 直线拟合基本理论已知数据点，分布大致为一条直线。作拟合直线，该直线不是通过所有的数据点，而是使偏差平方和为最小，其中每组数据与拟合曲线的偏差为 根据最小二乘原理，应该取和是有极小值，故和应满足下列条件： （3-8）即得如下正规方程： （3-9）3.2.3 多项式拟合基本理论有时所给数据点的分布并不一定近似地呈一条直线，这时仍用直线拟合显然是不合适的，可用多项式拟合7。对于给定的一组数据寻求次数不超过n (nm ) 的多项式，来拟合所给定的数据，与线性拟合类似，使偏差的平方和为最小，由于Q可以看作是关于 ( j=0,1,2, n)的多元函数，故上述拟合多项式的构造问题可归结为多元函数的极值问题。令，得，即有 （3-10）这是关于系数的线性方程组，通常称为正规方程组，有惟一解。3.2.4 最小二乘问题的数值结果导弹沿抛物线飞行，设导弹前进方向为轴，垂直高度为轴，建立直角坐标系。表2123450250500750100008151920设所求导弹飞行抛物线方程为计算可得： 则得到方程组：写成矩阵形式为解之得：平方误差：则所求抛物线方程为：再由y=0，得x=2044.233m即导弹着地时的水平距离约为2044m。应用MATLAB得到图1图1程序见附录。参考资料1 胡海岩. 机械振动与冲击M. 航空工业出版社, 2002.2 故海岩. 机械振动基础M. 北京航空航天大学出版社, 2005.3 季文美. 模态分析技术M. 科学出版社, 1985.4 郑兆昌. 机械振动（上册）M. 机械工业出版社, 1980.5 许君一. 方向控制最小二乘法理论M. 清华大学出版社, 2004.6 李德葆. 实验模态分析与应用M. 科学出版社, 2005.7