how can I plot this function?

1 次查看(过去 30 天)
Niloufar
Niloufar 2022-10-26
回答: David Hill 2022-10-27
My problem is how to plot this function
  5 个评论
Niloufar
Niloufar 2022-10-27
f(t) the bound for t is -5 to 5
Niloufar
Niloufar 2022-10-27
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

请先登录,再进行评论。

回答(1 个)

David Hill
David Hill 2022-10-27
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)

类别

Help CenterFile Exchange 中查找有关 Calculus 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by