-
大小: 42KB文件類(lèi)型: .zip金幣: 2下載: 0 次發(fā)布日期: 2023-07-20
- 語(yǔ)言: Matlab
- 標(biāo)簽: ieee118??節(jié)點(diǎn)數(shù)據(jù)??MATLAB??
資源簡(jiǎn)介
包含ieee118節(jié)點(diǎn)以及線(xiàn)路詳細(xì)數(shù)據(jù) 以及 MATLAB潮流計(jì)算程序

代碼片段和文件信息
%%??????運(yùn)用牛頓-拉夫遜法解IEEE118?節(jié)點(diǎn)潮流
%%%%%%???首先把原始文件分割成兩個(gè)獨(dú)立的純數(shù)據(jù)文件分別保存為Voltage.txt和Resis.txt
%%%%%%???計(jì)算n節(jié)點(diǎn)時(shí),當(dāng)Voltage.txt和Resis.txt準(zhǔn)備完畢后,只需將14改為n,28改為2*n,20改為線(xiàn)路條數(shù)即可
%%??讀取Voltage中的數(shù)據(jù)
clear?
load?Voltage.txt
Bus=Voltage(:1)?;?????????????????%設(shè)定節(jié)點(diǎn)數(shù)
Voltbase=Voltage(:2)?;????????????%設(shè)定電壓標(biāo)準(zhǔn)值
VType=Voltage(:5);????????????????%設(shè)定節(jié)點(diǎn)類(lèi)型
V=Voltage(:6)*100;????????????????????%節(jié)點(diǎn)電壓標(biāo)幺值
Angle=Voltage(:7);????????????????%電壓角
Pload=Voltage(:8);????????????????%負(fù)載有功
Qload=Voltage(:9);????????????????%負(fù)載無(wú)功
Pgen=Voltage(:10);????????????????%發(fā)電機(jī)發(fā)出有功
Qgen=Voltage(:11);????????????????%發(fā)電機(jī)發(fā)出無(wú)功
???????????????????????????????????%電壓設(shè)定點(diǎn)
Vset=Voltage(:13);????????????????%節(jié)點(diǎn)所接并聯(lián)電容器的電納
Qgmax=Voltage(:14);???????????????%Q的最大值
Qgmin=Voltage(:15);???????????????%Q的最小值
Qsh=Voltage(:17)*100;?????????????????%電容
baseMVA?=?100;?????????????????????%功率標(biāo)準(zhǔn)值
%%??讀取Resis中的數(shù)據(jù)
load?Resis.txt
II=Resis(:1);??????????????????
JJ=Resis(:2);??????????????????%兩端點(diǎn)
RType=Resis(:5);???????????????%線(xiàn)路類(lèi)型
R=Resis(:6);???????????????????%兩點(diǎn)間電阻
X=Resis(:7);???????????????????%兩點(diǎn)間電抗
Btwo=Resis(:8)/2;??????????????%線(xiàn)路兩端電納
K=Resis(:14);??????????????????%?變壓器變比
Kmax=Resis(:16);
Kmin=Resis(:17);?
Kset=Resis(:18);???????????????%變比要求
tic;
t=clock;
%%???開(kāi)始求導(dǎo)納?Y
y0=zeros(118)?;??????????????????%兩點(diǎn)間的電阻
yt=zeros(118)?;??????????????????%變壓器之間的導(dǎo)
Qsh=i*Qsh;??????????????????????%電容電抗標(biāo)幺值
K=1./K?;
for?x=1:179?
????switch?RType(x)?????????????%比較線(xiàn)路類(lèi)型??
????????case?2??????????????????%含變壓器的
????????????y0(II(x)JJ(x))=1./(K(x).*(R(x)+i*X(x)))?;???????????????????
????????????y0(JJ(x)II(x))=y0(II(x)JJ(x))????;???????????????
????????????yt(II(x)JJ(x))=(1-K(x))./((R(x)+i*X(x)).*K(x)^2);??%II端的電壓?
????????????yt(JJ(x)II(x))=(K(x)-1)./((R(x)+i*X(x)).*K(x));????%JJ端的電壓?
????????case?1?????????????????%不含變壓器的
????????????y0(II(x)JJ(x))=1/(R(x)+i*X(x));?????????????????
????????????y0(JJ(x)II(x))=1/(R(x)+i*X(x));??????????????????
????end
end
Y=zeros(118);??????????????????%導(dǎo)納矩陣
clear?x??
for?x=1:118
????Y(xx)=sum(y0(x:)+yt(x:))+i*Qsh(x)/baseMVA?;?????%求出對(duì)角導(dǎo)納(不含B/2)
end
Y=Y-y0;
clear?x
for?x=1:179
????Y(II(x)II(x))=Y(II(x)II(x))+i*Btwo(x)?;
????Y(JJ(x)JJ(x))=Y(JJ(x)JJ(x))+i*Btwo(x)?;??????????%加B
end
clear?x???????????????????????????%?Y導(dǎo)納矩陣結(jié)束
%%?????
U=zeros(1118);????????????????????????????%節(jié)點(diǎn)電壓
for?x=1:118
????U(x)=V(x)*exp(i*Angle(x)*pi/180)??;???%極坐標(biāo)電壓
end
clear?x
e=real(U);
f=imag(U);???????????%定義實(shí)部e跟虛部f?瓜分U
G=real(Y);
B=imag(Y);???????????%定義實(shí)部G跟虛部B?瓜分Y
D=ones(2361);???????%dQ?dP?dU^2?的矩陣
Ps=zeros(1118);?????
Qs=zeros(1118);?????????????????
N=0;
for?m=1:118
????Ps(m)=(Pgen(m)-Pload(m))/baseMVA;
????Qs(m)=(Qgen(m)-Qload(m))/baseMVA;
end
clear?m???????????????????????
while?max(abs(D))>0.00001????????????????????????%?開(kāi)始迭代過(guò)程求解
????for?m=1:118
????????switch?VType(m)?????????????????????%按照節(jié)點(diǎn)類(lèi)型分情況計(jì)算
????????????case?1
????????????????D(2*m-1)=Ps(m)-e(m)*sum(G(m:).*e-B(m:).*f)-f(m
?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件???????19411??2018-01-30?17:33??118節(jié)點(diǎn)電壓數(shù)據(jù).xlsx
?????文件???????27284??2018-01-21?10:20??118節(jié)點(diǎn)電阻數(shù)據(jù).xlsx
?????文件????????7271??2018-01-31?09:51??ieee118.m
評(píng)論
共有 條評(píng)論