线性代数在数学建模中的应用.docVIP

Word格式 完美整理 第六章 线性代数模型 §6.1 Matlab求解线性代数工具简介 　　１．矩阵的秩. 　　 rref?或?rrefmovie? 格式??R?=?rref(A)?????%用高斯—约当消元法和行主元法求?A 的行最简行矩阵R.?rrefmovie(A)??????%给出每一步化简的过程.? 　　２．方阵的行列式：det(A) 　　３．逆矩阵：inv(A) 　　指令inv(A)给出方阵A的逆矩阵，如果A不可逆，则inv(A)给出的矩阵的元素都是Inf． 　　利用初等变换也可以求出逆矩阵，构造n行2n列的矩阵(A E)，并进行行初等变换，当把A变为单位矩阵时，E就变成了A的逆矩阵．利用matlab命令rref可以求出矩阵的行简化阶梯形．输入命令： D=[A,eye(3)] D = 1 2 3 1 0 0 2 2 1 0 1 0 3 4 3 0 0 1 rref(D) ans = 1.0000 0 0 1.0000 3.0000 -2.0000 0 1.0000 0 -1.5000 -3.0000 2.5000 0 0 1.0000 1.0000 1.0000 -1.0000 　　线性方程组的求解是用矩阵除来完成的,，当且可逆时，给出唯一解．这时矩阵除相当于；当时，矩阵除给出方程的最小二乘解；当时，矩阵除给出方程的最小范数解． 例6.1：解方程组： 解：输入命令： a=[1 -1 1 2;1 1 -2 1;1 1 1 0;1 0 1 -1]; b=[1;1;2;1]; x=a\b x = 0.8333 0.7500 0.4167 0.2500 输入命令： z=inv(a)*b z = 0.8333 0.7500 0.4167 0.2500 例6.2：解方程组： 解：方程的个数和未知数不相等，用消去法，将增广矩阵化为行简化阶梯形，如果系数矩阵的秩不等于增广矩阵的秩，则方程组无解；如果系数矩阵的秩等于增广矩阵的秩，则方程组有解，方程组的解就是行简化阶梯形所对应的方程组的解．输入命令： a=[2 1 1 -1 -2 2;1 -1 2 1 -1 4;2 -3 4 3 -1 8]; rref(a) ans = 1 0 0 0 0 0 0 1 0 -1 -1 0 0 0 1 0 -1 2 由结果看出，,为自由未知量，方程组的解为： 　　　　　　 　　　　　　 　　　　　　 例６：解方程组： 　解：输入命令： a=[1 -1 -1 1;1 -1 1 -3;1 -1 0 -1;1 -1 -2 3]; rref(a) ans = 1 -1 0 -1 0 0 1 -2 0 0 0 0 0 0 0 0 由结果看出，,为自由未知量，方程组的解为： §6.3 交通流量模型 城市道路网中每条道路、每个交叉路口的车流量调查，是分析、评价及改善城市交通状况的基础。根据实际车流量信息可以设计流量控制方案，必要时设置单行线，以免大量车辆长时间拥堵。下面考虑一简化的交通流量统计问题。 问题描述：某城市有下图所示的交通图，每条道路都是单行线，需要调查每条道路每小时的车流量。图中的数字表示该条路段的车流数。如果每个交叉路口进入和离开的车数相等，整个图中进入和离开的车数相等。试解决以下问题： (1)建立确定每条道路流量的线性方程组； (2)分析哪些流量数据是多余的； (3)为了唯一确定未知流量, 需要增添哪几条道路的流量统计。 300 300 500 150 180 350 160 220 300 100 290 400 150 x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 图6.1 某城市单行线车流量 模型假设：根据题目要求，作以下假设：（1）全