一维梯度折射率介质光学特性的反演

在单个方波脉冲入射情况下,利用共轭梯度法对一维梯度折射率介质的折射率、吸收系数以及散射系数进行了反演.正问题采用间断有限元法求解,反问题的解则在正问题的基础上通过共轭梯度法得到.研究结果表明,利用单个方波脉冲入射情况下的时域半球反射率以及时域半球透射率作为测试值能够有效地反演一维梯度折射率介质的光学特性及其分布情况.     

吸收散射性三维矩形介质内辐射源项的反问题

提出了一种由边界出射辐射强度反演吸收散射性三维矩形介质内辐射源项分布的方法。该方法是在辐射传递方程离散坐标近似的基础上,用求目标函数极小值的共轭梯度法进行反演计算。通过对介质辐射特性、光学厚度等参数对反演精度影响的分析,结果表明,即使存在测量误差,本文所提出的方法可较精确地反演辐射源项。     

三维声波方程逆问题的共轭梯度法求解

考虑一个完整的三维声波方程的逆问题．通过构造一个表面声压偏差平方和形式的目标泛函，把声波方程的逆问题转化为一个控制声学特性参数分布使得目标泛函达到最小伍的优化问题．采用共轭梯度法来求解这个优化问题．通过引入一个对偶函数u（x，t），文中用微扰法求得了目标泛函梯度值的解析表达式，从而克服了以往用共轭梯度法求解偏微分方程控制的优化问题时计算目标活函梯度的困难，大大压缩了共轭梯度法计算目标泛函梯度的时间，而且提高了梯度值的计算精度．还进一步进行了反演声学特性参数三维分布的数值仿真计算．共轭梯度法完整解决了三维声波方程的逆问题．     

各向异性介质辐射特性参数联合反演

本文基于BP神经网络方法结合蒙特卡洛和BEER定律辐射传输模拟方法建立了联合反演各向异性散射介质的辐射特性参数模型。首先采用半球透射率结合半球反射率反演模型反演了各向同性介质的吸收系数和散射系数,在此基础上增加准直透射率,建立了联合反演各向异性介质的吸收系数、散射系数和散射不对称因子三参数联合反演模型。反演结果表明该模型能准确反演出介质辐射特性参数,具有实用意义。此外,为了检验测量误差对模型的反演准确性的影响,分别在不同程度测量误差情况下进行反演,结果显示测量误差对散射不对称因子反演值影响较大。     

管道内壁侵蚀形状识别的无网格法研究

本文使用直接配点无网格法结合共轭梯度法对管道内壁面侵蚀状况进行了反演识别。直接配点无网格法使用节点离散求解区域,采用移动最小二乘近似构造试函数,直接配点法构造线性方程组进行导热正问题求解;反演过程采用共轭梯度法使目标函数最小化,Akima三次样条插值将连续的几何边界反演问题转化为离散点几何位置的反演,并最终将这些离散点拟合成为光滑曲线。文中选择两个典型算例对数值方法进行验证,模拟结果表明使用直接配点无网格法结合共轭梯度法进行管道内壁几何边界识别具有较高精度。     

共轭梯度法求解双曲传热多宗量反问题

提出双曲传热反问题热物性参数和边界条件多宗量联合反演的一般数值求解模式,考虑了非均质和分布参数的影响,时域上采用时域精细算法进行离散,建立了便于敏度分析的有限元正演模型.由最小二乘原理建立反演模型,应用共轭梯度法进行求解.探讨了时间步长和测量误差对反演结果的影响,并进行了数值验证.     

基于频域辐射传输模型的光学成像研究

本文以超短脉冲激光照射参与性介质的光学成像为研究背景,分别构建了短脉冲激光在参与性介质内的频域辐射传输正问题模型和根据边界探测所得频域信号重建介质内部光学参数的逆问题模型。在瞬态辐射传输方程的基础上,利用傅里叶变换得到频域辐射传输方程,采用有限体积法求解频域传输方程,模拟超短脉冲激光在二维参与性介质内传输的过程,得到介质边界的出射频域辐射信号。选取共轭梯度法作为反演算法,采用伴随差分模型求解目标函数梯度,重建了二维非均匀参与性介质内不同位置内含物的光学参数分布。结果表明,基于频域辐射传输方程的伴随差分模型能够较为准确地反演多维参与性介质内的光学参数。     

采用共轭梯度法的管内壁温度导热反问题求解

建立了采用共轭梯度法的非稳态二维导热反问题数学模型,并利用FORTRAN语言编写了程序进行计算求解,由同心圆管外壁温度随时间的变化得到管内壁温度随时间的变化.计算结果表明,计算值与准实验值吻合,共轭梯度法能够准确地计算得到管内壁不同位置温度随时间的变化.     

共轭梯度法求解多宗量瞬态热传导反问题

提出瞬态热传导反问题中热物性参数和边界条件多宗量反演的一般数值求解模式,导出了敏度计算公式,并用共轭梯度法求解,探讨了测点数目、测量误差和变量初值对反演结果的影响,并进行了数值验证.     

基于荧光技术的时域近红外光学成像

