Scribd
Upload
423 views3 pages

Signal Processing in MATLAB

This document discusses using MATLAB for signal processing and Fourier transforms. It introduces the fft and ifft commands to compute the Fourier transform and inverse Fourier transform of a signal. The document provides an example of applying these commands to a sample signal to obtain the waveform and amplitude spectrum. It then discusses how fft and ifft can be used to filter noise from signals by removing components in the Fourier domain below a certain amplitude threshold. Assignments are provided asking the reader to experiment further with these Fourier transform concepts.

Uploaded by

Radu Vanca
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
Download as pdf or txt
423 views3 pages

Signal Processing in MATLAB

This document discusses using MATLAB for signal processing and Fourier transforms. It introduces the fft and ifft commands to compute the Fourier transform and inverse Fourier transform of a signal. The document provides an example of applying these commands to a sample signal to obtain the waveform and amplitude spectrum. It then discusses how fft and ifft can be used to filter noise from signals by removing components in the Fourier domain below a certain amplitude threshold. Assignments are provided asking the reader to experiment further with these Fourier transform concepts.

Uploaded by

Radu Vanca
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
Download as pdf or txt
Download as pdf or txt
You are on page 1/ 3

MCS 320

Introduction to Symbolic Computation

Spring 2007

MATLAB Lecture 7. Signal Processing in MATLAB


We have seen how to t data with polyt and how to design shapes with spline. In this lecture we cover another way to deal with approximate data, which is well suited for treating signals [1].

7.1 Fourier Transforms


Every signal can be written as a sum of sinusoids with dierent amplitudes and frequencies. The MATLAB command to compute the Fourier Transform and its inverse are respectively t and it, for example: >> >> >> >> >> x = rand(1,10); y = fft(x); iy = ifft(y); x2 = real(iy); norm(x-x2); % % % % % suppose 10 samples of a random signal Fourier transform of the signal inverse Fourier transform chop off tiny imaginary parts compare original with inverse of transformed

The t is the abbreviation of Fast Fourier Transform. This algorithm implements the discrete Fourier transform to transform data from time into the frequency domain. The study of this algorithm is normally covered in a good linear algebra course. First we give an example of the meaning of the Fourier transform before showing how Fourier transforms can be used to lter noise from signals.

7.2 Waveform and Amplitude Spectrum


Suppose we sample a signal during 4 seconds, at a sampling rate of 0.01: >> >> >> >> dt = 1/100; et = 4; t = 0:dt:et; y = 3*sin(4*2*pi*t) + 5*sin(2*2*pi*t); % % % % sampling rate end of the interval sampling range sample the signal

A natural plot is that of amplitude versus time: >> >> >> >> >> subplot(2,1,1); plot(t,y); grid on axis([0 et -8 8]); xlabel(Time (s)); ylabel(Amplitude); % % % % % first of two plots plot with grid adjust scaling time expressed in seconds amplitude as function of time

With the Fourier Transform we can visualize what characterizes this signal the most. From the Fourier transform we compute the amplitude spectrum: >> Y = fft(y); >> n = size(y,2)/2; >> amp_spec = abs(Y)/n; % compute Fourier transform % 2nd half are complex conjugates % absolute value and normalize

To visualize the amplitude spectrum, we execute the following commands >> >> >> >> >> subplot(2,1,2); freq = (0:79)/(2*n*dt); plot(freq,amp_spec(1:80)); grid on xlabel(Frequency (Hz)); ylabel(Amplitude); % % % % % second of two plots abscissa viewing window plot amplitude spectrum 1 Herz = number of cycles/second amplitude as function of frequency

On the amplitude spectrum we see two peaks: at 2 and 4. The location of the peaks occurs at the two frequencies in the signal. The heights of the peaks (5 and 3) are the amplitudes of the sines in the signal. Jan Verschelde, 23 April 2007 UIC, Dept of Math, Stat & CS MATLAB Lecture 7, page 1

MCS 320

Introduction to Symbolic Computation

Spring 2007

7.3 Filtering Noise from Signals


