MATLAB动态模拟地铁运行：从图论到动画的数学建模实践

1. 项目缘起：从一张静态图到动态旅程

最近在整理一些数学建模的旧项目时，翻到了一个挺有意思的案例：用MATLAB动态模拟武汉地铁2号线的运行。这个项目的初衷其实很简单，当时有学弟学妹在准备数学建模竞赛，题目涉及到城市交通网络的模拟与优化。他们手头有武汉地铁的线路数据，但不知道怎么把冷冰冰的站点坐标和连接关系，变成一个直观、能“动起来”的演示。跑来问我，有没有什么办法能快速做个可视化，最好还能模拟列车在不同区间的运行状态，比如速度变化、停站时间这些。

我一听，这不正是MATLAB的强项吗？用plot画静态线路图谁都会，但要让点在线段上移动，模拟出列车穿梭的效果，就需要用到动画和定时器的知识了。这个需求背后，其实是数学建模中一个非常核心的环节：将抽象的数学模型（如网络图、微分方程）转化为具象、可交互的可视化结果。一个好的动态模拟，不仅能验证模型逻辑是否正确，更能让评委或观众一眼看懂你的研究在解决什么问题，效果拔群。

所以，我花了一些时间，基于当时他们提供的粗略数据（主要是2号线各站点的经纬度坐标和连接顺序），写了一个基础的动态模拟框架。代码不算复杂，但涵盖了从数据预处理、图形绘制到动画生成的全流程。后来这个框架被多次复用，用于模拟不同假设下的客流、列车调度等场景，成了一个小“利器”。今天，我就把这个项目的核心思路、代码实现，以及其中几个容易踩坑的细节，完整地分享出来。无论你是想学习MATLAB动画制作，还是为数学建模寻找一个出彩的可视化方案，相信都能从中获得启发。

2. 核心思路拆解：动态模拟的“骨架”是什么

在做这个动态模拟之前，我们得先想清楚，所谓“地铁线路动态模拟”，到底要模拟哪些元素？它们之间的关系是什么？不能一上来就埋头写for循环。

我把整个系统分解为三个核心层：

2.1 数据层：静态网络的描述这是所有工作的基础。地铁线路本质上是一个图（Graph）。站点是节点（Node），轨道区间是边（Edge）。我们需要用数据把这个图描述清楚。

• 节点数据：每个站点的唯一标识（如站名或编号）、地理坐标（通常是经度、纬度，为了绘图方便，我们可能使用平面坐标x, y）。例如，汉口火车站可以是一个点 (x1, y1)，中山公园是另一个点 (x2, y2)。
• 边数据：描述哪些站点是相连的。通常用一个N×2的矩阵表示，每一行是一对相邻站点的编号。比如[1, 2; 2, 3; 3, 4]表示站点1-2相连，2-3相连，3-4相连，这样就构成了一条线。

对于武汉地铁2号线，我们首先需要获取这份拓扑数据。在实际项目中，这部分数据可能来自公开资料、GIS地图手动采集，或者赛题直接给出。有了它，我们才能画出那条静态的、蓝色的地铁线路图。

2.2 图形层：静态画面的绘制这一层的目标是把数据层描述的“图”，在MATLAB的图形窗口（Figure）里画出来。这一步相对直接：

1. 创建图形窗口：figure(‘Position‘, [100, 100, 800, 600])设置好大小。
2. 绘制站点（节点）：用scatter或plot命令，将每个站点的坐标画成一个小圆圈，并可以用text在旁边标注站名。
3. 绘制线路（边）：用plot命令，按照边数据连接对应的站点坐标，画出线段。为了美观，可以设置线宽、颜色（比如地铁2号线的标志色）。
4. 添加标题、坐标轴标签等修饰：title(‘武汉地铁2号线运行模拟‘)，xlabel(‘东-西方向‘)，ylabel(‘南-北方向‘)。

