Introduction

In this lab, we will look at the effect of filtering signals using a frequency domain implementation of an LTI system, i.e., multiplying the Fourier transform of the input signal with the frequency response of the system. In particular, we will filter sound signals, and investigate both low-pass and high-pass filters. Recall that a low-pass filter filters out high frequencies, allowing only the low frequencies to pass through. A high-pass filter does the opposite.

MATLAB Commands and Resources

All of the sounds for this lab can be downloaded from the Sound Resources page.

Transforming Signals to the Frequency Domain and Back

When working in MATLAB, the continuous-time Fourier transform cannot be done by the computer exactly, but a digital approximation is done instead. The approximation uses the discrete Fourier transform (more on that in EE 341). There are a couple important differences between the discrete Fourier transforms on the computer and the continuous Fourier transforms you are working with in class: finite frequency range and discrete frequency samples. The frequency range is related to the sampling frequency of the signal. In the example below, where we find the Fourier transform of the "fall" signal, the sampling frequency is Fs=8000, so the frequency range is [-4000,4000] Hz (or 2*pi times that for w in radians). The frequency resolution depends on the length of the signal (which is also the length of the frequency representation).

The MATLAB command for finding the Fourier transform of a signal is fft (for Fast Fourier Transform (FFT)). In this class, we only need the default version.

>> load fall     %load in the signal
>> x = fall;
>> X = fft(x);
           

>> N = length(x);
>> pfreq = [0:N/2]*Fs/N;
>> Xpos=X(1:N/2+1);
>> plot(pfreq,abs(Xpos));
>> figure;
>> freq = [-(N/2-1):N/2]*Fs/N;
>> plot(freq,abs(fftshift(X)));
            

Try looking at the frequency content of a few other signals. Note that the fall signal happens to have an even length, so N/2 is an integer. If the length is odd, you may have indexing problems, so it is easiest to just omit the last sample, as in x=x(1:length(x)-1);.

After you make modifications of a signal in the frequency domain, you typically want to get back to the time domain. The MATLAB command ifft will accomplish this task.

>> xnew = real(ifft(X));
            

Low-pass Filtering

An ideal low-pass filter eliminates high frequency components entirely, as in: HL ideal (ω)={ 1 | ω |≤B 0 | ω |>B } HL ideal (ω)={ 1 | ω |≤B 0 | ω |>B } MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadIeadaqhaaWcbaGaamitaaqaaiaadMgacaWGKbGaamyzaiaadggacaWGSbaaaOGaaiikaiabeM8a3jaacMcacqGH9aqpdaGadaqaauaabeqaciaaaeaacaaIXaaabaWaaqWaaeaacqaHjpWDaiaawEa7caGLiWoacqGHKjYOcaWGcbaabaGaaGimaaqaamaaemaabaGaeqyYdChacaGLhWUaayjcSdGaeyOpa4JaamOqaaaaaiaawUhacaGL9baaaaa@525A@ A real low-pass filter typically has low but non-zero values for | H L (ω) | | H L (ω) | MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaaemaabaGaamisamaaBaaaleaacaWGmbaabeaakiaacIcacqaHjpWDcaGGPaaacaGLhWUaayjcSdaaaa@3DFE@ at high frequencies, and a gradual (rather than an immediate) drop in magnitude as frequency increases. The simplest (and least effective) low-pass filter is given by (e.g. using an RC circuit): H L (ω)= α α+jω , α=cutoff frequency. H L (ω)= α α+jω , α=cutoff frequency. MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadIeadaWgaaWcbaGaamitaaqabaGccaGGOaGaeqyYdCNaaiykaiabg2da9maalaaabaGaeqySdegabaGaeqySdeMaey4kaSIaamOAaiabeM8a3baacaGGSaGaaeiiaiabeg7aHjabg2da9iaabogacaqG1bGaaeiDaiaab+gacaqGMbGaaeOzaiaabccacaqGMbGaaeOCaiaabwgacaqGXbGaaeyDaiaabwgacaqGUbGaae4yaiaabMhacaqGUaaaaa@5620@

