代码:
xd=2.26;
xq=2.26;
xd1=0.269; %xd’
xq1=0.269;%xq’
xd2=0.167; %xd’’
xq2=0.167; %xq’’
Eq0=1;
t=0;
Ta=2;
Td1=0.85; %Td’
Td2=0.1; %Td’’
w0=100*pi;
theta=0;
figure(1);
for theta=0:0.5:pi
t=[0:0.01:10];
ia=power(2,1/2)*Eq0*[(1/xd2-1/xd1)*exp(-t/Td2)+(1/xd1-1/xd)*exp(-t/Td1)+1/xd].*cos(theta+w0*t)-power(2,1/2)*Eq0*cos(theta)*exp(-t/Ta)/xd2;
plot(t,ia),xlabel('t'),ylabel('ia'),title('ia Graph');
hold on
end
figure(2);
for theta=0:0.5:pi
t=[0:0.01:10];
ib=power(2,1/2)*Eq0*[(1/xd2-1/xd1)*exp(-t/Td2)+(1/xd1-1/xd)*exp(-t/Td1)+1/xd].*cos(theta+w0*t-2*pi/3)-power(2,1/2)*Eq0*cos(theta-2*pi/3)*exp(-t/Ta)/xd2;
plot(t,ib),xlabel('t'),ylabel('ib'),title('ib Graph');
hold on
end
figure(3)
for theta=0:0.5:pi
t=[0:0.01:10];
ic=power(2,1/2)*Eq0*[(1/xd2-1/xd1)*exp(-t/Td2)+(1/xd1-1/xd)*exp(-t/Td1)+1/xd].*cos(theta+w0*t+2*pi/3)-power(2,1/2)*Eq0*cos(theta+2*pi/3)*exp(-t/Ta)/xd2;
plot(t,ic),xlabel('t'),ylabel('ic'),title('ic Graph');
hold on
end
图例是生成图片后在figure里加的