Práctica 2
Martha Galarraga
31/10/2020
Señales y Dinámica de Biosistemas
Práctica 2-Analisis Temporal de Señales de Flujo Respiratorio
Introducción
La hipercapnia es una enfermedad producida por el exceso de CO2 (dióxido de
carbono) en el torrente sanguíneo. Esta produce insuficiencia respiratoria por
falta de ventilación pulmonar, provocando que se inhale en exceso CO2 que
luego no se logra evacuar.
Para la presente practica se busca realizar un código/función en la herramienta
Matlab cuyo propósito es detectar el inicio de los ciclos respiratorios con el
objetivo de obtener el período respiratorio (ciclos/segundos), para
posteriormente calcular la frecuencia respiratoria, la cual nos permitirá evaluar
el efecto del estímulo ritmo/frecuencia respiratoria con datos experimentales
de un sujeto.
Objetivo
Medir frecuencia instantánea de cada ciclo
Desarrollo
A) Desarrollo de la función en Matlab
1. Establecer los parámetros de entrada y salida de la función a crear, en
nuestro caso usaremos la función frec_ins.
function [fr, T]= frec_ins(flujo, fs)
%Función que calcula frecuencia instantanea respiratoria
%
%[fr, T]=frec_inst(flujo,fs)
%
%Argumentos de entrada: Flujo
% fs=frecuencia de muestreo
%Argumentos de salida: fr=frecuencia instantanea
% T=ciclos
2. En el caso que el numero de argumentos de entrada es menor que 2, le
asignaremos una frecuencia de muestreo fs=100.
if nargin < 2
Martha Galarraga Página 1 | 5
Práctica 2
fs=100;
end
3. Realizamos el cálculo de del ciclo de cada flujo, en la imagen podemos
observar que cada ciclo de flujo estará definido por cada cruce de menos a
más.
Ciclo de flujo
Figura 1: Representación de puntos de un ciclo respiratorio
El cálculo estaría dado por:
F01 = flujo>0; %flujo tipo onda cuadrado: positiva=1,negativa=0
FD=diff(F01);
La onda cuadrada valdrá 0 para flujos negativos y 1 para flujos positivos.
Para corregir errores de transformación a onda cuadrada se realiza la
diferencia de entre los picos de la onda cuadrada.
Figura 2: Representación ciclo respiratorio en onda cuadrada
4. Cálculo de índices de picos positivos y periodo de cada ciclo. Se aplica filtro
para eliminar valores de artefactos no relacionados al flujo de muestreo.
Martha Galarraga Página 2 | 5
Práctica 2
ind=find(FD>0); %indices de picos positivos
Tind=ind/fs; %periodods donde ocurren los indices
Period=diff(Tind); %diferencia entre los tiempos/periodo de cada
ciclo
m=median(Period);
for i=1:length(Period)
if (Period(i)<0.7*m)|| (Period(i)>1.7*m)
Period(i)=0;
end
end
Nzeros=Period==0;
Period(Nzeros)= [];
fr=1./Period*60;
T=(1:length(fr))'; %ciclos respiratorios
B) Aplicación de la función frec_ins para nuestros datos.
1. Carga de datos
Figura 3: Señales de flujo respiratorio de un individuo
2. Se aplica la función frec_ins para cada flujo
>> [fr0,T0]=frec_ins(f0);
>> [fr1,T1]=frec_ins(f1);
>> [fr2,T2]=frec_ins(f2);
>> [fr3,T3]=frec_ins(f3);
>> [fr4,T4]=frec_ins(f4);
3. Se realiza la concatenación de los periodos
>> T1=T1+T0(end);
>> T2=T2+T1(end);
>> T3=T3+T2(end);
>> T4=T4+T3(end);
4. Graficación de los resultados
>> plot(T0,fr0,T1,fr1,T2,fr2,T3,fr3,T4,fr4)
>> axis([0 220 8 22])
>> xlabel('Ciclo respiratorio')
>> ylabel('Frecuencia respiratoria (respiración/min)')
Martha Galarraga Página 3 | 5
Práctica 2
>> title('Frecuencia respiratoria durante incremento de estimulo de
hipercapnia')
Figura 4: frecuencia instantanea y ciclo respiratorio datos
Interpretación de resultados
En la figura 4 se pude observar el aumento de la frecuencia respiratoria
a medida que pasa el tiempo.
f0 que representa la frecuencia respiratoria en reposo es casi similar a f1
la cual representa la frecuencia respiratoria en el primer nivel de hipercapnia,
de lo que podemos observar no hay mayores riesgos en este nivel por lo que se
podría decir que es un nivel seguro de exposición al CO2.
En f2 se nota una reducción de ciclos respiratorios que nos indicaría
cierta dificultad para respirar y un aumento de la frecuencia respiratoria en un
intento de oxigenar la sangre, a partir de f3 ya podemos observar mayor
cantidad de ciclos respiratorios e incremento de la frecuencia respiratoria
indicando que el sujeto realiza un mayor esfuerzo respiratorio pero al ser los
ciclos tan cortos no da tiempo a una correcta absorción de O2 por lo que la
eliminación de CO2 en la sangre sería muy reducido.
En f4 la frecuencia respiratoria es alta y al igual que en f3 los ciclos son
muy cortos para los niveles de CO2 absorbidos por lo que empezaría a afectar al
sujeto de muchas maneras, ya sea con mareos, desmayos, arritmia, etc.
Martha Galarraga Página 4 | 5
Práctica 2
Con estos datos podemos concluir que a mayores niveles de CO2 menor
será la capacidad respiratoria del sujeto lo cual podría causar falta de oxígeno
que se distribuya en su torrente sanguíneo pudiendo llegar a causar desmayos
Figura 4: frecuencia instantanea y ciclo respiratorio datos
y en peor de los casos la muerte.
Martha Galarraga Página 5 | 5
