harmonic.asv

Complex Wave Synthesizer

%Synthesis of complex wave
% Co-efficients and phases of sine wave
H=findobj(gcf,'Tag','b6');
b6=str2num(get(H,'String'));
H=findobj(gcf,'Tag','b5');
b5=str2num(get(H,'String')); 
H=findobj(gcf,'Tag','b4');
b4=str2num(get(H,'String'));
H=findobj(gcf,'Tag','b3');
b3=str2num(get(H,'String'));
H=findobj(gcf,'Tag','b2');
b2=str2num(get(H,'String'));
H=findobj(gcf,'Tag','b1');
b1=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pb6');
pb6=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pb5');
pb5=str2num(get(H,'String')); 
H=findobj(gcf,'Tag','pb4');
pb4=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pb3');
pb3=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pb2');
pb2=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pb1');
pb1=str2num(get(H,'String'));

%Co-efficients phases of cosine wave
H=findobj(gcf,'Tag','a6');
a6=str2num(get(H,'String'));
H=findobj(gcf,'Tag','a5');
a5=str2num(get(H,'String')); 
H=findobj(gcf,'Tag','a4');
a4=str2num(get(H,'String'));
H=findobj(gcf,'Tag','a3');
a3=str2num(get(H,'String'));
H=findobj(gcf,'Tag','a2');
a2=str2num(get(H,'String'));
H=findobj(gcf,'Tag','a1');
a1=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pa6');
pa6=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pa5');
pa5=str2num(get(H,'String')); 
H=findobj(gcf,'Tag','pa4');
pa4=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pa3');
pa3=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pa2');
pa2=str2num(get(H,'String'));
H=findobj(gcf,'Tag','pa1');
pa1=str2num(get(H,'String'));
% DC term
H=findobj(gcf,'Tag','a0');
a0=str2num(get(H,'String'));
%time 
H=findobj(gcf,'Tag','t0'); 
t0=str2num(get(H,'String'));
H=findobj(gcf,'Tag','t1'); 
t1=str2num(get(H,'String'));
H=findobj(gcf,'Tag','t2'); 
t2=str2num(get(H,'String'));
% frequency 
H=findobj(gcf,'Tag','f'); 
f=str2num(get(H,'String'));
 
t=t0:t1:t2;
va0=a0;
% Sine terms 
vb1  = b1*sin(2*pi*f*t+pb1/57.3);
vb2  = b2*sin(2*2*pi*f*t+pb2/57.3);
vb3  = b3*sin(3*2*pi*f*t+pb3/57.3);
vb4  = b4*sin(4*2*pi*f*t+pb4/57.3);
vb5  = b5*sin(5*2*pi*f*t+pb5/57.3);
vb6  = b6*sin(6*2*pi*f*t+pb6/57.3);
% Cosine terms
va1  = a1*cos(2*pi*f*t+pa1/57.3);
va2  = a2*cos(2*2*pi*f*t+pa2/57.3);
va3  = a3*cos(3*2*pi*f*t+pa3/57.3);
va4  = a4*cos(4*2*pi*f*t+pa4/57.3);
va5  = a5*cos(5*2*pi*f*t+pa5/57.3);
va6  = a6*cos(6*2*pi*f*t+pa6/57.3);

v=va0+va1+va2+va3+va4+va5+va6+vb1+vb2+vb3+vb4+vb5+vb6;
 
y0=a0+a1+a2+a3+a4+a5+a6+b1+b2+b3+b4+b5+b6;
y1=-y0;
plot(t,va0,'b',t,va1,'g',t,va2,'y',t,va3,'c',t,va4,t,va5,t,va6,t,vb1,'m',t,vb2,'v',t,vb3,t,vb4,t,vb5,t,vb6,t,v,'r-+')
grid on
xlabel('Time')
ylabel('Amplitude')