参与性介质内光学参数信息重建是目前辐射反问题的研究热点,广泛应用于生物医学成像等领域。荧光探针技术可以显著地提高生物体内健康组织与病变组织的对比度,在近红外光学成像领域受到了广泛的关注。本文利用扩散近似方程模拟激发激光和荧光在散射占优介质内的传输,求解介质边界的透反射信号。利用广义Gauss-Markov随机场模型构建正则化项,加入目标函数中以克服反问题的病态特性。采用共轭梯度法求解反问题,反演二维非均匀介质内的荧光产率和荧光寿命。结果表明,本文提出的算法能够很好地反演获得荧光产率和荧光寿命分布,所得图像能够清晰地分辨出介质内部结构。     

Solution of the inverse radiative load problem in a two-dimensional system

An inverse radiation analysis is presented for estimating the temperature and the heat load distributions of the heating surface from the temperature and the heat flux measurements of the heated object. The Monte Carlo method is employed to solve the direct radiation problem. The inverse radiation problem is solved using the conjugate gradient and singular value decomposition methods. The measured data are simulated by adding random errors to the exact solution of the direct problem. The effects of the measurement errors on the accuracy of the inverse analysis are investigated. The study shows that the heat load distribution of the heating surface can be estimated accurately for the exact and noisy data. And the conjugate gradient method is better than the singular value decomposition method since the former can obtain more accurate results if the measurement errors are the same.     

Simultaneous estimation of the temperature and heat rate distributions within the combustion region by a new inverse radiation analysis

A new inverse radiation analysis is presented for estimating the heat rate and temperature distributions in the combustion region from the information of the temperature and heat flux profiles of wall elements in the system. The Monte Carlo method is employed to solve the radiative heat transfer equation. The inverse radiation problem is posed as a minimization problem of the least squares criterion, which is solved by the conjugate gradient method. The performance of the present technique of inverse analysis is evaluated and the effects of the errors of the absorption coefficient, emissivity and convective heat transfer coefficient on the inverse analysis are investigated. The results show that the present technique is robust and yields accurate estimation even with noisy measurement.     

基于红外测温的试件内部缺陷的识别算法研究

内部缺陷形状位置的识别是传热逆问题的重要研究课题.本文对带有内部缺陷的试件建立了二维传热模型,通过有限元法求解了在一定加热条件下试件表面的温度分布;并通过共轭梯度法,提出了根据表面所测温度分布识别试件内部缺陷形状、位置的计算方法.数值实验证明了算法的有效性.     

An Effective Multigrid Preconditioned CG Algorithm for Millimeter Wave Scattering by an Infinite Plane Metallic Grating

An effective multigrid based preconditioned conjugate gradient method is developed to solve electromagnetic large matrix problem for millimeter wave scattering application. By using multigrid technique we restrict the large matrix equation to a relative smaller matrix and which can be solved rapidly. The solution is prolonged as the initial guess for the conjugate gradient (CG) method. Numerical results show that our developed method can reach five times improvement of computational complexity.     

Analysis of Millimeter Wave Scattering by an Electrically Large Metallic Grating Using Wavelet-Based Algebraic Multigrid Preconditioned CG Method

An effective wavelet based multigrid preconditioned conjugate gradient method is developed to solve electromagnetic large matrix problem for millimeter wave scattering application. By using wavelet transformation we restrict the large matrix equation to a relative smaller matrix and which can be solved rapidly. The solution is prolonged as the new improvement for the conjugate gradient (CG) method. Numerical results show that our developed wavelet based multigrid preconditioned CG method can reach large improvement of computational complexity. Due to the automaticity of wavelet transformation, this method is potential to be a block box solver without physical background.     

含噪声包裹相位图的加权最小二乘相位展开算法研究

二维相位展开广泛应用在精密光学测量、自适应光学、合成孔径雷达、图像处理等领域中。为处理含噪声包裹相位图,以预条件共轭斜量法求解权重最小二乘相位展开方程。引入非加权二维离散余弦变换求解泊松方程得到的最小二乘相位解作为共轭斜量法的初始解,从而加快了收敛速度,同时提出一种新质量图确定算法求解过程中的权重项。计算机模拟和试验表明算法计算速度快,能有效地消除传统路径积分法在处理信噪比低包裹相位图时的"拉线"现象,是一种有效的相位展开方法。     

Inverse estimation of the pulse parameters of a time-varying laser pulse to obtain desired temperature at the material surface

A proper selection of the pulse parameters is essential to achieve desired temperature at the material surface. This leads to obtain the required metallurgical changes in the surface vicinity when a time-varying laser pulse is used in a heating process such as surface modification. In this paper, the conjugate gradient method (CGM) for parameter estimation is successfully applied to estimate the unknown laser pulse parameters for those purposes during laser heating process. The determination of the pulse parameters is treated as a one-dimensional, transient, non-linear inverse heat conduction problem (IHCP). Based on a sensitivity analysis, the inverse problem is solved as an optimization problem comparing a desired temperature at the surface and a calculated one where the objective function is minimized by CGM. The method has been applied to a test case of a heating process on steel, appropriate pulse parameters and desired temperature distribution can also be returned.