Matrix Continued Fractions

A matrix continued fraction is defined and used for the approximation of a function known as a power series in 1/zwith matrix coefficientsp×q, or equivalently by a matrix of functions holomorphic at infinity. It is a generalization of P-fractions, and the sequence of convergents converges to the given function. These convergents have as denominators a matrix, the columns of which are orthogonal with respect to the linear matrix functional associated to . The case where the algorithm breaks off is characterized in terms of .     

Clifford algebra and the Lorentz transformation with (p, q) form

In this paper, the concepts of Lorentz inner product with (p, q) form, the Lorentz space and the Lorentz transformation with (p, q) form are given by using Clifford algebra. It is shown that L_m^p,q is the Lorentz transformation with (p, q) form, and the matrix equality relation of Minkowski space with (n − 1, 1) form is given. The examples are given to illustrate the corresponding results.     

Conical algorithm for the global minimization of linearly constrained decomposable concave minimization problems

In this paper, we are concerned with the linearly constrained global minimization of the sum of a concave function defined on ap-dimensional space and a linear function defined on aq-dimensional space, whereq may be much larger thanp. It is shown that a conical algorithm can be applied in a space of dimensionp + 1 that involves only linear programming subproblems in a space of dimensionp +q + 1. Some computational results are given.This research was accomplished while the second author was a Fellow of the Alexander von Humboldt Foundation, University of Trier, Trier, Germany.     

On IOM(q): The Incomplete Orthogonalization Method for Large Unsymmetric Linear Systems

The incomplete orthogonalization method (IOM(q)), a truncated version of the full orthogonalization method (FOM) proposed by Saad, has been used for solving large unsymmetric linear systems. However, no convergence analysis has been given. In this paper, IOM(q) is analysed in detail from a theoretical point of view. A number of important results are derived showing how the departure of the matrix A from symmetric affects the basis vectors generated by IOM(q), and some relationships between the residuals for IOM(q) and FOM are established. The results show that IOM(q) behaves much like FOM once the basis vectors generated by it are well conditioned. However, it is proved that IOM(q) may generate an ill-conditioned basis for a general unsymmetric matrix such that IOM(q) may fail to converge or at least cannot behave like FOM. Owing to the mathematical equivalence between IOM(q) and the truncated ORTHORES(q) developed by Young and Jea, insights are given into the convergence of the latter. A possible strategy is proposed for choosing the parameter q involved in IOM(q). Numerical experiments are reported to show convergence behaviour of IOM(q) and of its restarted version.     

On a q-gamma and a q-beta matrix functions

The main goal of this article is to extend the concept of q-special functions of complex variable to q-special matrix functions through the study of a q-gamma and a q-beta matrix function. The q-shifted factorial, q-gamma and q-beta matrix functions are defined and some of their properties are investigated.     

A Feynman integral in Lifshitz‐point and Lorentz‐violating theories in

We evaluate a 1‐loop, 2‐point, massless Feynman integral I_D,m(p,q) relevant for perturbative field theoretic calculations in strongly anisotropic d=D+m dimensional spaces given by the direct sum . Our results are valid in the whole convergence region of the integral for generic (noninteger) codimensions D and m. We obtain series expansions of I_D,m(p,q) in terms of powers of the variable X:=4p²/q⁴, where p=| p |, q=| q |, , , and in terms of generalised hypergeometric functions ₃F₂(−X), when X<1. These are subsequently analytically continued to the complementary region X≥1. The asymptotic expansion in inverse powers of X^1/2 is derived. The correctness of the results is supported by agreement with previously known special cases and extensive numerical calculations.     

On irreducible p,q‐representations of Lie algebras \documentclass{article}\usepackage{amssymb}\begin{document}\pagestyle{empty}$\mathcal {G} (0,1)$\end{document} and \documentclass{article}\usepackage{amssymb}\begin{document}\pagestyle{empty}$\mathcal {G} (0,0)$\end{document}

In this paper, the irreducible p, q‐representations of the Lie algebras $\mathcal {G}(0,1)$ and $\mathcal {G}(0,0)$ are discussed. We prove two theorems that classify certain irreducible p, q‐representations of these Lie algebras and construct their one variable models in terms of p, q‐derivative and dilation operators. As an application, we derive a p, q‐special function identity based on one such model.     

