An octave/MATLAB script to show the Fourier series of a string of impulse functions: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 10: Line 10:
for k=-M:M
for k=-M:M
v = 1;
v = 1;
f = f+v*exp(-j*2*pi*k*t/T);
f = f +v*exp(j*2*pi*k*t/T);
end
end
plot(t,f)
plot(t,f)
title(strcat('\Sigma_{k=-\infty}^\infty \delta(t-kT) \approx \Sigma_{k=-M}^M e^{-j2\pi nt/T}'
,'(M = ',num2str(2*M),' Terms, and T =',num2str(T),')'))
%title(strcat('Fourier Series of a Sum ofImpulse Functions with ',num2str(2*M),' Terms'))
xlabel('Time (S)')
</nowiki>
</nowiki>
[[File:Sampling_Function.jpg|800px]]

Revision as of 16:42, 2 November 2016

This checks this identity.

% This is a script to see if the Fourier series with unity coefficients is
% really a series of impulse functions.
M=100; % Number of terms 
T=1e-3;  % the period of the sampling function is 1 mS.
t=0:T/1000:10*T;

f=zeros(size(t));
for k=-M:M
  v = 1;
  f = f +v*exp(j*2*pi*k*t/T);
end
plot(t,f)
title(strcat('\Sigma_{k=-\infty}^\infty \delta(t-kT) \approx \Sigma_{k=-M}^M e^{-j2\pi nt/T}'
,'(M =  ',num2str(2*M),' Terms, and T =',num2str(T),')'))
%title(strcat('Fourier Series of a Sum ofImpulse Functions with ',num2str(2*M),' Terms'))
xlabel('Time (S)')

File:Sampling Function.jpg