线性代数在动力机械工程领域中的应用院系：能源动力学院 学号：201364180 姓名：董书金 在机械工程领域复杂线性方程组的数值求解是经常遇见的问题，而且机械工程中的一些 多解问题，例如机构转配构型，机器人机构树状解和设计方案的多解问题等，常常需要线性 代数中线性方程的一些理论求解。并且线性代数中的公式通用于能淬火硬化的各种碳素钢及 合金钢。实际上，这些方程可以当作是一种定量尺度，广泛用于设计或选择钢种、制定或修订 标准、控制熔炼成分等方面。现代飞行器外形设计，这个就需要先研究飞机表面的气流的过 程。把飞行器的外形分成若干大的部件，每个部件沿着其表面又用三维的细网格划分出许多 立方体，这些立方体包括了机身表面以及此表面内外的空气。对每个立方体列写出空气动力 学方程，其中包括了与它相邻的立方体的共同边界变量，这些方程通常都已经简化为线性方 程。结合高等代数广泛用于研制和提供能量转换机械，包括将热能、化学能、原子能、电能、 流体压力能和天然机械能转换为适合于应用的机械能的各种动力机械以及将机械能转换为 所需要的其他能量的能量变换机械。一、线性方程数据处理在理工科学习中的基础运用用于计算多元或者单元复杂结构极限（线性代数原理运用软件：MATLAB），解放人 工计算无法解决计算的复杂问题。1.1导数七 极值和积分弋Taylor&jG命令功能limitXtx 耳)计算 锂尹町计算悭Z3limit(flx,a right)计算竽侧極限血gleft)tl算单侧桃限lim8x注意：存片右楓.限不相等其计算原理用于求解导数，多元函数的偏导数，定积分，多重积分，进而解决实际 运用中计算不规则曲面，形体的面积，体积等，并用于航空器外壳，船舶形体量， 汽车制造，精密机械制造等工程设计中。1.描述n阶线性时不变（LTI ）连续系统的微分方程为dnydn_iydy, dmu , duaF aF + aa y = b+ bb u,n三m dtn2dtndtn+i1dtmm dtm+i我片右極.限有亍-不存=存时-M也TLAB的默认状心是求右极限*其中Pi，PJ，P是特征方程a,九n+a九n-i+.+ a九+ a,=0的根，它们可用roots(a)12n12nn+1语句求得。各系数C，cn由y及其各阶导数的初始值来确定。对此有C1+ 叨坨=yy = y(0)pn1C + pn1C HF pn1C1122n n卩心+ p2C2+ PnCn=Dy0(Dy0表示y的导数的初始值y(0)_ 11 1 -C 11y0 _写成矩阵形式为p1P 2pn C2=Dy丿0pn-1 *1pn-1p n-1nC1- nD n-1 y 0=D n-1 y0即 Vc = Y0，其解为C =V Y0式中C = C , C，C T； Y = y ,Dy ,，Dn-iy T1 2n0000p1Pn-1L1Pn-12pn-1nV为范德蒙矩阵，在MATLAB的特殊矩阵库中有vander函数可直接生成。 (参考文献：郭龙先，张毅敏，何建琼.高等代数M.北京：科学出版社，2011)用于求解多元线性，非线性方程或多解，方程的根，求和与级数求和，极值点，概 率论问题模拟分析等问题，进而运用于天气预测，流体力学分析，经济学最优解， 机械最优结构，空气动力学等多领域中。1在热传导的研究中，一个重要的问题是确定一块平板的稳态温度分布.根据定律， 只要测定一块矩形平板四周的温度就可以确定平板上各点的温度.图8 一块平板的温度分布图【模型准备】如图9所示的平板代表一条金属梁的截面.已知四周8个节点处的温度(单 位C),求中间4个点处的温度T1, T2, T3, T4.90T1T26080T3T45060.501图9一块平板的温度分布图 【模型假设】假设忽略垂直于该截面方向上的热传导, 并且每个节点的温度等于与它相邻的 四个节点温度的平均值.【模型建立】根据已知条件和上述假设,有如下线性方程组f 1T = (90 +100 + T + T)1 4 2 3T = -(80 + 60 + T + T )2 414T = -(80 + 60 + T + T )3 414T = -(50 + 50 + T + T)、4 423【模型求解】将上述线性方程组整理得4T - T - T = 190-T + 4T 3 - T = 140-T44T - T = 140 -1 - T - T3 + 4犷=100234在Mat lab命令窗口输入以下命令 A = 4,-1,-1,0;-1,4,0,-1;-1,0,4,-1;0,-1,-1,4; b = 190;140;140;100; x = Ab; xMatlab执行后得ans =82.9167 70.8333 70.8333 60.4167可见片=82.9167, T2 = 70.8333, T3 = 70.8333, T4 = 60.4167.(参考文献：陈怀琛, 高淑萍, 杨威, 工程线性代数 , 北京: 电子工业出版社 , 2007. 页码:15-16.)其计算原理用于平面解析几何分析，线性模拟计算，非线性多方向定向求解等复杂数学数据处理过程，进而用于平面设计，立体结构分析及结构优化，物理复杂受力求解，航天模拟，电路设计优化，发动机设计，能量转换优化等众多能源动力学领域中。1.设平移变换为(x, y) t (x+a, y+b)旋转变换(绕原点逆时针旋转e角度)为(x, y) t (xcos0 - ysinG, xsin0 + ycosG)放缩变换(沿x轴方向放大s倍，沿y轴方向放大t倍)为 (x, y) t (sx, ty)【模型求解】R2中的每个点(x, y)可以对应于R3中的(x, y, 1).它在xOy平面上方1单位 的平面上.我们称(x, y, 1)是(x, y)的齐次坐标.在齐次坐标下，平移变换(x, y) t (x+a, y+b)10 ar xr x + a 于是可以用矩阵乘积0 1 by=y + b(0 0 1 丿U丿L 1丿可以用齐次坐标写成(x, y, 1) t (x+a, y+b, 1).实现.旋转变换(x, y) t (xcos0 - ysin0, xsin0 + ycosG) 可以用齐次坐标写成(x, y, 1) t (xcos0 - ysin0, xsin0 + ycosG, 1).r cos 0- sin 00、r xrxcos0 - y sin0、于是可以用矩阵乘积sin 0 cos 00y=x sin 0 + y cos0实现L 0 0 1,L1,L 1 ,放缩变换(x, y) t (sx, ty) 可以用齐次坐标写成(x, y, 1) t (sx, ty, 1).xsxy=tyL1,实现.s 0于是可以用矩阵乘积0 t电路是电子元件的神经系统.0参数的计算是电路设计的重要环节.其依据来自两个方面：一是客观需要，二是物理学定律.图22 USB扩展板【模型准备】假设图23中的方框代表某类具有输入和输出终端的电路.用(v )1I i丿记录输入电压和输入电流(电压v以伏特为单位，电流i以安培为单位)，用(v )2I i丿v 2 y记录输出电压和输入电流.若(v )2I i丿v 2 y(v ).1I i丿称矩阵A为转移矩阵.输入终端v1电路输出终端v2图23 具有输入和输出终端的电子电路图图24给出了一个梯形网络，左边的电路称为串联电路，电阻为匕(单位：欧姆).右边的电路是并联电 路，电路R2.利用欧姆定理和楚列斯基定律，我们可以得到串联电路和并联电路的转移矩阵分别是一 1! AAA22:L3.v1R1V2 :V31i1i1 1 1 1(1R)(10 、011J和1/ R1 2就是图22中梯形网络的转移矩阵.于是，原问题转化为求R, R?的值使得1_1/ R 1 + R / R211_0.5【模型求解】由l_1/ R 1 + R / R211_0.5根据其中的前两个方程可得R1 = & 这就是说在图22中梯形网络中取RR=2 8=， 2R.22丿_R1_1/ R =_0.521 + R / R = 512= 8， R2 = 2代入上面的第三个方程确实能使等式成立.把R1= 2即为所求.串联电路并联电路图24 梯形网络(1 -8)设计一个梯形网络，其转移矩阵是l_o.55丿模型假设】假设导线的电阻为零.【模型建立】设人和A2分别是串联电路和并联电路的转移矩阵，则输入向量x先变换成A 再变换到 A2(A1x). 其中10、1 - R r 1-R 、一1/ R 1 丿、0 1=-1/ R1 + R 1R ,212A2A1 =2【模型分析】若要求的转移矩阵改为(_0.5_8 )4 丿， 则上面的梯形网络无法实现.因为这时对应的方程组分丿就土拯象钳堵進度殆方碣=& + QIR = -8i是|1/ R2 =5根据前两个方程依然得到R1 = & R2 = 2,但把R1 = & R2 = 2代入上第三个方程1 + R / R 二 4V 1 2 却不能使等式成立.（参考文献David C. Lay,线性代数及其应用，沈复兴，傅莺莺等译，北京：人民邮电出版社，2009. 页码：129-130.）二、线性代数学理论在能源动力工程中的综合运用航空飞行分析与设计叫二卜;J一“ 匚n二-：卜： 、* - a 丿 1 a J 匕 vj鱼hC45住巴OS# 竝*戒11曲 KbiuA + wf co$i7 j一 3 Keflsdtain/? + w 严 in 肝inQ 十(U 二G Jcoti H札血出+讥、妙 Y 氐 EcO* + iifiy ncs&in y f + nOsyf+ (Yf+ 41054.)作业：推导 气流坐标系 下的角速度 关系式%=讥+讥血