Application of MFCAV Riemann Solver to Maire's Cell-Centered Lagrangian Method

In this paper, we apply arbitrary Riemann solvers, which may not satisfy the Maire's requirement, to the Maire's node-based Lagrangian scheme developed in [P. H. Maire et al., SIAM J. Sci. Comput, 29 (2007), 1781-1824]. In particular, we apply the so-called Multi-Fluid Channel on Averaged Volume (MFCAV) Riemann solver and a Riemann solver that adaptively combines the MFCAV solver with other more dissipative Riemann solvers to the Maire's scheme. It is noted that neither of the two solvers satisfies the Maire's requirement. Numerical experiments are presented to demonstrate that the application of the two Riemann solvers is successful.     

基于风险管理的动态供应链超网络均衡模型

基于风险管理的动态供应链超网络均衡模型的研究有助于供应链超网络节点厂商在动态环境下优化其风险管理,降低风险损失,提升供应链网络在风险管理下的竞争优势。本文以三层供应链超网络为研究对象,采用风险发生概率和损失函数表达供应链超网络中节点厂商中断风险的特征,构建了基于风险管理的动态供应链超网络均衡模型。模型中包括三种类型的节点,产品生产商、零售商和需求市场,生产商考虑风险损失的情况下,基于动态变化的风险、需求和成本追求个体期望效益最优化。接着,通过进化变分不等式来表达动态供应链超网络风险管理下的均衡解,并采用投影动态系统求解进化变分不等式,通过数值算例验证方法的可靠性和合理性,通过投影动态系统解释某一个厂商趋近均衡解的过程。通过单一厂商趋近均衡解的过程,阐述其他厂商相应的最优决策。     

Perturbations of Drazin invertible operators

The necessary and sufficient conditions for the small norm perturbation of a Drazin invertible operator to be still Drazin invertible and the sufficient conditions for the finite rank perturbation of a Drazin invertible operator to be still Drazin invertible are established.     

基于概率加权共轭梯度算法的混凝土超声波层析成像

在混凝土层析成像中,为了提高反演的准确性和计算效率,针对共轭梯度算法提出一种加权算法——概率加权共轭梯度算法.新算法不同于常规的加权算法,权重是加在成像单元上而不是方程上.为取得较好的权重因子和较好的迭代初始值,采用IART算法的权重和迭代初始值的选取方法.模拟算例和混凝土试验均表明这种加权算法的可行性.     

一种新的单幅条纹图的相位解调方法

基于非载波条纹图的条纹灰度分布的极值准确定位,提出了一种新的由单幅条纹图解调相位方法.该方法首先准确定位每个条纹灰度的极值包括最大和最小,再基于极值图将条纹的强度值分布线性地变换到-1和1之间,最后利用arcos反余弦取代传统的atan2反正切算法求出相位,结合极值图就能快速准确地解调出含有封闭条纹的条纹图相位.本方法仅用一幅条纹图,不需要传统的条纹定级和对分数级条纹插值与拟合,自动判别条纹相位变化的波谷,准确获得条纹相位分布.这一方法使光测法研究动态和瞬态问题成为可能,而不需引入载波法调制条纹.文中通过实例说明了本方法的处理过程.     

THE APPROXIMATION PROBLEM ON THE CLOSED CONVEX CONE AND ITS NUMERICAL SOLUTION

In this paper, we study the approximation problem on the closed convexcone, and prove that there exists a unique solution of the approximation problem, thengive the algorithm to compute the unique solution.     

A note on multilinear singular integrals with rough kernel

In this paper, the author gives the weighted weak Lipschitz boundedness with power weight for roughmultilinear integral operators. A simple way is obtained that is closely linked with a class of rough frac-tional integral operators.     

调制度分析在等步长相移法相位展开中的应用

推导了两种常用等步长相移算法的调制度表达式,提出一种新的调制度分析方法。该方法用于等步长相移法中基于加权最小二乘法的相位展开,能够充分利用调制度信息,构造二值和小数权重,从而增强相位展开过程对多种干扰因素的免疫力。实验结果说明了该方法的有效性和实用性。最后比较了二值权重和小数权重在加权最小二乘法的相位展开中的性能表现。     

Truncation of small matrix elements based on the Euclidean norm for blocked data structures

Methods for the removal of small symmetric matrix elements based on the Euclidean norm of the error matrix are presented in this article. In large scale Hartree-Fock and Kohn-Sham calculations it is important to be able to enforce matrix sparsity while keeping errors under control. Truncation based on some unitary-invariant norm allows for control of errors in the occupied subspace as described in (Rubensson et al. J Math Phys 49, 032103). The Euclidean norm is unitary-invariant and does not grow intrinsically with system size and is thus suitable for error control in large scale calculations. The presented truncation schemes repetitively use the Lanczos method to compute the Euclidean norms of the error matrix candidates. Ritz value convergence patterns are utilized to reduce the total number of Lanczos iterations.     

Maximal multilinear Calderón-Zygmund operators with non-doubling measures

Under the assumption that μ is a non-doubling measure on ℝ ^d which only satisfies some growth condition, the authors prove that the maximal multilinear Calderón-Zygmund operator is bounded from (μ) × … × (μ) into L ^p (μ) for any p ₁, … p _m ∈ (1, ∞) and p with 1/p = 1/p ₁ + … + 1/p _m , and bounded from (μ) × … × (μ) into weak- ^{L p} (μ) if there exists any p _i = 1. Furthermore, the authors establish a weighted weak-type estimate for the maximal multilinear Calderón-Zygmund operator. The second (corresponding) author is supported by Tianyuan Fund of China (No. 10626055) and National Natural Science Foundation of China (No. 10701078).