ASN11 - Leakage demonstration in Octave: Difference between revisions
Jodi.Hodge (talk | contribs) No edit summary |
Jodi.Hodge (talk | contribs) No edit summary |
||
(4 intermediate revisions by the same user not shown) | |||
Line 1: | Line 1: | ||
[[Jodi Hodge|Back to my home page]] |
[[Jodi Hodge|Back to my home page]] |
||
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. |
|||
Leakage demonstration in Octave |
|||
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 |
Latest 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