FGMs稳态热传导分析的重心Lagrange插值配点法

提出一种求解二维功能梯度材料（FGMs）稳态热传导问题的重心Lagrange插值配点法.基于Chebyshev节点构造二维重心Lagrange插值函数及其偏导数,然后基于配点法将其直接代入FGMs热传导问题的控制方程和边界条件,得到系统离散方程.重心Lagrange插值配点法是一种真正的无网格方法,很好地融合了重心Lagrange插值和配点格式的优势,具有高效、稳定、高精度和易于数值实现的优点.采用重心Lagrange插值配点法分别对指数型、二次型和三角型FGMs热传导问题进行数值模拟.结果表明：该方法具有较高的计算效率和计算精度,对材料梯度参数的变化不敏感.可以进一步拓展到FGMs瞬态问题和FGMs的热力耦合分析.     

Estimation of mutual information objective function based on Fourier shift theorem: an application to eddy current distortion correction in diffusion tensor imaging

Diffusion tensor imaging requires correction of eddy current distortion in diffusion-weighted images. An effective retrospective correction approach is to transform a diffusion-weighted image to maximize the mutual information (MI) between the transformed diffusion-weighted image and the corresponding T2-weighted image. In the literature, either linear interpolation or partial volume interpolation is applied to estimate the MI objective function. However, these interpolation methods induce artifacts to the MI objective function, thus compromising correction results. In this work, the MI objective function is estimated based on interpolation using Fourier shift theorem. This method eliminates the artifacts incurred with the aforementioned interpolation methods. The algorithm is further improved by approximating pixel values using their nearest neighbors in the up-sampled spatial domain, resulting in dramatically increased computational efficiency without compromising the correction results. The effects of varying the number of quantization levels and using Parzen window filtering to smooth the MI objective function are also investigated to obtain optimized algorithm parameters. The diffusion tensor image quality after applying the proposed distortion correction method is significantly improved visually.     

基于平面问题的位移压力混合配点法

引入压力变量,将弹性力学控制方程表达为位移和压力的耦合偏微分方程组,采用重心插值近似未知量,利用重心插值微分矩阵得到平面问题控制方程的矩阵形式离散表达式.采用重心插值离散位移和应力边界条件,采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法求解过约束方程组,得到平面问题位移数值解.数值算例验证了所提方法的有效性和计算精度.     

浅海非均匀水平阵宽带声场信号无源孔径扩展

大孔径水平阵对于浅海低频水声物理实验研究和应用至关重要,然而受实际情况制约,通常使用的水平基阵孔径十分有限,通过孔径扩展处理来扩展基阵孔径是一条重要途径。传统的无源孔径扩展方法是建立在均匀线阵、匀速相对运动和相干窄带连续信号的基础上的,难以适用于非均匀阵以及非连续宽带声源的情况。针对这些问题,提出了非均匀阵宽带声场信号的无源孔径扩展方法。使用静止布设在海底的非均匀水平短阵接收运动声源重复发射的宽带声信号,先开展均匀空间插值,然后在阵元域和波束域进行宽带扩展孔径处理。仿真和实验结果表明,在水深约70 m的浅海波导环境中,纵向间隔27.5 m的2个阵元接收20~200 Hz宽带声场,其空间插值结果与真值的相关系数大于0.99,说明宽带声场插值方法的合理性。在水平非均匀、纵向孔径250 m的短阵接收声场无法分析简正波频散特征的情况下,仿真和实验数据经过宽带无源孔径扩展处理得到孔径大于1 km的均匀线阵的声场,能高分辨区分各阶简正波,证明所提方法是有效的。      

求解非线性偏微分方程的自适应小波精细积分法

以Burgers方程为例,提出了一种求解偏微分方程的自适应多层插值小波配置法,通过引入一种具有插值特性的拟Shannon小波并利用插值小波理论构造了多层自适应插值小波算子,从而在空间实现了偏微分方程的自适应离散.另外,精细时程积分方法和外推法的引入不但有助于提高求解速度和数值结果的精度,而且使时间积分步长的选取具有了自适应性.     

