Plotting The Amplitude for Fourier Series Coefficients

%% Problem 1

display('Problem 1')

syms n t

Wo = pi;

T = 2;

n = 1:10;

a0=(1/T)*(int(t,t,0,1))

an= (2/T)*(int(t*cos(n*Wo*t),0,1))

bn=(2/T)*(int(t*sin(n*Wo*t),0,1))

Answer 1

[Ameer Hamza]

Something like this

syms n t
Wo = pi;
T = 2;
n = 1:10;

a0 = (1/T)*int(t,t,0,1);
an = (2/T)*int(t*cos(n*Wo*t),0,1);
bn = (2/T)*int(t*sin(n*Wo*t),0,1);

subplot(2,1,1)
title('an');
stem(0:10, [a0 an])
subplot(2,1,2)
title('bn');
stem(1:10, bn)

[Salik Mallick]

Thank you, but how do I plot the actual function?

[Ameer Hamza]

Following plot actual function

syms n t
Wo = pi;
T = 2;
n = 1:10;

a0 = (1/T)*int(t,t,0,1);
an = (2/T)*int(t*cos(n*Wo*t),0,1);
bn = (2/T)*int(t*sin(n*Wo*t),0,1);

t = linspace(0, 1);
x = zeros(size(t));

x(:) = a0;
for i = 1:numel(an)
    x = x + an(i)*cos(Wo*i*t) + bn(i)*sin(Wo*i*t);
end

plot(t, x)