Block diagonalization and LU-equivalence of Hankel matrices

This article presents a new algorithm for obtaining a block diagonalization of Hankel matrices by means of truncated polynomial divisions, such that every block is a lower Hankel matrix. In fact, the algorithm generates a block LU-factorization of the matrix. Two applications of this algorithm are also presented. By the one hand, this algorithm yields an algebraic proof of Frobenius’ Theorem, which gives the signature of a real regular Hankel matrix by using the signs of its principal leading minors. On the other hand, the close relationship between Hankel matrices and linearly recurrent sequences leads to a comparison with the Berlekamp–Massey algorithm.     

A Generalized Rational Interpolation Problem and the Solution of the Welch–Berlekamp Key Equation

We show that an algorithm designed to solve the Welch–Berlekamp key equation may also be used to solve a more general problem, which can be regarded as a finite analogue of a generalized rational interpolation problem. As a consequence, we show that a single algorithm exists which can solve both Berlekamp's classical key equation (usually solved by the Berlekamp–Massey algorithm) and the Welch–Berlekamp key equation which arise in the decoding of Reed–Solomon codes.     

Berlekamp—Massey Algorithm, Continued Fractions, Pade Approximations, and Orthogonal Polynomials

The Berlekamp—Massey algorithm (further, the BMA) is interpreted as an algorithm for constructing Pade approximations to the Laurent series over an arbitrary field with singularity at infinity. It is shown that the BMA is an iterative procedure for constructing the sequence of polynomials orthogonal to the corresponding space of polynomials with respect to the inner product determined by the given series. The BMA is used to expand the exponential in continued fractions and calculate its Pade approximations.     

Fractional integral operators in the complex matrix variate case

A formal definition of fractional integrals in the complex matrix variate case is given here. This definition will encompass all the various fractional integral operators introduced by various authors in the real scalar and matrix cases. The new definition is introduced in terms of M-convolutions of products and ratios of matrices in the complex domain. Their connections to statistical distribution theory, Mellin convolutions, M-transforms and Mellin transform are pointed out. Some basic properties are given and a pathway extension of the new definition is also given. The pathway extension will provide a switching mechanism to move among three different families of functions.     

Products of idempotent matrices over integral domains

We say that a ring R has the idempotent matrices property if every square singular matrix over R is a product of idempotent matrices. It is known that every field, and more generally, every Euclidean domain has the idempotent matrices property. In this paper we show that not every integral domain has the idempotent matrices property and that if a projective free ring has the idempotent matrices property then it must be a Bezout domain. We also show that a principal ideal domain has the idempotent matrices property if and only if every fraction a/b with b≠0 has a finite continued fraction expansion. New proofs are also provided for the results that every field and every Euclidean domain have the idempotent matrices property.     

Refined convergents to the associated continued fractions for binary sequences

The relation between continued fractions and Berlekamp's algorithm was studied by some researchers. The latter is an iterative procedure proposed for decoding BCH codes. However, there remains an unanswered question wheter each of the iterative steps in the algorithm can be interpreted in terms of continued fractions. In this paper, we first introduce the so-called refined convergents to the continued fraction expansion of a binary sequence S, and then give a thorough answer to the question in the context of Massey's linear feedback shift register synthesis algorithm which is equivalent to that of Berlekamp, and at last we prove that there exists a one-to-one correspondence between then-th refined convergents and the lengthn segments.     

Riemann problem and singular integral equations with coefficients generated by piecewise constant functions

We study a Riemann boundary value problem with a shift into the interior of the domain. The problem has piecewise constant coefficients that take two values. We find conditions for the existence and uniqueness of a solution of the inhomogeneous problem and formulas for the number of linearly independent solutions of the homogeneous problem. We consider scalar singular integral operators with a shift and matrix characteristic operators whose coefficients are generated by piecewise constant functions and which have automorphic properties. For these operators, we find invertibility conditions.     

Wavelet-like bases for thin-wire integral equations in electromagnetics

In this paper, wavelets are used in solving, by the method of moments, a modified version of the thin-wire electric field integral equation, in frequency domain. The time domain electromagnetic quantities, are obtained by using the inverse discrete fast Fourier transform. The retarded scalar electric and vector magnetic potentials are employed in order to obtain the integral formulation. The discretized model generated by applying the direct method of moments via point-matching procedure, results in a linear system with a dense matrix which have to be solved for each frequency of the Fourier spectrum of the time domain impressed source. Therefore, orthogonal wavelet-like basis transform is used to sparsify the moment matrix. In particular, dyadic and M-band wavelet transforms have been adopted, so generating different sparse matrix structures. This leads to an efficient solution in solving the resulting sparse matrix equation. Moreover, a wavelet preconditioner is used to accelerate the convergence rate of the iterative solver employed. These numerical features are used in analyzing the transient behavior of a lightning protection system. In particular, the transient performance of the earth termination system of a lightning protection system or of the earth electrode of an electric power substation, during its operation is focused. The numerical results, obtained by running a complex structure, are discussed and the features of the used method are underlined.     

