The complex Jacobi iterative method for non-paraxial beam propagation in nonlinear optical waveguides

The recently introduced beam propagation method using complex Jacobi iteration adapted for modeling of non-paraxial beam propagation in nonlinear optical waveguides is presented in this paper. The beam propagation equation is based on our recently proposed modified Padé(1,1) approximant operator. The resulting approach is very efficient and well-suited for large structures with long propagation paths.     

A stable complex Jacobi iterative solution of 3D semivectorial wide-angle beam propagation using the iterated Crank–Nicholson method

An extension of the recently proposed three-dimensional (3D) wide-angle (WA) beam propagation method (BPM) using complex Jacobi iteration (CJI) taken into account polarization effects is presented. The resulting iterative BPM is faster than BPMs based on the traditional direct matrix inversion for waveguides with unchanging refractive index profiles during propagation direction. However, for varying refractive index waveguides the iterative method suffered from the fact that the iteration count between two successive cross-sections increases dramatically during the propagation direction. To overcome this problem, we propose the utility of the iterated Crank–Nicholson method. At each propagation step, the propagation equation is divided in multiple stages by the iterated Crank–Nicholson method and then each stage is recast in terms of a Helmholtz equation with source term, which is solved effectively by the complex Jacobi iterative method. The resulting approach is stable and well-suited for large structures with long propagation paths.     

迭代法重建平板光波导折射率分布

提出了一种用于重建渐变折射率平板光波导折射率分布的方法,这种方法通过数值法求解亥姆霍兹方程．从光波导的低阶模开始采用迭代法顺次对各个模式对应的“转折点”位置进行校正．并多次重复迭代以消除“转折点”外折射率分布的变化对相应模式“转折点”位置的影响。借助WKB近似原理说明了这种方法迭代过程中“转折点”坐标的收敛性。对折射率按照指数规律变化的平板光波导的折射率分布重建结果证明这种方法具有较高的精度。分析表明：这种方法克服了WKB法固有的缺陷,即忽略模式“转折点”外折射率分布对相应模式传输常鞋的影响和“转折点”处相移的不确定性,因而比逆WKB法具有更好的自洽性和精确性。     

A spectral collocation algorithm for two-point boundary value problem in fiber Raman amplifier equations

A novel algorithm implementing Chebyshev spectral collocation (pseudospectral) method in combination with Newton’s method is proposed for the nonlinear two-point boundary value problem (BVP) arising in solving propagation equations in fiber Raman amplifier. Moreover, an algorithm to train the known linear solution for use as a starting solution for the Newton iteration is proposed and successfully implemented. The exponential accuracy obtained by the proposed Chebyshev pseudospectral method is demonstrated on a case of the Raman propagation equations with strong nonlinearities. This is in contrast to algebraic accuracy obtained by typical solvers used in the literature. The resolving power and the efficiency of the underlying Chebyshev grid are demonstrated in comparison to a known BVP solver.     

任意色散周期结构的时域有限差分方法分析

利用辅助微分时域有限差分法求解了任意色散周期模型的电磁波传播问题。利用共轭复数对形式对任意色散媒质进行参数拟合,并将任意色散媒质的介电常数表示成公式形式,在FDTD迭代式中引入辅助微分方程,推导出了适用于多层任意色散模型的通用递推公式,分别求解了Debye、Drude与太阳能电池周期结构模型的电磁特性仿真问题。仿真结果表明:数值计算结果与CST商业软件仿真结果基本吻合,证明了所构建方法的有效性与普适性。     

Chirp-coefficient bisection iteration method for phase-intensity reconstruction of chirped pulses

A new method is proposed to systematically measure the phase-intensity information of chirped pulses that is based on the chirp-coefficient bisection iteration (CBI) concept. Through the CBI procedure with measured spectrum and temporal intensity profiles (or intensity autocorrelations), spectral and temporal amplitude–phase information can be rapidly retrieved. We experimentally verified that our method has high precision for nanosecond- and picosecond-level pulses and low precision for femtosecond-level pulses. Our proposed method does not require a sophisticated setup and has the advantage of accurate determination of temporal and spectral chirp coefficients with various orders. It also has lower cost, simple operation, in particular covers a wider measurement range than the main current methods. Moreover, the retrieved waveforms can reveal both the pulse shape and the actual intensity with spectral and temporal chirped coefficients of various orders, which can be directly used in various pulse propagation analyses such as chirped pulse amplification.     

The extended Fourier pseudospectral time-domain method for atmospheric sound propagation

An extended Fourier pseudospectral time-domain (PSTD) method is presented to model atmospheric sound propagation by solving the linearized Euler equations. In this method, evaluation of spatial derivatives is based on an eigenfunction expansion. Evaluation on a spatial grid requires only two spatial points per wavelength. Time iteration is done using a low-storage optimized six-stage Runge-Kutta method. This method is applied to two-dimensional non-moving media models, one with screens and one for an urban canyon, with generally high accuracy in both amplitude and phase. For a moving atmosphere, accurate results have been obtained in models with both a uniform and a logarithmic wind velocity profile over a rigid ground surface and in the presence of a screen. The method has also been validated for three-dimensional sound propagation over a screen. For that application, the developed method is in the order of 100 times faster than the second-order-accurate FDTD solution to the linearized Euler equations. The method is found to be well suited for atmospheric sound propagation simulations where effects of complex meteorology and straight rigid boundary surfaces are to be investigated.     

用串行迭代算法设计具有长焦深功能的衍射光学元件

提出一种串行迭代算法,将它用于控制波前传播进行衍射光学元件设计.用此算法设计出能使半导体激光束聚焦并具有长焦深功能的衍射光学元件.模拟计算结果表明了此设计方法的有效性.该衍射光学元件已用光刻和离子刻蚀技术制备并进行了光学实验研究,取得了与理论设计一致的实验结果.     

Wave propagation in a random medium: The solution of the m-nth moment equation with different wavenumbers

The moment equation with different wavenumbers and different transverse coordinates for wave propagation in a random medium is a linear differential equation. It often appears in the study of problems related to wave propagation in a random medium. The differential equation can be converted into an integral equation by using Green's functions and the integral equation can be solved by iteration. The moment equation is solved by the method of successive scatters, too. The solution of the moment equation is a Dyson expansion. The physical implication of the successive solution of the moment equation with different wavenumbers is explained.     

A spatial domain constraint iterative phase retrieval method for the propagation based X ray phase contrast imaging

Based on the propagation X ray phase contrast imaging theory, a spatial domain constraint iterative phase retrieval method is described in detail. This algorithm limits the object spatial domain according to the actual sample size firstly, and modify the image plane data with the actual test data, then the iteration can be terminated until iteration precision or the number of iterations meet the preset requirements. Finally the numerical simulation is made to evaluate the rapid phase retrieval algorithm performance, and a real column fiber material experiment is carried out using a micro focus X ray phase contrast imaging experiment platform, the phase distribution image of the column fiber is calculated out by this algorithm. The results show that this phase retrieval algorithm is effective, and the method has a potential stability and accuracy for X ray phase contrast imaging technology.     

Analytical solution for nonlinear wave propagation in shallow media using the variational iteration method

An analytical solution is presented for nonlinear surface wave propagation. A variational iteration method (VIM) was employed and time-dependent profiles of the surface elevation level and velocity obtained analytically for different initial conditions. It is shown that the VIM used here is a flexible and accurate approach and that it can rapidly converge to the same results obtained by the Adomian decomposition method.     

Multi-beam method to increase the stability of bridge deflection measurement

The main limitation to the accuracy of the detection of a bridge deflection using laser-based optical solutions is atmospheric turbulence because of the laser beam propagation in ground proximity. The multi-beam method is presented to increase the stability of bridge deflection measurement. It is based on the use of a four-beam optical system and a subpixel resolution algorithm for the measurement of the deflection of a laser beam that propagates through the system. To obtain accurate results, different algorithms for measuring the position of the deflected beam in different optical systems are tested and compared. Based on this comparison, the four-beam method based on the macropixels iteration centroid and four-beam optical system is selected, and an accuracy of 0.16 pixels is obtained by the determination of the beam position in our setup. The proposed method is adopted to detect the bridge deflection and an accuracy of 0.01 mm is gained when the scintillation index C_n² is 3×10⁻⁴ m^−2/3.     

基于近似模型的透平机械多目标优化

本文研究了采用多目标优化算法NSGA—II与前馈神经网络的耦合优化方法的收敛性问题。提出了一种改进的拥挤距离和coarse-to-fine的迭代策略,有效地解决了原耦合方法不收敛的问题。在此基础上,提出了基于数值模拟和耦合优化方法的多目标气动设计框架,并对转子37进行了三目标优化。优化结果表明,多目标优化结果的气动性能...     

水声正交频分复用异步多用户接入方法

水声信道传播时延大导致水声通信网中上行用户同步困难,因此需要接收端具备异步多用户检测能力。针对水声正交频分复用(OFDM)异步多用户系统提出一种基于滑动迭代的异步多用户接入方法,包括重叠截断、干扰消除、频域过采样、多用户信道估计与解码、以及干扰重构5个处理模块,并通过滑动迭代结构同时消除前向干扰和后向干扰影响。仿真结果表明,所提方法在不同异步到达时延,以及有、无信道估计误差下的误比特率均明显低于逐块解调异步接收算法。水池实验结果表明,滑动迭代异步接收算法在多径扩展18 ms、异步到达时延大于或等于循环前缀1.5倍时实现了两用户异步接入的可靠通信。      

耦合微扰简正波方法与孤立子内波

水平变化环境下声场简正波解的计算精度和效率取决于本地简正波的计算方法。提出一种完备的一阶微扰理论方法,并引入迭代算法,获得了本地简正波水平波数和本征函数的精确表达式。数值结果表明,改进后的微扰简正波方法得到的简正波水平波数和本征函数精度比前人方法更高,与KRAKENC计算结果吻合较好,而计算速度比KRAKENC快100倍。同时将微扰简正波方法与耦合简正波理论结合,应用到海水声速水平变化剧烈的孤立子内波群环境。数值结果表明,该方法计算得到的传播损失与COUPLE07在单次散射近似下的计算结果吻合较好,计算速度比COUPLE07快25倍,并将该方法在声场计算中的适用频率提高到了3 kHz。      

Spectral method for exploring patterns of diblock copolymers

A numerical method in Fourier-space is developed to solve the polymeric self-consistent field equations. The method does not require a priori symmetric information. More significantly, periodic structure can be adjusted automatically during the iteration process. In this article, we apply our method to AB linear diblock copolymer melt, thus reproduce all known stable phases, and reveal some meta-stable phases. It is worthy to point out that we also give an efficient strategy to estimating initial values for diblock copolymer system. Finally, by comparing with Matsen–Schick’s method, we show some advantages of our method.     

Radiation transfer in a three-dimensionally inhomogeneous stochastic medium

Within the framework of a small-angle iteration method the propagation of radiation in a three-dimensionally inhomogenous stochastic scattering medium is considered. The results of the analysis show a significant dependence of the average light field on the statistical structure of scattering parameters of the medium. Deceased. Institute of Applied Optics, National Academy of Sciences of Belarus, 11, Belynitskii-Birulya Str., Mogilev, 212793. Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 66, No. 2, pp. 234–240, March–April, 1999.     

Novel infrared image enhancement technology based on the frequency compensation approach

A novel infrared image enhancement method has been proposed in this paper. Our aim is to develop a detail enhancement method which is focused on the frequency feature of the image. The proposed method is following the most popular strategy of enhancing the infrared images nowadays, but concentrating on the frequency domain. Firstly, the original image is separated by a guided image filter into detail layer and the base layer. Quite unlike the traditional methods, we use the guided image filter to eliminate most of the noise and weak signal of the scenario. Then, by a designed iteration process, the higher frequency of the scenario will be calculated back and add to the detail layer. The noise will not be enhanced because the iteration is only focused on the leftover scenario frequency. We run many tests on the raw data captured by the 320 × 256 HgCdTe cooled thermal imager, and make a comparison between our approach with the previous method of bilateral filtering digital detail enhancement and guided image filtering digital detail enhancement. Figures and analytical data show that our method is better than the previous proposed researches. Our method could effectively process the infrared image with less noise and artifacts, which has potential applications in testing, manufacturing, chemical imaging, night vision, and surveillance security.     

不同非局域程度条件下空间光孤子的传输特性

光束在非局域非线性介质中传输由非局域非线性薛定谔方程描述.讨论了在不同非局域程度 条件下，空间光孤子的传输特性.提出了一个基于分步傅里叶算法数值求解孤子波形和分布 的迭代算法.假定介质的非线性响应函数为高斯型，得出了在不同非局域程度条件下空间光 孤子的数值解，并数值证明了它们的稳定性.结果表明，不论非局域程度如何，光束都能以 光孤子态在介质中稳定传输.光孤子的波形是从强非局域时的高斯型过渡到局域时的双曲正 割型，形成孤子的临界功率随非局域程度的减弱而减小，光孤子相位随距离线性增大，相位 的变化率随非局域程度的减弱而减小. 关键词： 非局域非线性薛定谔方程 空间光孤子 临界功率 相位     

Radiation Transfer in a Binary, Stochastic Medium with Strongly Fluctuating Parameters

The propagation of radiation in a three-dimensional, inhomogeneous, stochastic, scattering medium and in an equivalent homogeneous medium has been investigated by the small-angle iteration method. The statistical structure of the scattering parameters of a three-dimensional medium has been described within the framework of the two-dimensional Poisson process and the binary Markov process. The analytical data obtained point to the fact that the three-dimensional inhomogeneity of a stochastic medium significantly influences the radiation transfer in it.