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

。 电力系统稳定性分析作业一 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 35.161 35.282 35.523 35.882 36.356 36.943 37.642 38.450 拉 5 5 8 6 0 0 4 2 0 改 35.161 35.222 35.402 35.699 36.113 36.639 37.277 38.023 38.876 进 5 2 3 9 0 4 0 4 3 龙 35.161 35.221 35.401 35.698 36.111 36.637 37.274 38.020 38.873 格 5 9 6 9 6 6 7 7 1 时 0．09 0.10 0.11 0.12 0.13 0.14 0.15 0.16 0.17 间 欧 39.364 40.3835 41.504 42.636 43.777 44.923 46.070 37.642 38.450 拉 6 3 6 3 0 9 2 0 改 39.833 40 ．8918 42.005 43.125 44.250 45.375 46.499 47.618 48.730 进 2 1 6 1 7 5 7 5 龙 39.829 40.8875 42.000 43.120 44.244 45.369 46.492 47.611 48.722 格 5 2 0 0 0 3 0 4 （1）欧拉法 在 matlaB 中输入命令 [t,x,y,z]=euler(doty,doty2,doty3,0,5,0.1,0.01) 可得 t-w 曲线， t- δ曲线分别如下图所示。具体功角，角速度数据分别见文件 1.mat 和 2.mat 精选资料，欢迎下载 。 （2 ）欧拉改进法 在 matlab 命令窗口输入 [t,x,y,z]=reuler(doty,doty2,doty3,0,5,0.1,0.01) t-w 曲线， t- δ曲线分别如下图所示。具体功角，角速度数据分别见文件 3.mat 和 4.mat （3 ）龙格库塔法 在 matlab 命令窗口输入 [t,x,y,z]=kunta(doty,doty2,doty3,0,5,0.1,0.01) t-w 曲线， t- δ曲线分别如下图所示。具体功角，角速度数据分别见文件 5.mat 和 6.mat 精选资料，欢迎下载 。 2 运行 Runge-Kutta, 将参数阻尼 D 设置为 0.05 ，不断更改参数切