simulat.m

同步发电机空载三相两项短路情况下的各相波形以及励磁电流波形

x_d = 1.75581; x_ad = 1.56711; x_F = 1.73883; x_D = 1.63555; r_F = 0.00113; r_D = 0.00929; x_FD = 1.61805;
x_q = 1.75581; x_aq = 1.56711; x_H = 2.22911; x_Q = 1.64064; r_H = 0.01862; r_Q = 0.01788; r = 0.004;
DELT_t = 0.0005; t_max = 0.5 ;U = 1;sita_0 = 1.57;
T=0:DELT_t:1;
I=zeros(6,length(T));                      %%%I(:,j)表示第j点的值的存储
U_dq=zeros(2,length(T));
omg=1.0;
T_B=T*100*pi;
delt_t=DELT_t*100*pi;
U_dq0=[0;1];u_q0=U_dq0(2);
I(:,1)=[0;0;u_q0/x_ad;0;0;0];

C=[-(x_d+delt_t*r/2)    delt_t*omg*x_q/2               x_ad                 x_ad               -delt_t*omg*x_aq/2         -delt_t*omg*x_aq/2
    -delt_t*omg*x_d/2   -(x_q+delt_t*r/2)     delt_t*omg*x_ad/2       delt_t*omg*x_ad/2          x_aq                     x_aq
    -x_ad               0                       x_F+delt_t*r_F/2         x_FD                    0                         0
    -x_ad               0                           x_FD               x_D+delt_t*r_D/2          0                         0
    0                  -x_aq                        0                    0                      x_H+delt_t*r_H/2          x_aq
    0                  -x_aq                        0                    0                      x_aq                      x_Q+delt_t*r_Q/2];   %%%常矩阵C

j=length(T);
D=[delt_t/2 0 0 0 0 0
    0 delt_t/2 0 0 0 0]';
u_F=r_F*I(3,1);

for j=2:length(T)
    e1 = -(x_d-delt_t*r/2)*I(1,j-1) - delt_t*omg*x_q/2*I(2,j-1) + x_ad*I(3,j-1) + x_ad*I(4,j-1) + delt_t*omg*x_aq/2*I(5,j-1) + delt_t*omg*x_aq/2*I(6,j-1);
    e2 = delt_t*omg*x_d/2*I(1,j-1) - (x_q-delt_t*r/2)*I(2,j-1) - delt_t*omg*x_ad/2*I(3,j-1) - delt_t*omg*x_ad/2*I(4,j-1) + x_aq*I(5,j-1) + x_aq*I(6,j-1);
    e3 = -x_ad*I(1,j-1) + (x_F-delt_t*r_F/2)*I(3,j-1) + x_FD*I(4,j-1) + delt_t*u_F;
    e4 = -x_ad*I(1,j-1) + x_FD*I(3,j-1) + (x_D-delt_t*r_D/2)*I(4,j-1);
    e5 = -x_aq*I(2,j-1) + (x_H-delt_t*r_H/2)*I(5,j-1) +x_aq*I(6,j-1);
    e6 = -x_aq*I(2,j-1) + x_aq*I(5,j-1) + (x_Q+delt_t*r_Q/2)*I(6,j-1);
    E=[e1;e2;e3;e4;e5;e6];
    I(:,j)=C\E;
end
I_a=I(1,:).*cos(T_B)-I(2,:).*sin(T_B);
I_b=I(1,:).*cos(T_B-2*pi/3)-I(2,:).*sin(T_B-2*pi/3);
I_c=I(1,:).*cos(T_B+2*pi/3)-I(2,:).*sin(T_B+2*pi/3);
figure(1)
plot(T,I_a);
figure(2)
plot(T,I_b);
figure(3)
plot(T,I_c);
figure(4);
plot(T,I(3,:));