Filter Design Techniques

Similar documents
Quantitative Genomics and Genetics BTRY 4830/6830; PBSB

A general N-dimensional vector consists of N values. They can be arranged as a column or a row and can be real or complex.

Minimum Spanning Trees

(A) the function is an eigenfunction with eigenvalue Physical Chemistry (I) First Quiz

LINEAR SYSTEMS THEORY

Convergence Theorems for Two Iterative Methods. A stationary iterative method for solving the linear system: (1.1)

Fractions. Mathletics Instant Workbooks. Simplify. Copyright

Using SIP techniques to verify the trade-off between SNR and information capacity of a sigma delta modulator

Preview. Graph. Graph. Graph. Graph Representation. Graph Representation 12/3/2018. Graph Graph Representation Graph Search Algorithms

8. Linear Contracts under Risk Neutrality

TOPIC 5: INTEGRATION

INTEGRALS. Chapter 7. d dx. 7.1 Overview Let d dx F (x) = f (x). Then, we write f ( x)

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: Derivation of Ideal MHD Equation

ERDOS-SMARANDACHE NUMBERS. Sabin Tabirca* Tatiana Tabirca**

The Fourier Transform

Floating Point Number System -(1.3)

Floating Point Number System -(1.3)

Walk Like a Mathematician Learning Task:

Shortest Paths in Graphs. Paths in graphs. Shortest paths CS 445. Alon Efrat Slides courtesy of Erik Demaine and Carola Wenk

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

A Solution for multi-evaluator AHP

Multi-Section Coupled Line Couplers

Why the Junction Tree Algorithm? The Junction Tree Algorithm. Clique Potential Representation. Overview. Chris Williams 1.

INTERPOLATION(1) ELM1222 Numerical Analysis. ELM1222 Numerical Analysis Dr Muharrem Mercimek

Chapter Newton-Raphson Method of Solving a Nonlinear Equation

Lecture contents. Bloch theorem k-vector Brillouin zone Almost free-electron model Bands Effective mass Holes. NNSE 508 EM Lecture #9

Generalized Half Linear Canonical Transform And Its Properties

b.) v d =? Example 2 l = 50 m, D = 1.0 mm, E = 6 V, " = 1.72 #10 $8 % & m, and r = 0.5 % a.) R =? c.) V ab =? a.) R eq =?

Wave Phenomena Physics 15c

SAMPLE LITANY OF THE SAINTS E/G. Dadd9/F. Aadd9. cy. Christ, have. Lord, have mer cy. Christ, have A/E. Dadd9. Aadd9/C Bm E. 1. Ma ry and. mer cy.

Lecture 20: Minimum Spanning Trees (CLRS 23)

Lecture 3: Phasor notation, Transfer Functions. Context

Fundamentals of Continuum Mechanics. Seoul National University Graphics & Media Lab

1 Introduction to Modulo 7 Arithmetic

Improving Union. Implementation. Union-by-size Code. Union-by-Size Find Analysis. Path Compression! Improving Find find(e)

Lecture 11 Waves in Periodic Potentials Today: Questions you should be able to address after today s lecture:

Instructions for Section 1

September 27, Introduction to Ordinary Differential Equations. ME 501A Seminar in Engineering Analysis Page 1. Outline

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

1. Stefan-Boltzmann law states that the power emitted per unit area of the surface of a black

5/1/2018. Huffman Coding Trees. Huffman Coding Trees. Huffman Coding Trees. Huffman Coding Trees. Huffman Coding Trees. Huffman Coding Trees

Rate of Molecular Exchange Through the Membranes of Ionic Liquid Filled. Polymersomes Dispersed in Water

Vowel package manual

CONTINUITY AND DIFFERENTIABILITY

The Mathematics of Harmonic Oscillators

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

HIGHER ORDER DIFFERENTIAL EQUATIONS

FREE VIBRATION ANALYSIS OF FUNCTIONALLY GRADED BEAMS

Cycles and Simple Cycles. Paths and Simple Paths. Trees. Problem: There is No Completely Standard Terminology!

Spanning Tree. Preview. Minimum Spanning Tree. Minimum Spanning Tree. Minimum Spanning Tree. Minimum Spanning Tree 10/17/2017.

CS September 2018

