交通规划——容量限制增量加载计算机程序代码.doc

分配率K=10%Private Sub Command1_Click()Form2.HideForm1.ShowEnd SubPrivate Sub Command2_Click()Form2.HideForm3.ShowEnd SubDim s As StringDim t0(1 To 12) As DoubleDim f As DoublePrivate Sub Command1_Click()Dim v, c, e As IntegerDim t(9) As Doublev = InputBox(请输入路段的设计速度, , 50)c = InputBox(请输入路网外框路段的容量, , 750)Text1.Text = (Val(Text1.Text) / v) * (1 + 0.15 * (q / c) 4)Text2.Text = (Val(Text2.Text) / v) * (1 + 0.15 * (q / c) 4)Text3.Text = (Val(Text3.Text) / v) * (1 + 0.15 * (q / c) 4)Text4.Text = (Val(Text4.Text) / v) * (1 + 0.15 * (q / c) 4)Text5.Text = (Val(Text5.Text) / v) * (1 + 0.15 * (q / c) 4)Text6.Text = (Val(Text6.Text) / v) * (1 + 0.15 * (q / c) 4)Text7.Text = (Val(Text7.Text) / v) * (1 + 0.15 * (q / c) 4)Text8.Text = (Val(Text8.Text) / v) * (1 + 0.15 * (q / c) 4)Text9.Text = (Val(Text9.Text) / v) * (1 + 0.15 * (q / c) 4)Text10.Text = (Val(Text10.Text) / v) * (1 + 0.15 * (q / c) 4)Text11.Text = (Val(Text11.Text) / v) * (1 + 0.15 * (q / c) 4)Text12.Text = (Val(Text12.Text) / v) * (1 + 0.15 * (q / c) 4)t0(1) = Val(Text1.Text)t0(2) = Val(Text2.Text)t0(3) = Val(Text3.Text)t0(4) = Val(Text4.Text)t0(5) = Val(Text5.Text)t0(6) = Val(Text6.Text)t0(7) = Val(Text7.Text)t0(8) = Val(Text8.Text)t0(9) = Val(Text9.Text)t0(10) = Val(Text10.Text)t0(11) = Val(Text11.Text)t0(12) = Val(Text12.Text)End SubPrivate Sub Command10_Click()Text29.Text = : Text30.Text = : Text31.Text = : Text32.Text = Text33.Text = : Text34.Text = : Text35.Text = : Text36.Text = Text37.Text = : Text38.Text = : Text39.Text = : Text40.Text = Text41.Text = : Text42.Text = : Text43.Text = : Text44.Text = Text45.Text = : Text46.Text = : Text47.Text = : Text48.Text = Text49.Text = : Text50.Text = : Text51.Text = : Text52.Text = Dim X, Y, z, s, m, n, o As Integerz = InputBox(请输入第三次OD量分配率, , 0.15)n = 0.55T123 = Val(Text1.Text) + Val(Text2.Text)T14523 = Val(Text3.Text) + Val(Text6.Text) + Val(Text4.Text) + Val(Text2.Text)T14563 = Val(Text3.Text) + Val(Text6.Text) + Val(Text7.Text) + Val(Text5.Text)T12369 = Val(Text1.Text) + Val(Text2.Text) + Val(Text5.Text) + Val(Text10.Text)T12569 = Val(Text1.Text) + Val(Text4.Text) + Val(Text7.Text) + Val(Text10.Text)T12589 = Val(Text1.Text) + Val(Text4.Text) + Val(Text9.Text) + Val(Text12.Text)T14569 = Val(Text3.Text) + Val(Text6.Text) + Val(Text7.Text) + Val(Text10.Text)T14589 = Val(Text3.Text) + Val(Text6.Text) + Val(Text8.Text) + Val(Text12.Text)T14789 = Val(Text3.Text) + Val(Text8.Text) + Val(Text11.Text) + Val(Text12.Text)q = 750q1 = 750 A-BIf T123 T14523 And T123 T14563 Then If Val(Text14.Text * n) = q Then Text29.Text = Text14.Text * n Text31.Text = Text14.Text * n Else Text29.Text = q Text31.Text = q End IfEnd If A-CIf Val(Text3.Text) + Val(Text8.Text) Val(Text1.Text) + Val(Text6.Text) + Val(Text4.Text) + Val(Text8.Text) And Val(Text1.Text) + Val(Text2.Text) Val(Text1.Text) + Val(Text4.Text) + Val(Text9.Text) + Val(Text11.Text) Then If Val(Text15.Text * n) = q Then Text33.Text = Text15.Text * n Text43.Text = Text15.Text * n Else Text33.Text = q Text43.Text = q End IfEnd If 由于城市道路网四周网阻抗相同 C - A, B - A, B - D, D - B, c - D, D - c分配流量方法相同Text30.Text = IIf(Val(Text17.Text * n) = q, Val(Text17.Text * n), q)Text32.Text = IIf(Val(Text17.Text * n) = q, Val(Text17.Text * n), q)Text37.Text = IIf(Val(Text26.Text * n) = q, Val(Text26.Text * n), q)Text47.Text = IIf(Val(Text26.Text * n) = q, Val(Text26.Text * n), q)Text50.Text = IIf(Val(Text24.Text * n) = q, Val(Text24.Text * n), q)Text52.Text = IIf(Val(Text24.Text * n) = q, Val(Text24.Text * n), q)Text49.Text = IIf(Val(Text27.Text * n) = q, Val(Text27.Text * n), q)Text51.Text = IIf(Val(Text27.Text * n) = q, Val(Text27.Text * n), q)Text34.Text = IIf(Val(Text21.Text * n) = q, Val(Text21.Text * n), q)Text44.Text = IIf(Val(Text21.Text * n) = q, Val(Text21.Text * n), q)Text38.Text = IIf(Val(Text20.Text * n) = q, Val(Text20.Text * n), q)Text48.Text = IIf(Val(Text20.Text * n) = q, Val(Text20.Text * n), q) A-D流量分配f = IIf(T12369 T12569, T12369, T12569)f = IIf(f T12589, f, T12589)f = IIf(f T14569, f, T14569)f = IIf(f T14589, f, T14589)f = IIf(f T14789, f, T14789)Select Case f Case Is = T14569 temp = Val(Text33.Text) If q - Val(Text33.Text) Val(Text16.Text * n) Then Text33.Text = q Text29.Text = Val(Text29.Text) + temp + Val(Text16.Text * n) - q Text35.Text = Val(Text35.Text) + temp + Val(Text16.Text * n) - q Else Text33.Text = Val(Text33.Text) + Val(Text16.Text * n) Text35.Text = 0 End If Text39.Text = Val(Text33.Text) - temp Text41.Text = Val(Text33.Text) - temp Text48.Text = Val(Text48.Text) + Val(Text33.Text) - temp Text45.Text = Val(Text45.Text) + Val(Text35.Text) Text52.Text = Val(Text52.Text) + Val(Text35.Text)End Select B-C流量分配T36987 = Val(Text5.Text) + Val(Text10.Text) + Val(Text12.Text) + Val(Text11.Text)T32147 = Val(Text2.Text) + Val(Text1.Text) + Val(Text3.Text) + Val(Text8.Text)T32547 = Val(Text2.Text) + Val(Text4.Text) + Val(Text6.Text) + Val(Text8.Text)T32587 = Val(Text2.Text) + Val(Text4.Text) + Val(Text9.Text) + Val(Text11.Text)T36547 = Val(Text5.Text) + Val(Text7.Text) + Val(Text6.Text) + Val(Text8.Text)T36587 = Val(Text5.Text) + Val(Text7.Text) + Val(Text9.Text) + Val(Text11.Text)f = IIf(T36987 T32147, T36987, T32147)f = IIf(f T32547, f, T32547)f = IIf(f T32587, f, T32587)f = IIf(f T36547, f, T36547)f = IIf(f T36587, f, T36587)Select Case f Case T36547 temp1 = Val(Text38.Text) temp2 = Val(Text43.Text) temp9 = IIf(q - temp1 q - temp2, q - temp1, q - temp2) If temp9 Val(Text19.Text * n) Then Text38.Text = Val(Text38.Text) + temp9 Text32.Text = Val(Text32.Text) + Val(Text19.Text * n) - temp9 Text35.Text = Val(Text35.Text) + Val(Text19.Text * n) - temp9 Text42.Text = Val(Text42.Text) + temp9 Else Text38.Text = Val(Text38.Text) + Val(Text19.Text * n) End If If q - temp2 Val(Text19.Text * n) Then Text40.Text = q - temp2 Text43.Text = q Text45.Text = Val(Text45.Text) + Val(Text19.Text * n) - q + temp2 Text49.Text = Val(Text49.Text) + Val(Text19.Text * n) - q + temp2 Else Text40.Text = Val(Text40.Text) + Text19.Text * n Text43.Text = Val(Text19.Text * n) + Val(Text43.Text) End IfEnd Select A-D和D-A;B-C和C-B寻找最短路径方法一样并且路径也相同 C-B temp3 = q1 - Val(Text39.Text) temp4 = q - Val(Text44.Text) If temp3 temp4 Then temp5 = temp3ElseIf temp3 = temp4 Then temp5 = temp3Else temp5 = temp4End IfIf temp5 Val(Text22.Text * n) Then Text44.Text = Val(Text44.Text) + temp5 Text39.Text = Val(Text39.Text) + temp5 Text50.Text = Val(Text50.Text) + Val(Text22.Text * n) - temp5 Text46.Text = Val(Text46.Text) + Val(Text22.Text * n) - temp5Else Text44.Text = Val(Text44.Text) + Val(Text22.Text * n) Text39.Text = Val(Text39.Text) + Val(Text22.Text * n)End If temp6 = q1 - Val(Text41.Text) temp7 = q - Val(Text37.Text)If temp6 temp7 Then temp8 = temp6ElseIf temp6 = temp7 Then temp8 = temp6Else temp8 = temp7End IfIf temp8 Val(Text22.Text * n) Then Text41.Text = Val(Text41.Text) + temp8 Text37.Text = Val(Text37.Text) + temp8 Text36.Text = Val(Text36.Text) + Val(Text22.Text * n) - temp8 Text31.Text = Val(Text31.Text) + Val(Text22.Text * n) - temp8Else Text41.Text = Val(Text41.Text) + Val(Text22.Text * n) Text37.Text = Val(Text37.Text) + Val(Text22.Text * n)End If D-A temp10 = q1 - Val(Text42.Text) temp11 = q - Val(Text47.Text) temp12 = IIf(temp10 temp11, temp10, temp11)If temp12 Val(Text25.Text * n) Then Text47.Text = Val(Text47.Text) + temp12 Text42.Text = Val(Text42.Text) + temp12 Text51.Text = Val(Text51.Text) + Val(Text25.Text * n) - temp12 Text46.Text = Val(Text46.Text) + Val(Text25.Text * n) - temp12Else Text47.Text = Val(Text47.Text) + Val(Text25.Text * n) Text42.Text = Val(Text42.Text) + Val(Text25.Text * n)End If Text40.Text = Text39.Text Text34.Text = Text33.Text Text36.Text = Text35.Text Text30.Text = Text29.Text Text46.Text = Text45.Text Text42.Text = Text41.TextEnd SubPrivate Sub Command11_Click()v = 50c = 750e = 750Text1.Text = t0(1) * (1 + 0.15 * (Val(Text29.Text) / c) 4)Text2.Text = t0(2) * (1 + 0.15 * (Val(Text31.Text) / c) 4)Text3.Text = t0(3) * (1 + 0.15 * (Val(Text33.Text) / c) 4)Text4.Text = t0(4) * (1 + 0.15 * (Val(Text35.Text) / e) 4)Text5.Text = t0(5) * (1 + 0.15 * (Val(Text37.Text) / c) 4)Text6.Text = t0(6) * (1 + 0.15 * (Val(Text39.Text) / e) 4)Text7.Text = t0(7) * (1 + 0.15 * (Val(Text41.Text) / e) 4)Text8.Text = t0(8) * (1 + 0.15 * (Val(Text43.Text) / c) 4)Text9.Text = t0(9) * (1 + 0.15 * (Val(Text45.Text) / e) 4)Text10.Text = t0(10) * (1 + 0.15 * (Val(Text47.Text) / c) 4)Text11.Text = t0(11) * (1 + 0.15 * (Val(Text49.Text) / c) 4)Text12.Text = t0(12) * (1 + 0.15 * (Val(Text51.Text) / c) 4)End SubPrivate Sub Command12_Click()Text29.Text = : Text30.Text = : Text31.Text = : Text32.Text = Text33.Text = : Text34.Text = : Text35.Text = : Text36.Text = Text37.Text = : Text38.Text = : Text39.Text = : Text40.Text = Text41.Text = : Text42.Text = : Text43.Text = : Text44.Text = Text45.Text = : Text46.Text = : Text47.Text = : Text48.Text = Text49.Text = : Text50.Text = : Text51.Text = : Text52.Text = Dim X, Y, z, s, m, n, o As Integerz = InputBox(请输入第四次OD量分配率, , 0.1)n = 0.65T123 = Val(Text1.Text) + Val(Text2.Text)T14523 = Val(Text3.Text) + Val(Text6.Text) + Val(Text4.Text) + Val(Text2.Text)T14563 = Val(Text3.Text) + Val(Text6.Text) + Val(Text7.Text) + Val(Text5.Text)T12369 = Val(Text1.Text) + Val(Text2.Text) + Val(Text5.Text) + Val(Text10.Text)T12569 = Val(Text1.Text) + Val(Text4.Text) + Val(Text7.Text) + Val(Text10.Text)T12589 = Val(Text1.Text) + Val(Text4.Text) + Val(Text9.Text) + Val(Text12.Text)T14569 = Val(Text3.Text) + Val(Text6.Text) + Val(Text7.Text) + Val(Text10.Text)T14589 = Val(Text3.Text) + Val(Text6.Text) + Val(Text8.Text) + Val(Text12.Text)T14789 = Val(Text3.Text) + Val(Text6.Text) + Val(Text11.Text) + Val(Text12.Text)q = 750q1 = 750 A-BIf T123 T14523 And T123 T14563 Then If Val(Text14.Text * n) = q Then Text29.Text = Text14.Text * n Text31.Text = Text14.Text * n Else Text29.Text = q Text31.Text = q End IfEnd If A-CIf Val(Text3.Text) + Val(Text8.Text) Val(Text1.Text) + Val(Text6.Text) + Val(Text4.Text) + Val(Text8.Text) And Val(Text1.Text) + Val(Text2.Text) Val(Text1.Text) + Val(Text4.Text) + Val(Text9.Text) + Val(Text11.Text) Then If Val(Text15.Text * n) = q Then Text33.Text = Text15.Text * n Text43.Text = Text15.Text * n Else Text33.Text = q Text43.Text = q End IfEnd If 由于城市道路网四周网阻抗相同 C - A, B - A, B - D, D - B, c - D, D - c分配流量方法相同Text30.Text = IIf(Val(Text17.Text * n) = q, Val(Text17.Text * n), q)Text32.Text = IIf(Val(Text17.Text * n) = q, Val(Text17.Text * n), q)Text37.Text = IIf(Val(Text26.Text * n) = q, Val(Text26.Text * n), q)Text47.Text = IIf(Val(Text26.Text * n) = q, Val(Text26.Text * n), q)Text50.Text = IIf(Val(Text24.Text * n) = q, Val(Text24.Text * n), q)Text52.Text = IIf(Val(Text24.Text * n) = q, Val(Text24.Text * n), q)Text49.Text = IIf(Val(Text27.Text * n) = q, Val(Text27.Text * n), q)Text51.Text = IIf(Val(Text27.Text * n) = q, Val(Text27.Text * n), q)Text34.Text = IIf(Val(Text21.Text * n) = q, Val(Text21.Text * n), q)Text44.Text = IIf(Val(Text21.Text * n) = q, Val(Text21.Text * n), q)Text38.Text = IIf(Val(Text20.Text * n) = q, Val(Text20.Text * n), q)Text48.Text = IIf(Val(Text20.Text * n) = q, Val(Text20.Text * n), q) A-D流量分配f = IIf(T12369 T12569, T12369, T12569)f = IIf(f T12589, f, T12589)f = IIf(f T14569, f, T14569)f = IIf(f T14589, f, T14589)f = IIf(f T14789, f, T14789)Select Case f Case Is = T14569 temp = Val(Text33.Text) If q - Val(Text33.Text) Val(Text16.Text * n) Then Text33.Text = q Text29.Text = Val(Text29.Text) + temp + Val(Text16.Text * n) - q Text35.Text = Val(Text35.Text) + temp + Val(Text16.Text * n) - q Else Text33.Text = Val(Text33.Text) + Val(Text16.Text * n) Text35.Text = 0 End If Text39.Text = Val(Text33.Text) - temp Text41.Text = Val(Text33.Text) - temp Text48.Text = Val(Text48.Text) + Val(Text33.Text) - temp Text45.Text = Val(Text45.Text) + Val(Text35.Text) Text52.Text = Val(Text52.Text) + Val(Text35.Text)End Select B-C流量分配T36987 = Val(Text5.Text) + Val(Text10.Text) + Val(Text12.Text) + Val(Text11.Text)T32147 = Val(Text2.Text) + Val(Text1.Text) + Val(Text3.Text) + Val(Text8.Text)T32547 = Val(Text2.Text) + Val(Text4.Text) + Val(Text6.Text) + Val(Text8.Text)T32587 = Val(Text2.Text) + Val(Text4.Text) + Val(Text9.Text) + Val(Text11.Text)T36547 = Val(Text5.Text) + Val(Text7.Text) + Val(Text6.Text) + Val(Text8.Text)T36587 = Val(Text5.Text) + Val(Text7.Text) + Val(Text9.Text) + Val(Text11.Text)f = IIf(T36987 T32147, T36987, T32147)f = IIf(f T32547, f, T32547)f = IIf(f T32587, f, T32587)f = IIf(f T36547, f, T36547)f = IIf(f T36587, f, T36587)Select Case f Case T36547 temp1 = Val(Text38.Text) temp2 = Val(Text43.Text) temp9 = IIf(q - temp1 q - temp2, q - temp1, q - temp2) If temp9 Val(Text19.Text * n) Then Text38.Text = Val(Text38.Text) + temp9 Text32.Text = Val(Text32.Text) + Val(Text19.Text * n) - temp9 Text35.Text = Val(Text35.Text) + Val(Text19.Text * n) - temp9 Text42.Text = Val(Text42.Text) + temp9 Else Text38.Text = Val(Text38.Text) + Val(Text19.Text * n) End If If q - temp2 Val(Text19.Text * n) Then Text40.Text = q - temp2 Text43.Text = q Text45.Text = Val(Text45.Text) + Val(Text19.Text * n) - q + temp2 Text49.Text = Val(Text49.Text) + Val(Text19.Text * n) - q + temp2 Else Text40.Text = Val(Text40.Text) + Text19.Text * n Text43.Text = Val(Text19.Text * n) + Val(Text43.Text) End IfEnd Select A-D和D-A;B-C和C-B寻找最短路径方法一样并且路径也相同 C-B temp3 = q1 - Val(Text39.Text) temp4 = q - Val(Text44.Text) If temp3 temp4 Then temp5 = temp3ElseIf temp3 = temp4 Then temp5 = temp3Else temp5 = temp4End IfIf temp5 Val(Text22.Text * n) Then Text44.Text = Val(Text44.Text) + temp5 Text39.Text = Val(Text39.Text) + temp5 Text50.Text = Val(Text50.Text) + Val(Text22.Text * n) - temp5 Text46.Text = Val(Text46.Text) + Val(Text22.Text * n) - temp5Else Text44.Text = Val(Text44.Text) + Val(Text22.Text * n) Text39.Text = Val(Text39.Text) + Val(Text22.Text * n)End If temp6 = q1 - Val(Text41.Text) temp7 = q - Val(Text37.Text)If temp6 temp7 Then temp8 = temp6ElseIf temp6 = temp7 Then temp8 = temp6Else temp8 = temp7End IfIf temp8 Val(Text22.Text * n) Then Text41.Text = Val(Text41.Text) + temp8 Text37.Text = Val(Text37.Text) + temp8 Text36.Text = Val(Text36.Text) + Val(Text22.Text * n) - temp8 Text31.Text = Val(Text31.Text) + Val(Text22.Text * n) - temp8Else Text41.Text = Val(Text41.Text) + Val(Text22.Text * n) Text37.Text = Val(Text37.Text) + Val(Text22.Text * n)End If D-A temp10 = q1 - Val(Text42.Text) temp11 = q - Val(Text47.Text) temp12 = IIf(temp10 temp11, temp10, temp11)If temp12 Val(Text25.Text * n) Then Text47.Text = Val(Text47.Text) + temp12 Text42.Text = Val(Text42.Text) + temp12 Text51.Text = Val(Text51.Text) + Val(Text25.Text * n) - temp12 Text46.Text = Val(Text46.Text) + Val(Text25.Text * n) - temp12Else Text47.Text = Val(Text47.Text) + Val(Text25.Text * n) Text42.Text = Val(Text42.Text) + Val(Text25.Text * n)End If Text40.Text = Text39.Text Text34.Text = Text33.Text Text36.Text = Text35.Text Text30.Text = Text29.Text Text46.Text = Text45.Text Text42.Text = Text41.Text End SubPrivate Sub Command13_Click()v = 50c = 750e = 750Text1.Text = t0(1) * (1 + 0.15 * (Val(Text29.Text) / c) 4)Text2.Text = t0(2) * (1 + 0.15 * (Val(Text31.Text) / c) 4)Text3.Text = t0(3) * (1 + 0.15 * (Val(Text33.Text) / c) 4)Text4.Text = t0(4) * (1 + 0.15 * (Val(Text35.Text) / e) 4)Text5.Text = t0(5) * (1 + 0.15 * (Val(Text37.Text) / c) 4)Text6.Text = t0(6) * (1 + 0.15 * (Val(Text39.Text) / e) 4)Text7.Text = t0(7) * (1 + 0.15 * (Val(Text41.Text) / e) 4)Text8.Text = t0(8) * (1 + 0.15 * (Val(Text43.Text) / c) 4)Text9.Text = t0(9) * (1 + 0.15 * (Val(Text45.Text) / e) 4)Text10.Text = t0(10) * (1 + 0.15 * (Val(Text47.Text) / c) 4)Text11.Text = t0(11) * (1 + 0.15 * (Val(Text49.Text) / c) 4)Text12.Text = t0(12) * (1 + 0.15 * (Val(Text51.Text) / c) 4)End SubPrivate Sub Command14_Click()Text29.Text = : Text30.Text = : Text31.Text = : Text32.Text = Text33.Text = : Text34.Text = : Text35.Text = : Text36.Text = Text37.Text = : Text38.Text = : Text39.Text = : Text40.Text = Text41.Text = : Text42.Text = : Text43.Text = : Text44.Text = Text45.Text = : Text46.Text = : Text47.Text = : Text48.Text = Text49.Text = : Text50.Text = : Text51.Text = : Text52.Text = Dim X, Y, z, s, m, n, o As IntegerY = InputBox(请输入第五次OD量分配率, , 0.1)m = 0.75T123 = Val(Text1.Text) + Val(Text2.Text)T14523 = Val(Text3.Text) + Val(Text6.Text) + Val(Text4.Text) + Val(Text2.Text)T14563 = Val(Text3.Text) + Val(Text6.Text) + Val(Text7.Text) + Val(Text5.Text)T12369 = Val(Text1.Text) + Val(Text2.Text