DESIGN OF STABLE TWO-DIMENSIONAL IIR NOTCH FILTER USING ROOT MAP
|
|
- Rosanna Hopkins
- 5 years ago
- Views:
Transcription
1 8th European Sgnal Processng Conference EUSIPCO- Aalborg, Denmar, August -7, DESIGN OF STABLE TWO-DIMENSIONAL IIR NOTC FILTER USING ROOT MAP Soo-Chang Pe and Chen-Cheng Tseng Depart. of Electrcal Engneerng, Natonal Tawan Unversty Tape, Tawan Depart. of Computer and Communcaton Engneerng Natonal Kaohsung Frst Unversty of Sc. and Tech. Kaohsung, Tawan ABSTRACT In ths pape the two-dmensonal -D IIR notch flter desgn problem s presented. Frst, the lnear fractonal transformaton LFT n complex analyss s used to construct the desred root map of -D polynomal. Then, the numerator and denomnator of transfer functon of -D IIR notch flter can be obtaned from the LFT used n the constructon of root map. Ths desgn not only has closed-form transfer functon but also satsfes the bounded nput/bounded output BIBO stablty condton. Fnally, desgn example and mage applcaton are demonstrated to show the effectveness of the proposed approach.. INTRODUCTION In many sgnal processng applcatons, there s a need for a notch flter whch s charactered by a unt gan at all frequences except at the snusodal frequences where ther gan s ero. These applcatons nclude communcaton, control, mage processng and bomedcal engneerng etc. In one-dmensonal -D case, a typcal example s to cancel 5 or 6 power lne nterference n the recordng of electrocardograms ECGS []. In two-dmensonal -D case, two examples are to elmnate a -D snusodal nterference pattern supermposed on an mage [] and to reduce blocng artfact from DCT coded mage []. Thus, t s nterestng to desgn notch flter to remove the snusodal nterferences corrupted on a desred sgnal. The deal frequency response of -D notch flter s gven by, ± N N D, otherwse where N N s the prescrbed notch frequency. The problem s how to desgn -D flter, to approxmate D as well as possble. So fa the desgn of - D notch flter can be classfed nto two categores. One s FIR flter desgn, the other s IIR flter desgn. In the FIR case, the desgn methods nclude least squares approach [4], transformaton technques [5], sngular value decomposton approach [6] and constraned nonlnear optmaton method [7]. In the IIR case, the desgn approaches contan decomposton method [8] and outer product expanson [9]. So fa the -D adaptve IIR notch flters based on outer product expanson have been developed for removng snusodal nterference wth unnown or tme-varyng frequences [][]. On the other hand, the root map of -D polynomal has been a useful tool for checng the stablty of -D IIR dgtal flters []-[4]. Due to the success of ths tool n stablty checng, we wll use the root map to desgn -D notch flter n ths paper. The proposed desgn not only has closed-form transfer functon but also satsfes the BIBO stablty condton. Compared wth conventonal methods n [8][9], our method does not need to decompose the -D notch flter desgn nto several sub-flter desgns, so the proposed approach s easer to use than conventonal methods. Now, the desgn detals are descrbed n next sectons.. ROOT MAP AND LINEAR FRACTIONAL TRANSFORM In ths secton, the root map of -D polynomal s frst defned. Then, we study the relaton between root map and stablty of -D IIR flter. Fnally, the lnear fractonal transformaton s descrbed.. Defnton of Root Map Consder the -D polynomal wth order N, below: N N P, p,, N If P s regarded as a mappng between the plane and plane, the root map s defned as the mage of the unt crcle of one varable n the plane of other varable under ths mplct mappng P,. Now, let us descrbe the detals of these two maps below: Frst, the -D polynomal P, can be rearranged as the followng form: N N P, p, P[ ] [ ] The above expresson P can be nterpreted as -D polynomal n varable whose coeffcents are polynomal n varable. When s fxed, the P [ ] can be factored to yeld N roots. Thus, the frst root map Γ s defned as the loc of N roots of P [ ] when parameter traverses the unt crcle e. Second, the -D EURASIP, ISSN
2 polynomal P, can be rearranged as the followng form: N N P, p, 4 P[ ] [ ] where the expresson P can be nterpreted as -D polynomal n varable whose coeffcents are polynomal n varable. When s fxed, the P [ ] can be factored to yeld N roots. Thus, the second root map Γ s defned as the loc of N roots of P [ ] when parameter traverses the unt crcle e.. Stablty and Root Map For the -D frst quadrant IIR flter wth order M, M, N,, ts transfer functon s gven by N, M M N N q p,, Q, P, In [], Shans has shown that, s stable f and only f the followng condtons are true: P, f, 6a P, f, 6b It s obvous that the condton n Eq.6a mples that the root map Γ of the unt crcle of varable les nsde the unt crcle n the plane. And, the condton n Eq.6b means that the root map Γ of the unt crcle of les wthn the unt crcle n the plane. Thus, stablty condton of Shan's theorem wll be satsfed f two root maps le wthn unt crcle. So, we can control the stablty of -D IIR flter by controllng two root maps to le nsde the unt crcle.. Lnear Fractonal Transform In [5], the lnear fractonal transformaton LFT from complex varable to s defned by c + c 7 c + c4 where c, c, c and c are constants. It can be shown that 4 the LFT transforms crcles and lnes nto crcles and lnes. Moreove the Eq.7 can be rewrtten n the form c + c4 c c 8 Ths expresson s the same as P, wth order, f we choose p, c, p, c, 4 p, c, and p, c. Thus, we can use the propertes of LFT n complex analyss to construct the desred root maps of P, wth order,. In next 5 secton, the detals that desgn -D IIR notch flter usng root map and LFT wll be descrbed.. DESIGN OF -D IIR NOTC FILTER In ths secton, we wll summare some mportant root maps. Then, these root maps are used to desgn -D IIR notch flter. The frst mportant root map s descrbed by the followng fact: Fact : Choosng parameters ρ r r a, ρ b 9 r + r + and let -D polynomal wth order, be B, a + b + b + then we have the results: The two root maps Γ and Γ are the same. The two root maps are both crcles wth center at ρ and radus r. The equaltes B, ρ + r and B ρ + hold. Proof: Because B, s a specal LFT n complex varable textboo n [5], ths fact can be proved by usng technques n [5]. The detals are descrbed below. Snce the symmetry property B, B, s vald, two root maps are the same. So, we only need to fnd the root map Γ because Γ s the same as Γ. To fnd Γ, let complex varables and be represented by x + jy u + jv Substtutng Eq. nto B,, t yelds b + u x vy bu a vx + b + u y bv Snce s on the unt crcle, we have the constrant x + y. Based on ths constrant and Eq., we get b a a b u + v b b Substtutng Eq.9 nto Eq., we have u ρ + v r 4 Ths means that the root map Γ s a crcle wth center at ρ and radus r. Fnally, t s easy to show that the equaltes B, ρ + r and B ρ + hold by usng drect substtuton. Based on the results n Fact, we have the further results of root map below: Fact : Let -D polynomal wth order, be F Φ,, ρ, r,, [,, ] 5 where matrx Φ s 686
3 b cos b a cos + 6 Φ b cos ab cos + b cos b ab cos a wth parameters ρ r r a and ρ b, then we r + r + have the results: The root map Γ n plane conssts of two crcles wth center at ρe ± and rad r. The root map Γ n plane conssts of two crcles wth center at ρe ± and rad r. The equaltes F e, ρ + r e, ρ,, and e, ρ + r e, ρ,, hold. F ρ + r e, e, ρ,, ρ + r e, e, ρ,, hold. F 4 The equaltes and F Proof: The -D polynomal n Eq.5 can be factored as F,, ρ,, B, B, 7 where B, a e e + b e 8a + b e + B, a e e + b e 8b + b e + Thus, F,, ρ,, mples that B, or B,. Usng ths factoraton and the results n Fact, the proof can be shown easly. Now, let us use the Fact to construct the transfer functon of -D IIR notch flter wth order,,,. The desgn problem s how determne the flter coeffcents p, and q, n Eq.5 such that the frequency response e, e N N approxmates the deal response D n Eq. as well as possble f the notch frequency s gven. In ths pape the ey step to desgn IIR notch flter, s to choose numerator polynomal Q, and denomnator polynomal P, n Eq.5 as the form: Q, F,, ρ, r N N 9a P, F,, ρ, r N N 9b where parameters ρ, r, ρ, r need to be prescrbed. In ths pape we choose these four parameters below: ρ α, r α a ρ α, r α ε b where α s a number n the nterval,, and ε s a very small postve number. Based on ths choce, the root maps of numerator polynomal Q, and denomnator polynomal P, are shown n Fg.. Now, several remars are made below: From Eq.a, we have ρ r. Ths equalty maes the two root map crcles of numerator polynomal Q, be tangent to the unt crcle. Usng Fact and ρ r, t can be shown that N N N N Q e, e Q e, e Thus, there s ero gan at the notch frequency for the desgned IIR notch flter. The reason to choose ε as a postve number s to mae the root map crcles of denomnator polynomal P, le wthn unt crcle. Thus, the desgned IIR notch flter s guaranteed stable. In order to obtan unt gan at all frequences except at the notch frequency, the ε must be chosen as small as possble. When ε s very small, the root maps of P, wll be very close to ones of Q,. 4 Because we choose ρ ρ, the root map crcles of P, and Q, have the same center. Thus, we may refer to ths desgn as "concentrc crcle desgn". Fnally, substtutng Eq.9 and Eq. nto Eq.5, the transfer functon of the desgned -D IIR notch flter s gven by Q,, P, F,, α, α N N F,, α, α ε N N Usng Eq.5, the above transfer functon can be computed easly wthout requrng optmaton and complcated calculaton. In next secton, numercal examples wll be used to llustrate the effectveness of the proposed desgn. 4. DESIGN EXAMPLE AND APPLICATION In ths secton, one desgn example of proposed -D IIR notch flter s frst llustrated. Then, we compare the proposed desgn method wth conventonal outer product expanson method n [9]. Fnally, we use the -D IIR notch flter to remove snusodal nterference supermposed on an mage. Example : -D IIR Notch Flter Desgn In ths example, we wll study the performance of the desgned -D IIR notch flter. The notch frequency s chosen as N N.5π. And, the desgn parameters are α. and ε.. Fg.a show the resultant magntude response. It s clear that j e, e has unt gan at all frequences except at notch frequency.5π where gan s ero. So, the specfcaton s ftted well. oweve the detals of the notch are underneath the unt gan plane. In order to show the performance of the desgned flter bette Fg.b plots the loss j e, e. From ths result, we observe that the magntude response has small rpples n the vcnty of the notch. Fortunately, ths rpple can be reduced by reducng the parameter ε. To llustrate ths fact, Fg.cd show the magntude response and loss of notch flter when ε s reduced to. and other 687
4 desgn parameters are not changed. From ths result, t s clear that error rpple near notch frequency has been reduced. Moreove t s nterestng to study the shape of -db contour of notch n whch the magntude j e, e s equal to. Fg. shows the -db contours of the desgned -D IIR notch flters, for varous parameters α and ε when notch frequency s chosen as N N.5π. It s clear that the -db contour s symmetrc about the notch frequency.5π. And, the se of -db contour s reduced f parameters α and ε are reduced. Thus, we can use parameters α and ε to control the shape of -db contour of notch. Example : Comparson wth Conventonal -D IIR Notch Flter Desgn In ths example, we wll compare the proposed desgn method wth conventonal outer product expanson method n [9]. Gven the notch frequency and parameter N N BW, the transfer functon of conventonal -D IIR notch flter n [9] s gven by, b b a a where flters a a + b, 4 a + a b + a, 5 + b wth the coeffcents BW cos N tan a, a BW, 6 BW + tan + tan N π 4 N π + 4 sn b, 7 sn Fg.4 shows the magntude response j e, e and the loss j e, e of ths conventonal IIR notch flter for N N.5π and BW.5π. Comparng Fg. wth Fg.4, t can be seen that the performance of the proposed desgn method s smlar to one of the conventonal method n [9]. oweve the order of the proposed flter, n Eq. s,,, and the order of conventonal flter, n Eq. s,,,. So, the mplementaton complexty of proposed flter wll be less than the complexty of conventonal flter f the drect-form realaton s chosen. Example : Snusodal Interference Removal on Image In ths example, we wll use the proposed IIR notch flter to remove the snusodal nterference supermposed on an mage wth se 5 5. The mage shown n Fg.5a s the Lae mage corrupted by a snusodal nterference below: 5 sn.π m +.πn 8 Now, let us desgn a -D IIR notch flter wth notch frequency N N.π,.π, α. and ε. to remove ths snusodal nterference n spatal doman. The fltered mage s shown n Fg.5b. It s clear that the nterference has been removed by the proposed -D IIR notch flter. In the above, only the case of sngle snusod s studed. For the case of multple snusodal nterferences, the cascade notch flter can be used to remove the multple nterferences f the notch frequences are chosen as the frequences of snusods. 5. CONCLUSIONS In ths pape the root map has been presented to desgn a stable -D IIR notch flter. Frst, the lnear fractonal transformaton LFT n complex analyss s used to construct the desred root map of -D polynomal. Then, the numerator and denomnator of transfer functon of -D IIR notch flter can be obtaned from the LFT used n the constructon of root map. Fnally, desgn example and mage applcaton are demonstrated to show the effectveness of the proposed approach. oweve only IIR notch flter desgn s consdered here. Thus, t s nterestng to desgn other -D IIR flters by usng root map n the future. REFERENCES [] S.C. Pe and C.C. Tseng, "Elmnaton of AC nterference n electrocardogram usng IIR notch flter wth transent suppresson" IEEE Trans. on Bomedcal Engneerng, pp.8-, Nov [] R.C. Gonale and R.E. woods, Dgtal Image Processng, nd Edton, Prentce-all,. [] V.K. Srvastava and G.C. Ray, "Desgn of D-multple notch flter and ts applcaton n reducng blocng artfact from DCT coded mage," Proc. of the nd Annual EMBS Internatonal Conference, pp.89-8, July. [4] S.C. Pe and C.C. Tseng, "Two-dmensonal IIR and FIR dgtal notch flter desgn" Proc. of the ISCAS 99, pp.89-89, May 99. [5] V.L. Narayana Murthy and A. Mau "Desgn of some -D flter through the transformaton technque," IEE Proc. Vson, Image and Sgnal Processng, vol.4, pp.84-9, June 996. [6] S.C. Pe, W.S. Lu and C.C. Tseng, "Two-dmensonal FIR notch flter desgn usng sngular value decomposton," IEEE Trans. on Crcuts and Systems-I, pp.9-94, Mar [7] F. Wysoca-Schlla, "Desgn of equrpple -D lnear-phase FIR notch flters," EUROCON 7 The Internatonal Conference on "Computer as a Tool", pp.-5, Sept. 7. [8] S.C. Pe and C.C. Tseng, "Two dmensonal IIR dgtal notch flter desgn," IEEE Trans. on Crcuts and Systems-II, pp.7-, Mar [9] S.C. Pe, W.S. Lu and C.C. Tseng, "Analytcal twodmensonal IIR notch flter desgn usng outer product expanson," IEEE Trans. on Crcuts and Systems-II, pp , Sept [] T. namoto, N. Ieda, S. Nshmura and A. Do, "Desgn of two-dmensonal adaptve notch flters," Proc. of 5th Internatonal Conference on Sgnal Processng, pp.58-54, Aug.. [] S.C. Pe, C.L. Wu and J.J. Dng, "Smplfed structures for twodmensonal adaptve notch flters" Proc. of the ISCAS, pp.iv-46-iv-49, May. [] J. L. Shans, S. Tretel and J.. Justce, "Stablty and synthess of two-dmensonal recursve flters," IEEE Trans. on Audo Electroacoust., pp.5-8, June
5 [] D.E. Dudgeon and R.M. Mersereau, Multdmensonal Dgtal Sgnal Processng, Prentce-all, 984. [4] J.W. Woods, Multdmensonal Sgnal, Image, and Vdeo Processng and Codng, Academc Press, 6. [5] J.W. Brown and R.V. Churchll, Complex Varables and Applcatons, 7th Edton, McGraw-ll, 4. frequency x π.6 α. and ε. frequency x π.6 α. and ε frequency x π..6 frequency x π α. and ε.5 α. and ε. a b frequency x π.6. frequency x π frequency x π..6 frequency x π Fg. The -db contours of the desgned -D IIR notch flters, for varous parameters α and ε when notch frequency s chosen as.5π,.5 N N π. c d Fg. The root maps of numerator polynomal Q, and denomnator polynomal P, of -D IIR notch flter a Root map Γ of Q,. b Root map Γ of Q,. c Root map Γ of P,. d Root map Γ of P,. a b Fg.4 The magntude response e, e and loss j e, e of the conventonal -D IIR notch flter for N N.5π and BW.5 π. a b 4 4 c d Fg. The magntude response e, e and loss j e, e of the desgned IIR notch flter ab The results wth α. and ε.. cd The results wth α. and ε a b Fg.5 Example of snusodal nterference removal on an mage a The corrupted mage b The mage restored by usng -D IIR notch flter. 689
A New Design Approach for Recursive Diamond-Shaped Filters
A ew Desgn Approach for Recursve Damond-Shaped Flters RADU MATEI Faculty of Electroncs, Telecommuncatons and Informaton Technology Techncal Unversty Gh.Asach Bd. Carol I no., Ias 756 ROMAIA rmate@etc.tuas.ro
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More information290 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: FUNDAMENTAL THEORY AND APPLICATIONS, VOL. 45, NO. 3, MARCH H d (e j! ;e j!
9 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: FUNDAMENTAL THEORY AND APPLICATIONS, VOL. 45, NO. 3, MARCH 998 Transactons Brefs Two-Dmensonal FIR Notch Flter Desgn Usng Sngular Value Decomoston S.-C. Pe,
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationCOEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN
Int. J. Chem. Sc.: (4), 04, 645654 ISSN 097768X www.sadgurupublcatons.com COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN R. GOVINDARASU a, R. PARTHIBAN a and P. K. BHABA b* a Department
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationYong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )
Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationEE 508 Lecture 7. Degrees of Freedom The Approximation Problem
EE 508 Lecture 7 Degrees of Freedom The Approxmaton Problem vew from Last Tme Desgn Strategy Theorem: A crcut wth transfer functon T(s) can be obtaned from a crcut wth normalzed transfer functon T n (s
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationarxiv:cs.cv/ Jun 2000
Correlaton over Decomposed Sgnals: A Non-Lnear Approach to Fast and Effectve Sequences Comparson Lucano da Fontoura Costa arxv:cs.cv/0006040 28 Jun 2000 Cybernetc Vson Research Group IFSC Unversty of São
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationPop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing
Advanced Scence and Technology Letters, pp.164-168 http://dx.do.org/10.14257/astl.2013 Pop-Clc Nose Detecton Usng Inter-Frame Correlaton for Improved Portable Audtory Sensng Dong Yun Lee, Kwang Myung Jeon,
More informationFundamental loop-current method using virtual voltage sources technique for special cases
Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More information8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS
SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationAppendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis
A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationLecture 16 Statistical Analysis in Biomaterials Research (Part II)
3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan
More informationLectures on Multivariable Feedback Control
Lectures on Multvarable Feedback Control Al Karmpour Department of Electrcal Engneerng, Faculty of Engneerng, Ferdows Unversty of Mashhad June 200) Chapter 9: Quanttatve feedback theory Lecture Notes of
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationarxiv: v1 [math.co] 12 Sep 2014
arxv:1409.3707v1 [math.co] 12 Sep 2014 On the bnomal sums of Horadam sequence Nazmye Ylmaz and Necat Taskara Department of Mathematcs, Scence Faculty, Selcuk Unversty, 42075, Campus, Konya, Turkey March
More informationResearch Article Green s Theorem for Sign Data
Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationThe Minimum Universal Cost Flow in an Infeasible Flow Network
Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran
More informationRegularized Discriminant Analysis for Face Recognition
1 Regularzed Dscrmnant Analyss for Face Recognton Itz Pma, Mayer Aladem Department of Electrcal and Computer Engneerng, Ben-Guron Unversty of the Negev P.O.Box 653, Beer-Sheva, 845, Israel. Abstract Ths
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationIdentification of Linear Partial Difference Equations with Constant Coefficients
J. Basc. Appl. Sc. Res., 3(1)6-66, 213 213, TextRoad Publcaton ISSN 29-434 Journal of Basc and Appled Scentfc Research www.textroad.com Identfcaton of Lnear Partal Dfference Equatons wth Constant Coeffcents
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationSuppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl
RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com
More informationA Fast Computer Aided Design Method for Filters
2017 Asa-Pacfc Engneerng and Technology Conference (APETC 2017) ISBN: 978-1-60595-443-1 A Fast Computer Aded Desgn Method for Flters Gang L ABSTRACT *Ths paper presents a fast computer aded desgn method
More informationSINGLE OUTPUT DEPENDENT QUADRATIC OBSERVABILITY NORMAL FORM
SINGLE OUTPUT DEPENDENT QUADRATIC OBSERVABILITY NORMAL FORM G Zheng D Boutat JP Barbot INRIA Rhône-Alpes, Inovallée, 655 avenue de l Europe, Montbonnot Sant Martn, 38334 St Ismer Cedex, France LVR/ENSI,
More informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationPerron Vectors of an Irreducible Nonnegative Interval Matrix
Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of
More informationA new construction of 3-separable matrices via an improved decoding of Macula s construction
Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula
More informationCHAPTER 4 SPEECH ENHANCEMENT USING MULTI-BAND WIENER FILTER. In real environmental conditions the speech signal may be
55 CHAPTER 4 SPEECH ENHANCEMENT USING MULTI-BAND WIENER FILTER 4.1 Introducton In real envronmental condtons the speech sgnal may be supermposed by the envronmental nterference. In general, the spectrum
More informationExplicit constructions of all separable two-qubits density matrices and related problems for three-qubits systems
Explct constructons of all separable two-qubts densty matrces and related problems for three-qubts systems Y. en-ryeh and. Mann Physcs Department, Technon-Israel Insttute of Technology, Hafa 2000, Israel
More informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationChapter 7 Channel Capacity and Coding
Chapter 7 Channel Capacty and Codng Contents 7. Channel models and channel capacty 7.. Channel models Bnary symmetrc channel Dscrete memoryless channels Dscrete-nput, contnuous-output channel Waveform
More informationLOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin
Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationProjective change between two Special (α, β)- Finsler Metrics
Internatonal Journal of Trend n Research and Development, Volume 2(6), ISSN 2394-9333 www.jtrd.com Projectve change between two Specal (, β)- Fnsler Metrcs Gayathr.K 1 and Narasmhamurthy.S.K 2 1 Assstant
More informationSupplemental document
Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationEnergy Storage Elements: Capacitors and Inductors
CHAPTER 6 Energy Storage Elements: Capactors and Inductors To ths pont n our study of electronc crcuts, tme has not been mportant. The analyss and desgns we hae performed so far hae been statc, and all
More informationLecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.
prnceton u. sp 02 cos 598B: algorthms and complexty Lecture 20: Lft and Project, SDP Dualty Lecturer: Sanjeev Arora Scrbe:Yury Makarychev Today we wll study the Lft and Project method. Then we wll prove
More informationTHE HARTLEY TRANSFORM IN A FINITE FIELD
THE HARTLEY TRANSFORM IN A FINITE FIELD R. M. Campello de Souza H. M. de Olvera A. N. Kauffman CODEC - Grupo de Pesusas em Comuncações Departamento de Eletrônca e Sstemas - CTG - UFPE C.P. 78 57-97 Recfe
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationChapter 7 Channel Capacity and Coding
Wreless Informaton Transmsson System Lab. Chapter 7 Channel Capacty and Codng Insttute of Communcatons Engneerng atonal Sun Yat-sen Unversty Contents 7. Channel models and channel capacty 7.. Channel models
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationLimited Dependent Variables
Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages
More informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationThe Geometry of Logit and Probit
The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.
More informationThis model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:
1 Problem set #1 1.1. A one-band model on a square lattce Fg. 1 Consder a square lattce wth only nearest-neghbor hoppngs (as shown n the fgure above): H t, j a a j (1.1) where,j stands for nearest neghbors
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More information829. An adaptive method for inertia force identification in cantilever under moving mass
89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,
More informationAerodynamics. Finite Wings Lifting line theory Glauert s method
α ( y) l Γ( y) r ( y) V c( y) β b 4 V Glauert s method b ( y) + r dy dγ y y dy Soluton procedure that transforms the lftng lne ntegro-dfferental equaton nto a system of algebrac equatons - Restrcted to
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationTime-Varying Systems and Computations Lecture 6
Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy
More informationELE B7 Power Systems Engineering. Power Flow- Introduction
ELE B7 Power Systems Engneerng Power Flow- Introducton Introducton to Load Flow Analyss The power flow s the backbone of the power system operaton, analyss and desgn. It s necessary for plannng, operaton,
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More information), it produces a response (output function g (x)
Lnear Systems Revew Notes adapted from notes by Mchael Braun Typcally n electrcal engneerng, one s concerned wth functons of tme, such as a voltage waveform System descrpton s therefore defned n the domans
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationBOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS
BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all
More informationChapter 11: Simple Linear Regression and Correlation
Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationAdaptive sliding mode reliable excitation control design for power systems
Acta Technca 6, No. 3B/17, 593 6 c 17 Insttute of Thermomechancs CAS, v.v.. Adaptve sldng mode relable exctaton control desgn for power systems Xuetng Lu 1, 3, Yanchao Yan Abstract. In ths paper, the problem
More informationGEMINI GEneric Multimedia INdexIng
GEMINI GEnerc Multmeda INdexIng Last lecture, LSH http://www.mt.edu/~andon/lsh/ Is there another possble soluton? Do we need to perform ANN? 1 GEnerc Multmeda INdexIng dstance measure Sub-pattern Match
More informationADVANCED MACHINE LEARNING ADVANCED MACHINE LEARNING
1 ADVANCED ACHINE LEARNING ADVANCED ACHINE LEARNING Non-lnear regresson technques 2 ADVANCED ACHINE LEARNING Regresson: Prncple N ap N-dm. nput x to a contnuous output y. Learn a functon of the type: N
More informationPolynomials. 1 More properties of polynomials
Polynomals 1 More propertes of polynomals Recall that, for R a commutatve rng wth unty (as wth all rngs n ths course unless otherwse noted), we defne R[x] to be the set of expressons n =0 a x, where a
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More information8.6 The Complex Number System
8.6 The Complex Number System Earler n the chapter, we mentoned that we cannot have a negatve under a square root, snce the square of any postve or negatve number s always postve. In ths secton we want
More informationGraph Reconstruction by Permutations
Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer
More informationPerfect Competition and the Nash Bargaining Solution
Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange
More information