实验指导7 线性方程组的应用.docVIP

实验指导7 线性方程组的应用 下图给出了某城市部分单行街道的交通流量（每小时通过的车辆数） 图中有6个路口，已有9条街道记录了当天的平均车流量。另有7处的平均车流量未知，试利用每个路口的进出车流量相等关系推算这7处的平均车流量。 1.问题分析与数学模型 在图1中的任何一个路口（十字路口或丁字路口）处，都有车辆流进和流出。当一天结束后，流进的车辆数和流出的车辆数应该相等以达到平衡。在图中有的街道车流量有数据记录，而有的没有数据记录。我们可以理解为有数据记录的街道有专人（或设备）记录了当天的车流量情况，而没有记录的街道是由于人力不足（或设备的经费还没到位）造成的。为了填补空白，在没有数据记录的街道处假设车流量是未知数。在每一个路口处可根据进出的车流量相等关系，建立一个线性代数方程。图1中有六个路口，可建立含六个方程的线性方程组。问题的答案应该是在所列的线性方程组的通解中去寻找。将方程组写成矩阵向量形式为AX = b其中 X = [x1 x2 x3 x4 x5 x6 x7]T 显然，在这一方程组中，未知数个数多于方程的个数，是不定方程组。所以，当方程组的系数矩阵A的秩与增广矩阵[A b]的秩相等时，该问题有无穷多组解。由于图1中街道都是单行道，每一街道上的车流量只能是正数或者是零。故应在方程组的解集合中寻找非负解，如果方程组没有解或者没有非负解，则说明问题所给的数据有误。求解问题分三步：第一步，判断方程组是否有解；第二步，如果有解则求出方程组的通解；第三步，在通解中找非负特解。 2 . 程序和计算结果 在MATLAB环境中，首先输入方程组的系数矩阵A和方程组右端向量b A=[1 0 1 0 0 0 0；1 –1 0 1 0 0 0；0 1 0 0 –1 0 0；0 0 1 0 0 1 0；0 0 0 1 0 1 –1；0 0 0 0 1 0 -1] b = [700；200；200；500；0；-200] 然后用命令 r1=rank(A) r2=rank([A b])计算系数矩阵的秩r1，和增广矩阵[A b]的秩r2，得 r1= 5 r2 = 5 这说明系数矩阵和增广矩阵的秩相等。最后，用命令 rref([A b])将增广矩阵化为最简行阶梯形矩阵，得数据 ans =??? 1?? 0?? 0?? 0?? 0?? -1?? 0?? 200??? 0?? 1?? 0?? 0?? 0?? 0?? -1?? 0??? 0?? 0?? 1?? 0?? 0?? 1?? 0??? 500??? 0?? 0?? 0?? 1?? 0?? 1?? -1?? 0??? 0?? 0?? 0?? 0?? 1?? 0?? -1?? -200??? 0?? 0?? 0?? 0?? 0?? 0?? 0??? 0由此可确定对应的齐次方程组的基础解系以及非齐次方程组的通解。由于增广矩阵的秩为5，而方程组含有7个未知数x1，x2，x3，x4，x5，x6，x7，故方程组的通解中含有两个自由未知数。最后的数据结果是与增广矩阵等价的最简行阶梯形矩阵。由最简行阶梯形矩阵，可得与原方程组等价的简化后的方程组取x6，x7为自由未知数，直接可得原方程组的通解形式 3 . 由上面所得的方程组通解表达式，取适当的k1和k2使特解为非负数，即得一组满足问题条件的解。取k-1 = 0，k2 = 200，得 [ x1 x2 x3 x4 x5 x6 x7] = [200 200 500 200 0 0 200] 将对应数据填写，得下图 显然，这一问题的解是不唯一的。 ，