This low-pass filter can be implemented in MATLAB using what we know about the Fourier transform. Remember that multiplication in the Frequency domain equals convolution in the time domain. If our signal and filter are both in the frequency domain, we can simply multiply them to produce the result of the system. y(t)=x(t)∗h(t) Y(ω)=X(ω)H(ω) y(t)=x(t)∗h(t) Y(ω)=X(ω)H(ω) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOabaeqabaGaamyEaiaacIcacaWG0bGaaiykaiabg2da9iaadIhacaGGOaGaamiDaiaacMcacqGHxiIkcaWGObGaaiikaiaadshacaGGPaaabaGaamywaiaacIcacqaHjpWDcaGGPaGaeyypa0JaamiwaiaacIcacqaHjpWDcaGGPaGaamisaiaacIcacqaHjpWDcaGGPaaaaaa@4EBC@ Below is an example of using MATLAB to perform low-pass filtering on the input signal x with the FFT and the filter definition above.

The cutoff of the low-pass filter is defined by the constant a. The low-pass filter equation above defines the filter H in the frequency domain. Because the definition assumes the filter is centered around w = 0, the vector w is defined as such.

>> load fall     %load in the signal
>> x = fall;
>> X = fft(x);   % get the Fourier transform (uncentered)

>> N = length(X);
>> a = 100*2*pi;
>> w = (-N/2+1:(N/2));  % centered frequency vector
>> H = a ./ (a + i*w);       % centered version of H
>> plot(w*Fs/N,abs(H))
           

>> Hshift = fftshift(H);     % uncentered version of H
>> Y = X .* Hshift';         % filter the signal
           

>> y = real(ifft(Y));
>> sound(x, Fs)        % original sound
>> sound(y, Fs)        % low-pass-filtered sound
            

>> y = y * (max(abs(x))/max(abs(y)))
            

High-pass Filtering

An ideal high-pass filter eliminates low frequency components entirely, as in: H H ideal (ω)={ 0 | ω |<B 1 | ω |≥B } H H ideal (ω)={ 0 | ω |<B 1 | ω |≥B } MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadIeadaqhaaWcbaGaamisaaqaaiaadMgacaWGKbGaamyzaiaadggacaWGSbaaaOGaaiikaiabeM8a3jaacMcacqGH9aqpdaGadaqaauaabeqaciaaaeaacaaIWaaabaWaaqWaaeaacqaHjpWDaiaawEa7caGLiWoacqGH8aapcaWGcbaabaGaaGymaaqaamaaemaabaGaeqyYdChacaGLhWUaayjcSdGaeyyzImRaamOqaaaaaiaawUhacaGL9baaaaa@5263@ A real high-pass filter typically has low but non-zero values for | H L (ω) | | H L (ω) | MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaamaaemaabaGaamisamaaBaaaleaacaWGmbaabeaakiaacIcacqaHjpWDcaGGPaaacaGLhWUaayjcSdaaaa@3DFE@ at low frequencies, and a gradual (rather than an immediate) rise in magnitude as frequency increases. The simplest (and least effective) high-pass filter is given by (e.g. using an RC circuit): H H (ω)=1− H L (ω)=1− α α+jω , α=cutoff frequency. H H (ω)=1− H L (ω)=1− α α+jω , α=cutoff frequency. MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadIeadaWgaaWcbaGaamisaaqabaGccaGGOaGaeqyYdCNaaiykaiabg2da9iaaigdacqGHsislcaWGibWaaSbaaSqaaiaadYeaaeqaaOGaaiikaiabeM8a3jaacMcacqGH9aqpcaaIXaGaeyOeI0YaaSaaaeaacqaHXoqyaeaacqaHXoqycqGHRaWkcaWGQbGaeqyYdChaaiaacYcacaqGGaGaeqySdeMaeyypa0Jaae4yaiaabwhacaqG0bGaae4BaiaabAgacaqGMbGaaeiiaiaabAgacaqGYbGaaeyzaiaabghacaqG1bGaaeyzaiaab6gacaqGJbGaaeyEaiaab6caaaa@5F6C@ This filter can be implemented in the same way as the low pass filter above.

Sound Separation