昆明理工大学电力工程学院学生实验报告

昆明理工大学电力工程学院学生实验报告实验课程名称：《电力系统计算机辅助分析》
开课实验室：计算机中心408 2014年 12月 9日
实验一MATLAB软件的基本操作及程序的基本结构
实验内容：
1.验证实验
2.对于如下电路图，已知
()200cos(90)o
s
u t tω
=+
V，1
22
R=Ω
，
50
L mH
=，233
R=Ω
，
100
C F
μ
=
，自己编写程序，画出
()
c
u t
和
()
s
u t
的波形（四个周期）
1
R
2
R C
L
()
s
u t()
c
u t
+
-
+
-
程序：clear
mUsm=200;
aUsm=(pi/180)*(+90);
Usm=mUsm*exp(j*aUsm);
R1=22,L=0.05,R2=33,C=0.0001;
w=100*pi;
Z1=R1+(w*0.05*j);
Z2=(1/(1/R2+j*w*C));
2、将图2.1中的电压源分别改为2V、4V、6V、8V时，用因子表进行求解，体会该算法带来的益处；
3上机验证课本例17.6－17.7的计算结果
①例17.6
clear
clc
close
n=5;nl=5;
B=[1 2 0.03i 0 1.05 0;
2 3 0.08+0.3i 0.5i 1 0;
3 4 0.015i 0 1.05 1;
2 5 0.1+0.35i 0 1 0;
3 5 0.04+0.25i 0.5i 1 0];
X=[1 0;2 0;3 0;4 0;5 0];
Y=zeros(n);
for i=1:n
if X(i,2)~=0;
p=X(i,1); Y(p,p)=1./X(i,2);
end
end
for i=1:nl
if B(i,6)==0
p=B(i,1); q=B(i,2);
else
p=B(i,2); q=B(i,1);
end
Y(p,q)=Y(p,q)-1./(B(i,3)*B(i,5));
Y(q,p)=Y(p,q);
Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)^2)+B(i,4)./2; Y(p,p)=Y(p,p)+1./B(i,3)+B(i,4)./2;
end
disp('导纳矩阵Y=:');disp(Y)
for k=1:m+1
Z(m+1,k)=0;
end
end
end
end
disp('阻抗矩阵Z=');
disp(Z)
4.设计系统变更时修改导纳矩阵的函数（参考课本P.323），并调试通过，给出导纳矩阵修正函数的调用规则及其主要的实现代码。

n=4;nl=4;B=[1 2 0.08+0.4i 0 1 0;2 3 0.1+0.4i 0 1 0;
3 4 0.3i 0 1 0;1 3 0.12+0.5i 0 1 0];
X=[1 0;2 0;3 0;4 0]
Y=zeros(n);
for i=1:n
if X(i,2)~=0;
p=X(i,1);
otherwise
break
end
思考题：
1.计算短路电流的目的是什么？
答：选择电气设备,计算各种保护的定值,制定运行方案等
2.在计算短路电流时，是否用到了节点阻抗阵中的所有元素？据此，可以如何简化短路计算程序？
答：不需要，可以利用高斯消元法消去部分中间节点和联络节点，化简成等效网络再进行计算。

Js(2,1)=subs(J(2,1),{x1 x2},{Xn(1,1) Xn(2,1)});Js(2,2)=subs(J(2,2),{x1 x2},{Xn(1,1) Xn(2,1)});
Ja=-1*Js;
Xns=double(Ja\fn)
Xn=Xns+Xn
Xnsa(1)=abs(Xns(1))
Xnsa(2)=abs(Xns(2))
Xnsm=max(Xnsa)
t=t+1;
if Xnsm<a
b=0;
else b=1;
end
disp(Xn);
end
disp('迭代次数=')
disp(t);
disp('最终结果为')
disp(Xn);
2.上机验证课本例19.1和例19.2的计算结果。