In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied tho...In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.展开更多
We consider matrix integrable fifth-order mKdV equations via a kind of group reductions of the Ablowitz–Kaup–Newell–Segur matrix spectral problems. Based on properties of eigenvalue and adjoint eigenvalue problems,...We consider matrix integrable fifth-order mKdV equations via a kind of group reductions of the Ablowitz–Kaup–Newell–Segur matrix spectral problems. Based on properties of eigenvalue and adjoint eigenvalue problems, we solve the corresponding Riemann–Hilbert problems, where eigenvalues could equal adjoint eigenvalues, and construct their soliton solutions, when there are zero reflection coefficients. Illustrative examples of scalar and two-component integrable fifthorder mKdV equations are given.展开更多
This paper discusses the necessary and sufficient conditions for the existence of Hermite positive definite solutions of the quaternion matrix equation X<sup>m</sup>+ B*XB = C (m > 0) and its iterative ...This paper discusses the necessary and sufficient conditions for the existence of Hermite positive definite solutions of the quaternion matrix equation X<sup>m</sup>+ B*XB = C (m > 0) and its iterative solution method. According to the characteristics of the coefficient matrix, a corresponding algebraic equation system is ingeniously constructed, and by discussing the equation system’s solvability, the matrix equation’s existence interval is obtained. Based on the characteristics of the coefficient matrix, some necessary and sufficient conditions for the existence of Hermitian positive definite solutions of the matrix equation are derived. Then, the upper and lower bounds of the positive actual solutions are estimated by using matrix inequalities. Four iteration formats are constructed according to the given conditions and existence intervals, and their convergence is proven. The selection method for the initial matrix is also provided. Finally, using the complexification operator of quaternion matrices, an equivalent iteration on the complex field is established to solve the equation in the Matlab environment. Two numerical examples are used to test the effectiveness and feasibility of the given method. .展开更多
We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, ...We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.展开更多
The Geometrical Optics(GO)approach and the FAST Emissivity Model(FASTEM)are widely used to estimate the surface radiative components in atmospheric radiative transfer simulations,but their applications are limited in ...The Geometrical Optics(GO)approach and the FAST Emissivity Model(FASTEM)are widely used to estimate the surface radiative components in atmospheric radiative transfer simulations,but their applications are limited in specific conditions.In this study,a two-scale reflectivity model(TSRM)and a two-scale emissivity model(TSEM)are developed from the two-scale roughness theory.Unlike GO which only computes six non-zero elements in the reflectivity matrix,The TSRM includes 16 elements of Stokes reflectivity matrix which are important for improving radiative transfer simulation accuracy in a scattering atmosphere.It covers the frequency range from L-to W-bands.The dependences of all TSRM elements on zenith angle,wind speed,and frequency are derived and analyzed in details.For a set of downwelling radiances in microwave frequencies,the reflected upwelling brightness temperature(BTs)are calculated from both TSRM and GO and compared for analyzing their discrepancies.The TSRM not only includes the effects of GO but also accounts for the small-scale Bragg scattering effect in an order of several degrees in Kelvins in brightness temperature.Also,the third and fourth components of the Stokes vector can only be produced from the TSRM.For the emitted radiation,BT differences in vertical polarization between a TSEM and FASTEM are generally less than 5 K when the satellite zenith angle is less than 40°,whereas those for the horizontal component can be quite significant,greater than 20 K.展开更多
In this article, the generalized reflexive solution of matrix equations (AX = B, XC = D) is considered. With special properties of generalized reflexive matrices, the necessary and sufficient conditions for the solv...In this article, the generalized reflexive solution of matrix equations (AX = B, XC = D) is considered. With special properties of generalized reflexive matrices, the necessary and sufficient conditions for the solvability and the general expression of the solution are obtained. Moreover, the related optimal approximation problem to a given matrix over the solution set is solved.展开更多
Dehghan and Hajarian, [4], investigated the matrix equations A^TXB+B^TX^TA = C and A^TXB + B^TXA = C providing inequalities for the determinant of the solutions of these equations. In the same paper, the authors pre...Dehghan and Hajarian, [4], investigated the matrix equations A^TXB+B^TX^TA = C and A^TXB + B^TXA = C providing inequalities for the determinant of the solutions of these equations. In the same paper, the authors presented a lower bound for the product of the eigenvalues of the solutions to these matrix equations. Inspired by their work, we give some generalizations of Dehghan and Hajarian results. Using the theory of the numerical ranges, we present an inequality involving the trace of C when A, B, X are normal matrices satisfying A^T B = BA^T.展开更多
We derive necessary and sufficient conditions for the existence and an expression of the (anti)reflexive solution with respect to the nontrivial generalized reflection matrix P to the system of complex matrix equati...We derive necessary and sufficient conditions for the existence and an expression of the (anti)reflexive solution with respect to the nontrivial generalized reflection matrix P to the system of complex matrix equations AX = B and XC = D. The explicit solutions of the approximation problem min x∈Ф ||X - E||F was given, where E is a given complex matrix and Ф is the set of all reflexive (or antireflexive) solutions of the system mentioned above, and ||·|| is the Frobenius norm. Furthermore, it was pointed that some results in a recent paper are special cases of this paper.展开更多
Some new oscillation criteria are established for the second-order matrix differential system (r(t)Z'(t))' + p(t)Z'(t) + Q(t)F(Z'(t))G(Z(t)) = 0, t ≥ to 〉 0, are different from most known ...Some new oscillation criteria are established for the second-order matrix differential system (r(t)Z'(t))' + p(t)Z'(t) + Q(t)F(Z'(t))G(Z(t)) = 0, t ≥ to 〉 0, are different from most known ones in the sense that they are based on the information only on a sequence of subintervals of [to,∞), rather than on the whole half-line. The results weaken the condition of Q(t) and generalize some well-known results of Wong (1999) to nonlinear matrix differential equation.展开更多
In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermit...In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermite matrix polynomials,the orthogonality property and a Rodrigues' formula are given.展开更多
This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations in...This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations into a system of algebraic equations.Meanwhile,the error analysis is proven.Finally,the effectiveness of the approach is verified by two numerical examples.展开更多
In this paper, solutions to the generalized Sylvester matrix equations AX -XF = BY and MXN -X = TY with A, M ∈ R^n×n, B, T ∈ Rn×r, F, N ∈ R^p×p and the matrices N, F being in companion form, are est...In this paper, solutions to the generalized Sylvester matrix equations AX -XF = BY and MXN -X = TY with A, M ∈ R^n×n, B, T ∈ Rn×r, F, N ∈ R^p×p and the matrices N, F being in companion form, are established by a singular value decomposition of a matrix with dimensions n × (n + pr). The algorithm proposed in this paper for the euqation AX - XF = BY does not require the controllability of matrix pair (A, B) and the restriction that A, F do not have common eigenvalues. Since singular value decomposition is adopted, the algorithm is numerically stable and may provide great convenience to the computation of the solution to these equations, and can perform important functions in many design problems in control systems theory.展开更多
In this paper, the maximal and minimal ranks of the solution to a system of matrix equations over H, the real quaternion algebra, were derived. A previous known result could be regarded as a special case of the new re...In this paper, the maximal and minimal ranks of the solution to a system of matrix equations over H, the real quaternion algebra, were derived. A previous known result could be regarded as a special case of the new result.展开更多
A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. T...A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. This is the first time that the inverse of this remarkable matrix is determined directly and exactly. Thus, solving 1D Poisson equation becomes very accurate and extremely fast. This method is a very important tool for physics and engineering where the Poisson equation appears very often in the description of certain phenomena.展开更多
In this paper, two different methods are used to study the cyclic structure solution and the optimal approximation of the quaternion Stein equation AXB - X = F . Firstly, the matrix equation equivalent to the ta...In this paper, two different methods are used to study the cyclic structure solution and the optimal approximation of the quaternion Stein equation AXB - X = F . Firstly, the matrix equation equivalent to the target structure matrix is constructed by using the complex decomposition of the quaternion matrix, to obtain the necessary and sufficient conditions for the existence of the cyclic solution of the equation and the expression of the general solution. Secondly, the Stein equation is converted into the Sylvester equation by adding the necessary parameters, and the condition for the existence of a cyclic solution and the expression of the equation’s solution are then obtained by using the real decomposition of the quaternion matrix and the Kronecker product of the matrix. At the same time, under the condition that the solution set is non-empty, the optimal approximation solution to the given quaternion circulant matrix is obtained by using the property of Frobenius norm property. Numerical examples are given to verify the correctness of the theoretical results and the feasibility of the proposed method. .展开更多
In this paper, an improved gradient iterative (GI) algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared wi...In this paper, an improved gradient iterative (GI) algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.展开更多
Let P ∈ C^(n×n) be a Hermitian and {k + 1}-potent matrix, i.e., P^(k+1)= P = P~*,where(·)*~stands for the conjugate transpose of a matrix. A matrix X ∈ Cn×nis called{P, k + 1}-reflexive(anti-reflexive...Let P ∈ C^(n×n) be a Hermitian and {k + 1}-potent matrix, i.e., P^(k+1)= P = P~*,where(·)*~stands for the conjugate transpose of a matrix. A matrix X ∈ Cn×nis called{P, k + 1}-reflexive(anti-reflexive) if PXP = X(P XP =-X). The system of matrix equations AX = C, XB = D subject to {P, k + 1}-reflexive and anti-reflexive constraints are studied by converting into two simpler cases: k = 1 and k = 2, the least squares solution and the associated optimal approximation problem are also considered.展开更多
An interesting semi-analytic solution is given for the Helmholtz equation. This solution is obtained from a rigorous discussion of the regularity and the inversion of the tridiagonal symmetric matrix. Then, applicatio...An interesting semi-analytic solution is given for the Helmholtz equation. This solution is obtained from a rigorous discussion of the regularity and the inversion of the tridiagonal symmetric matrix. Then, applications are given, showing very good accuracy. This work provides also the analytical inverse of the skew-symmetric tridiagonal matrix.展开更多
We derive necessary and sufficient conditions for the existence of a Hermitian nonnegative-definite solution to the matrix equation AXB = C. Moreover, we derive a representation of a general Hermitian nonnegative-defi...We derive necessary and sufficient conditions for the existence of a Hermitian nonnegative-definite solution to the matrix equation AXB = C. Moreover, we derive a representation of a general Hermitian nonnegative-definite solution. We then apply our solution to two examples, including a comparison of our solution to a proposed solution by Zhang in [1] using an example problem given from [1]. Our solution demonstrates that the proposed general solution from Zhang in [1] is incorrect. We also give a second example in which we derive the general covariance structure so that two matrix quadratic forms are independent.展开更多
In this paper, we introduce a Hermite operational matrix collocation method for solving higher-order linear complex differential equations in rectangular or elliptic domains. We show that based on a linear algebra the...In this paper, we introduce a Hermite operational matrix collocation method for solving higher-order linear complex differential equations in rectangular or elliptic domains. We show that based on a linear algebra theorem, the use of different polynomials such as Hermite, Bessel and Taylor in polynomial collocation methods for solving differential equations leads to an equal solution, and the difference in the numerical results arises from the difference in the coefficient matrix of final linear systems of equations. Some numerical examples will also be given.展开更多
文摘In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.
基金supported in part by the National Natural Science Foundation of China (Grant Nos. 11975145, 11972291, and 51771083)the Ministry of Science and Technology of China (Grant No. G2021016032L)the Natural Science Foundation for Colleges and Universities in Jiangsu Province, China (Grant No. 17 KJB 110020)。
文摘We consider matrix integrable fifth-order mKdV equations via a kind of group reductions of the Ablowitz–Kaup–Newell–Segur matrix spectral problems. Based on properties of eigenvalue and adjoint eigenvalue problems, we solve the corresponding Riemann–Hilbert problems, where eigenvalues could equal adjoint eigenvalues, and construct their soliton solutions, when there are zero reflection coefficients. Illustrative examples of scalar and two-component integrable fifthorder mKdV equations are given.
文摘This paper discusses the necessary and sufficient conditions for the existence of Hermite positive definite solutions of the quaternion matrix equation X<sup>m</sup>+ B*XB = C (m > 0) and its iterative solution method. According to the characteristics of the coefficient matrix, a corresponding algebraic equation system is ingeniously constructed, and by discussing the equation system’s solvability, the matrix equation’s existence interval is obtained. Based on the characteristics of the coefficient matrix, some necessary and sufficient conditions for the existence of Hermitian positive definite solutions of the matrix equation are derived. Then, the upper and lower bounds of the positive actual solutions are estimated by using matrix inequalities. Four iteration formats are constructed according to the given conditions and existence intervals, and their convergence is proven. The selection method for the initial matrix is also provided. Finally, using the complexification operator of quaternion matrices, an equivalent iteration on the complex field is established to solve the equation in the Matlab environment. Two numerical examples are used to test the effectiveness and feasibility of the given method. .
基金Project supported by the National Natural Science Foundation of China(Grant No.11471262)the National Basic Research Program of China(Grant No.2012CB025904)the State Key Laboratory of Science and Engineering Computing and the Center for High Performance Computing of Northwestern Polytechnical University,China
文摘We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.
基金funded by the National Key Research and Development Program(Grant No.2022YFC3004200)the National Key Research and Development Program of China(Grant No.2021YFB3900400)+1 种基金Hunan Provincial Natural Science Foundation of China(Grant No.2021JC0009)the National Natural Science Foundation of China(Grant No.U2142212).
文摘The Geometrical Optics(GO)approach and the FAST Emissivity Model(FASTEM)are widely used to estimate the surface radiative components in atmospheric radiative transfer simulations,but their applications are limited in specific conditions.In this study,a two-scale reflectivity model(TSRM)and a two-scale emissivity model(TSEM)are developed from the two-scale roughness theory.Unlike GO which only computes six non-zero elements in the reflectivity matrix,The TSRM includes 16 elements of Stokes reflectivity matrix which are important for improving radiative transfer simulation accuracy in a scattering atmosphere.It covers the frequency range from L-to W-bands.The dependences of all TSRM elements on zenith angle,wind speed,and frequency are derived and analyzed in details.For a set of downwelling radiances in microwave frequencies,the reflected upwelling brightness temperature(BTs)are calculated from both TSRM and GO and compared for analyzing their discrepancies.The TSRM not only includes the effects of GO but also accounts for the small-scale Bragg scattering effect in an order of several degrees in Kelvins in brightness temperature.Also,the third and fourth components of the Stokes vector can only be produced from the TSRM.For the emitted radiation,BT differences in vertical polarization between a TSEM and FASTEM are generally less than 5 K when the satellite zenith angle is less than 40°,whereas those for the horizontal component can be quite significant,greater than 20 K.
基金supported by National Natural Science Foundation of China (10571047)and by Scientific Research Fund of Hunan Provincial Education Department of China Grant(06C235)+1 种基金by Central South University of Forestry and Technology (06Y017)by Specialized Research Fund for the Doctoral Program of Higher Education (20060532014)
文摘In this article, the generalized reflexive solution of matrix equations (AX = B, XC = D) is considered. With special properties of generalized reflexive matrices, the necessary and sufficient conditions for the solvability and the general expression of the solution are obtained. Moreover, the related optimal approximation problem to a given matrix over the solution set is solved.
基金partially supported by FCT(Portugal)with national funds through Centro de Matemática da Universidade de Trás-os-Montes e Alto Douro(PEst-OE/MAT/UI4080/2014)
文摘Dehghan and Hajarian, [4], investigated the matrix equations A^TXB+B^TX^TA = C and A^TXB + B^TXA = C providing inequalities for the determinant of the solutions of these equations. In the same paper, the authors presented a lower bound for the product of the eigenvalues of the solutions to these matrix equations. Inspired by their work, we give some generalizations of Dehghan and Hajarian results. Using the theory of the numerical ranges, we present an inequality involving the trace of C when A, B, X are normal matrices satisfying A^T B = BA^T.
基金supported by the National Natural Science Foundation of China (Grant No.60672160)
文摘We derive necessary and sufficient conditions for the existence and an expression of the (anti)reflexive solution with respect to the nontrivial generalized reflection matrix P to the system of complex matrix equations AX = B and XC = D. The explicit solutions of the approximation problem min x∈Ф ||X - E||F was given, where E is a given complex matrix and Ф is the set of all reflexive (or antireflexive) solutions of the system mentioned above, and ||·|| is the Frobenius norm. Furthermore, it was pointed that some results in a recent paper are special cases of this paper.
基金Supported by NECC and NSF of Shandong Proyilice,China(Y2005A06).
文摘Some new oscillation criteria are established for the second-order matrix differential system (r(t)Z'(t))' + p(t)Z'(t) + Q(t)F(Z'(t))G(Z(t)) = 0, t ≥ to 〉 0, are different from most known ones in the sense that they are based on the information only on a sequence of subintervals of [to,∞), rather than on the whole half-line. The results weaken the condition of Q(t) and generalize some well-known results of Wong (1999) to nonlinear matrix differential equation.
文摘In this paper we introduce the class of Hermite's matrix polynomials which appear as finite series solutions of second order matrix differential equations Y'-xAY'+BY=0.An explicit expression for the Hermite matrix polynomials,the orthogonality property and a Rodrigues' formula are given.
基金Supported by the NSF of Hubei Province(2022CFD042)。
文摘This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations into a system of algebraic equations.Meanwhile,the error analysis is proven.Finally,the effectiveness of the approach is verified by two numerical examples.
基金This work was supported by the Chinese Outstanding Youth Foundation(No.69925308)Program for Changjiang Scholars and Innovative ResearchTeam in University.
文摘In this paper, solutions to the generalized Sylvester matrix equations AX -XF = BY and MXN -X = TY with A, M ∈ R^n×n, B, T ∈ Rn×r, F, N ∈ R^p×p and the matrices N, F being in companion form, are established by a singular value decomposition of a matrix with dimensions n × (n + pr). The algorithm proposed in this paper for the euqation AX - XF = BY does not require the controllability of matrix pair (A, B) and the restriction that A, F do not have common eigenvalues. Since singular value decomposition is adopted, the algorithm is numerically stable and may provide great convenience to the computation of the solution to these equations, and can perform important functions in many design problems in control systems theory.
基金Project supported by the National Natural Science Foundation of China (Grant No.60672160)
文摘In this paper, the maximal and minimal ranks of the solution to a system of matrix equations over H, the real quaternion algebra, were derived. A previous known result could be regarded as a special case of the new result.
文摘A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. This is the first time that the inverse of this remarkable matrix is determined directly and exactly. Thus, solving 1D Poisson equation becomes very accurate and extremely fast. This method is a very important tool for physics and engineering where the Poisson equation appears very often in the description of certain phenomena.
文摘In this paper, two different methods are used to study the cyclic structure solution and the optimal approximation of the quaternion Stein equation AXB - X = F . Firstly, the matrix equation equivalent to the target structure matrix is constructed by using the complex decomposition of the quaternion matrix, to obtain the necessary and sufficient conditions for the existence of the cyclic solution of the equation and the expression of the general solution. Secondly, the Stein equation is converted into the Sylvester equation by adding the necessary parameters, and the condition for the existence of a cyclic solution and the expression of the equation’s solution are then obtained by using the real decomposition of the quaternion matrix and the Kronecker product of the matrix. At the same time, under the condition that the solution set is non-empty, the optimal approximation solution to the given quaternion circulant matrix is obtained by using the property of Frobenius norm property. Numerical examples are given to verify the correctness of the theoretical results and the feasibility of the proposed method. .
基金Project supported by the National Natural Science Foundation of China (Grant No.10271074), and the Special Funds for Major Specialities of Shanghai Education Commission (Grant No.J50101)
文摘In this paper, an improved gradient iterative (GI) algorithm for solving the Lyapunov matrix equations is studied. Convergence of the improved method for any initial value is proved with some conditions. Compared with the GI algorithm, the improved algorithm reduces computational cost and storage. Finally, the algorithm is tested with GI several numerical examples.
基金Supported by the Education Department Foundation of Hebei Province(QN2015218) Supported by the Natural Science Foundation of Hebei Province(A2015403050)
文摘Let P ∈ C^(n×n) be a Hermitian and {k + 1}-potent matrix, i.e., P^(k+1)= P = P~*,where(·)*~stands for the conjugate transpose of a matrix. A matrix X ∈ Cn×nis called{P, k + 1}-reflexive(anti-reflexive) if PXP = X(P XP =-X). The system of matrix equations AX = C, XB = D subject to {P, k + 1}-reflexive and anti-reflexive constraints are studied by converting into two simpler cases: k = 1 and k = 2, the least squares solution and the associated optimal approximation problem are also considered.
文摘An interesting semi-analytic solution is given for the Helmholtz equation. This solution is obtained from a rigorous discussion of the regularity and the inversion of the tridiagonal symmetric matrix. Then, applications are given, showing very good accuracy. This work provides also the analytical inverse of the skew-symmetric tridiagonal matrix.
文摘We derive necessary and sufficient conditions for the existence of a Hermitian nonnegative-definite solution to the matrix equation AXB = C. Moreover, we derive a representation of a general Hermitian nonnegative-definite solution. We then apply our solution to two examples, including a comparison of our solution to a proposed solution by Zhang in [1] using an example problem given from [1]. Our solution demonstrates that the proposed general solution from Zhang in [1] is incorrect. We also give a second example in which we derive the general covariance structure so that two matrix quadratic forms are independent.
文摘In this paper, we introduce a Hermite operational matrix collocation method for solving higher-order linear complex differential equations in rectangular or elliptic domains. We show that based on a linear algebra theorem, the use of different polynomials such as Hermite, Bessel and Taylor in polynomial collocation methods for solving differential equations leads to an equal solution, and the difference in the numerical results arises from the difference in the coefficient matrix of final linear systems of equations. Some numerical examples will also be given.