矩阵论研究报告..docxVIP

矩阵论在方程解耦及最小二乘法中的应用摘要：模态（也称为固有振动模态，或主模态）是多自由度线性系统的一种固有属性，可由系统的特征值（也称为固有值）与系统的特征矢量（也称为固有矢量，或者主振型）二者共同来表示的；它们分别从时空两个方面来刻画系统的振动特性。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型，其可以使得耦合方程组解耦。作用于一个n维自由度系统，可以转换到模态坐标下来解耦，确定在模态坐标下响应，然后通过线性变换得到物理坐标下的响应。惯常使用中，将线性定常系统振动微分方程组中的物理坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，以便求出系统的模态参数[1]。在科学实验和工程计算中，我们希望从给定的数据出发,构造一个近似函数，使数据点均在离曲线的上方或下方不远处，所求的曲线称为拟合曲线，它既能反映数据的总体分布，又不至于出现局部较大的波动，更能反映被逼近函数的特性，使求得的逼近函数与已知函数从总体上来说其偏差按某种方法度量达到最小，这就是最小二乘法。最小二乘法（又称最小平方法）是一种数学优化技术，它通过最小化误差的平方和寻找数据的最佳函数匹配，使这些求得的数据与实际数据之间误差的平方和为最小[2]，则需要范数的知识。关键字：模态，方程解耦，最小二乘一、引言数学中解耦是指使含有多个变量的数学方程变成能够用单个变量表示的方程组，即变量不再同时共同直接影响一个方程的结果，从而简化分析计算。通过适当的控制量的选取，坐标变换等手段将一个多变量系统化为多个独立的单变量系统的数学模型，即解除各个变量之间的耦合。对离散型函数(即数表形式的函数)考虑数据较多的情况.若将每个点都当作插值节点，则插值函数是一个次数很高的多项式，比较复杂，而且由于龙格振荡现象，这个高次的插值多项式可能并不接近原函数。最小二乘法在实际工程数据处理中应用广泛，在工程问题中，使用最小二乘法根据两个变量的几组实验数据可以轻松的找出这两个变量的函数关系的近似表达式,拟合成一条曲线来反映所给数据点总趋势[5]。二、预备知识2.1 坐标变换设是线性空间上的线性变换，（，…，）和（，…，）是的两组基，在这两个基下的表示矩阵分别为A，B则：（，…，）=（，…，）A ；（，…，）=（，…，）B设基变换公式为（，…，）=（，…，）C ，C为变换矩阵则 B=2.2 范数如果是数域上的线性空间，且对于的任以向量，对应于一个实数函数，它满足如下三个条件。非负性 当时；当时，；齐次性 ；三角不等式 ；则称为上的范数。可以证明对于向量的长度是一种范数，我们称为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）当时， 模态关于质量矩阵正交， 模态关于刚度矩阵正交当i＝j时 第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高度/国军情处分析得出该导弹沿抛物线轨道飞行。问题：预测该导弹在什么水平距离着地？3.2.2 直线拟合基本理论已知数据点，分布大致为一条直线。作拟合直线，该直线不是通过所有的数据点，而是使偏差平方和为最小，其中每组数据与拟合曲线的偏差为 根据最小二乘原理，应该取和是有极小值，故和应满足下列条件： （3-8）即得如下正规方程： （3-9）3.2.3 多项式拟合基本理论有时所给数据点的分布并不一定近似地呈一条直线，这时仍用直线拟合显然是不合适的，可用多项式拟合[7]。对于给定的一组数据寻求次数不超过n (nm ) 的多项式，来拟合所给定