Developing LSQR Method for Sylvester Quaternion Tensor Equations

Abstract

In this paper, we consider solving a least-squares problem to the generalized Sylvester quaternion tensor equation. From the properties of quaternions and the Kronecker product, a tensor form of the LSQR algorithm is proposed for solving this problem, the convergence analysis of which is then established. Numerical results are reported to illustrate the feasibility and validity of the proposed algorithm compared with the tensor form of the CGLS method, including when the algorithm is tested with some randomly generated data and color video restoration problems.

Keywords

Sylvester quaternion tensor equation / Tensor form / LSQR / Color video restoration / 15A69 / 15B33 / 65F10 / 94A08

References

[1] Benbow SJ. Solving generalized least-squares problems with LSQR. SIAM J. Matrix Anal. Appl., 1999, 21(1): 166-177

[2] Datta BN. Numerical Methods for Linear Control Systems, 2004, Oxford, Elsevier Science Press

[3] Fernandez, J.M., Schneeberger, W.A.: Quaternionic computing. arXiv:quant-ph/0307017v2 (2004)

[4] Fong DCL, Saunders M. LSMR: an iterative algorithm for sparse least-squares problems. SIAM J Sci. Comput., 2011, 33(5): 2950-2971

[5] Ginzberg P. Quaternion Matrices: Statistical Properties and Applications to Signal Processing and Wavelets, 2013, London, Imperial College

[6] Golub GH, Kahan W. Calculating a matrix’s singular values and pseudo-inverse. SIAM J. Numer. Anal., 1965, 2(2): 205-224

[7] Guide, M.E., Ichi, A.E., Beik, F.P.A.: Tensor GMRES, and Golub-Kahan bidiagonalization methods via the Einstein product with applications to image and video processing. arXiv:2005.07458 (2020)

[8] Hamilton WR. Elements of Quaternions, 1866, Harlow, Longmans, Green Company

[9] He ZH. The general solution to a system of coupled Sylvester-type quaternion tensor equations involving η\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\eta $$\end{document}-Hermicity. Bull. Iran. Math. Soc., 2019, 45: 1407-1430

[10] He ZH , Tian J, Zhao YF, Yu SW. Consistency and general solutions to some Sylvester-like quaternion matrix equations. Symmetry, 2022, 14(7): 1350

[11] He ZH, Wang XX, Zhao YF. Eigenvalues of quaternion tensors with applications to color video processing. J. Sci. Comput., 2023, 9411

[12] Hu JJ, Ke YF, Ma CF. Efficient iterative method for generalized Sylvester quaternion tensor equation. Comput. Appl. Math., 2023, 425237

[13] Huang Y, Jia ZX. Some results on the regularization of LSQR for large-scale discrete ill-posed problems. Sci. China Math., 2017, 60: 701-718

[14] Jia ZX. Some properties of LSQR for large sparse linear least squares problems. J. Syst. Sci. Complex., 2010, 23(4): 815-821

[15] Jiang TS, Zhao JL, Wei MS. A new technique of quaternion equality constrained least squares problem. J. Comput. Appl. Math., 2008, 216(2): 509-513

[16] Kolda TG, Bader BW. Tensor decompositions and applications. SIAM Rev., 2009, 51(3): 455-500

[17] Konstantinov M, Gu DW, Mehrmann VV. Perturbation Theory for Matrix Equations, 2003, Oxford, Elsevier Science Press

[18] Li SK, Huang TZ. LSQR iterative method for generalized coupled Sylvester matrix equations. Appl. Math. Model., 2012, 36(8): 3545-3554

[19] Paige CC, Saunders MA. LSQR: sparse linear equations and least squares problems. ACM Trans. Math. Soft., 1982, 8: 2

[20] Paige CC, Saunders MA. LSQR: an algorithm for sparse linear equations and sparse least squares. ACM Trans. Math. Soft., 1982, 8(1): 43-71

[21] Peng ZY. A matrix LSQR iterative method to solve matrix equation AXB=C\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$AXB=C$$\end{document}. Int. J. Comput. Math., 2010, 87(8): 1820-1830

[22] Qi, L.Q., Chen, H.B., Chen, Y.N.: Tensor Eigenvalues and Their Applications. Springer, Singapore (2018)

[23] Qin J, Wang QW. Solving a system of two-sided Sylvester-like quaternion tensor equations. Comput. Appl. Math., 2023, 42(5): 232

[24] Qin ZZ, Ming ZY, Zhang LP. Singular value decomposition of third order quaternion tensors. Appl. Math. Lett., 2022, 123 107597

[25] Reichel L, Sadok H, Zhang WH. Simple stopping criteria for the LSQR method applied to discrete ill-posed problems. Numer. Algorithms, 2020, 84: 1381-1395

[26] Rodman L. Topics in Quaternion Linear Algebra, 2014, Princeton, Princeton University Press

[27] Toutounian F, Karimi S. Global least squares method (Gl-LSQR) for solving general linear systems with several right-hand sides. Appl. Math. Comput., 2006, 178(2): 452-460

[28] Wang QW, Wang X, Zhang Y. A constraint system of coupled two-sided Sylvester-like quaternion tensor equations. Comput. Appl. Math., 2020, 39(4): 317

[29] Wang QW, Xu XJ, Duan XF. Least squares solution of the quaternion Sylvester tensor equation. Linear Multilinear Algebra, 2021, 691104-130

[30] Xie MY , Wang QW, He ZH, Saad MM. A system of Sylvester-type quaternion matrix equations with ten variables. Acta Math. Sin., 2022, 38(8): 1399-1420

[31] Xu XJ, Wang QW. Extending BiCG and BiCR methods to solve the Stein tensor equation. Comput. Math. Appl., 2019, 77(12): 3117-3127

[32] Yu SW, He ZH, Qin WL. Some systems of tensor equations under t-product and their applications. Filomat, 2021, 35(11): 3663-3677

[33] Zhang FZ. Quaternions and matrices of quaternions. Linear Algebra Appl., 1997, 251: 21-57

[34] Zhang XF, Ding W, Li T. Tensor form of GPBiCG algorithm for solving the generalized Sylvester quaternion tensor equations. J. Franklin Inst., 2023, 360(9): 5929-5946

[35] Zhang, X.F., Li, T., Ou, Y.G.: Iterative solutions of generalized Sylvester quaternion tensor equations. Linear and Multilinear Algebra (2023). https://doi.org/10.1080/03081087.2023.2176416

Funding

National Natural Science Foundation of China(11961018)

Natural Science Foundation of Hainan Province(122MS001)

RIGHTS & PERMISSIONS

Shanghai University