攀枝花学院本科毕业设计（论文） 结论
5）对用注入功率模拟线路开断的方法的思考
本文中用注入功率模拟线路开断，相对于变导纳矩阵的方法固然是中进步，但关于用注入功率模拟线路开断的方法是否能够在理论上做出真实模拟，提出疑问 ：一条线路开断将引起全网功率的分布的改变，而用注入功率模拟线路的开断只是对被模拟开断线路两侧节点进行了节点注入功率的修改，理论上并未对线路开断进行真实的模拟。对于该问题，在后续的研究中，初步拟定具体测试步骤如下：首先，在线路参数中删除一条线路的参数，即该线路不存在，即断开了，此时运行潮流程序对各个节点的电压幅值相角求解；然后，用修改注入功率模拟线路开断，两种情况下比对即可得出结论。
6）程序的适用性
关于程序的适用性上，此程序只限于5节点数据，是由于其初始化导纳矩阵维数的定义上的限制，如果对程序中的变量维数进行优化，维数用一具体变量替换，则该程序的适用性将不受节点个数的限制。此外由于时间和基础原因，程序运行结果不是友好界面，只是输出最终的运行数据和判断结果，即应用中至少应该是图形界面的。这是可以提高的。最后一点思考到的不足是，没有预留实时读取数据的接口函数。基础数据为手动输入数据，在真实开断模拟中，应实时跟踪采集电网数据，周期性的进行静态安全分，实时显示可能故障的支路，已供调度人员参考。
7）静态安全分析方法的优点
在本文的静态安全分析方法中，潮流迭代算法只用在第一次求系统正常运行状态参数上。随后的静态安全分程序则应用公式求出电压和相角的修正量，直接获得开断后节点的电压和相角。如果针对本文测试系统，在每一个断线预想故障状态下，都进行一次潮流计算的话，则会增加34次的迭代运算。因此，本文所研究的静态安全分析方法少了34次的迭代运算，这是在程序算法上的提高，运行效率提升很明显，也充分体现了本文所研究的静态安全分析方法的优点。
27
攀枝花学院本科毕业设计（论文） 参考文献
参 考 文 献
［1］KUNDUR. P． Power System Stability and Control [M]. 北京：中国电力出版社，2002. ［2］吴际舜．电力系统静态安全分析［M］．上海交通大学出版社，1985. ［3］王锡凡, 方万良, 杜正春．现代电力系统分析［M］．科学出版社，2003.
［4］吴天鹏，谢小介，彭斌. MATLAB电力系统设计与分析［M］. 北京:国防工业出版社,2004. ［5］王正林，刘明. 精通MATLAB7［M］.电子工业出版社，2006. ［6］何仰赞．电力系统分析［M］．武汉:武汉华中科技大学出版社，1985. ［7］于永源．电力系统分析［M］．北京：中国电力出版社，2007.
［8］陈珩．电力系统稳态分析(第三版)［M］. 北京：中国电力出版社，2007. ［9］李光琦. 电力系统暂态分析［M］．北京:中国电力出版社，1995.
［10］李文沅. 电力系统风险评估模型-方法和应用［M］．北京:科学出版社，2006. ［11］李文沅. 电力系统经济安全运行模型与方法［M］．重庆:重庆大学出版社，1989. ［12］周荣光．电力系统故障分析［M］．北京:北京清华大学出版社，1988. ［13］中华人民共和国国家经济贸易委员会. 电力系统安全稳定导则DL755-2001. ［14］许高达．电力系统之计算机应用［M］．新兴图书公司，1979. ［15］朱仁峰. 精通MATLAB7［M］.清华大学出版社，2006.
［16］Math Works. 产品和服务［EB/OL］. http://www.mathworks.cn/products/matlab/，2014-5-13. ［17］刘浏. 百度百科［EB/OL］．http:// baike.http://www.wodefanwen.com//view/10598.htm?fr=aladdin,2014-5-13.
28
攀枝花学院本科毕业设计（论文） 附录A:基础数据生成子程序
附录A:基础数据生成子程序
程序_基础数据生成子程序
function[p,q,voltage,ang,g,b0,b1,b2,y0,node_f,node_e]=basic_data();
%首节点 末节点 支路阻抗 对地导纳的1/2 变比
Z=[
1 2 0.04+j*0.25 j*0.25 1 3 0.1+j*0.35 0 2 3 0.08+j*0.3 j*0.25 4 2 j*0.015 0 5 3 j*0.03 0 ];
admittance=zeros(5); %初始化导纳矩阵 for m=1:5
p=Z(m,1); %p=1,1,2,4,5 q=Z(m,2); %q=2,3,3,2,3
admittance(p,q)=admittance(p,q)-1./(Z(m,3)*Z(m,5));
admittance(q,p)=admittance(p,q); %互导纳相等 %计算自导纳
admittance(q,q)=admittance(q,q)+1./(Z(m,3)*Z(m,5)^2)+Z(m,4); %是否添加对地导纳参数%计算高压侧自导纳 admittance(p,p)=admittance(p,p)+1./Z(m,3)+Z(m,4); %是否添加对地导纳参数%计算低压侧自导纳 end
conductance=real(admittance); %电导英文:conductance susceptance=imag(admittance); %电纳英文:susceptance b1=susceptance(1:4,1:4); %4*4
b1=inv(b1);%P-Q分解法中对应p中的电纳矩阵//仅用作求逆 b2=susceptance(1:3,1:3);
b2=inv(b2);%P-Q分解法中对应q的电纳矩阵 b0=susceptance;
29
1 1.05 1.05 1
1
攀枝花学院本科毕业设计（论文） 附录A:基础数据生成子程序
g=conductance;
var=[%发电机功率 负荷功率 电压初值 节点类型
0 1.6-j*0.8 1 2 0 2-j 1 2 0 3.7-j*1.3 1 2 5 0 1.05 3 0 0 1.05 1 ];
%/1平衡节点，2P-Q节点，3pv节点 voltage=zeros(5,1); %创建电压列向量 v=zeros(4,1);
ang =zeros(5,1); %创建相角列向量 voltage=var(:,3); %对电压列向量初始化
for m=1:5
s(m)=var(m,1)-var(m,2); end s=s'; p=real(s); q=imag(s);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算支路潮流所需的参数
y0 =Z(:,4); %取对地导纳矩阵 y0 =imag(y0); %对地导纳的共轭矩阵 node_f=zeros(5,1); node_e=zeros(5,1); node_f=Z(:,1); node_e=Z(:,2);
[voltage,ang]=P-Q(voltage,ang,p,q,g,b0,b1,b2); Return
30
攀枝花学院本科毕业设计（论文） 附录B:潮流计算子程序
附录B：潮流计算子程序
%,输入voltage,ang,p,q,g,b0,b1,b2,输出voltage,ang function [voltage,ang]=P-Q(voltage,ang,p,q,g,b0,b1,b2) %volatge,ang,为迭代计算的初值； iteration=1;
tol=1e-6; %%收敛容许误差 while(iteration<20) v=voltage(1:4,1); p_input=zeros(4,1);
for m=1:4 for n=1:5
p_input(m)=voltage(n)*(g(m,n)*cos(ang(m)-ang(n))+b0(m,n)*sin(ang(m)-ang(n)))+p_input(m);
end
p_input(m)=p_input(m)*voltage(m); end
p_correction=zeros(4,1); p_correction=p(1:4,1)-p_input;
ang_correction=zeros(4,1);
ang_correction=b1*p_correction./v./v;
ang(1:4,1)=ang(1:4,1)-ang_correction;
q_input=zeros(3,1); for m=1:3 for n=1:5
q_input(m,1)=voltage(n,1)*(g(m,n)*sin(ang(m)-ang(n))-b0(m,n)*cos(ang(m)-ang(n)))
31