We will see now how to use t and it to lter out the noise from signals. First we add random noise to our signal and compute the amplitude spectrum. >> >> >> >> >> noise = randn(1,size(y,2)); ey = y + noise; eY = fft(ey); n = size(ey,2)/2; amp_spec = abs(eY)/n; % % % % % random noise samples with noise Fourier transform of noisy signal use size for scaling compute amplitude spectrum

To interpret these calculations we make a plot of the waveform and amplitude spectrum: >> >> >> >> >> >> >> >> >> >> >> figure % plots in new window subplot(2,1,1); % first of two plots plot(t,ey); grid on % plot noisy signal with grid axis([0 et -8 8]); % scale axes for viewing xlabel(Time (s)); % time expressed in seconds ylabel(Amplitude); % amplitude as function of time subplot(2,1,2); % second of two plots freq = (0:79)/(2*n*dt); % abscissa viewing window plot(freq,amp_spec(1:80)); grid on % plot amplitude spectrum xlabel(Frequency (Hz)); % 1 Herz = number of cycles per second ylabel(Amplitude); % amplitude as function of frequency

On the rst plot we recognize the shape of the signal. In the plot of the amplitude spectrum, the peaks and their heights are the same as on the plot of the amplitude spectrum of the original signal. The wobbles we see around the peaks show that the amplitude of the noise is less than that of the original signal. We can visualize the output of the Fourier transforms: >> >> >> >> figure plot(Y/n,r+) hold on plot(eY/n,bx) % % % % new window for plot Fourier transform of original put more on same plot Fourier transform of noisy signal

Via the inverse Fourier transform, we lter out the noise. The command x rounds the elements of its argument to the nearest integers towards zero. For this example, we use x to set all elements in eY less than 100 to zero: >> fY = fix(eY/100)*100; >> ifY = ifft(fY); >> cy = real(ifY); % set numbers < 100 to zero % inverse Fourier transform of fixed data % remove imaginary parts

The vector cy contains the corrected samples. So, nally we plot this corrected signal: >> >> >> >> >> figure plot(t,cy); grid on axis([0 et -8 8]); xlabel(Time (s)); ylabel(Amplitude); % % % % % new window for plot plot corrected signal adjust scale for viewing time expressed in seconds amplitude as function of time

Here we ltered out noise of low amplitude. Note we can also remove noise of high frequency.

Jan Verschelde, 23 April 2007

UIC, Dept of Math, Stat & CS

MATLAB Lecture 7, page 2

MCS 320

Introduction to Symbolic Computation

Spring 2007

7.4 Assignments
1. Consider the following sequence of instructions: >> >> >> >> >> >> t = 0:0.1:10; y1 = sin(2*pi*t); y2 = sin(20*pi*t); plot(t,y1); hold on; plot(t,y2);

Why is the output of the second plot like this? Find a better range for t to plot sin(20*pi*t) right. Can you nd a good lower bound for the sampling interval in terms of the frequency? 2. Give the MATLAB commands to plot the amplitude spectrum for the signal
20

f (t) =
k=10

(20 k) sin(2kt).

In addition, plot the waveform spectrum of this signal. 3. Make a function to plot waveform and amplitude spectrum of a signal. The function has prototype: function specplot ( t, dt, et, y ) % % Opens a new figure window with two plots: % the waveform and amplitude spectrum of a signal. % % On entry : % t sampling range of the signal; % dt sampling rate; % et end of the range; % y samples of the signal over the range t. % So specplot computes the amplitude spectrum of the signal. For the abscissa viewing window you may take half of the range of t. Test your specplot with the signal of the previous assignment. 4. With t we can decompose a signal in low and high frequencies. Take the example signal from page 1. As noise we now add a sine of amplitude 4 and with frequency 50. Plot the waveform and amplitude spectrum of the new signal. Use t and it to remove this high frequency noise.

References
[1] S.D. Stearns and R.A. David. Signal Processing Algorithms in MATLAB. Prentice Hall, 1996.

Jan Verschelde, 23 April 2007

UIC, Dept of Math, Stat & CS

MATLAB Lecture 7, page 3

You might also like