当前位置:首页 >> 电力/水利 >>

天津大学-电力系统快速解耦法潮流分析及短路计算--matlab程序编程


电力系统快速解耦法潮流分析及短路计算 一.程序设计的基本思想: 由于电力系统潮流分析中要利用到矩阵运算,复数运算,故采用 matlab 编程。 二、潮流计算框图 开始 输入系统的支路参数和 PQ 参数 计算节点导纳矩阵 Yn 设置初值 U、θ 、P、Q 设置收敛条件及最大迭 代次数 Kmax=10 计算雅克比矩阵-B’与-B’’ Y 判断是否符合精度 N 求解修正方程得到Δ U 和Δ θ 输出潮流不收敛信息 对 U、θ 进行修正, U=U+Δ U,θ =θ +Δ θ 结束 K=K+1 输出结果 K>Kmax? Y 计算潮流分布 三、matlab 程序 N clear all %%计算节点导纳矩阵 network=[1 4 0 0.0576 0; 2 7 0 0.0625 0; 3 9 0 0.0586 0; 4 5 0.01 0.085 0.088; 4 6 0.017 0.092 0.079; 5 7 0.032 0.161 0.153; 6 9 0.039 0.17 0.179; 7 8 0.0085 0.072 0.0745; 8 9 0.0119 0.1008 0.1045]; Yn=zeros(9); [Height Width]=size(network); for i=1:Height m=network(i,1); n=network(i,2); Yn(m,m)=Yn(m,m)+1/(network(i,3)+1j*network(i,4))+1j*network(i,5); Yn(n,n)=Yn(n,n)+1/(network(i,3)+1j*network(i,4))+1j*network(i,5); Yn(m,n)=Yn(m,n)-1/(network(i,3)+1j*network(i,4)); Yn(n,m)=Yn(n,m)-1/(network(i,3)+1j*network(i,4)); end disp('节点导纳矩阵') Yn %%快速解耦法潮流计算 %B',B''矩阵计算 B=imag(Yn); B1=zeros(n-1); for i=1:n-1 for j=1:n-1 B1(i,j)=B(i+1,j+1); end end m=6; B2=zeros(m); for i=1:m for j=1:m B2(i,j)=B(n-m+i,n-m+j); end end B1 B2 %B矩阵n阶 %B'矩阵n-1阶 %B’’矩阵m阶 %%设置初值 U=[1.04,1.025,1.025,1,1,1,1,1,1]'; theta=[0,0,0,0,0,0,0,0,0]'; P=[1.63,0.85,0,-1.25,-0.9,0,-1,0]'; Q=[0,-0.5,-0.3,0,-0.35,0]'; delta=ones(14,1); delta2=ones(14,1); G=real(Yn); B=imag(Yn); k=0; Kmax=10; epsilon=1e-5; thetaij=zeros(8,8); while k<=Kmax&&max(abs(delta))>epsilon Utemp_1=zeros(9,1); Utemp_2=zeros(9,1); for i=1:9 for j=1:9 thetaij(i,j)=theta(i)-theta(j); Utemp_1(i)=Utemp_1(i)+U(j)*(G(i,j)*sin(thetaij(i,j))-B(i,j)*cos

相关文章:
更多相关标签: