电力系统建模与分析软件：MATLAB二次开发_（4）.电力系统动态分析.docx

PAGE1

PAGE1

电力系统动态分析

1.电力系统动态模型概述

电力系统动态分析是研究电力系统在各种扰动下（如短路故障、发电机跳闸、负荷变化等）的行为和性能。动态模型包括发电机模型、负荷模型、输电线路模型、励磁系统模型、调速器模型等。这些模型在MATLAB中可以通过Simulink进行仿真，也可以通过MATLAB脚本进行数值计算。

1.1发电机模型

发电机模型是电力系统动态分析的核心。常见的发电机模型包括同步发电机模型和异步发电机模型。同步发电机模型通常包括电磁暂态模型和机械暂态模型。

1.1.1同步发电机电磁暂态模型

同步发电机的电磁暂态模型描述了发电机在电磁暂态过程中的电压和电流关系。常用的模型有次暂态模型和暂态模型。

次暂态模型：

电压方程：

$$

V_d=E_d’-(X_d’-X_0)I_q

$$

$$

V_q=E_q’-(X_q’-X_0)I_d

$$

励磁方程：

$$

=-

$$

$$

=-

$$

暂态模型：

电压方程：

$$

V_d=E_d-(X_d-X_0)I_q

$$

$$

V_q=E_q-(X_q-X_0)I_d

$$

励磁方程：

$$

=-

$$

$$

=-

$$

1.1.2同步发电机机械暂态模型

同步发电机的机械暂态模型描述了发电机转子的运动特性。主要方程为：

转子运动方程：

$$

J=T_m-T_e-D(-1)

$$

其中，J是转子惯性常数，ω是转子角速度，Tm是机械转矩，Te是电磁转矩，D

1.1.3发电机模型在MATLAB中的实现

在MATLAB中，可以使用Simulink搭建同步发电机的动态模型。以下是一个简单的次暂态模型的实现示例：

%同步发电机次暂态模型

functiondydt=generator_model(t,y,params)

%解析参数

E_d_prime=y(1);

E_q_prime=y(2);

I_d=y(3);

I_q=y(4);

omega=y(5);

V_d=params.V_d;

V_q=params.V_q;

X_d_prime=params.X_d_prime;

X_q_prime=params.X_q_prime;

X_0=params.X_0;

tau_d_prime=params.tau_d_prime;

tau_q_prime=params.tau_q_prime;

J=params.J;

T_m=params.T_m;

D=params.D;

%电磁暂态方程

dE_d_prime_dt=-(E_d_prime-V_d+(X_d_prime-X_0)*I_q)/tau_d_prime;

dE_q_prime_dt=-(E_q_prime-V_q+(X_q_prime-X_0)*I_d)/tau_q_prime;

%电压方程

V_d=E_d_prime-(X_d_prime-X_0)*I_q;

V_q=E_q_prime-(X_q_prime-X_0)*I_d;

%电磁转矩

T_e=1.5*(E_d_prime*I_q-E_q_prime*I_d);

%机械暂态方程

domega_dt=(T_m-T_e-D*(omega-1))/J;

%返回状态变量的导数

dydt=[dE_d_prime_dt;dE_q_prime_dt;0;0;domega_dt];

end

%初始化参数

params=struct(V_d,1.0,V_q,0.0,X_d_prime,1.0,X_q_prime,1.0,X_0,0.1,...

tau_d_prime,0.05,tau_q_prime,0.1,J,1.0,T_m,1.0,D,0.1);

%初始状态

y0=[1.0;0.0;0.0;0.0;1.0];

%时间范围

tspan=[010];

%求解微分方程

[t,