2013级电力系统分析课程设计(3)

2019-01-27 17:11

1.818182

0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.818182 -2.651515 0.000000

0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.818182 0.000000 -2.651515 Zif如下 0.516426 0.535359 0.535359 0.197708 0.106525 0.057251 0.060176 0.051292 0.035171 0.035171

各节点电压为： U 1 = 0.316899 U 2 = 0.331537 U 3 = 0.331537 U 4 = 0.070494 U 5 = 0.000000 U 6 = 0.462564 U 7 = 0.435105 U 8 = 0.518503 U 9 = 0.669831 U 10 = 0.669831

各支路电流I为

0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000

0.000000 0.000000 0.000000 0.000000 -0.487930 0.000000 0.000000 0.000000 0.000000 0.000000

0.000000 0.000000 0.000000 0.487930 0.000000 3.116039 2.808605 0.000000 0.000000 0.000000

0.000000 0.000000 0.000000 0.000000 -3.116039 0.000000 -2.258322 0.000000 0.000000 0.000000

0.000000 0.000000 0.000000 0.000000 -2.808605 2.258322 0.000000 0.550282 0.000000 0.000000

0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -0.550282 0.000000 0.000000 0.000000

0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000

实验的结果比较理想，与实际所算出的结果只有微小的误差，已经足以忽略不计了。

六、设计体会

通过对这次电力系统分析课程设计，我对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练。刚开始进行的时候，感觉到根本无从下手，难度有点偏大，但是经过老师和同学的指导，我从图书馆借来相关的参考书来进行学习。Matlab虽然之前做电机课设的时候接触过，但是那次我用的是仿真功能，而这次用的是编程功能，在经过看书和实践操作后，我从人机界面的设计，到手动输入，到Y阵的生成，在到保存文件为TXT格式，每一步都要新的学习，每一步都是新的尝试。在后期的调试，将各个部分组成一个系统的程序，兼顾各个部分的兼容性，可以说实在一次次的尝试与失败中走过来的。在最后终于成功的那一刻，感觉之前的付出都是值得的。能够完成这次课程设计离不开老师和同学的帮助，我也会在以后的学习生活中更加努力。

七、参考文献

【1】电力系统分析（上册）何仰赞温增银，华中科技大学出版社，2002 【2】MATLAB GUI设计学习手册罗华飞，北京航空航天大学, 2011

【3】MATLAB 7.6从入门到精通张琨毕靖丛滨，电子工业出版社, 2009 【4】精通MATLAB 葛哲学，电子工业出版社, 2008

八、附录（主程序及其注释）

Sb=100; Vb=115;

n=str2num(get(handles.jiedian,'String')); f=str2num(get(handles.duanludian,'String')); p=str2num(get(handles.fdj,'String')); q2=str2num(get(handles.byq2,'String')); s=str2num(get(handles.xt,'String')); q3=str2num(get(handles.byq3,'String')); x=str2num(get(handles.xl,'String'));

%手动输入

F=zeros(10,10);

F=str2num(get(handles.jdz,'String')); G=zeros(10,10);

G=str2num(get(handles.fdjz,'String')); H=zeros(10,10);

H=str2num(get(handles.fdjs,'String')); % I=zeros(10,10);

% I=str2num(get(handles.fdjv,'String')); J=zeros(10,10);

J=str2num(get(handles.byq2z,'String')); K=zeros(10,10);

K=str2num(get(handles.byqz3,'String')); L=ones(10,10);

L=str2num(get(handles.xtz,'String')); global imagz; global z; global gz; global gs; global byqy; global byqy3; global sz;

imagz=zeros(n,n); z=zeros(n,n); gz=zeros(p,2); gs=zeros(p,1); byqy=zeros(q2,4); byqy3=zeros(q3,10); sz=zeros(s,2);

%标幺化线路阻抗 for nn=1:x i=F(nn,1); j=F(nn,2);

imagz(i,j)=F(nn,3)*F(nn,4)*Sb/Vb^2; imagz(j,i)=imagz(i,j); end

z=sqrt(-1)*imagz;

%标幺化发电机参数 for pp=1:p

gz(pp,1)=G(pp,1);

gz(pp,2)=G(pp,2)/H(pp,1)*Sb; end

%标幺化变压器参数 for qq=1:q2

byqy(qq,1)=J(qq,1); byqy(qq,2)=J(qq,2);

byqy(qq,4)=J(qq,4)/100*Sb/J(qq,3); end

%标幺化三绕组变压器参数 for qq3=1:q3

byqy3(qq3,1)=K(qq3,1); byqy3(qq3,2)=K(qq3,2); byqy3(qq3,3)=K(qq3,3); byqy3(qq3,4)=K(qq3,4);

byqy3(qq3,8)=1/2*(K(qq3,8)+K(qq3,9)-K(qq3,10))/100*Sb/K(qq3,5); byqy3(qq3,9)=1/2*(K(qq3,8)+K(qq3,10)-K(qq3,9))/100*Sb/K(qq3,6); byqy3(qq3,10)=1/2*(K(qq3,10)+K(qq3,9)-K(qq3,8))/100*Sb/K(qq3,7); end

%标幺化系统阻抗 for ss=1:s

sz(ss,1)=L(ss,1); sz(ss,2)=Sb/L(ss,2); end

Y=zeros(n,n); %节点导纳矩阵 V=zeros(n,1); %节点电压矩阵 I=zeros(n,n); %支路电流矩阵

for i=1:n %以下生成节点导纳矩阵 for j=1:n

if z(i,j)~=0

Y(i,i)=Y(i,i)+1/z(i,j); end end end

%考虑发电机参数影响 for pp=1:p

i=gz(pp,1);

Y(i,i)=Y(i,i)+1/sqrt(-1)/gz(pp,2); end

%考虑系统参数影响 for ss=1:s

i=sz(ss,1);

Y(i,i)=Y(i,i)+1/sqrt(-1)/sz(ss,2); end

for i=1:n %计算非对角线元素 for j=i+1:n

if z(i,j)==0 Y(i,j)=0;

else Y(i,j)=-1/z(i,j); end

Y(j,i)=Y(i,j); end end

%考虑双绕组变压器参数影响 for qq=1:q2 i=byqy(qq,1); j=byqy(qq,2);

Y(i,i)=Y(i,i)+1/sqrt(-1)/byqy(qq,4); Y(j,j)=Y(j,j)+1/sqrt(-1)/byqy(qq,4);

共4页:

2013级电力系统分析课程设计(3).doc 将本文的Word文档下载到电脑下载失败或者文档不完整，请联系客服人员解决！

下载这篇word文档