Power comparisons of tests of two multivariate hypotheses based on individual characteristic roots

Summary In this paper, power comparisons are made for tests of each of the following two hypotheses based on individual characteristic roots of a matrix arising in each case: (i) independence between ap-set and aq-set of variates in a (p+q)-variate normal population withp≦q and (ii) equality ofp-dimensional mean vectors ofl p-variate normal populations having a common covariance matrix. At first, a few lemmas are given which help to reduce the central distributions of the largest, smallest, second largest, and the second smallest roots in terms of incomplete beta functions or functions of them. Since the central distribution of the largest root has been discussed by Pillai earlier in several papers ([6], [8], [9], [11], [12], [13]) cdf’s of the three others in the central case are given. Further, the non-central distributions of the individual roots forp-3 are considered for the two hypotheses and that of the smaller root forp=2; that of the largest root forp=2 has been obtained by Pillai earlier, (Pillai [11], Pillai and Jayachandran [14]). The work of this author was supported by the National Science Foundation Grant No. GP-4600.     

Matrix transformations of lq（X） to lP（Y）

Let X and Y be Banach spaces,0 ＜ q ＜ +∞,1 ＜ p ＜ +∞.In this paper,we characterize matrix transformations of lq(X) to lp(Y).     

Relative Position,Sequences and Operators in Banach Spaces

A (p + q) × (p + q) matrix-valued inner function S in the unit disc ?? is called (p, q)-type Arov-inner if in the block partition . the p × p diagonal block S₁₁ and the q × q diagonal block S₂₂ are outer matrix-valued functions. A holomorphic p × q matrix-valued function f in ?? is called Arov-completable if there is a (p, q)-type Arov-inner function S such that S₁₂ = f Arov-completability of a given p × q Schur function f is characterized in terms of a (p + q)-variate stationary sequence (X_{n) ? Z}) in Hilbert space which is naturally associated with f. The necessary and sufficient condition for Arov-completability is an orthogonality condition for certain backward and forward innovation vectors generated by (X_{n) ? Z}.     

A Note on the Completeness of an Exponential Type Sequence

For any given coprime integers p and q greater than 1, in 1959, B. J. Birch proved that all sufficiently large integers can be expressed as a sum of pairwise distinct terms of the form paqb. As Davenport observed, Birch’s proof can be modified to show that the exponent b can be bounded in terms of p and q. In 2000, N. Hegyvari gave an effective version of this bound. The author improves this bound.     

A condition for arcs and MDS codes

There is polynomial function X _q in the entries of an m × m(q − 1) matrix over a field of prime characteristic p, where q = p ^h is a power of p, that has very similar properties to the determinant of a square matrix. It is invariant under multiplication on the left by a non-singular matrix, and under permutations of the columns. This gives a way to extend the invariant theory of sets of points in projective spaces of prime characteristic, to make visible hidden structure. There are connections with coding theory, permanents, and additive bases of vector spaces.     

Asymptotics of partition functions in a fermionic matrix model and of related q‐polynomials

In this paper, we study asymptotics of the thermal partition function of a model of quantum mechanical fermions with matrix‐like index structure and quartic interactions. This partition function is given explicitly by a Wronskian of the Stieltjes‐Wigert polynomials. Our asymptotic results involve the theta function and its derivatives. We also develop a new asymptotic method for general q‐polynomials.     

On the mutually non isomorphic ℓp(ℓq) spaces

We extend a result of Pe?czyński showing that {?_p(?_q): 1 ≤ p, q ≤ ∞} is a family of mutually non isomorphic Banach spaces. Some results on complemented subspaces of ?_p(?_q) are also given.     

Simultaneous Polynomial Approximation in the Weighted Chebyshev Norm

Simultaneous approximation means that a given sufficiently smooth function g:[-1, 1] and its derivatives up to order q are approximated by a single polynomial p and its derivatives. This paper deals with new error estimates (in a weighted norm with explicit constants) and corresponding algorithms in the most interesting cases q = 1 and q = 2. The described method is based on the close relationship between algebraic and trigonometric polynomial approximation.     

Lp-Computability

