How to Design a Fir Filter Using Matlab
AIM : Design of FIR filters of Low pass and high pass filter using Matlab commands
EQUIPMENTS:
Operating System - Windows XP
Constructor - Simulator
Software - CCStudio 3 & MATLAB 7.5
Description: Digital filters refers to the hard ware and software implementation of the mathematical algorithm which accepts a digital signal as input and produces another digital signal as output whose wave shape, amplitude and phase response has been modified in a specified manner. Digital filter play very important role in DSP. Compare with analog filters they are preferred in number of application due to following advantages.
? Truly linear phase response
? Better frequency response
? Filtered and unfiltered data remains saved for further use.
There are two type of digital filters.
? FIR (finite impulse response) filter
? IIR (infinite impulse response) filter
Description Of The Commands Used In FIR Filter Design
FIR1: FIR filters design using the window method.
B = FIR1(N,Wn) designs an N'th order low
pass FIR digital filter and returns the filter coefficients in length N+1 vector B. The cut-
off frequency Wn must be between 0 < Wn < 1.0, with 1.0 corresponding to half the
sample rate. The filter B is real and has linear phase. The normalized gain of the filter
at Wn is -6 dB.
B = FIR1(N,Wn,'high') designs an N'th order highpass filter.
You can also use B = FIR1(N,Wn,'low') to design a lowpass filter.
If Wn is a two-element vector, Wn = [W1 W2], FIR1 returns an order N bandpass filter with passband W1 < W < W2.
B = FIR1(N,Wn,'stop') is a bandstop filter if Wn = [W1 W2]. You can also specify If Wn is a multi-element vector, Wn = [W1 W2 W3 W4 W5 ... WN], FIR1 returns an order N multiband filter with bands 0 < W < W1, W1 < W < W2, ..., WN < W < 1.
B = FIR1(N,Wn,'DC-1') makes the first band a passband.
B = FIR1(N,Wn,'DC-0') makes the first band a stopband.
By default FIR1 uses a Hamming window. Other available windows, including Boxcar, Hann, Bartlett, Blackman, Kaiser and Chebwin can be specified with an optional trailing argument.
For example, B = FIR1(N,Wn,kaiser(N+1,4)) uses a Kaiser window with beta=4.
B = FIR1(N,Wn,'high', chebwin(N+1,R)) uses a Chebyshev window.
For filters with a gain other than zero at Fs/2, e.g., highpass and bandstop filters, N must
be even. Otherwise, N will be incremented by one. In this case the window length
should be specified as N+2.
By default, the filter is scaled so the center of the first pass band has magnitude exactly one after windowing. Use a trailing 'noscale' argument to prevent this scaling, e.g.
B = FIR1(N,Wn,'noscale')
B = FIR1(N,Wn,'high','noscale')
B = FIR1(N,Wn,wind,'noscale').
You can also specify the scaling explicitly, e.g. FIR1(N,Wn,'scale'), etc.
FREQZ Digital Filter Frequency Response.
[H,W] = FREQZ(B,A,N) returns the N-point complex frequency response vector H and
the N-point frequency vector W in radians/sample of the filter: given numerator and
denominator coefficients in vectors B and A. The frequency response is evaluated at N
points equally spaced around the upper half of the unit circle. If N isn't specified, it
defaults to 512.
[H,W] = FREQZ(B,A,N,'whole') uses N points around the whole unit circle.
H = FREQZ(B,A,W) returns the frequency response at frequencies designated in vector W, in radians/sample (normally between 0 and pi).
[H,F] = FREQZ(B,A,N,Fs) and [H,F] = FREQZ(B,A,N,'whole',Fs) return frequency vector F (in Hz), where Fs is the sampling frequency (in Hz).
H = FREQZ(B,A,F,Fs) returns the complex frequency response at the frequencies designated in vector F (in Hz), where Fs is the sampling frequency (in Hz).
[H,W,S] = FREQZ(...) or [H,F,S] = FREQZ(...) returns plotting information to be used with FREQZPLOT. S is a structure whose fields can be altered to obtain different frequency response plots. For more information see the help for FREQZPLOT.
FREQZ(B,A,...) with no output arguments plots the magnitude and unwrapped phase of the filter in the current figure window.
Designing A Low Pass Filter :
Program:
%Suppose the target is to pass all frequencies below 1200 Hz
fs=8000; % sampling frequency
n=50; % order of the filter
w=1200/ (fs/2);
b=fir1(n,w,'low'); % Zeros of the filter
freqz(b,1,128,8000); % Magnitude and Phase Plot of the filter figure(2)
[h,w]=freqz(b,1,128,8000);
subplot(2,1,1);
plot(w,abs(h));% Normalized Magnitude Plot
title('Normalized Magnitude Plot');
grid
subplot(2,1,2);
figure(2)
zplane(b,1); % the plot in lab
title('zplane');
Output:
Designing High Pass Filter:
%Now the target is to pass all frequencies above 1200 Hz
fs=8000;
n=50;
w=1200/ (fs/2);
b=fir1(n,w,'high');
freqz(b,1,128,8000); % this function plots the phase(degree)and magnitude in db
subplot(2,1,2)
figure(2)
[h,w]=freqz(b,1,128,8000);
plot(w,abs(h)); % Normalized Magnitude Plot
title('Magnitude Plot ');
grid
figure(3)
zplane(b,1); % this function plots fiq in zplane
Output:
Designing High Pass Filter:
fs=8000;
n=50;
w=1200/ (fs/2);
b=fir1(n,w,'high');
freqz(b,1,128,8000);
figure(2)
[h,w]=freqz(b,1,128,8000);
plot(w,abs(h)); % Normalized Magnitude Plot
grid
figure(3)
zplane(b,1);
Designing Band Pass Filter:
fs=8000;
n=40;
b=fir1(n,[1200/4000 1800/4000],'bandpass'); freqz(b,1,128,8000)
figure(2)
[h,w]=freqz(b,1,128,8000);
plot(w,abs(h)); % Normalized Magnitude Plot
grid
figure(3)
zplane(b,1);
Designing Notch Filter
fs=8000;
n=40;
b=fir1(n,[1500/4000 1550/4000],'stop'); freqz(b,1,128,8000)
figure(2)
[h,w]=freqz(b,1,128,8000);
plot(w,abs(h)); % Normalized Magnitude Plot
grid
figure(3)
zplane(b,1);
Designing Multi Band Filter
n=50;
w=[0.2 0.4 0.6];
b=fir1(n,w);
freqz(b,1,128,8000) figure(2)
[h,w]=freqz(b,1,128,8000);
plot(w,abs(h)); % Normalized Magnitude Plot
grid
figure(3)
zplane(b,1);
RESULT: The FIR low pass & high pass filter for given values is obtained. Hence the ideal and practical response of FIR filter are proved.
VIVA QUESTION:
- What is filter?
- What is FIR and IIR filter define, and distinguish between these two?
- What is window method? How you will design an FIR filter using window method?
- What are low-pass and band-pass filter and what is the difference between these two?
- What is the matlab command for Hamming window? Explain.
- What do you mean by built in function 'abs' and where it is used?
- Explain how the FIR filter are stable?
- Why is the impulse response "finite"?
- What does "FIR" mean?
- What are the advantages of FIR Filters (compared to IIR filters)?
- What are the disadvantages of FIR Filters (compared to IIR filters)?
- What terms are used in describing FIR filters?
- What is the delay of a linear-phase FIR?
- What is the Z transform of a FIR filter?
- What is the frequency response formula for a FIR filter?
- How Can I calculate the frequency response of a FIR using the Discrete Fourier Transform (DFT)?
- What is the DC gain of a FIR filter?
How to Design a Fir Filter Using Matlab
Source: https://vikramlearning.com/jntuh/notes/digital-signal-processing-lab/design-of-fir-filters-of-%C2%A0low-pass-and-high-pass-filter-%C2%A0using-matlab-commands/374
0 Response to "How to Design a Fir Filter Using Matlab"
Post a Comment