窄缝耦合的快速估算法

 基于小孔耦合理论和腔体格林函数,提出了长宽比大于10的窄缝耦合的快速算法,并与Micro-Stripes软件的计算结果进行比较,两者基本一致。研究了窄缝位置、数量及腔体尺寸对屏蔽效能的影响,结果表明：每增加一条相同尺寸、相同取向的窄缝可使腔体的腔体屏蔽效能下降约6 dB;窄缝长度不变的条件下,长宽比每增大1倍则屏蔽效能增加1 dB,并且增大腔体任一方向的尺寸都可以使屏蔽效能增大;除腔体的谐振频率与软件计算结果稍有差异之外,快速算法与软件的计算结果吻合很好,而快速算法的速度远大于软件计算速度,且适合于高频范围。     

基于电磁拓扑的多腔体屏蔽效能快速算法

提出基于电磁拓扑理论计算开孔多腔体屏蔽效能的快速方法.首先给出双腔体等效电路和电磁拓扑信号流图,并推导孔缝节点处的散射矩阵,给出拓扑网络的散射矩阵方程和传输矩阵方程,获得双腔体的广义Baum-Liu-Tesche(BLT)方程.在此基础上研究了开孔三腔体,包括串型级联三腔体和串并型混合级联三腔体的广义BLT方程.对于串型级联三腔体,其电磁拓扑网络和广义BLT方程在双腔体基础上直接扩展即可获得.而对于串并型混合级联三腔体,通过将位于三腔体公共面上的孔缝等效为三端口网络节点,并根据三端口网络散射参数定义推导获得该节点的散射矩阵,最终得到串并型混合级联三腔体的广义BLT方程.本文方法对双腔体的计算结果与文献结果和实验结果相符合,对3组不同类型和尺寸开孔腔的屏蔽效能的计算结果与时域有限差分法计算结果符合较好.该算法不仅效率高,通过对所有计算结果和实验结果的误差统计分析,表明该算法具有较高的计算准确度.     

开孔金属腔屏蔽效能快速算法适用性评估

对开孔金属腔屏蔽效能（SE）的四种快速算法传输线等效电路法、EMT法、ILCM法、Modal-MOM法的适用性进行评估,为不同实际工程问题选用恰当的算法提供了依据。为了解决传统SE结果描述方法不适用于快速算法适用性评估的问题,采用绝对误差和曲线相关系数这两个参数,提出SE随金属腔和孔缝尺寸单独变化的Contour图描述方法。在此基础上,利用864组不同形状、尺寸矩形金属腔的SE计算结果,分析并确定了四种SE快速算法满足不同计算精度时适用的腔体模式数。结果表明,在同等计算精度条件下,Modal-MOM算法可适用的腔体模式数最多,ILCM算法次之,而传输线等效电路法适用腔体模式数最少,这些结果为实际工程应用提供了指导。     

基于点和直线段对应的扩展正交迭代位姿估计算法

正交迭代算法是基于点特征,全局且快速收敛的位姿估计算法,是目前性能最优的实时位姿估计算法之一.对摄像机观测到的直线段提出了其不确定性描述的新方法:目标空间直线段误差,并把该误差融入到正交迭代算法的位姿计算过程中,形成扩展正交迭代算法.该算法可以同时利用点和直线段特征.算法先把点共面性方程表达为与正交迭代算法中点共线性方程一致的数学形式,再根据目标空间点共面性误差定义了目标空间直线段共面性误差.接着把目标空间点共线性误差和直线段共面性误差二者之和作为误差函数.最后推导出使该函数最小化的迭代求解过程.实验结果表明,本文算法是有效的,精确的.与正交迭代算法的计算结果对比,由于可以同时利用直线段特征,扩展正交迭代算法的位姿估计误差更低,抗噪声性能更强,算法更稳定.     

模拟短脉冲激光在介质中传输的扩散综合有限元法

