Fourier Series of a Squarewave

Question

Salik Mallick 2020-9-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/601738-fourier-series-of-a-squarewave

编辑： Ameer Hamza 2020-9-29

采纳的回答： Ameer Hamza

Why is this script not running for a Fourier Series of a Squarewave?

T=2

ao=1/2

an=0

bn= (-1/n*pi)*[ (-1)^n - 1 ])

%% Problem 1_Practice_Problem_17.1

display('Problem 1_Practice_Problem_17.1')

N = 10000;

Wo = pi;

t = -1:.001:11;

f = zeros(1,length(t));

for n = 1:2*n:N

an = [1-[(-1)^n]/n*pi;

f=f + an*cos(Wo*n*t);

end

% Plotting and plot formatting is done here

plot(t,f)

axis([-1 11 -1.1 1.1])

title(['Fourier Series Representation w/ ' num2str(N) ' Terms'],'FontSize',16)

xlabel('Time')

ylabel('Function Value')

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Ameer Hamza 2020-9-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/601738-fourier-series-of-a-squarewave#answer_502243

编辑：Ameer Hamza 2020-9-29

在 MATLAB Online 中打开

Check this code

ao=1/2;
an=0;
%% Problem 1_Practice_Problem_17.1
disp('Problem 1_Practice_Problem_17.1')
 
N = 10000;
Wo = pi;
t = -1:.001:11;
f = zeros(1,length(t)); 
for n = 1:2:N 
    bn = (1-(-1)^n)./(n*pi);
    f = f + an*cos(Wo*n*t) + bn*sin(Wo*n*t); 
end
% Plotting and plot formatting is done here
plot(t,f)
axis([-1 11 -1.1 1.1])
title(['Fourier Series Representation w/ ' num2str(N) ' Terms'],'FontSize',16)
xlabel('Time')
ylabel('Function Value')