电力系统 matlab程序,电力系统稳定性分析matlab程序

电力系统稳定性分析作业一

1

euler.m ,reuler.m, kunta.m分别为(1)中的欧拉法，改进欧拉法，龙格库塔法的主程序；doty.m,doty2.m,doty3.m均为(1)中子函数程序。Runge-Kutta.m为(2)和(3)的运行程序。

下表为三种方法的部分运行结果功角数据：

时

间

0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08

欧拉35.161

5

35.161

5

35.282

8

35.523

6

35.882

36.356

36.943

4

37.642

2

38.450

改进35.161

5

35.222

2

35.402

3

35.699

9

36.113

36.639

4

37.277

38.023

4

38.876

3

龙格35.161

5

35.221

9

35.401

6

35.698

9

36.111

6

36.637

6

37.274

7

38.020

7

38.873

1

时

间

0．09 0.10 0.11 0.12 0.13 0.14 0.15 0.16 0.17

欧拉39.364

6

40.3835 41.504

3

42.636

6

43.777

3

44.923

46.070

9

37.642

2

38.450

改进39.833

2

40．8918 42.005

1

43.125

6

44.250

1

45.375

7

46.499

5

47.618

7

48.730

5

龙格39.829

5

40.8875 42.000

2

43.120

44.244

45.369

46.492

3

47.611

48.722

4

(1)欧拉法

在matlaB中输入命令[t,x,y,z]=euler('doty','doty2','doty3',0,5,0.1,0.01)可得

t-w曲线，t-δ曲线分别如下图所示。具体功角，角速度数据分别见文件1.mat 和2.mat