Search

In this paper, we present a statistical approach to evaluate the relationship between variables observed in a two-factors experiment. We consider a three-level model with covariance structure Σ ⊗ Ψ 1 ⊗ Ψ 2 , where Σ is an arbitrary positive definite covariance matrix, and Ψ 1 and Ψ 2 are both correlation matrices with a compound symmetric structure corresponding to two different factors. The Rao's score test is used to test the hypotheses that observations grouped by one or two factors are uncorrelated. We analyze a fermentation process to illustrate the results. Supplementary materials accompanying this paper appear online. [ABSTRACT FROM AUTHOR]

Frieze patterns have been introduced by Coxeter [Acta Arith. 18 (1971), pp. 297–310] in the 1970's and have recently attracted renewed interest due to their close connection with Fomin-Zelevinsky's cluster algebras. Frieze patterns can be interpreted as assignments of values to the diagonals of a triangulated polygon satisfying certain conditions for crossing diagonals (Ptolemy relations). Weak friezes, as introduced by Çanakçı and Jørgensen [Adv. in Appl. Math. 131 (2021), Paper No. 102253], are generalizing this concept by allowing to glue dissected polygons so that the Ptolemy relations only have to be satisfied for crossings involving one of the gluing diagonals. To any frieze pattern one can associate a symmetric matrix using a triangular fundamental domain of the frieze pattern in the upper and lower half of the matrix and putting zeroes on the diagonal. Broline, Crowe and Isaacs [Geometriae Dedicata 3 (1974), pp. 171–176] have found a formula for the determinants of these matrices and their work has later been generalized in various directions by other authors. These frieze pattern determinants are the main focus of our paper. As our main result we show that this determinant behaves well with respect to gluing weak friezes: the determinant is the product of the determinants for the pieces glued, up to a scalar factor coming from the gluing diagonal. Then we give several applications of this result, showing that formulas from the literature, obtained by Broline-Crowe-Isaacs, Baur-Marsh [J. Combin. Theory Ser. A 119 (2012), pp. 1110–1122], Bessenrodt-Holm-Jørgensen [J. Combin. Theory Ser. A 123 (2014), pp. 30–42] and Maldonado [ Frieze matrices and infinite frieze patterns with coefficients , Preprint, arXiv: 2207.04120 , 2022] can all be obtained as consequences of our result. [ABSTRACT FROM AUTHOR]

The rational field Q is highly desired in many applications. Algorithms using the rational number field Q algebraic number fields use only integer arithmetics and are easy to implement. Therefore, studying and designing systems and expansions with coefficients in Q or algebraic number fields is particularly interesting. This paper discusses constructing quasi-tight framelets with symmetry over an algebraic field. Compared to tight framelets, quasi-tight framelets have very similar structures but much more flexibility in construction. Several recent papers have explored the structure of quasi-tight framelets. The construction of symmetric quasi-tight framelets directly applies the generalized spectral factorization of 2 × 2 matrices of Laurent polynomials with specific symmetry structures. We adequately formulate the latter problem and establish the necessary and sufficient conditions for such a factorization over a general subfield F of C , including algebraic number fields as particular cases. Our proofs of the main results are constructive and thus serve as a guideline for construction. We provide several examples to demonstrate our main results. [ABSTRACT FROM AUTHOR]

Matrix factorization is an inference problem that has acquired importance due to its vast range of applications that go from dictionary learning to recommendation systems and machine learning with deep networks. The study of its fundamental statistical limits represents a true challenge, and despite a decade-long history of efforts in the community, there is still no closed formula able to describe its optimal performances in the case where the rank of the matrix scales linearly with its size. In the present paper, we study this extensive rank problem, extending the alternative 'decimation' procedure that was recently introduced by Camilli and Mézard, and carry out a thorough study of its performance. Decimation aims at recovering one column/line of the factors at a time, by mapping the problem into a sequence of neural network models of associative memory at a tunable temperature. The main advantage of decimation is that its theoretical performance can be studied using statistical physics techniques. In this paper we provide the general analysis applying to a large class of compactly supported priors and we show that the replica symmetric free entropy of the neural network models takes a universal form in the low temperature limit. We then study the decimation phase diagrams in two concrete cases, the sparse Ising prior and the uniform prior, for which we show that matrix factorization is theoretically possible when the ration P / N of rank to variable is below a certain threshold. This suggest that a possible route to solving algorithmically the matrix factorization problem could be to find efficient decimation algorithms; in this respect, we show that simple simulated annealing, though being effective for limited signal sizes, is not scalable. [ABSTRACT FROM AUTHOR]

