Welcome to EDAboard.com

Welcome to our site! EDAboard.com is an international Electronics Discussion Forum focused on EDA software, circuits, schematics, books, theory, papers, asic, pld, 8051, DSP, Network, RF, Analog Design, PCB, Service Manuals... and a whole lot more! To participate you need to register. Registration is free. Click here to register now.

How to implement the filter give in Z transform using matlab

Status
Not open for further replies.

rrajbe

Newbie level 1
Joined
Jan 27, 2010
Messages
1
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Location
Tamilnadu/India
Activity points
1,300
Good morning to every one. . .

I am very new to this forum and not much expert in DSP field too. . .

I am trying to implement a IEEE paper named "A New Delayless Subband Adaptive Filtering Algorithm for Active Noise Control Systems" published in IEEE TRANSACTIONS ON AUDIO, SPEECH, AND LANGUAGE PROCESSING july 2009.


I have not done any such projects before and so i am strugling where to start and how to implement the given details.

This project is all about the following

1) A full-band filter that filters the input signal.
2) Decomposition of input and error signals into subbands.
3) Decimation in subbands.
4) Adaptive filters working in subbands.
5) A weight stacking method to combine all subbands weights into a full-band filter.


These are the details of a filter specification given in the paper. Could any one please help me to implement this in the matlab please?

x(n) is the input signal which is to be filtered into M subbands named

X'k(n)

where K ranges form 0 to M-1, using a analysis filter bank h(z) with a decimation factor D.

Filter h(z) is given by

h(z)=[ H0(z) ]
| H1(z) |
| . |
| . |
| . |
[ Hm-1(z)]

Where Hk(z) is the transfer function of Kth subband.


Given the low pass prototype filter for filter banks is

H0(z)=1+(Z)^-1+(Z)^-2+.......+(Z)^-M+1


The resulting filter bank is the simplest FIR filter perfect reconstruction filter bank which is made by


Hk(z)=H0(z.exp(-j*2*pi*k/M))

with a frequency response of

H0(exp(jw))={ M , when w=0;
{exp(-j*w(M-1)/2sin(wM/2)/sin(w/2)) , otherwise;


M is given about to be 16

D=M/4

For this proposed UDFT modulated filter banks, h(z) is defined by


h(z) = (1/M)F*[1 , (Z)^-1 , (Z)^-2 ,. . ., (Z)^-M+1] (transpose)

Where F is the DFT matix of order M. The central frequency of bandpass filters Hk(Z) are located at

wk=2*pi*k/M for 0 ≤ k ≤ M-1.

This filter can be realized by tapped delay line of length M followed by an inverse FFT block.
 

divyasagar

Newbie level 3
Joined
Jun 10, 2011
Messages
3
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,292
did u get a solution? that will be of great help to me as my problem is also similar
 

Status
Not open for further replies.

Similar threads

Part and Inventory Search

Welcome to EDABoard.com

Sponsor

Top