重心Lagrange插值配点法求解二维双曲电报方程

提出一种求解二维双曲电报方程的高精度重心Lagrange插值配点法.采用重心Lagrange插值构造包含时间和空间变量的近似函数.在给定Chebyshev-Gauss-Lobatto节点上,将多变量重心Lagrange插值近似函数代入双曲电报方程及其定解条件,得到离散代数方程组.包含狄里克雷和诺依曼边界条件的数值算例表明,本文方法程序实现方便并具有高精度,可应用于求解高维问题.     

Ising Model for Interpolation of Spatial Data on Regular Grids

We apply the Ising model with nearest-neighbor correlations (INNC) in the problem of interpolation of spatially correlated data on regular grids. The correlations are captured by short-range interactions between “Ising spins”. The INNC algorithm can be used with label data (classification) as well as discrete and continuous real-valued data (regression). In the regression problem, INNC approximates continuous variables by means of a user-specified number of classes. INNC predicts the class identity at unmeasured points by using the Monte Carlo simulation conditioned on the observed data (partial sample). The algorithm locally respects the sample values and globally aims to minimize the deviation between an energy measure of the partial sample and that of the entire grid. INNC is non-parametric and, thus, is suitable for non-Gaussian data. The method is found to be very competitive with respect to interpolation accuracy and computational efficiency compared to some standard methods. Thus, this method provides a useful tool for filling gaps in gridded data such as satellite images.     

改进S/B算法的近红外光谱模型转移

针对模型转移中S/B算法对于非线性问题的局限性,在传统S/B算法进行线性拟合、偏最小二乘法求参数的基础上加以改进,提出了引入变量的高次幂、使用Lagrange插值法与Newton插值法求待定系数和插值多项式来解决两组数据的非线性问题。为了验证改进算法的有效性,先对主机样品建模并分别预测主机和子机样品,然后通过实验数据和评价指标,筛选出最佳函数关系进行子机模型校正,并分别用改进的S/B算法和传统的S/B算法对子机未知样本进行预测。实验结果表明：直接用主机原模型对子机预测的值与真实值差距较大,利用改进的S/B算法(H-S/B)比传统的S/B算法预测值更接近真实值。改进的S/B算法提高了预测值的准确性,解决了传统S/B算法的非线性问题,实现了更好的模型转移效果,增强了网络化模型应用的通用性。     

脉冲激光测距时间间隔测量技术

在脉冲激光测距系统中,为提高时间间隔测量的精度,采用插值法进行时间间隔测量。在分析传统的数字计数法测量原理与误差的基础上,重点研究了插值法。其测量对象针对传统数字计数法中待测脉冲上升沿与下一个量化时钟脉冲上升沿之间的时间间隔。测量采用电容充放电技术,把时间间隔加入到一个时间扩展模块(通常为一只高精度的电容),实现时间上的放大,再对放大后的时间进行测量,可提高时间测量的精度。利用该方法得到的时间间隔测量精度可达到100 ps,对应于1.5 cm的测距精度。     

非均匀采样干涉数据光谱反演技术研究

干涉光谱成像仪获取的干涉数据是一种中间数据,需要进行光谱反演才能够得到目标光谱数据,傅里叶变换方法是常规的光谱反演方法。由于干涉数据中存在非均匀采样问题,若忽略光谱混叠,直接采用快速傅里叶变换会导致反演光谱的失真,难以满足实时处理需求。针对非均匀采样干涉数据的光谱反演需求,将插值及非均匀快速傅里叶变换(NUFFT)方法应用到光谱反演处理中,对过采样及部分欠采样情况下的非均匀采样干涉数据,提出了相应的光谱反演方法,并分析了方法的适用性。最后对过采样及部分欠采样情况下的光谱反演方法进行计算机仿真,过采样情况下采用NUFFT方法反演光谱的精度要明显高于插值方法,而部分欠采样情况下插值方法反演光谱的精度要明显高于NUFFT方法,并对欠采样造成的光谱混叠有一定的修正,验证了方法的有效性。     