完成这一步，你就得到了一张漂亮的、静态的武汉地铁2号线地图。但这还不够，它不会动。

2.3 动态层：让列车“跑”起来这是本项目最核心、也最具技巧性的部分。如何让代表列车的标记点，沿着我们画好的线路平滑移动？ 关键在于参数化路径和定时更新。

• 参数化路径：我们画出的线路是由多条直线段（站点到站点）首尾相连组成的。假设一列车要从站点A移动到站点B。A坐标是(xa, ya)，B是(xb, yb)。列车在t时刻的位置(x(t), y(t))，可以用一个比例参数k（0 ≤ k ≤ 1）来表示：x(t) = xa + k * (xb - xa)y(t) = ya + k * (yb - ya)当k=0时，车在A站；k=1时，车到达B站。k从0线性增加到1，车就从A匀速移动到B。如果线路有多段，我们需要记录当前车在哪一段，以及在该段上的k值。
• 定时更新：有了位置计算公式，我们需要一个机制来不断更新k的值，并重绘列车的位置。MATLAB提供了两种主流方法：
  • 循环暂停法：在for循环中计算新位置，用plot或scatter重画列车点（记得用hold on和hold off管理图形），然后用pause(0.05)暂停一小段时间（如0.05秒），制造动画效果。这种方法简单直观，适合初学者和快速演示。
  • 定时器（Timer）回调法：创建一个timer对象，设定一个固定的时间间隔（如0.05秒）。每次定时器触发，就执行一个我们写好的“回调函数”，在这个函数里更新列车位置并刷新图形。这种方法更专业、更高效，动画更流畅，且不阻塞MATLAB命令行的其他操作，适合构建更复杂的交互式模拟。

在这个项目中，为了清晰易懂，我主要采用第一种方法。但会在后面详细解释第二种方法的优势和改进方案。

3. 从零开始：手把手构建基础动态模拟

理论说清楚了，我们直接上代码。我会分模块解释，并附上完整的、可运行的MATLAB源码。你可以直接复制到MATLAB里运行。

3.1 第一步：准备线路数据（模拟数据）由于我们无法直接获取精确的官方坐标，这里根据武汉地铁2号线的大致走向，构造一组模拟的站点坐标。在实际应用中，你需要替换成真实数据。

% 假设的武汉地铁2号线部分站点坐标 (x, y) % 这里x和y可以理解为处理后的平面坐标单位，并非真实经纬度 stations = { ‘天河机场‘, [0, 10]; ‘航空总部‘, [2, 9.5]; ‘宋家岗‘, [4, 9]; ‘巨龙大道‘, [6, 8.5]; ‘盘龙城‘, [8, 8]; ‘宏图大道‘, [10, 7.5]; ‘常青城‘, [12, 7]; ‘金银潭‘, [14, 6.5]; ‘常青花园‘, [16, 6]; ‘长港路‘, [18, 5.5]; ‘汉口火车站‘, [20, 5]; ‘范湖‘, [22, 4.5]; ‘王家墩东‘, [24, 4]; ‘青年路‘, [26, 3.5]; ‘中山公园‘, [28, 3]; ‘循礼门‘, [30, 2.5]; ‘江汉路‘, [32, 2]; ‘积玉桥‘, [34, 1.5]; ‘螃蟹岬‘, [36, 1]; ‘小龟山‘, [38, 0.5]; ‘洪山广场‘, [40, 0]; ‘中南路‘, [42, -0.5]; ‘宝通寺‘, [44, -1]; ‘街道口‘, [46, -1.5]; ‘广埠屯‘, [48, -2]; ‘虎泉‘, [50, -2.5]; ‘杨家湾‘, [52, -3]; ‘光谷广场‘, [54, -3.5] }; num_stations = size(stations, 1); station_names = stations(:, 1); % 站名单元格数组 station_coords = cell2mat(stations(:, 2)); % 坐标矩阵，N行2列

