function?[?v?d?p?]?=?multi_driveway（?nlncfpdtnt?）
%??在某一特定車流密度下的（車流密度由fp決定）單、雙車道仿真模型
%??nc:車道數目（1或2），nl:車道長度——輸入參數
%??v:平均速度，d:換道次數（1000次）p:車流密度——輸出參數
%??dt:仿真步長時間，nt:仿真步長數目——輸入參數
%??fp:車道入口處新進入車輛的概率——輸入參數
%??test:
%??nl?=?400;fp?=?0.5;
%??nc?=?2;dt=0.01;nt=500;
???%構造元胞矩陣
???B=ones（2*nc+1nl+2）;
???%奇數行為不可行車道
???B（1:2:（2*nc+1）:）=1.2;
???%初始化仿真元胞狀態（1為無車，0為有車）
???bb=B（2:2:2*nc:）;bb（bb~=0）=1;B（2:2:2*nc:）=bb;B（2:2:2*ncend）=0;
???%顯示初始交通流圖
???figure（1）;
???H=imshow（B[]）;
???set（gcf‘position‘[241?132?560?420]）?;%241?132?560?420
???set（gcf‘doublebuffer‘‘on‘）;??%241
???title（‘cellular-automation?to?traffic?modeling‘‘color‘‘b‘）;
???%初始化化存儲元胞上車輛狀態的矩陣
???S（1:ncnl）?=?0;
???Q（1:nc1:2）?=?0;
???Acc（1:nc1:（nl+2））=0;
???%初始化換道頻率、平均速度、車流密度相關變量
???ad?=?0;
???av（1:nt）?=?0;
???ap（1:nt）?=?0;
???c?=?1;
???for?n?=?1:nt
??????A=B（2:2:2*nc:）;
??????%確定前n-2個車輛的狀態
??????S（::）?=?0;
??????S（A（:1:end-2）==0&A（:2:end-1）==1&A（:3:end）==1）=2;%加速的車
??????S（A（:1:e