The symmetric nonnegative inverse eigenvalue problem (SNIEP) is the problem of characterizing all possible spectra of entry-wise nonnegative symmetric matrices of given dimension. A list of real numbers is said to be symmetrically realizable if it is the spectrum of some nonnegative symmetric matrix. One of the most general sufficient conditions for realizability is the so-called C-realizability, which amounts to some kind of compensation between the positive and negative entries of the list of real numbers whose realizability one is trying to decide. A combinatorial characterization of C-realizable lists with zero trace was given in [11]. In this paper we make use of a recursive method for constructing simmetrically realizable lists due to Ellard and Šmigoc [3] to extend this combinatorial characterization of C-realizability to general lists with nonnegative trace. One consequence of this characterization is that the set of nonnegative C-realizable lists is a union of polyhedral cones whose faces are described by equations involving only linear combinations with coefficients 1 and −1 of the entries in the list. Another remarkable consequence is the monotonicity of C-realizability, i.e., the operation of increasing any positive entry of a C-realizable list preserves C-realizability. [ABSTRACT FROM AUTHOR]

In this paper, we analytically derive closed-form expressions for the tangency portfolio weights: the fully invested portfolio that maximizes the expected return over the risk-free rate, relative to the volatility of the portfolio return. We explicitly derive this portfolio from a range of underlying return models and show examples where it coincides with different well-known smart beta products. Specifically, we find the closed-form expression for the tangency portfolio weights for a return model with compound symmetric correlation matrix. We also deduce the tangency portfolio weights for the CAPM return model and illustrate in a case study that the estimated tangency portfolio weights may distinctly deviate from the market value weighted portfolio. Furthermore, we show that depending on the return model, the tangency portfolio weights may take a diverse set of shapes; from very diversified to highly concentrated portfolios. [ABSTRACT FROM AUTHOR]

Let G be a connected graph of order n, where n is a positive integer. A spanning subgraph F of G is called a path-factor if every component of F is a path of order at least 2. A P ≥ k -factor means a path-factor in which every component admits order at least k ( k ≥ 2 ). The distance matrix D (G) of G is an n × n real symmetric matrix whose (i, j)-entry is the distance between the vertices v i and v j . The distance signless Laplacian matrix Q (G) of G is defined by Q (G) = T r (G) + D (G) , where Tr(G) is the diagonal matrix of the vertex transmissions in G. The largest eigenvalue η 1 (G) of Q (G) is called the distance signless Laplacian spectral radius of G. In this paper, we aim to present a distance signless Laplacian spectral radius condition to guarantee the existence of a P ≥ 2 -factor in a graph and claim that the following statements are true: (i) G admits a P ≥ 2 -factor for n ≥ 4 and n ≠ 7 if η 1 (G) < θ (n) , where θ (n) is the largest root of the equation x 3 - (5 n - 3) x 2 + (8 n 2 - 23 n + 48) x - 4 n 3 + 22 n 2 - 74 n + 80 = 0 ; (ii) G admits a P ≥ 2 -factor for n = 7 if η 1 (G) < 25 + 161 2 . [ABSTRACT FROM AUTHOR]

Let r and s be positive integers such that r ⩾ 2. Let U and V be vector spaces over fields F and K , respectively, such that dim ⁡ U ⩾ 3 and F has at least r + 1 elements. In this paper, we characterize surjective maps ψ : S r U → S s V preserving adjacency in both directions on symmetric tensors of finite order, which generalizes Hua's fundamental theorem of geometry of symmetric matrices. We give examples showing the indispensability of the assumptions dim ⁡ U ⩾ 3 and the cardinality | F | ⩾ r + 1 in our result. [ABSTRACT FROM AUTHOR]

The parameterized block splitting (PBS) is a convergent and efficient iterative method to solve the large complex symmetric linear systems. In this paper, by using PBS iterative technique, the Newton equation is approximately solved, then we establish the modified Newton-PBS iterative method to solve the complex nonlinear systems whose Jacobian matrices are large, sparse, and complex symmetric. Subsequently, the local convergence analysis are explored under appropriate conditions. Ultimately, we apply the new method and several known methods to experimental numerical examples, and experimental results verify the superiority and efficiency of our new method. Especially, in terms of CPU time and iteration steps, our method is obviously better. [ABSTRACT FROM AUTHOR]