3.2 第二步：绘制静态线路图有了坐标，我们先把地铁线网画出来。

% 创建图形窗口 figure(‘Position‘, [100, 100, 900, 600]); hold on; % 保持图形，允许多次绘图叠加 grid on; % 显示网格，方便观察 box on; % 显示边框 axis equal; % 横纵坐标等比例，防止图形拉伸变形 % 1. 绘制线路（边）：用蓝色实线连接所有相邻站点 plot(station_coords(:,1), station_coords(:,2), ‘b-o‘, ... ‘LineWidth‘, 2, ... ‘MarkerSize‘, 8, ... ‘MarkerFaceColor‘, ‘w‘, ... ‘MarkerEdgeColor‘, ‘b‘); % ‘b-o‘ 表示蓝色(b)，线条和标记点用圆圈(o) % ‘MarkerFaceColor‘, ‘w‘ 将标记点内部填充为白色，更清晰 % 2. 高亮标注站点（节点） for i = 1:num_stations % 在站点坐标旁添加站名文本 text(station_coords(i,1)+0.3, station_coords(i,2), station_names{i}, ... ‘FontSize‘, 9, ‘Interpreter‘, ‘None‘); end % 3. 设置图形标题和标签 title(‘武汉地铁2号线线路图（静态）‘, ‘FontSize‘, 14, ‘FontWeight‘, ‘bold‘); xlabel(‘东-西方向 (模拟坐标)‘); ylabel(‘南-北方向 (模拟坐标)‘); xlim([min(station_coords(:,1))-2, max(station_coords(:,1))+2]); % 设置x轴范围，留点边距 ylim([min(station_coords(:,2))-2, max(station_coords(:,2))+2]); % 设置y轴范围 hold off;

运行这段代码，你应该能看到一条从“天河机场”蜿蜒至“光谷广场”的蓝色线路图，每个站点都有标注。这是我们的“舞台”。

3.3 第三步：实现单列车单向运行动画现在，让一列“车”从起点跑到终点。我们采用“循环暂停法”。

% 接续上面的代码，或者重新运行并hold on figure(‘Position‘, [100, 100, 900, 600]); hold on; grid on; box on; axis equal; % 绘制静态背景线路 plot(station_coords(:,1), station_coords(:,2), ‘b-o‘, ... ‘LineWidth‘, 2, ‘MarkerSize‘, 8, ‘MarkerFaceColor‘, ‘w‘, ‘MarkerEdgeColor‘, ‘b‘); for i = 1:num_stations text(station_coords(i,1)+0.3, station_coords(i,2), station_names{i}, ‘FontSize‘, 9); end title(‘武汉地铁2号线单列车运行动态模拟‘, ‘FontSize‘, 14, ‘FontWeight‘, ‘bold‘); xlabel(‘东-西方向 (模拟坐标)‘); ylabel(‘南-北方向 (模拟坐标)‘); xlim([min(station_coords(:,1))-2, max(station_coords(:,1))+2]); ylim([min(station_coords(:,2))-2, max(station_coords(:,2))+2]); % 动画参数设置 train_speed = 0.02; % 控制动画速度的参数，表示每帧k值的增量 % 这个值越小，车移动越慢；越大则越快。需要根据帧间隔调整。 % 初始化列车位置（从第一个站点开始） current_segment = 1; % 当前所在区间，即从站点1到站点2 k = 0; % 在当前区间上的位置比例，0表示在起点站 % 用红色五角星代表列车，先画在初始位置 train_plot = plot(station_coords(1,1), station_coords(1,2), ‘rp‘, ... ‘MarkerSize‘, 15, ‘MarkerFaceColor‘, ‘r‘); % ‘rp‘ 表示红色(r)，标记为五角星(p) % 动画主循环 while current_segment < num_stations % 当列车还未到达最后一个区间时 % 计算列车在当前区间上的新位置 start_idx = current_segment; end_idx = current_segment + 1; x_start = station_coords(start_idx, 1); y_start = station_coords(start_idx, 2); x_end = station_coords(end_idx, 1); y_end = station_coords(end_idx, 2); x_train = x_start + k * (x_end - x_start); y_train = y_start + k * (y_end - y_start); % 更新列车图形对象的位置 set(train_plot, ‘XData‘, x_train, ‘YData‘, y_train); % 增加比例k，模拟列车前进 k = k + train_speed; % 如果k超过1，说明已经到达当前区间的终点站 if k >= 1 % 列车“停靠”在终点站 set(train_plot, ‘XData‘, x_end, ‘YData‘, y_end); % 短暂停顿，模拟停靠时间 pause(0.5); % 切换到下一个区间，并重置k为0（从新区间的起点开始） current_segment = current_segment + 1; k = 0; % 如果已经到达最后一个站点，跳出循环 if current_segment >= num_stations break; end end % 刷新图形并暂停一小段时间，形成动画帧 drawnow; pause(0.05); % 控制动画帧率，约20帧/秒 end % 动画结束，列车到达终点站 text(station_coords(end,1), station_coords(end,2)-1, ‘列车已到达终点站: 光谷广场‘, ... ‘FontSize‘, 11, ‘Color‘, ‘r‘, ‘FontWeight‘, ‘bold‘, ‘HorizontalAlignment‘, ‘center‘); hold off;

