So as part of an assignment in a course I'm doing,the task was to explore aliasing.
Just to recap,
"Aliasing is when a signal is not sampled correctly(undersampled with respect to its Nyquist Rate),and this causes the original signal to appear like a lower frequency signal."
What this means is,the frequencies in the signal happen to assume the false identity of lower frequency signals,usually due to very low sampling rates.
This is what most textbooks will tell you.But I never quite felt 'okay' with this definition,till I did this MATLAB exercise.
It might appear trivial,but it has really helped me.Let me walk you through it.
So,to keep things simple,lets consider a 1800Hz tone,which we will sample at 44.1kHz(CD Audio) initially.
And to get a real idea of things,lets use wavplay command to listen to this tone.Nice and shrill.
In terms of the Frequency Spectrum,what we have is this.
|Frequency Spectrum for an 1800Hz tone sampled at 44.1kHz|
Now the ideal sampling rate according to Harry,is twice the maximum frequency at which the signal has non-zero energy.So that would be 3.6kHz in our case.Which means,anything below this,and we will start seeing the effects of aliasing,and as we go higher away from this(6kHz,8kHz and so on) this effect will die out.
Lets sample at 3kHz,which is a mild case of undersampling.Set Fs=3000.Now when you listen to it,you can observe that the sound is 'duller',or its lowered in frequency,even though the Tone generating function hasnt changed.(ToneFreq still 1800).What happened then? Aliasing. Lets look at the frequency Spectrum now,and this is what you will see:
|Frequency Spectrum for an 1800Hz tone 'mildly' under-sampled at 3kHz|
The picture in time domain is like this:
|1800Hz tone suffering from Aliasing due to mild under-sampling.|
This effect is further prominent at even lower sampling frequencies.Lets take 1.5kHz as our sampling frequency,just to make Harry cringe. :P Now we'd see the 1800Hz tone appear like a 1800-1500=300Hz tone.And thats exactly what our plots reflect:
|Frequency Spectrum for an 1800Hz tone badly under-sampled at 1.5kHz|
|1800Hz tone suffering from Aliasing due to aggressive under-sampling.|
Pretty cool,eh? Hope this helped you understand aliasing as much as it helped me. :)
I'm not explicitly uploading the MATLAB code,but drop a line in the comments if youre interested.
I was inspired by this other really cool article from the Mathworks blog - http://blogs.mathworks.com/steve/2010/03/03/aliasing-and-a-sampled-cosine-signal/