电力系统小干扰分析

电力系统发生小扰动时的小干扰稳定程序

%function [baseMVA bus gen branch areas gencost] = case39
%CASE39    Power flow data for 39 bus case.
%   Please see ‘help caseformat‘ for details on the case file format.
%
%   based on data from:
%
%       M. A. Pai Energy Function Analysis for Power System Stability
%       Kluwer Academic Publishers Boston 1989.

%   MATPOWER
%   $Id: case39.mv 1.6 2004/09/21 01:48:33 ray Exp $

%%-----  Power Flow Data  -----%%
%% system MVA base
baseMVA = 100;

%% 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.04	0       345	1	1.1	0.9
2	2	0	0	0	0	1	1.025	9.28	345	1	1.1	0.9
3	2	0	0	0	0	1	1.025	4.6648	345	1	1.1	0.9
4	1	0	0	0	0	1	1.0258	-2.2168	345	1	1.1	0.9
5	1	125	50	0	0	1	0.99563	-3.9888	345	1	1.1	0.9
6	1	90	30	0	0	1	1.0127	-3.6874	345	1	1.1	0.9
7	1	0	0	0	0	1	1.0258	3.7197	345	1	1.1	0.9
8	1	100	35	0	0	1	1.0159	0.72754	345	1	1.1	0.9
9	1	0	0	0	0	1	1.0324	1.9667	345	1	1.1	0.9

];

%% generator data
%	bus	Pg	Qg	Qmax	Qmin	Vg	mbase	status	Pmax	Pmin
gen0 = [
1	71.641	27.046	0	0	0	100	0	0	0
2	163	6.654	0	0	0	100 	0	0	0
3	85	-10.86	0	0	0	100	0	0	0
];
gen=[
% gen	bus	Tj	Ra	Xd	X‘d	Xq	X‘q	X‘‘d	T‘d0	T‘‘d0	T‘q0	T‘‘q0	D	xl
1	1	6.4     0	1.8     0.3     1.75	0.47	0.23	4.8     3.50E-02	1.5     7.00E-02	0	0.15
2	2	4.62	0	1.53	0.29	1.51	0.29	0.17	4       0.029       1.61    0.151       0.9	0.153
3	3	3.16	0	2.2     0.23	2       0.23	0.16	7.5     0.04        3.5     0.04        0.9	0.09
];


%% branch data
%	fbus	tbus	r	x	b	rateA	rateB	rateC	ratio	angle	status
branch = [
4	5	0.01	0.085	0.176	0	0	0	0	0	1
4	6	0.017	0.092	0.158	0	0	0	0	0	1
5	7	0.032	0.161	0.306	0	0	0	0	0	1
6	9	0.039	0.17	0.358	0	0	0	0	0	1
7	8	0.0085	0.072	0.149	0	0	0	0	0	1
8	9	0.0119	0.1008	0.209	0	0	0	0	0	1
1	4	0	0.0576	0	0	0	0	0	0	1
2	7	0	0.0625	0	0	0	0	0	0	1
3	9	0	0.0586	0	0	0	0	0	0	1

];

%%-----  OPF Data  -----%%
%% area data
areas = [
	1	1;
];

%% generator cost data
%	1	startup	shutdown	n	x0	y0	...	xn	yn
%	2	startup	shutdown	n	c（n-1）	...	c0
gencost = [
	2	0	0	3	0.01	0.3	0.2;
	2	0	0	3	0.01	0.3	0.2;
	2	0	0	3	0.01	0.3	0.2;
	2	0	0	3	0.01	0.3	0.2;
	2	0	0	3	0.01	0.3	0.2;
	2	0	0	3	0.01	0.3	0.2;
	2	0	0	3	0.01	0.3	0.2;
	2	0	0	3	0.01	0.3	0.2;
	2	0	0	3	0.006	0.3	0.2;
	2	0	0	3	0.006	0.3	0.2;
];
Exc = [ ...
%1Gen  2type 3VRmax 4VRmin 5KA 6TA 7KF 8TF 9no use 10TE 11TR 12Ae 13Be
1 2 5 -5 50 0.2 0.063 0.35 0.01 0.314 0.001 0.0039 1.555;
2 2 5 -5 50 0.2 0.063 0.35 0.01 0.314 0.001 0.0039 1.555;
3 2 5 -5 50 0.2 0.063 0.35 0.01 0.314 0.001 0.0039 1.555 ];

