矩阵论研究报告

1矩阵论在方程解耦及最小二乘法中的应用摘要：模态（也称为固有振动模态，或主模态）是多自由度线性系统的一种固有属性，可由系统的特征值（也称为固有值）与系统的特征矢量（也称为固有矢量，或者主振型）二者共同来表示的；它们分别从时空两个方面来刻画系统的振动特性。模态是机械结构的固有振动特性，每一个模态具有特定的固有频率、阻尼比和模态振型，其可以使得耦合方程组解耦。作用于一个n维自由度系统，可以转换到模态坐标下来解耦，确定在模态坐标下响应，然后通过线性变换得到物理坐标下的响应。惯常使用中，将线性定常系统振动微分方程组中的物理坐标变换为模态坐标，使方程组解耦，成为一组以模态坐标及模态参数描述的独立方程，以便求出系统的模态参数[1]。在科学实验和工程计算中，我们希望从给定的数据出发,构造一个近似函数，使数据点均在离曲线的上方或下方不远处，所求的曲线称为拟合曲线，它既能反映数据的总体分布，又不至于出现局部较大的波动，更能反映被逼近函数的特性，使求得的逼近函数与已知函数从总体上来说其偏差按某种方法度量达到最小，这就是最小二乘法。最小二乘法（又称最小平方法）是一种数学优化技术，它通过最小化误差的平方和寻找数据的最佳函数匹配，使这些求得的数据与实际数据之间误差的平方和为最小[2]，则需要范数的知识。关键字：模态，方程解耦，最小二乘一、引言数学中解耦是指使含有多个变量的数学方程变成能够用单个变量表示的方程组，即变量不再同时共同直接影响一个方程的结果，从而简化分析计算。通过适当的控制量的选取，坐标变换等手段将一个多变量系统化为多个独立的单变量系统的数学模型，即解除各个变量之间的耦合。对离散型函数(即数表形式的函数)考虑数据较多的情况.若将每个点都当作插值节点，则插值函数是一个次数很高的多项式，比较复杂，而且由于龙格振荡现象，这个高次的插值多项式可能并不接近原函数。最小二乘法在实际工程数据处理中应用广泛，在工程问题中，使用最小二乘法根据两个变量的几组实验数据可2以轻松的找出这两个变量的函数关系的近似表达式,拟合成一条曲线来反映所给数据点总趋势[5]。二、预备知识2.1坐标变换设是线性空间上的线性变换，（1，…，n）和（1，…，n）是的两组基，在这两个基下的表示矩阵分别为A，B则：（1，…，n）=（1，…，n）A；（1，…，n）=（1，…，n）B设基变换公式为（1，…，n）=（1，…，n）C，C为变换矩阵则B=1CAC2.2范数如果V是数域K上的线性空间，且对于V的任以向量，对应于一个实数函数，它满足如下三个条件。i.非负性当0时0；当0时，0；ii.齐次性,aaV；iii.三角不等式,,V；则称为V上的范数。可以证明对于向量12(,,,)n的长度22212n是一种范数，我们称为2-范数，记为2。三、坐标变换和2-范数在工程实践中运用3.1坐标变换在多自由度振动系统解耦中的运用33.1.1多自由度系统运动方程描述多自由度系统一般运动方程0KXXMnRX,M为质量矩阵，K为刚度矩阵，M、KnnR且都为正定矩阵。振动响应：)sin(tφX代入运动方程：0)(2φMKφ有非零解的充分必要条件：0MK2（3-1）得到特征多项式22(1)2110nnnnaaa得出每一阶固有圆频率ii=1,2…..n其中：为特征值（固有频率），φ为特征向量（模态）。推出2()()iiKM0φ（3-2）()iφ描述了系统做第i阶主振动时具有的振动形态，称为第i阶主振型，或第i阶模态。3.1.2模态关于质量矩阵和刚度矩阵的正交性由（3-2）式，得：转置后右乘()jφ()()iTjKφφ=2()()iTiiMφφ左乘()iTφ()()iTjKφφ=2()()iTijMφφ两式相减：（2i-2j）()()iTiMφφ=0恒成立（3-3）当ij时，ij()()0iTjMφφ模态关于质量矩阵正交，()()0iTjKφφ模态关于刚度矩阵正交当i＝j时()()iTipimMφφpim第i阶模态主质量，()()iTipikKφφpik第i阶模态主刚度，其中()iφ为第i阶主模态。()2()()2()iiijjjKMKMφφφφ43.1.3运动方程的解耦将()(1~)iinφ组成矩阵(1)()[]nΦφφnnR,Φ称为模态矩阵。分别对M、K处理，得以下等式(1)()(1)()[][]TnTnΦMΦMφφφφ(1)(1)()()[]TnnTMφφφφ(1)(1)()()[]TnnTMMφφφφ(1)(1)(1)()()(1)()()TTnnTnTnMMMMφφφφφφφφ100ppnmm（3-4）(1)()(1)()[]K[]TnTnKΦΦφφφφ(1)(1)()()[]TnnTKφφφφ(1)(1)()()[]TnnTKKφφφφ(1)(1)(1)()()(1)()()TTnnTnTnφKφφKφφKφφKφ100ppnkk（3-5）即11(,,)(,,)TppnpTppnpdiagmmdiagkkΦMΦMΦKΦKpM为主质量矩阵，pK为主刚度矩阵坐标变换：YXΦ则运动方程0KXXM变为0TTΦMΦYΦKΦY（3-6）0PPMYKY（3-7）物理空间0KXXM耦合主模态空间0piipiimyky解耦5在主坐标下成为n个独立的单自由度运动方程，可见实现了解耦，展开写即：11110ppmyky，第1阶固有频率111ppkm22220ppmyky第2阶固有频率222ppkm……第n阶固有频率pnnpnkm3.22-范数在最小二乘问题中的运用3.2.1实际问题描述一颗导弹从敌国发射，通过雷达我们观测到了它的飞行轨迹，具体有如下数据：表1水平距离/m02505007501000高度/m08151920我国军情处分析得出该导弹沿抛物线轨道飞行。问题：预测该导弹在什么水平距离着地？3.2.2直线拟合基本理论已知数据点miyxii,,2,1,,，分布大致为一条直线。作拟合直线xaaxy10)(，该直线不是通过所有的数据点iiyx,，而是使偏差平方和miiiyxaaaaF121010)(),(为最小，其中每组数据与拟合曲线的偏差为iiiiyxaayxy10)(mi,,2,1根据最小二乘原理，应该取0a和1a是),(10aaF有极小值，故和应满足下列条件：0pnnpnnmyky60)(2),(0)(2),(110110110010imiiimiiixyxaaaaaFyxaaaaaF（3-8）即得如下正规方程：miiimimiiimiimiiyxxaxayxama1110211110（3-9）3.2.3多项式拟合基本理论有时所给数据点的分布并不一定近似地呈一条直线，这时仍用直线拟合显然是不合适的，可用多项式拟合[7]。对于给定的一组数据,,1,2,,iixyim寻求次数不超过n(nm)的多项式，2012nnyaaxaxax来拟合所给定的数据，与线性拟合类似，使偏差的平方和210()mnjijiijQyax为最小，由于Q可以看作是关于(j=0,1,2,…,n)的多元函数，故上述拟合多项式的构造问题可归结为多元函数的极值问题。令0,0,1,2,,kQkna，得10()0,0,1,,mnjkijiiijyaxxkn，即有0121011201niniiniiniiinnnniiniiiamaxaxyaxaxaxxyaxaxaxxy（3-10）这是关于系数ja的线性方程组，通常称为正规方程组，有惟一解。3.2.4最小二乘问题的数值结果导弹沿抛物线飞行，设导弹前进方向为x轴，垂直高度为y轴，建立直角坐标系。7表2设所求导弹飞行抛物线方程为2012yaaxax计算可得：512500iix5211875000iix5311562500000iix541211.382812510iix5162iiy5143750iiixy52134937500iiixy则得到方程组：012012120125250018750006225001875000156250000043750187500015625000001.38281251034937500aaaaaaaaa写成矩阵形式为01122525001875000621562500000437502500187500034937500187500015625000001.382812510aaa解之得：0120.22857,0.039829,0.000019429aaa平方误差：2*220[y]0.457142882miiy则所求抛物线方程为：20.228570.0398290.000019429yxxi12345ix02505007501000iy081519208再由y=0，得x=2044.233m即导弹着地时的水平距离约为2044m。应用MATLAB得到图1图1程序见附录。参考资料[1]胡海岩.机械振动与冲击[M].航空工业出版社,2002.[2]故海岩.机械振动基础[M].北京航空航天大学出版社,2005.[3]季文美.模态分析技术[M].科学出版社,1985.9[4]郑兆昌.机械振动（上册）[M].机械工业出版社,1980.[5]许君一.方向控制最小二乘法理论[M].清华大学出版社,2004.[6]李德葆.实验模态分析与应用[M].科学出版社,2005.[7]刘炳涛.最小二乘法[M].测绘出版社,2001.附录clcclearx=[02505007501000];y=[08151920];p=polyfit(x,y,2);holdon;plot(x,y,'o')plot(x,polyval(p,x),'r');