线性代数在实际生活中地应用_第1页
线性代数在实际生活中地应用_第2页
线性代数在实际生活中地应用_第3页
已阅读5页,还剩8页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、线性代数在生活中的实际应用 制药工程学院 环境科学 苏雷 10204118 大学数学是自然科学的基本语言, 是应用模式探索现实世界物质运动机理的主 要手段。学习数学的意义不仅仅是学习一种专业的工具而已。 ;初等的数学 知识 学习线性代数数学建模 函数模型的建立及应用, 作为变化率的额倒数在几 何学、物理学、经济学中的应用,抛体运动的数学建模及其应用,最优化方法及 其在工程、经济、农业等领域中的应用,逻辑斯谛模型及其在人口预测、新产品 的推广与经济增长预测方面的应用, 网络流模型及其应用, 人口迁移模型及其应 用,常用概率模型及其应用,等等。线性代数中行列式 实质上是又一些竖直排列形成的数表按一

2、定的法则计算得 到的一个数。早在 1683 年与 1693 年,日本数学家关孝和与德国数学家莱布尼 茨就分别独立的提出了行列式的概念。 之后很长一段时间, 行列式主要应用与对 现行方程组的而研究。 大约一个半世纪后, 行列式逐步发展成为线性代数的一个 独立的理论分支。 1750 年瑞士数学家克莱姆也在他的论文中提出了利用行列式 求解线性方程组的著名法则克莱姆法则。 随后 1812 年,法国数学家柯西发 现了行列式在解析几何中的应用, 这一发现机器了人们对行列式的应用进行探索 的浓厚兴趣。如今,由于计算机和计算软件的发展, 在常见的高阶行列式计算中, 行列式的数值意义虽然不大, 但是行列式公式依

3、然可以给出构成行列式的数表的 重要信息。在线性代数的某些应用中,行列式的只是依然非常重要。 例如:有甲、乙、丙三种化肥,甲种化肥每千克含氮 70 克,磷 8 克,钾 2 克; 乙种、 化肥每千克含氮 64 克,磷 10 克,钾 0.6 克;丙种化肥每千克含氮 70 克,磷 5 克,钾 1.4 克若把此三种化肥混合,要求总重量 23 千克且含磷149克,钾30克,问三种化肥各需多少千克?解.设甲、乙、丙三种化肥 各需X1、X2、X3千克,依 解解题意得方程组Xi X2 X323,8X110 X2 5X3149,2 Xi °6X24X33°.此方程组的系数行列式D 27,又D1

4、81, D227, D38155由克莱姆法则,此方程组有唯一解:Xi 3;X2 5;X3 15.即甲乙丙三种化肥各需 3千克5千克15千克、矩阵实质上就是一张长方形的数表,无论是在日常生活中还是科学研究中,矩 阵是一种非常常见的数学现象。学校课表、成绩单、工厂里的生产进度表、车站 时刻表、价目表、故事中的证券价目表、科研领域中的数据分析表,它是表述或 处理大量的生活、生产与科研问题的有力的工具。矩阵的重要作用主要是它能把 头绪纷繁的十五按一定的规则清晰地展现出来,使我们不至于背一些表面看起来 杂乱无章的关系弄得晕头转向。塌还可以恰当的给出事物之间内在的联系, 并通 过矩阵的运算或变换来揭示事物

5、之间的内在联系。 它也是我们求解数学问题时候“数形结合”的途径。矩阵的运算是非常重要的内容。n 111nnn1n 11例:计算nnn11 n 1nnnnn 11nn1n 1解:nn111nnn2n 11111n 112n11n112nn1 11111 n 11nnn 11 1n 1nn(n 1)nn1nn(n 1)n2nnnn(n 1)n 11n n1 n 1n n1n1n在此例中,A A,所以A是幕等矩阵.矩阵的初等变化,矩阵的秩,初等矩阵,线性方程组的解。向量组的线性相关, 向量空间,向量组的秩,n维向量。这些都是线性代数的核心概念。线性代数在 应用上的重要性与计算机的计算性能成正比例增长

6、。而这一性能伴随着计算机软 硬件的不断创新提升,最终,计算机并行处理和大规模计算的迅猛发展将会吧计 算机科学与线性代数紧密的联系在一起并广泛应用于解决飞机制造,桥梁设计, 交通规划,石油勘探,经济管理等科学领域。线性模型比复杂的非线性模型更易 于用计算机进行计算。线性方程组应用广泛。主要有网络流模型,人口迁移模型, 基因问题,求血液的流率和血管分支点出的压强等等。 线性方程组的解法其中至 关重要的例如求解齐次线性方程组x-i 2x2X3X402x1 x2 2x3 2X40x1 x2 4x3 3x40ri22r5 一 34Xi 2x3即得与原方程组同解的方程组X2 2X3Q 0X4X45-34-

