基于MATLAB的单机无穷大系统功角稳定仿真

clear all;clc; Pm=[1.351 0 1.1024];%before fault,during fault,after fault Pmax % Siga=41.77*pi/180;%rotor angle % delt=0;%speed of angle % Num=1;%time counter % h=0.005;%step of time % t=0;%time %% 仿真功角随时间的变化过程 Tc=[0.07 0.086 0.087 0.1]; % Tcm=Tc(1); % [result_1]=Rotor_angle(Tcm); % plot(result_1(:,3),result_1(:,1),'g'); % hold on; % Tcm=Tc(2); % [result_1]=Rotor_angle(Tcm); % plot(result_1(:,3),result_1(:,1),'b'); % hold on; % Tcm=Tc(3); % [result_1]=Rotor_angle(Tcm); % plot(result_1(:,3),result_1(:,1),'r'); % hold on; Tcm=Tc(4); [result_1]=Rotor_angle(Tcm); % plot(result_1(20:800,3),result_1(20:800,1),'k');%故障恢复后的功角随时间的变化情况曲线 % hold on; %% 绘制故障恢复后，测量阻抗由于功角摇摆随时间变化的过程 x1=result_1(20:142,1); % Ea=1.1626*exp(x1*1i/180*pi); % Eb=0.90081; y1=1.1626*exp(x1*1i)/(1.1626*exp(x1*1i)-0.90081)*0.95i-0.45*1i; a=real(y1);b=imag(y1); plot(a,b); hold on