Implicit Shape Reconstruction of Unorganized Points Using PDE-Based Deformable 3D Manifolds

<正>In this work we consider the problem of shape reconstruction from an unorganized data set which has many important applications in medical imaging,scientific computing,reverse engineering and geometric modelling.The reconstructed surface is obtained by continuously deforming an initial surface following the Partial Differential Equation(PDE)-based diffusion model derived by a minimal volume-like variational formulation.The evolution is driven both by the distance from the data set and by the curvature analytically computed by it.The distance function is computed by implicit local interpolants defined in terms of radial basis functions.Space discretization of the PDE model is obtained by finite co-volume schemes and semi-implicit approach is used in time/scale.The use of a level set method for the numerical computation of the surface reconstruction allows us to handle complex geometry and even changing topology, without the need of user-interaction.Numerical examples demonstrate the ability of the proposed method to produce high quality reconstructions.Moreover,we show the effectiveness of the new approach to solve hole filling problems and Boolean operations between different data sets.     

Multivariate quasi-tight framelets with high balancing orders derived from any compactly supported refinable vector functions

Generalizing wavelets by adding desired redundancy and flexibility,framelets(i.e.,wavelet frames)are of interest and importance in many applications such as image processing and numerical algorithms.Several key properties of framelets are high vanishing moments for sparse multiscale representation,fast framelet transforms for numerical efficiency,and redundancy for robustness.However,it is a challenging problem to study and construct multivariate nonseparable framelets,mainly due to their intrinsic connections to factorization and syzygy modules of multivariate polynomial matrices.Moreover,all the known multivariate tight framelets derived from spline refinable scalar functions have only one vanishing moment,and framelets derived from refinable vector functions are barely studied yet in the literature.In this paper,we circumvent the above difficulties through the approach of quasi-tight framelets,which behave almost identically to tight framelets.Employing the popular oblique extension principle(OEP),from an arbitrary compactly supported M-refinable vector functionφwith multiplicity greater than one,we prove that we can always derive fromφa compactly supported multivariate quasi-tight framelet such that:(i)all the framelet generators have the highest possible order of vanishing moments;(ii)its associated fast framelet transform has the highest balancing order and is compact.For a refinable scalar functionφ(i.e.,its multiplicity is one),the above item(ii)often cannot be achieved intrinsically but we show that we can always construct a compactly supported OEP-based multivariate quasi-tight framelet derived fromφsatisfying item(i).We point out that constructing OEP-based quasi-tight framelets is closely related to the generalized spectral factorization of Hermitian trigonometric polynomial matrices.Our proof is critically built on a newly developed result on the normal form of a matrix-valued filter,which is of interest and importance in itself for greatly facilitating the study of refinable vector functions and multiwavelets/multiframelets.This paper provides a comprehensive investigation on OEP-based multivariate quasi-tight multiframelets and their associated framelet transforms with high balancing orders.This deepens our theoretical understanding of multivariate quasi-tight multiframelets and their associated fast multiframelet transforms.     

Singular integral operators on product domains along twisted surfaces

We introduce a class of singular integral operators on product domains along twisted surfaces.We prove that the operators are bounded on L^p provided that the kernels satisfy weak conditions.     

Poisson Stable Solutions for Stochastic Differential Equations with Lévy Noise

In this paper,we use a unified framework to study Poisson stable(including stationary,periodic,quasi-periodic,almost periodic,almost automorphic,Birkhoff recurrent,almost recurrent in the sense of Bebutov,Levitan almost periodic,pseudo-periodic,pseudo-recurrent and Poisson stable)solutions for semilinear stochastic differential equations driven by infinite dimensional L′evy noise with large jumps.Under suitable conditions on drift,diffusion and jump coefficients,we prove that there exist solutions which inherit the Poisson stability of coefficients.Further we show that these solutions are globally asymptotically stable in square-mean sense.Finally,we illustrate our theoretical results by several examples.     

手机套餐问题的一个数学模型

随着信息时代的到来,手机在人们日常工作、社交、经营等社会活动中的作用越来越重要.近年来我国通信业务量飞速增长,手机的功劳更是功不可没.手机资费问题也越来越受到人们的关注,并且对原有的各种资费方案越来越质疑.2007年1月以来上海、北京、广东等地相继推出的手机"套餐"琳琅满目,让人眼花缭乱,人们不能理性分辨手机"套餐"究竟优惠在何处.……     

Kobayashi’s and Teichmuller’s Metrics and Bers Complex Manifold Structure on Circle Diffeomorphisms