Gov = [ ...
%gen	type	ωref	R	Tmax	Tmin	T2	T1
 1 2 1 0.1 1.0 0.1 0.1 0.3;  
 2 2 1 0.1 1.0 0.1 0.1 0.3;
 3 2 1 0.1 1.0 0.1 0.1 0.3];
Pss = [ ...
%gen	type	KW	TW	TS1	TS2	TS3	TS4 TS5
    1    2    5   5    0.30    0.06    0.20    0.04    1e-4;
    ];

fault=[
    450.10.12;% fault 1
    460.10.12;% fault 2
    570.10.12;% fault 3
    690.10.12;% fault 4
    780.10.12;% fault 5
    890.10.12;% fault 6

    450.10.14;% fault 7
    460.10.14;% fault 8
    570.10.14;% fault 9
    690.10.14;% fault

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5829  2016-01-07 17:40  小干扰分析\case9.m

     文件       6128  2016-01-07 17:40  小干扰分析\compute_t0.m

     文件        197  2016-01-07 17:40  小干扰分析\convert_ID.m

     文件       3585  2016-01-07 17:40  小干扰分析\fm_call.m

     文件       6770  2016-01-07 17:40  小干扰分析\getJ.m

     文件       2352  2016-01-07 17:40  小干扰分析\getVI.m

     文件      90305  2016-01-07 17:40  小干扰分析\grid_simulation1.mat

     文件      92702  2016-01-07 17:40  小干扰分析\grid_simulation10.mat

     文件      97726  2016-01-07 17:40  小干扰分析\grid_simulation100.mat

     文件      97458  2016-01-07 17:40  小干扰分析\grid_simulation101.mat

     文件      96050  2016-01-07 17:40  小干扰分析\grid_simulation102.mat

     文件      96468  2016-01-07 17:40  小干扰分析\grid_simulation103.mat

     文件      95103  2016-01-07 17:40  小干扰分析\grid_simulation104.mat

     文件      98427  2016-01-07 17:40  小干扰分析\grid_simulation105.mat

     文件      97885  2016-01-07 17:40  小干扰分析\grid_simulation106.mat

     文件      97640  2016-01-07 17:40  小干扰分析\grid_simulation107.mat

     文件      96144  2016-01-07 17:40  小干扰分析\grid_simulation108.mat

     文件      92127  2016-01-07 17:40  小干扰分析\grid_simulation11.mat

     文件      91513  2016-01-07 17:40  小干扰分析\grid_simulation12.mat

     文件      92244  2016-01-07 17:40  小干扰分析\grid_simulation13.mat

     文件      91417  2016-01-07 17:40  小干扰分析\grid_simulation14.mat

     文件      93914  2016-01-07 17:40  小干扰分析\grid_simulation15.mat

     文件      93527  2016-01-07 17:40  小干扰分析\grid_simulation16.mat

     文件      92972  2016-01-07 17:40  小干扰分析\grid_simulation17.mat

     文件      92233  2016-01-07 17:40  小干扰分析\grid_simulation18.mat

     文件      92829  2016-01-07 17:40  小干扰分析\grid_simulation19.mat

     文件      89858  2016-01-07 17:40  小干扰分析\grid_simulation2.mat

     文件      91979  2016-01-07 17:40  小干扰分析\grid_simulation20.mat

     文件      94675  2016-01-07 17:40  小干扰分析\grid_simulation21.mat

     文件      94276  2016-01-07 17:40  小干扰分析\grid_simulation22.mat

............此处省略96个文件信息