25 lines
332 B
Matlab
25 lines
332 B
Matlab
|
|
[x, Fs] = audioread('wannabe_in_la.wav');
|
|
ft_x = fft(x);
|
|
|
|
Fc = 300;
|
|
p = Fs/(2*Fc);
|
|
|
|
H_f = zeros(length(ft_x), 1);
|
|
do = length(H_f)/(2*pi);
|
|
for k = 1:length(H_f)
|
|
if k < do*(pi/p)
|
|
H_f(k) = 1;
|
|
end
|
|
if k > do*(2*pi - pi/p)
|
|
H_f(k) = 1;
|
|
end
|
|
end
|
|
|
|
tic;
|
|
ft_x = fft(x);
|
|
|
|
ft_y = ft_x .* H_f;
|
|
|
|
y = ifft(ft_y, 'symmetric');
|
|
toc; |