高斯型空间光孤子相互作用的数值研究

 通过数值模拟的方法，对高斯孤子在对数型饱和非线性介质中的相互作用进行了研究，考查了两光束间的相对振幅和相对相位对其相互作用的影响。结果表明：高斯孤子之间的相互作用敏感地依赖于两光束间的相对振幅和相对相位。在不同的振幅差异范围内，光束间的主要作用交替地表现为相互排斥和相互吸引，并由于高斯孤子的不稳定性，导致了光束在碰撞后以一种尺寸周期性变化的呼吸模式传输。随着相对相位的增大，两光束间始终持续地表现出强烈的排斥作用，直到相对相位增加到一个2π周期之后。而且碰撞之后，光束也都以呼吸模式进行传输，其分离的角度越大，呼吸就越明显。     

一类相空间中的准几率分布函数系

定义了一类相空间中的准几率分布函数系，这个准几率分布函数系直接建立在具有更加广泛意义的量子相空间Schr?dinger方程解的基础之上，其中定义_α=αp-i?q和_α=(1-α)q+i?p.发现了两个有趣的关系.(1)建立的量子相空间Schr?dinger方程的解实际上是对函数φ(λ)exp［i(1-α)qp］做窗口Fourier变换.(2)这个窗口函数g(λ)起着选择窗口形式的作用，而且不同的窗口对应着不同的分布函数.当g(λ)是一个代表Gauss窗的Gauss函数的时候，准几率分布函数就是一个类似于Husimi的分布函数f^HL_α(q,p)；当g(λ)是一个表示椭圆的复函数时，准几率分布函数就是一个椭圆分布函数f^E_α(q,p)；再在g(λ)为复函数的基础上附加α=0，就可得到标准序分布函数f^S(q,p)、反标准序分布函数f^AS(q,p)和Wigner分布函数f^W(q,p)，此时g(λ)表示高度为1/12π?而长度为λ的矩形窗. 关键词： 窗口Fourier变换 相空间 Wigner分布函数     

Laguerre expansion on the Heisenberg group and Fourier-Bessel transform on C~n Dedicated to Professor Sheng GONG on the occasion of his 75th birthday

Given a principal value convolution on the Heisenberg group Hn = Cn×R, we study the relation between its Laguerre expansion and the Fourier-Bessel expansion of its limit on Cn. We also calculate the Dirichlet kernel for the Laguerre expansion on the group Hn.     

Semi‐analytical integration of the 8‐node plane element stiffness matrix using symbolic computation

The semi‐analytical integration of an 8‐node plane strain finite element stiffness matrix is presented in this work. The element is assumed to be super‐parametric, having straight sides. Before carrying out the integration, the integral expressions are classified into several groups, thus avoiding duplication of calculations. Symbolic manipulation and integration is used to obtain the basic formulae to evaluate the stiffness matrix. Then, the resulting expressions are postprocessed, optimized, and simplified in order to reduce the computation time. Maple symbolic‐manipulation software was used to generate the closed expressions and to develop the corresponding Fortran code. Comparisons between semi‐analytical integration and numerical integration were made. It was demonstrated that semi‐analytical integration required less CPU time than conventional numerical integration (using Gaussian‐Legendre quadrature) to obtain the stiffness matrix. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Numerical solution of the three‐dimensional parabolic equation with an integral condition

Developement of numerical methods for obtaining approximate solutions to the three dimensional diffusion equation with an integral condition will be carried out. The numerical techniques discussed are based on the fully explicit (1,7) finite difference technique and the fully implicit (7,1) finite difference method and the (7,7) Crank‐Nicolson type finite difference formula. The new developed methods are tested on a problem. Truncation error analysis and numerical examples are used to illustrate the accuracy of the new algorithms. The results of numerical testing show that the numerical methods based on the finite difference techniques discussed in the present article produce good results. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 193–202, 2002; DOI 10.1002/num.1040     

一类混合型微分差分方程的周期解

利用Fenchel变换,我们推出一类微分差分方程存在周期解等价于某泛函具有临界点,并求出方程具有周期解的充分条件.     

Construction algorithms for polynomial lattice rules for multivariate integration

We introduce a new construction algorithm for digital nets for integration in certain weighted tensor product Hilbert spaces. The first weighted Hilbert space we consider is based on Walsh functions. Dick and Pillichshammer calculated the worst-case error for integration using digital nets for this space. Here we extend this result to a special construction method for digital nets based on polynomials over finite fields. This result allows us to find polynomials which yield a small worst-case error by computer search. We prove an upper bound on the worst-case error for digital nets obtained by such a search algorithm which shows that the convergence rate is best possible and that strong tractability holds under some condition on the weights.

We extend the results for the weighted Hilbert space based on Walsh functions to weighted Sobolev spaces. In this case we use randomly digitally shifted digital nets. The construction principle is the same as before, only the worst-case error is slightly different. Again digital nets obtained from our search algorithm yield a worst-case error achieving the optimal rate of convergence and as before strong tractability holds under some condition on the weights. These results show that such a construction of digital nets yields the until now best known results of this kind and that our construction methods are comparable to the construction methods known for lattice rules.

We conclude the article with numerical results comparing the expected worst-case error for randomly digitally shifted digital nets with those for randomly shifted lattice rules.

     

On the counting function of the lattice profile of periodic sequences

The lattice profile analyzes the intrinsic structure of pseudorandom number sequences with applications in Monte Carlo methods and cryptology. In this paper, using the discrete Fourier transform for periodic sequences and the relation between the lattice profile and the linear complexity, we give general formulas for the expected value, variance, and counting function of the lattice profile of periodic sequences with fixed period. Moreover, we determine in a more explicit form the expected value, variance, and counting function of the lattice profile of periodic sequences for special values of the period.     

一种小波域音频信息隐藏方法

提出了一种基于量化的小波域音频隐藏算法，将保密语音隐藏到载体音频中．为提高隐藏重和保密语音传输的安全性，对保密语音进行了小波域压缩编码和m序列的扩频调制，生成待隐藏的比特序列；通过量化方法，将编码和调制后的保密语音隐藏到载体音频的小波系数中；保密语音的恢复过程不需要使用原始音频、仿真结果表明，隐藏有保密语音的载体音频听觉质量没有明显下降，提取的保密语音感知质量较好；该算法对重量化、加噪、低通滤波等攻击均有良好的鲁棒性．     

A parallel method for time discretization of parabolic equations based on Laplace transformation and quadrature

We consider the discretization in time of an inhomogeneous parabolicequation in a Banach space setting, using a representation ofthe solution as an integral along a smooth curve in the complexleft half-plane which, after transformation to a finite interval,is then evaluated to high accuracy by a quadrature rule. Thisreduces the problem to a finite set of elliptic equations withcomplex coefficients, which may be solved in parallel. The paperis a further development of earlier work by the authors, wherewe treated the homogeneous equation in a Hilbert space framework.Special attention is given here to the treatment of the forcingterm. The method is combined with finite-element discretizationin spatial variables.