Respuestas (1)

Geoff Hayes
Geoff Hayes el 17 de Oct. de 2015

0 votos

Abduladeam - in order to get the square wave, you will have to sum over the odd k (which is missing from your above equation). Also, you need to divide the 4 by (k*pi). Try the following code
freq = 1;
func = @(k,t)(4/(k*pi))*sin(2*pi*k*t*freq);
t = linspace(0,1,500);
k = 1;
% create a vector for k==1 across all t
fkt = func(k,t);
% now sum over the odd numbers
upperBound = 100;
for k=3:2:upperBound
fkt = fkt + func(k,t);
end
plot(t,fkt);
Adjust the freq to increase the frequency of the wave, and increase the upperBound in order to improve upon the square wave's approximation.

2 comentarios

Abduladeam abdullah
Abduladeam abdullah el 17 de Oct. de 2015
thanks but this is the main equation f(t)=(4/m*pi)sin(2*pi*m*t) m=1,3,5,7,9
and I need to plot the sum of harmonics against time for a time period between 0 and 10
f(t)=(4/pi)sin(2*pi*t)+(4/3*pi)sin(2*pi*3*t)+(4/5*pi)sin(2*pi*5*t)+(4/7*Pi)sin(2*pi*7*t)
Geoff Hayes
Geoff Hayes el 17 de Oct. de 2015
I'm not sure I understand your question then. Can you not use the above to code for your work?

Iniciar sesión para comentar.

Categorías

Más información sobre Fourier Analysis and Filtering en Centro de ayuda y File Exchange.

Etiquetas

Aún no se han introducido etiquetas.

Preguntada:

el 17 de Oct. de 2015

Comentada:

el 17 de Oct. de 2015

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by