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

昆明理工大学电力工程学院学生实验报告实验课程名称：《电力系统计算机辅助分析》

开课实验室：计算机中心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

b=0;

else b=1;

end

disp(Xn);

end

disp('迭代次数=')

disp(t);

disp('最终结果为')

disp(Xn);

2.上机验证课本例19.1和例19.2的计算结果