%This application is used to illustrate the harmonic
clc
clear all
vd=1;
n=1;
T=0:.0001:pi;
v=cos(n.*T./2);
v1=2.*sqrt(2)./pi.*vd./2.*cos(n.*T./2);
T=180/pi*T;            %convertingradiantodegree
plot(T,v)
title(' Harmonic Voltage ')
xlabel('Angle')
ylabel('voltage perunit')
grid on
hold on
n=3;                    %harmonic3
T=0:.0001:pi;
vn=abs(2.*sqrt(2)./pi.*vd./2./n.*cos(n.*T./2));
T=180/pi*T;            %convertingradiantodegree
plot(T,vn./v1,'r')
hold on
n=5;                         %harmonic5
T=0:.0001:pi;
vn=abs(2.*sqrt(2)./pi.*vd./2./n.*cos(n.*T./2));
T=180/pi*T;                 %convertingradiantodegree
plot(T,vn./v1,'g')
hold on
n=7;                        %harmonic7
T=0:.0001:pi;
vn=abs(2.*sqrt(2)./pi.*vd./2./n.*cos(n.*T./2));
T=180/pi*T;              %convertingradiantodegree
plot(T,vn./v1,'k')
hold on
n=9;                      %harmonic9
T=0:.0001:pi;
vn=abs(2.*sqrt(2)./pi.*vd./2./n.*cos(n.*T./2));
T=180/pi*T;                %convertingradiantodegree
plot(T,vn./v1,'y')
hold on
n=11;                      %harmonic11
T=0:.0001:pi;
vn=abs(2.*sqrt(2)./pi.*vd./2./n.*cos(n.*T./2));
T=180/pi*T;                %convertingradiantodegree
plot(T,vn./v1,'m')