Given a modulus of continuity ω,we consider the Teichmuller space TC1+ω as the space of all orientation-preserving circle diffeomorphisms whose derivatives are ω-continuous functions modulo the space of Mobius transformations preserving the unit disk.We study several distortion properties for diffeomorphisms and quasisymmetric homeomorphisms.Using these distortion properties,we give the Bers complex manifold structure on the Teichm(u| ")ller space TC^1+H as the union of over all0 <α≤1,which turns out to be the largest space in the Teichmuller space of C1 orientation-preserving circle diffeomorphisms on which we can assign such a structure.Furthermore,we prove that with the Bers complex manifold structure on TC^1+H ,Kobayashi’s metric and Teichmuller’s metric coincide.     

Difference of a Hauptmodul for Γ0(N ) and certain Gross-Zagier type CM value formulas

In this work, we show that the difference of a Hauptmodul for a genus zero group Γ₀(N) as a modular function on Y₀(N) × Y₀(N) is a Borcherds lift of type(2, 2). As applications, we derive the monster denominator formula like product expansions for these modular functions and certain Gross-Zagier type CM value formulas.     

Populations with individual variation in dispersal in heterogeneous environments: Dynamics and competition with simply diffusing populations

We consider a model for a population in a heterogeneous environment, with logistic-type local population dynamics, under the assumption that individuals can switch between two different nonzero rates of diffusion. Such switching behavior has been observed in some natural systems. We study how environmental heterogeneity and the rates of switching and diffusion affect the persistence of the population. The reactiondiffusion systems in the models can be cooperative at some population densities and competitive at others. The results extend our previous work on similar models in homogeneous environments. We also consider competition between two populations that are ecologically identical, but where one population diffuses at a fixed rate and the other switches between two different diffusion rates. The motivation for that is to gain insight into when switching might be advantageous versus diffusing at a fixed rate. This is a variation on the classical results for ecologically identical competitors with differing fixed diffusion rates, where it is well known that "the slower diffuser wins".     

A SPECTRAL METHOD FOR A WEAKLY SINGULAR VOLTERRA INTEGRO-DIFFERENTIAL EQUATION WITH PANTOGRAPH DELAY

In this paper,a Jacobi-collocation spectral method is developed for a Volterraintegro-differential equation with delay,which contains a weakly singular kernel.We use a function transformation and a variable transformation to change the equation into a new Volterra integral equation defined on the standard interval [-1,1],so that the Jacobi orthogonal polynomial theory can be applied conveniently.In order to obtain high order accuracy for the approximation,the integral term in the resulting equat...     

Solitary Wave Solutions of Some Nonlinear Physical Models Using Riccati Equation Approach

Riccati equation approach is used to look for exact travelling wave solutions of some nonlinear physical models.Solitary wave solutions are established for the modified KdV equation,the Boussinesq equation and the Zakharov-Kuznetsov equation.New generalized solitary wave solutions with some free parameters are derived.The obtained solutions,which includes some previously known solitary wave solutions and some new ones,are expressed by a composition of Riccati differential equation solutions followed by a polynomial.The employed approach,which is straightforward and concise,is expected to be further employed in obtaining new solitary wave solutions for nonlinear physical problems.     

基于Fuzzy等价关系的矩形聚类方法

证明了Boole矩阵与模糊矩阵等价的充分必要条件,总结出Fuzzy聚类的矩形方法.针对《中国大学评价》中的6所林业院校,应用矩形方法,直接从模糊关系矩阵逐步降低阈值,得到分类的矩形表.此方法操作简单,形象直观,分类合理.     

Maximum Likelihood Estimation of Hidden Markov Multivariate Normal Distribution Parameters

Hidden Markov model is widely used in statistical modeling of time, space and state transition data. The definition of hidden Markov multivariate normal distribution is given. The principle of using cluster analysis to determine the hidden state of observed variables is introduced. The maximum likelihood estimator of the unknown parameters in the model is derived. The simulated observation data set is used to test the estimation effect and stability of the method. The characteristic is simple classical statistical inference such as cluster analysis and maximum likelihood estimation. The method solves the parameter estimation problem of complex statistical models.     

Kolmogorov-Spieqel-Siveshinky方程大时间问题的Fourier拟谱逼近

该文讨论了Kolmogorov-Spieqel-Siveshinky方程的周期初值问题, 研究了半离散Fourier拟谱解的长时间行为, 证明了半离散系统的收敛性和整体吸引子的存在性. 构造了全离散的三层显式Fourier拟谱格式, 并证明了该格式的收敛性, 最后通过数值计算验证了格式的可信性. 数值结果表明: 该格式是长时间稳定并可取时间大步长. 作者模拟了方程的解在相空间的轨线, 得到了一些有意义的结论.     

Numerical solution of positive control problem via linear programming

The reachability problem for linear time-invariant discrete-time control systems with sign-restricted input is considered. The time-optimal control is constructed by an iterative procedure. Each step of the iteration is defined as a linear programming problem. This problem is solved by the simplex algorithm. The initial feasible solution for the simplex algorithm is provided by the preceding step of the iteration. The inversion of the basis matrix is reduced to a bordering procedure. The structural stability of the solution is investigated.     

Laplace方程有孔椭圆区域边值问题的解析解

本文利用相似流动替换方法 ,解决了中心有圆孔的椭园形区域上 Laplace方程第一类边值问题 ;采用分区域解法 ,给出了中心有椭园孔的椭园形区域上 Laplace方程第一类边值问题的解析通解 .这一结果在许多工程领域有重要应用 ,本文给出了油藏工程实例     

Exact solution of Marangoni convection in a binary fluid with throughflow and Soret effect

The onset of Marangoni convection in a non-reactive binary fluid layer in the presence of throughflow and Soret effect is determined. The bottom boundary of the fluid layer is assumed to be either conducting or insulating to temperature and solute concentration perturbations while the top boundary is free and insulating. The linear stability analysis is followed and an exact solution is obtained for the corresponding eigenvalue problem by assuming that stationary convection is exhibited at the neutral state. The contribution from the Soret effect is seen only when the throughflow is weak, but however for a wider range of upward throughflow when the bottom boundary is conducting. The instability gets advanced/delayed when the Soret parameter assumes negative/positive values. The results agree well with the existing results in the literature for some particular cases.     

B-Spline Collocation Method for Nonlinear Singularly-Perturbed Two-Point Boundary-Value Problems

A B-spline collocation method is presented for nonlinear singularly-perturbed boundary-value problems with mixed boundary conditions. The quasilinearization technique is used to linearize the original nonlinear singular perturbation problem into a sequence of linear singular perturbation problems. The B-spline collocation method on piecewise uniform mesh is derived for the linear case and is used to solve each linear singular perturbation problem obtained through quasilinearization. The fitted mesh technique is employed to generate a piecewise uniform mesh, condensed in the neighborhood of the boundary layers. The convergence analysis is given and the method is shown to have second-order uniform convergence. The stability of the B-spline collocation system is discussed. Numerical experiments are conducted to demonstrate the efficiency of the method.     

Patient-specific model of lung deformation using spatially dependent constitutive parameters

Breathing-induced spatially dependent lung deformation is predicted using patient-specific elastic properties with the contact–impact analysis model. The lung geometry is derived from 4D CT scan data of real patients. The spatially varying Young’s modulus for the patient is obtained from a previous study that used inverse deformation of the lung. The compact–impact analysis is implemented using the finite element method. The predicted lung deformation is compared with the results based on linear elasticity. The results are consistent with physiology, indicating large deformations near the diaphragm and smaller values at remote locations on the lobe. The effect of non-linearity of elastic property is most significant at the remote locations where the diaphragm-induced deformation is significantly attenuated.     

Stochastic optimal control of internal hierarchical labor markets

This paper develops an optimal control model for a graded manpower system where the demand for manpower is uncertain. The organization's objective is to minimize the discounted costs of operating the manpower system, including excess demand costs. The stock of workers in various grades can be adjusted in two ways. The first method is outside hiring flows, which is the usual control variable used in previous research. The second method is to control the transition rates between grades of the hierarchy, an instrument not previously studied. Incorporating the transition rates into the control variables creates time lags in the control process. The resulting problem is solved numerically using an approximation for the time-lagged control variables. The numerical example is based on the Air Force officer hierarchy. The model is used to examine such issues as the desirability of granting tenure to workers who are not promoted to the highest grade and the effects of length-of-service and demand uncertainty on manpower policy.     

Dynamic frictionless contact with adhesion

A model for the dynamic, adhesive, frictionless contact between a viscoelastic body and a deformable foundation is described. The adhesion process is modeled by a bonding field on the contact surface. The contact is described by a modified normal compliance condition. The tangential shear due to the bonding field is included. The problem is formulated as a coupled system of a variational equality for the displacements and a differential equation for the bonding field. The existence of a unique weak solution for the problem is established, together with a partial regularity result. The existence proof proceeds by construction of an appropriate mapping which is shown to be a contraction on a Hilbert space.