20 lines
513 B
Mathematica
20 lines
513 B
Mathematica
|
function output_y = DarellbandpassFilter(y,Fs,LOW,MED,HIGH)
|
||
|
|
Len = length(y);
|
||
|
|
F = Fs * (-Len/2 : (Len/2 - 1))/Len ;
|
||
|
|
Mod_Freq = fftshift(fft(y));
|
||
|
|
lenf = length(F);
|
||
|
|
output = 0 .* Mod_Freq; % zero array of len f
|
||
|
|
|
||
|
|
for n = 1:lenf
|
||
|
|
if ((LOW < abs(F(n))) && HIGH > abs(F(n)))
|
||
|
|
output(n) = 1;
|
||
|
|
else
|
||
|
|
output(n) = 0;
|
||
|
|
end
|
||
|
|
end
|
||
|
|
|
||
|
|
filtered_Mod_Freq = fftshift(Mod_Freq .* output);
|
||
|
|
|
||
|
|
output_y = real(ifft(filtered_Mod_Freq));
|
||
|
|
|
||
|
|
end
|