7、3Xi 2X3由此即得X22X35x4 ,%( x3,x4可任意取值).令x3C, x4g,把它写成通常的参数形式Xi2C2X22C253C2,43C2,x12X3Ci,X2X3CiC2X4C2,X40对系数矩阵A施行初等行变换:i22i a2rii22i3 ai2i22 03640ii43 ari0364r2( 3)0解:A14-302 2 02 10这些内容不仅在数方阵的特征值、特征向量理论及方阵的相似对角化的问题, 学本身的研究中具有重要的作用,在其他的许多科学领域中也有重要的应用。 例 如,在生物信息学中,人类基因的染色体图谱在进行 DNA序列对比是就用到了矩阵的相似,这个概念。线性代

8、数学习对数学建模十分必要。那么,为什么线性代数得到广泛运用,也就是说,为什么在实际的科学研究中解线性方程组是经常 的事,而并非解非线性方程组是经常的事呢?这是因为,大自然的许多现象恰好 是线性变化的。按照辩证唯物主义的观点,世间的一切事物都是在不断地运动着 的.所谓运动,从数学上描述,就是随时间而变化,因此,研究各个量随时间的变化率,即导数,与各个量的大小之间的关系,就是非常重要的以下为线性代数实际解决的应用问题:例1:基因间“距离”的表示在ABO血型的人们中,对各种群体的基因的频率进行了研究。如果我们把四 种等位基因Ai,A2,B,O区别开,有人报道了如下的相对频率,见表1.1表1.1基因的

9、相对频率爱斯基摩人f1i班图人f2i英国人f3i朝鲜人f4iA10.29140.10340.20900.2208A20.00000.08660.06960.0000B0.03160.12000.06120.2069O0.67700.69000.66020.5723合计1.0001.0001.0001.000问题 一个群体与另一群体的接近程度如何?换句话说, 就是要一个表示基因 的“距离”的合宜的量度。解 有人提出一种利用向量代数的方法。 首先,我们用单位向量来表示每一个群体。为此目的,我们取每一种频率的平方根,记Xki , fki 由于对这四种群4体的每一种有 fkii 141,所以我们得到x

10、2 1.这意味着下列四个向量的每个i 1都是单位向量.记x11x21x31x41x12x22x32x42,a2x13,a3x23,a4x33x43x14x24x34x44在四维空间中,这些向量的顶端都位于一个半径为1的球面上.现在用两个向量间的夹角来表示两个对应的群体间的“距离”似乎是合理的 如果我们ai和a2之间的夹角记为B,那么由于 | ai|=| a2|=1 ,再由内只公式,得cos a?0.53980.32160.00000.2943,a20.17780.34640.82280.8307故cosa1 a2 0.9187得23.2 °按同样的方式,我们可以得到表1.2.表1.2

11、基因间的“距离”爱斯基摩人班图人英国人朝鲜人爱斯基摩人0 °23.2 °16.4 °16.8 °班图人23.2 °0 °9.8 °20.4 °英国人16.4 °9.8 °0°19.6 °朝鲜人16.8 °20.4 °19.6 °0 °由表1.2可见,最小的基因“距离”是班图人和英国人之间的“距离”,而爱斯基摩人和班图人之间的基因“距离”最大.例2 :在医药领域也有着很重要的作用。例如:通过中成药药方配制问题,达至V理解向量组的线性相关性

12、、最大线性无关组向量的线性表示以及向量空间等线性代数的知识问题:某中药厂用9种中草药(A-I ),根据不同的比例配制成了 7种特效药,各用量成分见表1 (单位:克)1号成药2号成药3号成药4号成药5号成药6号成药7号成药A10214122038100B1201225356055C531105140D79255154735二012255336F255355355550G94172523925Hf651610103510821202620解:(1 )把每一种特效药看成一个九维列向量,分析7个列向量构成向量组的线性相关性。若向量组线性无关,则无法配制脱销的特效药;若向量组线性相关,并且能找到不含u3

13、,u6的一个最大线性无关组,则可以配制3号和6号药品。在Matlab窗口输入u1=10;12;5;7;0;25;9;6;8;u2=2;0;3;9;1;5;4;5;2;u3=14;12;11;25;2;35;17;16;12;u4=12;25;0;5;25;5;25;10;0;u5=20;35;5;15;5;35;2;10;0;u6=38;60;14;47;33;55;39;35;6;u7=100;55;0;35;6;50;25;10;20;U=u1,u2,u3,u4,u5,u6,u7U0,r=rref(U)计算结果为U0=r= 1 2 4 5 71 0 1 0 0 0 0从最简行阶梯型 U0

