資源簡(jiǎn)介
可以采用蒙特卡絡(luò)法對(duì)實(shí)際配電網(wǎng)例子進(jìn)行實(shí)驗(yàn)，驗(yàn)證了方法的正確性，10負(fù)荷點(diǎn)配電網(wǎng)蒙特卡洛法可靠性計(jì)算程序
代碼片段和文件信息
clc
disp（‘本配電網(wǎng)絡(luò)中共有線路（條）：‘）
Nline=10
disp（‘共有變壓器（臺(tái)）：‘）
Ntrans=10
disp（‘共有斷路器（臺(tái)）：‘）
Nbr=5
disp（‘配電網(wǎng)絡(luò)中共有元件數(shù)：‘）
N=Nline+Ntrans+Nbr
lth=[0.690.290.421.550.982.100.451.251.131.90];???%length?of?lines
load=[0.670.870.010.221.000.570.820.890.870.68];???%load?of?each?load?point
lpuser=[3820482022867437372];??%number?of?the?users
%故障率
disp（‘假設(shè)配電線路、變壓器、斷路器的故障率分別為：‘）
Kline=0.065
Ktrans=0.015
Kbr=0.006
%修復(fù)時(shí)間
disp（‘假設(shè)配電線路、變壓器、斷路器的修復(fù)時(shí)間分別為：‘）
Tline=5
Ttrans=200
Tbr=4
%初始化負(fù)荷點(diǎn)參數(shù)
lpgzcs=zeros（1Ntrans）;
lpgzsj=zeros（1Ntrans）;
lpgzl=zeros（1Ntrans）;
lpgzcxsj=zeros（1Ntrans）;
lpngzsj=zeros（1Ntrans）;
%初始化10個(gè)線路元件的故障率和修復(fù)時(shí)間
for?i=1:Nline
????yjgzl（i）=Kline*lth（i）;
????yjxfsj（i）=Tline;
end
%初始化10個(gè)變壓器元件的故障率和修復(fù)時(shí)間
for?i=（Nline+1）:（Nline+Ntrans）
????yjgzl（i）=Ktrans;
????yjxfsj（i）=Ttrans;
end
%初始化5個(gè)開關(guān)元件的故障率和修復(fù)時(shí)間
for?i=（Nline+Ntrans+1）:N
????yjgzl（i）=Kbr;
????yjxfsj（i）=Tbr;
end
%各個(gè)負(fù)荷點(diǎn)的供電路徑
lp1=[11121];
lp2=[121221];
lp3=[1231321];
lp4=[1234142122];
lp5=[12345152122];
lp6=[61623];
lp7=[671723];
lp8=[678182324];
lp9=[6789192324];
lp10=[610202325];
SUMMTTF=0;
while?SUMMTTF<=20000
????
????%獲得一對(duì)對(duì)偶數(shù)
????R1=rand（1N）;
????for?i=1:N
????????R2（i）=1-R1（i）;
????end
????
????%求各元件MTTF
????for?i=1:N
????????MTTF1（i）=-1/yjgzl（i）*log（R1（i））;
????end
????for?i=1:N
????????MTTF2（i）=-1/yjgzl（i）*log（R2（i））;
????end
????
????%確定故障元件，求取MTTR
????a1=MTTF1（1）;
????b1=1;
????for?i=1:N
????????if?a1>MTTF1（i）
????????????a1=MTTF1（i）;
????????????b1=i;
????????end
????end
????MTTR1=-yjxfsj（b1）*log（rand（1））;
????
????a2=MTTF2（1）;
????b2=1;
????for?i=1:N
????????if??a2>MTTF2（i）;
????????????a2=MTTF2（i）;
????????????b2=i;
????????end
????end
????MTTR2=-yjxfsj（b2）*log（rand（1））;
???
????SUMMTTF=SUMMTTF+a1+a2;
????
????%計(jì)算各個(gè)負(fù)荷點(diǎn)的故障次數(shù)與故障時(shí)間
????%計(jì)算負(fù)荷點(diǎn)1的故障次數(shù)與故障時(shí)間
????for?i=1:3
????????if?lp1（i）==b1
????????????lpgzcs（1）=lpgzcs（1）+1;
????????????lpgzsj（1）=lpgzsj（1）+MTTR1;
????????end
????end
????for?i=1:3
????????if?lp1（i）==b2
????????????lpgzcs（1）=lpgzcs（1）+1;
????????????lpgzsj（1）=lpgzsj（1）+MTTR2;
????????end
????end
??????
????%計(jì)算負(fù)荷點(diǎn)2的故障次數(shù)與故障時(shí)間
????for?i=1:4
????????if?lp2（i）==b1
????????????lpgzcs（2）=lpgzcs（2）+1;
????????????lpgzsj（2）=lpgzsj（2）+MTTR1;
????????end
????end
????for?i=1:4
????????if?lp2（i）==b2
????????????lpgzcs（2）=lpgzcs（2）+1;
????????????lpgzsj（2）=lpgzsj（2）+MTTR2;
????????end
????end
????
????%計(jì)算負(fù)荷點(diǎn)3的故障次數(shù)與故障時(shí)間
????for?i=1:5
????????if?lp3（i）==b1
????????????lpgzcs（3）=lpgzcs（3）+1;
????????????lpgzsj（3）=lpgzsj（3）+MTTR1;
????????end
????end
????for?i=1:5
????????if?lp3（i）==b2
????????????lpgzcs（3）=lpgzcs（3）+1;
????????????lpgzsj（3）=lpgzsj（3）+MTTR2;
????????end
????end
????
????%計(jì)算負(fù)荷點(diǎn)4的故障次數(shù)與故障時(shí)間
????for?i=1:7
????????if?lp4（i）==b1
????????????lpgzcs（4）=lpgzcs（4）+1;
????????????lpgzsj（4）=lpgzsj（4）+MTTR1;
????????end
????end
????for?i=1:7
????????if?lp4（i）==b2
???????????屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件???????6087??2017-11-27?20:45??配電系統(tǒng)可靠性計(jì)算\MATLAB_reliability.m
?????文件??????37888??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\Debug\reliability.exe
?????文件?????345708??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\Debug\reliability.ilk
?????文件?????461824??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\Debug\reliability.pdb
?????文件????2359296??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\ipch\reliability-853d85eb\reliability-7f0cc46.ipch
?????文件???????2390??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\cl.command.1.tlog
?????文件???????7166??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\CL.read.1.tlog
?????文件???????1628??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\CL.write.1.tlog
?????文件??????????2??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\li
?????文件???????5474??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\li
?????文件??????10832??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\li
?????文件???????3416??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\li
?????文件???????1400??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\mt.command.1.tlog
?????文件???????1358??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\mt.read.1.tlog
?????文件????????426??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\mt.write.1.tlog
?????文件???????2060??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\rc.command.1.tlog
?????文件???????1274??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\rc.read.1.tlog
?????文件???????1298??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\rc.write.1.tlog
?????文件????????406??2017-11-27?20:32??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\reliability.exe.em
?????文件????????472??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\reliability.exe.em
?????文件????????381??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\reliability.exe.intermediate.manifest
?????文件????????114??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\reliability.lastbuildstate
?????文件???????2834??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\reliability.log
?????文件??????26458??2017-11-27?21:04??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\reliability.obj
?????文件??????????0??2017-11-27?20:32??配電系統(tǒng)可靠性計(jì)算\reliability\reliability\Debug\reliability.write.1.tlog
............此處省略24個(gè)文件信息
評(píng)論
共有 條評(píng)論
川公網(wǎng)安備 51152502000135號(hào)