运行这段代码，你会看到一个红色的五角星，从天河机场站出发，沿着蓝色线路，一站一站地运行，最终停在光谷广场站。在每一站（线段端点）会有短暂的停顿。

注意：这里的train_speed和pause(0.05)共同决定了动画的视觉速度。train_speed是逻辑上的速度（每帧移动的路径比例），pause是帧间隔。你可以调整它们来获得更慢或更快的动画效果。pause(0.5)模拟了列车进站停靠的耗时。

4. 进阶与优化：从“能动”到“好用”

上面的基础版本已经实现了核心的动态效果，但离一个“好用”的模拟工具还有距离。在实际的数学建模项目中，我们往往需要更复杂的功能。下面介绍几个关键的进阶优化点。

4.1 多列车与双向运行模拟现实中，地铁线上不可能只有一列车。我们需要模拟多列车间隔发车、对向行驶。这需要对数据结构和动画逻辑进行升级。

核心思路是：为每一列车创建一个独立的状态对象，记录其当前所在区间、位置比例k、运行方向（上行/下行）和速度。在动画循环中，遍历所有列车，更新它们各自的状态和位置。

% 假设我们定义两列车 % 列车1：从起点（站1）开往终点（站N）， 方向为1（正向） % 列车2：从终点（站N）开往起点（站1）， 方向为-1（反向），并且延迟若干帧发车 trains = struct(); trains(1).segment = 1; % 当前区间起点索引 trains(1).k = 0; % 在区间上的位置 trains(1).direction = 1; % 运行方向：1为从低索引站到高索引站 trains(1).speed = 0.015; trains(1).plot_handle = []; % 预留，用于存储图形对象句柄 trains(1).departure_delay = 0; % 发车延迟帧数 trains(2).segment = num_stations - 1; % 反向车起始区间是倒数第二个区间 trains(2).k = 0; trains(2).direction = -1; % 运行方向：-1为从高索引站到低索引站 trains(2).speed = 0.015; trains(2).plot_handle = []; trains(2).departure_delay = 100; % 延迟100帧后发车 % 在图形窗口中初始化两列车的图形对象，用不同颜色和形状区分 figure(‘Position‘, [100, 100, 900, 600]); hold on; grid on; box on; axis equal; % 绘制静态背景（同上，省略...） % ... colors = {‘r‘, ‘g‘}; % 列车1红色，列车2绿色 markers = {‘p‘, ‘s‘}; % 五角星和方块 for i = 1:length(trains) if trains(i).departure_delay <= 0 % 如果无需延迟，立即初始化位置 [x_init, y_init] = getTrainPosition(trains(i), station_coords); trains(i).plot_handle = plot(x_init, y_init, [colors{i}, markers{i}], ... ‘MarkerSize‘, 12, ‘MarkerFaceColor‘, colors{i}); end % 如果需要延迟，plot_handle暂时为空，后续在循环中创建 end % 动画主循环 for frame = 1:1000 % 假设运行1000帧 for i = 1:length(trains) % 处理发车延迟 if trains(i).departure_delay > 0 trains(i).departure_delay = trains(i).departure_delay - 1; if trains(i).departure_delay == 0 % 延迟结束，创建该列车的图形对象 [x_init, y_init] = getTrainPosition(trains(i), station_coords); trains(i).plot_handle = plot(x_init, y_init, [colors{i}, markers{i}], ... ‘MarkerSize‘, 12, ‘MarkerFaceColor‘, colors{i}); end continue; % 如果还在延迟，跳过本次更新 end % 如果列车图形对象已创建，则更新其状态和位置 if ~isempty(trains(i).plot_handle) % 更新位置比例k trains(i).k = trains(i).k + trains(i).speed; % 判断是否到达区间端点 if trains(i).k >= 1 % 到达正向端点 trains(i).k = 0; trains(i).segment = trains(i).segment + trains(i).direction; % 检查是否越界（到达线路终点或起点） if trains(i).segment < 1 || trains(i).segment >= num_stations % 列车到达终点，可以设置为消失、折返或停止 % 这里简单设置为停止更新 set(trains(i).plot_handle, ‘Visible‘, ‘off‘); trains(i).plot_handle = []; continue; end % 模拟停靠 pause(0.02); elseif trains(i).k < 0 % 对于反向行驶的车，k可能从1递减 % 处理反向行驶的逻辑（略，原理类似） end % 计算新的坐标并更新图形 [x_new, y_new] = getTrainPosition(trains(i), station_coords); set(trains(i).plot_handle, ‘XData‘, x_new, ‘YData‘, y_new); end end drawnow; pause(0.03); % 控制整体动画帧率 end hold off; % 辅助函数：根据列车状态计算其当前坐标 function [x, y] = getTrainPosition(train, station_coords) if train.direction == 1 start_idx = train.segment; end_idx = train.segment + 1; else % direction == -1 start_idx = train.segment + 1; % 注意反向时起点和终点互换 end_idx = train.segment; end x_start = station_coords(start_idx, 1); y_start = station_coords(start_idx, 2); x_end = station_coords(end_idx, 1); y_end = station_coords(end_idx, 2); % 对于反向车，k的意义是“距离终点的比例”，计算时需注意 % 这里简化处理，统一用线性插值，方向由起终点决定 x = x_start + train.k * (x_end - x_start); y = y_start + train.k * (y_end - y_start); end

