数学建模案例解析分析线性代数建模案例解析(20例).doc

收藏

编号:2604185    类型:共享资源    大小:3.33MB    格式:DOC    上传时间:2020-04-23
10
金币
关 键 词:
数学 建模 案例 解析 分析 线性代数 20
资源描述:
^` 线性代数建模案例汇编 张小向 东南大学数学系 2012年6月 目 录 案例一. 交通网络流量分析问题 1 案例二. 配方问题 4 案例三. 投入产出问题 6 案例四. 平板的稳态温度分布问题 8 案例五. CT图像的代数重建问题 10 案例六. 平衡结构的梁受力计算 12 案例七. 化学方程式配平问题 14 案例八. 互付工资问题 16 案例九. 平衡价格问题 18 案例十. 电路设计问题 20 案例十一. 平面图形的几何变换 22 案例十二. 太空探测器轨道数据问题 24 案例十三. 应用矩阵编制Hill密码 25 案例十四. 显示器色彩制式转换问题 27 案例十五. 人员流动问题 29 案例十六. 金融公司支付基金的流动 31 案例十七. 选举问题 33 案例十八. 简单的种群增长问题 34 案例十九. 一阶常系数线性齐次微分方程组的求解 36 案例二十. 最值问题 38 附录 数学实验报告模板 39 这里收集了二十个容易理解的案例. 和各类数学建模竞赛的题目相比, 这些案例确实显得过于简单. 但如果学生能通过这些案例加深对线性代数基本概念、理论和方法的理解, 培养数学建模的意识, 那么我们初步的目的也就达到了. 案例一. 交通网络流量分析问题 城市道路网中每条道路、每个交叉路口的车流量调查,是分析、评价及改善城市交通状况的基础。根据实际车流量信息可以设计流量控制方案,必要时设置单行线,以免大量车辆长时间拥堵。 图1 某地交通实况 图2 某城市单行线示意图 【模型准备】 某城市单行线如下图所示, 其中的数字表示该路段每小时按箭头方向行驶的车流量(单位: 辆). 500 1 2 3 4 400 300 100 200 300 x1 x2 x3 x4 图3 某城市单行线车流量 (1) 建立确定每条道路流量的线性方程组. (2) 为了唯一确定未知流量, 还需要增添哪几条道路的流量统计? (3) 当x4 = 350时, 确定x1, x2, x3的值. (4) 若x4 = 200, 则单行线应该如何改动才合理? 【模型假设】 (1) 每条道路都是单行线. (2) 每个交叉路口进入和离开的车辆数目相等. 【模型建立】 根据图3和上述假设, 在①, ②, ③, ④四个路口进出车辆数目分别满足 500 = x1 + x2 ① 400 + x1 = x4 + 300 ② x2 + x3 = 100 + 200 ③ x4 = x3 + 300 ④ 【模型求解】根据上述等式可得如下线性方程组 其增广矩阵 (A, b) = 由此可得 即 . 为了唯一确定未知流量, 只要增添x4统计的值即可. 当x4 = 350时, 确定x1 = 250, x2 = 250, x3 = 50. 若x4 = 200, 则x1 = 100, x2 = 400, x3 = -100 < 0. 这表明单行线“③④”应该改为“③④”才合理. 【模型分析】(1) 由(A, b)的行最简形可见, 上述方程组中的最后一个方程是多余的. 这意味着最后一个方程中的数据“300”可以不用统计. (2) 由可得, , , 这就是说x1, x2, x3, x4这四个未知量中, 任意一个未知量的值统计出来之后都可以确定出其他三个未知量的值. 参考文献 陈怀琛, 高淑萍, 杨威, 工程线性代数, 北京: 电子工业出版社, 2007. 页码: 16-17. Matlab实验题 某城市有下图所示的交通图, 每条道路都是单行线, 需要调查每条道路每小时的车流量. 图中的数字表示该条路段的车流数. 如果每个交叉路口进入和离开的车数相等, 整个图中进入和离开的车数相等. 300 500 150 180 350 160 220 300 100 290 400 150 x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 图4 某城市单行线车流量 (1)建立确定每条道路流量的线性方程组. (2)分析哪些流量数据是多余的. (3)为了唯一确定未知流量, 需要增添哪几条道路的流量统计. 案例二. 配方问题 在化工、医药、日常膳食等方面都经常涉及到配方问题. 在不考虑各种成分之间可能发生某些化学反应时, 配方问题可以用向量和线性方程组来建模. 图5 日常膳食搭配 图6 几种常见的作料 【模型准备】一种佐料由四种原料A、B、C、D混合而成. 这种佐料现有两种规格, 这两种规格的佐料中, 四种原料的比例分别为2:3:1:1和1:2:1:2. 现在需要四种原料的比例为4:7:3:5的第三种规格的佐料. 问: 第三种规格的佐料能否由前两种规格的佐料按一定比例配制而成? 【模型假设】 (1) 假设四种原料混合在一起时不发生化学变化. (2) 假设四种原料的比例是按重量计算的. (3) 假设前两种规格的佐料分装成袋, 比如说第一种规格的佐料每袋净重7克(其中A、B、C、D四种原料分别为2克, 3克, 1克, 1克), 第二种规格的佐料每袋净重6克(其中A、B、C、D四种原料分别为1克, 2克, 1克, 2克). 【模型建立】 根据已知数据和上述假设, 可以进一步假设将x袋第一种规格的佐料与y袋第二种规格的佐料混合在一起, 得到的混合物中A、B、C、D四种原料分别为4克, 7克, 3克, 5克, 则有以下线性方程组 【模型求解】上述线性方程组的增广矩阵 (A, b) =, 可见 又因为第一种规格的佐料每袋净重7克, 第二种规格的佐料每袋净重6克, 所以第三种规格的佐料能由前两种规格的佐料按7:12的比例配制而成. 【模型分析】(1) 若令a1 = (2, 3, 1, 1)T, a2 = (1, 2, 1, 1)T, b = (4, 7, 5, 3)T, 则原问题等价于“线性方程组Ax = b是否有解”, 也等价于“b能否由a1, a2线性表示”. (2) 若四种原料的比例是按体积计算的, 则还要考虑混合前后体积的关系(未必是简单的叠加), 因而最好还是先根据具体情况将体积比转换为重量比, 然后再按上述方法处理. (3) 上面的模型假设中的第三个假设只是起到简化运算的作用. 如果直接设x克第一种规格的佐料与y克第二种规格的佐料混合得第三种规格的佐料, 则有下表 表1 混合后四种原料的含量 原料 佐料规格 A B C D 第一种 x x x x 第二种 y y y y 第三种 (x + y) (x + y) (x + y) (x + y) 因而有如下线性方程组 (*) 【模型检验】把x = 7, y = 12代入上述方程组(*), 则各等式都成立. 可见模型假设中的第三个假设不影响解的正确性. Matlab实验题 蛋白质、碳水化合物和脂肪是人体每日必须的三种营养, 但过量的脂肪摄入不利于健康.人们可以通过适量的运动来消耗多余的脂肪. 设三种食物(脱脂牛奶、大豆面粉、乳清)每100克中蛋白质、碳水化合物和脂肪的含量以及慢跑5分钟消耗蛋白质、碳水化合物和脂肪的量如下表. 表2 三种食物的营养成分和慢跑的消耗情况 营养 每100克食物所含营养(克) 慢跑5分钟消耗量(克) 每日需要的营养量(克) 牛奶 大豆面粉 乳清 蛋白质 36 51 13 10 33 碳水化合物 52 34 74 20 45 脂肪 10 7 1 15 3 问怎样安排饮食和运动才能实现每日的营养需求? 案例三. 投入产出问题 在研究多个经济部门之间的投入产出关系时, W. Leontief提出了投入产出模型. 这为经济学研究提供了强有力的手段. W. Leontief因此获得了1973年的Nobel经济学奖. 图7 三个经济部门 这里暂时只讨论一个简单的情形. 【模型准备】某地有一座煤矿, 一个发电厂和一条铁路. 经成本核算, 每生产价值1元钱的煤需消耗0.3元的电; 为了把这1元钱的煤运出去需花费0.2元的运费; 每生产1元的电需0.6元的煤作燃料; 为了运行电厂的辅助设备需消耗本身0.1元的电, 还需要花费0.1元的运费; 作为铁路局, 每提供1元运费的运输需消耗0.5元的煤, 辅助设备要消耗0.1元的电. 现煤矿接到外地6万元煤的订货, 电厂有10万元电的外地需求, 问: 煤矿和电厂各生产多少才能满足需求? 【模型假设】假设不考虑价格变动等其他因素. 【模型建立】设煤矿, 电厂, 铁路分别产出x元, y元, z元刚好满足需求. 则有下表 表3 消耗与产出情况 产出(1元) 产出 消耗 订单 煤 电 运 消耗 煤 0 0.6 0.5 x 0.6y + 0.5z 60000 电 0.3 0.1 0.1 y 0.3x + 0.1y + 0.1z 100000 运 0.2 0.1 0 z 0.2x + 0.1y 0 根据需求, 应该有 , 即 【模型求解】在Matlab命令窗口输入以下命令 >> A = [1,-0.6,-0.5;-0.3,0.9,-0.1;-0.2,-0.1,1]; b = [60000;100000;0]; >> x = A\b Matlab执行后得 x = 1.0e+005 * 1.9966 1.8415 0.5835 可见煤矿要生产1.9966105元的煤, 电厂要生产1.8415105元的电恰好满足需求. 【模型分析】令x =, A =, b =, 其中x称为总产值列向量, A称为消耗系数矩阵, b称为最终产品向量, 则 Ax == 根据需求, 应该有x - Ax = b, 即(E - A)x = b. 故x = (E - A)-1b. Matlab实验题 某乡镇有甲、乙、丙三个企业. 甲企业每生产1元的产品要消耗0.25元乙企业的产品和0.25元丙企业的产品. 乙企业每生产1元的产品要消耗0.65元甲企业的产品, 0.05元自产的产品和0.05元丙企业的产品. 丙企业每生产1元的产品要消耗0.5元甲企业的产品和0.1元乙企业的产品. 在一个生产周期内, 甲、乙、丙三个企业生产的产品价值分别为100万元, 120万元, 60万元, 同时各自的固定资产折旧分别为20万元, 5万元和5万元. (1) 求一个生产周期内这三个企业扣除消耗和折旧后的新创价值. (2) 如果这三个企业接到外来订单分别为50万元, 60万元, 40万元, 那么他们各生产多少才能满足需求? 案例四. 平板的稳态温度分布问题 在热传导的研究中, 一个重要的问题是确定一块平板的稳态温度分布. 根据…定律, 只要测定一块矩形平板四周的温度就可以确定平板上各点的温度. 图8 一块平板的温度分布图 【模型准备】如图9所示的平板代表一条金属梁的截面. 已知四周8个节点处的温度(单位C), 求中间4个点处的温度T1, T2, T3, T4. T1 T2 T3 T4 100 80 90 80 60 50 60 50 图9 一块平板的温度分布图 【模型假设】假设忽略垂直于该截面方向上的热传导, 并且每个节点的温度等于与它相邻的四个节点温度的平均值. 【模型建立】根据已知条件和上述假设, 有如下线性方程组 【模型求解】将上述线性方程组整理得 . 在Matlab命令窗口输入以下命令 >> A = [4,-1,-1,0;-1,4,0,-1;-1,0,4,-1;0,-1,-1,4]; b = [190;140;140;100]; >> x = A\b; x’ Matlab执行后得 ans = 82.9167 70.8333 70.8333 60.4167 可见T1 = 82.9167, T2 = 70.8333, T3 = 70.8333, T4 = 60.4167. 参考文献 陈怀琛, 高淑萍, 杨威, 工程线性代数, 北京: 电子工业出版社, 2007. 页码: 15-16. Matlab实验题 假定下图中的平板代表一条金属梁的截面, 并忽略垂直于该截面方向上的热传导. 已知平板内部有30个节点, 每个节点的温度近似等于与它相邻的四个节点温度的平均值. 设4条边界上的温度分别等于每位同学学号的后四位的5倍, 例如学号为16308209的同学计算本题时, 选择Tl = 40, Tu = 10, Tr = 0, Td = 45. Tu T1 T5 Tl Tl Td T2 T6 T7 T10 Tr Tr Tu T26 T30 Td T27 Tu Tr Td Tl 图10 一块平板的温度分布图 (1) 建立可以确定平板内节点温度的线性方程组. (2) 用Matlab软件求解该线性方程组. (3) 用Matlab中的函数mesh绘制三维平板温度分布图. 案例五. CT图像的代数重建问题 X射线透视可以得到3维对象在2维平面上的投影, CT则通过不同角度的X射线得到3维对象的多个2维投影, 并以此重建对象内部的3维图像. 代数重建方法就是从这些2维投影出发, 通过求解超定线性方程组, 获得对象内部3维图像的方法. 图11双层螺旋CT 图12 CT图像 这里我们考虑一个更简单的模型, 从2维图像的1维投影重建原先的2维图像. 一个长方形图像可以用一个横竖均匀划分的离散网格来覆盖, 每个网格对应一个像素, 它是该网格上各点像素的均值. 这样一个图像就可以用一个矩阵表示,其元素就是图像在一点的灰度值(黑白图像). 下面我们以33图像为例来说明. 表4 消耗与产出情况 33图像 各点的灰度值 水平方向上 的叠加值 x1 = 1 x2 = 0 x3 = 0 x1 + x2 + x3 = 1 x4 = 0 x5 = 0.5 x6 = 0.5 x4 + x5 + x6 = 1 x7 = 0.5 x8 = 0 x9 = 1 x7 + x8 + x9 = 1.5 竖直方向上的叠加值 x1 + x4 + x7 = 1.5 x2 + x5 + x8 = 0.5 x3 + x6 + x9 = 1.5 每个网格中的数字xi代表其灰度值, 范围在[0, 1]内. 0表示白色, 1表示黑色, 0.5表示灰色. 如果我们不知道网格中的数值, 只知道沿竖直方向和水平方向的叠加值, 为了确定网格中的灰度值, 可以建立线性方程组(含有6个方程, 9个未知数) 显然该方程组的解是不唯一的, 为了重建图像, 必须增加叠加值. 如我们增加从右上方到左下方的叠加值, 则方程组将增加5个方程 x1 = 1, x2 + x4 = 0, x3 + x5 + x7 = 1, x6 + x8 = 0.5, x9 = 1, 和上面的6个方程放在一起构成一个含有11个方程, 9个未知数的线性方程组. 【模型准备】设33图像中第一行3个点的灰度值依次为x1, x2, x3, 第二行3个点的灰度值依次为x4, x5, x6, 第三行3个点的灰度值依次为x7, x8, x9. 沿竖直方向的叠加值依次为1.5, 0.5, 1.5, 沿水平方向的叠加值依次为1, 1, 1.5, 沿右上方到左下方的叠加值依次为1, 0, 1, 0.5, 1. 确定x1, x2, …, x9的值. 【模型建立】由已知条件可得(含有11个方程, 9个未知数的)线性方程组 【模型求解】在Matlab命令窗口输入以下命令 >> A = [1,1,1,0,0,0,0,0,0;0,0,0,1,1,1,0,0,0;0,0,0,0,0,0,1,1,1; 1,0,0,1,0,0,1,0,0;0,1,0,0,1,0,0,1,0;0,0,1,0,0,1,0,0,1; 1,0,0,0,0,0,0,0,0;0,1,0,1,0,0,0,0,0;0,0,1,0,1,0,1,0,0; 0,0,0,0,0,1,0,1,0;0,0,0,0,0,0,0,0,1]; >> b = [1;1;1.5;1.5;0.5;1.5;1;0;1;0.5;1]; >> x = A\b; x’ Matlab执行后得 Warning: Rank deficient, rank = 8 tol = 4.2305e-015. ans = 1.0000 0.0000 0 -0.0000 0.5000 0.5000 0.5000 -0.0000 1.0000 可见上述方程组的解不唯一. 其中的一个特解为 x1 = 1, x2 = 0, x3 = 0, x4 = 0, x5 = 0.5, x6 = 0.5, x7 = 0.5, x8 = 0, x9 = 1. 【模型分析】上述结果表明, 仅有三个方向上的叠加值还不够.可以再增加从左上方到右下方的叠加值. 在实际情况下, 由于测量误差, 上述线性方程组可能是超定的. 这时可以将超定方程组的近似解作为重建的图像数据. Matlab实验题 给定一个33图像的2个方向上的灰度叠加值: 沿左上方到右下方的灰度叠加值依次为0.8, 1.2, 1.7, 0.2, 0.3; 沿右上方到左下方的灰度叠加值依次为0.6, 0.2, 1.6, 1.2, 0.6. (1) 建立可以确定网格数据的线性方程组, 并用Matlab求解. (2) 将网格数据乘以256, 再取整, 用Matlab绘制该灰度图像. 案例六. 平衡结构的梁受力计算 在桥梁、房顶、铁塔等建筑结构中, 涉及到各种各样的梁. 对这些梁进行受力分析是设计师、工程师经常做的事情. 图13埃菲尔铁塔全景 图14 埃菲尔铁塔局部 下面以双杆系统的受力分析为例, 说明如何研究梁上各铰接点处的受力情况. 【模型准备】在图15所示的双杆系统中, 已知杆1重G1 = 200牛顿, 长L1 = 2米, 与水平方向的夹角为q1 = p/6, 杆2重G2 = 100牛顿, 长L2 = 米, 与水平方向的夹角为q2 = p/4. 三个铰接点A, B, C所在平面垂直于水平面. 求杆1, 杆2在铰接点处所受到的力. A B C 杆1 杆2 p/6 p/4 图15双杆系统 【模型假设】假设两杆都是均匀的. 在铰接点处的受力情况如图16所示. 【模型建立】对于杆1: 水平方向受到的合力为零, 故N1 = N3, 竖直方向受到的合力为零, 故N2 + N4 = G1, 以点A为支点的合力矩为零, 故(L1sinq1)N3 + (L1cosq1)N4 = (L1cosq1)G1. A B C 杆1 杆2 C N1 N2 N4 N3 N7 N8 N5 N6 G1 G2 图16 两杆受力情况 对于杆2类似地有 N5 = N7, N6 = N8 + G2, (L2sinq2)N7 = (L2cosq2)N8 + (L2cosq2)G2. 此外还有N3 = N7, N4 = N8. 于是将上述8个等式联立起来得到关于N1, N2, …, N8的线性方程组: 【模型求解】在Matlab命令窗口输入以下命令 >> G1=200; L1=2; theta1=pi/6; G2=100; L2=sqrt(2); theta2=pi/4; >> A = [1,0,-1,0,0,0,0,0;0,1,0,1,0,0,0,0; 0,0,L1*sin(theta1),L1*cos(theta1),0,0,0,0;0,0,0,0,1,0,-1,0; 0,0,0,0,0,1,0,-1;0,0,0,0,0,0,L2*sin(theta2),-L2*cos(theta2); 0,0,1,0,0,0,-1,0;0,0,0,1,0,0,0,-1]; >> b = [0;G1;0.5*L1*cos(theta1)*G1;0;G2;0.5*L2*cos(theta2)*G2;0;0]; >> x = A\b; x’ Matlab执行后得 ans = 95.0962 154.9038 95.0962 45.0962 95.0962 145.0962 95.0962 45.0962 【模型分析】最后的结果没有出现负值, 说明图16中假设的各个力的方向与事实一致. 如果结果中出现负值, 则说明该力的方向与假设的方向相反. 参考文献 陈怀琛, 高淑萍, 杨威, 工程线性代数, 北京: 电子工业出版社, 2007. 页码: 157- 158. Matlab实验题 有一个平面结构如下所示, 有13条梁(图中标号的线段)和8个铰接点(图中标号的圈)联结在一起. 其中1号铰接点完全固定, 8号铰接点竖直方向固定, 并在2号, 5号和6号铰接点上, 分别有图示的10吨, 15吨和20吨的负载. 在静平衡的条件下,任何一个铰接点上水平和竖直方向受力都是平衡的. 已知每条斜梁的角度都是45. (1) 列出由各铰接点处受力平衡方程构成的线性方程组. (2) 用Matlab软件求解该线性方程组, 确定每条梁受力情况. 图17 一个平面结构的梁 案例七. 化学方程式配平问题 在用化学方法处理污水过程中, 有时会涉及到复杂的化学反应. 这些反应的化学方程式是分析计算和工艺设计的重要依据. 在定性地检测出反应物和生成物之后,可以通过求解线性方程组配平化学方程式. 图18 污水处理 【模型准备】某厂废水中含KCN, 其浓度为650mg/L. 现用氯氧化法处理, 发生如下反应: KCN + 2KOH + Cl2 = KOCN + 2KCl + H2O. 投入过量液氯, 可将氰酸盐进一步氧化为氮气. 请配平下列化学方程式: KOCN + KOH + Cl2 === CO2 + N2 + KCl + H2O. (注: 题目摘自福建省厦门外国语学校2008-2009学年高三第三次月考化学试卷) 【模型建立】设 x1KOCN + x2KOH + x3Cl2 === x4CO2 + x5N2 + x6KCl + x7H2O, 则 , 即 【模型求解】在Matlab命令窗口输入以下命令 >> A = [1,1,0,0,0,-1,0;1,1,0,-2,0,0,-1;1,0,0,-1,0,0,0; 1,0,0,0,-2,0,0;0,1,0,0,0,0,-2;0,0,2,0,0,-1,0]; >> x = null(A,’r’); format rat, x’ Matlab执行后得 ans = 1 2 3/2 1 1/2 3 1 可见上述齐次线性方程组的通解为 x = k(1, 2, 3/2, 1, 1/2, 3, 1)T. 取k = 2得x = (2, 4, 3, 2, 1, 6, 2)T. 可见配平后的化学方程式如下 2KOCN + 4KOH + 3Cl2 === 2CO2 + N2 + 6KCl + 2H2O. 【模型分析】利用线性方程组配平化学方程式是一种待定系数法. 关键是根据化学方程式两边所涉及到的各种元素的量相等的原则列出方程. 所得到的齐次线性方程组Ax = q中所含方程的个数等于化学方程式中元素的种数s, 未知数的个数就是化学方程式中的项数n. 当r(A) = n -1时, Ax = q的基础解系中含有1个(线性无关的)解向量. 这时在通解中取常数k为各分量分母的最小公倍数即可. 例如本例中 1, 2, 3/2, 1, 1/2, 3, 1 分母的最小公倍数为2, 故取k = 2. 当r(A) n -2时, Ax = q的基础解系中含有2个以上的线性无关的解向量. 这时可以根据化学方程式中元素的化合价的上升与下降的情况, 在原线性方程组中添加新的方程. 参考文献 陈怀琛, 高淑萍, 杨威, 工程线性代数, 北京: 电子工业出版社, 2007. 页码: 84-85. Matlab实验题 配平下列反应式 (1) FeS + KMnO4 + H2SO4 —— K2SO4 + MnSO4 + Fe2(SO4)3 + H2O + S↓ (2) Al2(SO4)3 + Na2CO3 + H2O —— Al(OH)3↓+ CO2↑+ Na2SO4 案例八. 互付工资问题 互付工资问题是多方合作相互提供劳动过程中产生的. 比如农忙季节, 多户农民组成互助组, 共同完成各户的耕、种、收等农活. 又如木工, 电工, 油漆工等组成互助组, 共同完成各家的装潢工作. 由于不同工种的劳动量有所不同, 为了均衡各方的利益, 就要计算互付工资的标准. 图19 农忙互助 图20 装修互助 【模型准备】现有一个木工, 电工, 油漆工. 相互装修他们的房子, 他们有如下协议: (1) 每人工作10天(包括在自己家的日子), (2) 每人的日工资一般的市价在60~80元之间, (3) 日工资数应使每人的总收入和总支出相等. 表5 工作天数 在谁家 工人 木工 电工 油漆工 木工家 2 1 6 电工家 4 5 1 油漆工家 4 4 3 求每人的日工资. 【模型假设】假设每人每天工作时间长度相同. 无论谁在谁家干活都按正常情况工作, 既不偷懒, 也不加班. 【模型建立】设木工, 电工, 油漆工的日工资分别为x, y, z元, 则由下表 表6 各家应付工资和各人应得收入 在谁家 工人 木工 电工 油漆工 各家应付工资 木工家 2x 1y 6z 2x + y + 6z 电工家 4x 5y 1z 4x + 5y + z 油漆工家 4x 4y 3z 4x + 4y + 3z 各人应得收入 10x 10y 10z 可得 , 即 【模型求解】在Matlab命令窗口输入以下命令 >> A = [-8,1,6;4,-5,1;4,4,-7]; >> x = null(A,’r’); format rat, x’ Matlab执行后得 ans = 31/36 8/9 1 可见上述齐次线性方程组的通解为x = k(31/36, 8/9, 1)T. 因而根据“每人的日工资一般的市价在60~80元之间”可知 60 k > A = [1,-0.4,-0.6;-0.6,0.9,-0.2;-0.4,-0.5,0.8]; >> x = null(A,’r’); format short, x’ Matlab执行后得 ans = 0.9394 0.8485 1.0000 可见上述齐次线性方程组的通解为 x = k(0.9394, 0.8485, 1)T. 这就是说, 如果煤炭、电力、钢铁行业每年总产出的价格分别0.9394亿元, 0.8485亿元, 1亿元, 那么每个行业的投入与产出都相等. 【模型分析】实际上, 一个比较完整的经济系统不可能只涉及三个行业, 因此需要统计更多的行业间的分配数据. Matlab实验题 假设一个经济系统由煤炭、石油、电力、钢铁、机械制造、运输行业组成, 每个行业的产出在各个行业中的分配如下表所示: 表8 行业产出分配表 产出分配 购买者 煤炭 石油 电力 钢铁 制造 运输 0 0 0.2 0.1 0.2 0.2 煤炭 0 0 0.1 0.1 0.2 0.1 石油 0.5 0.1 0.1 0.2 0.1 0.1 电力 0.4 0.1 0.2 0 0.1 0.4 钢铁 0 0.1 0.3 0.6 0 0.2 制造 0.1 0.7 0.1 0 0.4 0 运输 每一列中的元素表示占该行业总产出的比例. 求使得每个行业的投入与产出都相等的平衡价格. 参考文献 David C. Lay, 线性代数及其应用, 沈复兴, 傅莺莺等译, 北京: 人民邮电出版社, 2009. 页码: 49-50. 案例十. 电路设计问题 电路是电子元件的神经系统. 参数的计算是电路设计的重要环节. 其依据来自两个方面: 一是客观需要, 二是物理学定律. 图22 USB扩展板 【模型准备】假设图23中的方框代表某类具有输入和输出终端的电路. 用记录输入电压和输入电流(电压v以伏特为单位, 电流i以安培为单位), 用记录输出电压和输入电流. 若= A, 则称矩阵A为转移矩阵. 输入终端v1 输出终端v2 i1 i2 电路 图23 具有输入和输出终端的电子电路图 图24给出了一个梯形网络, 左边的电路称为串联电路, 电阻为R1(单位: 欧姆). 右边的电路是并联电路, 电路R2. 利用欧姆定理和楚列斯基定律, 我们可以得到串联电路和并联电路的转移矩阵分别是 和 v1 v2 i1 i2 R1 v3 i2 i3 R2 串联电路 并联电路 图24 梯形网络 设计一个梯形网络, 其转移矩阵是. 【模型假设】假设导线的电阻为零. 【模型建立】设A1和A2分别是串联电路和并联电路的转移矩阵, 则输入向量x先变换成A1x, 再变换到A2(A1x). 其中 A2A1 == 就是图22中梯形网络的转移矩阵. 于是, 原问题转化为求R1, R2的值使得=. 【模型求解】由=可得. 根据其中的前两个方程可得R1 = 8, R2 = 2. 把R1 = 8, R2 = 2代入上面的第三个方程确实能使等式成立. 这就是说在图22中梯形网络中取R1 = 8, R2 = 2即为所求. 【模型分析】若要求的转移矩阵改为, 则上面的梯形网络无法实现. 因为这时对应的方程组是. 根据前两个方程依然得到R1 = 8, R2 = 2, 但把R1 = 8, R2 = 2代入上第三个方程却不能使等式成立. 参考文献 David C. Lay, 线性代数及其应用, 沈复兴, 傅莺莺等译, 北京: 人民邮电出版社, 2009. 页码: 129-130. 练习题 根据基尔霍夫回路电路定律(各节点处流入和流出的电流强度的代数和为零, 各回路中各支路的电压降之和为零), 列出下图所示电路中电流i1, i2, i3所满足的线性方程组, 并用矩阵形式表示: E1 E2 R1 R2 R3 R4 R5 i1 i2 i3 ① ② ③ 图25 简单的回路 案例十一. 平面图形的几何变换 随着计算机科学技术的发展, 计算机图形学的应用领域越来越广, 如仿真设计、效果图制作、动画片制作、电子游戏开发等. 图26 计算机图形学的广泛应用 图形的几何变换, 包括图形的平移、旋转、放缩等, 是计算机图形学中经常遇到的问题. 这里暂时只讨论平面图形的几何变换. 【模型准备】平面图形的旋转和放缩都很容易用矩阵乘法实现, 但是图形的平移并不是线性运算, 不能直接用矩阵乘法表示. 现在要求用一种方法使平移、旋转、放缩能统一用矩阵乘法来实现. 【模型假设】设平移变换为 (x, y) (x+a, y+b) 旋转变换(绕原点逆时针旋转q角度)为 (x, y) (xcosq - ysinq, xsinq + ycosq) 放缩变换(沿x轴方向放大s倍, 沿y轴方向放大t倍)为 (x, y) (sx, ty) 【模型求解】R2中的每个点(x, y)可以对应于R3中的(x, y, 1). 它在xOy平面上方1单位的平面上. 我们称(x, y, 1)是(x, y)的齐次坐标. 在齐次坐标下, 平移变换 (x, y) (x+a, y+b) 可以用齐次坐标写成 (x, y, 1) (x+a, y+b, 1). 于是可以用矩阵乘积=实现. 旋转变换 (x, y) (xcosq - ysinq, xsinq + ycosq) 可以用齐次坐标写成 (x, y, 1) (xcosq - ysinq, xsinq + ycosq, 1).
展开阅读全文
提示  淘文阁 - 分享文档赚钱的网站所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
关于本文
本文标题:数学建模案例解析分析线性代数建模案例解析(20例).doc
链接地址:https://www.taowenge.com/p-2604185.html
关于淘文阁 - 版权申诉 - 用户使用规则 - 积分规则 - 联系我们

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号 © 2020-2023 www.taowenge.com 淘文阁 

收起
展开