如何用matlab实现无限大功率电源供电系统的三相短路分析

暂态上机计算教学指导书（适用于Matlab）
一、基础知识
1． 为完成本次设计，需要掌握的Matlab的基础知识有：数组的创建和计算，循环语句，条件选择语句，条件判断语句，绘制二维图的语句；
2． 提供了Matlab的学习参考《精通Matlab6.5》光盘版；
3． Matlab是区分大小写的；
4． 本次设计需要学习和掌握的专业知识：《电力系统暂态分析》相关知识；
5． matlab 在三角函数计算时默认是弧度，而不是度；
6． Matlab 数组（矩阵）的下标从1开始；
7． Matlab 默认小数点后保留4位；
8． 如果不想显示变量的数值，则在该行末尾加“;”符号。

三、程序设计指导
提供了两个程序，（1）不对称短路，名字是buduichen.m；（2）稳定计算的程序，名字是wending.m。所有程序里“%”表示该行为注释行。
1． 不对称短路程序

function buduichen
%数据来源于教材《电力系统暂态分析》P129例（5-1）
%给学生提供的不对称短路故障的部分程序。
%要求的工作（1）自动形成导纳阵；（2）完善短路类型，两相相间短路和两相短路接地；（3）求短路后发电机的电压
clear %清除所有变量的值

%输入支路各序阻抗，ZZ1表示支路正序阻抗，ZZ2表示支路负序阻抗，ZZ0表示支路零序阻抗
ZZ1(1,2)=j*0.1;
ZZ1(1,3)=j*0.1;
ZZ1(2,3)=j*0.1;
ZZ2(1,2)=j*0.1;
ZZ2(2,3)=j*0.1;
ZZ2(1,3)=j*0.1;
ZZ0(1,2)=j*0.2;
ZZ0(2,3)=j*0.2;
ZZ0(1,3)=j*0.2;
%输入节点导纳矩阵，Y1表示不计发电机节点的正序网络节点导纳阵，Y2表示不计发电机节点的负序网络节点导纳阵，Y0表示不计发电机节点的零序网络节点导纳阵，
Y1=[-j*26.6266 j*10 j*10