这个版本引入了direction（方向）和departure_delay（发车延迟），使得模拟更加真实。你可以扩展这个结构，加入列车ID、当前载客量、运行状态（正常、延误）等属性，构建一个更复杂的仿真系统。

4.2 使用Timer对象实现更流畅的动画for循环加pause的方法简单，但有一个缺点：在循环执行期间，MATLAB的命令行窗口会被“阻塞”，你无法进行其他交互操作。对于需要长时间运行或与其他组件（如GUI按钮）联动的模拟，更好的选择是使用timer对象。

timer的工作原理是设定一个周期性的定时任务，MATLAB在后台执行，不干扰前台。修改上面的多列车动画主循环部分：

% 创建图形窗口和绘制静态背景（同上，省略...） % 初始化列车状态数组 trains（同上，省略...） % 创建并配置定时器 t = timer; t.Period = 0.03; % 定时器触发间隔，单位秒 t.ExecutionMode = ‘fixedRate‘; % 固定速率执行 t.TimerFcn = @(~, ~) updateAnimation(trains, station_coords); % 回调函数 % 注意：这里需要将trains和station_coords通过其他方式（如appdata或嵌套函数）传递给回调函数 % 下面用嵌套函数的方式演示 % 在图形窗口的UserData或应用数据中存储状态 hFig = gcf; appData.trains = trains; appData.station_coords = station_coords; guidata(hFig, appData); % 存储数据 % 定义嵌套的回调函数（实际写代码时，这部分应放在主函数内部） function updateAnimation(~, ~) % 从图形窗口获取数据 appData = guidata(hFig); trains = appData.trains; station_coords = appData.station_coords; for i = 1:length(trains) % 更新列车状态的逻辑（与之前for循环内的逻辑几乎相同） % ... (省略具体更新代码) % 计算新坐标 [x_new, y_new] = getTrainPosition(trains(i), station_coords); % 更新图形 if ~isempty(trains(i).plot_handle) && ishandle(trains(i).plot_handle) set(trains(i).plot_handle, ‘XData‘, x_new, ‘YData‘, y_new); end end % 将更新后的状态存回 appData.trains = trains; guidata(hFig, appData); drawnow; % 刷新图形 end % 启动定时器 start(t); % 此时，命令行可以继续输入命令，动画在后台运行。 % 需要停止时，使用 stop(t); delete(t);