建立有限元模型,通过求解瞬态辐射传输方程模拟短脉冲激光在半透明介质中的传输.针对散射占优性半透明介质内辐射传输求解效率较差的问题,采用扩散综合加速迭代算法,提高计算效率,缩短计算时间.结果表明:采用精确解析式描述脉冲激光散射源项的求解策略可以获得准确的计算结果,精确地模拟快速变化的波前,不会产生数值扩散和数值振荡.此外,扩散综合迭代算法的计算时间仅为源项迭代的50%~60%.     

腔体内传输线耦合的电磁仿真软件与传输线方程的混合算法

研究了一种有限积分法软件与传输线方程相结合的混合算法,用于解决复杂电磁环境下屏蔽腔体内传输线的电磁耦合问题。利用有限积分法软件实现屏蔽腔体的建模,仿真得到腔体内部空间电磁场分布,并设置电场探针提取出传输线的激励场。利用传输线方程建立腔体内传输线的耦合模型,将得到的传输线激励场引入到传输线方程作为等效分布电压和电流源。利用时域有限差分(FDTD)格式对传输线方程进行离散,从而迭代求解出传输线终端负载上的电压和电流响应。通过与文献以及传统数值算法的计算结果进行对比,验证混合算法的正确性。研究表明,该混合算法在模拟电大尺寸腔体内传输线的电磁耦合方面,具有较高的精度和计算效率。     

任意形状金属腔体内电磁脉冲的三维计算

提出了模拟任意形状腔体中的内电磁脉冲的三维直角坐标系时域有限差分(FDTD)算法。该算法采用FDTD共形网格技术模拟任意形状腔体的边界,可以解决腔体内非对称的边界问题。推导了射线斜入射的差分方程,进行了三维数值计算,并采用直角坐标系FDTD算法和柱坐标系FDTD算法计算了射线斜入射圆柱腔体产生的内电磁脉冲,二者吻合很好,验证了直角坐标系FDTD算法正确性。     

水平层状横向同性地层中频率测深资料全参数快速迭代反演

提出水平层状横向同性地层中频率测深资料的全参数快速迭代反演算法,以便从测量资料中同时确定各个地层的横向、 纵向电导率及层界面深度.首先,利用水平层状介质中电磁场并矢Green函数在频率波数域中的解析解和Sommerfeld积分的 快速计算技术确定正演响应.然后,利用摄动理论和Fourier逆变换公式,研究建立一套快速求解全参数Fréchet导数的 有效算法,并利用规范化处理和奇异值分解技术给出迭代反演过程.最后,给出数值试验证明反演理论的有效性以及反演算法 的抗噪声能力.     

一种内置条状金属板的双层金属腔体屏蔽效能的理论模型

针对复杂屏蔽腔体往往是由多个空间构成的实际情况, 本文构建了内置条状金属板的双层金属腔体物理解析模型, 将外层腔体的近场电磁干扰等效为电偶极子, 基于Bethe小孔耦合理论并利用推广的腔体格林函数推导了内腔体的电磁场分布的近似表达式. 利用该解析模型计算分析了条状金属板的位置和方向对屏蔽效能的影响. 通过计算结果与全波仿真软件CST仿真结果的对比, 证实了本文所建理论模型的有效性, 为复杂腔体屏蔽效能的快速计算提供了理论参考.     

基于强度传输方程的相位检索

为了从直接测量的强度图像中计算出相位信息,研究了基于强度传输方程的相位检索技术.设计了一个实用的相位检索系统,包含经典傅里叶相位检索和基于整体变分的相位检索算法.通过测量聚焦平面和前后散焦平面光波场的空间强度,求解二阶微分方程得到相位.与迭代相位检索技术相比克服了其迭代不确定性、收敛速度慢等缺点.实验表明,利用强度传输方程检索相位的方法可以有效快速地从强度图像中计算出相位|整体变分算法相比傅里叶算法,在恢复相位的同时可以保持良好的边缘.     

An iterative method for improving the quality of reconstruction of a three-dimensional surface

A complex image with constraints imposed on the amplitude and phase image components is processed using the Gerchberg iterative algorithm for the first time. The use of the Gerchberg iterative algorithm makes it possible to improve the quality of a three-dimensional surface profile reconstructed by the previously proposed method that is based on the multiangle projection of fringes and the joint processing of the obtained images by Fourier synthesis.     

