三相短路短路电流计算(matlab程序)
基于MATLAB短路电流的计算及仿真
基于MATLAB短路电流的计算及仿真短路电流是指在电力系统中，当发生短路时，电流可以达到的最大值。
短路电流的计算和仿真对于电力系统的设计和运行具有重要的意义。
本文将介绍基于MATLAB的短路电流计算和仿真方法。
短路电流的计算可以通过解析法和数值法两种方法进行。
解析法是通过解闭合方程组得到短路电流的解析解，而数值法则是通过迭代计算来逼近短路电流的数值解。
首先，我们将介绍基于解析法的短路电流计算方法。
在电力系统中，短路电流通常可以用下式表示：Isc = U / Z其中，Isc为短路电流，U为系统电压，Z为系统的等值阻抗。
为了计算短路电流，我们首先需要确定系统的等值阻抗。
等值阻抗通常可以通过系统的参数和网络拓扑来计算。
在MATLAB中，可以使用电力系统仿真工具箱（Power System Toolbox）来计算等值阻抗。
具体的步骤如下：1. 创建一个空的系统模型：使用“new_system”函数创建一个新的系统模型。
2. 导入系统参数：使用“psat”函数将系统参数导入系统模型中。
3. 创建系统配置：使用“psscfg”函数创建一个系统配置，并设置相关参数。
4. 计算等值阻抗：使用“psadeqz”函数计算等值阻抗，并保存结果。
5.计算短路电流：根据上述公式计算短路电流。
以上就是基于解析法的短路电流计算的主要步骤。
通过这种方法，我们可以快速准确地计算出系统的短路电流。
接下来，我们将介绍基于数值法的短路电流计算方法。
数值法通常采用迭代过程来逼近短路电流的数值解。
在MATLAB中，可以使用电力系统仿真工具箱（Power System Toolbox）中的短路计算功能来进行数值计算。
具体的步骤如下：1. 创建一个空的系统模型：使用“new_system”函数创建一个新的系统模型。
2. 导入系统参数：使用“psat”函数将系统参数导入系统模型中。
3. 创建系统配置：使用“psscfg”函数创建一个系统配置，并设置相关参数。
电力系统三相短路的编程计算
电力系统分析课程专题报告学生姓名：班级：学号：指导教师：所在单位：提交日期：评分电力系统三相短路的编程计算某某某（学院， ）摘要：在电力系统中，三相短路故障造成的危害是最大的，发生的几率也最高，故短路计算对电力系统的稳定运行具有十分重要的意义。
作为电力系统三大计算之一，计算三相短路故障发生时的短路电流、各节点电压、各支路电流是短路计算的基本内容。
在电力系统短路电流的工程计算中，由于快速继电保护的应用，最重要的是计算短路电流基频交流分量的初始值，即次暂态电流I ''。
在给定电源电势时，实际上就是稳态交流电路的求解。
本文基于教材例3-2应用MATLAB 编程计算三相短路故障的电流电压情况，并于例3-3进行进一步的验证和完善。
关键词：三相短路；MATLAB中图分类号：TM 713 文献标识码：A0 引言在电力系统的四种短路类型中，三相短路是其中最严重的，其短路电流可达数万安以至十几万安，随之产生的热效应和电动力效应将使电气设备遭受严重破坏。
因此，计算短路电流的主要应用目的是电力系统设计中的电气设备选择，短路计算已成为电力系统运行分析、设计计算的重要环节。
实际电力系统短路电流交流分量初始值的计算，小型系统可以手算，而对于结构复杂的大型系统，短路电流计算量较大，用计算机进行辅助计算成为大势所趋。
1 解析法求三相短路电流1.1 参数说明(1) 为了元件参数标幺值计算方便，取基准容量B S 为A MV ⋅60，可设任意值，但必须唯一值参与计算。
(2) 取基准电压B U 为平均额定电压AV U ,基于例3-2中系统的额定电压等级有kv 10、kv 110，平均额定电压分别为kv 115、kv 5.10，平均额定电压与线路额定电压相差5%，为简化计算，故取平均额定电压。
(3) I ''为次暂态短路电流有效值，短路电流周期分量的初值等于时间0=t 时的有效值。
满足产生最大短路电流的三个条件下的最大次暂态短路电流作为计算依据。
MATLAB应用基础短路计算仿真
南京工程学院MATLAB应用基础(作业)题目题目八短路计算仿真：work1.mdl课程名称 MATLAB应用基础院（系、部、中心）电力工程学院专业电力系统及其自动化班级 K电力122学生姓名于湘唐(42) 张译天(43)赵晨(44) 周博梵(45)朱高斐(46)无穷大功率电源供电系统发生三相短路故障MATLAB 模拟1.1 实例在图1-1-1所示的网络中，当降压变电所10.5kV 母线上发生了三相短路时，可将系统视为无限大容量电源，试求此时短路点的冲击电流imp i ，短路电流的最大有效值imp I 和短路功率kt S 。
图1-1-11.2 理论计算过程解 取MV A 100B =S 、n av B •=U U ，已知km 401Ω=.x 。
首先计算各元件参数的标幺值电抗()525020100100510100N B k 1..S S %U X *=⨯==292037100104022n av B 12..U S l x X *=⨯⨯==•()192231001007100%N B k 43..S S U X X **=⨯=⨯==取1=*E 作成等值网络如图1-2-1所示。
图1-2-1等值网络图短路回路的等值电抗为91211922129205250....X *=⨯++=Σ短路电流周期分量的有效值为5230912111..X I **===∑ω ()kA 88251031005230B ...I I I *=⨯⨯==ωω若取冲击系数81imp .K =，则冲击电流为()kA 347882552281imp ...I .i =⨯=⨯=ω短路电流的最大有效值为()kA 384882521521imp ...I .I =⨯==ω短路功率为()MV A 3521005230B k t ..S I S *=⨯==ω2.1 无穷大功率电源供电系统仿真模型构建假设无穷大功率电源供电系统如图2-1-1所示，在0.02s 时刻变压器低压母线发生三相短路故障，仿真其短路电流周期分量幅值和冲击电流的大小。
基于matlab的电力系统故障分析与仿真(毕业论文毕业设计)[管理资料]
![基于matlab的电力系统故障分析与仿真(毕业论文毕业设计)[管理资料]](https://img.taocdn.com/s3/m/fc5da7eebed5b9f3f80f1c30.png)
毕业设计（论文）基于MATLAB的电力系统故障分析与仿真学号：姓名：专业：电气工程及其自动化系别：指导教师：二〇一三年六月毕业设计（论文）基于MATLAB的电力系统故障分析与仿真学号：姓名：专业：电气工程及其自动化系别：指导教师：二〇一三年六月北京交通大学毕业设计（论文）成绩评议题目：基于MATLAB的电力系统故障分析与仿真系别：专业：电气工程及其自动化姓名：学号：指导教师建议成绩：84评阅教师建议成绩：86答辩小组建议成绩：82总成绩：84答辩委员会主席签字：年月日北京交通大学毕业设计（论文）任务书北京交通大学毕业设计（论文）开题报告北京交通大学毕业设计（论文）指导教师评阅意见北京交通大学毕业设计（论文）评阅教师评阅意见北京交通大学毕业设计（论文）答辩小组评议意见毕业设计（论文）诚信声明本人声明所呈交的毕业设计（论文），是本人在指导教师的指导下，独立进行研究工作所取得的成果，除了文中特别加以标注和致谢中所罗列的内容以外，毕业设计（论文）中不包含其他人已经发表或撰写过的研究成果，也不包含为获得北京交通大学或其他教育机构的学位或证书而使用过的材料。
申请学位论文与资料若有不实之处，本人承担一切相关责任。
本人签名：日期：毕业设计（论文）使用授权书本人完全了解北京交通大学有关保管、使用论文的规定，其中包括：①学校有权保管、并向有关部门送交学位论文的原件与复印件；②学校可以采用影印、缩印或其它复制手段复制并保存论文；③学校可允许论文被查阅或借阅；④学校可以学术交流为目的，复制赠送和交换学位论文；⑤学校可以公布学位论文的全部或部分内容。
本人签名：日期：摘要本设计分析了电力系统短路故障的电气特征，并利用Matlab/Simulink软件对其进行仿真，进一步研究短路故障的特点。
通过算例对电力系统短路故障进行分析计算。
然后运用Matlab/Simulink对算例进行电力系统短路故障仿真，得出仿真结果。
并将电力系统短路故障的分析计算结果与Matlab仿真的分析结果进行比较，从而得出结论。
短路电流计算公式
短路电流计算公式
（1）三相短路电流计算：起始短路电流周期分量有效值/KA
In =1.05Un /√3 /Z∑=1.05Uφ/√R2∑+X2∑
式中Un、U/φ——网络标称电压（线电压、相电压），V。
Z∑、R2∑、X2∑——计算电路总阻抗、总电阻、总电抗，主要为系统、变压器、母线、及线路阻抗，mΩ。
（2）单相接地故障电流及单相短路电流计算由序网分析可知，单相接地故障电流及单相短路电流可由下式求得：
In =3Uφ/|Z1∑+Z2∑+Z0∑|= √3Un/(R1∑+R2∑+R0∑)2+(X1∑+X2∑+X0∑)2
式中：Un、Uφ——网络标称电压（线电压、相电压），V。
Z1∑、Z2∑、Z0∑计算电路正序、负序及零序总阻抗，mΩ。
R1∑、R2∑、R0∑计算电路正序、负序及零序总电阻，mΩ。
X1∑+X2∑+X0∑计算电路正序、负序及零序总电抗，mΩ。
负序阻抗与正序阻抗相等。
零序阻抗为相线零序阻抗与3倍保护线/中性线的零序阻抗之和。
由于配电变压器一般均采用Dyh或Yyh联结，故在计算时无需考虑变压器及高压侧的零序阻抗。
对于用过阻抗接地的TT系统，该阻抗应按3倍计入计算电路的零序阻抗。
可通过计算出相保（相线与保护线PE、PEN)/相零（相线与中性线N)回路阻抗的方法直接求取单相接地故障电流及单相短路电流/KA，此时
In =Un√3/Z = Uφ/√R2+X2
式中Un、Uφ——网络标称电压（线电压、相电压），V。
Z、R、X——相保/相零回路阻抗、电阻、电抗，mΩ。
基于MATLAB的牛顿拉夫逊迭代法计算潮流(附加短路计算)
这个程序可以适用于三机九节点系统（参数见主程序），本来是想编一个通用各种结构的程序的，但是因为鄙人比较懒，老师留作业时候又没说要通用，就没改完。
惭愧啊。
不过大同小异啦。
有兴趣的慢慢改吧。
使用方法：按后文中给出的代码建立.m文件放在一个文件夹里面。
先运行Powerflow_main.m计算算例系统的潮流；然后运行ShortcircuitCalc.m计算算例系统三相短路电流；程序说明详见各.m文件注释部分，写的已经很详细了，慢慢看吧。
Powerflow_main.m文件代码如下：clear%牛顿拉夫逊迭代法计算潮流format short %规定参数数据显示精度%节点参数矩阵%第一列为节点编号%第二列表示有功注入P%第三列表示无功注入Q%第四列表示电压幅值U%第五列表示电压角度θ%第六列表示发电机x′%第七列表示发电机E′%第八列表示节点类型（2表示平衡节点，1表示PV节点，0表示PQ节点）Node_p=[ 1, 0, 0, 1.04 , 0, 0.3, 1.137, 2;2, 1.63, 0, 1.025, 0, 0.3, 1.211, 1;3, 0.85, 0, 1.025, 0, 0.3, 1.047, 1;4, 0, 0, 1.0, 0, 0, 0, 0;5, -1.25, -0.5, 1.0, 0, 0, 0, 0;6, -0.9, -0.3, 1.0, 0, 0, 0, 0;7, 0, 0, 1.0, 0, 0, 0, 0;8, -1, -0.35, 1.0, 0, 0, 0, 0;9, 0, 0, 1.0, 0, 0, 0, 0];count_s=0;countPV=0;for k=1:size(Node_p,1)if Node_p(k,8)==1countPV=countPV+1;else if Node_p(k,8)==2count_s=count_s+1;end;end;end;countPV;count_s;countPQ=size(Node_p,1)-1-countPV;%显示节点参数disp('节点参数如下：')disp(Node_p)%支路参数%第一列为首节点，第二列为末节点，第三列表示R，第四列表示X，第五列表示B/2 %第六列表示支路类型（1为变比为1的变压器元件；2为输电线元件;0为接地支路）Branch_p =[ 1, 4, 0 , 0.0576, 0 , 1;2, 7, 0 , 0.0625, 0 , 1;3, 9, 0 , 0.0586, 0 , 1;4, 5, 0.01 , 0.085 , 0.088 , 2;4, 6, 0.017 , 0.092 , 0.079 , 2;5, 7, 0.032 , 0.161 , 0.153 , 2;6, 9, 0.039 , 0.17 , 0.179 , 2;7, 8, 0.0085, 0.072 , 0.0745, 2;8, 9, 0.0119, 0.1008, 0.1045, 2];%显示支路参数disp('支路参数如下：')disp(Branch_p)%设置节点初值U=Node_p(:,4);e_ang=Node_p(:,5);P=Node_p(:,2);Q=Node_p(:,3);save data.matformat long%计算结果数据显示精度%显示节点导纳矩阵admi();disp('节点导纳矩阵Y');sparseYKmax=10; %设置最大迭代次数kaccuracy=10^-7;%设置迭代精度k=0;%迭代次数初始化为零for k1=1:Kmax[dP,dQ,y]=getY(U,e_ang);if max(abs(y))<accuracybreak;end;J=jacob(U,e_ang,dP,dQ);x=-inv(J)*y;de_ang=[0;x(1:8)];dU=x(9:14);u1=U(2:9)*dU.';u=diag(u1);U(4:9)=U(4:9)+u;e_ang=e_ang+de_ang;k=k+1;end;e_ang=e_ang/pi*180;save result.matdisp('迭代次数：')kdisp('4号节点至9号节点电压幅值如下：') disp(U(4:9))disp('2号节点至9号节点电压相角如下：') disp(e_ang(2:9))admi.m文件代码如下：%节点导纳矩阵的形成format long %规定数据格式NI=size(Branch_p,1);k_t=1;Y=zeros(NI);for m1=1:NI;I=Branch_p(m1,1);J=Branch_p(m1,2);R=Branch_p(m1,3);X=Branch_p(m1,4);b=Branch_p(m1,5);Style=Branch_p(m1,6);if Style==1 %判断为变压器元件Y(I,I)=Y(I,I)+1/(R+1j*X);Y(J,J)=Y(J,J)+1/(R+1j*X)/k_t/k_t;Y(I,J)=Y(I,J)-1/(R+1j*X)/k_t;Y(J,I)=Y(J,I)-1/(R+1j*X)/k_t;else if Style==0 %判断为母线接地支路元件Y(I,J)=Y(I,J)+1/(R+1j*X);else %判断为输电线元件Y(I,I)=Y(I,I)+1j*b+1/(R+1j*X);Y(J,J)=Y(J,J)+1j*b+1/(R+1j*X);Y(I,J)=Y(I,J)-1/(R+1j*X);Y(J,I)=Y(J,I)-1/(R+1j*X);end;end;end;Y;G=real(Y);B=imag(Y);sparseY=sparse(Y);save data.mat;getY.m文件代码如下：（线性方程组常写作AX=Y形式，故此处命名为get_Y）%计算ΔP和ΔQ的函数function [dP,dQ,y]=getY(U,e_ang)load data P Q G B Node_p countPV count_s;dP=zeros(size(Node_p,1),1); sum1=zeros(size(Node_p,1),1);for i=1:size(Node_p,1)for j=1:size(Node_p,1)sum1(i)=sum1(i)+U(j)*(G(i,j)*cos(e_ang(i)-e_ang(j))+B(i,j)*sin(e_ang(i)-e_ang(j)));end;dP(i)=P(i)-U(i)*sum1(i);end;dQ=zeros(size(Node_p,1),1);sum2=zeros(size(Node_p,1),1);for i=1:size(Node_p,1)for j=1:size(Node_p,1)sum2(i)=sum2(i)+U(j)*(G(i,j)*sin(e_ang(i)-e_ang(j))-B(i,j)*cos(e_ang(i)-e_ang(j)));end;dQ(i)=Q(i)-U(i)*sum2(i);end;y=[dP((count_s+1):9);dQ((countPV+count_s+1):9)]; %拼接ΔP和ΔQ构成方程-J*x=y的向量yjacob.m文件的代码如下：%形成雅克比矩阵的函数function J=jacob(U,e_ang,dP,dQ)load data B G Q P Node_p countPV count_s;size_Y=size(Node_p,1);%H矩阵H1=zeros(size_Y);for i=1:size_Yfor j=1:size_Yif j==i;H1(i,j)= U(i)*U(i)*B(i,j)+Q(i)-dQ(i);elseH1(i,j)=-U(i)*U(j)*(G(i,j)*sin(e_ang(i)-e_ang(j))-B(i,j)*cos(e_ang(i)-e_ang(j)));end;endendH=H1(2:size_Y,2:size_Y);%N矩阵N1=zeros(size_Y);for i=1:size_Yfor j=1:size_Yif j==i;N1(i,j)= -U(i)*U(i)*G(i,j)-(P(i)-dP(i));elseN1(i,j)=-U(i)*U(j)*(G(i,j)*cos(e_ang(i)-e_ang(j))+B(i,j)*sin(e_ang(i)-e_ang(j)));end;end;end;N=N1(2:size_Y,(countPV+count_s+1):size_Y);%M矩阵M1=zeros(size_Y);for i=1:size_Yfor j=1:size_Yif j==i;M1(i,j)=U(i)*U(i)*G(i,j)-(P(j)-dP(j));elseM1(i,j)=U(i)*U(j)*(G(i,j)*cos(e_ang(i)-e_ang(j))+B(i,j)*sin(e_ang(i)-e_ang(j)));end;end;end;M=M1((countPV+count_s+1):size_Y,2:size_Y);%L矩阵L1=zeros(size_Y);for i=1:size_Yfor j=1:size_Yif j==i;L1(i,j)= U(i)*U(i)*B(i,j)-(Q(i)-dQ(i));elseL1(i,j)= -U(i)*U(j)*(G(i,j)*sin(e_ang(i)-e_ang(j))-B(i,j)*cos(e_ang(i)-e_ang(j)));end;endendL=L1((countPV+count_s+1):size_Y,(countPV+count_s+1):size_Y);J=[H N;M L];%拼接构成雅克比矩阵ShortcircuitCalc.m文件的代码如下：clear%三相短路计算format long;%对YN进行修正，形成包括发电机内阻抗和负荷阻抗的节点导纳矩阵re_admi();Z=inv(rY);%计算节点阻抗矩阵disp('4节点发生金属短路')f=4;%短路点为4节点%输出节点阻抗矩阵的短路点所在列disp('节点阻抗矩阵的短路点所在列Z(:,f)=');Z(:,f)%短路电流If计算load result U e_ang;e_angf=e_ang(f);Uf=U(f)*(cos(e_angf/180*pi)+1j*sin(e_angf/180*pi));Zff=Z(f,f);zf=0;If=Uf/(Zff+zf);i_angf=angle(If)*180/pi;If=abs(If);%短路电流计算结果显示disp('短路电流幅值:')Ifdisp('短路电流相角(单位为°)')i_angf%短路时各节点电压计算U_k=U;for x1=1:size(Node_p,1)U_k(x1)=U(x1)*(cos(e_ang(x1)*pi/180)+1j*sin(e_ang(x1)*pi/180))-Z(x1,f)*If*(cos(i_angf*pi/180)+ 1j*sin(i_angf*pi/180));end;%Uk为短路时各节点电压幅值Uk=abs(U_k);%uk_ang为短路时各节点电压相角（单位为°）uk_ang=angle(U_k)*180/pi;%输出短路时各节点电压disp('短路时1-9节点电压：')disp('幅值：')Ukdisp('相角：(单位为°)')uk_ang%计算短路时各支路电流%输出矩阵初始化%第一列为首节点i%第二列为末节点j%第三列为Ii幅值%第四列为Ii相角(单位为°)%第五列为Ij幅值%第六列为Ij相角（单位为°）Ik=[1 4 1 0 1 0;2 7 1 0 1 0;3 9 1 0 1 0;4 5 1 0 1 0;4 6 1 0 1 0;5 7 1 0 1 0;6 9 1 0 1 0;7 8 1 0 1 0;8 9 1 0 1 0];load data Branch_p k_t;%其中变压器变比为k_tNI=size(Branch_p,1);for m1=1:NI;I=Branch_p(m1,1);J=Branch_p(m1,2);R=Branch_p(m1,3);X=Branch_p(m1,4);b=Branch_p(m1,5);Style=Branch_p(m1,6);if Style~=1 %判断为非变压器支路Ik(m1,3)=(U_k(I)-U_k(J))/(R+1j*X)+U_k(I)*1j*b;Ik(m1,4)=angle(Ik(m1,3))*180/pi;Ik(m1,3)=abs(Ik(m1,3));Ik(m1,5)=(U_k(J)-U_k(I))/(R+1j*X)+U_k(J)*1j*b;Ik(m1,6)=angle(Ik(m1,5))*180/pi;Ik(m1,5)=abs(Ik(m1,5));else %否则为变压器支路Ik(m1,3)=(U_k(I)-U_k(J))/(R+1j*X)/k_t;Ik(m1,4)=angle(Ik(m1,3))*180/pi;Ik(m1,3)=abs(Ik(m1,3));Ik(m1,5)=(U_k(J)-U_k(I))/(R+1j*X)/k_t;Ik(m1,6)=angle(Ik(m1,5))*180/pi;Ik(m1,5)=abs(Ik(m1,5));endend%输出短路时各节点电流disp('第一列为首节点i')disp('第二列为末节点j')disp('第三列为Ii幅值')disp('第四列为Ii相角(单位为°)')disp('第五列为Ij幅值')disp('第六列为Ij相角（单位为°）')format shortIksave result_k%形成包括发电机内阻抗、负荷阻抗的节点导纳矩阵format long %规定数据显示精度load data P Q Y Node_p ;load result U;rY=Y;for n1=1:size(Node_p,1)Xd=Node_p(n1,6);Ed=Node_p(n1,7);style=Node_p(n1,8);if style~=0 %判断是否为发电机节点rY(n1,n1)=rY(n1,n1)+1/(1j*Xd);%YN中与发电机节点对应的对角线元素增加发电机导纳else %其他节点(包含负荷节点）rY(n1,n1)=rY(n1,n1)+(-P(n1)+1j*Q(n1))/(U(n1)*U(n1));%与负荷节点对应的对角线元素增加负荷导纳%对于非负荷节点Y矩阵元素不做修改但仍满足上式end;end;rY;。
基于MATLAB的短路电流计算程序编制
．
１ ＭＡ Ｌ Ｂ 介 ＴＡ 简
Ｍ Ｔ Ａ 是 英 文 Ｍ ｔｘＬｂｒｔｙ ＡＬＢ ａ ｉ ａｏ ｏ 的缩 写 ， ｒ ａｒ 自 １ ８ 年 美 国Ｍａ Ｗａｋ 司 推 向市 场 以来 ， 到 广泛 ９４ ｔ ｒ公 ｈ 得
Ｖ ’ ｆ 一 ， （） ３
式 中： 是 故障节点ｆ 的自电抗 。
这 ３ 方程 中包 含 了６ 个 个未 知量 。因此需 要根 据 不对 称 短路 的具 体 边 界条 件 写 出另￣ ３ 边界 条 件 ｂ个 才能求 解 ｌ ３ ｌ 。 ２２１ 单相 （相 ） ．． ａ 接地 短路 单相 接地 短路 的边界 条 件为
假 定 系统 中的节点 经 过 渡 阻抗 发 生短 路 。 这 个 过 渡 阻抗 不 参 与形 成 网络 的节 点 阻抗 矩 阵 。保
编程技术与数据库技术较好地统一起来 。本文使用 Ｍ Ｔ Ａ 软件实现短路电流计算建模及程序编制 。 ＡＬＢ
一 ０舞 ◇～ ◇
持故障处 的边界条件不变 ， 把网络的原有部分同故 障支路分开 。 容易得 出 ， 对于正常状态的网络而言 ， 发生短路相 当于在故 障节点蹭 加 了一个注入 电流
便地调用有关Ｆ ｒａ 、语言的子程序 。 ｏｒ Ｃ ｔｎ 另外 ， 它不需 定义数组的维数 ， 并给出矩阵函数 、 特殊矩阵专门的 库 函数 ， 使之在求解诸如信号处理 、 建模 、 系统辨识 、 控 制 、 化 等领 域 的问题 时更加 简捷 、 效 和方便 。 优 高
分 析 、 真 发展 的必然 趋 势 。同时 仿 随着 数据 库技 术
维普资讯
基 于
《短 路 电流 计 算 震 序编 制 》 《 》
matlab直流短路电流计算仿真
在MATLAB中进行直流短路电流计算的仿真，可以按照以下步骤进行：
1. 导入系统参数：使用“psat”函数将系统参数导入系统模型中。
这些参数通常包括电源电压、系统阻抗、短路位置和短路类型等。
2. 创建系统模型：使用MATLAB中的Simulink模块库，根据需要创建一个新的系统模型。
可以选择无穷大功率电源供电系统或有限大功率电源供电系统。
3. 计算等值阻抗：使用“psadeqz”函数计算等值阻抗，并保存结果。
等值阻抗是计算短路电流的关键参数，它反映了系统的阻抗特性。
4. 计算短路电流：根据上述公式计算短路电流。
通常，短路电流的计算公式是基于欧姆定律和基尔霍夫定律的。
5. 进行仿真模拟：在Simulink环境中，设置仿真时间、仿真算法等参数，然后运行仿真模拟。
观察仿真结果，包括短路电流的波形、幅值和持续时间等。
6. 分析仿真结果：对仿真结果进行深入分析，包括计算电流的周期性变化、冲击电流的大小等。
这些数据可以为实际系统的设计和优化提供参考。
需要注意的是，在进行直流短路电流计算的仿真时，需要考虑系统的实际情况和各种不确定性因素。
例如，电源电压的波动、负载的变化和线路阻抗等因素都可能对仿真结果产生影响。
因此，在进行仿
真时需要进行合理的假设和简化，同时对仿真结果进行谨慎的分析和解释。
三相短路电流计算
短路电流次暂态值 I ′′ = Id = 9.16 = 2.92kA X*KΣ 3.14
短路电流冲击值 ish = 2K sh I ′′ = 2.55 × 2.92 = 7.45kA
次暂态短路功率 S ′′ = 3U av I ′′ = 3 × 6.3 × 2.92 = 31.86MVA
4.4.2 有限容量电源系统的三相短路电流计算
1．有效值的计算---运算曲线法 有限容量电源系统发生三相短路后，其母线电压不再保持恒定，短路电流周期分量也随 之发生变化。如果已知短路后某一时刻发电机的电势，则短路电流周期分量相应时刻的有效 值可按下式求取
I pt =
Et 3X KΣ
（4.4.9）
式中： Et ——短路后 t 时刻的发电机的电势；
（4.4.4）
式中： I d ——计算点所在电压级的基准电流， I d =
Sd 3U d
X *KΣ ——短路回路总电抗的标幺值， X *KΣ = X KΣ
3I d Ud
于是可得：
2.短路电流冲击值
I ′′ =I ∞=
Ip
=
Id X *KΣ
（4.4.5）
由式（4.2.9）和（4.2.10）可求得短路电流冲击值和短路冲击电流有效值
X KΣ ——短路回路总电抗。
但是同步发电机突然短路时，电势随时间变化的规律是很复杂的，用上式计算比较困难。 电力部门根据国产同步发电机参数和容量配置情况，用概率统计的方法分别制定了汽轮发电 机和水轮发电机的短路电流运算曲线。利用运算曲线可以方便地查出三相短路电流周期分量
的有效值 I*pt ，因此在实际工程计算中，通常采用“运算曲线”来求解三相短路电流周期分
Id1 =
Sd = 3U av
电力系统短路故障分析的MATLAB辅助程序设计短路计算程序
电力系统短路故障分析的MATLAB辅助程序设计短路计算程序电力系统短路故障分析是电力系统设计和运行过程中非常重要的一环。
短路故障会导致电力系统各个部分的电压、电流和功率的突然变化，对设备的保护和稳定运行产生不利影响甚至引起事故。
因此，进行短路计算和故障分析非常必要。
MATLAB是一种功能强大的数值计算和数据可视化工具，对于电力系统短路计算和故障分析也可以发挥重要的作用。
下面将介绍如何使用MATLAB设计一个简单的电力系统短路计算程序。
首先，我们需要建立一个电力系统的模型。
电力系统可以用图模型表示，其中节点表示发电机、变压器、负荷等设备，边表示导线、变压器等电力连接。
我们可以使用MATLAB中的图模型工具箱创建电力系统模型，并且设置各个节点和边的属性，例如电压、电流、阻抗等。
然后，我们需要编写短路计算程序。
短路计算可以分为对称故障和不对称故障两种情况。
对称故障是指短路故障发生在电力系统的正常运行条件下，例如三相短路。
不对称故障是指短路故障发生在电力系统的不正常运行条件下，例如单相接地短路。
对于对称故障，我们可以使用节点电流法进行计算。
首先，应用基尔霍夫电流定律，根据电压和阻抗计算电流。
然后，根据节点电流方程和电流方程计算电流分布。
最后，根据电流分布计算短路电流和故障点的电压。
对于不对称故障，我们可以使用仿真方法进行计算。
首先，需要设置故障位置和故障类型，例如A相到地短路。
然后，根据故障位置和类型修改节点和边的参数，例如将故障位置的阻抗设置为零。
最后，使用数值方法求解电力系统的动态响应，得到短路电流和故障点的电压。
在MATLAB中，可以使用矩阵运算和数值求解函数实现短路计算。
例如，可以使用矩阵乘法和矩阵求逆函数计算节点电流和电流分布。
可以使用ODE求解器求解动态响应方程。
可以使用MATLAB的绘图函数绘制电力系统的电流分布和故障点的电压。
总结起来，电力系统短路故障分析的MATLAB辅助程序设计涉及建立电力系统模型、编写短路计算程序并使用MATLAB的数值计算和数据可视化工具进行计算和分析。
三相短路电流的计算
技
业的贡献，有必要从专业、技能、经验、学历以及年龄等多方面进行考
术
量。
方
面
➢ 软件技术领域最好根据研发人员所擅长的技术能力来划分不同的专业技
术组。
➢ 培养研发人员对新技术的学习掌握，并创新应用。
5.1.2 软件生产函数
提高软件生产率的主要举措有：
➢ 建立良好的人才激励机制，通过物质、精神奖励以及职务晋升制度提高
X *ca
X *dK
S N Sd
（4.4.11）
南京师范大学电气工程系
有限容量电源系统三相短路电流的计算
注意：
通过运算曲线法求得的I*pt 是以向短路点直接提供短路 电流的发电机总容量为基准功率所对应的三相短路电流 周期分量在t时刻的有效值的标幺值。因此， 所求短路 后t时刻的短路电流周期分量有效值为
③生产函数具有非负性，总产出必 须是正值，且总产量是生产要素组 合的结果，单一要素的投入是不能 获得产出的
5.1.1 生产函数
3. 规模报酬 规模报酬又称为规模经济，主要研究所有要素都同比例增加时，产量时如何发生变化的。
通常在规模较小时，增加投入可以促进规模报酬递增，而当规模达到一定 程度时，则会出现规模报酬递减的情况。
Ip
U av 3X K
无穷大电源系统三相短路电流的计算
2.短路电流冲击值
由式（4.2.9）和（4.2.10）可求得短路电流冲 击值和短路冲击电流有效值
ish 2Ksh I I sh I 1 2(Ksh 1)2
（4.4.6） （4.4.7）
无穷大电源系统三相短路电流的计算
3.三相短路功率
4.4 三相短路电流的计算
无穷大容量电源系统的三相短路电流计算 有限容量电源系统的三相短路电流计算
matlab基于节点阻抗矩阵的三相短路计算
matlab基于节点阻抗矩阵的三相短路计算MATLAB基于节点阻抗矩阵的三相短路计算三相短路是电力系统中最常见的故障类型之一，也是最严重的一种故障，其产生的电流会对设备造成故障、损坏电力设备，甚至会导致火灾等事故。
因此，对电力系统进行三相短路计算及分析非常必要，MATLAB是一款通用的工具软件，可用于电力系统的短路计算中，本文就基于节点阻抗矩阵介绍MATLAB的三相短路计算。
一、节点阻抗矩阵节点阻抗矩阵是一种直观、简单、易于理解的方法，其基本思想是将电力系统中每个节点的短路电流计算单独列成一个向量，向量中每个元素都代表着该节点与其他节点之间的电流响应系数。
节点阻抗矩阵根据电力系统的拓扑结构所形成，元素值来自于两个节点之间的阻抗和电导之和。
二、节点阻抗矩阵的计算1、定义节点位置和电力设备参数在MATLAB中，首先需要定义电力系统中所有节点的位置参数（x，y）、所有支路的编号、阻抗参数值等，定义方法如下：node_position=[1 1; 1 1.5; 1.5 1; 1.5 1.5]; %四个节点的位置line=[1 2 0.02 0.04; 1 3 0.01 0.03; 2 4 0.03 0.06; 3 40.02 0.04];%四条线路的起始节点、结束节点、电阻、电抗2、构建节点阻抗矩阵根据节点位置和电力设备参数，可以通过以下语句构建节点阻抗矩阵：[~,Y_node]=Admittance_Matrix(line,4); %采用自定义函数计算导纳矩阵Z_node=inv(Y_node);%计算节点阻抗矩阵其中，Admittance_Matrix函数是一个自定义函数，用于求取系统的导纳矩阵，在导纳矩阵中，对称元素等于发电机到负载电工的导纳，非对称元素等于接线点的线路阻抗之和，函数的具体实现方法可以查阅MATLAB帮助文档。
三、三相短路计算有了节点阻抗矩阵，就可以进行三相短路计算，MATLAB中可以通过以下步骤进行计算：1、定义故障考虑的节点编号和负荷首先需要确定故障考虑的节点，即需要计算的节点，有多少个节点，就需要计算多少次，在节点矩阵中的位置，判定方法如下：fault_node=2;%节点2作为故障节点bus=[1 0; 2 10+4*j; 3 10+j; 4 15];%分别代表总发电机、故障节点、非故障节点1和非故障节点2的编号和电阻2、计算故障前的稳态电压在进行三相短路计算之前，需要先计算故障前的稳态电压值，详细计算方式可以参考MATLAB帮助文档或其他相关资料，计算发电机电动势和短路电流，具体计算方法如下：[num_bus,~]=size(bus);V_node=zeros(num_bus,num_bus);for i=1:num_busfor j=1:num_busV_node(i,j)=(bus(i,2)-bus(j,2))/Z_node(i,j);endend%计算电动势和短路电流Es=bus(1,2)-V_node(1,2)*Z_node(1,2);Fault_Pre_Curr=(Es-bus(fault_node,2))/Z_node(1,fault_node);3、计算短路电流和故障后电流得到故障前稳态电压后，可以根据一定的公式计算短路电流和故障后电流，具体计算方法如下：Z_fault=0.03+0.02*j;%故障阻抗I_fault=bus(fault_node,2)/(Z_node(fault_node,fault_node) +Z_fault);I_fault_phase=I_fault/(3^0.5);4、计算故障后电压最后，可以根据故障后电流计算出故障时间的电压，公式如下：V_fault_node=bus(fault_node,2)-I_fault*Z_node(fault_node,fault_node);总结本文简要介绍了MATLAB基于节点阻抗矩阵的三相短路计算方法，通过以上的步骤，可以较为准确地计算出电力系统中三相短路的电流和电压，为电力系统的安全稳定运行提供了重要保障。
基于Matlab的电力系统短路电流分析与仿真
ｓ ｔ ｙｓｅｍ ．
Ｋｅ ｒ ｓ Ｍ ａ ｌ ｂ Ｓ ｍｕ ａ ｉ ｎ； ｉ ｐ ｗｅ Ｓ ｓ ｅ ｓ Ｓ ｏ ｔｃ ｒ ｕ ｔｃ ｒ ｅ ｔ ａ ｃ ｌ ｔｏ ｙ ｗｏ ｄ ： ｔａ ； ｉ ｌ ｔｏ Ｓ ｍ ｏ ｒ ｙ ｔ ｍ ： ｈ ｒ ． ｉｃ ｉ ｕ ｒ ｎ ｌ ｕ ａ ｉ ｎ ｃ
ｌ 引言
值 或 百 分值 给 出 的 ，而 在 进 行 电力 系 统 计 算 时 ，必 须取 统 一 的 基 准 值 ， 所 以要 将 原 来 的 以本 身 额 定 值 为基 准 值 的 阻抗 标 么 值 换 算 到 统 一 的基 准 值 。 一般 先选 定 基 准 容 量Ｓ和 基 准 电 上 ｕ。 而在 工 程 设 计 中通 常 取 Ｓ ＩＯ Ｖ・ ，Ｕ＝ Ｋ。 Ｏ Ｍ Ａ ｕ，其 中 ｕ为 ： 短 路 点 的短 路 计 算 电压 。 冈此 ：
十 分 严 重 的 ， 有 必 要 进 行 短 路 计 算 。而 短 路 电 流计 算 是 一
项复杂而又责任重大 的工作 ，这些 工作历来 南技术人员手
工 计 算 完 成 ，误 差 比较 难 免 ， 因 而一 个 完 整 的 工 作 必 需 要
基准电流为：
基准电抗为：
瓦ｄ Ｓ ： Ｓ ｄ
Ｓ ｏ ｔｃｒ ｕ ｔＣｕ ｒ ｎｔ ａｙ ｉ ｎｄＳｉ ｕｌｔｏ ｓ ｄ ｏ ａ ｌｂ Ｐｒ ｇ ａ ｓｔ ｈ ｅ ｅＰ ｗｅ ｙ ｔｍ ｆＳ ｏ ｔｃｒ ｕ ｔ ｈ ｒ ｉｃ ｉ ｒ ｅ Ａｎ ｌｓｓａ ｍ ａ ｉ ｎＢａ ｅ ｎ Ｍ ｔａ ｏ ｒ ｍ ｏ Ａｃ ｉｖ ｏ ｒＳ ｓｅ ｏ ｈ ｒ —ｉｃ ｉ Ｃｕ ｒ ｎ ａ ｙ ｉ ｎｄ Ｓ ｍ ｕ ａ ｉ ｎ ｒｅ ｔ ｎ Ａ ｌ ｓｓ ａ ｉ ｌ ｔｏ
基于MatlabSimulink的三段电流继电保护在教学及实验中的应用
基于Matlab/Si mul i nk 的三段电流继电保护在教学及实验中的应用鄢Appl i c at i on of Three —stage Current Relay Protecti o& Based on Matlab/Si mul i &k in Teach i &g and Exper i me&t 王海军（南京铁道职业技术学院，江苏南京210031）张紫莹（中国铁路上海局集团有限公司南京供电段，江苏南京210015）摘要：电流保护是继保专业的教学基础,一般都采用实验手段对其进行验证。
为了提高实验手段和效率，通过实例进行 三段式电流保护的计算，利用MATLAB/Si mul i nk 建立对应的仿真模型，并进行相应的动作仿真和线路故障分析，验证了三段式电流保护的行为特性和模型的正确性,给其他实际工况对象奠定仿真基础遥关键词：电流保护；Matlab/Si m ul i n k ；保护模型；仿真分析;教学实验Abstract :Current protect i ons the teach i n g bas i s of relay protect ion,and it is generally ver i f i e d by exper i ment.ln order to improve the e xper i m ental methods and eff i c i e ncy,the art i c le uses examples of three-stage current protect i o n to ver i f y,thecorrespond i n g s i m ulat i o n model is establ i s hed w i t h MATLAB/Si m ul ink .I t is correct that The character i s t i c s and model of the three-stage current were ver i f i e d,and s imulat i o n foundat i o n of the other s i m ulat i on is establ i s hed 。
Matlab-Three-PhaseSource（三相电源）
Matlab-Three-PhaseSource（三相电源）内部R-L 感性阻抗的三相电源的实现库电源模块库电源描述三相电源模块执⾏带有内部R L 阻抗的平衡的三相电压来源。
三电压来源在中与中性点连接成Y 或者其他⽅式。
你可以通过规定电源短路容量或者是直接规定电源内部的电阻和电感R 和L 或者是间接地设置X 与R ⽐率。
对话框和参数参数表线电压有效值内部线电压单位：VA 相初相⾓指定内部电压的相⾓为A 相相⾓，⽤度数表⽰。
三相电压的相序为正序。
这样，B 相和C 相的电压分别落后A 相120°和240°。
频率三相电源的频率单位：Hz内部连接⽅式三相电源的内部连接。
改变电源连接，电源模块图标就跟着改变。
选择其中⼀个连接⽅式：通过短路电流容量指定阻抗选择它表明你想通过短路电流容量和X/R 的⽐值来指定阻抗值。
基准电压下三相短路电流容量基准电压下三相电感短路功率单位：VA ，⽤来计算内部电感值。
只有选中通过短路电流容量指定阻抗选项时，才能够使⽤该参数。
在给定基准相电压V base （单位：V ）和电感三相短路功率Psc （单位：VA ）时，内部电感L （单位：H ）可由下式得到：212base sc V L P fπ=? 基准电压相基准电压，单位：V ，⽤来指定三相短路电流容量。
基准电压值通常就是名义电压源的⼤⼩。
只有选中通过短路电流容量指定阻抗选项时，才能够使⽤该参数。
X/R ⽐值在标称电源频率或者X/R ⽐值内部电源阻抗品质因数下的X/R ⽐值。
只有选中通过短路电流容量指定阻抗选项时，才能够使⽤该参数。
在给定内部电阻R （单位：Ω）和电感三相短路功率X （单位：Ω）时，X/R ⽐值可由下式得到：2(/)/X fL R X R X Rπ==只有不选中通过短路电流容量指定阻抗选项时，才能够使⽤该参数。
电源电阻单位：Ω电感只有不选中通过短路电流容量指定阻抗选项时，才能够使⽤该参数。
电感单位：H潮流计算表表中的这些参数是 Powergui模块潮流计算的⼯具，这些潮流计算参数只能给模型进⾏初始化，对块模型和仿真参数没有影响。
基于MATLAB的短路电流计算程序编制
言， 而且 几乎 与数学 表达式 相 同 ， 言 中基 本 元素是 语
节 点 电压 的故 障分 量 。两者 的叠 加就 等于发 生短 路 后 节点 ｉ 的实 际 电压 ， ： 即 Ｖ ＝Ｖ ’ Ｉｆ ｉ ｉ 一Ｚｆ Ｉ （） ２
序编制 。
关 键 词 短路 电流计算
ＭＡ Ｌ Ｂ 短路故 障 ＴＡ
在 电力 系统 的设计 、 行 的各个 部 门 ， 运 无论 是在 设 计时 需要选 择 断 路器 等 电气设 备 ， 还是 在 电 力 系 统 运行 时要 进行 继 电保 护 整 定 值计 算 ， 要 事先 进 都
以直接 调用 。 同时通 过建 立 Ｍｅ调 文件 的形式 可 以 方便地 调 用有 关用 Ｆ ｒａ 、 言的子程 序 。另 外 ， ｏ ｒｎ Ｃ语 ｔ 它不 需定 义数 组 的维数 ， 给 出矩阵 函数 、 殊矩 阵 并 特
统分 析 、 仿真发 展 的必然 趋势 。 同时 ， 随着数 据库 技
术 的发展 ， 数据 的管理更 加有 序 、 便 。这 为 电力应 方
保持 故障处 的边界条件不 变 ， 网络 的原有部 分 同故 把 障支 路分开 。容得 得 出， 于正 常状 态 的网 络而 言 ， 对 发生 短路相 当 于在 故 障节 点 ｆ 加 了一 个 注入 电流 增
（） ２ 式也 适用 于故 障节 点 ｆ于是 有 ，
Ｖ ＝Ｖ 一ｚ Ｉ ｆ ｆ’ ｆ （） ３
矩阵， 可提供各种矩阵的运算和操作 ， 且具有符号计 算 、 学 和文字 统一 处理 ， 数 离线 和在 线计算 等 功能 。 与传 统 的 Ｂ ｓ 、 ａｉ Ｃ语 言 及 Ｆ ｒａ ｃ ｏｔ ｎ语 言 相 比 ， ｒ Ｍ ＴＡ Ａ Ｔ ＡＢ简 介
基于MATLAB短路电流的计算与仿真之欧阳术创编
摘要短路故障也称为横向故障，是破坏电力系统正常安全稳定运行的主要原因。
一般分为四种常见的短路情况，分别是三相短路，两相短路、单相接地短路及两相接地短路。
利用传统的数学计算分析电力系统故障情况不但工作量大而且危险性高。
课题中利用了MATLAB软件对电力系统短路故障进行仿真模拟，这样既可以不偏离实际情况又可以减少实验的危险性。
首先建立一个空白的模型，选择Simulink元件库中本身具有的电力系统模块库的模块，分别建立无穷大功率电源供电系统系统和有限大功率电源供电系统，得出四种短路情况分别发生时的图形，然后对两种不同供电系统发生三相短路的进行数学计算得出周期电流的幅值大小和冲击电流的大小。
通过相关的系统仿真与计算，结果表明利用MATLAB软件对电力系统进行仿真，方便快捷，同时也具有很高的可靠性。
关键词：短路故障；MATLAB；电力系统ABSTRACTShort circuit, also known as lateral fault, is the main reason for the destruction of normal power system safe and stable operation. Generally divided into four common short-circuit conditions, which are three-phase short, two-phase short circuit, single-phase ground short circuit and two ground fault. Calculation of power system failure not only a heavy workload and high-risk use of traditional mathematics. The use of MATLAB software subject to short circuit power system simulation, so that both can not deviate from the actual situation and can reduce the risk of the experiment. First create a blank model, select Simulink component library itself has a power system module library modules, namely the establishment of the infinite power supply system power supply systems and limited system, draw graphics were four short-circuit occurs, Then two different phase short-circuit of the power supply system mathematically calculated the amplitude of the current cycle and theimpact of the current size. Related systems through simulation and calculation results show that the use of MATLAB software for power system simulation, convenient, but also has high reliability.Keywordsshort-circuit fault; MATLAB; power system目录1 绪论11.1 研究的背景与意义11.2 设计的主要工作12 电力系统短路故障简介32.1 短路的基本概念32.2 短路的危害42.3 短路计算的目的53 仿真软件63.1 MATLAB的简介63.2 Simulink在仿真中的应用74 无穷大功率电源短路故障计算与仿真104.1 三相短路104.2 两相接地短路194.3 单相短路接地214.4 两相短路235 有限大功率电源短路故障计算与仿真265.1 三相短路故障265.2 两相接地短路315.3 单相短路接地325.4 两相短路336 结论35参考文献36致谢391 绪论1.1 研究的背景与意义随着社会的不断发展，电力已经渗入到了人们的学习、生活等各个方面，只有保证不间断地供电才能使工厂正常生产、人们正常生活。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
B node3 X2
H 1 0.27758
Y(node2(aa),node2(aa))=Y(node2(aa),node2(aa))+1/(R(aa)+j*X1(aa));
end
[m3,m4]=size(node3);
for aa=1:m4
Y(node3(aa),node3(aa))= Y(node3(aa),node3(aa))+1/(j*X2(aa));
H 5 0.1041
end
fclose(fid);
n=max(max(node1),max(node2)); %求取结点个数
%% 形成导纳矩阵（初始化）
Y=zeros(n,n);
[m1,m2]=size(node1);
fprintf('程序计算的冲击电流Im: %10.4fKA , 手算结果： 225.7195KA \n\n',1.85*sqrt(2)*abs(1.0664/Z(4,4))*100/sqrt(3)/10.5);
下面是程序读入的txt文件，自己新建一个即可(格式模仿bpa格式)
L 1 2 0 0.0328
L 2 3 0 0.0189
L 3 4 0 0.1167
%% 三门水电站三相短路程序
%% 电气一班 刘永康 0805010114
%% 短路数据读入部分（读入数据为自定应格式）
clear
clc
format long
fid=fopen('C:\Users\john\Desktop\三门水电站程序\三门水电站设计短路部分\短路读入数据.txt','r');
L 3 5 0 0.0583
L 4 6 0 0.075
L 5 6 0 0.0375
fprintf('程序计算的冲击电流Im: %10.4fKA ， 手算结果： 15.7184KA \n\n',1.85*sqrt(2)*abs(1.0118/Z(3,3))*100/sqrt(3)/230);
disp('中压侧：');
fprintf('程序计算的短路电流I": %10.4fKA ， 手算结果： 7.9154KA\n',abs(1.0610/Z(6,6))*100/sqrt(3)/115);
end
while tline(1)=='H'
node3=[node3 eval(tline(5))];
X2=[X2 eval(tline(12:18))];
e=fgetl(fid);
for aa=1:m2
Y(node1(aa),node2(aa))=Y(node1(aa),node2(aa))-1/(R(aa)+j*X1(aa));
Y(node2(aa),node1(aa))=Y(node1(aa),node2(aa));
Y(node1(aa),node1(aa))=Y(node1(aa),node1(aa))+1/(R(aa)+j*X1(aa));
end
%% 计算三相短路电流并输出
%% 有潮流计算结果可知：高压侧电压为1.0118，中压侧电压为1.0610，低压侧电压为1.0664
Z=inv(Y);
disp('高压侧：');
fprintf('程序计算的短路电流I": %10.4fKA ， 手算结果： 6.0079KA\n',abs(1.0118/Z(3,3))*100/sqrt(3)/230);
fprintf('程序计算的冲击电流Im: %10.4fKA , 手算结果： 20.709KA \n\n',1.85*sqrt(2)*abs(1.0610/Z(6,6))*100/sqrt(3)/115);
disp('低压侧：');
fprintf('程序计算的短路电流I": %10.4fKA ， 手算结果： 86.2745KA\n',abs(1.0664/Z(4,4))*100/sqrt(3)/10.5);
node2=[node2 eval(tline(13))];
R=[R eval(tline(19))];
X1=[X1 eval(tline(25:30))];
tline=fgetl(fid);
end
while tline(1)~='H'
tline=fgetl(fid);
H 1 0.4496
H 2 0.0756
H 4 0.20825
1234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890
A node1 node2 R X1
tline=fgetl(fid);
node1=[];node2=[];node3=[];R=[];X1=[]; X2=[]; %初始化，或者后面加个变量a，用node1（a）,就不用初始定义了
while tline(1)~='L'
tline=fgetl(fid);
end
while tline(1)=='L'
node1=[node1 eval(tline(5))]; %平衡结点为1，PV结点为2，PQ结点为3