In this paper, we advocate a new technique to determine inner bounds for the extreme eigenvalues of real symmetric matrices. Our method involves the matrix elements and compares favourably with existing methods. We also show how the bounds can be optimized. [ABSTRACT FROM AUTHOR]

Large-scale symmetric and asymmetric matrices have emerged in predicting the relationship between genes and diseases. The emergence of large-scale matrices increases the computational complexity of the problem. Therefore, using low-rank matrices instead of original symmetric and asymmetric matrices can greatly reduce computational complexity. In this paper, we propose an approximation conjugate gradient method for solving the low-rank matrix recovery problem, i.e., the low-rank matrix is obtained to replace the original symmetric and asymmetric matrices such that the approximation error is the smallest. The conjugate gradient search direction is given through matrix addition and matrix multiplication. The new conjugate gradient update parameter is given by the F-norm of matrix and the trace inner product of matrices. The conjugate gradient generated by the algorithm avoids SVD decomposition. The backtracking linear search is used so that the approximation conjugate gradient direction is computed only once, which ensures that the objective function decreases monotonically. The global convergence and local superlinear convergence of the algorithm are given. The numerical results are reported and show the effectiveness of the algorithm. [ABSTRACT FROM AUTHOR]

In this paper, we generalize the Pascal matrices presented in [9] as a functional for several variables. The entries of these matrices are defined on the set R ⊆ Z ≥ 0 n with multidimensional binomial coefficients along with several variables. In the following, we will investigate the algebraic properties of these matrices in terms of products, powers, and inverses of the matrices and also some factorization formulas using this method. [ABSTRACT FROM AUTHOR]

In this paper, we deal with the construction of symmetric matrix whose corresponding graph is connected and unicyclic using some pre-assigned spectral data. Spectral data for the problem consist of the smallest and the largest eigenvalues of each leading principal submatrices. Inverse eigenvalue problem (IEP) with this set of spectral data is generally known as the extremal IEP. We use a standard scheme of labeling the vertices of the graph, which helps in getting a simple relation between the characteristic polynomials of each leading principal submatrix. Sufficient condition for the existence of the solution is obtained. The proof is constructive, hence provides an algorithmic procedure for finding the required matrix. Furthermore, we provide the condition under which the same problem is solvable when two particular entries of the required matrix satisfy a linear relation. [ABSTRACT FROM AUTHOR]

In this paper, we study parameter deformations of matrix valued orthogonal polynomials. These deformations are built on the use of certain matrix valued operators which are symmetric with respect to the matrix valued inner product defined by the orthogonality weight. We show that the recurrence coefficients associated with these operators satisfy generalizations of the non-Abelian lattice equations. We provide a Lax pair formulation for these equations, and an example of deformed Hermite-type matrix valued polynomials is discussed in detail. [ABSTRACT FROM AUTHOR]

Let G be a simple connected graph on n vertices, and let G σ be an orientation of G with skew adjacency matrix S (G σ). Let d i be the degree of the vertex v i in G. The skew Randić matrix of G σ is the n × n real skew symmetric matrix R S (G σ) = [ (R S) i j ] , where (R S) i j = − (R S) j i = (d i d j) − 1 2 if (v i , v j) is an arc of G σ , and (R S) i j = (R S) j i = 0 otherwise. The skew spectral radius ρ S (G σ) and the skew Randić spectral radius ρ R S (G σ) of G σ are defined as the spectral radius of S (G σ) and R S (G σ) respectively. In this paper we give upper bounds for the skew spectral radius and skew Randić spectral radius of general random oriented graphs. [ABSTRACT FROM AUTHOR]

An efficient iteration method is provided in this paper for solving a class of nonlinear systems whose Jacobian matrices are complex and symmetric. The modified Newton-NDSS method is developed and applied to the class of nonlinear systems by adopting the modified Newton method as the outer solver and a new double-step splitting (NDSS) iteration scheme as the inner solver. Additionally, we theoretically analyze the local convergent properties of the new method under the weaker Hölder conditions. Lastly, the new method is compared numerically with some existing ones and the numerical experiments solving the nonlinear equations demonstrate the superiority of the Newton-NDSS method. [ABSTRACT FROM AUTHOR]

