%FS=FRECUENCIA DE MUESTREO
[y fs]=audioread('[Link]'); %señal del audio
y1=y(:,1);
sound(y1,fs); %SONIDO DEL AUDIO
figure(1) % nombre DE LA grafica de la señal de
audio
plot(y1) %grafica señal del audio
pause on %FUNCION PAUSA ACTIVADA
pause(3) % 3 SEGUNDOS DE ESPERA
fo = 4000; % frecuencia de 4khz
n = 0:length(y1)-1;
x = sin(2*pi*(fo/fs)*n); %tono de 4khz
sound(x,fs) %reproduce sonido con la frecuencia de
muestreo
figure(2) % nombre grafica de la señal de tono de
4khz
plot(x) % grafica de la señal de tono de 4khz
pause on %FUNCION PAUSA ACTIVADA
pause(3) % 3 SEGUNDOS DE ESPERA
y2=y1+x'; % señal de audio + frecuencia de 4khz
sound(y2,fs) %reproduce sonido de la señal
combinada con la frecuencia de muestreo
figure(3) % NOMBRE grafica de la señal de audio +
señal de tono 4khz
plot(y2) %grafica de la señal de audio + señal
de tono 4khz
pause on %FUNCION PAUSA ACTIVADA
pause(3) % 3 SEGUNDOS DE ESPERA
%Filtro notch
F1=3880/(fs/2) %frecuencia que quiero eliminar
F2=0.50 %Factor de calidad del filtro
[b,a] = iirnotch(F1,F2); %variables del filtro
filtrada=filter(b,a,y2);
sound(filtrada,fs) %sonido de la señal filtrada
figure (4) % figura de la señal filtrada
plot(filtrada); %grafica de la señal filtrada no al
100%