matlab交通问题求解,交通流量问题MATLAB求解.doc

交通流量问题MATLAB求解

交通流量问题 MATLAB

下图给出了某城市部分单行街道的交通流量(每小时通过的车辆数)

图 1 中有 6 个路口，已有 9 条街道记录了当天的平均车流量。另

有 7 处的平均车流量未知，试利用每个路口的进出车流量相等关系推

算这 7 处的平均车流量。

问题分析与数学模型

在图 1 中的任何一个路口(十字路口或丁字路口)处，都有车辆

流进和流出。当一天结束后，流进的车辆数和流出的车辆数应该相等

以达到平衡。在图中有的街道车流量有数据记录，而有的没有数据记

录。我们可以理解为有数据记录的街道有专人(或设备)记录了当天

的车流量情况，而没有记录的街道是由于人力不足(或设备的经费还

没到位)造成的。为了填补空白，在没有数据记录的街道处假设车流

量是未知数。在每一个路口处可根据进出的车流量相等关系，建立一

个线性代数方程。图 1 中有六个路口，可建立含六个方程的线性方程

组。问题的答案应该是在所列的线性方程组

? ? ?

x x 700

1 3

? ? ? ?

x x x 200

?

1 2 4

? ? ?

?x x 200

2 5

? ? ?

x x 500

?

3 6

? ? ? ?

x x x 0

?

4 6 7

?x ? x ? ?200

?

5 7

的通解中去寻找。将方程组写成矩阵向量形式为 AX = 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

? ?，x ??x ?

1

? ?

x

? ?

2

? ?

? ?

?x ?

? ?

72 . 程序和计算结果

在 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

这说明系数矩阵和增广矩阵的秩相等。最后，用命令

[R,ID]=rref([A b])

将增广矩阵化为最简行阶梯形矩阵，得数据

R =

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

ID = 1 2 3 4 5

由此可确定对应的齐次方程组的基础解系以及非齐次方程组的

通解。取 x6 , x7 为自由未知数，直接可得原方程组的通解形式

? ? ? ? ? ? ? ?

x

1 0 200

1

? ? ? ? ? ? ? ?

x 0 1 0

? ? ? ? ? ? ? ?

2

? ? ?? ? ? ? ? ?

x 1 0 500

? ? ? ? ? ? ? ?

3

?x ? c ? ? c ? ? ? ?

? ?1 ? 1 ? 0

4 1 2

? ? ? ? ? ? ?? ?

x 0 1 200

? ? ? ? ? ? ? ?

5

? ? ? ? ? ? ? ?

x 1 0 0

6

? ? ? ? ? ? ? ?

0 1 0

x ? ? ? ? ? ?

? ?

7

3 .问题解答

由上面所得的方程组通解表达式，取适当的c 和

1c 使特解为非

2负数，即得一组满足问题条件的解。例如，取 1 0, 2 200

c ? c ? ，得

x ? [200 200 500 200 0 0 200]T

将对应数据填写，得下图

显然，这一问题的解是不唯一的。