A necessary and sufficient condition is obtained for the generalized eigenfunction systems of 2 ×2 operator matrices to be a block Schauder basis of some Hilbert space, which offers a mathematical foundation of s...A necessary and sufficient condition is obtained for the generalized eigenfunction systems of 2 ×2 operator matrices to be a block Schauder basis of some Hilbert space, which offers a mathematical foundation of solving symplectic elasticity problems by using the method of separation of variables. Moreover, the theoretical result is applied to two plane elasticity problems via the separable Hamiltonian systems.展开更多
A block representation of the BLU factorization for block tridiagonal matrices is presented. Some properties on the factors obtained in the course of the factorization are studied. Simpler expressions for errors incur...A block representation of the BLU factorization for block tridiagonal matrices is presented. Some properties on the factors obtained in the course of the factorization are studied. Simpler expressions for errors incurred at the process of the factorization for block tridiagonal matrices are considered.展开更多
In this paper, we present a unified approach to decomposing a special class of block tridiagonal matrices <i>K</i> (<i>α</i> ,<i>β</i> ) into block diagonal matrices using similar...In this paper, we present a unified approach to decomposing a special class of block tridiagonal matrices <i>K</i> (<i>α</i> ,<i>β</i> ) into block diagonal matrices using similarity transformations. The matrices <i>K</i> (<i>α</i> ,<i>β</i> )∈ <i>R</i><sup><i>pq</i>× <i>pq</i></sup> are of the form <i>K</i> (<i>α</i> ,<i>β</i> = block-tridiag[<i>β B</i>,<i>A</i>,<i>α B</i>] for three special pairs of (<i>α</i> ,<i>β</i> ): <i>K</i> (1,1), <i>K</i> (1,2) and <i>K</i> (2,2) , where the matrices <i>A</i> and <i>B</i>, <i>A</i>, <i>B</i>∈ <i>R</i><sup><i>p</i>× <i>q</i></sup> , are general square matrices. The decomposed block diagonal matrices <img src="Edit_00717830-3b3b-4856-8ecd-a9db983fef19.png" width="15" height="15" alt="" />(<i>α</i> ,<i>β</i> ) for the three cases are all of the form: <img src="Edit_71ffcd27-6acc-4922-b5e2-f4be15b9b8dc.png" width="15" height="15" alt="" />(<i>α</i> ,<i>β</i> ) = <i>D</i><sub>1</sub> (<i>α</i> ,<i>β</i> ) ⊕ <i>D</i><sub>2</sub> (<i>α</i> ,<i>β</i> ) ⊕---⊕ <i>D</i><sub>q</sub> (<i>α</i> ,<i>β</i> ) , where <i>D<sub>k</sub></i> (<i>α</i> ,<i>β</i> ) = <i>A</i>+ 2cos ( <i>θ<sub>k</sub></i> (<i>α</i> ,<i>β</i> )) <i>B</i>, in which <i>θ<sub>k</sub></i> (<i>α</i> ,<i>β</i> ) , k = 1,2, --- q , depend on the values of <i>α</i> and <i>β</i>. Our decomposition method is closely related to the classical fast Poisson solver using Fourier analysis. Unlike the fast Poisson solver, our approach decomposes <i>K</i> (<i>α</i> ,<i>β</i> ) into <i>q</i> diagonal blocks, instead of <i>p</i> blocks. Furthermore, our proposed approach does not require matrices <i>A</i> and <i>B</i> to be symmetric and commute, and employs only the eigenvectors of the tridiagonal matrix <i>T</i> (<i>α</i> ,<i>β</i> ) = tridiag[<i>β b</i>, <i>a</i>,<i>αb</i>] in a block form, where <i>a</i> and <i>b</i> are scalars. The transformation matrices, their inverses, and the explicit form of the decomposed block diagonal matrices are derived in this paper. Numerical examples and experiments are also presented to demonstrate the validity and usefulness of the approach. Due to the decoupled nature of the decomposed matrices, this approach lends itself to parallel and distributed computations for solving both linear systems and eigenvalue problems using multiprocessors.展开更多
A fast algorithm FBTQ is presented which computes the QR factorization a block-Toeplitz matrix A (A∈R) in O(mns3) multiplications. We prove that the QR decomposition of A and the inverse Cholesky decomposition can be...A fast algorithm FBTQ is presented which computes the QR factorization a block-Toeplitz matrix A (A∈R) in O(mns3) multiplications. We prove that the QR decomposition of A and the inverse Cholesky decomposition can be computed in parallel using the sametransformation.We also prove that some kind of Toeplltz-block matrices can he transformed into the corresponding block-Toeplitz matrices.展开更多
A new collapse model of the trapdoors,three-dimensional rectangular trapdoor(3DRT),is presented for ground surface collapse.Undrained stability of 3DRT is examined with the upper bound method of plasticity limit analy...A new collapse model of the trapdoors,three-dimensional rectangular trapdoor(3DRT),is presented for ground surface collapse.Undrained stability of 3DRT is examined with the upper bound method of plasticity limit analysis theory.The soil where the trapdoors are located is assumed to be a perfectly plastic model with a Tresca yield criterion.Block analysis technique is employed to investigate the collapse of 3DRT.The model is divided into five different block types and added up to ten rigid blocks.According to the law of conservation of energy,the critical stability ratios of 3DRT are obtained through a search proceeding.The results of upper bound solution for 3DRT are given,and three trapdoor models with depth various are discussed during the application in the stability analysis of square trapdoors.The critical stability ratios can be used in the design of underground excavation and support force.展开更多
In this article, the computation of μ-values known as Structured Singular Values SSV for the companion matrices is presented. The comparison of lower bounds with the well-known MATLAB routine mussv is investigated. T...In this article, the computation of μ-values known as Structured Singular Values SSV for the companion matrices is presented. The comparison of lower bounds with the well-known MATLAB routine mussv is investigated. The Structured Singular Values provides important tools to analyze the stability and instability analysis of closed loop time invariant systems in the linear control theory as well as in structured eigenvalue perturbation theory.展开更多
Objective:Few studies have been conducted to establish animal models of left bundle branch block by using three-dimensional mapping systems.This research was aimed at creating a canine left bundle branch block model b...Objective:Few studies have been conducted to establish animal models of left bundle branch block by using three-dimensional mapping systems.This research was aimed at creating a canine left bundle branch block model by using a three-dimensional mapping system.Materials and Methods:We used a three-dimensional mapping system to map and ablate the left bundle branch in beagles.Results:Ten canines underwent radiofrequency ablation,among which left bundle branch block was successfully es-tablished in eight,one experienced ventricular fibrillation,and one developed third-degree atrioventricular block.The maximum HV interval measured within the left ventricle was 29.00±2.93 ms,and the LBP-V interval at the ablation site was 20.63±2.77 ms.The LBP-V interval at the ablation target was 71.08%of the maximum HV interval.Conclusion:This three-dimensional mapping system is a reliable and effective guide for ablation of the left bundle branch in dogs.展开更多
Rock joints infilled with sediments can strongly influence the strength of rock mass. As infilled joints often exist under unsaturated condition, this study investigated the influence of matric suction of infill on th...Rock joints infilled with sediments can strongly influence the strength of rock mass. As infilled joints often exist under unsaturated condition, this study investigated the influence of matric suction of infill on the overall joint shear strength. A novel technique that allows direct measurement of matric suction of infill using high capacity tensiometers(HCTs) during direct shear of infilled joints under constant normal stiffness(CNS) is described. The CNS apparatus was modified to accommodate the HCT and the procedure is explained in detail. Joint specimens were simulated by gypsum plaster using threedimensional(3D) printed surface moulds, and filled with kaolin and sand mixture prepared at different water contents. Shear behaviours of both planar infilled joints and rough joints having joint roughness coefficients(JRCs) of 8-10 and 18-20 with the ratios of infill thickness to asperity height(t/a)equal to 0.5 were investigated. Matric suction shows predominantly unimodal behaviour during shearing of both planar and rough joints, which is closely associated with the variation of unloading rate and volumetric changes of the infill material. As expected, two-peak behaviour was observed for the rough joints and both peaks increased with the increase of infill matric suction. The results suggest that the contribution of matric suction of infill on the joint peak normalised shear stress is relatively independent of the joint roughness.展开更多
Block matrices associated with discrete Trigonometric transforms （DTT＇s） arise in the mathematical modelling of several applications of wave propagation theory including discretizations of scatterers and radiators ...Block matrices associated with discrete Trigonometric transforms （DTT＇s） arise in the mathematical modelling of several applications of wave propagation theory including discretizations of scatterers and radiators with the Method of Moments, the Boundary Element Method, and the Method of Auxiliary Sources. The DTT＇s are represented by the Fourier, Hartley, Cosine, and Sine matrices, which are unitary and offer simultaneous diagonalizations of specific matrix algebras. The main tool for the investigation of the aforementioned wave applications is the efficient inversion of such types of block matrices. To this direction, in this paper we develop an efficient algorithm for the inversion of matrices with U-diagonalizable blocks （U a fixed unitary matrix） by utilizing the U- diagonalization of each block and subsequently a similarity transformation procedure. We determine the developed method＇s computational complexity and point out its high efficiency compared to standard inversion techniques. An implementation of the algorithm in Matlab is given. Several numerical results are presented demonstrating the CPU-time efficiency and accuracy for ill-conditioned matrices of the method. The investigated matrices stem from real-world wave propagation applications.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 11361034 and 11371185)the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20111501110001)the Natural Science Foundation of Inner Mongolia, China (Grant Nos. 2012MS0105 and 2013ZD01 )
文摘A necessary and sufficient condition is obtained for the generalized eigenfunction systems of 2 ×2 operator matrices to be a block Schauder basis of some Hilbert space, which offers a mathematical foundation of solving symplectic elasticity problems by using the method of separation of variables. Moreover, the theoretical result is applied to two plane elasticity problems via the separable Hamiltonian systems.
文摘A block representation of the BLU factorization for block tridiagonal matrices is presented. Some properties on the factors obtained in the course of the factorization are studied. Simpler expressions for errors incurred at the process of the factorization for block tridiagonal matrices are considered.
文摘In this paper, we present a unified approach to decomposing a special class of block tridiagonal matrices <i>K</i> (<i>α</i> ,<i>β</i> ) into block diagonal matrices using similarity transformations. The matrices <i>K</i> (<i>α</i> ,<i>β</i> )∈ <i>R</i><sup><i>pq</i>× <i>pq</i></sup> are of the form <i>K</i> (<i>α</i> ,<i>β</i> = block-tridiag[<i>β B</i>,<i>A</i>,<i>α B</i>] for three special pairs of (<i>α</i> ,<i>β</i> ): <i>K</i> (1,1), <i>K</i> (1,2) and <i>K</i> (2,2) , where the matrices <i>A</i> and <i>B</i>, <i>A</i>, <i>B</i>∈ <i>R</i><sup><i>p</i>× <i>q</i></sup> , are general square matrices. The decomposed block diagonal matrices <img src="Edit_00717830-3b3b-4856-8ecd-a9db983fef19.png" width="15" height="15" alt="" />(<i>α</i> ,<i>β</i> ) for the three cases are all of the form: <img src="Edit_71ffcd27-6acc-4922-b5e2-f4be15b9b8dc.png" width="15" height="15" alt="" />(<i>α</i> ,<i>β</i> ) = <i>D</i><sub>1</sub> (<i>α</i> ,<i>β</i> ) ⊕ <i>D</i><sub>2</sub> (<i>α</i> ,<i>β</i> ) ⊕---⊕ <i>D</i><sub>q</sub> (<i>α</i> ,<i>β</i> ) , where <i>D<sub>k</sub></i> (<i>α</i> ,<i>β</i> ) = <i>A</i>+ 2cos ( <i>θ<sub>k</sub></i> (<i>α</i> ,<i>β</i> )) <i>B</i>, in which <i>θ<sub>k</sub></i> (<i>α</i> ,<i>β</i> ) , k = 1,2, --- q , depend on the values of <i>α</i> and <i>β</i>. Our decomposition method is closely related to the classical fast Poisson solver using Fourier analysis. Unlike the fast Poisson solver, our approach decomposes <i>K</i> (<i>α</i> ,<i>β</i> ) into <i>q</i> diagonal blocks, instead of <i>p</i> blocks. Furthermore, our proposed approach does not require matrices <i>A</i> and <i>B</i> to be symmetric and commute, and employs only the eigenvectors of the tridiagonal matrix <i>T</i> (<i>α</i> ,<i>β</i> ) = tridiag[<i>β b</i>, <i>a</i>,<i>αb</i>] in a block form, where <i>a</i> and <i>b</i> are scalars. The transformation matrices, their inverses, and the explicit form of the decomposed block diagonal matrices are derived in this paper. Numerical examples and experiments are also presented to demonstrate the validity and usefulness of the approach. Due to the decoupled nature of the decomposed matrices, this approach lends itself to parallel and distributed computations for solving both linear systems and eigenvalue problems using multiprocessors.
文摘A fast algorithm FBTQ is presented which computes the QR factorization a block-Toeplitz matrix A (A∈R) in O(mns3) multiplications. We prove that the QR decomposition of A and the inverse Cholesky decomposition can be computed in parallel using the sametransformation.We also prove that some kind of Toeplltz-block matrices can he transformed into the corresponding block-Toeplitz matrices.
基金the Fundamental Research Funds for the Provincial Universities,China(No.702/000007020303)。
文摘A new collapse model of the trapdoors,three-dimensional rectangular trapdoor(3DRT),is presented for ground surface collapse.Undrained stability of 3DRT is examined with the upper bound method of plasticity limit analysis theory.The soil where the trapdoors are located is assumed to be a perfectly plastic model with a Tresca yield criterion.Block analysis technique is employed to investigate the collapse of 3DRT.The model is divided into five different block types and added up to ten rigid blocks.According to the law of conservation of energy,the critical stability ratios of 3DRT are obtained through a search proceeding.The results of upper bound solution for 3DRT are given,and three trapdoor models with depth various are discussed during the application in the stability analysis of square trapdoors.The critical stability ratios can be used in the design of underground excavation and support force.
文摘In this article, the computation of μ-values known as Structured Singular Values SSV for the companion matrices is presented. The comparison of lower bounds with the well-known MATLAB routine mussv is investigated. The Structured Singular Values provides important tools to analyze the stability and instability analysis of closed loop time invariant systems in the linear control theory as well as in structured eigenvalue perturbation theory.
基金This work was supported by the National Science Foundation for Young Researchers of China(grant Nos:82000315,82000325 and 82100325).
文摘Objective:Few studies have been conducted to establish animal models of left bundle branch block by using three-dimensional mapping systems.This research was aimed at creating a canine left bundle branch block model by using a three-dimensional mapping system.Materials and Methods:We used a three-dimensional mapping system to map and ablate the left bundle branch in beagles.Results:Ten canines underwent radiofrequency ablation,among which left bundle branch block was successfully es-tablished in eight,one experienced ventricular fibrillation,and one developed third-degree atrioventricular block.The maximum HV interval measured within the left ventricle was 29.00±2.93 ms,and the LBP-V interval at the ablation site was 20.63±2.77 ms.The LBP-V interval at the ablation target was 71.08%of the maximum HV interval.Conclusion:This three-dimensional mapping system is a reliable and effective guide for ablation of the left bundle branch in dogs.
基金The financial support provided by the China Scholarship Council (Grant No. 201406420027)
文摘Rock joints infilled with sediments can strongly influence the strength of rock mass. As infilled joints often exist under unsaturated condition, this study investigated the influence of matric suction of infill on the overall joint shear strength. A novel technique that allows direct measurement of matric suction of infill using high capacity tensiometers(HCTs) during direct shear of infilled joints under constant normal stiffness(CNS) is described. The CNS apparatus was modified to accommodate the HCT and the procedure is explained in detail. Joint specimens were simulated by gypsum plaster using threedimensional(3D) printed surface moulds, and filled with kaolin and sand mixture prepared at different water contents. Shear behaviours of both planar infilled joints and rough joints having joint roughness coefficients(JRCs) of 8-10 and 18-20 with the ratios of infill thickness to asperity height(t/a)equal to 0.5 were investigated. Matric suction shows predominantly unimodal behaviour during shearing of both planar and rough joints, which is closely associated with the variation of unloading rate and volumetric changes of the infill material. As expected, two-peak behaviour was observed for the rough joints and both peaks increased with the increase of infill matric suction. The results suggest that the contribution of matric suction of infill on the joint peak normalised shear stress is relatively independent of the joint roughness.
文摘Block matrices associated with discrete Trigonometric transforms （DTT＇s） arise in the mathematical modelling of several applications of wave propagation theory including discretizations of scatterers and radiators with the Method of Moments, the Boundary Element Method, and the Method of Auxiliary Sources. The DTT＇s are represented by the Fourier, Hartley, Cosine, and Sine matrices, which are unitary and offer simultaneous diagonalizations of specific matrix algebras. The main tool for the investigation of the aforementioned wave applications is the efficient inversion of such types of block matrices. To this direction, in this paper we develop an efficient algorithm for the inversion of matrices with U-diagonalizable blocks （U a fixed unitary matrix） by utilizing the U- diagonalization of each block and subsequently a similarity transformation procedure. We determine the developed method＇s computational complexity and point out its high efficiency compared to standard inversion techniques. An implementation of the algorithm in Matlab is given. Several numerical results are presented demonstrating the CPU-time efficiency and accuracy for ill-conditioned matrices of the method. The investigated matrices stem from real-world wave propagation applications.
