CODE Matlab Fourier Programme
clc;
clear all;
close all;
[x,Fs]=audioread('audio BTL.mp4'); %Tai file am thanh
X = fft(x); %Bien doi fourier khong nam o trung tam
N =length(x); %Xac dinh kich thuoc cua x
a = 2000*2*pi; %tan so cat
w0 =(-N/2+1:(N/2)); % Vector tan so trung tam
w = w0.*Fs/N; % LAY MOT TAN SO MAU
H = 1-(a./(a + (1i*w))); %H nam o trung tam
Hshift = fftshift(H);
Y = X .*Hshift' ; % loc tin hieu, chuyen vi de nhan 2 ma tran
y = real(ifft(Y)); % Bien doi fourier ngich
subplot(4,1,1)
plot(w,abs(fftshift(x)))
title('Am thanh ban dau - mien thoi gian');
subplot(4,1,2)
plot(w,abs(fftshift(X)))
title('Am thanh ban dau - mien tan so');
subplot(4,1,3)
plot(w,abs(fftshift(Y)))
title('Am thanh khu nhieu - mien ttan so');
subplot(4,1,4)
plot(w,abs(fftshift(y)))
title('Am thanh khu nhieu - mien thoi gian');