Quantum Mechanics & Spectroscopy Prof. Jason Goodpaster. Problem Set #2 ANSWER KEY (5 questions, 10 points)

I. The Connection between Spectroscopy and Quantum Mechanics

Ch 1.2: Solutions of Some Differential Equations

(2) If we multiplied a row of B by λ, then the value is also multiplied by λ(here lambda could be 0). namely

Case Study VI Answers PHA 5127 Fall 2006

SAMPLE CSc 340 EXAM QUESTIONS WITH SOLUTIONS: part 2

Section 3: Antiderivatives of Formulas

??? Dynamic Causal Modelling for M/EEG. Electroencephalography (EEG) Dynamic Causal Modelling. M/EEG analysis at sensor level. time.

1- I. M. ALGHROUZ: A New Approach To Fractional Derivatives, J. AOU, V. 10, (2007), pp

CSE303 - Introduction to the Theory of Computing Sample Solutions for Exercises on Finite Automata

Chapter 5 Transient Analysis

10. The Discrete-Time Fourier Transform (DTFT)

ME 522 PRINCIPLES OF ROBOTICS. FIRST MIDTERM EXAMINATION April 19, M. Kemal Özgören

b. How many ternary words of length 23 with eight 0 s, nine 1 s and six 2 s?

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

UNIT # 08 (PART - I)

Grand Canonical Ensemble

ECE COMBINATIONAL BUILDING BLOCKS - INVEST 13 DECODERS AND ENCODERS

In which direction do compass needles always align? Why?

Dennis Bricker, 2001 Dept of Industrial Engineering The University of Iowa. MDP: Taxi page 1

The University of Sydney MATH 2009

Planar Upward Drawings

, each of which is a tree, and whose roots r 1. , respectively, are children of r. Data Structures & File Management

1.9 Cartesian Tensors

Introduction to Numerical Integration Part II

A Note on Estimability in Linear Models

Discrete Shells Simulation

Partially Observable Systems. 1 Partially Observable Markov Decision Process (POMDP) Formalism

Weighted Graphs. Weighted graphs may be either directed or undirected.

2. Finite Impulse Response Filters (FIR)

FSA. CmSc 365 Theory of Computation. Finite State Automata and Regular Expressions (Chapter 2, Section 2.3) ALPHABET operations: U, concatenation, *

3.4 Properties of the Stress Tensor

Economics 600: August, 2007 Dynamic Part: Problem Set 5. Problems on Differential Equations and Continuous Time Optimization

More Foundations. Undirected Graphs. Degree. A Theorem. Graphs, Products, & Relations

Garnir Polynomial and their Properties