Nowadays, deep representations have gained significant attention due to their outstanding performance in a wide range of tasks. However, the interpretability of deep representations in specific applications poses a significant challenge. For instances where the generated quantity matrices exhibit symmetry, this paper introduces a variant of deep matrix factorization (deep MF) called deep L p smooth symmetric non-negative matrix factorization (DSSNMF), which aims to improve the extraction of clustering structures inherent in complex hierarchical and graphical representations in high-dimensional datasets by improving the sparsity of the factor matrices. We successfully applied DSSNMF to synthetic datasets as well as datasets related to post-traumatic stress disorder (PTSD) to extract several hierarchical communities. Specifically, we identified non-disjoint communities within the partial correlation networks of PTSD psychiatric symptoms, resulting in highly meaningful clinical interpretations. Numerical experiments demonstrate the promising applications of DSSNMF in fields like network analysis and medicine. [ABSTRACT FROM AUTHOR]

In this paper we consider n coupled Schrödinger equations − d 2 y d t 2 + u (ω t) y = E y , y ∈ R n , where E = diag ( λ 1 2 , ... , λ n 2 ) is a diagonal matrix, u(ωt) is a real analytic quasi-periodic symmetric matrix. If the basic frequencies ω = (1, α), where α is irrational, it is proved that for most of sufficiently large λj, j = 1, ..., n, all the solutions of n coupled Schrödinger equations are bounded. Furthermore, if the basic frequencies satisfy that 0 ≤ β(α) < r, where β (α) = lim sup n → ∞ ln q n + 1 q n , qn is the sequence of denominations of the best rational approximations for α ∈ R \ Q , r is the initial analytic radius, we obtain the existence of n pairs of conjugate quasi-periodic solutions for most of sufficiently large λj, j = 1, ..., n. [ABSTRACT FROM AUTHOR]

An n-by-n real symmetric matrix $ H=[h_{i,j}] $ H = [ h i , j ] is said to be a Hankel matrix if $ h_{i,j}=h_{i-1,j+1} $ h i , j = h i − 1 , j + 1 , for each $ i=2,\ldots,n $ i = 2 , ... , n and $ j=1,\ldots,n-1 $ j = 1 , ... , n − 1. The symmetric nonnegative inverse eigenvalue problem (SNIEP) asks when a list $ \Lambda =\{\lambda _{1},\ldots,\lambda _{n}\} $ Λ = { λ 1 , ... , λ n } of real numbers is the spectrum of an n-by-n symmetric nonnegative matrix H. In this paper, we search for conditions on the list $ \Lambda =\{\lambda _{1},\ldots,\lambda _{n}\} $ Λ = { λ 1 , ... , λ n } for the matrix H to be Hankel. For n = 3, sufficient conditions are established. In particular, a necessary and sufficient condition is obtained if Λ is a list of three nonnegative numbers. Also, if $ \sum _{i=1}^{n}\lambda _{i}=0 $ ∑ i = 1 n λ i = 0 , we give conditions for realizability by a Hankel matrix. Finally, we present a special type of list that can serve as the spectrum of a Hankel nonnegative matrix with positive trace. Several of our results are constructive and provide a Hankel realizing matrix. [ABSTRACT FROM AUTHOR]

The spectral radius ρ (A) is the maximum absolute value of the eigenvalues of a matrix A. In this paper, we establish some relationship between the spectral radius of a symmetric matrix and its principal submatrices, i.e., if A is partitioned as a 2 × 2 block matrix A = ( 0 A 12 A 21 A 22 ) , then ρ (A) 2 ≤ ρ 2 2 + θ ⁎ , where θ ⁎ is the largest real root of the equation μ 2 = (x − ν) 2 (ρ 2 2 + x) and ρ 2 = ρ (A 22) , μ = ρ (A 12 A 22 A 21) , ν = ρ (A 12 A 21). Furthermore, the results are used to obtain several upper bounds of the spectral radius of graphs, which strengthen or improve some known results. [ABSTRACT FROM AUTHOR]