使用timer后，动画的流畅度通常更好，且程序响应性更高。这对于构建交互式仿真界面至关重要。

4.3 集成真实数据与交互功能一个完整的模拟项目，最终需要和真实数据结合，并可能提供一些交互功能。

• 集成真实坐标：将模拟的station_coords替换为从百度地图/高德地图API获取的经纬度坐标，或者从公开的GIS数据中提取的平面坐标。注意，经纬度是球面坐标，直接绘图会导致变形。通常需要将其转换为某种平面投影坐标（如Web墨卡托投影），或者简单地进行线性缩放和平移，使其适应图形窗口。
• 添加交互控件：使用MATLAB的GUIDE或更现代的App Designer，可以创建图形用户界面（GUI）。在界面中添加滑块（uislider）控制模拟速度，按钮（uibutton）控制开始/暂停/重置，下拉菜单（uidropdown）选择显示哪条线路或哪种列车，甚至文本框（uieditfield）输入发车间隔等参数。这样，你的模拟就从一个脚本升级为一个可交互的工具。
• 数据输出与记录：在模拟过程中，除了图形化显示，还可以记录每列车在每个时间点的位置、速度、所在区间等信息到矩阵或表格中。模拟结束后，可以将这些数据保存为.mat或.csv文件，用于后续的分析，比如计算平均旅行时间、区间通过能力等，这正是数学建模定量分析的部分。

5. 实战踩坑与性能调优心得

在实现和多次使用这个模拟框架的过程中，我积累了一些宝贵的经验，这些在官方文档里通常不会细说。

5.1 图形刷新与性能瓶颈当站点数量很多（比如模拟整个城市地铁网络），或者同时运行的列车数量很大（比如几十列）时，动画可能会变得卡顿。这是因为在每一帧中，MATLAB都需要重绘大量的图形对象（train_plot）。

• 优化技巧1：减少不必要的绘图对象。对于固定的背景（如线路、站点、标签），只在初始化时绘制一次，之后不再改动。使用hold on后，它们会一直保留。动画循环中只更新移动对象（列车）的属性（XData,YData），这比每次都重新调用plot要高效得多。上面的代码已经采用了这种方法。
• 优化技巧2：使用drawnow limitrate或drawnow expose。drawnow命令强制MATLAB刷新图形。drawnow limitrate会限制刷新频率，防止程序过度消耗CPU，在动画场景下通常够用且更高效。drawnow expose只更新需要更新的图形区域，效率更高。可以尝试替换简单的drawnow。
• 优化技巧3：简化图形细节。如果还是卡，可以考虑暂时关闭网格（grid off）、将站点标记从圆圈(‘o‘)改为点(‘.‘)、减少文本标签的字体大小或数量。在调试阶段追求流畅，最终展示时再开启美化选项。

