群控电梯matlab仿真代码

以下是一个简单的群控电梯MATLAB仿真代码示例：

clear all;
close all;

% 定义电梯参数
num_floors = 10; % 楼层数
num_elevators = 3; % 电梯数量
capacity = 15; % 电梯容量
elevator_speed = 1; % 电梯速度

% 定义人流量
arrival_rate = 1/3; % 到达速率
departure_rate = 1/5; % 离开速率

% 初始化电梯位置和状态
elevator_positions = zeros(1, num_elevators);
elevator_directions = ones(1, num_elevators);

% 初始化楼层状态
floors = struct('up_button', zeros(1, num_floors), 'down_button', zeros(1, num_floors));

% 初始化人员信息
people = struct('start_floor', [], 'end_floor', [], 'elevator', []);

% 初始化仿真时间和时间步长
sim_time = 1000;
dt = 0.1;

% 开始模拟
for t = 0:dt:sim_time
    
    % 处理人员运动和电梯状态
    for i = 1:length(people)
        if isempty(people(i).elevator) % 如果人不在电梯里
            if elevator_positions(1) == people(i).start_floor && elevator_directions(1) == 1 % 如果电梯在起点且向上运行
                if sum([people.elevator] == 1) < capacity % 如果电梯还有容量
                    people(i).elevator = 1; % 加入电梯
                end
            elseif elevator_positions(1) == people(i).start_floor && elevator_directions(1) == -1 % 如果电梯在起点且向下运行
                if sum([people.elevator] == 1) < capacity % 如果电梯还有容量
                    people(i).elevator = 1; % 加入电梯
                end
            elseif floors(people(i).start_floor).up_button == 1 && elevator_directions(1) == 1 % 如果电梯向上运动并且有人在该层按上行按钮
                if sum([people.elevator] == 1) < capacity % 如果电梯还有容量
                    people(i).elevator = 1; % 加入电梯
                    floors(people(i).start_floor).up_button = 0; % 关闭该层上行按钮
                end
            elseif floors(people(i).start_floor).down_button == 1 && elevator_directions(1) == -1 % 如果电梯向下运动并且有人在该层按下行按钮
                if sum([people.elevator] == 1) < capacity % 如果电梯还有容量
                    people(i).elevator = 1; % 加入电梯
                    floors(people(i).start_floor).down_button = 0; % 关闭该层下行按钮
                end
            end
        else % 如果人在电梯里
            if elevator_positions(people(i).elevator) == people(i).end_floor % 如果到达目的地
                people(i).elevator = []; % 离开电梯
            end
        end
    end
    
    % 处理电梯位置和方向
    for i = 1:num_elevators
        if elevator_positions(i) == num_floors % 如果电梯在顶层
            elevator_directions(i) = -1; % 改变方向
        elseif elevator_positions(i) == 1 % 如果电梯在底层
            elevator_directions(i) = 1; % 改变方向
        end
        
        if elevator_directions(i) == 1 % 向上运动
            elevator_positions(i) = elevator_positions(i) + elevator_speed*dt;
        else % 向下运动
            elevator_positions(i) = elevator_positions(i) - elevator_speed*dt;
        end
    end
    
    % 处理人员到达和离开
    if rand() < arrival_rate*dt % 如果有人到达
        start_floor = randi(num_floors); % 随机选择起点楼层
        end_floor = randi(num_floors); % 随机选择目的地楼层
        while end_floor == start_floor % 如果目的地与起点一样，则重新随机选择
            end_floor = randi(num_floors);
        end
        people(end+1) = struct('start_floor', start_floor, 'end_floor', end_floor, 'elevator', []); % 添加新的人员
        if end_floor > start_floor % 如果目的地在起点上方
            floors(start_floor).up_button = 1; % 按上行按钮
        else % 如果目的地在起点下方
            floors(start_floor).down_button = 1; % 按下行按钮
        end
    end
    
    if rand() < departure_rate*dt && ~isempty(people) % 如果有人离开
        idx = randi(length(people)); % 随机选择一个人
        people(idx) = []; % 删除该人员
    end
    
    % 可视化
    clf;
    hold on;
    plot(0, 0, 'ko', 'MarkerFaceColor', 'k');
    plot(num_floors+1, 0, 'ko', 'MarkerFaceColor', 'k');
    for i = 1:num_floors
        text(-0.5, i, num2str(num_floors-i+1));
        line([0,num_floors+1], [i,i]);
        if floors(i).up_button == 1
            text(0.5, i, '↑');
        end
        if floors(i).down_button == 1
            text(0.5, i, '↓');
        end
    end
    for i = 1:num_elevators
        rectangle('Position',[i-0.5,elevator_positions(i)-0.5,1,1],'FaceColor','b');
    end
    for i = 1:length(people)
        if isempty(people(i).elevator)
            rectangle('Position',[0.5,people(i).start_floor-0.5,0.5,0.5],'FaceColor','r');
        else
            rectangle('Position',[people(i).elevator-0.5,elevator_positions(people(i).elevator)-0.5,0.5,0.5],'FaceColor','r');
        end
    end
    axis([-1,num_elevators+1,-1,num_floors+1]);
    axis equal;
    drawnow;
    
end

此代码模拟了一个有三台电梯和十层楼的建筑，人员到达和离开速率是固定的，并且在每层楼都有上行和下行按钮。电梯按照最简单的规则运行：到达顶层或底层时改变方向。当电梯到达某层楼时，如果有人在该层按下行或上行按钮，则电梯会停下并让这些人进入。如果电梯还有容量，则在电梯起点处进入电梯的人也会被让进去。电梯内的人直到到达目的地才离开电梯。