function fangzhen 
tspan = [0  1];
y0 = [0;0;17.8876/0.13;0;0];
[t,y]=ode23tb('tongbufadianji',tspan,y0);
subplot(5,1,1)
plot(t,y(:,1))
xlabel('Time[s]')
ylabel('Id[A]')

hold on 
subplot(5,1,2)
plot(t,y(:,2))
xlabel('Time[s]')
ylabel('Iq[A]')

hold on
subplot(5,1,3)
plot(t,y(:,3))
xlabel('Time[s]')
ylabel('Ifd[A]')

hold on
subplot(5,1,4)
plot(t,y(:,4))
xlabel('Time[s]')
ylabel('Ikd[A]')

hold on
subplot(5,1,5)
plot(t,y(:,5))
hold on
xlabel('Time[s]')
ylabel('Ikq[A]')