# how can I plot this function?

3 views (last 30 days)
Niloufar on 26 Oct 2022
Answered: David Hill on 27 Oct 2022
My problem is how to plot this function
Niloufar on 27 Oct 2022
f(t) the bound for t is -5 to 5
Niloufar on 27 Oct 2022
close all;clear;clc; t = linspace(-5,5,1000); %ylim([-0.01 0.01]); N = 1000;
syms n t; an = (100*n*pi*(-1).^(n+1).*(1000-121*n.*n*pi*pi))/((1000-121*n.*n*pi*pi)^2 + (1320*n*pi)^2)*(50*n*pi); bn = (100*n*pi*(-1)^(n+1)*(1320*n*pi))/((1000-121*n.*n*pi*pi)^2 + (1320*n*pi)^2)*(50*n*pi); f = symsum(an.*cos(100*n*pi*t) + bn.*sin(100*n*pi*t),n,1,Inf); fplot(f);
I have done it with Symsum but it doesn't work.what is the problem

David Hill on 27 Oct 2022
t = linspace(-5,5,1000);
N = 1000;n=1:N;
an = (100*n*pi.*(-1).^(n+1).*(1000-121*n.^2*pi^2))/(((1000-121*n.^2*pi^2).^2 + (1320*n*pi).^2).*(50*n*pi));
bn = (100*n*pi.*(-1).^(n+1).*(1320*n*pi))./(((1000-121*n.^2*pi^2).^2 + (1320*n*pi).^2).*(50*n*pi));
f = @(t)sum(an.*cos(100*n*pi*t) + bn.*sin(100*n*pi*t));
for k=1:length(t)
F(k)=f(t(k));
end
plot(t,F)