example3_4_2.m

来自「经典《信号与系统》教程的matlab例程,对深入理解信号与系统相关概念有很大帮助」· M 代码 · 共 10 行

% Gibbs
t=-3:6/1000:3;
N=input('input the Number of harmonics N=');
c0=0.5;
w0=pi;
xN=c0*ones(1,length(t));
for k=1:2:N
    xN=xN-j*exp(j*k*pi*t)/(k*pi)+j*exp(-j*k*pi*t)/(k*pi);
end
plot(t,xN)