蒙特卡洛随机潮流与和密度估计Matlab程序

程序实现了基于蒙特卡洛的含源配电网随机潮流计算，并通过平滑核密度函数重构了状态变量的分布情况。

function mpc = IEEE33
%% MATPOWER Case Format : Version 2
mpc.version = ‘2‘;

%%-----  Power Flow Data  -----%%
%% system MVA base
mpc.baseMVA = 10;

%% bus data
%	bus_i	type	Pd	Qd	Gs	Bs	area	Vm	Va	baseKV	zone	Vmax	Vmin
mpc.bus = [
	1	3	0	0	0       0	2	1.05	0	12.66	1	1.07	0.90;
	2	1	0.10	0.06	0	0	2	1.0 0	12.66	1	1.07	0.90;
	3	1	0.09	0.04	0	0	2	1.0	0	12.66	1	1.07	0.90;
	4	1	0.12	0.08	0	0	1	1.0	0	12.66	1	1.07	0.90;
	5	1	0.06	0.03	0	0	1	1.0	0	12.66	1	1.07	0.90;
	6	1	0.06	0.02	0	0	1	1.0	0	12.66	1	1.07	0.90;
	7	1	0.20	0.10	0	0	1	1.0	0	12.66	1	1.07	0.90;
	8	1	0.20	0.10	0	0	1	1.0	0	12.66	1	1.07	0.90;
	9	1	0.06	0.02	0	0	1	1.0	0	12.66	1	1.07	0.90;
	10	1	0.06	0.02	0	0	1	1.0	0	12.66	1	1.07	0.90;
	11	1	0.045	0.03	0	0	1	1.0	0	12.66	1	1.07	0.90;
	12	1	0.06	0.035	0	0	1	1.0	0	12.66	1	1.07	0.90;
	13	1	0.06	0.035	0	0	1	1.0	0	12.66	1	1.07	0.90;
	14	1	0.12	0.08	0	0	1	1.0	0	12.66	1	1.07	0.90;
	15	1	0.06	0.01	0	0	3	1.0	0	12.66	1	1.07	0.90;
	16	1	0.06	0.02	0	0	3	1.0	0	12.66	1	1.07	0.90;
	17	1	0.06	0.02	0	0	2	1.0	0	12.66	1	1.07	0.90;
	18	1	0.09	0.04	0	0	2	1.0	0	12.66	1	1.07	0.90;
	19	1	0.09	0.04	0	0	3	1.0	0	12.66	1	1.07	0.90;
	20	1	0.09	0.04	0	0	3	1.0	0	12.66	1	1.07	0.90;
	21	1	0.09	0.04	0	0	3	1.0	0	12.66	1	1.07	0.90;
	22	1	0.09	0.04	0	0	3	1.0	0	12.66	1	1.07	0.90;
	23	1	0.09	0.05	0	0	3	1.0	0	12.66	1	1.07	0.90;
	24	1	0.42	0.20	0	0	3	1.0	0	12.66	1	1.07	0.90;
	25	1	0.42	0.20	0	0	2	1.0	0	12.66	1	1.07	0.90;
	26	1	0.06	0.025	0	0	2	1.0	0	12.66	1	1.07	0.90;
	27	1	0.06	0.025	0	0	2	1.0	0	12.66	1	1.07	0.90;
	28	1	0.06	0.020	0	0	3	1.0	0	12.66	1	1.07	0.90;
	29	1	0.12	0.070	0	0	3	1.0	0	12.66	1	1.07	0.90;
	30	1	0.20	0.60	0	0	2	1.0	0	12.66	1	1.07	0.90;
	31	1	0.15	0.07	0	0	1	1.0	0	12.66	1	1.07	0.90;
	32	1	0.21	0.10	0	0	1	1.0	0	12.66	1	1.07	0.90;
	33	1	0.06	0.04	0	0	3	1.0	0	12.66	1	1.07	0.90;
];

