交通运输学Matlab画图
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
用matlab画交通流密速关系
为了作图更明显,这里分别假设畅行速度Sf的数值为1000,阻塞密度Dj的数值为500,这可能和实际情况不符,但是具体运用是可以对程序中的Sf 值和Dj值进行更改。
一、密度、流速和流量关系图及各种假设情况下的图。
图中:
左上:V=DS;
右上:1931年Green Shields 假设S=Sf(1-D/Dj);
左下:1959年Greenberg 假设S=Sf*ln(Dj/D);
右下:1961年Underwood 假设。
Matlab的程序如下:
Sf=1000;
Dj=500;
subplot(2,2,1);
for D=0:100:Dj;
for S=0:10:Sf;
V=D*S;
plot3(S,D,V,'+r');
hold on;
end;
end;
for S=0:100:Sf;
for D=0:10:Dj;
V=D*S;
plot3(S,D,V,'+r');
hold on;
end;
end;
Zlabel('流量','FontWeight','bold')
Xlabel('流速','FontWeight','bold')
Ylabel('密度','FontWeight','bold')
title('\fontname{隶书}空间图形')
subplot(2,2,2)
D=0:10:Dj;
S=Sf*(1-D/Dj);
plot(S,D,'b*');
Zlabel('流量','FontWeight','bold')
Xlabel('流速','FontWeight','bold')
Ylabel('密度','FontWeight','bold')
title('\fontname{隶书} 1934年投影到XY平面') view(-37.5,80)
subplot(2,2,3)
Sf=1000;
Dj=500;
D=0:10:Dj;
S=Sf*log(Dj./D);
plot(S,D,'kh');
Zlabel('流量','FontWeight','bold')
Xlabel('流速','FontWeight','bold')
Ylabel('密度','FontWeight','bold')
title('\fontname{隶书} 1959年投影到XY平面')
view(-37.5,80)
subplot(2,2,4)
Sf=1000;
Dj=500;
D=0:10:Dj;
S=Sf*exp(-D/Dj);
plot(S,D,'gp');
Zlabel('流量','FontWeight','bold')
Xlabel('流速','FontWeight','bold')
Ylabel('密度','FontWeight','bold')
title('\fontname{隶书} 1961年投影到XY平面')
view(-37.5,80)
二、Green Shields假设涉及到的图形。
图中:
左上:Green Shields假设图形;
右上:Green Shields假设图形在XOY平面上的投影;左下:Green Shields假设图形在YOZ平面上的投影;右下:Green Shields假设图形在ZOX平面上的投影。
Matlab的程序如下:
Sf=1000;
Dj=500;
subplot(2,2,1)
D=0:10:Dj;
S=Sf*(1-D/Dj);
V=S.*D;
plot3(S,D,V,'rh');
Zlabel('流量','FontWeight','bold')
Xlabel('流速','FontWeight','bold')
Ylabel('密度','FontWeight','bold')
title('\fontname{隶书} Green Shields假设')
subplot(2,2,2)
D=0:10:Dj;
S=Sf*(1-D/Dj);
plot(S,D,'g*');
Zlabel('流量','FontWeight','bold')
Xlabel('流速','FontWeight','bold')
Ylabel('密度','FontWeight','bold')
title('\fontname{隶书} Green Shields假设XOY面投影')
view(-37.5,80)
subplot(2,2,3)
D=0:10:Dj;
S=Sf*(1-D/Dj);
V=S.*D;
plot3(S,D,V,'kd');
Zlabel('流量','FontWeight','bold')
Xlabel('流速','FontWeight','bold')
Ylabel('密度','FontWeight','bold')
title('\fontname{隶书} Green Shields假设YOZ面投影')
view(90,0)
subplot(2,2,4)
D=0:10:Dj;
S=Sf*(1-D/Dj);
V=S.*D;
plot3(S,D,V,'bs');
Zlabel('流量','FontWeight','bold')
Xlabel('流速','FontWeight','bold')
Ylabel('密度','FontWeight','bold')
title('\fontname{隶书} Green Shields假设ZOX面投影') view(360,0)。