The study of regular incidence structures, such as projective planes and symmetric block designs, is a well established topic in discrete mathematics. The work of Bruck, Ryser and Chowla in the mid-twentieth century applied the Hasse-Minkowski local-global theory for quadratic forms to derive non-existence results for certain design parameters. Several combinatorialists have provided alternative proofs of this result, replacing conceptual arguments with algorithmic ones. In this paper, we show that the methods required are purely linear-algebraic and are no more difficult conceptually than the theory of the Jordan Canonical Form. Computationally, they are rather easier. We conclude with some classical and recent applications to design theory, including a novel application to the decomposition of incidence matrices of symmetric designs. [ABSTRACT FROM AUTHOR]

This paper studies a priori and regularity estimates of Evans-Krylov type in Hölder spaces for fully nonlinear uniformly elliptic and parabolic equations of second order when the operator fails to be concave or convex in the space of symmetric matrices. In particular, it is assumed that either the level sets are convex or the operator is concave, convex or close to a linear function near infinity. As a byproduct, these results imply polynomial Liouville theorems for entire solutions of elliptic equations and for ancient solutions to parabolic problems. [ABSTRACT FROM AUTHOR]

Let S be the polynomial ring on the space of non-square generic matrices or the space of odd-sized skew-symmetric matrices, and let I be the determinantal ideal of maximal minors or Pf the ideal of sub-maximal Pfaffians, respectively. Using desingularizations and representation theory of the general linear group we expand upon work of Raicu–Weyman–Witt [Adv. Math. 250 (2014), pp. 596–610] to determine the S-module structures of Ext^j_S(S/I^t, S) and Ext^j_S(S/Pf^t, S), from which we get the degrees of generators of these Ext modules. As a consequence, via graded local duality we answer a question of Wenliang Zhang [J. Pure Appl. Algebra 225 (2021), Paper No. 106789] on the socle degrees of local cohomology modules of the form H^j_\mathfrak {m}(S/I^t). [ABSTRACT FROM AUTHOR]

( R , S )-(skew) symmetric matrices have numerous applications in civil engineering, information theory, numerical analysis, etc. In this paper, we deal with the ( R , S )-(skew) symmetric solutions to the quaternion matrix equation A X B = C . We use a real representation A τ to obtain the necessary and sufficient conditions for A X B = C to have ( R , S )-(skew) symmetric solutions and derive the solutions when it is consistent. We also derive the least-squares ( R , S )-(skew) symmetric solution to the above matrix equation. [ABSTRACT FROM AUTHOR]

The paper is a survey and concerns with infinite symmetric block Jacobi matrices J with m × m-matrix entries. We discuss several results on general block Jacobi matrices to be either self-adjoint or have maximal as well as intermediate deficiency indices. We also discuss several conditions for J to have discrete spectrum. [ABSTRACT FROM AUTHOR]

In this paper we are concerned with plane wave discontinuous Galerkin (PWDG) methods for time-harmonic Maxwell equations in three-dimensional anisotropic media, for which the coefficients of the equations are piecewise constant symmetric matrices, where each constant symmetric matrix is defined on a medium (subdomain). By using suitable scaling transformations and coordinate (complex) transformations on every subdomain, the original Maxwell equation in anisotropic media is transformed into a Maxwell equation in isotropic media occupying a union domain of specific subdomains of complex Euclidean space. Based on these transformations, we define anisotropic plane wave basis functions and discretize the considered Maxwell equations by PWDG method with the proposed plane wave basis functions. We derive error estimates of the resulting approximate solutions, and further introduce a practically feasible local hp-refinement algorithm, which substantially improves accuracies of the approximate solutions. Numerical results indicate that the approximate solutions generated by the proposed PWDG methods possess high accuracy for the case of strong discontinuity media. [ABSTRACT FROM AUTHOR]

In this paper, we consider a coupled Brusselator model of chemical reactions, for which no symmetry for the coupling matrices is assumed. We show that the model can undergo a Hopf bifurcation, and consequently periodic solutions can arise when the dispersal rates are large. Moreover, the effect of the coupling matrices on the Hopf bifurcation value is considered for a special case. [ABSTRACT FROM AUTHOR]

