資源簡介
可以采用蒙特卡絡法對實際配電網例子進行實驗，驗證了方法的正確性，10負荷點配電網蒙特卡洛法可靠性計算程序
代碼片段和文件信息
clc
disp（‘本配電網絡中共有線路（條）：‘）
Nline=10
disp（‘共有變壓器（臺）：‘）
Ntrans=10
disp（‘共有斷路器（臺）：‘）
Nbr=5
disp（‘配電網絡中共有元件數：‘）
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（‘假設配電線路、變壓器、斷路器的故障率分別為：‘）
Kline=0.065
Ktrans=0.015
Kbr=0.006
%修復時間
disp（‘假設配電線路、變壓器、斷路器的修復時間分別為：‘）
Tline=5
Ttrans=200
Tbr=4
%初始化負荷點參數
lpgzcs=zeros（1Ntrans）;
lpgzsj=zeros（1Ntrans）;
lpgzl=zeros（1Ntrans）;
lpgzcxsj=zeros（1Ntrans）;
lpngzsj=zeros（1Ntrans）;
%初始化10個線路元件的故障率和修復時間
for?i=1:Nline
????yjgzl（i）=Kline*lth（i）;
????yjxfsj（i）=Tline;
end
%初始化10個變壓器元件的故障率和修復時間
for?i=（Nline+1）:（Nline+Ntrans）
????yjgzl（i）=Ktrans;
????yjxfsj（i）=Ttrans;
end
%初始化5個開關元件的故障率和修復時間
for?i=（Nline+Ntrans+1）:N
????yjgzl（i）=Kbr;
????yjxfsj（i）=Tbr;
end
%各個負荷點的供電路徑
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
????
????%獲得一對對偶數
????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;
????
????%計算各個負荷點的故障次數與故障時間
????%計算負荷點1的故障次數與故障時間
????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
??????
????%計算負荷點2的故障次數與故障時間
????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
????
????%計算負荷點3的故障次數與故障時間
????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
????
????%計算負荷點4的故障次數與故障時間
????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
???????????屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????6087??2017-11-27?20:45??配電系統可靠性計算\MATLAB_reliability.m
?????文件??????37888??2017-11-27?21:04??配電系統可靠性計算\reliability\Debug\reliability.exe
?????文件?????345708??2017-11-27?21:04??配電系統可靠性計算\reliability\Debug\reliability.ilk
?????文件?????461824??2017-11-27?21:04??配電系統可靠性計算\reliability\Debug\reliability.pdb
?????文件????2359296??2017-11-27?21:04??配電系統可靠性計算\reliability\ipch\reliability-853d85eb\reliability-7f0cc46.ipch
?????文件???????2390??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\cl.command.1.tlog
?????文件???????7166??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\CL.read.1.tlog
?????文件???????1628??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\CL.write.1.tlog
?????文件??????????2??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\li
?????文件??????????2??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\li
?????文件???????5474??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\li
?????文件??????10832??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\li
?????文件???????3416??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\li
?????文件???????1400??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\mt.command.1.tlog
?????文件???????1358??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\mt.read.1.tlog
?????文件????????426??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\mt.write.1.tlog
?????文件???????2060??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\rc.command.1.tlog
?????文件???????1274??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\rc.read.1.tlog
?????文件???????1298??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\rc.write.1.tlog
?????文件????????406??2017-11-27?20:32??配電系統可靠性計算\reliability\reliability\Debug\reliability.exe.em
?????文件????????472??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\reliability.exe.em
?????文件????????381??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\reliability.exe.intermediate.manifest
?????文件????????114??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\reliability.lastbuildstate
?????文件???????2834??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\reliability.log
?????文件??????26458??2017-11-27?21:04??配電系統可靠性計算\reliability\reliability\Debug\reliability.obj
?????文件??????????0??2017-11-27?20:32??配電系統可靠性計算\reliability\reliability\Debug\reliability.write.1.tlog
............此處省略24個文件信息
- 上一篇：abaqus軟件對隧道開挖過程的模擬
- 下一篇：長沙市POI數據.zip
評論
共有 條評論
川公網安備 51152502000135號