Scribd
Upload
127 views

Sampling Theorem

Procesamiento de señales biológicas, adquisición de señales con Matlab, Base de datos de Physionet.

Uploaded by

Israx Theron
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
127 views

Sampling Theorem

Procesamiento de señales biológicas, adquisición de señales con Matlab, Base de datos de Physionet.

Uploaded by

Israx Theron
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 33

2-9-2014

PROCESSING OF BIOSIGNALS
3BM3

PRACTICE 1. - SIGNAL ACQUISITION AND SAMPLING


THEOREM

DR. TOVAR CORONA BLANCA

ISRAEL MUOZ GARCA


UPIITA-IPN

INTRODUCTION
The sampling frequency defines the number of samples per second taken from a continuous signal
to make a discrete signal. The common notation for sampling frequency is fs.
The Nyquist-Shannon sampling theorem states that perfect reconstruction of a signal is possible
when the sampling frequency is greater than twice the maximum frequency of the signal being
sampled. If lower sampling rates are used, the original signals information may not be completely
recoverable from the sampled signal.

GENERAL PURPOSE
Apply the sampling frequency in order to convert, in this case, bio-signals from continuous signal to digital or discrete signal and
observe what would happen if we choose a sampling frequency lower than the maximum frequency.

PROCEDURE
PART I
I took bio-signals from the data base of PhysioNet. Each bio-signal has its frequency provided by PhysioNet and this is the code in
Matlab to obtain the sampling frequency and show the signal.
% in this case is a signal of an ECG at a frequency of 250Hz
fo =250;
fs = 2*fo; %sampling frequency
t1 = 0:1/fs: 10-(1/fs);
x = load('rec_1m.mat'); %File taken from physionet
figure,plot(t1,val),xlabel('Time(sec)'),ylabel('Amplitude(mV)'),title('ECG')
% EMG at a frequency of 2000Hz
fo =2000;
fs = 2*fo; %sampling frequency
t1 = 0:1/fs: 12.715-(1/fs);
x = load('emg_healthym.mat');
figure,plot(t1,val,'r'),xlabel('Time(sec)'),ylabel('Amplitude(mV)'),title('EMG')

% apnea ECG at a frequency of 50Hz


fo =50;
fs = 2*fo; %sampling frequency
t1 = 0:1/fs: 60-(1/fs);
x = load('a02m.mat');
figure,plot(t1,val,'k'),xlabel('Time(sec)'),ylabel('Amplitude(mV)'),title('Apnea ECG')

Part II
In this part of the practice we are going to show a sine wave at a frequency of 10 KHz and 1 KHz and we are going to compare the
signals increasing the frequency every 50 Hz.
Code in Matlab:
clc,clear all,close all;
fs=10000;
fs1=1000;
t=0:1/fs:0.05;
t1=0:1/fs1:0.05;
for i=50:50:1250 %increase the frequency every 50 Hz
x=sin(2*pi*i*t);
x2=sin(2*pi*i*t1);
figure,subplot(2,1,1),plot(t,x),xlabel('time(s)'),ylabel('Amplitude(V)'),title('10KHz'),
subplot(2,1,2),plot(t1,x2,'g'),xlabel('time(s)'),ylabel('Amplitude(V)'),title('1KHz')
end

RESULTS
PART I

PART II
Frequency at 50 Hz

100 Hz

150Hz

200Hz

250Hz

300Hz

350Hz

400Hz

450Hz

500Hz

550Hz

600Hz

650Hz

700Hz

750Hz

800Hz

850Hz

900Hz

950Hz

1000Hz

1050Hz

1100Hz

1150Hz

1200Hz

1250Hz

CONCLUSION
Increasing the sampling frequency increases the number of data points acquired in a given time period. Often, a fast sampling frequency provides
a better representation of the original signal than a slower sampling frequency. In part II we observe that an aliased signal provides a poor
representation of the analog signal because the sampling frequency is lower than the original signal and sometimes decompose it.

BIBLIOGRAPHY
http://www.princeton.edu/~achaney/tmve/wiki100k/docs/Sampling_rate.html
http://zone.ni.com/reference/en-XX/help/370051V-01/cvi/libref/analysisconcepts/aliasing/
http://www.physionet.org/physiobank/

You might also like