培训营销型网站建设,怎样在微信上开店卖东西,东莞南城做网站,昆明设计网站建设#x1f4a5;#x1f4a5;#x1f49e;#x1f49e;欢迎来到本博客❤️❤️#x1f4a5;#x1f4a5; #x1f3c6;博主优势#xff1a;#x1f31e;#x1f31e;#x1f31e;博客内容尽量做到思维缜密#xff0c;逻辑清晰#xff0c;为了方便读者。 ⛳️座右铭欢迎来到本博客❤️❤️ 博主优势博客内容尽量做到思维缜密逻辑清晰为了方便读者。 ⛳️座右铭行百里者半于九十。 本文目录如下 目录 1 概述 2 运行结果 2.1 风能 2.2 IEEE6和IEEE118 2.3 IEEE6 2.4 IEEE118 3 参考文献 4 Matlab代码及数据实现 1 概述
机会约束规划(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 运行结果
2.1 风能 部分代码 figure yyaxis left bar(x,y,FaceColor,[0.3010 0.7450 0.9330],EdgeColor,[0 0.4470 0.7410]);box off xlim([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 right fpfplot((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和IEEE118 x[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 IEEE6 2.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 guarantee n_dv4*N*T; Numceil(n_dv/(epsilon*c_level))-1; % generating wind scenario indexceil(T*Num*(1-epsilon)); d_reald_f; virtual_bbbb; [gm] gen_samples(error_data); for k1:w_num locw_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; % end rng(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); end
function[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
function[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
3 参考文献 部分理论来源于网络如有侵权请联系删除。 [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.
4 Matlab代码及数据实现
