- 微信
- 微博
  
  分享文章到微博
- 复制链接
  
  复制链接到剪贴板

【数学建模】基于matlab模拟十字路口车辆通行【含Matlab源码 1106期】

海神之光发表于 2022/05/29 02:13:04 2022/05/29

【摘要】一、获取代码方式获取代码方式1：完整代码已上传我的资源：【数学建模】基于matlab模拟十字路口车辆通行【含Matlab源码 1106期】获取代码方式2：通过订阅紫极神光博客付费专栏，凭支付凭...

一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【数学建模】基于matlab模拟十字路口车辆通行【含Matlab源码 1106期】

获取代码方式2：
通过订阅紫极神光博客付费专栏，凭支付凭证，私信博主，可获得此代码。

备注：
订阅紫极神光博客付费专栏，可免费获得1份代码（有效期为订阅日起，三天内有效）；

二、部分源代码

clear;clc;
width=30    %道路宽
rl=5*width/8;     %左转虚拟车道
rr=1*width/8;     %右转虚拟车道
drr=pi/2*rr;
dstr=width-sqrt(rl^2-(width/8)^2);   %直行距离
dl=(acos(width/8/rl))*rl;      %左拐距离
ck=1.8                           ;     %车宽
cl=3.8                           ;     %车长
ddstr=sqrt((rl+ck/2)^2-(width/8+ck/2)^2)-sqrt(rl^2-(width/8)^2);
ddl=(acos(width/8/rl)-acos((width/8+ck/2)/(rl+ck/2)))*rl;
sl=dl;
sstr=dstr;
onetime=160;
dmin=2;

azb=[];
bzb=[];

alcar=[];                           %行车矩阵
arcar=[];
astrcar=[];
blcar=[];
brcar=[];
bstrcar=[];
cjjlal=[];
cjjlar=[];
cjjlastr=[];%车间距
cjjlbl=[];
cjjlbr=[];
cjjlbstr=[];%车间距
red=30;
yel=10;
val0=4      ;                        %a左行车速度（此处设未冲突速度）
var0=5       ;                         %a右行车速度（此处设未冲突速度）
vastr0=6     ;                         %a直行车速度（此处设未冲突速度）
vbl0=val0       ;                        %b左行车速度（此处设未冲突速度）
vbr0=var0      ;                         %b右行车速度（此处设未冲突速度）
vbstr0= vastr0    ;                         %b直行车速度（此处设未冲突速度）

val=[];
val(1)=val0    ;
var=[];
var(1)=var0  ;
vastr=[]  ;
vastr(1)=vastr0;
vbl=[]  ;
vbl(1)=vbl0;
vbr=[];
vbr(1)=vbr0    ;
vbstr=[]
vbstr(1)=vbstr0  ;
dt=0.1       ;                       %时间步长
alcar(1)=val*dt;
arcar(1)=vbr(1)*dt;
astrcar(1)=vastr(1)*dt;

blcar(1)=vbl(1)*dt;
brcar(1)=vbr(1)*dt;
bstrcar(1)=vbstr(1)*dt;

T=dt;
sumal=1;
sumar=1;
sumastr=1;
sumbl=1;
sumbr=1;
sumbstr=1;

kal=1;
kar=1;
kastr=1;
kbl=1;
kbr=1;
kbstr=1;

dertal=0;
dertas=0;

%外部条件参数===========================================================================================================
yctime=0.5;%加速延迟时间
jsxs=0.7;%减速系数
ajia=8;%最大加速度
ajian=8;%最大减速度
azxgl=1;
ayzxgl=1;
bzxgl=1;
byzxgl=1;
pro=0.02;
sigma=1;
n=20;%画图中的模拟车！！！！！！！！！！！！！！！！！！！！！
ttt=0.000000001;