In this paper we investigate conditions for L^p-computability which are in accordance with the classical Grzegorczyk notion of computability for a continuous function. For a given computable real number p ≥ 1 and a compact computable rectangle I ? ?^q, we show that an L^p function f ∈ L^p(I) is L^P-computable if and only if (i) f is sequentially computable as a linear functional and (ii) the L^p-modulus function of f is effectively continuous at the origin of ?^q.     

Continued fraction expansions of matrix eigenvectors

We examine various properties of the continued fraction expansions of matrix eigenvector slopes of matrices from the SL(2, ℤ) group. We calculate the average period length, maximum period length, average period sum, maximum period sum, and the distributions of 1s, 2s, and 3s in the periods versus the radius of the ball within which the matrices are located. We also prove that the periods of continued fraction expansions from the real irrational roots of x ²+px+q=0 are always palindromes.      

Approximate method for eigensensitivity analysis of a defective matrix

Based on the exact modal expansion method, an arbitrary high-order approximate method is developed for calculating the second-order eigenvalue derivatives and the first-order eigenvector derivatives of a defective matrix. The numerical example shows the validity of the method. If the different eigenvalues μ(1),…,μ(q) of the matrix are arranged so that |μ(1)|≤?≤|μ(q)| and satisfy the condition that |μ(q₁)|<|μ(q₁+1)| for some q₁<q, and if the approximate method only uses the left and right principal eigenvectors associated with μ(1),…,μ(q₁), then associated with μ(h)(h≤q₁) the errors of the eigenvalue and eigenvector derivatives by the pth-order approximate method are nearly proportional to |μ(h)/μ(q₁+1)|^p+1.     

Identification of Refined ARMA Echelon Form Models for Multivariate Time Series

In the present article, we are interested in the identification of canonical ARMA echelon form models represented in a “refined” form. An identification procedure for such models is given by Tsay (J. Time Ser. Anal.10(1989), 357-372). This procedure is based on the theory of canonical analysis. We propose an alternative procedure which does not rely on this theory. We show initially that an examination of the linear dependency structure of the rows of the Hankel matrix of correlations, with originkin (i.e., with correlation at lagkin position (1, 1)), allows us not only to identify the Kronecker indicesn₁, …, n_d, whenk=1, but also to determine the autoregressive ordersp₁, …, p_d, as well as the moving average ordersq₁, …, q_dof the ARMA echelon form model by settingk>1 andk<1, respectively. Successive test procedures for the identification of the structural parametersn_i,p_i, andq_iare then presented. We show, under the corresponding null hypotheses, that the test statistics employed asymptotically follow chi-square distributions. Furthermore, under the alternative hypothesis, these statistics are unbounded in probability and are of the formNδ{1+o_p(1)}, whereδis a positive constant andNdenotes the number of observations. Finally, the behaviour of the proposed identification procedure is illustrated with a simulated series from a given ARMA model.     

Solutions in Sobolev spaces of vector refinement equations with a general dilation matrix

In this paper, we present a necessary and sufficient condition for the existence of solutions in a Sobolev space W_p^k(ℝ^s) (1≤p≤∞) to a vector refinement equation with a general dilation matrix. The criterion is constructive and can be implemented. Rate of convergence of vector cascade algorithms in a Sobolev space W_p^k(ℝ^s) will be investigated. When the dilation matrix is isotropic, a characterization will be given for the L_p (1≤p≤∞) critical smoothness exponent of a refinable function vector without the assumption of stability on the refinable function vector. As a consequence, we show that if a compactly supported function vector φ∈L_p(ℝ^s) (φ∈C(ℝ^s) when p=∞) satisfies a refinement equation with a finitely supported matrix mask, then all the components of φ must belong to a Lipschitz space Lip(ν,L_p(ℝ^s)) for some ν>0. This paper generalizes the results in R.Q. Jia, K.S. Lau and D.X. Zhou (J. Fourier Anal. Appl. 7 (2001) 143–167) in the univariate setting to the multivariate setting. Dedicated to Professor Charles A. Micchelli on the occasion of his 60th birthday Mathematics subject classifications (2000) 42C20, 41A25, 39B12. Research was supported in part by the Natural Sciences and Engineering Research Council of Canada (NSERC Canada) under Grant G121210654.