【机会约束、鲁棒优化】具有排放感知型经济调度中机会约束和鲁棒优化研究【IEEE6节点、IEEE118节点算例】（Matlab代码实现）
欢迎来到本博客❤️❤️博主优势博客内容尽量做到思维缜密逻辑清晰为了方便读者。完整资源、论文复现、期刊合作、论文辅导及科研仿真定制事宜点击本文完整资源下载⛳️座右铭行百里者半于九十。⛳️赠与读者做科研涉及到一个深在的思想系统需要科研者逻辑缜密踏实认真但是不能只是努力很多时候借力比努力更重要然后还要有仰望星空的创新点和启发点。建议读者按目录次序逐一浏览免得骤然跌入幽暗的迷宫找不到来时的路它不足为你揭示全部问题的答案但若能解答你胸中升起的一朵朵疑云也未尝不会酿成晚霞斑斓的别一番景致万一它给你带来了一场精神世界的苦雨那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。或许雨过云收神驰的天地更清朗.......第一部分——内容介绍机会约束规划(chance constrained programming简称CCP)也称为概率约束规划,最初由Charnes和 Cooper[1在1958年提出,是针对约束条件中含有随机参数并在一定概率条件满足下的最优化问题.近几十年,机会约束广泛应用于参数不确定的优化和决策问题中.比如,网络优化2]、供应链管理3、生产计划[4和水流箱管理5等关于CCP的发展情形,感兴趣的读者可参见综述文献[6].机会约束规划的一般形式为:常规机会约束假设已知随机变量概率分布Р的准确信息.然而,在许多实际应用中随机变量∈ 难以精确预测,仅能够获得其部分信息(如矩和样本信息等).满足这些概率信息的分布构从伏)方法结合.建立在分布不确定集合最环情优L好分布鲁棒优化问题优化思想和随机优化方法结合,建立在分布不确定鲁棒优化(distributionally robust optimization,简称结合机会约束即为分布鲁棒机会约束优化问题(distributionally robust chance constrainedoptimization program,简称DRCC).分布鲁棒机会约束优化的一般形式为:调度者存在多元化的风险偏好导致对“最恶劣场景”的评估产生差异。如图 3 所示日前预测总偏差由偏离时长和对应时刻的偏离度决定。也就是说控制调度周期内不确定变量的总体偏离度就能够度量调度者风险偏好。但是在给定的总偏离度条件下学界普遍认为“最恶劣场景”发生在不确定域边界这在线性系统中必然成立但是在非线性系统中“最恶劣场景”可能发生在不确定域内。为此本文通过有限离散化不确定域生成梯次偏离的不确定域边界扩充“恶劣场景集”。在此基础上控制不确定变量的总体偏离值表征调度者的保守程度。考虑风电不确定性的电力系统调度模型中的部分输入量Pw和Pg具有不确定性通过传统数学规划方法得到的优化调度方案可能无法实现系统的经济性最优甚至可能违反系统的安全运行约束。针对传统的数学规划方法难以消除不确定参量对优化模型的负面影响Soyster提出了鲁棒优化方法。RO在优化模型建立之初就对其中的不确定参量进行分析明确不确定参量的波动区间。RO旨在得到这样一个优化方案:在不确定参量可能的取值范围内优化模型的所有约束条件都能被满足且使得优化模型目标函数能取得最恶劣情况下的最优值。因此RO所得的优化结果对不确定参量的波动具有免疫能力并且不失良好性能。对包含不确定参量的优化问题而言RO提出了一种“劣甲选优的处理刀八关注点在于最恶劣情况下模型的最优解仅需知道不确定参量波动的边界信息即可进行决策计算效率高需要了解的不确定参量信息少优化方案的抗干扰能力强因此具有很高的实用价值。RO的提出和发展弥补了传统优化方法的不足己成为处理含有不确定参量优化问题的重要工具。传统的数学优化模型一般形式如下:第二部分——运行结果2.1 风能部分代码figureyyaxis leftbar(x,y,FaceColor,[0.3010 0.7450 0.9330],EdgeColor,[0 0.4470 0.7410]);box offxlim([mu-3*sigma,mu3*sigma])xlabel(Error,FontSize,13.2,FontName,Times New Roman,FontWeight,Bold)ylabel(Number of samples,FontSize,13.2,FontName,Times New Roman,FontWeight,Bold)% a2axes;yyaxis rightfpfplot((x)normpdf(x,mu,sigma),[mu-3*sigma,mu3*sigma]);xlim([mu-3*sigma,mu3*sigma])fp.LineStyle -.;fp.LineWidth 2;fp.Color [0.9290 0.6940 0.1250];%set(a2,box,off,yaxislocation,right,color,none)% set(gca,xticklabel,[]);ylabel(Probability density,FontSize,13.2,FontName,Times New Roman,FontWeight,Bold)set(gcf,Position,[100 100 400 500]);set(gca,FontSize,12,FontName,Times New Roman)ax gca;ax.YAxis(1).Exponent 4;2.2 IEEE6和IEEE118x[1,2,3,4];bBarPlotBreak(x,y_time_record,15,20,Line,0.99);x_label{CC (Gaussian),DRO (Binomial),Clairvoyant,SO (Scenario)};set(gca,XTickLabel, x_label)legend(6 bus,118 bus)ylabel(time (s),FontSize,13.2,FontName,Times New Roman,FontWeight,Bold)set(gca,FontSize,12,FontName,Times New Roman)b(1).FaceColor[0.9290, 0.6940, 0.1250];b(2).FaceColor[0.4660, 0.6740, 0.1880];2.3 IEEE62.4 IEEE118部分代码function[fval_avg,x_avg]gen_SSO_average_performance(c_level,epsilon,T,N,M,bb,d_f,p,q,H,fmax,gmin,gmax,w_loc,w_num,error_data,ramp_rate,DR,UR)% scenarios to be generated to meet the probability guaranteen_dv4*N*T;Numceil(n_dv/(epsilon*c_level))-1;% generating wind scenarioindexceil(T*Num*(1-epsilon));d_reald_f;virtual_bbbb;[gm] gen_samples(error_data);for k1:w_numlocw_loc(k);% random_numzeros(T*Num,1);% for i1:Num% rng(default);% [random_T,~]random(gm,T);% random_num((i-1)*T1:i*T)random_T;% endrng(default);random_numrandom(gm,T*Num);random_numsort(random_num);virtual_bb(loc)max([bb(loc)-random_num(index),0]);end[standard_delta,real_bb] gen_standard_delta(bb,virtual_bb,w_loc);umin-real_bb*ramp_rate;umaxreal_bb*ramp_rate;[x_avg,fval_avg]MinC(T,N,M,real_bb,d_real,p,q,H,fmax,gmin,gmax,umin,umax,DR,UR);endfunction[wind_error] gen_wind_data(real_capacity_scale,filename)%filename.\data\WindGenTotalLoadYTD_2020.xls;input_data xlsread(filename, 1, B25:C52428);input_scalemean(input_data(:,1));wind_datainput_data./input_scale*real_capacity_scale;wind_errorwind_data(:,2)-wind_data(:,1);% error_normnormalize(error_data,scale);% wind_errorerror_norm*ratio;endfunction[wind_error] gen_wind_data(real_capacity_scale,filename)%filename.\data\WindGenTotalLoadYTD_2020.xls;input_data xlsread(filename, 1, B25:C52428);input_scalemean(input_data(:,1));wind_datainput_data./input_scale*real_capacity_scale;wind_errorwind_data(:,2)-wind_data(:,1);% error_normnormalize(error_data,scale);% wind_errorerror_norm*ratio;end第三部分——参考文献文章中一些内容引自网络会注明出处或引用为参考文献难免有未尽之处如有不妥请随时联系删除。(文章内容仅供参考具体效果以运行结果为准)[1]耿晓路,童小娇.分布鲁棒机会约束优化问题的研究[J].运筹学学报,2020,24(01):115-130.DOI:10.15960/j.cnki.issn.1007-6093.2020.01.009.[2]随权,林湘宁,童宁,李雪松,汪致洵,胡志豪,李正天,孙士茼.基于改进两阶段鲁棒优化的主动配电网经济调度[J].中国电机工程学报,2020,40(07):2166-21792396.DOI:10.13334/j.0258-8013.pcsee.182259.[3]程凤璐. 在线经济调度的鲁棒优化方法研究[D].山东大学,2015.[4]王晨曦. 含大规模风电的电力系统鲁棒优化调度研究[D].华南理工大学,2019.DOI:10.27151/d.cnki.ghnlu.2019.002201.[5]许书伟,吴文传,朱涛,王珍意.机会约束随机动态经济调度的凸松弛迭代求解法[J].电力系统自动化,2020,44(17):43-51.第四部分——本文完整资源下载资料获取更多粉丝福利MATLAB|Simulink|Python|数据|文档等完整资源获取本文完整资源下载