A Series Approach to Stochastic Differential Equations with Infinite Dimensional Noise

This paper considers semilinear stochastic differential equations in Hilbert spaces with Lipschitz nonlinearities and with the noise terms driven by sequences of independent scalar Wiener processes (Brownian motions). The interpretation of such equations requires a stochastic integral. By means of a series of Itô integrals, an elementary and direct construction of a Hilbert space valued stochastic integral with respect to a sequence of independent scalar Wiener processes is given. As an application, existence and strong and weak uniqueness for the stochastic differential equation are shown by exploiting the series construction of the integral.     

Maximally movable Riemannian spaces with torsion

We prove that, among all Riemannian spaces of constant curvature, only three-dimensional spaces have torsion which is invariant under the group of motions. The torsion tensor in these spaces is covariantly constant and determines the torsion form. The ratio of the integral of this form over a bounded domain to its volume is a constant determining the torsion of the space. We introduce the notions of volume torsion and scalar torsion.     

DFT calculation for the {2}-inverse of a polynomial matrix with prescribed image and kernel

In this paper, we first give the finite algorithm for generalized inverse of a matrix A over an integral domain, and, based on it and the discrete Fourier transform, present an algorithm for calculating {2}-inverses of a polynomial matrix with prescribed image and kernel. And the algorithm is implemented in the Mathematica programming language and expands the algorithms in [13].     

A normal crack in an elastic half-space with stress-free surface

The problem of an elastic half-space with stress-free surface and a crack of arbitrary shape with prescribed displacements or tractions is reduced to an equivalent system of integral equations on the crack. For a pressurized crack in a plane perpendicular to the free surface, a scalar integral equation is derived. In properly chosen function spaces, unique solvability of the integral equation and regularity of solutions for regular data are proven.     

MATRIX ALGORITHMS AND ERROR FORMULA FOR BIVARIATE THIELE-TYPE RECTANGULAR MATRIX VALUED RATIONAL INTERPOLATION

A new method for the construction of bivariate matrix valued rational interpolants (BGIRI) on a rectangular grid is presented in [6]. The rational interpolants are of Thiele-type continued fraction form with scalar denominator. The generalized inverse introduced by [3]is gen-eralized to rectangular matrix case in this paper. An exact error formula for interpolation is ob-tained, which is an extension in matrix form of bivariate scalar and vector valued rational interpola-tion discussed by Siemaszko[l2] and by Gu Chuangqing [7] respectively. By defining row and col-umn-transformation in the sense of the partial inverted differences for matrices, two type matrix algorithms are established to construct corresponding two different BGIRI, which hold for the vec-tor case and the scalar case.     

Method of Integral Equations for the Three-Dimensional Problem of Wave Reflection from an Irregular Surface

The problem on the reflection of the field of a plane H-polarized three-dimensional electromagnetic wave from a perfectly conducting interface between media which contains a local perfectly conducting inhomogeneity is considered. To construct a numerical algorithm, the boundary value problem for the system of Maxwell equations in an infinite domain with irregular boundary is reduced to a system of singular integral equations, which is solved by the approximation–collocation method. The elements of the resulting complex matrix are calculated by a specially developed algorithm. The solution of the system of singular integral equations is used to obtain an integral representation for the reflected electromagnetic field and computational formulas for the directional diagram of the reflected electromagnetic field in the far region.     

Functional integral representations of the Pauli-Fierz model with spin 1/2

A Feynman-Kac-type formula for a Lévy and an infinite-dimensional Gaussian random process associated with a quantized radiation field is derived. In particular, a functional integral representation of e−tHPF generated by the Pauli-Fierz Hamiltonian with spin 1/2 in non-relativistic quantum electrodynamics is constructed. When no external potential is applied HPF turns translation-invariant and it is decomposed as a direct integral . The functional integral representation of e−tHPF(P) is also given. Although all these Hamiltonians include spin, nevertheless the kernels obtained for the path measures are scalar rather than matrix expressions. As an application of the functional integral representations energy comparison inequalities are derived.     

Characteristic polynomials of symmetric matrices III: some counterexamples

When is a monic polynomial the characteristic polynomial of a symmetric matrix over an integral domain D? Known necessary conditions are shown to be insufficient when D is the field of 2-adic numbers and when D is the rational integers. The latter counterexamples lead to totally real cubic extensions of the rationals whose difierents are not narrowly equivalent to squares. Furthermorex³-4x+1 is the characteristic polynomial of a rational symmetric matrix and is the characteristic polynomial of an integral symmetric p-adic matrix for every prime p, but is not the characteristic polynomial of a rational integral symmetric matrix.