Phase-only stereoscopic hologram calculation based on Gerchberg–Saxton iterative algorithm

A phase-only computer-generated holography(CGH) calculation method for stereoscopic holography is proposed in this paper.The two-dimensional(2D) perspective projection views of the three-dimensional(3D) object are generated by the computer graphics rendering techniques.Based on these views,a phase-only hologram is calculated by using the Gerchberg–Saxton(GS) iterative algorithm.Comparing with the non-iterative algorithm in the conventional stereoscopic holography,the proposed method improves the holographic image quality,especially for the phase-only hologram encoded from the complex distribution.Both simulation and optical experiment results demonstrate that our proposed method can give higher quality reconstruction comparing with the traditional method.     

Three-dimensional electromagnetic nonlinear inversion in layered media by a hybrid diagonal tensor approximation: Stabilized biconjugate gradient fast Fourier transform method

This paper presents an efficient three-dimensional nonlinear electromagnetic inversion method in a multilayered medium for radar applications where the object size is comparable to the wavelength. In the first step of this two-step inversion algorithm, the diagonal tensor approximation is used in the Born iterative method. The solution of this approximate inversion is used as an initial guess for the second step in which further inversion is carried out using a distorted Born iterative method. Since the aim of the second step is to improve the accuracy of the inversion, a full-wave solver, the stabilized biconjugate-gradient fast Fourier transform algorithm, is used for forward modelling. The conjugate-gradient method is applied at each inversion iteration to minimize the functional cost. The usage of an iterative solver based on the FFT algorithm and the developed recursive matrix method combined with an interpolation technique to evaluate the layered medium Green's functions rapidly, makes this method highly efficient. An inversion problem with 32 768 complex unknowns can be solved with 1% relative error by using a simple personal computer. Several numerical experiments for arbitrarily located source and receiver arrays are presented to show the high efficiency and accuracy of the proposed method.     

Three-dimensional electromagnetic nonlinear inversion in layered media by a hybrid diagonal tensor approximation: Stabilized biconjugate gradient fast Fourier transform method

This paper presents an efficient three-dimensional nonlinear electromagnetic inversion method in a multilayered medium for radar applications where the object size is comparable to the wavelength. In the first step of this two-step inversion algorithm, the diagonal tensor approximation is used in the Born iterative method. The solution of this approximate inversion is used as an initial guess for the second step in which further inversion is carried out using a distorted Born iterative method. Since the aim of the second step is to improve the accuracy of the inversion, a full-wave solver, the stabilized biconjugate-gradient fast Fourier transform algorithm, is used for forward modelling. The conjugate-gradient method is applied at each inversion iteration to minimize the functional cost. The usage of an iterative solver based on the FFT algorithm and the developed recursive matrix method combined with an interpolation technique to evaluate the layered medium Green's functions rapidly, makes this method highly efficient. An inversion problem with 32 768 complex unknowns can be solved with 1% relative error by using a simple personal computer. Several numerical experiments for arbitrarily located source and receiver arrays are presented to show the high efficiency and accuracy of the proposed method.     

Accurate complex modulation by the iterative spatial cross-modulation method

A method is proposed to realize accurate spatial complex modulation based on the spatial cross-modulation method(SCMM) via a phase-only spatial light modulator. The conventional SCMM cannot achieve high quality reconstruction, especially when the diffusion ratio is small. We propose an iterative algorithm in the calculation of a computer-generated hologram to implement accurate complex modulation. It enables us to generate a high quality reconstruction under a small diffusion ratio. The feasibility of the method is verified by both a numerical simulation and an optical experiment.     

有限口径非周期亚波长衍射聚焦偏转器的严格矢量优化

有限口径非周期亚波长衍射光学器件的亚波长结构特性使得传统标量衍射理论失效,必须借助严格电磁理论进行器件的设计。提出一种循环平面波谱算法(IterativePlaneWaveSpectrumAlgorithm,IPWSA)作为器件优化算法,利用二维时域有限差分方法作为严格电磁计算模式,根据Farn方法对器件面型进行亚波长处理,设计出有限口径非周期具有亚波长结构的衍射聚焦偏转器件,偏转角为24.8°,该器件可与激光器集成构成微型集成光学系统,也可以用在红外量子阱系统中。详细分析了算法原理,设计实例验证了算法的优化效果。     

