【正文】 謝開貴，中壓配電網(wǎng)可靠性的模糊評估，重慶大學(xué)學(xué)報 (自然科學(xué)版 )2021， 29(2)，4549. 【 10】 Billin ， Allan R N， Reliability Evaluation of Power System. Boston， Pistman， 1984. 【 11】別朝紅，王錫凡，配電系統(tǒng)的可靠性分析，中國電力， 1997， 30(5)， 1013. 【 12】張焰，電網(wǎng)規(guī)劃中的模糊可靠性評估方法，中國電機(jī)工程學(xué)報， 2021， 20(11)， 7780. 【 13】劉 莉，陳學(xué)允，基于模糊遺傳算法的配電網(wǎng)絡(luò)重構(gòu) [J]，中國電機(jī)工程學(xué)報， 2021， 20(2)， 6669. 21 22 附錄： 81 80 79 78 76 75 74 73 72 71 70 69 68 67 66 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 38 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 6 5 4 3 2 圖 5－ 1 某地區(qū)配網(wǎng)系統(tǒng)圖 the distribution system drawing 23 145 137 138 113 33 129 28 128 36 127 75 126 35 74 73 31 125 32 124 123 72 122 26 121 71 24 120 30 119 27 117 22 21 68 116 115 23 114 16 19 118 20 15 29 70 69 25 112 111 67 17 66 18 110 108 14 106 11 105 13 12 104 103 10 102 101 61 51 47 50 44 79 40 139 76 41 42 54 130 131 38 132 53 43 134 136 135 78 140 34 133 46 144 80 143 49 142 48 141 45 146 147 148 52 55 56 57 149 150 151 59 58 152 153 154 81 60 156 155 158 63 157 62 64 107 109 2 The topology after simplification 圖5－2 簡化后網(wǎng)絡(luò)拓?fù)鋱D 表示開關(guān)閉合 表示開關(guān)斷開 24 序號 節(jié)點(diǎn)（ From） 節(jié)點(diǎn)（ To） R（Ω） X（Ω） B（μ F） 1 1 101 2 101 10 3 10 102 4 102 66 5 66 103 6 103 13 7 66 104 8 104 12 9 1 105 10 105 11 11 11 106 12 106 14 13 1 107 14 107 14 15 1 108 16 108 18 17 18 109 18 109 13 19 1 110 20 110 17 21 17 67 22 67 111 23 111 12 24 67 112 25 112 15 26 1 70 27 70 25 28 25 69 29 69 113 30 113 15 31 1 16 32 16 118 33 118 19 34 1 114 35 114 23 36 23 115 37 115 22 38 22 116 39 116 68 40 68 21 41 21 117 42 127 20 43 1 28 44 28 129 45 129 33 46 1 119 47 119 30 48 30 120 49 120 24 50 24 19 51 1 71 52 71 123 53 123 31 25 54 71 121 55 121 26 56 26 122 57 122 72 58 72 27 59 72 124 60 124 32 61 1 125 62 125 73 63 73 31 64 73 74 65 74 32 66 74 35 67 35 126 68 126 75 69 75 33 70 75 127 71 127 36 72 36 128 73 128 29 74 1 130 75 130 38 76 38 131 77 131 43 78 1 45 79 45 146 80 146 47 81 1 134 82 134 78 83 78 135 84 135 42 85 42 136 86 136 46 87 78 41 88 41 76 89 76 138 90 138 40 91 76 137 92 137 44 93 1 139 94 139 34 95 34 140 96 140 40 97 1 133 98 133 46 99 1 141 100 141 48 101 48 142 102 142 79 103 79 44 104 79 49 105 49 143 106 143 80 107 80 47 26 108 80 144 109 144 50 110 50 145 111 145 51 112 1 57 113 57 149 114 149 58 115 1 152 116 152 60 117 1 54 118 54 53 119 53 132 120 132 43 121 1 56 122 56 147 123 147 55 124 55 52 125 52 148 126 148 51 127 1 150 128 150 59 129 59 151 130 151 58 131 1 62 132 62 157 133 157 63 134 63 158 135 158 64 136 1 153 137 153 81 138 81 154 139 154 61 140 61 155 141 155 60 142 81 156 143 156 64 27 程序原文： function [MVAbase, bus, gen, branch, success, et] = runpf(casename, mpopt, fname, solvedcase) %RUNPF Runs a power flow. % % [baseMVA, bus, gen, branch, success, et] = ... % runpf(casename, mpopt, fname, solvedcase) % % Runs a power flow (full AC Newton39。s method by default) and optionally % returns the solved values in the data matrices, the objective function % value, a flag which is true if the algorithm was successful in finding a % solution, and the elapsed time in seconds. All input arguments are % optional. If casename is provided it specifies the name of the input % data file or struct (see also 39。help caseformat39。 and 39。help loadcase39。) % containing the power flow data. The default value is 39。case939。. If the % mpopt is provided it overrides the default MATPOWER options vector and % can be used to specify the solution algorithm and output options among % other things (see 39。help mpoption39。 for details). If the 3rd argument is % given the pretty printed output will be appended to the file whose name % is given in fname. If solvedcase is specified the solved case will be % written to a case file in MATPOWER format with the specified name. If % solvedcase ends with 39。.mat39。 it saves the case as a MATfile otherwise it % saves it as an Mfile. % MATPOWER % $Id: ,v 2021/08/23 20:59:36 ray Exp $ % by Ray Zimmerman, PSERC Cornell % Copyright (c) 19962021 by Power System Engineering Research Center (PSERC) 28 % See for more info. %% initialize %% define named indices into bus, gen, branch matrices a10(1)=000000。 for i=2:64 digits(6) a10(i)=a10(i1)+1。 end a2=dec2bin(a10)。 k=[a2]。 for i=1:64 digits(6) baseMVA=10。 %% bus data % bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin bus = [ 1 3 0 0 0 0 1 1 0 10 1 。 10 1 0 0 0 0 1 1 0 10 1 。 11 1 0 0 0 0 1 1 0 10 1 。 12 1 0 0 0 0 1 1 0 10 1 。 13 1 0 0