while T<=onetime
    T=T+dt;

    lor=mod(floor(T/(red+yel)),2);
    goy=mod(T,red+yel);
    %虚拟交警================================================================================================================
    if kal~=0
        cjjlal(1)=0;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%555
    end
    if kar~=0
        cjjlar(1)=0;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%555
    end
    if kastr~=0
        cjjlastr(1)=0;%%%%%%%%%%%%%%%%%%%%%%5
    end
    if kbl~=0
        cjjlbl(1)=0;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%555
    end
    if kbr~=0
        cjjlbr(1)=0;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%555
    end
    if kbstr~=0
        cjjlbstr(1)=0;%%%%%%%%%%%%%%%%%%%%%%5

    end
    if kbl~=0&&kastr~=0
        pbl=dl-blcar(1)-sl;
        pastr=dstr-astrcar(1)-sstr;
        if pbl<=0&&pastr<=0
            ttsl=(dl-blcar(1)-ddl)/vbl(1)+dertal;
            ttsstr=(dstr-astrcar(1)-ddstr)/vastr(1)+dertas;
            if(ttsstr<=ttsl)
                if  vastr(1)<vastr0
                    vastr(1)=min(vastr(1)+ajia*dt,vastr0);
                end
                xvbl=[(dl-blcar(1)-ddl)/((dstr-astrcar(1)+ddstr+cl)/vastr(1)),vbl0];

                vbl(1)=min(xvbl);
                jishi1=dt;
            elseif(ttsstr>ttsl)
                if vbl(1)<vbl0
                    vbl(1)=min(vbl(1)+ajia*dt,vbl0);
                end
                xvastr=[(dstr-astrcar(1)-ddstr)/((dl-blcar(1)+ddl+cl)/vbl(1)),vastr0];

                vastr(1)=min(xvastr);
                jishi1=dt;
            end
        end
    end

    if kal~=0
        cjjlal(1)=0;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%555
    end
    if kbstr~=0
        cjjlbstr(1)=0;%%%%%%%%%%%%%%%%%%%%%%5
    end
    if kal~=0&&kbstr~=0
        pal=dl-alcar(1)-sl;
        pbstr=dstr-bstrcar(1)-sstr;
        if pal<=0&&pbstr<=0
            ttsl=(dl-alcar(1)-ddl)/val(1)+dertal;
            ttsstr=(dstr-bstrcar(1)-ddstr)/vbstr(1)+dertas;
            if(ttsstr<=ttsl)
                if vbstr(1)<vbstr0
                    vbstr(1)=min(vbstr(1)+ajia*dt,vbstr0);
                end
                xval=[(dl-alcar(1)-ddl)/((dstr-bstrcar(1)+ddstr+cl)/vbstr(1)),val0];

                val(1)=min(xval);
                jishi1=dt;
            elseif(ttsstr>ttsl)
                if val(1)<val0
                    val(1)=min(val(1)+ajia*dt,val0);
                end
                xvbstr=[(dstr-bstrcar(1)-ddstr)/((dl-alcar(1)+ddl+cl)/val(1)),vbstr0];

                vbstr(1)=min(xvbstr);
                jishi1=dt;
            end
        end
    end

  
 
  1
  2
  3
  4
  5
  6
  7
  8
  9
  10
  11
  12
  13
  14
  15
  16
  17
  18
  19
  20
  21
  22
  23
  24
  25
  26
  27
  28
  29
  30
  31
  32
  33
  34
  35
  36
  37
  38
  39
  40
  41
  42
  43
  44
  45
  46
  47
  48
  49
  50
  51
  52
  53
  54
  55
  56
  57
  58
  59
  60
  61
  62
  63
  64
  65
  66
  67
  68
  69
  70
  71
  72
  73
  74
  75
  76
  77
  78
  79
  80
  81
  82
  83
  84
  85
  86
  87
  88
  89
  90
  91
  92
  93
  94
  95
  96
  97
  98
  99
  100
  101
  102
  103
  104
  105
  106
  107
  108
  109
  110
  111
  112
  113
  114
  115
  116
  117
  118
  119
  120
  121
  122
  123
  124
  125
  126
  127
  128
  129
  130
  131
  132
  133
  134
  135
  136
  137
  138
  139
  140
  141
  142
  143
  144
  145
  146
  147
  148
  149
  150
  151
  152
  153
  154
  155
  156
  157
  158
  159
  160
  161
  162
  163
  164
  165
  166
  167
  168
  169
  170
  171
  172
  173
  174
  175
  176
  177
  178
  179
  180
  181
  182
  183
  184
  185
  186
  187
  188
  189
  190
  191
  192
  193
  194
  195