Conditional Source-term Estimation using dynamic ensemble selection and parallel iterative solution

A modified version of the Least-Square QR-factorisation (LSQR) algorithm has been implemented in conjunction with Conditional Source-term Estimation (CSE) for lean, turbulent premixed methane–air combustion via Large Eddy Simulation (LES). The iterative solver can reduce computational times by an order of magnitude during the inversion phase of CSE in comparison with the conventional LU-decomposition method. The advantages of iterative and parallel iterative solvers become more prominent as the size of the system increases. The ensemble selection procedure for computing averages within localised regions of the simulation domain has also been updated to a dynamic routine. This allows for more flexible and efficient allocation of computational resources along with reduced input from the user, especially for complex geometries. Preliminary LES calculations have shown that the implementation of an iterative solver and a dynamic ensemble selection algorithm will reduce computational times significantly with negligible error contribution for one-condition CSE, which is applicable to purely premixed or non-premixed turbulent combustion problems. In addition, these algorithms provide the foundation for exceptional computational cost savings for the inversion in two-condition CSE, or Doubly Conditional Source-term Estimation (DCSE), which has shown promise for predicting partially-premixed combustion. Parallel computation of the inverse solution is particularly beneficial to DCSE as the computational cost of the inversion process is considerably larger than in one-condition CSE.     

Vector fuzzy control iterative algorithm for the design of sub-wavelength diffractive optical elements for beam shaping

The vector fuzzy control iterative algorithm (VFCIA) is proposed for the design of phase-only sub-wavelength diffractive optical elements (SWDOEs) for beam shaping. The vector diffraction model put forward by Mansuripur is applied to relate the field distributions between the SWDOE plane and the output plane. Fuzzy control theory is used to decide the constraint method for each iterative process of the algorithm. We have designed a SWDOE that transforms a circular flat-top beam to a square irradiance pattern. Computer design results show that the SWDOE designed by the VFCIA can produce better results than the vector iterative algorithm (VIA). And the finite difference time-domain method (FDTD), a rigorous electromagnetic analysis technique, is used to analyze the designed SWDOE for further confirming the validity of the proposed method.     

一种基于角谱理论的改进型相位恢复迭代算法

进行星间激光通信的光学发射天线光束整形器设计时, 首要解决的问题是根据输入光场及理想的输出光场, 确定整形器的相位分布, 其核心就是相位恢复. 基于角谱传播理论, 在传统 Gerchberg-Saxton (G-S)迭代算法的基础上, 提出了一种幅度梯度加成迭代算法, 给出了算法的详细流程与分析. 与G-S相比, 新算法利用迭代过程, 构建光场幅度反馈回路, 利用梯度搜索最佳迭代路径, 两者的联合作用加速其迭代收敛进程. 数值仿真表明, 新算法的单位迭代次数所引起迭代误差下降的速度是G-S算法的1.7倍, 其收敛速度明显优于G-S算法; 对不同的随机初始相位, 新算法都能进行有效迭代, 表现出适应性强, 且收敛一致性好的优点. 幅度梯度加成迭代算法为复杂光场的高效相位恢复提供了一种新思路, 为设计各种衍射光学元件提供了技术支持. 关键词： 相位恢复 迭代算法 角谱理论 光通信     

Iterative retrieval of one-dimensional x ray wave field using a single intensity measurement

The problem of retrieving a complex function from the modulus of its Fourier transform has non-unique solutions in one dimension. Therefore iterative phase retrieval methods cannot in general be confidently applied to one-dimensional problems, due to the presence of ambiguities. We present a method for a posteriori reduction of the ambiguities based on the correlation analysis of the solution of a large number of runs of an iterative phase retrieval algorithm with different random starting phases. The method is applied to experimentally measured diffraction patterns from an x ray waveguide illuminated by hard x rays. We demonstrate the possibility of retrieving the complex wave field at the exit face of the waveguide and compare the result with theoretical prediction.