ASN11 - Leakage demonstration in Octave: Difference between revisions
Jodi.Hodge (talk | contribs) No edit summary |
Jodi.Hodge (talk | contribs) No edit summary |
||
Line 9: | Line 9: | ||
fs=1500 |
fs=1500 |
||
t=[1:20]/fs; %is multiple of the period 1/fs |
|||
t=[1:20]/fs; |
|||
y=sin(2*pi*t); |
y=sin(2*pi*t); |
||
x= dft(y,t); %function x is continous |
|||
x= dft(y); |
|||
Example octave script with leakage |
Example octave script with leakage |
||
Line 19: | Line 19: | ||
fs=1500 |
fs=1500 |
||
t=[1:20]/500; %does not include whole periods |
|||
t=[1:20]/fs; |
|||
y=sin(2*pi*t); |
y=sin(2*pi*t); |
||
x= dft(y,t); %function x is not continuous and demonstrates leakage |
|||
x= dft(y); |
Revision as of 15:21, 21 December 2009
Leakage occurs when the Discrete Fourier Transform (DFT) is applied to a signal whose range is not a multiple of the total period. The DFT function assumes that the signal is periodic it when it repeats the function should be periodic. If the sample range is not a multiple of the period them when the DFT repeat the function will have disctonuity at the restarting point.
Example octave script no leakage
fs=1500
t=[1:20]/fs; %is multiple of the period 1/fs
y=sin(2*pi*t);
x= dft(y,t); %function x is continous
Example octave script with leakage
fs=1500
t=[1:20]/500; %does not include whole periods
y=sin(2*pi*t);
x= dft(y,t); %function x is not continuous and demonstrates leakage