Given a family of nearly commuting symmetric matrices, we consider the task of computing an orthogonal matrix that nearly diagonalizes every matrix in the family. In this paper, we propose and analyze randomized joint diagonalization (RJD) for performing this task. RJD applies a standard eigenvalue solver to random linear combinations of the matrices. Unlike existing optimization-based methods, RJD is simple to implement and leverages existing high-quality linear algebra software packages. Our main novel contribution is to prove robust recovery: Given a family that is e-near to a commuting family, RJD jointly diagonalizes this family, with high probability, up to an error of norm O (epsilon). We also discuss how the algorithm can be further improved by deflation techniques and demonstrate its state-of-the-art performance by numerical experiments with synthetic and real-world data. [ABSTRACT FROM AUTHOR]

This paper presents a stability-preserving model reduction approach for a vibro-acoustic finite element model including poroelastic materials. Most of the research on these systems in the past was conducted in the frequency domain and there were less focus on the stability properties. However, with the increasing of interest in time-domain auralization and virtual sensing, stability-preserving model order reduction is becoming essential for efficient time-domain simulations. The original finite element models for such systems are already well established but not extended to the time domain because of its high computational demand. Therefore, this paper proposes a method to generate stable reduced-order models. The solid displacement–total displacement ( u s − u t ) formulation of the Biot's model is used to describe the poroelastic media. This formulation leads to a set of positive and symmetric system matrices with passive transfer functions. Applying the Kalman–Yakubovich–Popov lemma and the Cholesky/LDL decomposition, this formulation is modified to satisfy the previously proven stability-preserving conditions under one-sided model order reduction. Furthermore, the coupling conditions between the poroelastic media, air, and structure are also investigated. It is finally shown that the stability-preserving property is kept for the coupled models with appropriate variables. The proposed method is verified by several numerical simulations. • A time-domain formulation of Biot's equations is derived. • The time-domain Biot's model is stability-preserving under one-sided model order reduction methods. • The coupling models between the proposed time-domain Biot's model and vibro-acoustic systems are still stability-preserving. • The reduced order model enables efficient time-domain simulations of the vibro-acoustic system with the Biot's model. [ABSTRACT FROM AUTHOR]

A comprehensive approach to the spectrum characterization (derivation of eigenvalues and the corresponding multiplicities) for non-normalized, symmetric discrete trigonometric transforms (DTT) is presented in the paper. Eight types of the DTT are analyzed. New explicit analytic expressions for the eigenvalues, together with their multiplicities, for the cases of three DTT (DCT (1) , DCT (5) , and DST (8) ), are the main contribution of this paper. Moreover, the presented theory is supplemented by new, original derivations for the closed-form expressions of the square and the trace of all symmetric DTT matrices including existing results and the new ones. • New analytic expressions for the eigenvalues are derived for DCT1, DCT5, and DST • Closed-form expressions of the trace square of all symmetric DTT matrices are derived. [ABSTRACT FROM AUTHOR]

This paper addresses the problem of maximizing the effectiveness of the ranking produced by information retrieval or recommender systems and at the same time maximizing two fairnesses, that of the group and that of the individual. The context of this paper is therefore that of access to information carried out by users, who aim to satisfy their own information needs, to documents produced by authors and curators who aim to be exposed in a fair manner, i.e. without discriminating between groups nor individuals. The paper describes a general method based on the spectral decomposition of mixtures of symmetric matrices, each of which represents a variable to be maximized, and experiments conducted with the use of a test collection. The method described in this paper has explained if and how the trade-offs between effectiveness, group fairness and individual fairness manifest themselves. The experimental results show that maintaining an acceptable level of effectiveness and fairness at the same time is feasible and (b) the trade-offs exist but the order of magnitude of the variations depends on the measure of effectiveness used and therefore by what the user's model of access to information is as well as on the fairness measures and therefore on how authors or editors should be exposed. • Modern information access systems should balance fairness and effectiveness. • A single eigensystem achieves simultaneous maximization. • Fairness, effectiveness, and access measures are crucial in trade-offs. [ABSTRACT FROM AUTHOR]

In this paper we treat statistical inference for a wavelet estimator of curves of symmetric positive definite (SPD) using the log-Euclidean distance. This estimator preserves positive-definiteness and enjoys permutation-equivariance, which is particularly relevant for covariance matrices. Our second-generation wavelet estimator is based on average-interpolation (AI) and allows the same powerful properties, including fast algorithms, known from nonparametric curve estimation with wavelets in standard Euclidean set-ups. The core of our work is the proposition of confidence sets for our AI wavelet estimator in a non-Euclidean geometry. We derive asymptotic normality of this estimator, including explicit expressions of its asymptotic variance. This opens the door for constructing asymptotic confidence regions which we compare with our proposed bootstrap scheme for inference. Detailed numerical simulations confirm the appropriateness of our suggested inference schemes. • Confidence regions for curves in the space of symmetric positive-definite matrices • Based on second-generation wavelet estimators of the curves of SPD matrices • Asymptotic normality of estimators including derivation of their asymptotic variance • Wild bootstrap confidence regions shown to be valid via derived asymptotic normality [ABSTRACT FROM AUTHOR]

