i want to plot a square wave

3 次查看（过去 30 天）
Abduladeam abdullah
Abduladeam abdullah 2015-10-17
评论： Star Strider 2015-10-17
the equation is f(t)=(4/k*pi)sin(2*pi*k*t) and k is odd number

请先登录，再回答此问题。

回答（1 个）

Geoff Hayes
Geoff Hayes 2015-10-17
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 个评论
Abduladeam abdullah
Abduladeam abdullah 2015-10-17
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 2015-10-17
I'm not sure I understand your question then. Can you not use the above to code for your work?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

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

Start Hunting!

Translated by