We first give the style spectral decomposition of a special skew circulant matrix and then get the style decomposition of arbitrary skew circulant matrix by making use of the Kronecker products between the elements of first row in skew circulant and the special skew circulant . Besides that, we obtain the singular value of skew circulant matrix as well. Finally, we deal with the optimal backward perturbation analysis for the linear system with skew circulant coefficient matrix on the base of its style spectral decomposition. 1. Introduction A skew circulant matrix with the first row is a square matrix of the form denoted by . Skew circulant matrices have important applications in various disciplines including image processing, signal processing, solving Toeplitz matrix problems, and preconditioner. The skew circulant matrices are considered as preconditioners for linear-multistep-formulae (LMF-) based ordinary differential equations (ODEs) codes; Hermitian and skew-Hermitian Toeplitz systems are considered in [1–4]. Lyness and S?revik [5] employed a skew circulant matrix to construct -dimensional lattice rules. Spectral decompositions of skew circulant and skew left circulant matrices are discussed in [6]. Akhondi and Toutounian [7] presented a new iteration method for the numerical solution of Hermitian positive definite Toeplitz systems of linear equations. Narasimha [8] believed that the linear convolution required in block filtering can be decomposed into a sum of skew-circulant convolutions and such convolutions can be realized efficiently with half-length complex transforms when the signals are real. Liu and Vaidyanathan [9] presented a new family of normal form state-space structures, the method used allows people to synthesize in normal form, most IIR transfer functions, and the state transition matrices involved are either circulant or skew circulant matrices. Vaidyanathan and Pal [10] examined a case where two arrays are generated by matrices that are adjugates of each other; in this case, it is possible to obtain a dense rectangular tiling of the 2 frequency plane from a pair of coarse 2 DFT filter banks; the special case where the adjugate pairs are generated by skew circulant matrices has some advantages, which are examined in detail. An additional convolution-multiplication property for the skew-circulant convolution operation , where is a skew-circulant matrix; besides, skew-circulant convolution is the underlying form of convolution in half of the 40 cases of symmetric convolution, and the convolution is an extension of a result Vernet's
References
[1]
D. Bertaccini and M. K. Ng, “Skew-circulant preconditioners for systems of LMF-based ODE codes,” Numerical Analysis and Its Applications, vol. 1988, pp. 93–101, 2001.
[2]
R. H. Chan and X.-Q. Jin, “Circulant and skew-circulant preconditioners for skew-hermitian type Toeplitz systems,” BIT Numerical Mathematics, vol. 31, no. 4, pp. 632–646, 1991.
[3]
R. H. Chan and K.-P. Ng, “Toeplitz preconditioners for Hermitian Toeplitz systems,” Linear Algebra and Its Applications, vol. 190, pp. 181–208, 1993.
[4]
T. Huclke, “Circulant and skew-circulant matrices for solving Toeplitz matrix problems,” SIAM Journal on Matrix Analysis and Applications, vol. 13, pp. 767–777, 1992.
[5]
J. N. Lyness and T. S?revik, “Four-dimensional lattice rules generated by skew-circulant matrices,” Mathematics of Computation, vol. 73, no. 245, pp. 279–295, 2004.
[6]
H. Karner, J. Schneid, and C. W. Ueberhuber, “Spectral decomposition of real circulant matrices,” Linear Algebra and Its Applications, vol. 367, pp. 301–311, 2003.
[7]
N. Akhondi and F. Toutounian, “Accelerated circulant and skew circulant splitting methods for Hermitian positive definite to eplitz systems,” Advances in Numerical Analysis, vol. 2012, Article ID 973407, 17 pages, 2012.
[8]
M. J. Narasimha, “Linear convolution using skew-cyclic convolutions,” IEEE Signal Processing Letters, vol. 14, no. 3, pp. 173–176, 2007.
[9]
V. C. Liu and P. P. Vaidyanathan, “Circulant and skew-circulant matrices as new normal-form realization of IIR digital filters,” IEEE transactions on circuits and systems, vol. 35, no. 6, pp. 625–635, 1988.
[10]
P. P. Vaidyanathan and P. Pal, “Adjugate pairs of sparse arrays for sampling two dimensional signals,” in Proceedings of the 36th IEEE International Conference on Acoustics, Speech, and Signal Processing, ICASSP 2011, pp. 3936–3939, May 2011.
[11]
J. L. Vernet, “Real signals fast Fourier transform. Storage capacity and step number reduction by means of an odd discrete Fourier transform,” Proceedings of the IEEE, vol. 59, no. 10, pp. 1531–1532, 1971.
[12]
T. M. Foltz and B. M. Welsh, “Symmetric convolution of asymmetric multidimensional sequences using discrete trigonometric transforms,” IEEE Transactions on Image Processing, vol. 8, no. 5, pp. 640–651, 1999.
[13]
X. G. Liu and X. X. Guo, “On optimal backward perturbation analysis for the linear system with block cyclic coefficient matrix,” Numerical Mathematics, vol. 12, no. 2, pp. 162–172, 2003.
[14]
J.-G. Sun and Z. Sun, “Optimal backward perturbation bounds for underdetermined systems,” SIAM Journal on Matrix Analysis and Applications, vol. 18, no. 2, pp. 393–402, 1997.
[15]
J. L. Rigal and J. Gaches, “On the compatibility of a given solution with the data of a linear system,” Journal of the ACM, vol. 14, pp. 543–548, 1967.
[16]
H.-J. Wittsack, A. M. Wohlschl?ger, E. K. Ritzl et al., “CT-perfusion imaging of the human brain: advanced deconvolution analysis using circulant singular value decomposition,” Computerized Medical Imaging and Graphics, vol. 32, no. 1, pp. 67–77, 2008.