奇异源项问题的重心插值数值解

点源热传导问题和集中力作用梁变形问题的数学模型中,源项为奇异的Delta函数.采用数值稳定性好的重心型插值近似未知函数,利用Delta函数与Heaviside函数的导数关系以及Delta函数的积分筛选性,建立求解含有奇异源项问题的重心插值配点法和重心插值Galerkin法.通过数值算例比较两个方法的有效性和计算精度.     

Spectral collocation method with a flexible angular discretization scheme for radiative transfer in multi-layer graded index medium

The spectral collocation method (SCM) is employed to solve the radiative transfer in multi-layer semitransparent medium with graded index. A new flexible angular discretization scheme is employed to discretize the solid angle domain freely to overcome the limit of the number of discrete radiative direction when adopting traditional S_N discrete ordinate scheme. Three radial basis function interpolation approaches, named as multi-quadric (MQ), inverse multi-quadric (IMQ) and inverse quadratic (IQ) interpolation, are employed to couple the radiative intensity at the interface between two adjacent layers and numerical experiments show that MQ interpolation has the highest accuracy and best stability. Variable radiative transfer problems in double-layer semitransparent media with different thermophysical properties are investigated and the influence of these thermophysical properties on the radiative transfer procedure in double-layer semitransparent media is also analyzed. All the simulated results show that the present SCM with the new angular discretization scheme can predict the radiative transfer in multi-layer semitransparent medium with graded index efficiently and accurately.     

Nonuniform sampled angular spectrum method by using trigonometric interpolation

The angular spectrum method(ASM) is a popular numerical approach for scalar diffraction calculations. However,traditional ASM has an inherent problem in that nonuniform sampling is precluded. In an attempt to address this limitation,an improved trigonometric interpolation ASM(TIASM) is proposed, in which the fast Fourier transform(FFT) is replaced by a trigonometric interpolation. The results show that TIASM is more suitable to situations in which the source field has a simple and strong frequency contrast, irrespective of whether the original phase distribution is a plane wave or a Fresnel zone plate phase distribution.     

Multisoliton Solutions of the (2 1)-Dimensional KdV Equation

Using the extension homogeneous balance method,we have obtained some new special types of soliton solutions of the (2 1)-dimensional KdV equation.Starting from the homogeneous balance method,one can obtain a nonlinear transformation to simple (2 1)-dimensional KdV equation into a linear partial differential equation and two bilinear partial differential equations.Usually,one can obtain only a kind of soliton-like solutions.In this letter,we find further some special types of the multisoliton solutions from the linear and bilinear partial differential equations.``     

M尺度函数及其在图像缩放处理中的应用

针对图像缩小和放大处理中经典方法存在的问题,提出了M尺度插值法,推广了小波理论中的双尺度关系.对构造M尺度函数的方法和M尺度函数的逼近性质进行了研究,以满足3尺度关系的样条函数为例给出了3尺度函数插值算法,该算法与经典的最临近插值法、双线性插值法和双三次插值法相比,不仅保持了较高的清晰度,而且图像边缘较平滑,提高了图像的质量,实践证明了该算法的有效性.本方法还可以推广到图像处理的其他应用中.     

Multisoliton Solutions of the (2+1)-Dimensional KdV Equation

Using the extension homogeneous balance method,we have obtained some new special types of soliton solutions of the (2+1)-dimensional KdV equation.Starting from the homogeneous balance method,one can obtain a nonlinear transformation to simple (2+1)-dimensional KdV equation into a linear partial differential equation and two bilinear partial differential equations.Usually,one can obtain only a kind of soliton-like solutions.In this letter,we find further some special types of the multisoliton solutions from the linear and bilinear partial differential equations.     

Transient heat conduction analysis using the MLPG method and modified precise time step integration method

