Algorithms for solving inverse problems in the optics of layered media based on comparing the extrema of spectral characteristics

A new algorithm for determining the optical parameters of deposited multilayer optical coatings based on comparing the positions of the extrema of the optical characteristics of multilayer optical coatings is proposed. Two versions of this algorithm are compared. Using a series of numerical simulation experiments, the advantage of one of these versions is demonstrated. It is shown that this version decreases the influence of systematic errors in the spectral data.     

Numerical solution of the problem of image formation in two planes using the integral fresnel formula

The paper considers numerical solution of the problem of constructing a plane phase optical element which forms prespecified images in two focal planes when illuminated by a plane wave front. The propagation of radiation from the optical element to the focal planes obeys the integral Fresnel formula. The phase function defining the optical element, is from the space of feasible phase functions satisfying both computational and technological constraints. The computer algorithm is based on efficient numerical solution of the inverse directional diagram problem. The phase function is constructed using the projector on the space of feasible phase functions. The efficiency of the phase optical element produced by the proposed procedure is demonstrated by solving a prototype problem involving simulation of the direct operator. The corresponding optical element was manufactured by electron-beam lithography. Physical experiments with a laser source illuminating the element confirmed the findings of the numerical experiment. Translated from Prikladnaya Matematika i Informatika, No. 1, pp. 86–89, 1999.     

Second kind integral equation formulation for the mode calculation of optical waveguides

We present a second kind integral equation (SKIE) formulation for calculating the electromagnetic modes of optical waveguides, where the unknowns are only on material interfaces. The resulting numerical algorithm can handle optical waveguides with a large number of inclusions of arbitrary irregular cross section. It is capable of finding the bound, leaky, and complex modes for optical fibers and waveguides including photonic crystal fibers (PCF), dielectric fibers and waveguides. Most importantly, the formulation is well conditioned even in the case of nonsmooth geometries. Our method is highly accurate and thus can be used to calculate the propagation loss of the electromagnetic modes accurately, which provides the photonics industry a reliable tool for the design of more compact and efficient photonic devices. We illustrate and validate the performance of our method through extensive numerical studies and by comparison with semi-analytical results and previously published results.     

A unified framework for constructing globally convergent algorithms for multidimensional coefficient inverse problems

We present a unified framework for constructing the globally convergent algorithms for a broad class of multidimensional coefficient inverse problems arising in natural science and industry. Based on the convexification approach, the unified framework substantiates the numerical solution of the corresponding problem of nonconvex optimization. A globally convergent iterative algorithm for an inverse problem of diffuse optical mammography is constructed. It utilizes the contraction property of a nonlinear operator resulting from applying the convexification approach. The effectiveness of this algorithm is demonstrated in computational experiments.     

一种寻找星图光斑质心的高精度算法及在成像畸变矫正中的应用

首先讨论星图识别中的质心寻找问题,提出了一种改进的匹配算法.使用高斯分布对光斑内像素点的强度分布进行近似模拟, 并采用最小二乘法中的法方程法对光斑图像质心坐标进行优化计算.算例表明该算法可以得到更为精确的结果, 在噪声小于最大像素值的$5\%$时,质心坐标计算精度高于0.04个像素, 优于以往的匹配算法的结果.最后讨论了算法在望远镜光学成像畸变矫正中的应用, 设计了矫正算法, 并通过算例验证了算法的合理性.     

A Staged Continuous Tabu Search Algorithm for the Global Optimization and its Applications to the Design of Fiber Bragg Gratings

A novel staged continuous Tabu search (SCTS) algorithm is proposed for solving global optimization problems of multi-minima functions with multi-variables. The proposed method comprises three stages that are based on the continuous Tabu search (CTS) algorithm with different neighbor-search strategies, with each devoting to one task. The method searches for the global optimum thoroughly and efficiently over the space of solutions compared to a single process of CTS. The effectiveness of the proposed SCTS algorithm is evaluated using a set of benchmark multimodal functions whose global and local minima are known. The numerical test results obtained indicate that the proposed method is more efficient than an improved genetic algorithm published previously. The method is also applied to the optimization of fiber grating design for optical communication systems. Compared with two other well-known algorithms, namely, genetic algorithm (GA) and simulated annealing (SA), the proposed method performs better in the optimization of the fiber grating design.     

