Structured condition numbers and statistical condition estimation for the LDU factorization

In this article, we consider the structured condition numbers for LDU, factorization by using the modified matrix-vector approach and the differential calculus, which can be represented by sets of parameters. By setting the specific norms and weight parameters, we present the expressions of the structured normwise, mixed, componentwise condition numbers and the corresponding results for unstructured ones. In addition, we investigate the statistical estimation of condition numbers of LDU factorization using the probabilistic spectral norm estimator and the small-sample statistical condition estimation method, and devise three algorithms. Finally, we compare the structured condition numbers with the corresponding unstructured ones in numerical experiments.

ON THE APPROXIMATE COMPUTATION OF EXTREME EIGENVALUES AND THE CONDITION NUMBER OF NONSINGULAR MATRICES

From the formulas of the conjugate gradient, a similarity between a symmetric positive definite (SPD) matrix A and a tridiagonal matrix B is obtained. The elements of the matrix B are determined by the parameters of the conjugate gradient. The computation of eigenvalues of A is then reduced to the case of the tridiagonal matrix B. The approximation of extreme eigenvalues of A can be obtained as a 'by-product' in the computation of the conjugate gradient if a computational cost of O(s) arithmetic operations is added, where s is the number of iterations This computational cost is negligible compared with the conjugate gradient. If the matrix A is not SPD, the approximation of the condition number of A can be obtained from the computation of the conjugate gradient on AT A. Numerical results show that this is a convenient and highly efficient method for computing extreme eigenvalues and the condition number of nonsingular matrices.

Geometric Condition of 3UPS-S Parallel Mechanism in Singular Configuration

The existing researches on singularity of parallel mechanism are mostly limited to the property and regularity of singularity locus and there is no further research into the geometric relationship between uncontrolled kinematic screw and parallel mechanism in singularity. A 3UPS-S parallel mechanism is presented which fulfils 3-DOF in rotation. The regularity of nutation angle singularity is analyzed based on the Jacobian matrix, and the singularity surface of 3UPS-S parallel mechanisms is obtained. By applying the concept of reciprocal product in screw theory, the singular kinematic screw is derived when 3UPS-S parallel mechanism is in singularity. The geometric relationship between singular kinematic screw and singular configuration of 3UPS-S parallel mechanism is investigated by using programs in MATLAB. It is revealed that there are two kinds of situation. Firstly, the three limbs of 3UPS-S parallel mechanism intersect the singular kinematic screw in space simultaneously; Secondly, two limbs cross the singular kinematic screw while the third limb parallels with that screw. It is concluded that the nutation angle singularity of 3UPS-S parallel mechanism belongs to the singular linear complexes. This paper sheds light into and clarifies the geometric relationship between singular kinematic screw and singular configuration of 3UPS-S parallel mechanism.

Novel method based on ant colony opti mization for solving ill-conditioned linear systems of equations

A novel method based on ant colony optimization （ACO）, algorithm for solving the ill-conditioned linear systems of equations is proposed. ACO is a parallelized bionic optimization algorithm which is inspired from the behavior of real ants. ACO algorithm is first introduced, a kind of positive feedback mechanism is adopted in ACO. Then, the solu- tion problem of linear systems of equations was reformulated as an unconstrained optimization problem for solution by an ACID algorithm. Finally, the ACID with other traditional methods is applied to solve a kind of multi-dimensional Hilbert ill-conditioned linear equations. The numerical results demonstrate that ACO is effective, robust and recommendable in solving ill-conditioned linear systems of equations.

On Strong Law of Large Numbers for Random Sequence

这笔记被奉献介绍有条件地统治的随机的变量的一个新概念。在下面合适限制条件，为任意的连续随机的变量的大数字的一条一般强壮的法律被获得。

On the Preconditioning Properties of RHSS Preconditioner for Saddle-Point Linear Systems

In this paper,for the regularized Hermitian and skew-Hermitian splitting(RHSS)preconditioner introduced by Bai and Benzi(BIT Numer Math 57:287–311,2017)for the solution of saddle-point linear systems,we analyze the spectral properties of the preconditioned matrix when the regularization matrix is a special Hermitian positive semidefinite matrix which depends on certain parameters.We accurately describe the numbers of eigenvalues clustered at(0,0)and(2,0),if the iteration parameter is close to 0.An estimate about the condition number of the corresponding eigenvector matrix,which partly determines the convergence rate of the RHSS-preconditioned Krylov subspace method,is also studied in this work.

2D Minimum Compliance Topology Optimization Based on a Region Partitioning Strategy

