Design of Infinite Impulse Response (IIR) filters with almost linear phase characteristics via Hankel-norm approximation
|
|
- Christal Henry
- 5 years ago
- Views:
Transcription
1 Design of Infinite Impulse Response IIR filters with almost linear phase characteristics via Hankel-norm approximation GD Halikias, IM Jaimoukha and E Milonidis Abstract We propose a method for designing Infinite Impulse Response IIR filters with specified magnitude frequency-response tolerances and approximately linear phase characteristics The design method consists of two steps First, a Finite Impulse Response FIR filter with linear phase is designed using standard optimisation techniques eg linear programming In the second stage, the designed FIR filter is approximated by a stable IIR filter using a Hankel-norm approximation method This is based on a recent technique for approximating discrete-time descriptor systems and requires only standard linear algebraic routines, while avoiding altogether the solution of two matrix Lyapunov equations which is computationally expensive A detailed solution is developed in state-space form, using only the Markov parameters of the FIR filter Further, a-priori bounds on the magnitude and phase errors are obtained which may be used to select the lowest-order IIR filter order which satisfies the specified design tolerances The efectiveness of the method is illustrated with a numerical example Keywords: FIR filters, linear phase, Hankel-norm approximation, descriptor systems 2 Introduction Finite Impulse Response FIR filters are widely used in many digital signal processing applications, especially in the field of communications, because they can exhibit linear phase characteristics Unfortunately, in many cases the order of such filters is prohibitively high for practical implementation In general, the number of delay elements and multipliers needed for an FIR design tends to be much higher compared to similar Infinite Impulse Response IIR implementations This is especially true for filters with sharp cut-off characteristics which have a long impulse response [6, 9] It is therefore natural to ask whether a linear-phase FIR filter Hz can be approximated by a loworder IIR filter, without degrading significantly its magnitude and phase characteristics In this paper we apply Hankel-norm approximation techniques to matrix FIR filters Our technique is based on a recent result involving approximations of discrete-time descriptor systems [, 2] This allows us to treat systems with poles at the origin and has been applied successfully in the context of mixed H /H 2 optimisation problems [5] Our results apply both to the γ-suboptimal and the strictly optimal Control Engineering Research Centre, School of Engineering and Mathematical Sciences, City University, Northampton Square, London ECV HB, UK Department of Electrical and Electronic Engineering, Imperial College of Science, Technology and Medicine, Exhibition Road, London SW7 2BT
2 problem, although in the later case the resulting state-space realisation is non-minimal Using an allpass matrix dilation technique, a state-space parametrisation of the complete family of solutions can be obtained, in the form of a linear fractional map of the ball of unstable contractions [, 2, 4] This can be subsequently reformulated using only the matrix Markov parameters of Hz When the results are specialised to the scalar case, magnitude and phase error bounds of the approximation error can be obtained in terms of the Hankel singular values of Hz, which allows the a-priori determination of the IIR filter order satisfying specified magnitude ripple specifications and an acceptable phase deviation from linearity The effectiveness of the method is illustrated via a numerical example, involving the approximation of a high-order linear-phase FIR filter, designed using linear programming techniques [8] Finally, extensions and further applications of the method are briefly discussed 3 Notation Most of the notation used is standard and is reproduced here for completeness D denotes the open unit disc D = {ζ C : ζ < }, with D and D its closure and boundary respectively L 2 D denotes the Hilbert space of all matrix-valued functions F defined on the unit circle such that 2π trace[f e jθ F e jθ ]dθ < where denotes the complex conjugate transpose of a matrix The corresponding inner product of two L 2 D functions F and G of compatible dimensions is given as: < F, G >= 2π trace[f e jθ Ge jθ ]dθ 2π H 2 D and H 2 D are the closed subspaces of L 2 consisting of all functions analytic in C \ D and D, respectively L D is the space of all uniformly-bounded matrix-valued functions in D, ie all functions defined on the unit circle whose norm: F = sup σ[f e jθ ] θ [,2π is finite Here σ denotes the largest singular value of a matrix H D and H D denote the closed subspaces of L D consisting of all functions analytic in C \ D and D, respectively, while H,k D is the set of all functions in L D with no more than k poles in D Spaces of real-rational functions will be indicated by the suffix R before the corresponding space symbol The unit ball of H D is the set BH D = {U H D : U } If G L D, then the Hankel operator with symbol G is defined as: Γ G : L 2 D H 2 D, Γ G = Π + M G H 2 D in which M G denotes the multiplication operator M G : L 2 D L 2 D, M G f = Gf and Π + denotes the orthogonal projection Π + : L 2 D H 2 D If H L D has dimensions p +p 2 m +m 2 and U RL D has dimensions m 2 p 2, we define the lower linear fractional map of H and U as F l H, U = H + H 2 UI H 22 U H 2, provided the indicated inverse exists If U is a set of m 2 p 2 matrix functions, then F l H, U denotes the set {F l H, U : U U} 2
3 4 FIR filters with linear phase response The fact that FIR filters can have a linear phase response has been used extensively in many digital signal processing applications A filter with a nonlinear phase characteristic causes distortion to its input signal, since the various frequency components in the signal will be delayed, in general, by amounts which are not proportional to their frequency, thus altering their mutual harmonic relationships A distortion of this form is undesirable in many practical digital signal processing applications such as music, video, data transmission and biomedicine and can be avoided by using filters with linear phase over the frequency range of interest [9] Suppose that the unit pulse response of an FIR filter is given by {h, h, h2,, hn} For this filter to have a linear phase response, say θω = αω for α constant, it must have an impulse response with positive symmetry, ie hn = hn n, where n varies as n =,, 2,, N 2 for N odd and n =,, 2,, N 2 for N even Such a filter has identical phase and group delay and these are independent of frequency, ie T p = θω ω = T g = dθω dω An FIR filter whose impulse response has negative symmetry, ie hn = hn n, has an affine phase-response characteristic of the form θω = β αω with α and β constant and its group delay T g = α is independent of ω The simple relationship between the impulse response of an FIR filter and its frequency response has been exploited to design filters of this type via a variety of optimisation methods [3], [7], [8], [] One of the earliest and simplest approaches is [8], in which a linear programming procedure is proposed for designing linear phase FIR filters with specified magnitude-response characteristics low-pass, high-pass, etc As an example, consider the filter Hz = h + hz + + hp z p + hpz p + hp + z p+ + + h2pz 2p in which positive symmetry of the impulse response is enforced around the central coefficient hp by setting h = h2p, h = h2p,, hp = hp + The frequency response of the filter may be written as: He jω = e jωp Mω = α where Mω is a real linear function of the filter s coefficients: Mω = 2 cospω 2 cosp ω h h hp Suppose now that we want to satisfy the following specifications: a δ He jω + δ for all frequencies in the pass-band ω [, ω p ], b He jω δ for all frequencies in the stopband ω [ω s, π] where ω p < ω s, and c minimise the ripple δ subject to constraints a and b Specifications a and b can enforced by discretising the pass-band and stop-band frequency intervals using n frequencies, say, ie = ω < ω 2 < < ω n = ω p 3
4 and ω s = ω n+ < ω n+2 < < ω 2n = π and enforcing the specifications via the inequalities: δ Mω i + δ, i =, 2,, n and δ Mω i δ, i = n +, n + 2,, 2n 2 The optimisation problem: min δ subject to and 2, can now be formulated as a linear programme in the standard form: min c x st Ax b where x = h h hp δ, c =, b is a 4n-dimensional vector and A is a 4n p + 2- dimensional matrix This can be solved efficiently using standard techniques, eg the simplex algorithm or interior point methods 5 Hankel-norm approximation of FIR filters In this section we propose a Hankel-norm method for approximating FIR filters via lower-order IIR filters The method is based on a recent result involving model-reduction of descriptor discrete-time systems [], [2] The main advantage of Hankel-norm approximation methods over other model-reduction techniques is that they offer tight bounds on the infinity-norm of the approximation error; in particular, it is shown in [4] how to construct k-th order approximations Xz RH D of Hz RH D with degxz = k < deghz = n, such that: n Hz + Xz σ i Γ H 3 i=k+ where σ i Γ H denotes the i-th singular value of Γ H, indexed in non-increasing order of magnitude This inequality can be used to determine a-priori the order k of the low-order system Xz which satisfies magnitude error specifications Theorem 2 below parametrises all solutions Xz to the Hankel-norm approximation problem Γ H + Γ X γ, in which Hz is the matrix FIR filter Hz = H + H z + + H n z n and Xz is a matrix IIR filter of degree degxz k The parametrisation is given in descriptor form and hence applies both in the sub-optimal case σ k+ Γ H < γ < σ k Γ H and the optimal case γ = σ k+ Γ H Before stating this theorem, however, the following preliminary result is needed: Theorem : Let Hz = Hz H = H z + H 2 z H n z n with H i R p l for i =, 2,, n Then: 4
5 The singular values of Γ H, σ i Γ H indexed in decreasing order of magnitude are the singular values of the Hankel matrix: H H 2 H n H n H 2 H 3 H n H 3 H 4 R = H n H n H n 2 A state-space realisation of Hz is given by Hz = CzI A B with: I p A = I p R np np I p B = H n H n H 2 H R np l C = I p R p np Further, the realisation is output balanced, ie the unique solution of the Lyapunov equation Q = A QA + C C is given by Q = I np 3 The controllability grammian of the realisation in part 2, ie the unique solution of the Lyapunov equation P = AP A + BB, can be factored as P = R 2 R 2 where: H n H R 2 = n H n Equivalently P can be block-partitioned as: and σ 2 i Γ H = λ i P for each i H H 2 H n mini,j P ij j=,2,n i=,2,n where P ij = H n i+k H n j+k k= Proof: Straightforward and hence omitted The following Theorem gives a complete parametrization of all k th order γ-suboptimal Hankelnorm approximations in descriptor form All optimal approximations can also be obtained by setting γ = σ k+ Γ H Theorem 2: Let Hz be the matrix FIR filter defined in Theorem Let γ satisfy σ k+ Γ H γ < σ k Γ H Then all Xz H,k D such that Hz + Xz γ are generated via the lower linear fractional transformation: X = {F l S a, U : γu BH D} 5
6 where: with and E S = B S = S a = D S + C S ze S A S B S I np I np γ 2 I np P, A S = X A CX P, C S = B B C CX B γi p D s = γi l where A, B, C and P are defined in Theorem, and X = I np +A, A = P γ 2 XA X, B = γxx C and C = γb X Proof: Minor adaptation of result in [], [2] The proof of Theorem 2 follows readily by specialising the results of [] and [2] to our case, using the state-space description in Theorem Note that the generator of all Hankel norm approximations S a z is given in descriptor form; this can be converted into a state-space form, if required, using a standard procedure In particular, it is always possible to obtain a state-space realisation of S a z of order 2np RankA ; if A is non-singular, we obtain a state-space description of order np Note also, that although Theorem 2 is still valid in the optimal case γ = σ k+ Γ H, the resulting realisation is non-minimal In this case, a minimal realisation can be obtained in closed form via a singular perturbation argument see [2] for details It is clear from Theorem 2 that the solution of Lyapunov equations is completely avoided in our framework; in fact the generator S a z is completely defined by the Markov parameters of Hz and γ We make this dependence explicit in Theorem 4 below, which derives a state-space realization of the generator of all γ-suboptimal Hankel-norm approximations We first reformulate the parametrization of all Hankel-norm approximations in a state-space rather than descriptor setting: Theorem 3: Let Hz be the matrix FIR filter defined in Theorem Let γ satisfy σ k+ Γ H γ < σ k Γ H and assume that A is non-singular Then all Xz H,k D such that Hz+Xz γ are generated via the lower linear fractional transformation: where: with X = {F l H a, U : γu BH D} H a = D H + C H zi A H B H A H = I pn γ 2 I pn P A X, B H = γ 2 I n P A B B C H = CX P C A X, D H = CX B CX P A B γi CX P A B γi C A B C A B 6
7 where all variables are defined in Theorem and Theorem 2 Proof: Follows by the standard procedure of transforming descriptor realisations to state-space form The following Theorem shows that in the state-space realization given in Theorem 3 has a special structure In particular, the corresponding A matrix of the realization is in companion form, which allows us to obtain an analytic expression for the transfer function of the so-called central approximation corresponding to Uz = For simplicity the results are stated only for the scalar case Theorem 4: Let all variables be defined as in Theorems -3 above, set p = m = and assume that A is non-singular and that σ k+ Γ h < γ < σ k+ Γ h sub-optimal case Then the realization of H a z = A H, B H, C H, D H defined in Theorem 3 is of the form: a a 2 a 3 a n a n b b 2 h n A H =, B H = h n 2 h C H = c c 2 c 3 c n c n, D H = d 22 where hz = h z + h 2 z h n z n and in which a i = γ2 θ i + θ j γ 2 i =, 2,, n, a n = γ2 θ n θ γ 2 θ γ 2 n γ b = h n + γ 2 θ i h n i+, b 2 = θ γ 2 θ θ j = i= n i ˆPij, j =, 2,, n, i= ˆP = P γ 2 I n In particular, the transfer function of the central approximation Uz = is given as: Xz = h z n + h 2 h a z n h n h a n 2 z + b h a n z n a z n a 2 z n 2 a n z a n Proof: First note that since σ k+ Γ h < γ < σ k+ Γ h, P γ 2 I n is nonsingular It is also easy to see that: X =, X = 7 n+ n n+
8 Thus, P γ 2 I n A = γ 2 XA X I n = γ 2 n := uv and hence, using the matrix inversion lemma, A = P γ 2 I n uv = P γ 2 I n + λp γ 2 I n uv P γ 2 I n where we have defined Setting λ = v P γ 2 I n u = v ˆP u θ = θ θ 2 θ n = v P γ 2 I n shows that and hence A H = I n γ 2 I n P A X = γ 2 I n P A = I n λuθ λγ 2 θ + θ 2 λγ 2 θ 2 + θ 3 λγ 2 θ n + θ n λγ 2 θ n which is of the required form on noting that λ = Next consider B H Using the fact that γ 2 I n P A = θ u = γ 2 = θ γ 2 n i= i ˆP i + λγ 2 θ λγ 2 θ 2 λγ 2 θ n λγ 2 θ n the first column of B H can be written as: γ 2 I n P A B = h n + γ2 γ 2 θ n i= θ ih n i+ h n h 2 h 8
9 Similarly, since B = γxx C = γ n = γ the second column of B H can be written as: + λγ 2 θ λγ 2 θ 2 λγ 2 θ n λγ 2 θ n γ 2 I n P A B = γ = γ γ 2 θ as required Next consider the first row of C H ; this is given by: CX P A X = CX P P γ 2 XA X X = CX P γ 2 XA X + γ 2 XA X P γ 2 XA X X = C γ 2 CA X P γ 2 XA X X = C since CA = Finally, consider D H Using a similar argument to the one above, its, element is: Similarly, D H, = CX B CX P P γ 2 XA X B = CX B CX B γ 2 CA X P γ 2 XA X B = D H,2 = γ γcx P P γ 2 XA X XX C = γ γcx C γ 3 CA X P γ 2 XA X XX C = γ γcc = Finally consider the 2, element of D H, D H,2 = γ C A B We first show that A B is a zero column vector, except from its first element which is equal to h n Write A = P γ 2 XA X = x X 2 where x is a column vector and set where φ is a scalar Then: A B = φ φ 2 φ x + X 2 φ 2 = B It is easy that the first column of XA X is zero, and hence x is the first column of P which is given by x = h n h n h h n 9
10 Thus x is a multiple of B, from which it follows that φ = h n and φ 2 =, so that A B = h n Thus, D H,2 = γ γb X A B = γ γh n h n h n h n n+ = The transfer function of the central approximation follows by routine calculations since A H companion form is in In the remaining part of this section we derive a bound on the phase of the approximation error Note first, that in the scalar case, the solution to the Hankel-norm approximation problem is, in general, unique only in the optimal case The so-called central solution is obtained by setting Uz = Having obtained the central approximation, one next needs to remove its anti-stable component The extraction of the stable projection can be performed either by factoring the denominator polynomial of Xz, pz = z n a z n a 2 z n 2 a n z a n which is guaranteed to have k roots strictly inside the unit circle and n k roots strictly outside the unit circle Alternatively, the stable projection can be extracted from the state-space realization of Xz by transforming the system to block-schur form using an appropriate orthogonal state-space transformation and ordering the eigenvalues of A H in ascending order of magnitude; decoupling of the stable and anti-stable parts then requires the solution of a matrix Sylvester equation Reference [4] shows that having solved the k-th order Hankel-norm approximation problem, it is always possible to choose a constant term x so that the approximation error satisfies the bound: n hz + xz + x σ i Γ h 4 i=k+ Since this bound applies uniformly in frequency, it can be used to give an immediate bound on the approximation phase error: Theorem 5: Let ˆxz = xz + x be a k-th order optimal Hankel norm approximation of the scalar FIR filter hz such that 4 holds Then, if φω = arghe jω + argˆxe jω, we have that: sinφω n i=k+ σ iγ h he jω 5 for every ω [, π In particular, if the frequency interval [ω ω 2 ] lies in the filter s passband in which δ he jω + δ, then n i=k+ sinφω σ iγ h 6 δ for every ω [ω ω 2 ]
11 Proof: Straightforward and hence omitted The phase error bound given in Theorem 5 can be used to select the minimum order of approximation k consistent with a worst-case phase-error specification; If hz is a linear-phase FIR filter, then the RHS of 5 or 6 quantifies the deviation in the phase of the IIR approximation of hz from linearity If the magnitude tolerances resulting from the optimization procedure are tight, the specifications need to be tightened to account for the approximation error 6 Example In this section some of the results presented in the paper are illustrated by means of a computer example First, a linear phase FIR filter of order n = 2 was designed using the linear programming procedure outlined in section 3 The design specifications were defined as: ω p = rad/sample, ω s = 5 rads/sample and the frequency intervals [, ω p ] and [ω s, π] were each discretised using 5 equally-spaced frequencies The linear programme was then set up and solved using Matlab s function linprogm The minimum ripple was obtained as δ = 232; the first optimal impulse response coefficients of the resulting filter hz are tabulated below the last coefficients are symmetric and are not included: i hi i hi i hi Table : Impulse response coefficients The magnitude frequency response of the filter is shown in Figure below Note that the response in the passband and the stopband lies within the desired bounds ± δ Next the Hankel-norm approximation method outlined in Theorem 2 was applied to the FIR filter hz The Hankel singular values of hz are shown in Table 2 Parameter γ was chosen as γ = 3 resulting in a 7-th order sub-optimal Hankel-norm approximation The seven stable poles of the approximant were obtained as: 7467, 284 ± 852j, 4977 ± 6347j and 6886 ± 3363j i σ i i σ i i σ i i σ i Table 2: Hankel singular values
12 4 Optimal FIR filter and bounds 2 8 Gain angular frequency rads/sample Figure : Magnitude frequency response The magnitude response of hz and its IIR approximation is shown in Figure 2 It can be seen that the IIR filter has a slightly larger ripple in the pass-band Figure 3 shows the impulse responses of the two filters 4 FIR/IIR magnitude response 2 Magnitude linear angular frequency rads/sample Figure 2: FIR/IIR magnitude frequency response Finally, Figure 4 shows the phase error in the passband ω rad/sample arising from the approximation ie the deviation of the phase of the IIR filter from linearity, while Figure 5 shows the pole/zero pattern of the IIR filter there is an additional zero at z = 5652 not shown in the figure 2
13 4 FIR and IIR impulse response Amplitude index Figure 3: FIR/IIR impulse response 7 Conclusions In this work we have presented a systematic way for designing low-order IIR filters with approximately linear phase response characteristics The method relies on two steps: a First a linear-phase FIR filter is designed using linear programming or an alternative optimisation technique; b In step 2, the FIR filter is approximated by a low-order IIR filter using Hankel-norm approximation methods The method is computationally efficient and relies on the specialization of recent results for approximating discrete-time descriptor systems A-priori bounds are obtained for the magnitude and phase approximation error which can assist the designer to choose an approximation order consistent with the design specifications A low-order example has illustrated the effectiveness of the method There are a number of issues related to the proposed technique that we intend to pursue in the future These include: Derivation of tighter error bounds than those applying in general for the Hankel-norm approach This seems possible given the special structure of FIR filters all poles lie at the origin Investigation of the sensitivity properties of the proposed design under finite-precision implementation Generalisation of the method to other types of approximation, eg relative-error and general frequency weighted approximations Extension of the method to more general application domains in the fields of Signal Processing, Systems Identification and Control References [] M Al-Husari, IM Jaimoukha and DJN Limebeer 993, A descriptor approach for the solution of the one-block discrete distance problem, Proc IFAC World Congress, Sydney, Australia 3
14 2 Phase error in passband 5 Phase deg angular frequency rads/sample Figure 4: Phase error in passband deg [2] M Al-Husari, IM Jaimoukha and DJN Limebeer 997, A descriptor solution to a class of discrete distance problems, IEEE Trans Auto Control, 42: [3] X Chen and TW Parks 987, Design of FIR filters in the complex domain, IEEE Trans Acoust Speech and Signal Processing, 35:44-53 [4] K Glover 984, All optimal Hankel-norm approximations of linear multivariable systems and their L error bounds, Int J Control, 39:5-93 [5] GD Halikias, IM Jaimoukha and DA Wilson 997, A numerical solution to the matrix H /H 2 optimal control problem, Int J Robust and Nonlinear Control, 77:7-726 [6] AV Oppenheim and RW Schafer 975, Digital Signal Processing, Prentice Hall, Inc,, Englewood Cliffs, NJ [7] K Preuss 989, On the design of FIR filters by complex Chebychev approximations, IEEE Trans Acoust Speech and Signal Processing, 37:72-72 [8] LR Rabiner 972, Linear Program Design of Finite Impulse Response FIR Digital Filters, IEEE Tran Audio and Electroacoustics, AU24: [9] H Stark and FB Tuteur 979, Modern Electrical Communications, Prentice Hall, Inc,, Englewood Cliffs, NJ [] SP Wu, S Boyd and L Vandenberghe 997, FIR filter design via Spectral factorization and Convex Optimization, Applied Computational Control, Signal and Comunications, Biswa Datta ed, Birkhauser 4
15 5 Pole/zero plot 5 Imaginary axis Real axis Figure 5: Pole/zero pattern of IIR filter 5
DESIGN OF INFINITE IMPULSE RESPONSE (IIR) FILTERS WITH ALMOST LINEAR PHASE CHARACTERISTICS
DESIGN OF INFINITE IMPULSE RESPONSE (IIR) FILTERS WITH ALMOST LINEAR PHASE CHARACTERISTICS GD Halikias and IM Jaimoukha Control Engineering Research Centre, School of Engineering and Mathematical Sciences,
More informationPS403 - Digital Signal processing
PS403 - Digital Signal processing 6. DSP - Recursive (IIR) Digital Filters Key Text: Digital Signal Processing with Computer Applications (2 nd Ed.) Paul A Lynn and Wolfgang Fuerst, (Publisher: John Wiley
More informationVIBRATION CONTROL OF CIVIL ENGINEERING STRUCTURES VIA LINEAR PROGRAMMING
4 th World Conference on Structural Control and Monitoring 4WCSCM-65 VIBRATION CONTROL OF CIVIL ENGINEERING STRUCTURES VIA LINEAR PROGRAMMING P. Rentzos, G.D. Halikias and K.S. Virdi School of Engineering
More informationQuadrature-Mirror Filter Bank
Quadrature-Mirror Filter Bank In many applications, a discrete-time signal x[n] is split into a number of subband signals { v k [ n]} by means of an analysis filter bank The subband signals are then processed
More informationStatic Output Feedback Stabilisation with H Performance for a Class of Plants
Static Output Feedback Stabilisation with H Performance for a Class of Plants E. Prempain and I. Postlethwaite Control and Instrumentation Research, Department of Engineering, University of Leicester,
More informationOn the Frequency-Domain Properties of Savitzky-Golay Filters
On the Frequency-Domain Properties of Savitzky-Golay Filters Ronald W Schafer HP Laboratories HPL-2-9 Keyword(s): Savitzky-Golay filter, least-squares polynomial approximation, smoothing Abstract: This
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #24 Tuesday, November 4, 2003 6.8 IIR Filter Design Properties of IIR Filters: IIR filters may be unstable Causal IIR filters with rational system
More informationButterworth Filter Properties
OpenStax-CNX module: m693 Butterworth Filter Properties C. Sidney Burrus This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 3. This section develops the properties
More informationDigital Signal Processing Lecture 9 - Design of Digital Filters - FIR
Digital Signal Processing - Design of Digital Filters - FIR Electrical Engineering and Computer Science University of Tennessee, Knoxville November 3, 2015 Overview 1 2 3 4 Roadmap Introduction Discrete-time
More informationAn LMI Approach to the Control of a Compact Disc Player. Marco Dettori SC Solutions Inc. Santa Clara, California
An LMI Approach to the Control of a Compact Disc Player Marco Dettori SC Solutions Inc. Santa Clara, California IEEE SCV Control Systems Society Santa Clara University March 15, 2001 Overview of my Ph.D.
More informationComputational Methods for Feedback Control in Damped Gyroscopic Second-order Systems 1
Computational Methods for Feedback Control in Damped Gyroscopic Second-order Systems 1 B. N. Datta, IEEE Fellow 2 D. R. Sarkissian 3 Abstract Two new computationally viable algorithms are proposed for
More informationDESIGN OF QUANTIZED FIR FILTER USING COMPENSATING ZEROS
DESIGN OF QUANTIZED FIR FILTER USING COMPENSATING ZEROS Nivedita Yadav, O.P. Singh, Ashish Dixit Department of Electronics and Communication Engineering, Amity University, Lucknow Campus, Lucknow, (India)
More informationFilter Design Problem
Filter Design Problem Design of frequency-selective filters usually starts with a specification of their frequency response function. Practical filters have passband and stopband ripples, while exhibiting
More informationCONVEX OPTIMIZATION OVER POSITIVE POLYNOMIALS AND FILTER DESIGN. Y. Genin, Y. Hachez, Yu. Nesterov, P. Van Dooren
CONVEX OPTIMIZATION OVER POSITIVE POLYNOMIALS AND FILTER DESIGN Y. Genin, Y. Hachez, Yu. Nesterov, P. Van Dooren CESAME, Université catholique de Louvain Bâtiment Euler, Avenue G. Lemaître 4-6 B-1348 Louvain-la-Neuve,
More informationResearch Article Design of One-Dimensional Linear Phase Digital IIR Filters Using Orthogonal Polynomials
International Scholarly Research Network ISRN Signal Processing Volume, Article ID 8776, 7 pages doi:.54//8776 Research Article Design of One-Dimensional Linear Phase Digital IIR Filters Using Orthogonal
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #21 Friday, October 24, 2003 Types of causal FIR (generalized) linear-phase filters: Type I: Symmetric impulse response: with order M an even
More informationDIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS. 3.6 Design of Digital Filter using Digital to Digital
DIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS Contents: 3.1 Introduction IIR Filters 3.2 Transformation Function Derivation 3.3 Review of Analog IIR Filters 3.3.1 Butterworth
More informationMultimedia Signals and Systems - Audio and Video. Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2
Multimedia Signals and Systems - Audio and Video Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2 Kunio Takaya Electrical and Computer Engineering University of Saskatchewan December
More informationDesign of Stable IIR filters with prescribed flatness and approximately linear phase
Design of Stable IIR filters with prescribed flatness and approximately linear phase YASUNORI SUGITA Nagaoka University of Technology Dept. of Electrical Engineering Nagaoka city, Niigata-pref., JAPAN
More informationAn LQ R weight selection approach to the discrete generalized H 2 control problem
INT. J. CONTROL, 1998, VOL. 71, NO. 1, 93± 11 An LQ R weight selection approach to the discrete generalized H 2 control problem D. A. WILSON², M. A. NEKOUI² and G. D. HALIKIAS² It is known that a generalized
More informationOptimal Design of Real and Complex Minimum Phase Digital FIR Filters
Optimal Design of Real and Complex Minimum Phase Digital FIR Filters Niranjan Damera-Venkata and Brian L. Evans Embedded Signal Processing Laboratory Dept. of Electrical and Computer Engineering The University
More informationEFFICIENT REMEZ ALGORITHMS FOR THE DESIGN OF NONRECURSIVE FILTERS
EFFICIENT REMEZ ALGORITHMS FOR THE DESIGN OF NONRECURSIVE FILTERS Copyright 2003- Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 24, 2007 Frame # 1 Slide # 1 A. Antoniou EFFICIENT
More informationStability Condition in Terms of the Pole Locations
Stability Condition in Terms of the Pole Locations A causal LTI digital filter is BIBO stable if and only if its impulse response h[n] is absolutely summable, i.e., 1 = S h [ n] < n= We now develop a stability
More information( ) John A. Quinn Lecture. ESE 531: Digital Signal Processing. Lecture Outline. Frequency Response of LTI System. Example: Zero on Real Axis
John A. Quinn Lecture ESE 531: Digital Signal Processing Lec 15: March 21, 2017 Review, Generalized Linear Phase Systems Penn ESE 531 Spring 2017 Khanna Lecture Outline!!! 2 Frequency Response of LTI System
More informationMINIMUM PEAK IMPULSE FIR FILTER DESIGN
MINIMUM PEAK IMPULSE FIR FILTER DESIGN CHRISTINE S. LAW AND JON DATTORRO Abstract. Saturation pulses rf(t) are essential to many imaging applications. Criteria for desirable saturation profile are flat
More informationHankel Optimal Model Reduction 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Hankel Optimal Model Reduction 1 This lecture covers both the theory and
More informationFilter Analysis and Design
Filter Analysis and Design Butterworth Filters Butterworth filters have a transfer function whose squared magnitude has the form H a ( jω ) 2 = 1 ( ) 2n. 1+ ω / ω c * M. J. Roberts - All Rights Reserved
More informationEE 521: Instrumentation and Measurements
Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA November 1, 2009 1 / 27 1 The z-transform 2 Linear Time-Invariant System 3 Filter Design IIR Filters FIR Filters
More informationTitle without the persistently exciting c. works must be obtained from the IEE
Title Exact convergence analysis of adapt without the persistently exciting c Author(s) Sakai, H; Yang, JM; Oka, T Citation IEEE TRANSACTIONS ON SIGNAL 55(5): 2077-2083 PROCESS Issue Date 2007-05 URL http://hdl.handle.net/2433/50544
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationA DESIGN OF FIR FILTERS WITH VARIABLE NOTCHES CONSIDERING REDUCTION METHOD OF POLYNOMIAL COEFFICIENTS FOR REAL-TIME SIGNAL PROCESSING
International Journal of Innovative Computing, Information and Control ICIC International c 23 ISSN 349-498 Volume 9, Number 9, September 23 pp. 3527 3536 A DESIGN OF FIR FILTERS WITH VARIABLE NOTCHES
More informationDenis ARZELIER arzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15
More informationAnalysis Of Ill-Conditioning Of Multi-Channel Deconvolution Problems
Analysis Of Ill-Conditioning Of Multi-Channel Deconvolution Problems Ole Kirkeby, Per Rubak, and Angelo Farina * Department of Communication Techology, Fredrik Bajers Vej 7, Aalborg University, DK-9220
More informationMaximally Flat Lowpass Digital Differentiators
Maximally Flat Lowpass Digital Differentiators Ivan W. Selesnick August 3, 00 Electrical Engineering, Polytechnic University 6 Metrotech Center, Brooklyn, NY 0 selesi@taco.poly.edu tel: 78 60-36 fax: 78
More informationOptimal Polynomial Control for Discrete-Time Systems
1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should
More informationLike bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform.
Inversion of the z-transform Focus on rational z-transform of z 1. Apply partial fraction expansion. Like bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform. Let X(z)
More informationFilter structures ELEC-E5410
Filter structures ELEC-E5410 Contents FIR filter basics Ideal impulse responses Polyphase decomposition Fractional delay by polyphase structure Nyquist filters Half-band filters Gibbs phenomenon Discrete-time
More informationThe Approximation Problem
EE 508 Lecture 3 The Approximation Problem Classical Approximating Functions - Thompson and Bessel Approximations Review from Last Time Elliptic Filters Can be thought of as an extension of the CC approach
More informationEE 313 Linear Signals & Systems (Fall 2018) Solution Set for Homework #7 on Infinite Impulse Response (IIR) Filters CORRECTED
EE 33 Linear Signals & Systems (Fall 208) Solution Set for Homework #7 on Infinite Impulse Response (IIR) Filters CORRECTED By: Mr. Houshang Salimian and Prof. Brian L. Evans Prolog for the Solution Set.
More informationDigital Signal Processing 2/ Advanced Digital Signal Processing, Audio/Video Signal Processing Lecture 10, Frequency Warping, Example
Digital Signal Processing 2/ Advanced Digital Signal Processing, Audio/Video Signal Processing Lecture 10, Gerald Schuller, TU Ilmenau Frequency Warping, Example Example: Design a warped low pass filter
More informationUse: Analysis of systems, simple convolution, shorthand for e jw, stability. Motivation easier to write. Or X(z) = Z {x(n)}
1 VI. Z Transform Ch 24 Use: Analysis of systems, simple convolution, shorthand for e jw, stability. A. Definition: X(z) = x(n) z z - transforms Motivation easier to write Or Note if X(z) = Z {x(n)} z
More informationPart 4: IIR Filters Optimization Approach. Tutorial ISCAS 2007
Part 4: IIR Filters Optimization Approach Tutorial ISCAS 2007 Copyright 2007 Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 24, 2007 Frame # 1 Slide # 1 A. Antoniou Part4: IIR Filters
More informationCANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM
CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics, Universiteit Maastricht, P.O. Box 616, 6200 MD Maastricht,
More informationPoles and Zeros in z-plane
M58 Mixed Signal Processors page of 6 Poles and Zeros in z-plane z-plane Response of discrete-time system (i.e. digital filter at a particular frequency ω is determined by the distance between its poles
More informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More informationUsing fractional delay to control the magnitudes and phases of integrators and differentiators
Using fractional delay to control the magnitudes and phases of integrators and differentiators M.A. Al-Alaoui Abstract: The use of fractional delay to control the magnitudes and phases of integrators and
More informationH Optimal Nonparametric Density Estimation from Quantized Samples
H Optimal Nonparametric Density Estimation from Quantized Samples M. Nagahara 1, K. I. Sato 2, and Y. Yamamoto 1 1 Graduate School of Informatics, Kyoto University, 2 Graduate School of Economics, Kyoto
More informationKatholieke Universiteit Leuven Department of Computer Science
A convex optimization method to solve a filter design problem Thanh Hieu Le Marc Van Barel Report TW 599, August 2011 Katholieke Universiteit Leuven Department of Computer Science Celestijnenlaan 200A
More informationDiscrete-time Symmetric/Antisymmetric FIR Filter Design
Discrete-time Symmetric/Antisymmetric FIR Filter Design Presenter: Dr. Bingo Wing-Kuen Ling Center for Digital Signal Processing Research, Department of Electronic Engineering, King s College London. Collaborators
More informationUNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 9th, 011 Examination hours: 14.30 18.30 This problem set
More informationMinimax Design of Complex-Coefficient FIR Filters with Low Group Delay
Minimax Design of Complex-Coefficient FIR Filters with Low Group Delay Wu-Sheng Lu Takao Hinamoto Dept. of Elec. and Comp. Engineering Graduate School of Engineering University of Victoria Hiroshima University
More informationA METHOD OF ADAPTATION BETWEEN STEEPEST- DESCENT AND NEWTON S ALGORITHM FOR MULTI- CHANNEL ACTIVE CONTROL OF TONAL NOISE AND VIBRATION
A METHOD OF ADAPTATION BETWEEN STEEPEST- DESCENT AND NEWTON S ALGORITHM FOR MULTI- CHANNEL ACTIVE CONTROL OF TONAL NOISE AND VIBRATION Jordan Cheer and Stephen Daley Institute of Sound and Vibration Research,
More informationTopic 7: Filter types and structures
ELEN E481: Digital Signal Processing Topic 7: Filter tpes and structures 1. Some filter tpes 2. Minimum and maximum phase 3. Filter implementation structures 1 1. Some Filter Tpes We have seen the basics
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time
More informationThe Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation
The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse
More informationLinear Matrix Inequality (LMI)
Linear Matrix Inequality (LMI) A linear matrix inequality is an expression of the form where F (x) F 0 + x 1 F 1 + + x m F m > 0 (1) x = (x 1,, x m ) R m, F 0,, F m are real symmetric matrices, and the
More informationPrediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate
www.scichina.com info.scichina.com www.springerlin.com Prediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate WEI Chen & CHEN ZongJi School of Automation
More informationECE503: Digital Signal Processing Lecture 5
ECE53: Digital Signal Processing Lecture 5 D. Richard Brown III WPI 3-February-22 WPI D. Richard Brown III 3-February-22 / 32 Lecture 5 Topics. Magnitude and phase characterization of transfer functions
More informationLecture 7 Discrete Systems
Lecture 7 Discrete Systems EE 52: Instrumentation and Measurements Lecture Notes Update on November, 29 Aly El-Osery, Electrical Engineering Dept., New Mexico Tech 7. Contents The z-transform 2 Linear
More informationFREQUENCY DOMAIN IDENTIFICATION TOOLBOX
NASA Technical Memorandum 109039 FREQUENCY DOMAIN IDENTIFICATION TOOLBOX Lucas G. Horta and Jer-Nan Juang NASA Langley Research Center Hampton, Virginia 23681-0001 Chung-Wen Chen North Carolina State University
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationNew Design of Orthogonal Filter Banks Using the Cayley Transform
New Design of Orthogonal Filter Banks Using the Cayley Transform Jianping Zhou, Minh N. Do and Jelena Kovačević Department of Electrical and Computer Engineering University of Illinois at Urbana-Champaign,
More informationAnalysis of Finite Wordlength Effects
Analysis of Finite Wordlength Effects Ideally, the system parameters along with the signal variables have infinite precision taing any value between and In practice, they can tae only discrete values within
More informationIterative Learning Control Analysis and Design I
Iterative Learning Control Analysis and Design I Electronics and Computer Science University of Southampton Southampton, SO17 1BJ, UK etar@ecs.soton.ac.uk http://www.ecs.soton.ac.uk/ Contents Basics Representations
More informationTransform Representation of Signals
C H A P T E R 3 Transform Representation of Signals and LTI Systems As you have seen in your prior studies of signals and systems, and as emphasized in the review in Chapter 2, transforms play a central
More informationPartial Eigenvalue Assignment in Linear Systems: Existence, Uniqueness and Numerical Solution
Partial Eigenvalue Assignment in Linear Systems: Existence, Uniqueness and Numerical Solution Biswa N. Datta, IEEE Fellow Department of Mathematics Northern Illinois University DeKalb, IL, 60115 USA e-mail:
More informationSignals and Systems. Lecture 11 Wednesday 22 nd November 2017 DR TANIA STATHAKI
Signals and Systems Lecture 11 Wednesday 22 nd November 2017 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE LONDON Effect on poles and zeros on frequency response
More informationAnalog and Digital Filter Design
Analog and Digital Filter Design by Jens Hee http://jenshee.dk October 208 Change log 28. september 208. Document started.. october 208. Figures added. 6. october 208. Bilinear transform chapter extended.
More informationJ-SPECTRAL FACTORIZATION
J-SPECTRAL FACTORIZATION of Regular Para-Hermitian Transfer Matrices Qing-Chang Zhong zhongqc@ieee.org School of Electronics University of Glamorgan United Kingdom Outline Notations and definitions Regular
More informationDigital Signal Processing Lecture 8 - Filter Design - IIR
Digital Signal Processing - Filter Design - IIR Electrical Engineering and Computer Science University of Tennessee, Knoxville October 20, 2015 Overview 1 2 3 4 5 6 Roadmap Discrete-time signals and systems
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #20 Wednesday, October 22, 2003 6.4 The Phase Response and Distortionless Transmission In most filter applications, the magnitude response H(e
More informationT ulation and demodulation. It has applications not only
IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT. VOL 39. NO 6. DECEMBER 1990 847 Optimal FIR and IIR Hilbert Transformer Design Via LS and Minimax Fitting ISTVAN KOLLAR, MEMBER, IEEE, RIK PINTELON,
More informationDiscrete Time Systems
Discrete Time Systems Valentina Hubeika, Jan Černocký DCGM FIT BUT Brno, {ihubeika,cernocky}@fit.vutbr.cz 1 LTI systems In this course, we work only with linear and time-invariant systems. We talked about
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
More informationRank-one LMIs and Lyapunov's Inequality. Gjerrit Meinsma 4. Abstract. We describe a new proof of the well-known Lyapunov's matrix inequality about
Rank-one LMIs and Lyapunov's Inequality Didier Henrion 1;; Gjerrit Meinsma Abstract We describe a new proof of the well-known Lyapunov's matrix inequality about the location of the eigenvalues of a matrix
More informationAn Fir-Filter Example: Hanning Filter
An Fir-Filter Example: Hanning Filter Josef Goette Bern University of Applied Sciences, Biel Institute of Human Centered Engineering - microlab Josef.Goette@bfh.ch February 7, 2018 Contents 1 Mathematical
More informationIIR digital filter design for low pass filter based on impulse invariance and bilinear transformation methods using butterworth analog filter
IIR digital filter design for low pass filter based on impulse invariance and bilinear transformation methods using butterworth analog filter Nasser M. Abbasi May 5, 0 compiled on hursday January, 07 at
More informationThe Discrete-Time Fourier
Chapter 3 The Discrete-Time Fourier Transform 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 3-1-1 Continuous-Time Fourier Transform Definition The CTFT of
More informationTransform Analysis of Linear Time-Invariant Systems
Transform Analysis of Linear Time-Invariant Systems Discrete-Time Signal Processing Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-Sen University Kaohsiung, Taiwan ROC Transform
More informationUsing Hankel structured low-rank approximation for sparse signal recovery
Using Hankel structured low-rank approximation for sparse signal recovery Ivan Markovsky 1 and Pier Luigi Dragotti 2 Department ELEC Vrije Universiteit Brussel (VUB) Pleinlaan 2, Building K, B-1050 Brussels,
More informationFast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma
Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent
More informationAPPLIED SIGNAL PROCESSING
APPLIED SIGNAL PROCESSING DIGITAL FILTERS Digital filters are discrete-time linear systems { x[n] } G { y[n] } Impulse response: y[n] = h[0]x[n] + h[1]x[n 1] + 2 DIGITAL FILTER TYPES FIR (Finite Impulse
More informationAPPROXIMATE GREATEST COMMON DIVISOR OF POLYNOMIALS AND THE STRUCTURED SINGULAR VALUE
PPROXIMTE GRETEST OMMON DIVISOR OF POLYNOMILS ND THE STRUTURED SINGULR VLUE G Halikias, S Fatouros N Karcanias ontrol Engineering Research entre, School of Engineering Mathematical Sciences, ity University,
More informationREAL TIME DIGITAL SIGNAL PROCESSING
REAL TIME DIGITAL SIGNAL PROCESSING www.electron.frba.utn.edu.ar/dplab Digital Filters FIR and IIR. Design parameters. Implementation types. Constraints. Filters: General classification Filters: General
More informationDSP. Chapter-3 : Filter Design. Marc Moonen. Dept. E.E./ESAT-STADIUS, KU Leuven
DSP Chapter-3 : Filter Design Marc Moonen Dept. E.E./ESAT-STADIUS, KU Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/stadius/ Filter Design Process Step-1 : Define filter specs Pass-band, stop-band,
More information-Digital Signal Processing- FIR Filter Design. Lecture May-16
-Digital Signal Processing- FIR Filter Design Lecture-17 24-May-16 FIR Filter Design! FIR filters can also be designed from a frequency response specification.! The equivalent sampled impulse response
More informationTransform analysis of LTI systems Oppenheim and Schafer, Second edition pp For LTI systems we can write
Transform analysis of LTI systems Oppenheim and Schafer, Second edition pp. 4 9. For LTI systems we can write yœn D xœn hœn D X kd xœkhœn Alternatively, this relationship can be expressed in the z-transform
More informationFIR BAND-PASS DIGITAL DIFFERENTIATORS WITH FLAT PASSBAND AND EQUIRIPPLE STOPBAND CHARACTERISTICS. T. Yoshida, Y. Sugiura, N.
FIR BAND-PASS DIGITAL DIFFERENTIATORS WITH FLAT PASSBAND AND EQUIRIPPLE STOPBAND CHARACTERISTICS T. Yoshida, Y. Sugiura, N. Aikawa Tokyo University of Science Faculty of Industrial Science and Technology
More informationVISA: Versatile Impulse Structure Approximation for Time-Domain Linear Macromodeling
VISA: Versatile Impulse Structure Approximation for Time-Domain Linear Macromodeling Chi-Un Lei and Ngai Wong Dept. of Electrical and Electronic Engineering The University of Hong Kong Presented by C.U.
More informationModel reduction via tangential interpolation
Model reduction via tangential interpolation K. Gallivan, A. Vandendorpe and P. Van Dooren May 14, 2002 1 Introduction Although most of the theory presented in this paper holds for both continuous-time
More informationReview of Discrete-Time System
Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationAn Optimum Fitting Algorithm for Generation of Reduced-Order Models
An Optimum Fitting Algorithm for Generation of Reduced-Order Models M.M. Gourary 1, S.G. Rusakov 1, S.L. Ulyanov 1, M.M. Zharov 1, B.J. Mulvaney 2 1 IPPM, Russian Academy of Sciences, Moscow 1523, e-mail:
More informationThe Approximation Problem
EE 508 Lecture The Approximation Problem Classical Approximating Functions - Elliptic Approximations - Thompson and Bessel Approximations Review from Last Time Chebyshev Approximations T Type II Chebyshev
More informationSYNTHESIS OF BIRECIPROCAL WAVE DIGITAL FILTERS WITH EQUIRIPPLE AMPLITUDE AND PHASE
SYNTHESIS OF BIRECIPROCAL WAVE DIGITAL FILTERS WITH EQUIRIPPLE AMPLITUDE AND PHASE M. Yaseen Dept. of Electrical and Electronic Eng., University of Assiut Assiut, Egypt Tel: 088-336488 Fax: 088-33553 E-Mail
More informationModel Reduction for Unstable Systems
Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model
More information1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0
4 INUTES. If, ω, ω, -----, ω 9 are the th roots of unity, then ( + ω) ( + ω ) ----- ( + ω 9 ) is B) D) 5. i If - i = a + ib, then a =, b = B) a =, b = a =, b = D) a =, b= 3. Find the integral values for
More informationSpectral Theorem for Self-adjoint Linear Operators
Notes for the undergraduate lecture by David Adams. (These are the notes I would write if I was teaching a course on this topic. I have included more material than I will cover in the 45 minute lecture;
More informationLMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System
LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System Gholamreza Khademi, Haniyeh Mohammadi, and Maryam Dehghani School of Electrical and Computer Engineering Shiraz
More informationInterval solutions for interval algebraic equations
Mathematics and Computers in Simulation 66 (2004) 207 217 Interval solutions for interval algebraic equations B.T. Polyak, S.A. Nazin Institute of Control Sciences, Russian Academy of Sciences, 65 Profsoyuznaya
More information