作者您好，你这个单向三车道的元胞自动机怎么写代码出它的密度流量图呢

Question

不知怎么出流量密度图
%% 单向3车道主函数
%% 注意：在本模型中车辆位置信息由车道和列两个维度唯一确定
function []= main(lane_length,car_rate,v_max,time_max,time_span,p_slowdown,p_changelane)

car_number = fix(1+(3*lane_length-1)*car_rate); %按车辆占有率算出的车辆数
%% 创建空间
space = zeros(3,lane_length);%元胞空间
car = struct('v',zeros(1,car_number),'m',zeros(1,car_number),'n',zeros(1,car_number));% 车辆信息结构体从左到右为速度，车道，列
%% 随机生成初始车辆信息
[space,car] = initialize(space,car,car_number,lane_length,v_max);

%% 显示初始仿真图
figure('name','单向3车道模型','position',[241 132 560 420],'doublebuffer','on');%窗口名称，位置,双缓存)
space = -1space;
H = imshow(space,[]);
title('单向3车道模型','color','red');
space = -1space;
%% 开始仿真
for time=1:time_max
%% 换道阶段
[car,space] = change_lane(space,car,car_number,p_changelane,v_max,lane_length);
%% 跟驰阶段
[space,car] = NaSch(space,car,v_max,lane_length,p_slowdown,car_number);
%% 显示仿真图
space = -1space;
set(H,'CData',space);
pause(time_span);
space = -1space;
end
end