Missing Observations in Spatial Models: A Spectral EM Algorithm

This article presents a model-based thin-plate smoothing method for optimal signal extraction and interpolation of missing data in spatial datasets. The method is based on a spectral EM algorithm where the two steps can be carried out in the frequency domain. In essence, the approach allows both dimensions to be treated separately from each other effectively rendering a likelihood that is easy to evaluate. As a result the algorithm is computationally inexpensive, in terms of both memory size and computing time, while allowing us to obtain an analytic expression for the asymptotic variance of the signal-to-noise ratio with which to construct confidence intervals of the missing data. Some numerical Monte Carlo simulations and a real data example using remotely sensed global aerosol optical thickness data illustrate the results given. Supplemental materials (Matlab computer code and dataset) are available online.     

基于绳索作动器的大型太空望远镜桁架结构的振动主动控制

薄膜衍射是一种新型的太空望远镜的成像方式,它具有轻质、易折叠与展开、光学成像精度高等许多优点,是当今太空望远镜技术的研究热点.该文针对一类薄膜衍射太空望远镜桁架结构的振动主动控制进行了研究,提出了一种基于绳索作动器的振动主动控制策略.首先建立了望远镜桁架结构的动力学模型,然后采用粒子群优化算法研究了绳索作动器的优化布置...     

Continuation and stability analysis for Bloch waves of the Gross-Pitaevskii equation

We describe a two-parameter continuation algorithm for computing Bloch waves of Bose-Einstein condensates (BEC) in optical lattices which is governed by the Gross-Pitaevskii equation (GPE). The Fourier collocation method and fourth-order Adini’s elements with penalty are used to discretize the GPE. We propose two different approaches so that the two-parameter continuation algorithm can be modified to compute closed tubes at the four corners of the Bloch band. We also study linear stability analysis for the GPE. We show that all the discrete steady-state solutions are numerically neutrally stable. Numerical results show that the four edges of the Bloch waves are surrounded by closed loops if the coefficient of the cubic nonlinear term is greater than that of the periodic potential. Moreover, closed tubes at the four corners of the Bloch band are obtained. The numerical results display superfluidity of BEC.     

Stabilization of the motions of mechanical systems in prescribed phase-space manifolds

A method for constructing a mathematical model of the dynamics of a mechanical system is proposed. An algorithm is constructed for determining the expressions for the control forces and the components of the constraint reactions. A modification is made to the dynamic equations which enables one to solve the problem of stabilizing the constraints and which ensures the required accuracy in the numerical solution of the corresponding system of differential-algebraic equations describing the constraints imposed on the system, its kinematics and dynamics. By virtue of well-known dynamic analogies, the proposed method can be used to investigate the dynamics of different physical systems. The problem of modelling the dynamics of an adaptive optical system with two degrees of freedom is considered.     

An alternative to the Bathe algorithm

This paper presents a new composite sub-steps algorithm for solving reliable numerical responses in structural dynamics. The newly developed algorithm is a two sub-steps, second-order accurate and unconditionally stable implicit algorithm with the same numerical properties as the Bathe algorithm. The detailed analysis of the stability and numerical accuracy is presented for the new algorithm, which shows that its numerical characteristics are identical to those of the Bathe algorithm. Hence, the new sub-steps scheme could be considered as an alternative to the Bathe algorithm. Meanwhile, the new algorithm possesses the following properties: (a) it produces the same accurate solutions as the Bathe algorithm for solving linear and nonlinear problems; (b) it does not involve any artificial parameters and additional variables, such as the Lagrange multipliers; (c) The identical effective stiffness matrices can be obtained inside two sub-steps; (d) it is a self-starting algorithm. Some numerical experiments are given to show the superiority of the new algorithm and the Bathe algorithm over the dissipative CH-α algorithm and the non-dissipative trapezoidal rule.     

L~1空间中第二类Fredholm积分方程数值解法探究

在L¹空间中对第二类Fredholm积分方程提出了一种求其数值解的算法,证明了算法的收敛性,并给出相应的误差估计.数值算例进一步验证了算法的合理性.     

A smoothing-type algorithm for solving nonlinear complementarity problems with a non-monotone line search

The smoothing-type algorithm has been successfully applied to solve various optimization problems. In general, the smoothing-type algorithm is designed based on some monotone line search. However, in order to achieve better numerical results, the non-monotone line search technique has been used in the numerical computations of some smoothing-type algorithms. In this paper, we propose a smoothing-type algorithm for solving the nonlinear complementarity problem with a non-monotone line search. We show that the proposed algorithm is globally and locally superlinearly convergent under suitable assumptions. The preliminary numerical results are also reported.     

Multiscale computational method for nonstationary integrated heat transfer problem in periodic porous materials

This article discusses multiscale analysis and numerical algorithm for the nonstationary integrated heat transfer problem with rapidly oscillating coefficients. The multiscale asymptotic expansion of the solution for this kind of problems is presented first. Then, error estimates of the multiscale approximate solution are derived, and a numerical algorithm based on the multiscale method for temperature field is introduced. Finally, using some numerical models, we verify the validity and relevancy of the proposed algorithm. The numerical results show that the algorithm is effective to predict the heat transfer performance of porous materials, and support the convergence theorem reported in this article. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 510–530, 2016     

A residual algorithm for finding a fixed point of a nonexpansive mapping

This paper presents a weak convergence residual algorithm for finding a fixed point of a nonexpansive mapping in a real Hilbert space. To study the numerical behavior of the algorithm it is included an extensive series of numerical experiments. Our computational experiments show that the new algorithm is computationally efficient.     

An integrated solution and analysis of bioluminescence tomography and diffuse optical tomography

While diffuse optical tomography (DOT) has been studied for years, bioluminescence tomography (BLT) is emerging as a promising optical molecular imaging tool. These two modalities have different goals. DOT is for reconstruction of optical parameters of a medium such as a breast from surface measurements induced by external sources. BLT is for reconstruction of a bioluminescent source distribution in a medium such as a mouse from surface measurements induced by internal bioluminescent sources. However, an important pre-requisite for BLT reconstruction is the knowledge on the distribution of optical parameters within the medium, which is the output of DOT. In this paper, we propose a mathematical model integrating BLT and DOT at the fundamental level; that is, performing the two types of reconstructions simultaneously instead of doing them sequentially. The model is introduced through minimizing the difference between predicted quantities and boundary measurements, as well as incorporating regularization terms. Then, we show the solution existence, introduce numerical schemes and prove convergence of the numerical solution. We also present numerical results to illustrate the utility of our approach.     

A reliable algorithm for solving tenth-order boundary value problems

In this paper we present an efficient numerical algorithm for solving linear and nonlinear boundary value problems with two-point boundary conditions of tenth-order. The differential transform method is applied to construct the numerical solutions. The proposed algorithm avoids the complexity provided by other numerical approaches. Several illustrative examples are given to demonstrate the effectiveness of the present algorithm.      

A Dynamical Method for Solving the Obstacle Problem

In this paper, we consider the unilateral obstacle problem, trying to find the numerical solution and coincidence set. We construct an equivalent format of the original problem and propose a method with a second-order in time dissipative system for solving the equivalent format. Several numerical examples are given to illustrate the effectiveness and stability of the proposed algorithm. Convergence speed comparisons with existent numerical algorithm are also provided and our algorithm is fast.     

Bakhvalov-Shishkin网格上求解边界层问题的差分进化算法

该文在Bakhvalov-Shishkin网格上求解具有左边界层或右边界层的对流扩散方程,并采用差分进化算法对Bakhvalov-Shishkin网格中的参数进行优化,获得了该网格上具有最优精度的数值解.对三个算例进行了数值模拟,数值结果表明:采用差分进化算法求解具有较高的计算精度和收敛性,特别是边界层的数值解精度明显...     

API stepsize control for the numerical solution of ordinary differential equations

A control-theoretic approach is used to design a new automatic stepsize control algorithm for the numerical integration of ODE's. The new control algorithm is more robust at little extra expense. Its improved performance is particularly evident when the stepsize is limited by numerical stability. Comparative numerical tests are presented.