This paper presents an extended sequential element rejection and admission(SERA)topology optimizationmethod with a region partitioning strategy.Based on the partitioning of a design domain into solid regions and weak regions,the proposed optimizationmethod sequentially implements finite element analysis(FEA)in these regions.After standard FEA in the solid regions,the boundary displacement of the weak regions is constrained using the numerical solution of the solid regions as Dirichlet boundary conditions.This treatment can alleviate the negative effect of the material interpolation model of the topology optimization method in the weak regions,such as the condition number of the structural global stiffness matrix.For optimization,in which the forward problem requires nonlinear structural analysis,a linear solver can be applied in weak regions to avoid numerical singularities caused by the over-deformedmesh.To enhance the robustness of the proposedmethod,the nonmanifold point and island are identified and handled separately.The performance of the proposed method is verified by three 2D minimum compliance examples.

新型旋叶式汽车空调压缩机受力分析

本文针对为汽车空调系统而开发的新型旋叶式压缩机的结构特点，建立了受力分析的数学模型。并且在此基础上开发出了新型旋叶式汽车空调压缩机分析软件。应用这一软件，详细分析了新型旋叶式汽车空调压缩机中各种作用力的变化规律及其对压缩机性能的影响。同时详细分析了旋叶式压缩机中叶片数这一重要参数对压缩机受力的影响。

弱不连续问题高阶有限元离散系统的GAMG法

弱不连续问题(如含夹杂问题)是固体力学计算中的一类重要问题。高阶有限元方法由于其具有更好的逼近效果,是确保数值解在界面保持较高精度的计算方法之一。但与线性元相比,高阶单元需要更多的计算机存储单元,具有更高的计算复杂性。本文利用两水平算法的思想,将高阶有限元离散系统化归于线性元离散系统的求解,为弱不连续问题高阶有限元离散系统设计了一种新的基于几何与分析信息的代数多重网格(GAMG)法,并应用于圆形求解域含单夹杂问题的高阶有限元离散系统的求解。数值试验结果表明,相比于常用GAMG法,新方法的迭代次数基本不依赖于问题规模、单元阶次以及杨氏模量的间断性,CPU计算时间得到明显改善,具有更好的计算效率和鲁棒性,可大大提高弱不连续问题有限元分析的整体效率。

窗式空调器平行流冷凝器流程和扁管的优化分析

选择百叶窗肋片及小通道内流体传热流动的经验关联式,采用效能-传热单元数法建立了平行流冷凝器的稳态仿真模型。对不同流程数和扁管数方案的冷凝热和制冷剂压降进行了模拟分析,得出了扁管总数为39时的优化方案。

一类数列收敛的充分条件

分析了用于定义Euler-Mascheroni常数的数列，并在此基础上给出了其收敛的两个充分条件而非必要条件，这些条件较已有的条件更为简洁。

Novel method of improving the alignment accuracy of SINS on revolving mounting base

In the process of initial alignment for a strapdown inertial navigation system （SINS） on a stationary base, the east gyro drift rate is an important factor affecting the alignment accuracy of the azimuth misalignment angle. When the Kalman filtering algorithm is adopted in initial alignment, it yields a constant error in the estimation of the azimuth misalignment angle because the east gyro drift rate cannot be estimated. To improve the alignment accuracy, a novel alignment method on revolving mounting base is proposed. The Kalman filtering algorithm of extending the measured values is studied. The theory of spectral condition number is utilized to analyze the degrees of observability of states. Simulation results show that the estimation accuracy of the azimuth misalignment angle is greatly improved through revolving mounting base, and the proposed method is efficient in initial alignment for a medium accurate SINS.

5mm管房间空调器流路的优化设计

将目前房间空调系统中室外机和室内机常用的直径7mm及以上的内螺纹铜管改为5mm内螺纹管,可以达到节省材料与提高传热效率的效果。本文提出了空调系统换热管由7mm及以上铜管改为5mm铜管的设计方法和改进措施,如调整换热器纵向和横向换热管数目,制冷剂流路的分路数和压缩机更换等。通过一个分体式房间空调系统设计算例,具体说明将空调系统由7mm及以上管径改为5mm管径的具体的改进措施。设计结果表明在换热器具有相同换热性能的条件下,室外机由9mm管改为5mm管时,需同时增加换热器纵向和横向换热管数目,并调整分路数;室内机由7mm管改为5mm管室内机时,仅需要改变分路数。同时,采用高效压缩机,可进一步提高5mm管空调系统的性能。

Diffusion tensor imaging with multiple diffusion-weighted gradient directions

