Scribd
Upload
46 views2 pages

Function: (Out) BPF (f1, f2, In)

This document describes a function called bpf that designs and implements a band-pass filter. It takes in lower and upper cutoff frequencies (f1 and f2) to design a Chebyshev band-pass filter. It can optionally take an input signal and return the band-pass filtered output. If no output is provided, it will plot the magnitude and phase response of the designed filter.

Uploaded by

fikrirozi
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
46 views2 pages

Function: (Out) BPF (f1, f2, In)

This document describes a function called bpf that designs and implements a band-pass filter. It takes in lower and upper cutoff frequencies (f1 and f2) to design a Chebyshev band-pass filter. It can optionally take an input signal and return the band-pass filtered output. If no output is provided, it will plot the magnitude and phase response of the designed filter.

Uploaded by

fikrirozi
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

function

[out] = bpf(f1,f2,in)

% BPF ......... Band-pass filter design and filtering of the input sequence.
%
%
BPF(F1,F2) specifies a Chebychev band-pass filter with lower cut-off
%
frequency = min(F1,F2) and upper cut-off frequency = max(F1,F2).
%
Here, 0 < F1, F2 < Fs/2, where Fs is the sampling frequency.
%
Magnitude and phase responses are plotted.
%
%
[Y] = BPF(F1,F2,X) returns vector Y obtained through band-pass
%
filtering the input sequence X.
%
%
See also LPF, FREQZ and FILTER.
%
%
%
%
%
%
%

AUTHORS : M. Zeytinoglu & N. W. Ma


Department of Electrical & Computer Engineering
Ryerson Polytechnic University
Toronto, Ontario, CANADA
DATE
: August 1991.
VERSION : 1.0

%===========================================================================
% Modifications history:
% ---------------------%
o
Added "checking" 11.30.1992 MZ
%
o
Tested (and modified) under MATLAB 4.0/4.1 08.16.1993 MZ
%
o
Modified under OCTAVE 2.0.14 2000.08.12
%===========================================================================
global
global
global
global
global
global

START_OK;
SAMPLING_CONSTANT;
SAMPLING_FREQ;
BINARY_DATA_RATE;
BELL;
WARNING;

check;
if ((nargin < 2) | (nargin > 3))
error(eval('eval(BELL),eval(WARNING),help bpf'));
return;
end
%---------------------------------------------------------------------------%
set the sampling frequency fs, lower cut-off frequency fl, and
%
the upper cut-off frequency fu.
%---------------------------------------------------------------------------fs = SAMPLING_FREQ;
fl = min(f1,f2);
fu = max(f1,f2);

% Default sampling frequency.

%---------------------------------------------------------------------------%
check the consistency of fl, fu and fs.

%---------------------------------------------------------------------------if ( fl < 0 )
error('Lower cut-off frequency must be positive')
end
if ( fu >= fs/2 )
fprintf('Upper cut-off frequency must be less than %6.2f
[kHz].\n',SAMPLING_FREQ/2000);
error('');
end
%---------------------------------------------------------------------------%
determine filter coefficients
%---------------------------------------------------------------------------passband
ripple
filt_order
[Bc,Ac]

=
=
=
=

[fl/(fs/2) fu/(fs/2)];
% Passband specification
.1;
% Allowable ripple, in decibels
8;
cheby1(filt_order, ripple, passband);

%---------------------------------------------------------------------------%
Output routines
%
===============
%
If no output arguments, plot the magnitude response, otherwise
%
compute the channel output vector "y"
%---------------------------------------------------------------------------if (nargin == 2)
fpts
f
Hc
mag
phase

=
=
=
=
=

256;
fs/(2*fpts) * (0:fpts-1);
freqz(Bc,Ac,fpts);
abs(Hc);
angle(Hc);

title('BPF magnitude response'),


xlabel('Frequency (Hz)'),
ylabel('Mag [dB]'), grid;
axis([0 fs/2 -100 0]);
subplot(211), plot( f, 20*log10(mag) );

...

title('BPF phase response'),


xlabel('Frequency (Hz)'),
ylabel('Phase'), grid;
axis;
subplot(212), plot( f, unwrap(phase) );

...

else
out = filter(Bc,Ac,in);
end

...
...
...

...
...
...

You might also like