該樓層疑似違規(guī)已被系統(tǒng)折疊?隱藏此樓查看此樓

n=input('Please input n\n'); %n表示系統(tǒng)的節(jié)點(diǎn)數(shù)

d0=input('Please input d0\n'); %d0表示系統(tǒng)的誤差

m=input('Please inpit m\n'); %m表示系統(tǒng)的pq節(jié)點(diǎn)的個(gè)數(shù)

p0=ones(n-1,1);

q0=ones(m,1);

d=1;

dpv=ones(n-1,1);

dqv=ones(m,1);

dp=zeros(n-1,1);

dq=zeros(m,1);

for a=1:n

for b=1:n

G(a,b)=real(Y(a,b));

B(a,b)=imag(Y(a,b));

end

end %把Y矩陣的實(shí)部和虛部分別導(dǎo)出

for a=1:n-1

for b=1:n-1

B1(a,b)=B(a,b);

end

end %表示出B1

for a=1:m

for b=1:m

B2(a,b)=B(a,b);

end

end %表示出B2

while(d>d0)

for a=1:n-1

p0(a)=0;

for b=1:n

p0(a)=p0(a)+v(a)*v(b)*(G(a,b)*cos(w(a)-w(b))+B(a,b)*sin(w(a)-w(b)));

end

end %利用初值求得相應(yīng)的功率

for a=1:m

q0(a)=0;

for b=1:n

q0(a)=q0(a)+v(a)*v(b)*(G(a,b)*sin(w(a)-w(b))-B(a,b)*cos(w(a)-w(b)));

end

end %利用初值求得相應(yīng)的無功

for a=1:n-1

dp(a)=p(a)-p0(a);

end %求dp

for a=1:m

dq(a)=q(a)-q0(a);

end %求dq

for a=1:n-1

dpv(a)=dp(a)/v(a);

end %利用初值求出dp/v

for a=1:m

dqv(a)=dq(a)/v(a);

end %利用初值求出dq/v

dw=-1*B1\dpv;

for a=1:n-1

dw(a)=dw(a)/v(a);

end %求出dw的修正值

dv=-1*B2\dqv; %求出v的修正值

d=max(max(abs(dw)),max(abs(dv)));%誤差衡量

for a=1:n-1

w(a)=w(a)+dw(a);

end %求修正后的ｗ

for a=1:m

v(a)=v(a)+dv(a);

end %求修正后的v

end

這個(gè)是基于PQ分解法解潮流方程，我是按照自己的習(xí)慣編寫的==。在運(yùn)行本程序前，需要自行建立p(有功)、q(無功)、Y(節(jié)點(diǎn)導(dǎo)納矩陣)、w(相角)、v(電壓幅值)矩陣。

其中，節(jié)點(diǎn)編號(hào)請(qǐng)按照PQ節(jié)點(diǎn)、PV節(jié)點(diǎn)、平衡節(jié)點(diǎn)的次序編號(hào)。為了獲得更好的收斂，可以將w統(tǒng)一設(shè)置為0，v統(tǒng)一設(shè)置為1。PQ解耦法適用于大電力系統(tǒng)，對(duì)于小型配電網(wǎng)等電力系統(tǒng)可能會(huì)存在較大的計(jì)算誤差。

總結(jié)

以上是生活随笔為你收集整理的matlab电力系统潮流计算,大神们，求个电力系统潮流计算的matlab程序。的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。