function [a,b] = cascadeIIRNotch(notchFrequencies, alpha, beta)

% y = IIRNotchFilter(soundfilename, zeroCoefficients, poleCoefficients, alpha)

%

% Computes the coefficients for a cascaded IIR digital notch filter. Each

% cascaded IIR notch filter is a second order filter. Thus, for each notch

% frequency specified a second order filter is added to the cascaded system.

%

% notchFrequencies A vector specifying the digital frequencies to be notched.

% alpha The magnitude of the poles

% beta The magnitude of the zeros

%

% a Coefficients of the numerator polynomial

% b Coefficient of the denominator polynomial

%

%

% Cascaded Digital Notch Filter Implementation

% Dr. Ratliff

% ECE505

%Calculate locations of each pole and zero, noting that each pole-zero

%also has a matching complex conjugates. All I am doing here is using the

%convolution command to multiply out the transfer function polynomials to

%obtain the coefficients.

w = notchFrequencies;

for k=1:length(w)

if(k==1)

Z = conv([ 1 -beta*exp(1i*w(k))], [ 1 -beta*exp(-1i*w(k))]);

P = conv([ 1 -alpha*exp(1i*w(k))], [ 1 -alpha*exp(-1i*w(k))]);

else

Z = conv(Z, [ 1 -beta*exp(1i*w(k))]);

Z = conv(Z, [ 1 -beta*exp(-1i*w(k))]);

P = conv(P, [ 1 -alpha*exp(1i*w(k))]);

P = conv(P, [ 1 -alpha*exp(-1i*w(k))]);

end

end

%Take the real part of the coefficients for use in the difference equation and display them

a = real(Z); b = real(P);

disp(sprintf(‘Cascaded Notch Filter Numerator (Zero) Coefficients’));

disp(a);

disp(sprintf(‘Cascaded Notch Filter Denominator (Pole) Coefficients’));

disp(b);ECE505 Digital Signal Processing

Laboratory Project – IIR Notch Filter Design

Dr. Ratliff

1 Project Definition:

In this project, you are given a digital audio clip that has been corrupted by sinusoidal tones resulting from

pressing the numeric keys on a touch-tone phone. The goal is to remove the tones from the signal by designing

cascaded digital notch filters. This project is designed to make use of the mathematical tools we have learned

this semester by applying them to a practical DSP problem.

2 Problem Description:

The Touch-Tone phone system, developed during the 1960’s by AT&T, was created to eliminate the need for

rotary-dialing telephone systems. This is accomplished using Dual-Tone Multi-Frequency (DTMF) signaling

and has resulted in the familiar tones we hear when dialing a telephone. As the name implies, each tone is

generated by introducing a sinusoidal wave at two distinct frequencies. As such, the DTMF telephone keypad

is laid out in a 4 × 4 matrix of push buttons in which each row represents the low frequency component and

each column represents the high frequency component of the DTMF signal. Pressing a key sends a combination

of the row and column frequencies, as described in the figure below:

For example, pressing the ‘1’ key produces a superposition of tones at 697 and 1209 hertz (Hz). When these

keys are pressed, they are mixed with the audio signal from the phone microphone. When such an audio stream

is digitized and recorded, it may be desirable to later remove these tones from the audio stream. One way to

do this is to design a group of filters such that each respective filter removes the content at a single frequency.

Such filters are called digital notch filters. A block diagram for such a system can be represented by:

In this case each hk[n], for k = 1, 2, . . . , 8, represents a digital notch filter for one of the numeric tone frequencies

from the DTMF keypad. Thus, if such an ideal filter bank is employed, it would indeed remove the tones resulting

from any key-press. Notice that the total impulse response of such a filter bank can be obtained as

h[n] = h1[n] ∗ h2[n] ∗ . . . ∗ h8[n]. (1)

1

In this project we are going to design second order IIR notch filters to remove each sinusoidal frequency term.

Thus, we expect to have two poles and two zeros for each filter. For notch filters of this type the design is straight-

forward. For each frequency we want to filter, we simply place a pole and zero along a line corresponding to

the digital frequency we wish to notch, along with their respective complex conjugate pairs, as shown in the

pole-zero diagram below.

Notice that the zero is placed upon the unit circle with a magnitude of 1, i.e., znk

= e±jθk , and each pole is

placed at the same angle, θk, only at a magnitude of α, i.e., zdk = αe±jθk , where 0 < α < 1, znk denotes system zeros and zdk denote system poles. Here θk corresponds directly to the digital frequency of the sinusoid to be notched. Later we

#### Why Choose Us

- 100% non-plagiarized Papers
- 24/7 /365 Service Available
- Affordable Prices
- Any Paper, Urgency, and Subject
- Will complete your papers in 6 hours
- On-time Delivery
- Money-back and Privacy guarantees
- Unlimited Amendments upon request
- Satisfaction guarantee

#### How it Works

- Click on the “Place Order” tab at the top menu or “Order Now” icon at the bottom and a new page will appear with an order form to be filled.
- Fill in your paper’s requirements in the "
**PAPER DETAILS**" section. - Fill in your paper’s academic level, deadline, and the required number of pages from the drop-down menus.
- Click “
**CREATE ACCOUNT & SIGN IN**” to enter your registration details and get an account with us for record-keeping and then, click on “PROCEED TO CHECKOUT” at the bottom of the page. - From there, the payment sections will show, follow the guided payment process and your order will be available for our writing team to work on it.