线性代数实践matlab4.ppt

7.6.4 弹性梁的柔度矩阵 ? 设简支梁如图 7.3 所示，在梁的三个位置分别施加 力 f1 ， f2 和 f3 后，在该处产生的综合变形为图示的 y1 ， y2 和 y3 ，通常称为挠度。根据虎克定律，在 材料未失去弹性的范围内，力与它引起的变形呈 线性关系，可以写出： ? 矩阵中的元素 d 为单位力 f 引起的挠度，它愈大， 表明这个梁愈柔软。 1 11 12 13 1 2 21 22 23 2 3 31 32 33 3 y d d d f y d d d f y Df y d d d f ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 数字实例 ? 设柔度矩阵 （ 1 ）在 1 ， 2 ， 3 处施加的力为 30 ， 50 和 20 试求出其挠度。 （ 2 ）要在 3 处产生 0.4 挠度，其他两处为零，求应加的力。 ? 程序 ag764 ? D ? 0.001*[5,2,1;2,4,3;1,3,6] % 输入柔度矩阵 ? f ? [30;50;20], y ? D*f （ 排齐） % 给定力，求挠度 ? y1 ? [0;0;0.4] % 给定挠度， ? K ? inv(D), f1 ? K*y1 % 求刚度矩阵 , 求力 .005 .002 .001 .002 .004 .003 .001 .003 .006 D ? ? ? ? ? ? ? ? ? ? ? 梁的刚度矩阵计算 ? 柔度矩阵的逆就是刚度矩阵 K ， K ? D ? 1 ， ? 其中 1 11 12 13 1 2 21 22 23 2 3 31 32 33 3 f k k k y f k k k y f k k k y ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 254.2373 152.5424 33.8983 ( ) 152.5424 491.5254 220.3390 33.8983 220.3390 271.1864 K inv D ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 7.6.5 网络和图 ? 图为 1 ， 2 ， 3 ， 4 四个城市之间的空运航线，用有向图表 示。则该图可以用下列航路矩阵表示： ? 经过一次转机（也就是坐两次航班）能到达的城市，可 以由邻接矩阵的平方 A2 ? A1^2 来求得。 0 0 1 1 1 0 0 0 1 0 1 0 0 1 0 1 0 A ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 0 0 1 1 0 0 1 1 1 1 1 0 1 0 0 0 1 0 0 0 0 0 1 1 2 1 1 0 1 0 0 0 1 0 0 1 0 0 0 1 0 1 0 1 0 1 0 0 1 1 1 A A A ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 第 8 章 用向量空间解方程组 ? 8.1 向量和向量空间 二维空间 R 2 中的向量用两个沿列向的元素表示 – u=[2;4]; v=[3;-1]; – plot([2,3],[4, ? 1],x);hold on – % 若用中的子程序 drawvec ， – drawvec(u);hold on – drawvec(v,g);hold off 1 1 2 2 2 3 , , 4 1 u v u v u v ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 二维向量张成的空间 ? 平面上的任何一点 [w1;w2] 是不是一定能用 u 和 v 的线性组 合来实现？即是不是一定能找到一组常数 [c1,c2] ，使得 ? c1,c2 取所有可能的值，得到的 w 的集合就是 u 和 v 张成的子 空间，在所给的 u 和 v 下，它是一个平面。 ? 若 u 和 v 两个向量的各元素成简单的比例关系，合成的向量 只能在一根直线上，不可能张