电力系统三相短路仿真是很多电气工程专业课程设计和毕业设计里的必做项目。我在论坛上看到最多的提问是为什么我的仿真波形没有任何暂态冲击为什么短路电流曲线只有一条正弦曲线为什么把故障模块的切换时间填了模型运行后却报错这些问题暴露了一个共同点大部分教程只教“怎么搭建模型”没有讲清楚“三相短路暂态过程到底在 Simulink 中怎么体现”“仿真结果应该长成什么样”“如何用理论手算来验证仿真”。所以本文从物理原理讲到仿真实现从参数设置讲到结果验证把所有关键点串成一条线。先说结论用 Matlab Simulink 做三相短路仿真真正的难点不是建模型而是理解短路瞬间电流为什么不能突变、非周期分量如何衰减、冲击电流大小由什么决定以及如何把仿真波形和这些理论对应起来。只要理解了这些换个参数、换条线路你也能很快搭出可信的仿真模型。文章会包含三相短路基础概念、Simulink/Simscape Electrical 环境准备、完整建模步骤、可复制的初始化脚本和分析脚本、运行结果验证方法以及常见报错排查表。1. 这篇文章真正要解决的问题在电力系统分析课程中三相短路是必须掌握的内容。但很多同学做仿真时习惯按照教程把模块拖出来、参数照着填、运行完截一张 Scope 波形图交作业就算结束。这种做法丢失了最核心的部分——用仿真去验证理论模型而不是用仿真替代理论。这里有一个容易被忽略的大前提三相短路是电力系统中最严重的故障类型。虽然在实际电网中单相接地短路发生的概率更高但三相短路电流最大、对系统的暂态稳定影响最严重所以设备选型、保护整定、母线动热稳定校验都以三相短路作为计算工况。这也是为什么 MATLAB 课程设计、电力系统自动化专业课设、研究生复试里三相短路仿真出现频率这么高。本文要解决的几个具体问题第一短路电流为什么不是简单地由“电压除以电阻”得到。短路瞬间由于电感中磁链不能突变电流里会出现一个逐渐衰减的非周期分量使得短路后第一个周波的电流峰值明显高于稳态值。这个物理过程如果理解不深仿真波形就看不懂。第二Simulink 里的模块参数怎么设置。很多新手用 Three-Phase Fault 模块以为只要把 Fault 类型勾选为 A、B、C 三相就完事了却没有搞懂 Transition Time、接地电阻、故障点位置都会直接影响仿真结果。第三仿真结果怎么验证。不能仅仅说“波形看起来是对的”要能用短路理论公式对仿真数值做校核能解释冲击电流、稳态短路电流、衰减时间常数分别对应波形上的哪一段。读完这篇文章你可以独立完成一个无穷大电源系统三相短路的 Simulink 仿真并具备把仿真结果与手算理论值对比排查的能力。2. 三相短路暂态的核心概念与原理2.1 为什么三相短路最严重电力系统短路故障有四种常见类型三相短路、两相短路、两相接地短路、单相接地短路。从故障概率看单相接地故障最多但从故障严重程度看三相短路最严重。原因很简单三相短路时短路点的阻抗是对称的不存在零序分量故障电流最大。对设备而言最大冲击电流决定了断路器的开断能力对系统而言最严重故障下的稳定性决定了运行方式。故障类型 | 对称性 | 短路电流幅值 | 危害侧重点 三相短路 | 对称 | 最大 | 设备动热稳定、系统失稳 两相短路 | 不对称 | 约为三相的 0.87倍 | 负序电流引起发热 两相接地短路 | 不对称 | 视系统接地方式而定 | 零序保护动作 单相接地短路 | 不对称 | 可能较小 | 非故障相电压升高因此电力系统课程设计、继电保护整定计算、电气设备选型中都是先按三相短路计算最大短路电流再根据不同保护的配置原则考虑其他故障类型。这个背景决定了我们仿真建模时也要优先把三相短路模型做对。2.2 短路瞬间发生了什么短路发生前系统处于正常的正弦稳态。短路发生的瞬间相当于回路中突然并联了一个很小的阻抗电源电动势仍然存在但回路阻抗突然大幅减小。由于输电线路、变压器、发电机绕组都有电感电感中的电流不能突变所以短路电流不是从短路瞬间的负载电流直接跳到稳态短路电流而要经历一个过渡过程这就是“暂态”。如果从电路角度理解短路后的回路可以用一个 R-L 串联电路来近似。电感上的电流满足一阶微分方程。这个方程的解由两部分组成一个是与电源同频率的周期分量也叫稳态短路分量另一个是按指数规律衰减的非周期分量也就是直流分量。非周期分量的初始值大小取决于短路发生的瞬间电压的相位。如果短路发生在电压过零附近非周期分量初始值最大短路电流的第一个峰值最剧烈如果短路发生在电压峰值附近非周期分量初始值可能很小波形接近于纯正弦。这一点在仿真中非常重要也是很多人调 Fault 模块时间后发现波形冲击时大时小的原因。2.3 短路电流的关键指标三相短路暂态电流的经典表达式可以简化为i_k(t) I_pm * cos(ωt θ - φ_k) - I_pm * cos(θ - φ_k) * e^(-t/T_a)其中第一项是周期分量频率为工频 50Hz幅值保持不变第二项是非周期分量按时间常数 T_a 衰减T_a L / Rθ 是短路发生时刻电源电压的相位角φ_k 是短路回路的阻抗角。这里有两个必须记住的指标。第一个是冲击电流。短路后约半个周波50Hz 系统约 10ms时非周期分量尚未完全衰减与周期分量方向相同叠加成最大瞬时值称为冲击电流。工程上写成i_sh √2 * I * K_sh其中 I 是次暂态短路电流有效值K_sh 是冲击系数一般取 1.6 到 1.85取决于短路点电气距离的 R/X 比值。R/X 越小非周期分量衰减越慢冲击系数越大。第二个是稳态短路电流。当非周期分量衰减完毕短路电流进入稳定状态只剩基波周期分量。这个值就是电力系统短路计算中的稳定短路电流有效值一般用 I_inf 表示。2.4 容易混淆的几个“电流”很多初学者把瞬时值、有效值、冲击值搞混导致仿真结果对不上。这里做一个简单区分瞬时值是每一时刻电流的实际大小Scope 里看到的波形就是瞬时值。有效值是对一个周期做均方根得到的等效值常用电流表测量的就是有效值。冲击值是短路后第一个周波内的最大瞬时值不是有效值而是瞬时值的峰值。区别体现在如果仿真波形中冲击峰值为 10kA不代表短路电流有效值是 10kA要换算成有效值后才能和短路计算手册中的数值对比。3. Simulink 环境准备与建模前知识点3.1 需要安装的组件做电力系统仿真除了 MATLAB 基础模块和 Simulink 之外还需要安装电气仿真专用工具箱。在较新的版本中这个工具箱叫 Simscape Electrical在较早的版本中叫 SimPowerSystems 或 Simscape Power Systems。不管叫什么它在 Simulink Library Browser 中的路径一般是Simscape Electrical Specialized Power Systems如何确认自己有没有安装这个工具箱在 MATLAB 命令行窗口输入ver(simscape) ver(simscapeelectrical)如果返回版本信息说明已经安装。如果提示找不到该工具箱需要单独安装或更新 MATLAB 版本。一般情况下MATLAB 官网安装包中勾选 Simscape Electrical 即可版本以实际安装为准本文演示的模块路径在主流版本中均适用。3.2 Powergui 模块的作用从 Specialized Power Systems 库中拖入的第一个模块通常是 Powergui。很多新手不知道它的作用直接跳过不添加然后发现仿真报错或者波形异常。Powergui 是电力系统模型仿真环境的图形界面用来设置仿真类型、求解器参数、执行潮流计算、FFT 分析等。可以这样说没有 PowerguiSpecialized Power Systems 库中的模块无法正常工作。Powergui 的 Simulation types 有三种Continuous、Discrete、Phasor。Continuous连续仿真适合暂态过程分析也是本文推荐的方式。Discrete离散仿真需要设置采样时间适合实时仿真或与数字控制器联合仿真。Phasor相量仿真只适合稳态分析不能用于研究短路暂态波形。做三相短路暂态仿真时必须选择 Continuous否则无法看到非周期分量衰减过程。3.3 常用模块速查表这里先给出后面建模会用到的模块清单。记住这些模块在库中的路径可以少走很多弯路。模块名称 | 库路径 | 功能说明 Three-Phase Source | Simscape Electrical Specialized Power Systems | 三相电压源模拟电源 Three-Phase V-I Measurement | 同上的 Fundamental Blocks Measurements | 测量三相电压和电流 Three-Phase PI Section Line | 同上的 Fundamental Blocks Elements | 输电线路模型 Three-Phase Fault | 同上的 Fundamental Blocks Elements | 三相短路故障模块 Three-Phase Series RLC Load | 同上的 Fundamental Blocks Elements | 三相负载 Three-Phase Breaker | 同上的 Fundamental Blocks Elements | 三相断路器 Scope | Simulink Sinks | 波形显示 To Workspace | Simulink Sinks | 数据导出到工作区 Powergui | 同上的 Powergui | 仿真分析与求解器配置4. 三相短路 Simulink 建模核心流程4.1 新建模型与模型命名在 MATLAB 工具栏点击 Simulink选择 Blank Model保存为英文文件名比如three_phase_fault_demo.slx。为什么不建议用中文模型名因为在调用脚本、批量仿真、生成代码时中文路径和中文模型名容易出现兼容性问题。这是一个小的工程习惯养成了能省很多麻烦。4.2 搭建主电路拓扑示例系统采用一个简化的无穷大系统结构三相电源 → 测量模块 → 输电线路 → 故障点 → 负载。故障点设置在线路末端也就是负载的母线处。结构并不复杂但足够覆盖短路暂态仿真的所有要点。搭建顺序建议从左到右。先放置 Three-Phase Source然后放 Three-Phase V-I Measurement再放 Three-Phase PI Section Line然后是 Three-Phase Fault最后接一个 Three-Phase Series RLC Load。Powergui 单独放在模型中任意空白位置即可它与主电路之间没有连线关系。4.3 配置三相电源双击 Three-Phase Source主要配置项有Phase-to-phase voltage线电压有效值本文示例用 110kVFrequency50HzInternal connection选择 Yg星形接地这是电力系统中最常见的电源接法Internal Resistance、Internal Inductance电源内阻抗。电源内阻抗不能设置为 0否则短路电流理论上会趋于无穷大仿真会报数值错误。内阻抗的取值可以参考系统短路容量来折算示例中取一个很小的值只是为了限制数值溢出实际工程中应由系统等值电路计算。4.4 配置输电线路输电线路用 Three-Phase PI Section Line 模块。双击后可以设置Frequency used for RLC specification50HzLength线路长度示例为 50kmPositive- and zero-sequence resistances、inductances、capacitances正序和零序参数。这里有一个必须注意的点做三相短路仿真时正序参数决定短路电流大小。如果只做三相短路可以不关心零序参数但如果后续要扩展成单相接地短路仿真零序参数就必须认真填写。线路参数可以用仿真所在地区的典型架空线参数也可以从电力系统教材附录中查找关键是把单位换算正确电阻单位是 Ω/km电感单位是 H/km电容单位是 F/km。4.5 配置三相故障模块Three-Phase Fault 是本次仿真的核心模块。双击打开后关键配置项是Fault type三相短路时勾选 A、B、C 三相Fault resistance短路点的接地电阻金属性短路可以填一个很小的值比如 0.001ΩTransition Time切换时间以向量形式填写例如 [0.05 0.4]表示 0.05s 时发生短路0.4s 时故障切除。Transition Time 非常容易理解错。它不是“故障开始时间”而是“状态切换时间点”。故障模块在初始状态是打开的不短路到达第一个时间点后切换到短路状态到第二个时间点后再切回来。如果只填一个时间点比如 [0.05]模块不会自动恢复。这个细节是很多报错和波形异常的来源。4.6 配置测量、显示和结果存储在电源出口处串联一个 Three-Phase V-I Measurement它有两个输出端口一个是测量到的电压一个是电流。将电流输出连接到 Scope就可以观察到三相短路电流波形。但 Scope 中的波形只能肉眼看不能量化分析。为了让数据进入 MATLAB 工作区更推荐的做法是从测量模块的电流输出端再引一路到 To Workspace 模块。To Workspace 模块的 Variable Name 设置为IabcSave format 选择 Timeseries这样仿真结束后工作区中会出现一个包含三相电流时间序列的变量方便后续做 FFT、求有效值、算冲击系数。4.7 设置仿真参数在 Simulink 工具栏点击 Simulation Model Configuration Parameters求解器类型选择 Variable-step求解器建议用 ode23tb相对容差设置为 1e-4 或更小。为什么推荐 ode23tb电力电子和电力系统暂态过程往往是刚性问题短路瞬间电流变化率很大ode45 在部分情况下会花很长时间甚至不收敛ode23tb 对这类问题更稳定。Stop Time 设置为 0.6s这样故障前有 0.05s 的稳定运行段故障后有 0.35s 的暂态衰减段能完整观察非周期分量的变化趋势。5. 完整示例无穷大系统三相短路仿真5.1 示例系统参数这个示例的目标是搭建一个 110kV、50Hz 的无穷大电源系统在 0.05s 时于线路末端发生三相金属性短路0.4s 时切除故障观察短路电流的暂态过程。参数取值如下表实际项目请按真实系统替换参数项 | 数值 电源线电压有效值 | 110 kV 电源频率 | 50 Hz 电源内电阻 | 0.01 Ω 电源内电感 | 0.0001 H 线路长度 | 50 km 线路正序电阻 | 0.1153 Ω/km 线路正序电感 | 0.00105 H/km 线路正序电容 | 11.33e-9 F/km 三相负载有功功率 | 10 MW 三相负载无功功率 | 1 Mvar 故障电阻 | 0.001 Ω 故障切换时刻 | [0.05 0.4] 总仿真时间 | 0.6 s5.2 初始化脚本推荐把系统参数写在一个独立的脚本中而不是直接填在模块界面里。这样做的好处是一次定义多处引用后期改参数时不需要逐个模块去改而且不容易出错。新建一个init_short_circuit.m内容如下% init_short_circuit.m % 三相短路仿真系统初始化脚本 % 在 Simulink 模型的 PreLoadFcn 回调中调用 freq 50; Vll 110e3; % 线电压有效值单位 V Vph Vll / sqrt(3); % 相电压有效值单位 V % 电源内阻抗 Rs 0.01; % 电源内电阻单位 Ω Ls 0.0001; % 电源内电感单位 H % 线路正序参数按总长度计算 line_length 50; % 线路长度单位 km R_line 0.1153 * line_length; % 线路总电阻单位 Ω L_line 0.00105 * line_length; % 线路总电感单位 H C_line 11.33e-9 * line_length; % 线路总电容单位 F % 负载参数 P_load 10e6; % 有功功率单位 W Q_load 1e6; % 无功功率单位 Var % 由 R/L 估算非周期分量衰减时间常数 Ta (Ls L_line) / (Rs R_line); fprintf(系统初始化完成\n); fprintf( 相电压有效值 %.2f V\n, Vph); fprintf( 线路总电阻 %.4f Ω\n, R_line); fprintf( 线路总电感 %.5f H\n, L_line); fprintf( 非周期分量时间常数 ≈ %.3f ms\n, Ta * 1000);然后在模型窗口的 Model Properties Callbacks PreLoadFcn 中填入init_short_circuit。之后每次打开模型工作区都会自动得到这批参数。5.3 模块参数与变量关联为了让初始化脚本中的参数真正生效模块参数栏中不要写数字而是写变量名。例如Three-Phase Source 的 Internal Resistance 填RsInternal Inductance 填LsThree-Phase PI Section Line 的 Length 填line_length电阻、电感、电容也可以引用R_line、L_line、C_lineThree-Phase Series RLC Load 的 Active Power 填P_loadInductive Reactive Power 填Q_load。这样做以后改参数时只需要改init_short_circuit.m中的数值再重新运行一次脚本模型就会用新参数仿真。这是一个很实用的建模习惯不仅适用于短路仿真也适用于变压器、电机等任何带参数的电力系统模型。5.4 运行仿真在模型搭建完成参数配置好之后点击 Run 或在命令行执行sim(three_phase_fault_demo)如果 To Workspace 模块的变量名为IabcStop Time 设置为 0.6s仿真结束后工作区会自动出现Iabc变量。可以用下面的命令快速查看plot(Iabc.time, Iabc.signals.values)第一列是时间第二列起是 A、B、C 三相电流瞬时值。如果波形显示正常就得到了短路暂态的原始数据。6. 运行结果与效果验证6.1 预期波形特征正常情况下的短路电流波形应该有三个特征故障发生前电流是负载电流幅值不大。0.05s 瞬间电流突然增大并且第一个周波的峰值明显高于后续稳态峰值。随后非周期分量逐渐衰减电流波形趋于关于横轴对称的标准正弦波。如果切换时间设置为 [0.05 0.4]那么在 0.4s 时电流会突然下降之后带负载恢复正常。三相电流的暂态过程并不完全相同。由于短路瞬间各相电压相位不同非周期分量的初始值不同所以某一相可能冲击很明显另一相冲击较小。这是正常现象不要误以为模型设置错误。6.2 用理论值验算仿真值短路回路的近似总阻抗为R_sum Rs R_line 0.01 0.1153 * 50 5.775 Ω X_sum 2 * pi * 50 * (Ls L_line) 2 * pi * 50 * 0.0526 ≈ 16.53 Ω Z_sum sqrt(R_sum^2 X_sum^2) ≈ 17.51 Ω稳态短路电流理论值约等于相电压有效值除以短路阻抗幅值也就是 63.51kV / 17.51Ω ≈ 3.63kA。在仿真中读取故障后 0.2s 到 0.4s 之间的电流数据做均方根运算结果应该在 3.6kA 左右。如果偏差超过 10%就要检查线路参数、故障电阻或仿真步长设置。6.3 FFT 分析脚本均方根法求出的有效值包含非周期分量残留更严格的做法是用 FFT 提取基波分量。在实际应用中短路电流周期分量的有效值是校验设备热稳定的重要依据。新建analyze_short_circuit.m% analyze_short_circuit.m % 读取 Iabc 数据计算三相短路稳态电流的有效值 % 使用前请先运行仿真确保工作区存在 Iabc 变量 i_a Iabc.signals.values(:, 1); t Iabc.time; % 选取故障后 0.25s 到 0.4s 的数据段此时非周期分量已衰减 mask (t 0.25) (t 0.4); data i_a(mask); t_data t(mask); % 默认采样间隔 dt median(diff(t_data)); fs 1 / dt; % 去直流分量后做 FFT data_clean data - mean(data); N length(data_clean); spectrum fft(data_clean .* hann(N)); freq_axis (0:N-1) / N * fs; % 找 50Hz 附近的幅值 [~, idx] min(abs(freq_axis - 50)); amp_50 2 * abs(spectrum(idx)) / N; irms_50 amp_50 / sqrt(2); fprintf(稳态短路电流基波有效值 %.1f A\n, irms_50); fprintf(相对理论值误差 %.2f%%\n, abs(irms_50 - 3.63e3) / 3.63e3 * 100);如果 FFT 计算出的基波有效值和理论值对得上说明模型搭建和参数设置是可信的。这里有一个提醒FFT 对采样均匀性敏感如果 To Workspace 模块采用 Variable-step 输出数据点可能不均匀。建议在 To Workspace 模块的 Sample Time 中填一个固定步长比如 1e-4或者在信号进入 To Workspace 之前串联一个 Zero-Order Hold 模块。6.4 如何判断成功判断仿真是否成功不能只看 Scope 有没有波形建议按顺序检查三件事一是波形形态。短路瞬间电流是否明显增大第一个峰值是否高于后续峰值非周期分量是否在几百毫秒内衰减完毕。二是数值量级。稳态短路电流是否与理论手算结果接近误差是否在可接受范围。三是故障切除现象。0.4s 后电流是否恢复正常负载电流如果能观察到故障切除后的暂态过程说明 Fault 模块的切换时间设置正确。如果波形和数值都符合本次仿真才算真正完成。7. 常见问题与排查方法下面是电力系统 Simulink 三相短路仿真中最常见的问题按出现频率从高到低排列。问题现象 | 可能原因 | 排查方式 | 解决方案 模型运行报错未添加 Powergui | Specialized Power Systems 模块缺少仿真环境 | 检查模型中是否有 Powergui | 从库中拖入 Powergui仿真类型选 Continuous 波形中没有短路冲击 | Fault 模块的 Transition Time 设置错误 | 双击 Fault 模块确认向量时间点 | 正确设置切换时间如 [0.05 0.4] 短路电流数值为 0 或很小 | 测量模块接线错误或 Fault 类型未勾选三相 | 检查 V-I Measurement 的电流信号是否接入 Scope | 将 Fault 类型勾选为 A、B、C且接地电阻设为极小值 仿真非常慢或长时间不收敛 | 求解器不适合刚性电气暂态过程 | 查看求解器类型和容差设置 | 改用 ode23tb相对容差设为 1e-4 波形出现明显高频振荡 | Powergui 选择的是 Discrete步长过大 | 查看 Powergui 的 Simulation type | 改用 Continuous或减小离散采样时间 To Workspace 中保存的变量无法做 FFT | 数据点不等间隔 | 查看数据时间序列的 diff 是否恒定 | 在 To Workspace 中设置固定采样时间这里重点说一下仿真不收敛的问题。在短路瞬间电流变化率很大如果求解器步长控制不好会出现数值振荡或者计算时间极长。解决办法是使用 ode23tb、相对容差缩小到 1e-5、限制最大步长不超过 1e-4。但要注意最大步长设置过小会大幅增加仿真时间在仿真精度和速度之间要取得平衡。8. 最佳实践与工程建议8.1 参数变量化前面提到了将模块参数绑定到初始化脚本中的变量。这个做法的收益在实践中比想象中大得多。当你做电网参数敏感性分析、多故障场景批量仿真时只需要在循环里修改初始化脚本中的变量再调用sim命令就可以自动完成几十组工况的仿真不需要手动在模块界面上改参数。这也是 Matlab/Simulink 平台相比纯图形化仿真软件的核心优势之一。8.2 数据导出与分析Scope 波形适合人眼观察但不适合后期输出论文插图或生成计算报告。建议所有关键信号都通过 To Workspace 导出然后编写统一的绘图脚本来出图。绘图时注意设置线宽、字体、坐标轴标签这样得到的图片可以直接用于课程设计报告。对于短路电流波形建议同时画出三相并且用xline(0.05)、xline(0.4)在故障发生和切除时刻画竖线读者一眼能看出故障区间。8.3 仿真步长与数据采样在选择仿真步长时有一个常见误区不是步长越小越好。步长过小会显著增加仿真时间步长过大又会丢失短路瞬间的细节。对于 50Hz 工频系统建议至少保证每个工频周期有 50 到 100 个采样点。如果使用固定步长1e-4s 左右可以满足大多数三相短路暂态分析的精度需求。如果做 FFT固定步长是必要条件这一点务必记住。8.4 模型规范化与团队协作如果是课程设计小组或工程项目多人协作模型命名、模块命名、信号命名都要规范化。模块名不要用Subsystem1、Gain2这种默认名建议用能表达功能的名称比如Fault_At_Line_End、Load_Bus_Current。这一步有利于代码审查、仿真调试和后期保护逻辑设计。模型文件中不要放多余的中文注释变量尤其是在跨平台使用时避免编码问题。8.5 安全提醒每次谈到电力系统短路仿真都值得强调一次安全边界真实电力系统短路试验属于破坏性试验只能由电网调度和检修专业人员在高电压试验场或有严格安全措施的继电保护测试系统中进行普通电气工程师和学生绝不能在任何场合模仿“制造短路”。Simulink/Matlab 仿真的价值之一就是用数值计算替我们完成危险又昂贵的物理试验因此不要因为仿真中轻易实现了短路故障就低估真实系统的危险性。9. 总结与后续学习方向三相短路 Simulink 仿真看似是电力系统分析课程的一个小作业但它串联起了电路暂态分析、电力系统对称分量法、电气设备校验、继电保护整定等多个知识点。真正跑通模型、验证完数据之后你会发现自己读短路计算手册时不再觉得公式抽象了。如果你已经把本文的示例模型跑通了下一步可以从三个方向深入第一个方向是不对称短路。把 Fault 模块改为单相接地故障或两相短路改用对称分量法分析对比不同故障类型下短路电流的差异。这时候线路的零序参数会变得很关键需要补充完整。第二个方向是暂态稳定分析。在无穷大电源模型中加入同步发电机模块观察短路期间发电机功角的变化曲线分析临界切除时间。这是电力系统暂态稳定课程设计的经典内容。第三个方向是继电保护仿真。在故障点两端加电流互感器和距离保护模型验证保护动作逻辑和跳闸时间这会让你的仿真成果更接近工程应用。建议把你搭建好的模型做一次完整参数整理并且把验证脚本和模型一起保存。这些代码和参数在你之后做课程设计、毕业论文甚至实际项目时都可以当作可复用的仿真基础模块。