Diffusion tensor MRI （DT-MRI or DTI） is emerging as an important non-invasive technology for elucidating intemal brain structures. It has recently been utilized to diagnose a series of diseases that affect the integrity of neural systems to provide a basis for neuroregenerative studies. Results from the present study suggested that neural tissue is reconstructed with multiple diffusion-weighted gradient directions DTI, which varies from traditional imaging methods that utilize 6 gradient directions. Simultaneously, the diffusion tensor matrix is obtained by multiple linear regressions from an equation of echo signal intensity. The condition number value and standard deviation of fractional anisotropy for each scheme can be used to evaluate image quality. Results demonstrated that increasing gradient direction to some extent resulted in improved effects. Therefore, the traditional 6 and 15 directions should not be considered optimal scan protocols for clinical DTI application. In a scheme with 20 directions, the condition number and standard deviation of fractional anisotropy of the encoding gradients matrix were significantly reduced, and resulted in more clearly and accurately displayed neural tissue. Results demonstrated that the scheme with 20 diffusion gradient directions provided better accuracy of structural renderings and could be an optimal scan protocol for clinical DTI application.

A complete solution of an improved universal 3D coordinate similarity transformation model

When linearizing three-dimensional(3 D)coordinate similarity transformation model with large rotations,we usually encounter the ill-posed normal matrix which may aggravate the instability of solutions.To alleviate the problem,a series of conversions are contributed to the 3 D coordinate similarity transformation model in this paper.We deduced a complete solution for the 3 D coordinate similarity transformation at any rotation with the nonlinear adjustment methodology,which involves the errors of the common and the non-common points.Furthermore,as the large condition number of the normal matrix resulted in an intractable form,we introduced the bary-centralization technique and a surrogate process for deterministic element of the normal matrix,and proved its benefit for alleviating the condition number.The experimental results show that our approach can obtain the smaller condition number to stabilize the convergence of the interested parameters.Especially,our approach can be implemented for considering the errors of the common and the non-common points,thus the accuracy of the transformed coordinates improves.

Domain Decomposition for Wavelet Single Layer on Geometries with Patches

We focus on the single layer formulation which provides an integral equation of the first kind that is very badly conditioned. The condition number of the unpreconditioned system increases exponentially with the multiscale levels. A remedy utilizing overlapping domain decompositions applied to the Boundary Element Method by means of wavelets is examined. The width of the overlapping of the subdomains plays an important role in the estimation of the eigenvalues as well as the condition number of the additive domain decomposition operator. We examine the convergence analysis of the domain decomposition method which depends on the wavelet levels and on the size of the subdomain overlaps. Our theoretical results related to the additive Schwarz method are corroborated by numerical outputs.

Stability of GM(1,1) power model on vector transformation

The morbidity problem of the GM（1,1） power model in parameter identification is discussed by using multiple and rotation transformation of vectors. Firstly we consider the morbidity problem of the special matrix and prove that the condition number of the coefficient matrix is determined by the ratio of lengths and the included angle of the column vector, which could be adjusted by multiple and rotation transformation to turn the matrix to a well-conditioned one. Then partition the corresponding matrix of the GM（1,1） power model in accordance with the column vector and regulate the matrix to a well-conditioned one by multiple and rotation transformation of vectors, which completely solve the instability problem of the GM（1,1） power model. Numerical results show that vector transformation is a new method in studying the stability problem of the GM（1,1） power model.

Novel high-performance element in the electromagnetic finite-element method——node-edge element

It is known in the computational electromagnetics （CEM） that the node element has a relative wellconditioned matrix, but suffers from the spurious solution problem; whereas the edge element has no spurious solutions, but usually produces an ill-conditioned matrix. Particularly, when the mesh is over dense, the iterative solution of the matrix equation from edge element converges very slowly. Based on the node element and edge element, a node-edge element is presented, which has no spurious solutions and better-conditioned matrix. Numerical experiments demonstrate that the proposed node-edge element is more efficient than now-widely used edge element.

Capacity analysis of high-rank line-of-sight MIMO channels

Matrix perturbation theory is utilized to investigate high-rank line of sight multiple input multiple output channels in a microwave relay system. The upper and lower bounds of channel capacity are derived based on space time block codes technique and singular values decomposition. A useful constraint for designing LOS MIMO channels is developed by the use of the condition number of the MIMO channel matrix. The theoretical analysis of channel capacity is confirmed by the simulation. The results show that the proposed method is able to give a physical explanation of the high-rank LOS MIMO channel matrix characteristics.

ON THE SENSITIVITY OF THE LDU FACTORIZATION

In this paper,we give the sensitivity analyses by two approaches for L,D,U in factorization A=LDU of general perturbations in A which sufficiently small in norm. By the matrix vector equation approach,we derive the sharp condition number for L,D and U factors .By the matrix equation approach we derive corresponding condition estimates. When A is a symmetric matrix,the corresponding results can be obtained for LDL T factorization.