%% generator data
%	bus	Pg	Qg	Qmax	Qmin	Vg	mbase	status	Pmax	Pmin	Pc1	Pc2	Qc1min	Qc1max	Qc2min	Qc2max	ramp_agc	ramp_10	ramp_30	ramp_q	apf
mpc.gen = [
	1	0	0	500	-500	1.05	10	1	1000	-1000	0	0	0	0	0	0	0	0	0	0	0;
    17	0	0	500	-500	1.00	10	1	1000	-1000	0	0	0	0	0	0	0	0	0	0	0;
];

%% branch data
%	fbus	tbus	r	x	b	rateA	rateB	rateC	ratio	angle	status	angmin	angmax
mpc.branch = [
	1	2	0.00575	0.00293	0	600	600	600	0	0	1	-360	360;
	2	3	0.03076	0.01567	0	900	900 900	0	0	1	-360	360;
	3	4	0.02284	0.01163	0	500	500	500	0	0	1	-360	360;
	4	5	0.02378	0.01211	0	500	500	500	0	0	1	-360	360;
	5	6	0.05110	0.04411	0	500	500	500	0	0	1	-360	360;
	6	7	0.01168	0.03861	0	500	500	500	0	0	1	-360	360;
	7	8	0.04439	0.01467	0	600	600	600	0	0	1	-360	360;
	8	9	0.06426	0.04617	0	500	500	500	0	0	1	-360	360;
	9	10	0.06514	0.04617	0	900	900	900	0	0	1	-360	360;
	10  11	0.01227	0.00406	0	900	900	900	0	0	1	-360	360;
	11	12	0.02336	0.00772	0	900	900	900	0	0	1	-360	360;
	12	13	0.09159	0.07206	0	480	480	480	0	0	1	-360	360;
	13	14	0.03379	0.04448	0	900	900	900	0	0	1	-360	360;
	14	15	0.03687	0.03282	0	900	900	900	0	0	1	-360	360;
	15	16	0.04656	0.03400	0	900	900	900	0	0	

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3935  2017-11-21 20:52  蒙特卡洛随机潮流与核密度估计Matlab程序\IEEE33.m
     文件        6529  2017-11-21 21:15  蒙特卡洛随机潮流与核密度估计Matlab程序\MCPPF_ksdensity.m
     目录           0  2017-11-21 21:16  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\
     目录           0  2017-11-21 21:16  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\
     文件        4745  2010-06-09 22:56  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\add_constraints.m
     文件        5694  2010-04-27 03:45  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\add_costs.m
     文件        3076  2010-04-27 03:45  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\add_vars.m
     文件        4368  2011-12-02 01:04  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\build_cost_params.m
     文件        4061  2011-12-02 01:04  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\compute_cost.m
     文件        3959  2010-04-27 03:45  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\display.m
     文件        1657  2010-04-27 03:45  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\get.m
     文件        2740  2011-12-02 01:03  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\get_cost_params.m
     文件        3012  2010-06-09 22:56  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\get_idx.m
     文件        1450  2010-04-27 03:45  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\get_mpc.m
     文件        2405  2010-06-09 22:56  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\getN.m
     文件        2221  2010-04-27 03:45  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\getv.m
     文件        2945  2010-04-27 03:45  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\linear_constraints.m
     文件        8294  2010-06-09 22:56  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\opf_model.m
     文件        2149  2010-05-05 02:44  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\@opf_model\userdata.m
     文件        6175  2010-04-27 03:45  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\add_userfcn.m
     文件        2998  2010-04-27 03:45  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\bustypes.m
     文件        7587  2010-03-12 01:06  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case_ieee30.m
     文件        8768  2015-12-02 14:32  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case_ieee33.m
     文件       32633  2010-03-11 02:08  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case118.m
     文件        4521  2010-03-11 02:08  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case14.m
     文件      340794  2010-03-11 02:08  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case2383wp.m
     文件        8708  2010-03-11 02:08  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case24_ieee_rts.m
     文件      717636  2010-03-11 02:08  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case2736sp.m
     文件      715220  2010-03-11 02:08  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case2737sop.m
     文件      732898  2010-03-11 02:08  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case2746wop.m
     文件      736274  2010-03-11 02:08  蒙特卡洛随机潮流与核密度估计Matlab程序\Modified_matpower\case2746wp.m
............此处省略248个文件信息