The meshless local Petrov–Galerkin (MLPG) method in conjunction with the modified precise time step integration method in the time domain is proposed for transient heat conduction analysis in this paper. The MLPG method is often referred to as a truly meshless method because it requires no elements or background cells for either field interpolation or background integration. Local weak forms are developed using weighted residual method locally from the partial differential equation of transient heat conduction. In order to simplify the treatment of essential boundary conditions, the natural neighbour interpolation (NNI) is employed for the construction of trial functions. Moreover, the three-node triangular FEM shape functions are taken as test functions to reduce the order of integrands involved in domain integrals. The semi-discrete heat conduction equation is solved numerically with modified precise time step integration method in the time domain. The availability and accuracy of the present method for transient heat conduction analysis are tested through numerical examples.     

场方法的改进及其在积分Riemann-Cartan空间运动方程中的应用

和Hamilton-Jacobi方法类似,Vujanovi?场方法把求解常微分方程组特解的问题转化为寻找一个一阶拟线性偏微分方程(基本偏微分方程)完全解的问题,但Vujanovi?场方法依赖于求出基本偏微分方程的完全解,而这通常是困难的,这就极大地限制了场方法的应用.本文将求解常微分方程组特解的Vujanovi?场方法改进为寻找动力学系统运动方程第一积分的场方法,并将这种方法应用于一阶线性非完整约束系统Riemann-Cartan位形空间运动方程的积分问题中.改进后的场方法指出,只要找到基本偏微分方程的包含m(m≤ n,n为基本偏微分方程中自变量的数目)个任意常数的解,就可以由此找到系统m个第一积分.特殊情况下,如果能够求出基本偏微分方程的完全解(完全解是m=n时的特例),那么就可以由此找到≤系统全部第一积分,从而完全确定系统的运动.Vujanovi?场方法等价于这种特殊情况.     

Alternating Direction Implicit Orthogonal Spline Collocation on Non-Rectangular Regions

The alternating direction implicit (ADI) method is a highly efficient technique for solving multi-dimensional time dependent initial-boundary value problems on rectangles. When the ADI technique is coupled with orthogonal spline collocation (OSC) for discretization in space, we not only obtain the global solution efficiently, but the discretization error with respect to space variables can be of an arbitrarily high order. In [2], we used a Crank Nicolson ADI OSC method for solving general nonlinear parabolic problems with Robin's boundary conditions on rectangular polygons and demonstrated numerically the accuracy in various norms. A natural question that arises is: Does this method have an extension to non-rectangular regions? In this paper, we present a simple idea of how the ADI OSC technique can be extended to some such regions. Our approach depends on the transfer of Dirichlet boundary conditions in the solution of a two-point boundary value problem (TPBVP). We illustrate our idea for the solution of the heat equation on the unit disc using piecewise Hermite cubics.     

Texton-based super-resolution for achieving high spatiotemporal resolution in hybrid camera system

Many super-resolution methods have been proposed to enhance the spatial resolution of images by using iteration and multiple input images. In a previous paper, we proposed the example-based super-resolution method to enhance an image through pixel-based texton substitution to reduce the computational cost. In this method, however, we only considered the enhancement of a texture image. In this study, we modified this texton substitution method for a hybrid camera to reduce the required bandwidth of a high-resolution video camera. We applied our algorithm to pairs of high- and low-spatiotemporal-resolution videos, which were synthesized to simulate a hybrid camera. The result showed that the fine detail of the low-resolution video can be reproduced compared with bicubic interpolation and the required bandwidth could be reduced to about 1/5 in a video camera. It was also shown that the peak signal-to-noise ratios (PSNRs) of the images improved by about 6 dB in a trained frame and by 1.0–1.5 dB in a test frame, as determined by comparison with the processed image using bicubic interpolation, and the average PSNRs were higher than those obtained by the well-known Freeman’s patch-based super-resolution method. Compared with that of the Freeman’s patch-based super-resolution method, the computational time of our method was reduced to almost 1/10.