, between the vertical lines x a and x b. Given a demand curve, having price as a function of quantity, p f (x) at height k is the curve f ( x,

Errata for Second Edition, First Printing

Single Correct Type. cos z + k, then the value of k equals. dx = 2 dz. (a) 1 (b) 0 (c)1 (d) 2 (code-v2t3paq10) l (c) ( l ) x.

The Hyperelastic material is examined in this section.

12/3/12. Outline. Part 10. Graphs. Circuits. Euler paths/circuits. Euler s bridge problem (Bridges of Konigsberg Problem)

FOURIER SERIES. Series expansions are a ubiquitous tool of science and engineering. The kinds of

5/9/13. Part 10. Graphs. Outline. Circuits. Introduction Terminology Implementing Graphs

Having a glimpse of some of the possibilities for solutions of linear systems, we move to methods of finding these solutions. The basic idea we shall

Limits Indeterminate Forms and L Hospital s Rule

Chapter 6 Student Lecture Notes 6-1

8. INVERSE Z-TRANSFORM

The Schur-Cohn Algorithm

VIBRATION CONTROL OF SMART COMPOSITE SHELL STRUCTURES USING GENETIC ALGORITHM

TIME MINIMIZING TRANSPORTATION PROBLEM WITH FRACTIONAL BOTTLENECK OBJECTIVE FUNCTION

Transcription:

Fltr Dsgn chnqus Fltr Fltr s systm tht psss crtn frquncy componnts n totlly rcts ll othrs Stgs of th sgn fltr Spcfcton of th sr proprts of th systm ppromton of th spcfcton usng cusl scrt-tm systm Rlzton of th systm Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 1

Rvw of scrttm systms Frquncy rspons : proc : pro for rl mpuls rspons h[] Mgntu rspons H s vn functon Phs rspons s o functon mpl : H 1.5-4 - 4 5 yqust frquncy..., 1, 1,1, 1,1,... -5-4 - 4 Fltr Dsgn-FIR cwlu@twns..nctu.u.tw

Rvw of scrttm systms `Populr frquncy rsponss for fltr sgn : low-pss P hgh-pss HP bn-pss BP bn-stop mult-bn Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 3

Rvw of scrttm systms FIR fltrs fnt mpuls rspons: H z B z 1 b + b z +... 1 z + b z Movng vrg fltrs M fltrs pols t th orgn z hnc gurnt stblty zros zros of Bz, ll zro fltrs corrspons to ffrnc quton y[ ] b u [ ] + b1u [ 1] +... + b u [ mpuls rspons h[ ] b, h[1] b1,..., h[ ] b, h[ 1] + ],... Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 4

nr Phs FIR Fltrs on-cusl zro-phs fltrs : mpl: symmtrc mpuls rspons h[-],.h[-1],h[],h[1],...,h[] h[]h[-], 1.. frquncy rspons s H +. h[ ]... cos -.. rl-vlu zro-phs trnsfr functon - cusl mplmntton by ntroucng group ly Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 5

nr Phs FIR Fltrs Cusl lnr-phs fltrs non-cusl zro-phs + ly mpl: symmtrc mpuls rspons & vn h[],h[1],.,h[] vn h[]h[-],.. frquncy rspons s H. h[ ]... cos.. cusl mplmntton of zro-phs fltr, by ntroucng group ly z z Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 6

nr Phs FIR Fltrs yp-1 yp- yp-3 yp-4 vn +1o vn +1o symmtrc symmtrc nt-symmtrc nt-symmtrc h[]h[-] h[]h[-] h[]-h[-] h[]-h[-] / cos / cos cos / sn 1 cos /. sn cos P/HP/BP P/BP zro t zro t, zro t HP Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 7

nr Phs FIR Fltrs ffcnt rct-form rlzton. mpl: u[] + Δ Δ Δ Δ Δ Δ Δ Δ + + + + bo b1 b b3 b4 y[] + + + + PS: IIR fltrs cn EVER hv lnr-phs proprty! Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 8

Fltr Spcfcton E: ow-pss 1. 1 δ P 1.8 Pssbn Rppl P S.6 Pssbn Cutoff -> <- Stopbn Cutoff.4. Stopbn Rppl δ S.5 1 1.5.5 3 Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 9

Fltr Dsgn Problm Dsgn of fltrs s problm of functon ppromton For FIR fltr, t mpls polynoml ppromton For IIR fltr, t mpls ppromton by rtonl functon of z Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 1

Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 11 Fltr Dsgn by Optmzton I Wght st Squrs Dsgn : slct on of th bsc forms tht yl lnr phs.g. yp-1 spcfy sr frquncy rspons P,HP,BP, optmzton crtron s whr s wghtng functon cos / / H / H 444 4 3 444 14,...,,...,,..., mn mn F W H H W + + W

Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 1 Fltr Dsgn by Optmzton ths s quvlnt to stnr Qurtc Optmzton problm [ ] [ ]... cos... cos 1... } mn{ 1,..., + μ μ c c W p c c W Q p Q F 4 4 4 8 44 64 7 p Q OP 1

Fltr Dsgn by Optmzton Empl: ow-pss sgn 1. Pssbn Rppl 1 1, <, S P pss - bn stop - bn optmzton functon s.8.6.4. Pssbn Cutoff -> <- Stopbn Cutoff Stopbn Rppl F.. P +.5 1 1.5.5 3,..., 1 + γ.... W... 144 443 pss - bn S 1443 4 stop - bn Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 13

Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 14 Fltr Dsgn by Optmzton smplr problm s obtn by rplcng th F.. by whr th w sr st of n smpl frquncs h qurtc optmzton problm s thn quvlnt to lst-squrs problm +++ : smpl --- : unprctbl bhvor n btwn smpl frquncs. c W W F :,..., { { { } mn{ mn...... + b b b b c c W b S 1 Compr to p.1

Fltr Dsgn by Optmzton thn ll ths s oftn supplmnt wth tonl constrnts Empl: ow-pss P sgn contnu pss-bn rppl control : 1 δ P, < P δ P s pss - bn rppl stop-bn rppl control : δs, S δs s stop - bn rppl Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 15

Fltr Dsgn by Optmzton Empl: ow-pss P sgn contnu rlstc wy to mplmnt ths constrnts, s to mpos th constrnts only on st of smpl frquncs, P1 P,..., Pm n th pss-bn n, S1 S,..., Sn n th stop-bn h rsultng optmzton problm s : mnmz : F,..., [... ] subct to pss-bn constrnts P bp stop-bn constrnts S bs `Qurtc nr Progrmmng problm... 1 Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 16

Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 17 Fltr Dsgn by Optmzton II `Mnm Dsgn : slct on of th bsc forms tht yl lnr phs.g. yp-1 spcfy sr frquncy rspons P,HP,BP, optmzton crtron s whr s wghtng functon ]cos [ / / H / H m mn m mn,...,,..., W H H W W

Fltr Dsgn by Optmzton Concluson: I wght lst squrs sgn II mnm sgn prov gnrl `frmwor, procurs to trnslt fltr sgn problms nto stnr optmzton problms In prctc n n ttboos: mphss on spcfc -hoc procurs : - fltr sgn bs on wnows - qu-rppl sgn Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 18

Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 19 Fltr Dsgn usng Wnows Empl : ow-pss fltr sgn l low-pss fltr s hnc l tm-omn mpuls rspons s trunct h [] to +1 smpls : group ly to turn nto cusl fltr < < < 1 C C H H h c c α sn... 1 ] [. + < < othrws / / ] [ ] [ h h on-cusl n nfntly long

Fltr Dsgn usng Wnows Empl : ow-pss fltr sgn contnu not : t cn b shown tht tm-omn truncton corrspons to solvng wght lst-squrs optmzton problm wth th gvn H, n wghtng functon truncton corrspons to pplyng rctngulr wnow : h[ ] h [ ] w[ ] w[ ] 1 W 1 / < < othrws +++: smpl procur lso for HP,BP, --- : truncton n th tm-omn rsults n Gbbs ffct n th frquncy omn,.. lrg rppl n pss-bn n stop-bn, whch cnnot b ruc by ncrsng th fltr orr. / Fltr Dsgn-FIR cwlu@twns..nctu.u.tw

Fltr Dsgn usng Wnows Rmy : pply wnows othr thn rctngulr wnow: tm-omn multplcton wth wnow functon w[] corrspons to frquncy omn convoluton wth Wz : h[ ] h [ ] w[ ] H z H z W z cnt wnows : Hn, Hmmng, Blcmn, Ksr,. s ttboos wnow choc/sgn tr-off btwn s-lob lvls fn p pss-/stop-bn rppl n wth mn-lob fns trnston bnwth Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 1

Wnowng Effct Gbbs phnomnon Fltr Dsgn-FIR cwlu@twns..nctu.u.tw

Wnowng 3

Equrppl Dsgn Strtng pont s mnm crtron,.g. mn mw mn,...,,..., m E Bs on thory of Chbyshv ppromton n th ltrnton thorm, whch roughly stts tht th optml sr such tht th m mmum wght ppromton rror s obtn t + trml frquncs m E E for 1,.., + tht hnc wll hbt th sm mmum rppl qurppl Itrtv procur for computng trml frquncs, tc. Rmz chng lgorthm, Prs-McCllln lgorthm Vry flbl, tc., vlbl n mny softwr pcgs Dtls omtt hr s ttboos Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 4

Softwr FIR Fltr sgn bunntly vlbl n commrcl softwr Mtlb: bfr1n,wn,typ,wnow, wnow lnr-phs FIR sgn, n s fltr orr, Wn fns bn-gs, typ s `hgh,`stop, bfrn,f,m,wnow, wnow FIR sgn bs on nvrs Fourr trnsform wth frquncy ponts f n corrsponng mgntu rspons m brmzn,f,m, qurppl lnr-phs FIR sgn wth Prs-McCllln Rmz chng lgorthm Fltr Dsgn-FIR cwlu@twns..nctu.u.tw 5