電力系統(tǒng)穩(wěn)態(tài)分析基于matlab的潮流計算快速分解法-在線瀏覽
2025-03-02 14:14本頁面
【正文】 ，有 （16）對于一般的電網(wǎng)，可能有較復(fù)雜的結(jié)構(gòu)。下面將證明，對于樹形電網(wǎng)或所有支路的比值都相同的環(huán)形網(wǎng)絡(luò)，與相等。這說明對樹形電網(wǎng)，就是XB型快速分解法中的陣。但和相比，在的零元素處，相應(yīng)的元素近似等于零；在的非零元素處，相應(yīng)的元素近似和的非零元素相等。圖1 四節(jié)點(diǎn)電力系統(tǒng)以圖1所示的四節(jié)點(diǎn)系統(tǒng)為例，圖中給出了支路阻抗。以上推導(dǎo)過程中，只在對有功功率偏差量和無功功率偏差量的計算處做了線性化近似，既沒有rx的假設(shè)，也沒有引用PQ解耦的假設(shè)。以100MVA、220KV為基準(zhǔn)的節(jié)點(diǎn)初始標(biāo)幺值和網(wǎng)絡(luò)參數(shù)見表2。圖2 電力系統(tǒng)模型 節(jié)點(diǎn)12345PQ表1 節(jié)點(diǎn)負(fù)荷值支路rxb變比12130．000232435表2 線路參數(shù)4．程序流程圖程序開始輸入節(jié)點(diǎn)和支路數(shù)據(jù)數(shù)據(jù)建立節(jié)點(diǎn)導(dǎo)納矩陣形成矩陣和置迭代次數(shù)有功與無功的迭代偏差均大于計算無功不平衡量計算相角修正量，求得計算有功不平衡量計算電壓修正量，求得計算平衡節(jié)點(diǎn)功率有功迭代偏差小于無功迭代偏差小于無功迭代偏差是否小于？有功迭代偏差是否小于？？輸出數(shù)據(jù)是否是是否否5．MatLab程序清單本文采用兩種程序來仿真電力系統(tǒng)的潮流，兩個程序得到的仿真結(jié)果是相同的。%節(jié)點(diǎn)數(shù)x[y]=5。%誤差要求point=[1 1 0 。 1 1 0 。 2 0 0 ]。 1 3 0 0。 2 4 0 。]。%TYPE為節(jié)點(diǎn)類型矩陣U=zeros(x,1)。%a為節(jié)點(diǎn)電壓相角矩陣P=zeros(x,1)。%Q為節(jié)點(diǎn)無功功率I=zeros(y,1)。%J為終止節(jié)點(diǎn)編號矩陣Rij=zeros(y,1)。%X為線路電抗Zij=Rij+j*Xij。%Y為n階節(jié)點(diǎn)導(dǎo)納方陣G=zeros(x)。%B為n階節(jié)點(diǎn)電納方陣B0=zeros(y,1)。 %K為ij支路y*1階變壓器變比,若K=0表示無變壓器,K=1則為標(biāo)準(zhǔn)變比,K不等于1為非標(biāo)準(zhǔn)變比%矩陣賦初值：TYPE=point(:,1)。O=point(:,3)。Q=point(:,5)。J=branch(:,2)。Xij=branch(:,4)。B0=branch(:,5)。%求節(jié)點(diǎn)導(dǎo)納矩陣Yfor m=1:y %求Y中非對角元的元素Yij Y(I(m),J(m))=1/Zij(m)。end。 % Yii為Yij加上線路對地電導(dǎo)的一半乘j end。 end。%求B39。%求Y的虛部,節(jié)點(diǎn)電納矩陣 B11=zeros(x,x)。 B11(J(m),I(m))=1/(Xij(m))。 for m=1:x for n=1:y if(I(n)==m|J(n)==m) B11(m,m)=B11(m,m)1/(Xij(n))。 end。 p=find(TYPE(:,1)==3)。%平衡節(jié)點(diǎn)編號對應(yīng)行置空 B11(p,:)=[]。 B1=inv(B1)。39。%找出非PQ節(jié)點(diǎn)的編號 B22=B。%非PQ節(jié)點(diǎn)編號對應(yīng)行置空 B22(npq,:)=[]。 B2=inv(B2)。 %k為迭代次數(shù) ep=0。 %計算U不平衡量dQi的收斂標(biāo)志(0表示不收斂,1表示收斂) np=find(TYPE(:,1)3)。%dPi為x*1階矩陣,x即為節(jié)點(diǎn)數(shù) pq=find(TYPE(:,1)==1)。 pqn=pqn(1)。%dQi為pqn*1階矩陣while((~eq|~ep)amp。 for m=1:(x1)%求dPi sum1=0。 end dPi(m)=P(np(m))sum1。 %Up為中間變量 Up(p)=[]。%求得除平衡節(jié)點(diǎn)外的電壓列向量 dO=(B1*(dPi./Unp))./Unp。 end max1=abs(dPi(1)/U(np(1)))。 end end if max1=e %如果最大值滿足要求,則ep置為1,表示收斂 ep=1。 for n=1:x sum2=sum2+U(pq(m))*U(n)*(G(pq(m),n)*sin(O(pq(m))
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1