Muy buen video, solo una duda. En el primer código de pasa bajas ¿dónde o cómo coloco la frecuencia de core? Tengo una señal de menos de 300Hz y quiero eliminar el ruido ¿se hace en Wn? si es así cómo sé qué valor poner porque dicen que es de 0 a 1.
Saludos, ese error indica que no tienes el archivo de la señal a filtrar ECG.txt en la carpeta de acceso principal (PATH de MATLAB). La autora del video dejó un link en la descripción donde está el archivo.
hola buenos dias tengo el siguiente programa para procesar un audio con filtro fir2,al audio de musica le puse ruido para quitarselo con el filtro fir pero no lo quitarle del todo el rudo y q se vea en las dos graficas el antes de filtrala y el despues del filtrado,me puedes ayudar? clear all close all clc % IMPORTACION DE AUDIO [a, fs] = audioread("AC-DC.wav"); d = length(a)/fs; % duracion de la cancion a_m = 0.5*(a(:,1)+a(:,2)).'; % audio en un solo canal t = linspace(0,d,length(a_m)); %AGREGAR RUIDO a_r = a_m + randn(1,length(t))*0.05; %FILTRO FIR n =1500; f = [0 0.23 0.23 1]; %[0 0.6 0.6 1]; m = [1 1 0 0]; b1 = fir2(n,f,m); %Crea filtro FIR2 a_F= filter(b1,1,a_r); %señal de audio filtrada figure(2) subplot(3,1,1) plot(t,a_r ) hold on subplot(3,1,1) plot(t,a_F,'r') legend("señal de audio","señal Filtrada"); title(" señal de audio y señal filtrada "); xlabel("tiempo [s]"); ylabel("Amplitud"); [h1,w1] = freqz(b1,1,fs); figure(2) subplot(3,1,2) plot(f,m,w1/pi,abs(h1)) legend("ideal","fir2"); xlabel("\omega/ \pi"); %ESPECTRO %FRECUENCIA _TRANFORMADA FFT PARA PASAR DEL DOMINIO DEL TIEMPO A FRECUENCIA A_r = fftshift(fft(a_r)); %transformada rapida de furier de nuestro audio ffthshift para dejar el espectro unido f = linspace(-fs/2,fs/2,length(A_r)); mag_A_r = abs(A_r); A_F= fftshift(fft(a_F)); %transformada rapida de furier de nuestro audio ffthshift para dejar el espectro unido f = linspace(-fs/2,fs/2,length(A_F)); mag_A_F = abs(A_F); figure(2); subplot(3,1,3) plot(f,mag_A_r/max(mag_A_r)); hold on subplot(3,1,3) plot(f,mag_A_F/max(mag_A_F),'r') legend("Espectro de audio con ruido","Espectro de audio con ruido Filtrado"); title(" espectro de audio y Espectro filtrado"); xlabel("Frecuencia [Hz]"); ylabel("Amplitud"); ax = gca; % se usa para dejar x10^3 ax.XAxis.Exponent = 3; sound(a_F ,fs); pause(d+1); sound(a_r,fs)
Hola, muy buen video. Gracias por compartir.
Hola. Se agradece mucho la explicación y haber dejado los archivos para descargar. Saludos.
Excelente vídeo, muchas gracias por compartir tu conocimiento, ¿cuando subes el siguiente vídeo?
Gracias! En el enlace de la descripción está el link para descargar los archivos.
Excelente explicación, gracias por tu gran aporte y enseñanza, ¿Dictas clases virtuales?
Muchas gracias 👍
Excelente!
bueas, de donde sale el txt? o cual seria la senal en ese caso?
Muy buen video, solo una duda.
En el primer código de pasa bajas ¿dónde o cómo coloco la frecuencia de core?
Tengo una señal de menos de 300Hz y quiero eliminar el ruido ¿se hace en Wn? si es así cómo sé qué valor poner porque dicen que es de 0 a 1.
Genial. ¡Qué buenos videos!
Por cierto, ¿qué estudiaste? 🤔
Veo que eres muy conocedora en el tema de procesamiento digital de señales
Ingeniería Biomédica, Maestría en Bioelectrónica y Doctorado en curso en Ingeniería Aplicada
@@mundotecnologico6659 Quede asombrado con tu respuesta! Mi admiración completa!!! Gracias su información me a ayudado mucho!
¡@@mundotecnologico6659, felicitaciones! . Eres un ejemplo para muchas personas que con esfuerzo se logran objetivos complicados 🙂
Excelente vídeo, pero no comprendo bien cómo calculaste los valores de f y m para el segundo filtro (fir2)?
Como agrego la paqueteria EGC , para que me ejute el comando load , lo intente hacer en MATLAB pero no me corre
no complila me dice
Error using load
Unable to read file 'ECG.txt'. No such file or directory.
Saludos, ese error indica que no tienes el archivo de la señal a filtrar ECG.txt en la carpeta de acceso principal (PATH de MATLAB). La autora del video dejó un link en la descripción donde está el archivo.
hola buenos dias tengo el siguiente programa para procesar un audio con filtro fir2,al audio de musica le puse ruido para quitarselo con el filtro fir pero no lo quitarle del todo el rudo y q se vea en las dos graficas el antes de filtrala y el despues del filtrado,me puedes ayudar?
clear all
close all
clc
% IMPORTACION DE AUDIO
[a, fs] = audioread("AC-DC.wav");
d = length(a)/fs; % duracion de la cancion
a_m = 0.5*(a(:,1)+a(:,2)).'; % audio en un solo canal
t = linspace(0,d,length(a_m));
%AGREGAR RUIDO
a_r = a_m + randn(1,length(t))*0.05;
%FILTRO FIR
n =1500;
f = [0 0.23 0.23 1]; %[0 0.6 0.6 1];
m = [1 1 0 0];
b1 = fir2(n,f,m); %Crea filtro FIR2
a_F= filter(b1,1,a_r); %señal de audio filtrada
figure(2)
subplot(3,1,1)
plot(t,a_r )
hold on
subplot(3,1,1)
plot(t,a_F,'r')
legend("señal de audio","señal Filtrada");
title(" señal de audio y señal filtrada ");
xlabel("tiempo [s]");
ylabel("Amplitud");
[h1,w1] = freqz(b1,1,fs);
figure(2)
subplot(3,1,2)
plot(f,m,w1/pi,abs(h1))
legend("ideal","fir2");
xlabel("\omega/ \pi");
%ESPECTRO
%FRECUENCIA _TRANFORMADA FFT PARA PASAR DEL DOMINIO DEL TIEMPO A FRECUENCIA
A_r = fftshift(fft(a_r)); %transformada rapida de furier de nuestro audio ffthshift para dejar el espectro unido
f = linspace(-fs/2,fs/2,length(A_r));
mag_A_r = abs(A_r);
A_F= fftshift(fft(a_F)); %transformada rapida de furier de nuestro audio ffthshift para dejar el espectro unido
f = linspace(-fs/2,fs/2,length(A_F));
mag_A_F = abs(A_F);
figure(2);
subplot(3,1,3)
plot(f,mag_A_r/max(mag_A_r));
hold on
subplot(3,1,3)
plot(f,mag_A_F/max(mag_A_F),'r')
legend("Espectro de audio con ruido","Espectro de audio con ruido Filtrado");
title(" espectro de audio y Espectro filtrado");
xlabel("Frecuencia [Hz]");
ylabel("Amplitud");
ax = gca; % se usa para dejar x10^3
ax.XAxis.Exponent = 3;
sound(a_F ,fs);
pause(d+1);
sound(a_r,fs)
la técnica de filtro de Gabor no lineal tienes algún video donde explique como utilizarlo. gracias
como podemos obtener los coeficientes del filtro FIR con respuesta en la frecuencia? ahi los utilizamos pero como puedo verlos en la pantalla?
Se guardan en las variables a y b. Sería cuestión de graficarlas
y ese .txt de donde se saca?
La autora del video dejó un link en la descripción donde está el archivo ECG.txt, tambien puedes entrar a Physionet
Podrías proporcionarme un nombre para referenciar tus videos?
María Corzo