PMP考点怎么选？看完这篇不纠结！
2026/7/10 8:37:57
目录
1.设计题目
2.设计要求
3.基础资料（用于仿真程序验证）
4.设计内容
三、设计内容
1.仿真程序流程图
3.仿真程序使用说明
四、设计总结
1.设计过程中遇到的问题及解决方法
2.设计体会
3.对设计的建议
五、参考文献
1.设计题目
两相短路故障短路电流仿真程序设计
2.设计要求
采用计算机算法计算故障点不对称故障电流、短路功率等参数。需绘制计算机算法的流程图，并基于流程图编写相应计算程序。
3.基础资料（用于仿真程序验证）
电力系统图1所示，输电线路、电力变压器、系统、发电机及负载的参数均标注于图中，系统图f点发生不对称短路（b、c两相短路），分析与计算短路电流与短路功率（不对称故障的短路功率
为故障电流有效值）。变压器连接组标号均为YNd11。各元件三序参数可参考后面附表。电动机三序参数可参考教材表10.1，或者按照祝淑萍的书135页略去不计。距离较远，电动机电势对短路电流的影响可忽略。本题复合序网只有正序和负序，故计算元件参数时，可只计算正序和负序参数。关于线路的几何均距，参考后面表格自己选择即可。
4.设计内容
（1）设计两相短路故障计算程序流程图。
（2）编写基于节点导纳矩阵的系统两相短路相关参数计算程序，该程序应实现以下主要功能：
①基于网络各序标幺值等值参数，采用导纳矩阵求逆法求取系统的正、负序节点阻抗矩阵；
②针对任意给定的故障点，计算其在发生任意两相短路时的故障点正、负序电流、故障点短路电流以及短路功率。
③在故障点发生两相短路的条件下，计算任意指定支路的各相电流和任意指定节点的各相电压。
（3）以3基础资料所述系统为例，以100 MVA为基准容量，采用平均额定电压作为基准电压，计算各元件的正、负序标幺值等值参数，据此绘制标注各节点序号的系统正、负序等值电路(请自行对各节点进行编号)，采用所编写程序对其仿真分析，验证程序正确性。
三、设计内容
1.仿真程序流程图
2.仿真程序源代码
%%节点阻抗矩阵函数 % B Z 支路参数矩阵 n1*6 % Y 节点导纳矩阵 function [Z]=T3(B,X) nl=length(B(:,1)); n=length(X(:,1)); Y=zeros(n,n); for i=1:nl if(B(i,6)==1) L=B(i,1); R=B(i,2); else L=B(i,2); R=B(i,1); end Y(L,L)=Y(L,L)+1/B(i,3)*B(i,5)^2+B(i,4)/2; Y(R,R)=Y(R,R)+1/B(i,3)+B(i,4)/2; Y(L,R)=Y(L,R)-1/B(i,3)*B(i,5); Y(R,L)=Y(L,R); end for i=1:n Y(i,i)=Y(i,i)+X(i,2); end Z=inv(Y); end B=[1 2 0.058i 0 0.9565 1; 2 3 0.329i 0 1 0; 3 4 0.1665i 0 0.9565 0; 3 6 0.161i 0 1 0; 6 5 1.05i 0 1.00435 0]; X1=[1 1/0.204i; 2 0; 3 1/6.321i; 4 1/0.208i; 5 1/8.5i; 6 0]; X2=[1 1/0.204i; 2 0; 3 1/3.321i; 4 1/0.208i; 5 1/8.5i; 6 0]; %正序阻抗矩阵 Z1=T3(B,X1); %负序阻抗矩阵 Z2=T3(B,X2); %正负序电流 I_1=1.12/(Z1(3,3)+Z2(3,3)) I_2=-I_1 %短路电流标幺值 有名值 If_1=abs(I_1)*sqrt(3) If=If_1*100/sqrt(3)/115*1000 %短路功率 P=100*If_1 %输入节点位置 i=3; j=2; %%各节点电压 %正序 U_0=1.12; U_i_1=U_0-Z1(i,3)*I_1 Zjk_1=0.0887i; U_j_1=U_0-Z1(j,3)*I_1 %负序 U_i_2=-Z2(i,3)*I_2 Zjk_2=0.0862i; U_j_2=-Z2(j,3)*I_2 %%各支路电流正序 %正序 I_ij_1=(U_i_1-U_j_1)/0.329i %负序 I_ij_2=(U_i_2-U_j_2)/0.329i %%各相电流 a=-0.5+1i*sqrt(3)/2; I_a=I_ij_1+I_ij_2 I_b=a^2*I_ij_1+a*I_ij_2 I_c=a*I_ij_1+a^2*I_ij_2 %%各相电压 U_a=U_i_1+U_i_2 U_b=a^2*U_i_1+a*U_i_2 U_c=a*U_i_1+a^2*U_i_2 %%判断支路是否存在变压器 %t是判断条件 t=0; % 角度转弧度 theta_pos = deg2rad(+30); % 正序：+30°（YNd11） theta_neg = deg2rad(-30); % 负序：-30° if(t==1) disp('支路存在变压器，从高压侧到低压侧'); %正序参数相位偏移 U_i_1=U_i_1*exp(1i*theta_pos) U_j_1=U_j_1*exp(1i*theta_pos) I_ij_1=I_ij_1*exp(1i*theta_pos) %负序参数相位偏移 U_i_2=U_i_2*exp(1i*theta_neg) U_j_2=U_j_2*exp(1i*theta_neg) I_ij_2=I_ij_2*exp(1i*theta_neg) %%各相电流 a=-0.5+0.866i; I_a=I_ij_1+I_ij_2 I_b=a^2*I_ij_1+a*I_ij_2 I_c=a*I_ij_1+a^2*I_ij_2 %%各相电压 U_a=U_i_1+U_i_2 U_b=a^2*U_i_1+a*U_i_2 U_c=a*U_i_1+a^2*U_i_2 end3.仿真程序使用说明
(一)、程序功能
该MATLAB仿真程序用于电力系统不对称BC两相短路故障计算，依托自定义T3函数实现导纳矩阵求逆法求解序网节点阻抗矩阵。程序先录入电网支路、节点对地序参数，分别搭建正序网络、负序网络节点导纳矩阵，对导纳矩阵求逆得到正负序节点阻抗矩阵；再输入故障节点编号，根据两相短路故障边界条件，计算故障点正序、负序电流标幺值，推导故障短路电流、短路容量；接着遍历全部网络节点，计算所有节点的正序、负序序电压，提取任意支路的正负序支路电流；最后利用对称分量变换公式，将序电压、序电流转换为实际A、B、C三相电压与三相电流，完整实现含变压器、发电机、线路的多节点电网不对称短路数值仿真，可用于课程设计故障数值验算与电网故障特性分析。
(二)、输入参数、输入矩阵及各元素含义
1.支路参数矩阵B
该矩阵为核心拓扑输入矩阵，每一行对应电网一条电气支路，固定六列参数。第一列代表支路起始节点编号，第二列代表支路终止节点编号，用来确定电网各元件的拓扑连接关系；第三列是支路串联电抗的复数标幺值，本次仿真无电阻仅保留虚数电抗；第四列为线路对地充电电纳，本次模型不考虑电容，全部填充0；第五列是变压器变比标幺值，普通输电线路无变压器，统一填1，变压器支路填入折算后的标准变比；第六列为接地标志数字，数值1代表该支路连接接地元件，数值0代表支路两端均不接地。
2.节点对地正序导纳矩阵X1
矩阵为两列多行结构，用来录入每个节点对地的正序电气参数。每行第一列是节点编号，第二列是该节点对地元件的正序导纳标幺值，发电机节点填入发电机正序电抗倒数，无接地元件的线路、负荷节点直接填0，用于函数生成完整正序导纳矩阵。
3.节点对地负序导纳矩阵X2
矩阵格式与X1完全一致，仅第二列为节点对地负序导纳标幺值，发电机、电动机正负序电抗数值不同，对应节点导纳数值会发生变化，无源节点依旧填充0，用于构建负序网络导纳矩阵。
4.单独输入常量参数
包含系统基准容量、基准电压、发电机空载电动势标幺值、故障节点编号，以及对称分量变换计算所需的复数算子a；若网络含YNd11变压器，还需输入相位校正系数，用于修正变压器带来的序分量相位偏移。
(三)、输出参数、输出矩阵及各元素含义
1.正序节点阻抗矩阵Z1、负序节点阻抗矩阵Z2
二者均为和网络节点数量一致的方阵，矩阵行号、列号均对应电网节点编号。矩阵对角线上的元素Z(k,k)为第k个节点的序自阻抗，物理含义是仅在k节点通入单位电流时，k节点自身产生的序电压；非对角线元素Z(m,n)为m节点与n节点之间的序互阻抗，代表n节点通入单位电流时m节点产生的序电压，完整反映网络各节点之间的电气耦合关系。
2.故障点电气量输出参数
输出故障点正序电流、负序电流复数标幺值，BC两相短路工况下二者幅值相等；输出故障短路电流标幺值与有名值、短路容量标幺值，直观反映短路故障严重程度。
3.节点、支路序电气量输出参数
输出全网每个节点对应的正序电压、负序电压复数标幺值；输出任意指定支路流过的正序支路电流、负序支路电流复数标幺值，复数同时包含电气量幅值与相位信息。
4.三相电气量输出参数
依托对称分量公式换算得到A相、B相、C相三相电压、三相电流复数，可直接提取幅值与相位，直观体现不对称短路带来的三相电压跌落、三相电流不平衡问题；存在变压器时，额外输出相位校正后的修正版三相电压、三相电流，消除变压器绕组相位带来的计算误差。
(四)、自定义函数T3输入输出补充说明
函数输入：第一个输入B为支路参数矩阵，第二个输入X为对应序网的节点对地导纳矩阵（X1或X2）；
函数输出：Z为对应序网完整节点阻抗方阵，矩阵行列对应系统各节点编号。
函数内部逻辑：先搭建节点导纳矩阵，逐条录入支路的自导、互导，叠加对地电容与发电机对地导纳，最后对导纳矩阵求逆得到节点阻抗矩阵输出。
4.仿真程序验证
四、设计总结
1.设计过程中遇到的问题及解决方法
（1）、在绘制标注各节点序号的系统正负序等值电路时，有问题。
解决办法：通过和老师及交流和查询课本例题找到错误。
（2）、仿真时，节点导纳矩阵等于零，写入的数据有误，矩阵X第二列，接地导纳值输入有问题。
解决办法：灵活运用Chatgpt 、deepseek和Gemini等等以及和老师交流查询课本资料找到错误并修改。
（3）、计算相关变压器支路时，没考虑到相位偏移。
解决办法：通过查询教材上公式和资料得知正负序电压电流经过星三角形变压器的相位偏移情况。
2.设计体会
本次两相短路故障仿真程序设计，是我第一次完整经历“理论推导→参数计算→程序编写→结果校验”的电力系统工程实践，收获远超预期。
最初我以为只是把课本公式敲进MATLAB，实际却踩了不少“隐形坑”。比如绘制正、负序等值电路时，我一开始忽略了变压器YNd11接线对相位的影响，导致支路电流计算结果与理论值有偏差；后来在搭建节点导纳矩阵时，又误将节点对地导纳填为0，导致矩阵奇异、程序报错，反复排查后才发现是发电机次暂态电抗的导纳值漏写。这些问题让我深刻意识到：电力系统计算容不得半点马虎，一个参数的疏漏就会导致整个结果失效。
最让我有成就感的是对“仿真结果与理论边界条件差异”的分析。一开始我发现短路点a相电压接近0，和课本中“两相短路a相电压不为零”的结论矛盾，一度以为是程序错误。后来通过查阅《电力系统故障分析》，才明白课本结论是“理想不接地系统”的假设，而实际工程中110kV系统多为中性点接地，短路点电位会被钳位到地电位。这次“纠错”让我真正打通了“理论-仿真-工程”的逻辑链条：仿真不是对理论的简单复刻，而是要反映真实系统的复杂性。
此外，我也体会到模块化编程的优势。把导纳矩阵构建、阻抗矩阵求逆、序分量合成拆分成独立模块后，调试效率提升了不少，也更理解了计算机算法在电力系统大规模计算中的必要性——如果用手算这个6节点系统，至少需要半天，而程序几秒就能出结果。这次设计不仅巩固了不对称故障分析的知识，更培养了我的工程严谨性和问题解决能力，为后续的毕业设计和工程实践打下了坚实基础。
3.对设计的建议
（1）、 增加“理论-仿真差异”的引导分析
本次设计中“短路点a相电压接近0”的问题困扰了我很久，后来才明白是理想模型与实际系统的差异。建议在教学中专门增设一节“仿真结果合理性分析”，引导学生区分“教科书理想假设”和“工程实际场景”，比如讲解中性点接地方式对短路电压的影响，避免学生把“数值差异”简单归为“程序错误”，培养工程思维。
（2）、 优化参数查询的便利性
本次设计需要查询线路电抗、变压器阻抗等参数，附录中的表格排版较乱，查找效率低。建议将常用参数整理成电子表格，按电压等级、设备类型分类索引，或提供在线参数查询工具，减少学生在非核心任务上的时间消耗。
（3）、增加“变压器相位偏移”的专项练习
YNd11变压器的±30°相位偏移是不对称故障计算的重点和难点，本次设计中我差点遗漏这一环节。建议在设计任务中明确要求对比“考虑相位偏移”和“不考虑相位偏移”的结果差异，加深学生对变压器接线组别影响的理解。
（4）、开放“扩展故障类型”的选做任务
本次设计仅要求实现两相短路仿真，对学有余力的同学可以开放单相接地、两相接地短路等扩展任务，鼓励对比不同故障类型的序网连接方式和计算结果，进一步提升综合能力。
五、参考文献
[1] 何仰赞,温增银.电力系统分析[M].4版.武汉:华中科技大学出版社,2016.
[2] 陈珩.电力系统稳态分析[M].4版.北京:中国电力出版社,2018.
[3] 李光琦.电力系统暂态分析[M].3版.北京:中国电力出版社,2012.
[6] 刘万顺.电力系统故障分析[M].3版.北京:中国电力出版社,2010.
[7] DL/T 5429-2009 电力系统设计技术规程[Z].北京:中国电力出版社,2009.
[9] 电力工业部电力规划设计总院.电力系统设计手册[Z].北京:中国电力出版社,1998.
[10] 孟祥萍，高燕.电力系统分析.3版. 北京:高等教育出版社，2021.5