基于干扰对消的红外焦平面非均匀性校正算法

红外焦平面器件的非均匀性产生机理复杂,难以准确拟合探测元响应曲线。提出了一种基于相关干扰抵消的非均匀性校正算法,以预先采集到的一帧黑体面源图像做为自适应干扰对消器的参考输入图像,自适应滤波器由参考输入图像迭代计算出待校正红外图像的空间噪声的最佳估计,实现从空间噪声中提取真实图像信号。自适应滤波算法采用变步长最小均方误差算法,减少了算法的运算量,提高了算法的收敛速度。理论分析以及针对实际红外图像的仿真结果表明,提出的算法校正效果好,收敛速度快,更易于工程实现。     

改进GSO与二维OTSU融合的红外图像多阈值分割方法

针对二维最大累间方差（Nobuyuki OTSU,OTSU）图像分割算法在电气设备故障诊断与定位中,其红外图像的多阈值分割中存在的耗时多、分割精度低、误分割等不足,造成故障区域欠分割或者过分割的问题,提出一种改进的萤火虫算法（glowworm swarm optimization,GSO）与二维OTSU的融合算法来提高电气设备红外图像多阈值分割的实时性与准确度。寻优过程中,将局部寻优扩展到全局寻优,并引入非线性递减步长及新的移动策略对GSO进行优化改进。实验结果表明:该融合算法在分割结果上较二维OTSU及未改进GSO与二维OTSU融合算法更能准确分割运行电气设备图像异常区域,分割速度分别提高19倍、1.28倍,为红外图像早期故障的有效识别与定位奠定基础。     

胞映射搜索最优控制路径的新策略

在胞映射求解最优控制问题的现有方法的基础上，提出了一种利用胞映射搜索最优控制路径的新策略。该策略基于负步长逆向数值积分技巧，采用了一种新的搜索过程，使得求解问题的效率大为提高。同现有的方法相比，新策略节省了几倍，甚至十几倍的时间。而且计算结果可靠，有很强的工程应用价值，为求解高维的最优控制问题可能会提供一定的基础。     

SR-CT技术乘型迭代算法的研究

算法为SR-CT(synchrotron radiation computed tomography)技术中的重要组成部分,其中乘型迭代算法为该技术的一种有效算法。本文对乘型迭代算法进行了研究,并对该算法中的一些重要参数,如迭代步长、迭代次数及初始解进行了优化分析,给出了这些参数对重建图像及运算时间的影响关系以及这些参数间的相互影响关系,同时给出了获取较优初始解的方法以及迭代步长和迭代次数的最佳取值范围。     

求解伪单调变分不等式的两种投影算法

本文提出了两种求解伪单调变分不等式的定步长的投影算法．这与Solodov ＆ Tseng(1996)和He(1997)的变步长策略不同．我们证明了算法的全局收敛性，并且还在一定条件下证明了算法的Q-线性收敛性．     

环形等离子体中电子温度梯度不稳定性的粒子模拟

用粒子模拟方法求解描述环形等离子体中电子温度梯度静电模的回旋动力学方程.方程采用圆磁通面的轴对称环形几何系统，考虑了有限拉摩半径、电子渡越频率k_∥v _∥以及环形漂移(曲率和磁场梯度)运动w_D(v²_⊥,v ²_∥,θ)的效应.简述了粒子模拟的基本方法.采用了四阶变步长积分格式，使计算省时、简便.讨论了模的基本特征 , 并且给出了临界梯度对电子温度与离子温度之 关键词： 电子温度梯度不稳定性 粒子模拟 变步长积分格式 临界梯度     

控制棒下落与流体流动的耦合状态方程及其保辛算法

针对核反应堆内控制棒下落问题,提出了描述控制棒下落与流体流动的耦合非线性状态方程。该状态方程对于落棒过程内不同的流体状态,具有统一的表达形式,可以很方便地处理不同工况下的落棒问题。为高效分析落棒过程,准确捕捉落棒过程内流动状态的突变,并保证时程积分的数值稳定,提出了一种基于时间步长自适应的保辛算法。数值算例表明,提出的数值模型可以采用较大的时间步长精确计算控制棒在下落过程中的位移、速度、加速度、落棒时间等关键数据,计算结果与商业软件所得结果高度吻合。     

结合广义Armijo步长搜索的一类新的三项共轭梯度算法及其收敛特征

In this paper, we consider the convergence properties of a new class of three terms conjugate gradient methods with generalized Armijo step size rule for minimizing a continuously differentiable function f on R^π without assuming that the sequence {xk} of iterates is bounded. We prove that the limit infimum of ‖↓△f(xk)‖ is Zero. Moreover, we prove that, when f(x) is pseudo-convex (quasi-convex) function, this new method has strong convergence results: either xk→x* and x* is a minimizer (stationary point); or ‖xk‖→∞, arg min{f(x) :x∈R^n} =φ, and.f(xk) ↓ inf(f(x) : x∈R^n}. Combining FR, PR, HS methods with our new method, FR, PR, HS methods are modified to have global convergence property.Numerical result show that the new algorithms are efficient by comparing with FR,PR, HS conjugate gradient methods with Armijo step size rule.     

等离子体中Fokker-Planck方程的有限元模拟

利用有限元方法设计了一套相对简单明了的求解Fokker-Planck方程的方案.这个方案不必严格限制计算格点的步长和时间步长,就可以确保分布函数的非负性和粒子数的守恒.通过一维程序模拟,进一步证实了这个方案的可靠性.对于多维问题的分析和一维问题完全一样,所以非常容易将其推广到多维问题.     

计算周期解时出现刚性问题的一种处理方法

给出计算周期解时出现刚性问题的一种具体处理方法,采用MATLAB软件中求解刚性初值问题的解题器ODE15s,同时在处理问题的过程中,为便于计算和减少工作量使用变步长的Gear方法,给出相应的计算公式,并给出数值实例.