二重数值积分公式的高精度对偶公式及其应用

给出了计算二重积分的Simpson公式与两点高斯公式的对偶公式的构造过程,得到与之对应的高精度对偶修正解,提高了二重数值积分公式的计算精度,同时给出了二重积分的一种估值方法.最后,应用于几个典型的数值算例,计算结果表明：对偶修正解比对应的数值积分公式及其对偶公式的解有更高的计算精度和更快的收敛速度.     

奇异微分方程的广义Sturm--Liouville边值问题正解的存在性

By employing the fixed point theorem of cone expansion and compression of norm type, we investigate the existence of positive solutions of generalized Sturm-Liouville boundary value problems for a nonlinear singular differential equation with a parameter. Some sufficient conditions for the existence of positive solutions are established. In the last section, an example is presented to illustrate the applications of our main results.     

Physalis method for heterogeneous mixtures of dielectrics and conductors: Accurately simulating one million particles using a PC

Prosperetti’s seminal Physalis method, an Immersed Boundary/spectral method, had been used extensively to investigate fluid flows with suspended solid particles. Its underlying idea of creating a cage and using a spectral general analytical solution around a discontinuity in a surrounding field as a computational mechanism to enable the accommodation of physical and geometric discontinuities is a general concept, and can be applied to other problems of importance to physics, mechanics, and chemistry. In this paper we provide a foundation for the application of this approach to the determination of the distribution of electric charge in heterogeneous mixtures of dielectrics and conductors. The proposed Physalis method is remarkably accurate and efficient. In the method, a spectral analytical solution is used to tackle the discontinuity and thus the discontinuous boundary conditions at the interface of two media are satisfied exactly. Owing to the hybrid finite difference and spectral schemes, the method is spectrally accurate if the modes are not sufficiently resolved, while higher than second-order accurate if the modes are sufficiently resolved, for the solved potential field. Because of the features of the analytical solutions, the derivative quantities of importance, such as electric field, charge distribution, and force, have the same order of accuracy as the solved potential field during postprocessing. This is an important advantage of the Physalis method over other numerical methods involving interpolation, differentiation, and integration during postprocessing, which may significantly degrade the accuracy of the derivative quantities of importance. The analytical solutions enable the user to use relatively few mesh points to accurately represent the regions of discontinuity. In addition, the spectral convergence and a linear relationship between the cost of computer memory/computation and particle numbers results in a very efficient method. In the present paper, the accuracy of the method is numerically investigated by example computations using one dielectric particle, one isolated conductor particle, one conductor particle connected to an external source with imposed voltage, and two conductor/dielectric particles with strong interactions. The efficiency of the method is demonstrated with one million particles, which suggests that the method can be used for many important engineering applications of broad interest.     

宽量程高时间分辨扫描变像管

 为改善大偏转角度的偏转像质,设计了一种高时间分辨、大扫描长度的扫描变像管电子光学和行波偏转系统。对变像管后加速技术进行了实验研究,结合图像外增强技术,解决了大面积扫描图像的探测灵敏度问题,达到了同时扩展扫描变像管的时间量程和提高时间分辨力的目的,研制出了宽时间量程、高时间分辨力的扫描变像管。初步实验表明,在量程为5 ns时,时间分辨力达到5 ps,也就是时间分辨力与时间量程之比达到1/1 000。     

一种新的自动聚焦算法的研究

聚焦评价函数和极点搜索策略是自动聚焦的重要模块,直接影响聚焦的效果.为同时满足自动聚焦速度和精度上的要求,提出了一种新的自动聚焦算法.该算法综合了快速改进的灰度差分法和高精度小波变换法的优点,根据各阶段聚焦评价函数的梯度选择合适的评价函数和镜头移动步长,并且实现了从大步长快速粗扫到小步长精确细调的过渡.实验表明,与传统...     

标准镜的高精度柔性支撑镜框的优化设计

为了实现干涉仪标准镜中光学元件的高精度定位,设计了一种柔性支撑镜框,研究了该结构的力学模型、结构参数、定位精度和透镜变形。首先,根据材料力学原理将柔性镜框等效为一个弹簧系统;根据力学方程和几何关系,建立了透镜中心位置与柔性结构的挠度之间的二元方程。然后,分析了安装位置、温度、结构参数对透镜位置以及作用力的影响。最后,应...     

低振幅中心对称光折变暗和灰空间孤子的时间特性

为了研究中心对称光折变晶体中低振幅暗和灰空间孤子的时间特性,推导出了含时间参量的空间电荷场和光波的动态演化方程.采用数值方法,得到了中心对称光折变晶体中低振幅空间暗和灰孤子强度包络和强度半峰全宽的时间演化特性.结果表明,初始阶段形成的宽度较宽的孤子,其宽度随时间单调递减到一个最小值直至稳态孤子的形成;在相同的演化时间内...     

Efficient and stable perfectly matched layer for CEM

An efficient unsplit perfectly matched layer for numerical simulation of electromagnetic waves in unbounded domains is derived via a complex change of variables. In order to surround a Cartesian grid with the PML, the time-dependent PML requires only one (scalar) auxiliary variable in two space dimensions and six (scalar) auxiliary variables in three space dimensions. It is therefore cheap and straightforward to implement. We use Fourier and energy methods to prove the stability of the PML. We extend the stability result to a semi-discrete PML approximated by central finite differences of arbitrary order of accuracy and to a fully discrete problem for the ‘Leap-Frog’ schemes. This makes precise the usefulness of the derived PML model for longtime simulations. Numerical experiments are presented, illustrating the accuracy and stability of the PML.     

Probabilistic solution of random homogeneous linear second-order difference equations

This paper deals with the computation of the first probability density function of the solution of random homogeneous linear second-order difference equations by the Random Variable Transformation method. This approach allows us to generalize the classical solution obtained in the deterministic scenario. Several illustrative examples are provided.     

基于Beta函数及其偏导数的广义积分的高精度快速计算

考虑Beta函数偏导数的计算以及与此相关的广义积分的高精度快速计算问题.首先将Beta函数B(x,y)的定义扩展到整个复平面上,并建立了在整个复平面上Beta函数B(x,y)的偏导数的递推公式.对许多广义积分我们给出Beta函数偏导数的表示形式,因而利用Beta函数的偏导数计算这些广义积分.数值计算表明,算法无论从计算精度还是计算速度,远好于数值积分.另外,得到了B_(p,q)(x,y)存在闭形式的条件,并给出一些广义积分的闭形式.