clear,clc,clf;
Vm=30000*sqrt(2);
w=2*pi*50;
T=1/50;
dt=1e-4;
k=0;
for t=0:dt:50*T
k=k+1;
tt(k)=t;
Va(k)=Vm*sin(w*t);
Vb(k)=Vm*sin(w*t-(2*pi/3));
Vc(k)=Vm*sin(w*t+(2*pi/3));
end;
%grid on;
%plot(Va);
%hold on;
%plot(Vb);
%hold on;
%plot(Vc);
k=0;
for t=(2*pi/(3*w)):dt:50*T
k=k+1;
Vd(k)=(2/3)*(Va(k)*sin(w*t)+Vb(k)*sin(w*t-2*pi/3)+Vc(k)*sin(w*t+(2*pi/3)));
Vq(k)=(2/3)*(Va(k)*cos(w*t)+Vb(k)*cos(w*t-2*pi/3)+Vc(k)*cos(w*t+2*pi/3));
Vo(k)=(1/3)*(Va(k)+Vb(k)+Vc(k));
end
subplot(3,1,1);
grid on;
plot(Vd);
subplot(3,1,2);
hold on;
plot(Vq);
subplot(3,1,3);
hold on;
plot(Vo);