5.2 时间步长与速度的匹配train_speed（逻辑速度）和pause（帧间隔）共同决定了视觉速度。这里有个常见的坑：如果pause的时间太短，而train_speed太大，列车可能会“跳站”。因为在一帧内，k的增加可能直接超过了1，导致判断if k >= 1时，列车已经“飞”过了站点，停靠逻辑可能出错。

经验公式：确保train_speed * (1 / pause_time)不会过大。例如，pause(0.05)对应约20帧/秒，那么train_speed最好设置在0.01~0.05之间，这样列车至少需要20~100帧才能走完一个区间，移动看起来是平滑的。更稳健的做法是，在更新k后，增加一个边界检查：if k > 1, k = 1; end，确保不会超限。

5.3 处理折返与交路在实际地铁运营中，列车到达终点站后需要折返。在我们的简单模型里，当current_segment到达最后一个区间末端时，可以让direction乘以-1，并将current_segment设为倒数第二个区间，k设为0（或1，取决于方向定义），从而实现掉头。这需要仔细处理区间索引和k的增减逻辑，容易出错。建议单独写一个updateTrainStatus函数来封装这些状态转移逻辑，使主循环更清晰。

5.4 坐标变换与比例尺如果你用的是真实的经纬度（如[114.25, 30.62]），直接绘图会发现图形被压扁或集中在一个小点。因为经度和纬度的数值范围不同（武汉的经度差和纬度差大约为0.2和0.1），而且绝对值很大。我们需要进行归一化和缩放。

% 假设 lons 和 lats 是原始的经纬度数组 lons = [114.20, 114.21, 114.22, ...]; lats = [30.65, 30.64, 30.63, ...]; % 方法1：简单线性缩放至绘图窗口大小 x_coords = (lons - min(lons)) / (max(lons) - min(lons)) * 50; % 缩放到0~50范围 y_coords = (lats - min(lats)) / (max(lats) - min(lats)) * 30; % 缩放到0~30范围 station_coords = [x_coords‘, y_coords‘]; % 得到用于绘图的坐标 % 方法2：保持纵横比（Aspect Ratio） range_lon = max(lons) - min(lons); range_lat = max(lats) - min(lats); scale = 50 / max(range_lon, range_lat); % 以50为基准尺度 x_coords = (lons - min(lons)) * scale; y_coords = (lats - min(lats)) * scale; station_coords = [x_coords‘, y_coords‘]; % 这样处理后，图形的纵横比能反映实际地理形状的相对比例。

5.5 代码的模块化与复用最初的脚本可能把所有代码写在一个文件里。随着功能增加（多列车、交互、数据记录），代码会变得冗长难维护。好的做法是：

• 将数据加载和预处理写成独立的函数，如loadStationData.m。
• 将列车状态更新逻辑封装成函数，如updateTrainPosition.m。
• 将绘图初始化部分也写成函数，如plotMetroNetwork.m。
• 主脚本只负责调用这些函数和控制流程。这样不仅清晰，也方便你将这个“地铁模拟引擎”应用到其他城市或线路，只需更换数据加载函数即可。

这个基于MATLAB的武汉地铁2号线动态模拟项目，从最初一个简单的动画想法，逐步扩展为一个包含多列车调度、可集成真实数据、具备一定交互能力的轻量级仿真工具。它完美地展示了如何将数学建模中的图论概念，通过编程转化为直观的动态可视化，无论是用于竞赛演示、课程作业还是学术研究，都是一个非常实用的起点。希望这份详细的拆解和代码，能帮你打开思路，制作出属于你自己的、更精彩的地铁交通模拟。