14、中可以看0 1 2 0 0 3 0出,R ( U) =5 ,向量组线性0 0 0 1 0 1 0相关,一个最大无关组为0 0 0 0 1 1 0u1, u2, u4, u5, u7,0 0 0 0 0 0 1u3=u1+2u2四个零行u6=3u2+u4+u5 故可以配制新药2)三种新药用 v1 ,v2,v3 表示,问题化为 v1,v2 ,v3 能否由 u1-u7 线性表示,若能表示,则可配制;否则,不能配制。令 U=u1,u2,u3,u4,u5,u6,u7,v1,v2,v3U0,r=rref(U)由 U0 的最后三列可以看出结果计 算 结 果 r=1,2,4,5,7,10 则 可 以 看 出

15、v1=u1+3u2+2u4v2=3u1+4u2+2u4+u7v3 不能被线性表示,所以无法配制例 3 :化学方程的配平:( x1) C3 H 8 (x2 )O2 (x3)CO2 (x4)H2O确定 x1,x2,x3,x4 ,使两边原子数相等称为配平,方程为30x1 8 x2 00210x3 0 x4 221写成矩阵方程308002100 - 221x1x2x3x4Ax 0例 4 :卫星上用三种可见光和四种红外光进行摄像,对每一个区域,可以获得七张遥感图象。 利用多通道的遥感图可以获取尽可能多的地面信息, 因为各种地貌、作物和气象特征可能对不同波段的光敏感。 而在实用上应该寻找每一个地方的主因素

16、, 成为一张实用的图象。 每一个象素上有七个数据, 形成一个多元的变量数组,在其中合成并求取主因素的问题,就与线性代数中要讨论的特征值问题有关。例5 :用逆阵进行保密编译码在英文中有一种对消息进行保密的措施,就是把英文字母用一个整数来表示。然后传送这组整数。这种方法是很容易根据数字出现的频率来破译,例如出现频率特别高的数字,很可能对应于字母 E。可以用乘以矩阵A的方法来进一步加密。假如 A是一个行列式等于土 1的整 数矩阵,则A 1的元素也必定是整数。而经过这样变换过的消息,同样两个 字母对应的数字不同,所以就较难破译。接收方只要将这个消息乘以 A 1就可以复原。例6 : Euler的四面体问

17、题问题 如何用四面体的六条棱长去表示它的体积?这个问题是由Euler (欧拉)提出的解 建立如图2.1所示坐标系,设A, B, C三点的坐标分别为(a1,b1,C1),( a2,b2,C2)和(a3,b3,c),并设四面体 O-ABC的六条棱长分别为l ,m, n, p,q,r.由立体几何知道,该四面体的体积 V等于以向量OA,OB,OC1组成右手系时,以它们为棱的平行六面体的体积V6的6 .而aibiCiV6 OA OBOCa2b2C2a3b3C3aibiCi6V a2 b2 C2.将上式平方,得a3 b3 C3ab|Ciab|Cia 2b?C2a:p C2a3b3C3*3bjC336V22

18、aib22Cia a2bib2CiC2aia3b1b3a a:bib2G C22a2b;2C2a2a3b2b3aia3bib3C2C3a2a3b2b3C2C3al b3GC3C2C3203于是得根据向量的数量积的坐标表示,有OA OA2ai.22biCi,OA OBaa2bib2g,OAOCaia3bib3CiC3,OBOBa| b2 22C2OBOCa2a3db3C2C3,OCOCafb2 gOAOAOAOBOA36V2OAOAOBOCOBOBOBOCOBOCOCOCOC(2.1 )由余弦定理,可行OA OBq cos2P2同理OAOCOBOCr2 l22将以上各式代入(2.1)式,236V

19、2P2q22 2 2p r m222 222 2pq nPr m2222PPr l222 2pr l2 r2n2(2.2)于是这就是Euler的四面体体积公式 例7 : 一块形状为四面体的花岗岩巨石,量得六条棱长分别为l=10m,m=15m,n=12m,p=14m,q=13m,r=11m.则P2q2n2 110.5,-p2 r2m2p2 r2 l246,95.222代入(2.1)式,得196110.54636V2110.5169951369829.75.4695121于是V2 38050.82639(195m3)2.即花岗岩巨石的体积约为195m 3。古埃及的金字塔形状为四面体,因而可通 过测量其六条棱长去计算金字塔的体积。线性代数(Lin ear Algebra )是数学的一个分支,它的研究对象是向量,向量 空间(或称线性空间),线性变换和有限维的线性方程组。向量

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论