This paper investigates the synchronization control of neural networks (NNs) considering parameter uncertainties by utilizing an improved asymmetric Lyapunov–Krasovskii functional (ALKF) method. Firstly, an uncertain NNs model with discrete and distributed delays is developed. To synchronize the master–slave system, a memory sampled-data controller is designed. Secondly, an improved ALKF is constructed, in which the positive-definite and symmetric restrictions of matrices are relaxed. Furthermore, improved synchronization criteria are established by linear matrix inequalities (LMIs), which are based on the ALKF method and the integral inequality technique. Finally, two simulation examples are conducted to demonstrate the feasibility and superiority of the established stability conditions in reaching a maximum sampling period. [ABSTRACT FROM AUTHOR]

Optimal averaged Gauss quadrature rules provide estimates for the quadrature error in Gauss rules, as well as estimates for the error incurred when approximating matrix functionals of the form u T f (A) v with a large matrix A ∈ R N × N by low-rank approximations that are obtained by applying a few steps of the symmetric or nonsymmetric Lanczos processes to A ; here u , v ∈ R N are vectors. The latter process is used when the measure associated with the Gauss quadrature rule has support in the complex plane. The symmetric Lanczos process yields a real tridiagonal matrix, whose entries determine the recursion coefficients of the monic orthogonal polynomials associated with the measure, while the nonsymmetric Lanczos process determines a nonsymmetric tridiagonal matrix, whose entries are recursion coefficients for a pair of sets of bi-orthogonal polynomials. Recently, it has been shown, by applying the results of Peherstorfer, that optimal averaged Gauss quadrature rules, which are associated with a nonnegative measure with support on the real axis, can be expressed as a weighted sum of two quadrature rules. This decomposition allows faster evaluation of optimal averaged Gauss quadrature rules than the previously available representation. The present paper provides a new self-contained proof of this decomposition that is based on linear algebra techniques. Moreover, these techniques are generalized to determine a decomposition of the optimal averaged quadrature rules that are associated with the tridiagonal matrices determined by the nonsymmetric Lanczos process. Also, the splitting of complex symmetric tridiagonal matrices is discussed. The new splittings allow faster evaluation of optimal averaged Gauss quadrature rules than the previously available representations. Computational aspects are discussed. [ABSTRACT FROM AUTHOR]

A linear code with complementary dual (or an LCD code) is defined to be a linear code which intersects its dual code trivially. Let I be an identity matrix and T be a Toeplitz matrix of the same order over a finite field. A Double Toeplitz code (or a DT code) is a linear code generated by a generator matrix of the form (I , T). In 2021, Shi et al. obtained necessary and sufficient conditions for a Double Toeplitz code to be LCD when T is symmetric and tridiagonal. In this paper, by using a result on factoring Dickson polynomials over finite fields, we determine when a Double Toeplitz code is LCD for T being a skew symmetric and tridiagonal matrix. In addition, using a concatenation, we construct LCD codes with arbitrary minimum distance from DT codes over extension fields, provided the length of which is increased if necessary. [ABSTRACT FROM AUTHOR]

Multivariate time-series have become abundant in recent years, as many data-acquisition systems record information through multiple sensors simultaneously. In this paper, we assume the variables pertain to some geometry and present an operator-based approach for spatiotemporal analysis. Our approach combines three components that are often considered separately: (i) manifold learning for building operators representing the geometry of the variables, (ii) Riemannian geometry of symmetric positive-definite matrices for multiscale composition of operators corresponding to different time samples, and (iii) spectral analysis of the composite operators for extracting different dynamic modes. We propose a method that is analogous to the classical wavelet analysis, which we term Riemannian multi-resolution analysis (RMRA). We provide some theoretical results on the spectral analysis of the composite operators, and we demonstrate the proposed method on simulations and on real data. [ABSTRACT FROM AUTHOR]