《MATLAB/Simulink电力系统建模与仿真》第2版p125
简单系统,短路故障(II段)发电机摇摆曲线δ-t(转子运动方程——两个一阶非线性常微分方程):
已知切除故障时间(III段),求δ-t曲线判断系统稳定性
教材常用分段计算法和常微分方程数值解法——改进欧拉法,但是由于MATLAB进行常微分方程求解算法中没有改进欧拉法,编程时采用Runge-Kutta龙格-库塔法。
%****发电机转子摇摆曲线的非线性一阶微分方程的Runge-Kutta龙格-库塔法的m函数***
求解微分方程程序
%求解发电机转子摇摆曲线的微分方程
建立simulink仿真模型: