Minimization of Frequency-Weighted l 2 -Sensitivity Subject to l 2 -Scaling Constraints for Two-Dimensional State-Space Digital Filters
|
|
- Sheena Daniels
- 5 years ago
- Views:
Transcription
1 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ 007 Minimization of Frequency-Weighted l -Sensitivity Subject to l -Scaling Constraints for Two-Dimensional State-Space Digital Filters Taao Hinamoto, Toru Oumi, Osemehian I. Omoifo and Wu-Sheng Lu Abstract This paper investigates the problem of frequencyweighted l -sensitivity minimization subject to l -scaling constraints for two-dimensional (-D) state-space digital filters described by the Roesser model. It is shown that the Fornasini- Marchesini second model can be imbedded in the Roesser model. Two iterative methods are developed to solve the constrained optimization problem encountered. The first iterative method introduces a Lagrange function and optimizes it using some matrix-theoretic techniques and an efficient bisection method. The second iterative method converts the problem into an unconstrained optimization formulation by using linear-algebraic techniques and solves it by applying an efficient quasi-newton algorithm. The optimal filter structure with minimum frequencyweighted l -sensitivity and no overflow is then synthesized by an appropriate coordinate transformation. Case studies are presented to demonstrate the validity and effectiveness of the proposed techniques. Index Terms -D digital filters, Roesser s model, Fornasini- Marchesini s second model, frequency-weighted l -sensitivity minimization, l -scaling constraints, no overflow, Lagrange function, bisection method, quasi-newton method I. INTRODUCTION It is well nown that there exist infinitely many minimal state-space realizations for a given transfer function, and some inherent properties such as controllability, observability, stability, etc. are invariant within these realizations. However, performance measures such as coefficient sensitivity, output roundoff noise, overflow oscillations, etc. may be significantly varying among the realizations. Consider a transfer function with coefficients of infinite accuracy, which meets certain design specifications including stability. When the transfer function is implemented by a state-space model with a finite binary representation, truncation or rounding of the state-space model is required to satisfy the finite word length (FWL) constraints. As a result, the characteristics of the stable filter might be so altered that the filter may become unstable. This motivates the study of the coefficient sensitivity minimization problem. To date, several techniques have been reported for synthesizing the state-space descriptions with minimum coefficient sensitivity. The techniques can be divided into two Manuscript received mm dd, 007; revised mm dd, 007. T. Hinamoto, T. Oumi and O. I. Omoifo are with the Graduate School of Engineering, Hiroshima University, Higashi-Hiroshima , Japan. ( {hinamoto,oumi,osei}@hiroshima-u.ac.jp, Phone: , Fax: ) W.-S. Lu is with the Depertment of Electrical and Computer Engineering, University of Victoria, Victoria, B.C, Canada, V8W 3P6. ( wslu@ece.uvic.ca, Phone: , Fax: ) main classes: those for l /l -mixed sensitivity minimization -5 and those for l -sensitivity minimization 6-. In 6-0, it has been argued that the sensitivity measure based on a sole l -norm is more natural and reasonable relative to the l /l -mixed sensitivity minimization. For -D statespace digital filters, the l /l -mixed sensitivity minimization problem -7 and l -sensitivity minimization problem 0,7-0 have also been investigated. It has been realized that solutions for frequency-weighted sensitivity minimization would be of practical use as these solutions allow to emphasize or de-emphasize the filter s sensitivity in certain frequency regions of interest. Synthesis procedures of the optimal FWL -D filter structures that minimize the frequency-weighted sensitivity measure have been considered 5-8. However, the minimization methods proposed in the above wor do not impose constraints on the scaling of the design variables. As a result, elimination of overflow cannot be ensured. More recently, the minimization problem of l -sensitivity subject to l -scaling constraints has been explored for -D and a class of -D state-space digital filters -3. It is well nown that the use of scaling constraints can be beneficial for suppressing overflow 4,5. However, frequency-weighted sensitivity measure has not yet been considered in -3. In this paper, we investigate the problem of minimizing a frequency-weighted l -sensitivity measure subject to l -scaling constraints for -D state-space digital filters described by the Roesser local state-space (LSS) model 6. We then proceed by introducing an expression for evaluating the frequencyweighted l -sensitivity, and formulating the minimization problem for the frequency-weighted l -sensitivity measure subject to l -scaling constraints. Next, two iterative methods are developed for solving the constrained optimization problem. The first iterative method introduces a Lagrange function, and maes use of some matrix-theoretic techniques and an efficient bisection method. The second iterative method relies on a technique that converts the constrained optimization problem into an unconstrained optimization formulation and utilizes an efficient quasi-newton method with closed-form formula for gradient evaluation. Finally, case studies are presented to demonstrate the validity and effectiveness of the proposed techniques. One of the contributions made in this paper is to show that either the Fornasini-Marchesini (FM) second LSS model 7 or its transposed-structure model,8 can be imbedded in the Roesser model as a special case. This justifies the use of the Roesser model in our studies. Another
2 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ 007 contribution is that a bisection method is applied to obtain the Lagrange multipliers, which maes it possible to attain considerably faster convergence than the algorithm reported in. Moreover, unlie and, the present paper investigates a frequency-weighted l -sensitivity measure under l -scaling constraints. Although this extension is technically manageable, to the best of our nowledge, this is the first time a frequency-weighted l -sensitivity measure under l -scaling constraints is addressed in the l /l framewor for state-space digital filters. Throughout the paper, I n stands for the identity matrix of dimension n n, is used to denote the direct sum of matrices, the transpose (conjugate transpose) of a matrix A is indicated by A T (A ), and the trace and ith diagonal element of a square matrix A are denoted by tra and (A) ii, respectively. II. PROBLEM FORMULATION. System Models Consider a stable, separately locally controllable and separately locally observable LSS model for -D recursive digital filters x h (i +,j) A A x h (i, j) b + u(i, j) x v (i, j +) A 3 A 4 x v (i, j) b y(i, j) x h (i, j) c c + du(i, j) x v (i, j) () which was originally proposed by Roesser 6,9, x h (i, j) is an m horizontal state vector, x v (i, j) is an n vertical state vector, u(i, j) is a scalar input, y(i, j) is a scalar output, and A, A, A 3, A 4, b, b, c, c,andd are m m, m n, n m, n n, m, n, m, n, and real constant matrices, respectively. A bloc diagram of the LSS model in (39) is shown in Fig.. The Fig.. The bloc diagram of the Roesser LSS model. transfer function of the LSS model in (39) is given by H(z,z )c(z A) b + d () A, b, andc are (m + n) (m + n), (m + n), and (m + n) real constant matrices defined by A A b A, b, c c c, A 3 A 4 b respectively, and Z z I m z I n. For the sae of simplicity, the LSS model in (39) is denoted hereafter by (A, b, c,d) m,n. Alternatively, an LSS model for a class of -D recursive digital filters can be described by,8 x(i +,j+) A A x(i, j +) y(i, j) c c x(i +,j) (3) b + u(i, j) d x(i, j) is an N local state vector, u(i, j) is a scalar input, y(i, j) is a scalar output, and A, A, b, c, c,andd are N N, N N, N, N, N, and real constant matrices, respectively. The transfer function of the LSS model in (3) is given by D(z,z )(z c + z c ) If we define (I n z A z A ) b + d. (4) x(i, j +)x h (i, j), x(i +,j)x v (i, j) (5) then the LSS model in (3) can then be imbedded in that of (39) as a special case as follows: x h (i +,j) A A x h (i, j) b + u(i, j) x v (i, j +) A A x v (i, j) b y(i, j) x h (i, j) c c + du(i, j) x v (i, j) (6) m n N. It is noted that D(z,z ) T can be viewed as a transfer function of the FM second LSS model 7, which reveals that the LSS model of D(z,z ) T can be realized by a transposed structure of that in (6). Therefore, the LSS model in (39) is more general than either the LSS model in (3) or the FM second LSS model 7 (and vice versa with the same dimension, but increased number of coefficients 7). In addition, we note that the technique reported in has merely treated the l -sensitivity minimization problem for the LSS model in (3) subject to l -scaling constraints. Recall that the total numbers of the coefficients in (39) and (3) are (m + n) +(m + n)+ and N +3N +, respectively. This means that the LSS model in (39) has less number of the coefficients than that of (3) when their local state vectors possess the same dimension, i.e., m + n N. Under these circumstances, it is worthwhile to consider the more general problem of minimizing the frequency-weighted l -sensitivity for the LSS model in (39) subject to l -scaling constraints, because its solutions will mae it possible to emphasize or deemphasize the filter s sensitivity in certain frequency regions of interest, and because the LSS model in (39) owns less number of the coefficients in case m + n N.. A Frequency-Weighted l -Sensitivity Measure Suppose that the LSS model in (39) is implemented by FWL fixed-point arithmetic with a B-bit fractional representation, and is realized with coefficient matrices
3 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ Ã A +ΔA, b b +Δb (7) c c +Δc, d d +Δd ΔA, Δb, Δc, andδd stand for the quantization errors of the coefficient matrices. The transfer function of the FWL realization is then expressed as H(z,z ) c(z Ã) b + d. (8) Let {p i,i,,,m} be the set of the ideal parameters of a realization and let { p i,i,,,m} be its FWL version p i p i +Δp i with Δp i indicating the corresponding parameter perturbation. If all Δp i are sufficiently small in magnitude, then the first-order approximation of the Taylor series expansion yields ΔH(z,z ) H(z,z ) H(z,z ) M H(z,z ) Δp i. p i i Obviously, smaller H(z,z )/ p i for i,,,m yield smaller transfer function error ΔH(z,z ). For a fixed-point implementation of B bits, the parameter perturbations can be considered to be independent random-variables uniformly distributed within the range B, B. Under these circumstances, a measure of the transfer function error can statistically be defined as σ ΔH (πj) z z (9) E ΔH(z,z ) dz dz z z (0) E( ) denotes the ensemble-average operation. Since Δp i s are uniformly-distributed independent random variables, it follows that E ΔH(z,z ) M i H(z,z ) p i σ E(Δp i ) B. Thus, if the measure S o (πj) z M z i H(z,z ) p i σ () dz dz () z z is minimized then the minimum variance σδh can be attained because of the relation σδh S oσ. The measure in () is referred to as an l -sensitivity measure. The frequency-weighted l -sensitivity of the LSS model in (39) is defined as follows. Definition :LetX be an m n real matrix and let f(x) be a scalar complex function of X, that is differentiable with respect to all the entries of X. The sensitivity function of f(x) with respect to X is then defined as 5 S X f(x) X, (S X ) ij f(x) (3) x ij x ij denotes the (i, j)th entry of matrix X. Definition : In order to tae into account the sensitivity of the transfer function in a specified frequency band, or even at some discrete frequency points, the weighted sensitivity functions are defined as 5,7 δh(z,z ) W A (z,z ) H(z,z ) δa A δh(z,z ) W B (z,z ) H(z,z ) δb b δh(z,z ) δc T W C (z,z ) H(z,z ) c T (4) W A (z,z ), W B (z,z ), and W C (z,z ) are scalar, stable, causal functions of the complex variables z and z. Notice that δ in (4) is not meant to be a derivative operator, but rather a notation for defining the weighted parameter sensitivity. Definition 3 :LetX(z,z ) be an m n complex matrix valued function of the complex variables z and z.thel norm of X(z,z ) is then defined by X(z,z ) ( tr (πj) Γ X(z,z )X (z,z ) dz dz Γ z z ) (5) j and Γ i {z i : z i } for i,. From () and Definitions -3, the overall frequencyweighted l -sensitivity measure for the LSS model in (39) can be defined as S δh(z,z ) δa + δh(z,z ) δb + δh(z,z ) δc T W A (z,z )F (z,z )G(z,z ) T + W B (z,z )G T (z,z ) + W C(z,z )F (z,z ) (6) F (z,z )(Z A) b, G(z,z )c (Z A). It follows that the frequency-weighted l -sensitivity measure in (6) can be written as S trm A +trw B +trk C (7) M A, W B,andK C are obtained by the following general expression: X (πj) Y (z,z )Y (z,z ) dz dz Γ z z Γ with Y (z,z )W A (z,z )F (z,z )G(z,z ) T for X M A, Y (z,z )W B (z,z )G T (z,z ) for X W B,and Y (z,z )W C (z,z )F (z,z ) for X K C..3 Problem Formulation Define a state-space coordinate transformation by 6,9 x h (i, j) T 0 x h (i, j) (8) x v (i, j) 0 T 4 x v (i, j)
4 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ T and T 4 are m m and n n nonsingular matrices, respectively. New realizations can then be characterized as (A, b, c,d) m,n with A T AT, b T b, c ct (9) T T T 4. For a new realization, the frequencyweighted l -sensitivity measure in (7) is changed to S(P )trm A (P )P +trw B P +trk C P trn A (P )P +trw B P +trk C P (0) with P TT T P P 4 and M A (P ) (πj) Y (z,z )P Y (z,z ) dz dz Γ z z N A (P ) (πj) Γ Γ Y (z,z )PY(z,z ) dz dz Γ z z Y (z,z )W A (z,z )F (z,z )G(z,z ) T. If l -scaling constraints are imposed on the horizontal and vertical state vectors x h (i, j) and x v (i, j), we require that 30 (K ) ξξ (T K T T ) ξξ for ξ,,,m (K 4 ) ζζ (T 4 K 4T T 4 ) ζζ for ζ,,,n () K (πj) F (z,z )F (z,z ) dz dz Γ Γ z z K K K 3 K 4 is the local controllability Gramian for the LSS model in (39) with an m m submatrix K and an n n submatrix K 4 along its diagonal 9. Thus, the l -scaling constrained frequency-weighted l - sensitivity minimization problem can be formulated as follows: Given matrices A, b, and c, obtain a bloc-diagonal nonsingular matrix T T T 4 which minimizes S(P ) in (0) subject to l -scaling constraints in (). III. PROBLEM SOLUTION 3. A Constrained Optimization Method Solving the optimization problem formulated above consists of several steps. First, we relax the problem of minimizing S(P ) in (0) subject to l -scaling constraints in () into the problem minimize S(P ) in (0) subject to trk P m and trk 4P 4 n. () If trk P m(trk 4P 4 n) is satisfied, then an m m (n n) orthogonal matrix U (U 4 ) matrix can always be constructed so that T P / U (T 4 P / 4 U 4 ) satisfies l -scaling constraints in (). This justifies the relaxation made in (). In order to solve problem () for P P P 4,wedefine the Lagrange function of the problem as J(P,λ,λ 4 ) trm (P )P +trm 4 (P )P 4 +trw B P + trw 4B P 4 +trk C P +trk 4CP 4 (3) + λ (trk P m)+λ 4 (trk 4 P 4 n) λ and λ 4 are the Lagrange multipliers, and M (P ) M (P ) M A (P ) M 3 (P ) M 4 (P ) W B W B KC K C W B, K C W 3B W 4B K 3C K 4C with an m m submatrix and an n n submatrix along its diagonal for each matrix. Using the formula for evaluating matrix gradient 3, p.75 tr(mx) M T X (4) tr(mx ) (X MX ) T X we compute J(P,λ,λ 4 ) P J(P,λ,λ 4 ) P 4 M (P ) P N (P )P + W B P K C P λ P K P M 4 (P ) P 4 N 4 (P )P 4 + W 4B P 4 K 4C P 4 λ 4 P 4 K 4 P 4 (5) N (P ) N (P ) N A (P ) N 3 (P ) N 4 (P ) with an m m submatrix N (P ) and an n n submatrix N 4 (P ) along its diagonal. Setting J(P,λ,λ 4 )/ P 0 and J(P,λ,λ 4 )/ P 4 0, it follows that P F (P )P G (P,λ ) (6) P 4 F 4 (P )P 4 G 4 (P,λ 4 ) F (P )M (P )+W B, F 4 (P )M 4 (P )+W 4B G (P,λ )N (P )+K C + λ K G 4 (P,λ 4 )N 4 (P )+K 4C + λ 4 K 4. The equations in (6) are highly nonlinear with respect to P and P 4. Namely, P i F i (P )P i for i, 4 has a rational type R/D of nonlinearity, the degree of the nonlinearity of R is m + n + and the degree of the nonlinearity of D is m+n, while G i (P,λ i ) for i, 4 depends on P linearly. An effective approach for solving these equations is to relax them into the following recursive second-order matrix equations: P (+) F (P () )P (+) G (P (),λ (+) ) (7) P (+) 4 F 4 (P () )P (+) 4 G 4 (P (),λ (+) 4 )
5 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ with initial condition P (0) P (0) P (0) 4 I m+n. Noting that PWP M has the unique solution 5 P W W MW W (8) W > 0 and M 0 are symmetric, the unique solutions P (+) and P (+) 4 of (7) are found to be P (+) F (P () ) F (P () ) G (P (),λ (+) )F (P () ) F (P () ) P (+) 4 F 4 (P () ) F 4 (P () ) G 4 (P (),λ (+) 4 )F 4 (P () ) F 4 (P () ). (9) Here, the Lagrange multipliers λ (+) and λ (+) 4 can be efficiently obtained using a bisection method 3 so that f (λ (+) () )m tr K G () (λ (+) ) 0 f 4 (λ (+) () 4 ) n tr K 4 G () 4 (λ (+) 4 ) 0 are satisfied K () F (P () ) K F (P () ) K () 4 F 4 (P () ) K 4 F 4 (P () ) (30) G () (λ (+) )F (P () ) G (P (),λ (+) )F (P () ) G () 4 (λ(+) 4 )F 4 (P () ) G 4 (P (),λ (+) 4 )F 4 (P () ). The iteration process continues until J(P (),λ (+),λ (+) 4 ) J(P ( ),λ (),λ() 4 ) <ε (3) for a prescribed tolerance ε>0. If the iteration is terminated at step, thenp () is claimed to be a solution point. It is noted that a straightforward extension of the iterative algorithm reported in for updating λ and λ 4 does not wor well in this case, we do not require trkp m+n, but both trk P m and trk 4P 4 n. In addition, the bisection method offers an exponential convergence rate of (/) L L is the number of iterations used. As such, accurate solutions of the Lagrange multipliers λ (+) and λ (+) 4 in (30) can be identified with just a few iterations. In our simulation studies, the bisection method was found considerably faster than the iterative algorithm proposed in. All the Gramians can be evaluated by truncating the corresponding infinite summations, see Appendix I for details. 3. An Unconstrained Optimization Method By defining ˆT ˆT ˆT 4 (T T 4 ) T (K K 4 ), (3) it follows that T i K i T i T T ˆT i ˆT i for i, 4. (33) Thus, a convenient way to eliminate the l -scaling constraints in () is to choose ˆT and ˆT 4 as ˆT t t, t t,, t m t m (34) ˆT 4 t 4 t 4, t 4 t 4,, t 4n t 4n t i is an m vector for i,,,m and t 4j is an n vector for j,,,n. With (34), all the T diagonal elements of ˆT ˆT T and ˆT 4 ˆT 4 are found to be unity. Taing (3) and (34) into account, we conclude that the coordinate transformation matrix T of the form T T T 4 (K K 4 ) ( ˆT ˆT 4 ) T (35) automatically satisfies the l -scaling constraints in (). By substituting (35) into S(P ) in (0), the frequency-weighted l -sensitivity measure can be expressed as J o (x) tr ˆT ˆM A ( ˆP ) ˆT T +tr ˆT Ŵ B ˆT T +tr ˆT T ˆKC ˆT (36) with x (t T, tt,, tt m, tt 4, tt 4,, tt 4n )T ˆM A ( ˆP ) (πj) ˆP ˆT T ˆT and Γ Ŷ (z,z ) ˆP Ŷ (z,z ) dz dz Γ z z Ŷ (z,z )(K K 4 ) Y (z,z )(K K 4 ) Y (z,z )W A (z,z )F (z,z )G(z,z ) T Ŵ B (K K 4 ) W B (K K 4 ) ˆK C (K K 4 ) K C (K K 4 ). This shows that the problem of obtaining the bloc-diagonal nonsingular matrix T T T 4 which minimizes S(P ) in (0) subject to the l -scaling constraints in () can be converted into an unconstrained optimization problem of obtaining an (m + n ) vector x which minimizes J o (x) in (36). Applying a quasi-newton algorithm to minimize J o (x) in (36), in the th iteration the most recent point x is updated to point x + as 33 d S J o (x ), S + S + S 0 I m +n, x + x + α d (37) ( + γ T S γ γ T δ α arg min J o(x + αd ) α ) δ δ T γ T δ δ γ T S +S γ δ T γ T δ δ x + x γ J o (x + ) J o (x ). Here, J o (x) is the gradient of J o (x) with respect to x, and S is a positive-definite approximation of the inverse Hessian matrix of J o (x). The algorithm starts with a trivial initial point
6 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ x 0 obtained from an initial assignment ˆT I m+n, and this iteration process continues until J o (x + ) J o (x ) <ε (38) ε>0 is a prescribed tolerance. The implementation of (37) requires the gradient of J o (x), which can be efficiently evaluated using closed-form expressions, see Appendix II for details. IV. CASE STUDIES In this section, we examine the proposed algorithms by applying them to a recursive -D state-space digital filter. In addition, the algorithms are applied to a -D filter utilized in and the results are compared with those obtained by the method of. Example : Consider a -D stable recursive digital filter realization (A o, b o, c o,d), 30, p A o b o T c o d After carrying out the l -scaling for the above realization with a diagonal coordinate matrix T o diag{9.3364, ,.0650, }, we obtain the -D state-space digital filter (A, b, c,d), characterized by A b T c d A T o A o T o, b T o b o,andc c o T o.the -D state-space digital filter (A, b, c,d), now satisfies the l -scaling constraints. The frequency-weighting functions used in this example were given by a -D nonrecursive low-pass digital filter with the following unit-sample response 34, p. 895: w A (i, j) w B (i, j) w C (i, j) exp {(i 4) +(j 4) } for (0, 0) (i, j) (0, 0), and zero else. Using (A.) and (A.) with truncation (0, 0) (i, j) (50, 50) to evaluate the Gramians K C, W B, M A and K, it was found that K C W B M A K In what follows, we evaluate the frequency-weighted l - sensitivity of the -D state-space digital filter (A, b, c,d),. Later on, this sensitivity measure will be used as a benchmar in the examination of the performance of the proposed algorithms. Using (7), the frequency-weighted l -sensitivity of the LSS model (A, b, c,d), was found to be S As will be seen next, the proposed algorithms are able to deduce an equivalent state-space realization with much reduced frequency-weighted l -sensitivity. 4. Application of the Lagrange method Choosing P (0) P (0) P (0) 4 I 4 in (9) as an initial estimate and a tolerance ε 0 8 in (3) as well as in the bisection method, it too the Lagrange-based algorithm 0 iterations to converge to the solution P opt or equivalently, T opt The minimized frequency-weighted l -sensitivity measure in (3) corresponding to the above solution was found to be J(P opt,λ,λ 4 ) with λ and λ , and the optimal state-space filter structure (A, b, c,d), (that
7 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ minimizes (0) subject to the l -scaling constraints in ()) was synthesized by substituting matrix T opt into (9) as A b T c The profile of the frequency-weighted l -sensitivity measure J(P,λ,λ 4 ) and the profiles of the Lagrange multipliers λ and λ 4 for the first 0 iterations are shown in Figs. and 3, respectively. 5 x 0 4 or equivalently, T opt and the minimized frequency-weighted l -sensitivity was found to be J o ( ˆT opt ) We see that the results obtained by this method are identical to those obtained by the Lagrange-based method. The profile of the l -sensitivity measure J o ( ˆT ) during the first 6 iterations is shown in Fig x (P (), λ (+), λ 4 (+) ) J Fig.. Profile of J(P,λ,λ 4 ) during the first 0 iterations. J o (x ) Fig. 4. Profile of J o ( ˆT ) during the first 6 iterations. λ (), λ4 ().5 x λ λ Fig. 3. Profiles of λ and λ 4 during the first 0 iterations. 4. Application of the quasi-newton method By choosing ˆT I I (therefore T (K K 4 ) / in (3)) as an initial estimate and a tolerance ε 0 8 in (38), the quasi-newton algorithm too 6 iterations to converge to the solution ˆT opt We now explain the effectiveness of the optimal realization that minimizes the frequency-weighted l -sensitivity subject to l -scaling constraints. The magnitude response of the original realization (A, b, c,d), is shown in Fig. 5, the maximum value and the l -norm (which corresponds to a square root of the summation of squared values at 0 0 sampling points) were and 33.50, respectively. When all coefficients in the original realization were rounded to power-of-two representation with 8 bits after binary point, the magnitude-response deviation between the original realization and that with rounded coefficients is shown in Fig. 6, the maximum deviation and the l -norm were 4.44 and.5633, respectively. Alternatively, when all coefficients in the optimal realization were rounded in the same manner, the magnitude-response deviation between the original realization and the optimal one with rounded coefficients is shown in Fig. 7, the maximum deviation and the l -norm were and.3380, respectively. From these figures and data, it is observed that the coefficient sensitivity of the optimal realization is considerably lower than that of the original realization. Example : Let a -D stable recursive digital filter realization be specified by (A o, b o, c o,d) 4,4 A o A o A o b o, b o, c o c o A o A o b o c o
8 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ Fig. 5. Magnitude response of the original realization. Fig. 7. Magnitude-response deviation between the original realization and the optimal one with rounded coefficients. and then applying the same frequency-weighted functions as in Example to the resulting realization (A, b, c,d) 4,4,the frequency-weighted l -sensitivity in (7) was found to be S Fig. 6. Magnitude-response deviation between the original realization and that with rounded coefficients. with A o A o b o bo T c o c o d This -D filter was obtained by imbedding the LSS model of Example in into the Roesser LSS model. By performing the l -scaling for the above realization with a diagonal coordinate matrix T o diag{.00000,.00000, , ,.00000,.00000, , } with truncation (0, 0) (i, j) (00, 00) in (A.) and (A.). 4.3 Application of the Lagrange method Choosing P (0) P (0) P (0) 4 I 8 in (9) as an initial estimate and a tolerance ε 0 8 in (3) as well as in the bisection method, it too the Lagrange-based algorithm 04 iterations to converge to the solution P opt or equivalently, T opt The minimized frequency-weighted l -sensitivity measure in (3) corresponding to the above solution was found to be J(P opt,λ,λ 4 ) with λ and λ The profile of the frequency-weighted l -sensitivity measure J(P,λ,λ 4 ) and the profiles of the Lagrange multipliers λ and λ 4 for the first 04 iterations are shown in Figs. 8 and 9, respectively.
9 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ J ( P () (+) (+), λ, λ 4 ) The minimized frequency-weighted l -sensitivity in (36) was found to be J o ( ˆT opt ) which is identical to the minimum value of the frequencyweighted l -sensitivity measure, obtained by the Lagrangebased method. The profile of the l -sensitivity measure J o ( ˆT ) during the first 54 iterations is shown in Fig Fig. 8. Profile of J(P,λ,λ 4 ) during the first 04 iterations. J o (x ) , λ 4 () λ () Fig. 9. Profiles of λ and λ 4 during the first 04 iterations. 4.4 Application of the quasi-newton method By choosing ˆT ˆT ˆT 4 I 8 in (37) as an initial estimate and a tolerance ε 0 8 in (38), the quasi-newton algorithm too 54 iterations to converge to ˆT opt or equivalently, T opt λ λ Fig. 0. Profile of J o ( ˆT ) during the first 54 iterations. To compare the technique reported in with the proposed ones, the optimal realization derived in Example of was imbedded in the Roesser LSS model by using (6). Then the frequency-weighted l -sensitivity of the resulting imbedded model was computed from (7) as S It is observed that this value is.9 times greater than the minimized frequency-weighted l -sensitivity obtained by the proposed techniques. Moreover, the proposed Lagrangebased method attains considerably faster convergence than that reported in, which required 000 iterations for its convergence. V. CONCLUSION We have investigated the problem of minimizing the frequency-weighted l -sensitivity subject to l -scaling constraints for -D state-space digital filters described by the Roesser LSS model. It has been shown that the FM second LSS model can be imbedded in the Roesser LSS model as a special case. Two iterative methods have been developed to solve the problem at hand. The first iterative method is based on the introduction of a Lagrange function and maes use of an efficient bisection method. In our simulation studies, the bisection method was found to be considerably faster than the iterative method proposed in. The second iterative method relies on the conversion of the constrained optimization problem into an unconstrained optimization formulation and utilizes an efficient quasi-newton algorithm. The optimal state-space realization with minimum frequency-weighted l - sensitivity and no overflow has then been constructed by applying an appropriate coordinate transformation. Our computer simulation results have demonstrated the validity and effectiveness of the proposed techniques.
10 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ x h (i +,j+) A A x h (i, j +) x v (i +,j+) 0 0 x v (i, j +) 0 0 x h (i +,j) + A 3 A 4 x v (i +,j) b 0 + u(i, j +)+ u(i +,j) 0 b y(i, j) x h (i, j) c c + du(i, j) x v (i, j) (39) APPENDIX I COMPUTATIONS OF GRAMIANS The matrices K C, W B,andM A can be computed using A (,0) K C W B M A i0 j0 i0 j0 A A 0 0 i0 j0 f C (i, j)f C T (i, j) gb T (i, j)g B(i, j) H A T (i, j)h A(i, j), A (0,) 0 0 A 3 A 4 (A.) A (0,0) I m+n, A ( i,j ) 0 (i ), A (i, j ) 0 (j ) A (i,j) A (,0) A (i,j) + A (0,) A (i,j ) A (i,j) A (,0) + A (i,j ) A (0,), (i, j) > (0, 0) + A (i,j ) 0 b f (i, j) A (i,j) b 0 + ca (i,j ) I n g (i, j) ca (i,j) Im H(i, j) f (, r)g (i, j r) f C (i, j) g B (i, j) H A (i, j) (0,0) (,r)<(i,j) (0,0) (,r)<(i,j) (0,0) (,r)<(i,j) (0,0) (,r)<(i,j) w C (, r)f (i, j r) w B (, r)g (i, j r) w A (, r)h (i, j r) with partial ordering for integer pairs (i, j) 6, p., and w A (, r), w B (, r), andw C (, r) denoting the unit-sample responses of frequency-weighting functions W A (z,z ), W B (z,z ),andw C (z,z ), respectively. The local controllability Gramian K and the other Gramians M A (P ), N A (P ) and ˆM A ( ˆP ) can be computed using K f (i, j)f T (i, j) M A (P ) N A (P ) ˆM A ( ˆP ) i0 j0 i0 j0 i0 j0 i0 j0 HA T (i, j) P H A (i, j) H A (i, j) PH T A(i, j) Ĥ T A(i, j) ˆP Ĥ A (i, j) Ĥ A (i, j) (K K 4 ) HA (i, j)(k K 4 ). APPENDIX II GRADIENT EVALUATION OF J o (x) J o ( ˆT ) t ij (A.) lim Δ 0 J o ( ˆT ij ) J o ( ˆT ) Δ β β +β 3 β 4 (A.3) ˆT ij is the matrix obtained from ˆT ˆT ˆT 4 with a perturbed (i, j)th component, which is given by 35, p. 655 with ˆT ij ˆT + Δ ˆTg ij e T ˆT j Δe T ˆTg, ˆT ij ˆT Δg ij e T j j ij { } tj g ij / t ij t j t j 3 (t ijt j t j e i ) β e T j ˆM A ( ˆT ) ˆTg ij β e T ˆT T j Ĥ A (p, q) ˆT T ˆT Ĥ T A (p, q) g ij p0 q0 β 3 e T j ˆT Ŵ B ˆT T ˆTgij, β 4 e T j ˆT T ˆKc g ij {t, t,, t m+n } {t, t,, t m } {t 4, t 4,, t 4n }. REFERENCES L. Thiele, Design of sensitivity and round-off noise optimal state-space discrete systems, Int. J. Circuit Theory Appl., vol., pp , Jan V. Tavsanoglu and L. Thiele, Optimal design of state-space digital filters by simultaneous minimization of sensitivity and roundoff noise, IEEE Trans. Circuits Syst., vol. CAS-3, pp , Oct L. Thiele, On the sensitivity of linear state-space systems, IEEE Trans. Circuits Syst., vol. CAS-33, pp , May M. Iwatsui, M. Kawamata and T. Higuchi, Statistical sensitivity and minimum sensitivity structures with fewer coefficients in discrete time linear systems, IEEE Trans. Circuits Syst., vol. 37, pp. 7-80, Jan G. Li, B. D. O. Anderson, M. Gevers and J. E. Perins, Optimal FWL design of state-space digital systems with weighted sensitivity minimization and sparseness consideration, IEEE Trans. Circuits Syst. I, vol. 39, pp , May W.-Y. Yan and J. B. Moore, On L -sensitivity minimization of linear state-space systems, IEEE Trans. Circuits Syst. I, vol. 39, pp , Aug. 99.
11 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. YY, ZZZ G. Li and M. Gevers, Optimal synthetic FWL design of state-space digital filters, in Proc. 99 IEEE Int. Conf. Acoust., Speech, Signal Processing, vol. 4, pp M. Gevers and G. Li, Parameterizations in Control, Estimation and Filtering Problems: Accuracy Aspects, Springer-Verlag, C. Xiao, Improved L -sensitivity for state-space digital system, IEEE Trans. Signal Processing, vol. 45, pp , Apr T. Hinamoto, S. Yooyama, T. Inoue, W. Zeng and W.-S. Lu, Analysis and minimization of L -sensitivity for linear systems and twodimensional state-space filters using general controllability and observability Gramians, IEEE Trans. Circuits Syst. I, vol. 49, pp , Sept. 00. S. Yamai, M. Abe and M. Kawamata, A closed form solution to L - sensitivity minimization of second-order state-space digital filters, in Proc. 006 IEEE Int. Symp. Circuits Syst., pp M. Kawamata, T. Lin and T. Higuchi, Minimization of sensitivity of - D state-space digital filters and its relation to -D balanced realizations, in Proc. 987 IEEE Int. Symp. Circuits Syst., pp T. Hinamoto, T. Hamanaa and S. Maeawa, Synthesis of -D state-space digital filters with low sensitivity based on the Fornasini- Marchesini model, IEEE Trans. Acoust., Speech, Signal Processing, vol. ASSP-38, pp , Sept T. Hinamoto, T. Taao and M. Muneyasu, Synthesis of -D separabledenominator digital filters with low sensitivity, J. Franlin Institute, vol. 39, pp , T. Hinamoto and T. Taao, Synthesis of -D state-space filter structures with low frequency-weighted sensitivity, IEEE Trans. Circuits Syst. II, vol. 39, pp , Sept T. Hinamoto and T. Taao, Minimization of frequency-weighting sensitivity in -D systems based on the Fornasini-Marchesini second model, in 99 IEEE Int. Conf. Acoust., Speech, Signal Processing, pp T. Hinamoto, Y. Zempo, Y. Nishino and W.-S. Lu, An analytical approach for the synthesis of two-dimensional state-space filter structures with minimum weighted sensitivity, IEEE Trans. Circuits Syst. I, vol. 46, pp. 7-83, Oct G. Li, On frequency weighted minimal L sensitivity of -D systems using Fornasini-Marchesini LSS model, IEEE Trans. Circuits Syst. I, vol. 44, pp , July G. Li, Two-dimensional system optimal realizations with L -sensitivity minimization, IEEE Trans. Signal Processing, vol. 46, pp , Mar T. Hinamoto and Y. Sugie, L -sensitivity analysis and minimization of -D separable-denominator state-space digital filters, IEEE Trans. Signal Processing, vol. 50, pp , Dec. 00. T. Hinamoto, H. Ohnishi and W.-S. Lu, Minimization of L -sensitivity for state-space digital filters subject to L -dynamic-range scaling constraints, IEEE Trans. Circuits Syst.-II, vol. 5, pp , Oct T. Hinamoto, K. Iwata and W.-S. Lu, L -sensitivity Minimization of one- and two-dimensional state-space digital filters subject to L -scaling constraints, IEEE Trans. Signal Processing, vol. 54, pp , May S. Yamai, M. Abe and M. Kawamata, A novel approach to L - sensitivity minimization of digital filters subject to L -scaling constraints, in Proc. 006 IEEE Int. Symp. Circuits Syst., pp C. T. Mullis and R. A. Roberts, Synthesis of minimum roundoff noise fixed-point digital filters, IEEE Trans. Circuits Syst., vol. CAS-3, pp , Sept S. Y. Hwang, Minimum uncorrelated unit noise in state-space digital filtering, IEEE Trans. Acoust., Speech, Signal Processing, vol. ASSP- 5, pp. 73-8, Aug R. P. Roesser, A discrete state-space model for linear image processing, IEEE Trans. Automat. Contr., vol. AC-0, pp. -0, Feb E. Fornasini and G. Marchesini, Doubly-indexed dynamical systems: State-space models and structural properties, Math Syst. Theory, vol., pp. 59-7, T. Hinamoto, A novel local state-space model for -D digital filters and its properties, in Proc. 00 IEEE Int. Symp. Circuits Syst., vol., pp S. Kung, B. C. Levy, M. Morf and T. Kailath, New results in -D systems theory, Part II: -D state-space model Realization and the notions of controllability, observability, and minimality, Proc. IEEE, vol. 65, pp , June W.-S. Lu and A. Antoniou, Synthesis of -D state-space fixed-point digital filter structures with minimum roundoff noise, IEEE Trans. Circuits Syst., vol. CAS-33, pp , Oct L. L. Scharf, Statistical Signal Processing, Reading, MA: Addison- Wesley, H. Togawa, Handboo of Numerical Methods, Saiensu-sha, Toyo, R. Fletcher, Practical Methods of Optimization, nd ed., Wiley, New Yor, S. A. H. Aly and M. M. Fahmy, Spatial-domain design of twodimensional recursive digital filters, IEEE Trans. Circuits Syst., vol. CAS-7, pp , Oct T. Kailath, Linear System, Englewood Cliffs, N.J.: Prentice Hall, 980.
Two-Dimensional State-Space Digital Filters with Minimum Frequency-Weighted L 2 -Sensitivity
wo-dimensional State-Space Digital Filters with Minimum Frequency-Weighted L -Sensitivity under L -Scaling Constraints. Hinamoto,. Oumi, O. I. Omoifo W.-S. Lu Graduate School of Engineering Electrical
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 informationWe are IntechOpen, the first native scientific publisher of Open Access books. International authors and editors. Our authors are among the TOP 1%
We are IntechOpen, the first native scientific publisher of Open Access books 3,35 18, 1.7 M Open access books available International authors and editors Downloads Our authors are among the 151 Countries
More informationDURING THE last two decades, many authors have
614 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: FUNDAMENTAL THEORY AND APPLICATIONS, VOL 44, NO 7, JULY 1997 Stability the Lyapunov Equation for -Dimensional Digital Systems Chengshan Xiao, David J Hill,
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 informationKNOWN approaches for improving the performance of
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 58, NO. 8, AUGUST 2011 537 Robust Quasi-Newton Adaptive Filtering Algorithms Md. Zulfiquar Ali Bhotto, Student Member, IEEE, and Andreas
More informationON THE REALIZATION OF 2D LATTICE-LADDER DISCRETE FILTERS
Journal of Circuits Systems and Computers Vol. 3 No. 5 (2004) 5 c World Scientific Publishing Company ON THE REALIZATION OF 2D LATTICE-LADDER DISCRETE FILTERS GEORGE E. ANTONIOU Department of Computer
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 informationDIGITAL filters capable of changing their frequency response
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL 46, NO 8, AUGUST 1999 1035 An Improved Weighted Least-Squares Design Variable Fractional Delay FIR Filters Wu-Sheng
More informationSECTION: CONTINUOUS OPTIMISATION LECTURE 4: QUASI-NEWTON METHODS
SECTION: CONTINUOUS OPTIMISATION LECTURE 4: QUASI-NEWTON METHODS HONOUR SCHOOL OF MATHEMATICS, OXFORD UNIVERSITY HILARY TERM 2005, DR RAPHAEL HAUSER 1. The Quasi-Newton Idea. In this lecture we will discuss
More informationClosed-Form Design of Maximally Flat IIR Half-Band Filters
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 49, NO. 6, JUNE 2002 409 Closed-Form Design of Maximally Flat IIR Half-B Filters Xi Zhang, Senior Member, IEEE,
More informationMATH 4211/6211 Optimization Quasi-Newton Method
MATH 4211/6211 Optimization Quasi-Newton Method Xiaojing Ye Department of Mathematics & Statistics Georgia State University Xiaojing Ye, Math & Stat, Georgia State University 0 Quasi-Newton Method Motivation:
More informationLMI based Stability criteria for 2-D PSV system described by FM-2 Model
Vol-4 Issue-018 LMI based Stability criteria for -D PSV system described by FM- Model Prashant K Shah Department of Electronics Engineering SVNIT, pks@eced.svnit.ac.in Abstract Stability analysis is the
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 informationComputing the pth Roots of a Matrix. with Repeated Eigenvalues
Applied Mathematical Sciences, Vol. 5, 2011, no. 53, 2645-2661 Computing the pth Roots of a Matrix with Repeated Eigenvalues Amir Sadeghi 1, Ahmad Izani Md. Ismail and Azhana Ahmad School of Mathematical
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 informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon s method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
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 informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationQUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY
QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY D.D. Olesky 1 Department of Computer Science University of Victoria Victoria, B.C. V8W 3P6 Michael Tsatsomeros Department of Mathematics
More informationECE133A Applied Numerical Computing Additional Lecture Notes
Winter Quarter 2018 ECE133A Applied Numerical Computing Additional Lecture Notes L. Vandenberghe ii Contents 1 LU factorization 1 1.1 Definition................................. 1 1.2 Nonsingular sets
More informationSensitivity Analysis of Coupled Resonator Filters
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 47, NO. 10, OCTOBER 2000 1017 Sensitivity Analysis of Coupled Resonator Filters Smain Amari, Member, IEEE Abstract
More informationEfficient and Accurate Rectangular Window Subspace Tracking
Efficient and Accurate Rectangular Window Subspace Tracking Timothy M. Toolan and Donald W. Tufts Dept. of Electrical Engineering, University of Rhode Island, Kingston, RI 88 USA toolan@ele.uri.edu, tufts@ele.uri.edu
More informationDiagonal matrix solutions of a discrete-time Lyapunov inequality
Diagonal matrix solutions of a discrete-time Lyapunov inequality Harald K. Wimmer Mathematisches Institut Universität Würzburg D-97074 Würzburg, Germany February 3, 1997 Abstract Diagonal solutions of
More informationTime-Varying Systems and Computations Lecture 3
Time-Varying Systems and Computations Lecture 3 Klaus Diepold November 202 Linear Time-Varying Systems State-Space System Model We aim to derive the matrix containing the time-varying impulse responses
More informationImproved Quasi-Newton Adaptive-Filtering Algorithm Md Zulfiquar Ali Bhotto, Student Member, IEEE, and Andreas Antoniou, Life Fellow, IEEE
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: REGULAR PAPERS, VOL. 57, NO. 8, AUGUST 2010 2109 Improved Quasi-Newton Adaptive-Filtering Algorithm Md Zulfiquar Ali Bhotto, Student Member, IEEE, and Andreas
More informationA NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1
A NOVEL OPTIMAL PROBABILITY DENSITY FUNCTION TRACKING FILTER DESIGN 1 Jinglin Zhou Hong Wang, Donghua Zhou Department of Automation, Tsinghua University, Beijing 100084, P. R. China Control Systems Centre,
More informationControl for stability and Positivity of 2-D linear discrete-time systems
Manuscript received Nov. 2, 27; revised Dec. 2, 27 Control for stability and Positivity of 2-D linear discrete-time systems MOHAMMED ALFIDI and ABDELAZIZ HMAMED LESSI, Département de Physique Faculté des
More informationOn Solving Large Algebraic. Riccati Matrix Equations
International Mathematical Forum, 5, 2010, no. 33, 1637-1644 On Solving Large Algebraic Riccati Matrix Equations Amer Kaabi Department of Basic Science Khoramshahr Marine Science and Technology University
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
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 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 informationBobby Hunt, Mariappan S. Nadar, Paul Keller, Eric VonColln, and Anupam Goyal III. ASSOCIATIVE RECALL BY A POLYNOMIAL MAPPING
Synthesis of a Nonrecurrent Associative Memory Model Based on a Nonlinear Transformation in the Spectral Domain p. 1 Bobby Hunt, Mariappan S. Nadar, Paul Keller, Eric VonColln, Anupam Goyal Abstract -
More informationADAPTIVE signal processing algorithms (ASPA s) are
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 46, NO. 12, DECEMBER 1998 3315 Locally Optimum Adaptive Signal Processing Algorithms George V. Moustakides Abstract We propose a new analytic method for comparing
More informationDeconvolution Filtering of 2-D Digital Systems
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 50, NO. 9, SEPTEMBER 2002 2319 H Deconvolution Filtering of 2-D Digital Systems Lihua Xie, Senior Member, IEEE, Chunling Du, Cishen Zhang, and Yeng Chai Soh
More information/$ IEEE
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 55, NO. 9, SEPTEMBER 2008 937 Analytical Stability Condition of the Latency Insertion Method for Nonuniform GLC Circuits Subramanian N.
More informationJim Lambers MAT 610 Summer Session Lecture 2 Notes
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 2 Notes These notes correspond to Sections 2.2-2.4 in the text. Vector Norms Given vectors x and y of length one, which are simply scalars x and y, the
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 informationUTILIZING PRIOR KNOWLEDGE IN ROBUST OPTIMAL EXPERIMENT DESIGN. EE & CS, The University of Newcastle, Australia EE, Technion, Israel.
UTILIZING PRIOR KNOWLEDGE IN ROBUST OPTIMAL EXPERIMENT DESIGN Graham C. Goodwin James S. Welsh Arie Feuer Milan Depich EE & CS, The University of Newcastle, Australia 38. EE, Technion, Israel. Abstract:
More informationCOVARIANCE CALCULATION FOR FLOATING POINT STATE SPACE REALIZATIONS
opyright 2002 IFA 15th riennial World ongress, Barcelona, Spain OVARIANE ALULAION FOR FLOAING POIN SAE SPAE REALIZAIONS Sangho Ko,1 Robert R Bitmead,1 Department of Mechanical and Aerospace Engineering,
More informationOptimum Sampling Vectors for Wiener Filter Noise Reduction
58 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 50, NO. 1, JANUARY 2002 Optimum Sampling Vectors for Wiener Filter Noise Reduction Yukihiko Yamashita, Member, IEEE Absact Sampling is a very important and
More informationNP-hardness of the stable matrix in unit interval family problem in discrete time
Systems & Control Letters 42 21 261 265 www.elsevier.com/locate/sysconle NP-hardness of the stable matrix in unit interval family problem in discrete time Alejandra Mercado, K.J. Ray Liu Electrical and
More informationSplit Algorithms and ZW-Factorization for Toeplitz and Toeplitz-plus-Hankel Matrices
Split Algorithms and ZW-Factorization for Toeplitz and Toeplitz-plus-Hanel Matrices Georg Heinig Dept of Math& CompSci Kuwait University POBox 5969 Safat 13060, Kuwait Karla Rost Faultät für Mathemati
More information3 Finite Wordlength Effects
Bomar, B.W. Finite Wordlength Effects Digital Signal Processing Handbook Ed. Vijay K. Madisetti and Douglas B. Williams Boca Raton: CRC Press LLC, 1999 c 1999byCRCPressLLC 3 Finite Wordlength Effects Bruce
More information1 Introduction 198; Dugard et al, 198; Dugard et al, 198) A delay matrix in such a lower triangular form is called an interactor matrix, and almost co
Multivariable Receding-Horizon Predictive Control for Adaptive Applications Tae-Woong Yoon and C M Chow y Department of Electrical Engineering, Korea University 1, -a, Anam-dong, Sungbu-u, Seoul 1-1, Korea
More informationMatrix Arithmetic. a 11 a. A + B = + a m1 a mn. + b. a 11 + b 11 a 1n + b 1n = a m1. b m1 b mn. and scalar multiplication for matrices via.
Matrix Arithmetic There is an arithmetic for matrices that can be viewed as extending the arithmetic we have developed for vectors to the more general setting of rectangular arrays: if A and B are m n
More informationON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION
Annales Univ. Sci. Budapest., Sect. Comp. 33 (2010) 273-284 ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION L. László (Budapest, Hungary) Dedicated to Professor Ferenc Schipp on his 70th
More informationx i2 j i2 ε + ε i2 ε ji1 ε j i1
T-S Fuzzy Model with Linear Rule Consequence and PDC Controller: A Universal Framewor for Nonlinear Control Systems Hua O. Wang Jing Li David Niemann Laboratory for Intelligent and Nonlinear Control (LINC)
More informationNonlinear Discrete-Time Observer Design with Linearizable Error Dynamics
622 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 4, APRIL 2003 Nonlinear Discrete-Time Observer Design with Linearizable Error Dynamics MingQing Xiao, Nikolaos Kazantzis, Costas Kravaris, Arthur
More information17 Solution of Nonlinear Systems
17 Solution of Nonlinear Systems We now discuss the solution of systems of nonlinear equations. An important ingredient will be the multivariate Taylor theorem. Theorem 17.1 Let D = {x 1, x 2,..., x m
More informationTowards Global Design of Orthogonal Filter Banks and Wavelets
Towards Global Design of Orthogonal Filter Banks and Wavelets Jie Yan and Wu-Sheng Lu Department of Electrical and Computer Engineering University of Victoria Victoria, BC, Canada V8W 3P6 jyan@ece.uvic.ca,
More informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
More informationMaximum Likelihood Diffusive Source Localization Based on Binary Observations
Maximum Lielihood Diffusive Source Localization Based on Binary Observations Yoav Levinboo and an F. Wong Wireless Information Networing Group, University of Florida Gainesville, Florida 32611-6130, USA
More informationCovariance Calculation for Floating Point. State Space Realizations
Covariance Calculation for Floating Point 1 State Space Realizations Sangho Ko, Robert R. Bitmead Abstract This paper provides a new method for analyzing floating point roundoff error for digital filters
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationFundamentals of Engineering Analysis (650163)
Philadelphia University Faculty of Engineering Communications and Electronics Engineering Fundamentals of Engineering Analysis (6563) Part Dr. Omar R Daoud Matrices: Introduction DEFINITION A matrix is
More informationMinimum Repair Bandwidth for Exact Regeneration in Distributed Storage
1 Minimum Repair andwidth for Exact Regeneration in Distributed Storage Vivec R Cadambe, Syed A Jafar, Hamed Malei Electrical Engineering and Computer Science University of California Irvine, Irvine, California,
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationThe Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems
668 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL 49, NO 10, OCTOBER 2002 The Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems Lei Zhang, Quan
More informationIN neural-network training, the most well-known online
IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 10, NO. 1, JANUARY 1999 161 On the Kalman Filtering Method in Neural-Network Training and Pruning John Sum, Chi-sing Leung, Gilbert H. Young, and Wing-kay Kan
More informationComputing running DCTs and DSTs based on their second-order shift properties
University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering Information Sciences 000 Computing running DCTs DSTs based on their second-order shift properties
More informationTransients in Reconfigurable Signal Processing Channels
936 IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT, VOL. 50, NO. 4, AUGUST 2001 Transients in Reconfigurable Signal Processing Channels Tamás Kovácsházy, Member, IEEE, Gábor Péceli, Fellow, IEEE,
More informationEE-597 Notes Quantization
EE-597 Notes Quantization Phil Schniter June, 4 Quantization Given a continuous-time and continuous-amplitude signal (t, processing and storage by modern digital hardware requires discretization in both
More informationAdaptive State Feedback Nash Strategies for Linear Quadratic Discrete-Time Games
Adaptive State Feedbac Nash Strategies for Linear Quadratic Discrete-Time Games Dan Shen and Jose B. Cruz, Jr. Intelligent Automation Inc., Rocville, MD 2858 USA (email: dshen@i-a-i.com). The Ohio State
More informationMinimal positive realizations of transfer functions with nonnegative multiple poles
1 Minimal positive realizations of transfer functions with nonnegative multiple poles Béla Nagy Máté Matolcsi Béla Nagy is Professor at the Mathematics Department of Technical University, Budapest, e-mail:
More informationJ. Liang School of Automation & Information Engineering Xi an University of Technology, China
Progress In Electromagnetics Research C, Vol. 18, 245 255, 211 A NOVEL DIAGONAL LOADING METHOD FOR ROBUST ADAPTIVE BEAMFORMING W. Wang and R. Wu Tianjin Key Lab for Advanced Signal Processing Civil Aviation
More informationPAPER Closed Form Expressions of Balanced Realizations of Second-Order Analog Filters
565 PAPER Closed Form Expressions of Balanced Realizations of Second-Order Analog Filters Shunsuke YAMAKI a), Memer, Masahide ABE ), Senior Memer, and Masayuki KAWAMATA c), Fellow SUMMARY This paper derives
More informationQuasi-Newton methods for minimization
Quasi-Newton methods for minimization Lectures for PHD course on Numerical optimization Enrico Bertolazzi DIMS Universitá di Trento November 21 December 14, 2011 Quasi-Newton methods for minimization 1
More informationarxiv: v1 [cs.sy] 2 Apr 2019
On the Existence of a Fixed Spectrum for a Multi-channel Linear System: A Matroid Theory Approach F Liu 1 and A S Morse 1 arxiv:190401499v1 [cssy] 2 Apr 2019 Abstract Conditions for the existence of a
More informationFREQUENCY-WEIGHTED MODEL REDUCTION METHOD WITH ERROR BOUNDS FOR 2-D SEPARABLE DENOMINATOR DISCRETE SYSTEMS
INERNAIONAL JOURNAL OF INFORMAION AND SYSEMS SCIENCES Volume 1, Number 2, Pages 105 119 c 2005 Institute for Scientific Computing and Information FREQUENCY-WEIGHED MODEL REDUCION MEHOD WIH ERROR BOUNDS
More informationNON-LINEAR CONTROL OF OUTPUT PROBABILITY DENSITY FUNCTION FOR LINEAR ARMAX SYSTEMS
Control 4, University of Bath, UK, September 4 ID-83 NON-LINEAR CONTROL OF OUTPUT PROBABILITY DENSITY FUNCTION FOR LINEAR ARMAX SYSTEMS H. Yue, H. Wang Control Systems Centre, University of Manchester
More informationBisimilar Finite Abstractions of Interconnected Systems
Bisimilar Finite Abstractions of Interconnected Systems Yuichi Tazaki and Jun-ichi Imura Tokyo Institute of Technology, Ōokayama 2-12-1, Meguro, Tokyo, Japan {tazaki,imura}@cyb.mei.titech.ac.jp http://www.cyb.mei.titech.ac.jp
More informationTHE LARGE number of components in a power system
IEEE TRANSACTIONS ON POWER SYSTEMS, VOL. 19, NO. 2, MAY 2004 857 Order Reduction of the Dynamic Model of a Linear Weakly Periodic System Part I: General Methodology Abner Ramirez, Member, IEEE, Adam Semlyen,
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,
More informationStructured Condition Numbers of Symmetric Algebraic Riccati Equations
Proceedings of the 2 nd International Conference of Control Dynamic Systems and Robotics Ottawa Ontario Canada May 7-8 2015 Paper No. 183 Structured Condition Numbers of Symmetric Algebraic Riccati Equations
More informationMATHEMATICAL ENGINEERING TECHNICAL REPORTS
MATHEMATICAL ENGINEERING TECHNICAL REPORTS Combinatorial Relaxation Algorithm for the Entire Sequence of the Maximum Degree of Minors in Mixed Polynomial Matrices Shun SATO (Communicated by Taayasu MATSUO)
More informationLapped Unimodular Transform and Its Factorization
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL 50, NO 11, NOVEMBER 2002 2695 Lapped Unimodular Transform and Its Factorization See-May Phoong, Member, IEEE, and Yuan-Pei Lin, Member, IEEE Abstract Two types
More informationH 2 -optimal model reduction of MIMO systems
H 2 -optimal model reduction of MIMO systems P. Van Dooren K. A. Gallivan P.-A. Absil Abstract We consider the problem of approximating a p m rational transfer function Hs of high degree by another p m
More informationThe Steepest Descent Algorithm for Unconstrained Optimization
The Steepest Descent Algorithm for Unconstrained Optimization Robert M. Freund February, 2014 c 2014 Massachusetts Institute of Technology. All rights reserved. 1 1 Steepest Descent Algorithm The problem
More informationAlgebraic Algorithm for 2D Stability Test Based on a Lyapunov Equation. Abstract
Algebraic Algorithm for 2D Stability Test Based on a Lyapunov Equation Minoru Yamada Li Xu Osami Saito Abstract Some improvements have been proposed for the algorithm of Agathoklis such that 2D stability
More informationA STATE-SPACE APPROACH FOR THE ANALYSIS OF WAVE AND DIFFUSION FIELDS
ICASSP 2015 A STATE-SPACE APPROACH FOR THE ANALYSIS OF WAVE AND DIFFUSION FIELDS Stefano Maranò Donat Fäh Hans-Andrea Loeliger ETH Zurich, Swiss Seismological Service, 8092 Zürich ETH Zurich, Dept. Information
More informationTaylor series based nite dierence approximations of higher-degree derivatives
Journal of Computational and Applied Mathematics 54 (3) 5 4 www.elsevier.com/locate/cam Taylor series based nite dierence approximations of higher-degree derivatives Ishtiaq Rasool Khan a;b;, Ryoji Ohba
More informationON THE QR ITERATIONS OF REAL MATRICES
Unspecified Journal Volume, Number, Pages S????-????(XX- ON THE QR ITERATIONS OF REAL MATRICES HUAJUN HUANG AND TIN-YAU TAM Abstract. We answer a question of D. Serre on the QR iterations of a real matrix
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 School of Mathematical Sciences Projet APICS Universiteit
More informationDesign of Projection Matrix for Compressive Sensing by Nonsmooth Optimization
Design of Proection Matrix for Compressive Sensing by Nonsmooth Optimization W.-S. Lu T. Hinamoto Dept. of Electrical & Computer Engineering Graduate School of Engineering University of Victoria Hiroshima
More informationMathematical Constraint on Functions with Continuous Second Partial Derivatives
1 Mathematical Constraint on Functions with Continuous Second Partial Derivatives J.D. Franson Physics Department, University of Maryland, Baltimore County, Baltimore, MD 15 Abstract A new integral identity
More informationLeast Squares Approximation
Chapter 6 Least Squares Approximation As we saw in Chapter 5 we can interpret radial basis function interpolation as a constrained optimization problem. We now take this point of view again, but start
More informationInverse of a Square Matrix. For an N N square matrix A, the inverse of A, 1
Inverse of a Square Matrix For an N N square matrix A, the inverse of A, 1 A, exists if and only if A is of full rank, i.e., if and only if no column of A is a linear combination 1 of the others. A is
More informationPerformance Comparison of Two Implementations of the Leaky. LMS Adaptive Filter. Scott C. Douglas. University of Utah. Salt Lake City, Utah 84112
Performance Comparison of Two Implementations of the Leaky LMS Adaptive Filter Scott C. Douglas Department of Electrical Engineering University of Utah Salt Lake City, Utah 8411 Abstract{ The leaky LMS
More informationAssignment 10. Arfken Show that Stirling s formula is an asymptotic expansion. The remainder term is. B 2n 2n(2n 1) x1 2n.
Assignment Arfken 5.. Show that Stirling s formula is an asymptotic expansion. The remainder term is R N (x nn+ for some N. The condition for an asymptotic series, lim x xn R N lim x nn+ B n n(n x n B
More informationNumerical Methods I Solving Square Linear Systems: GEM and LU factorization
Numerical Methods I Solving Square Linear Systems: GEM and LU factorization Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 18th,
More informationDesign of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear Continuous Systems
Systems Science and Applied Mathematics Vol. 1 No. 3 2016 pp. 29-37 http://www.aiscience.org/journal/ssam Design of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear
More informationReview of Basic Concepts in Linear Algebra
Review of Basic Concepts in Linear Algebra Grady B Wright Department of Mathematics Boise State University September 7, 2017 Math 565 Linear Algebra Review September 7, 2017 1 / 40 Numerical Linear Algebra
More informationFilter-Generating Systems
24 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 47, NO. 3, MARCH 2000 Filter-Generating Systems Saed Samadi, Akinori Nishihara, and Hiroshi Iwakura Abstract
More informationMAT 610: Numerical Linear Algebra. James V. Lambers
MAT 610: Numerical Linear Algebra James V Lambers January 16, 2017 2 Contents 1 Matrix Multiplication Problems 7 11 Introduction 7 111 Systems of Linear Equations 7 112 The Eigenvalue Problem 8 12 Basic
More informationAlgorithms for Constrained Optimization
1 / 42 Algorithms for Constrained Optimization ME598/494 Lecture Max Yi Ren Department of Mechanical Engineering, Arizona State University April 19, 2015 2 / 42 Outline 1. Convergence 2. Sequential quadratic
More informationAdaptive beamforming for uniform linear arrays with unknown mutual coupling. IEEE Antennas and Wireless Propagation Letters.
Title Adaptive beamforming for uniform linear arrays with unknown mutual coupling Author(s) Liao, B; Chan, SC Citation IEEE Antennas And Wireless Propagation Letters, 2012, v. 11, p. 464-467 Issued Date
More informationMathematical preliminaries and error analysis
Mathematical preliminaries and error analysis Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan September 12, 2015 Outline 1 Round-off errors and computer arithmetic
More informationLinear Algebra. Linear Equations and Matrices. Copyright 2005, W.R. Winfrey
Copyright 2005, W.R. Winfrey Topics Preliminaries Systems of Linear Equations Matrices Algebraic Properties of Matrix Operations Special Types of Matrices and Partitioned Matrices Matrix Transformations
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More information