我不能运行此代码,请修复它

3视图(30天)

显示旧的评论

Fajri Mardiansyah 2018年5月22日

0
链接

这个问题直接联系

https://de.mathworks.com/matlabcentral/answers/402018-i-can-t-run-this-code-please-fix-it

评论道: 为副总经理 2022年8月6日

                         如果真正的
                        
                         %这段代码从4号线开始
                        
                         结束
                        
                         N = 33;%的组件数量
                        
                         λ= (0.3979;0.8209;0.7666;0.1206;0.8577;…
                        
                         0.2586;0.7049;0.2429;0.3521;0.8118;…
                        
                         0.1850;0.5135;0.3920;0.7194;0.6776;…
                        
                         0.2065;0.3197;0.3226;0.8566;0.3348;…
                        
                         0.8950;0.8833;0.6400;0.0412;0.0518;…
                        
                         0.9458;0.2257;0.7303;0.2191;0.0101;…
                        
                         0.7205;0.1289;0.1327);
                        
                         MTTR = (0.3546;0.6970;0.8490;0.8724;0.0411;…
                        
                         0.2098;0.7382;0.7379;0.1978;0.4534;…
                        
                         0.2299;0.0704;0.3979;0.8555;0.6809;…
                        
                         0.2954;0.8536;0.7195;0.3405;0.0495;…
                        
                         0.0174;0.7846;0.2554;0.6597;0.8496;…
                        
                         0.2965;0.2238;0.0066;0.0684;0.4306;…
                        
                         0.7953;0.7759;0.9673);
                        
                         权力55 = [0 0 34 16 0 0 0 5 0 0 4 40 7 45 146 82…
                        
                         67 9 405 45 83 0 0 4 0 32 0 28 2 0 450);
                        
                         用户=装天花板(功率/ 2)1;
                        
                         % %失败历史
                        
                         时间= 1000;%年
                        
                         interuption = 0;
                        
                         outage_time = 0;
                        
                         maxT = 0;
                        
                         为i = 1: N
                        
                         downT(我),upT (i)] = failure_history2(λ(i), MTTR(我),持续时间);
                        
                         坏蛋= max (upT{我}(1:结束);
                        
                         如果(maxT < cur)
                        
                         maxT =坏蛋;
                        
                         结束
                        
                         interuption =中断+长度(downT{我})1;
                        
                         outage_time = outage_time +总和(upT{我}-downT{我});
                        
                         %图;
                        
                         线(1:长度(downT{我}),downT{我});
                        
                         线(1:长度(downT{我}),upT {},“颜色”,“r”,“标记”,“。”);
                        
                         结束
                        
                         %指标计算
                        
                         average_interuption = interuption /持续时间;
                        
                         average_outage_time = outage_time /持续时间;
                        
                         customers_num =总和(用户);
                        
                         total_power =总和(权力);
                        
                         SAIFI = average_interuption * customers_num / customers_num;
                        
                         SAIDI = average_outage_time * customers_num / customers_num;
                        
                         CAIDI = SAIDI / SAIFI;
                        
                         ASAI = (8760 - average_outage_time customers_num * * customers_num) / (customers_num * 8760);
                        
                         外加厚= average_outage_time * total_power;%千瓦*小时
                        
                         disp (“- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”);
                        
                         disp (“步骤1”);
                        
                         disp ([“SAIFI”num2str (SAIFI)]);
                        
                         disp ([“说”num2str(说)]);
                        
                         disp ([“CAIDI”num2str (CAIDI)]);
                        
                         disp ([“ASAI”num2str (ASAI)]);
                        
                         disp ([外加厚的num2str(外加厚)]);

19日的评论
显示隐藏 18年长的评论

为副总经理 2021年9月27日

编辑:为副总经理 2021年9月27日

就其价值而言,这就是从这里:

https://www.scirp.org/journal/PaperInformation.aspx?PaperID=78287

我编辑它修复的部分网站格式打破东西。它仍然无法运行没有缺失的功能。

“行”不是代码行号,但电气配电线路;同样,“酒吧”是母线,或上下文中的变量命名、开关。

至于如何重新安排这个在其他节点来模拟,我不知道。我必须用我的头在网络模型知道/元素的行,和我必须开发课程及研究文章。

                                清晰;
                               
                                N = 33;%的组件数量
                               
                                λ= (0.3979;0.8209;0.7666;0.1206;0.8577;…
                               
                                0.2586;0.7049;0.2429;0.3521;0.8118;…
                               
                                0.1850;0.5135;0.3920;0.7194;0.6776;…
                               
                                0.2065;0.3197;0.3226;0.8566;0.3348;…
                               
                                0.8950;0.8833;0.6400;0.0412;0.0518;…
                               
                                0.9458;0.2257;0.7303;0.2191;0.0101;…
                               
                                0.7205;0.1289;0.1327);
                               
                                MTTR = (0.3546;0.6970;0.8490;0.8724;0.0411;…
                               
                                0.2098;0.7382;0.7379;0.1978;0.4534;…
                               
                                0.2299;0.0704;0.3979;0.8555;0.6809;…
                               
                                0.2954;0.8536;0.7195;0.3405;0.0495;…
                               
                                0.0174;0.7846;0.2554;0.6597;0.8496;…
                               
                                0.2965;0.2238;0.0066;0.0684;0.4306;…
                               
                                0.7953;0.7759;0.9673);
                               
                                功率= 55 [0 0 34 16 0 0 0 5 0 0 4 40 7 45 15 146 82 67 9 405 45 83 0 0 4 0 32 0 28 2 0 450);
                               
                                用户=装天花板(功率/ 2)1;
                               
                                % %失败历史
                               
                                时间= 1000;%年
                               
                                interuption = 0;
                               
                                outage_time = 0;
                               
                                maxT = 0;
                               
                                为i = 1: N
                               
                                我[downT {}, upT{我}]= failure_history2(λ(i), MTTR(我),持续时间);
                               
                                坏蛋= max (upT{我}(1:结束);
                               
                                如果(maxT < cur)
                               
                                maxT =坏蛋;
                               
                                结束
                               
                                interuption =中断+长度(downT{我})1;
                               
                                outage_time = outage_time +总和(upT{我}-downT{我});
                               
                                %图;
                               
                                线(1:长度(downT{我}),downT{我});
                               
                                线(1:长度(downT{我}),upT {},“颜色”,“r”,“标记”,“。”);
                               
                                结束
                               
                                %指标计算
                               
                                average_interuption = interuption /持续时间;
                               
                                average_outage_time = outage_time /持续时间;
                               
                                customers_num =总和(用户);
                               
                                total_power =总和(权力);
                               
                                SAIFI = average_interuption * customers_num / customers_num;
                               
                                SAIDI = average_outage_time * customers_num / customers_num;
                               
                                CAIDI = SAIDI / SAIFI;
                               
                                ASAI = (8760 - average_outage_time customers_num * * customers_num) / (customers_num * 8760);
                               
                                外加厚= average_outage_time * total_power;%千瓦*小时
                               
                                disp (“- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”);
                               
                                disp (“步骤1”);
                               
                                disp ([“SAIFI”num2str (SAIFI)]);
                               
                                disp ([“说”num2str(说)]);
                               
                                disp ([“CAIDI”num2str (CAIDI)]);
                               
                                disp ([“ASAI”num2str (ASAI)]);
                               
                                disp ([外加厚的num2str(外加厚)]);
                               
                                % %步骤2开关
                               
                                switch_834_lines =(20 21日22日31日32);
                               
                                switch_832_lines =[17 18 19 30 23日24日25日switch_834_lines);
                               
                                switch_834_elements = [19 20 21 30 31];
                               
                                switch_832_elements = (switch_834_elements 18 22日23日24日25日29日32);
                               
                                interaptionXcustomer_number = 0;
                               
                                interaptionXcustomer_duration = 0;
                               
                                unservedEnergy = 0;
                               
                                cur_t = 0;
                               
                                而cur_t < maxT
                               
                                %找到下一个失败的时间
                               
                                min_t =(时间+ 1)* 24 * 365;
                               
                                fail_line_number = 0;
                               
                                为i = 1: N
                               
                                清晰的indx;
                               
                                indx =找到(downT{我}> cur_t);
                               
                                如果(~ isempty (indx))
                               
                                indx = indx (1);
                               
                                如果我(downT {} (indx) < min_t)
                               
                                min_t = downT{我}(indx);
                               
                                fail_line_number =我;
                               
                                n = indx;
                               
                                cur_outage_time = upT{我}(indx)——downT{我}(indx);
                               
                                结束
                               
                                结束
                               
                                结束
                               
                                cur_t = upT {fail_line_number} (n);
                               
                                in_switch_834 = ~ isempty(找到(switch_834_lines = = fail_line_number));
                               
                                customer_834 =总和(用户(switch_834_elements));
                               
                                power_834 =总和(权力(switch_834_elements));
                               
                                in_switch_832 = ~ isempty(找到(switch_832_elements = = fail_line_number));
                               
                                customer_832 =总和(用户(switch_832_elements));
                               
                                power_832 =总和(权力(switch_832_elements));
                               
                                % interaption计算
                               
                                如果(~ in_switch_832)
                               
                                interaptionXcustomer_number = interaptionXcustomer_number + customers_num;
                               
                                interaptionXcustomer_duration = interaptionXcustomer_duration + customers_num * cur_outage_time;
                               
                                unservedEnergy = unservedEnergy + total_power * cur_outage_time;
                               
                                其他的
                               
                                如果(in_switch_834)
                               
                                interaptionXcustomer_number = interaptionXcustomer_number + customer_834;
                               
                                interaptionXcustomer_duration = interaptionXcustomer_duration + customer_834 * cur_outage_time;
                               
                                unservedEnergy = unservedEnergy + power_834 * cur_outage_time;
                               
                                其他的
                               
                                interaptionXcustomer_number = interaptionXcustomer_number + customer_832;
                               
                                interaptionXcustomer_duration = interaptionXcustomer_duration + customer_832 * cur_outage_time;
                               
                                unservedEnergy = unservedEnergy + power_832 * cur_outage_time;
                               
                                结束
                               
                                结束
                               
                                结束
                               
                                %指标计算
                               
                                SAIFI = interaptionXcustomer_number / customers_num /持续时间;
                               
                                SAIDI = interaptionXcustomer_duration / customers_num /持续时间;
                               
                                CAIDI = SAIDI / SAIFI;
                               
                                ASAI = (customers_num * 8760 - interaptionXcustomer_duration /持续时间)/ (customers_num * 8760);
                               
                                外加厚= unservedEnergy /持续时间;%千瓦*小时
                               
                                disp (“- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”);
                               
                                disp (“步骤2”);
                               
                                disp ([“SAIFI”num2str (SAIFI)]);
                               
                                disp ([“说”num2str(说)]);
                               
                                disp ([“CAIDI”num2str (CAIDI)]);
                               
                                disp ([“ASAI”num2str (ASAI)]);
                               
                                disp ([外加厚的num2str(外加厚)]);
                               
                                % %步骤3开关从步骤2和DGs
                               
                                switch_834_lines =(20 21日22日31日32);
                               
                                switch_832_lines =[17 18 19 30 23日24日25日switch_834_lines);
                               
                                switch_834_elements = [19 20 21 30 31];
                               
                                switch_832_elements = (switch_834_elements 18 22日23日24日25日29日32);
                               
                                interaptionXcustomer_number = 0;
                               
                                interaptionXcustomer_duration = 0;
                               
                                unservedEnergy = 0;
                               
                                cur_t = 0;
                               
                                而cur_t < maxT
                               
                                %找到下一个失败的时间
                               
                                min_t =(时间+ 1)* 24 * 365;
                               
                                fail_line_number = 0;
                               
                                为i = 1: N
                               
                                清晰的indx;
                               
                                indx =找到(downT{我}> cur_t);
                               
                                如果(~ isempty (indx))
                               
                                indx = indx (1);
                               
                                如果我(downT {} (indx) < min_t)
                               
                                min_t = downT{我}(indx);
                               
                                fail_line_number =我;
                               
                                n = indx;
                               
                                cur_outage_time = upT{我}(indx)——downT{我}(indx);
                               
                                结束
                               
                                结束
                               
                                结束
                               
                                cur_t = upT {fail_line_number} (n);
                               
                                in_switch_834 = ~ isempty(找到(switch_834_lines = = fail_line_number));
                               
                                customer_834 =总和(用户(switch_834_elements));
                               
                                power_834 =总和(权力(switch_834_elements));
                               
                                in_switch_832 = ~ isempty(找到(switch_832_elements = = fail_line_number));
                               
                                customer_832 =总和(用户(switch_832_elements));
                               
                                power_832 =总和(权力(switch_832_elements));
                               
                                customer_800 = customers_num-customer_832;
                               
                                power_800 = total_power-power_832;
                               
                                % interaption计算
                               
                                如果(~ in_switch_832)
                               
                                interaptionXcustomer_number = interaptionXcustomer_number + customer_800;
                               
                                interaptionXcustomer_duration = interaptionXcustomer_duration + customer_800 * cur_outage_time;
                               
                                unservedEnergy = unservedEnergy + power_832 * cur_outage_time;
                               
                                其他的
                               
                                如果(in_switch_834)
                               
                                interaptionXcustomer_number = interaptionXcustomer_number + customer_834;
                               
                                interaptionXcustomer_duration = interaptionXcustomer_duration + customer_834 * cur_outage_time;
                               
                                unservedEnergy = unservedEnergy + power_834 * cur_outage_time;
                               
                                其他的
                               
                                interaptionXcustomer_number = interaptionXcustomer_number + customer_832;
                               
                                interaptionXcustomer_duration = interaptionXcustomer_duration + customer_832 * cur_outage_time;
                               
                                unservedEnergy = unservedEnergy + power_832 * cur_outage_time;
                               
                                结束
                               
                                结束
                               
                                结束
                               
                                %指标计算
                               
                                SAIFI = interaptionXcustomer_number / customers_num /持续时间;
                               
                                SAIDI = interaptionXcustomer_duration / customers_num /持续时间;
                               
                                CAIDI = SAIDI / SAIFI;
                               
                                ASAI = (customers_num * 8760 - interaptionXcustomer_duration /持续时间)/ (customers_num * 8760);
                               
                                外加厚= unservedEnergy /持续时间;%千瓦*小时
                               
                                disp (“- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -”);
                               
                                disp (“第三步”);
                               
                                disp ([“SAIFI”num2str (SAIFI)]);
                               
                                disp ([“说”num2str(说)]);
                               
                                disp ([“CAIDI”num2str (CAIDI)]);
                               
                                disp ([“ASAI”num2str (ASAI)]);
                               
                                disp ([外加厚的num2str(外加厚)]);