基于分形的捷联惯组历次测试数据混合插值算法

针对捷联惯组历次测试数据小样本、非等间隔、非线性的问题,提出了一种基于分形插值的三次混合插值算法.通过第一次分形插值保证原始测试数据的变化趋势;通过第二次样条函数插值保证了插值的准确性,实现原始测试数据等间隔化;通过第三次分段线性插值,扩大样本容量,同时保证了原有测试序列的统计特性不变.实例分析表明,该算法很好的实现了对捷联惯组历次测试数据的等间隔处理和样本容量扩大,为捷联惯组历次测试数据小样本建模分析提供良好的基础.     

Multi-dimensional self-affine fractal interpolation model in tensor form

Iterated Function System (IFS) models have been explored to represent discrete sequences where the attractor of an IFS is self-affine either in R ² or R ³ (R is the set of real numbers). In this paper, the self-affine IFS model is extended from R ³ to R ⁿ (n is an integer and greater than 3), which is called the multi-dimensional self-affine fractal interpolation model. This new model is presented by introducing the defined parameter “mapping partial derivative”. A constrained inverse algorithm is given for the identification of the model parameters. The values of this new model depend continuously on all of the variables. That is, the function is determined by the coefficients of the possibly multi-dimensional affine maps. So the new model is presented as much more general and significant. Moreover, the multi-dimensional self-affine fractal interpolation model in tensor form is more terse than in the usual matrix form.     

On star product fractal surfaces and their dimensions

Ｇｅｎｅｒａｌｌｙｓｐｅａｋｉｎｇ,ｆｒａｃｔｕｒｅｓｕｒｆａｃｅｓｏｆｍａｔｅｒｉａｌｓａｒｅｒｏｕｇｈ,ｉｒｒｅｇｕｌａｒａｎｄｒａｎｄｏｍ．Ｔｈｅｆｒａｃｔｕｒｅａｎｄｗｅａｋｎｅｓｓｉｎｒｏｃｋ,ｃｏｎｃｒｅｔｅ,ｃｅｒａｍｉｃｓａｎｄｍｅｔａｌｓｉｇｎｉｆｉｃａｎｔｌｙａｆｆｅｃｔｔｈｅｄｅｆｏｒｍａｔｉｏｎ,ｓｔｒｅｎｇｔｈ,ａｎｄｃｏｎｄｕｃｔｉｖｉｔｙｏｆｔｈｅｓｅｍａｔｅｒｉａｌｓ．Ｍａｎｙｐａｐｅｒ（ｓｕｃｈａｓ［１～３］）ｆｏｃｕｓｅｄｏｎｔｈｅｆｒａｃｔｕｒｅｓｕｒｆａ…     

三维自然单元法插值形函数的导数计算

根据Voronoi胞的几何性质,获得了积分点的二阶Voronoi胞顶点的表达式,并对各邻近结点相关的顶点进行排序以使其生成的二阶Voronoi胞切割面为凸多边形,从而获得各切割凸多边形面域的面积表达式;最后,基于复合函数链式求导法则,获得了三维自然单元法non-Sibson插值形函数导数的显式格式。相比Lasserre算法,该方法具有直观、便于编程且计算量小的特点。悬臂梁的算例结果进一步说明了该方法的可靠性,证实了文献[2,7,8]关于自然单元法具有比有限元中常应变单元更高的精度,理论上和双线性单元的精度同阶的结论。     

用Voronoi图进行新型自然邻居插值的几何学方法与特性

新的基于Voronoi图的Natural Neighbour插值是自然单元法的数学基础，也是一种新型的几何插值方法，具有与其他传统常用插值不同的构造方法，并表现出一定的优越性。本文介绍了基于Natural Noighbour关系的Sibson插值和non-Sibsonian插值，并与有限元法和无单元法所用的插值方法，就离散插值方案和网格总体特性、形函数支撑域、本征边界条件、空间维数扩展与计算工作量等诸问题进行了比较分析。     

Numerical simulation of ground water flow via a new approach to the local radial point interpolation meshless method

A novel approach to local radial point interpolation meshless (LRPIM) method is introduced to investigate the influence of leakage on tidal response in a coastal leaky confined aquifer system, based on a local weighted residual method with the Heaviside step function as the weighting function over a local sub-domain. The present approach is a truly meshless method based only on a number of randomly located nodes. In this approach, neither global background integration mesh nor domain integration is needed. Radial basis functions (RBFs) interpolation is employed in shape function and its derivatives construction for evaluating the local weak form integrals. Due to satisfaction of kronecker delta property in RBF interpolation, no special treatment is needed to impose the essential boundary conditions. In order to obtain the optimum parameters, shape parameters of multiquadrics (MQ)-RBF are tuned and studied. The leakage has a significant impact on the tidal behaviour of the confined aquifer. The numerical results of this research indicate that both tidal amplitude of groundwater head in the aquifer and the distance over which the aquifer can be disturbed by the tide are considerably reduced by leakage. The novelty of the approach is the use of a local Heaviside weight function in the LRPIM which does not need local domain integration and only integrations on the boundary of the local domains are needed. Therefore, in this research a new local Heaviside weight function has been proposed. Numerical results are presented and compared with the results of analytical solution. It is observed that the obtained results agreed very well with the results of analytical solution. The numerical results show that the use of a local Heaviside weight function in the LRPIM is highly accurate, fast and robust. It is also noticed that this novel meshless approach using MQ radial basis is very stable.     

二维Logistic映射的分岔与分形

理论分析了二维Logistic映射的分岔,并采用相图、分岔图、功率谱、Lyapunov指数和分维数计算的方法,揭示出：二维Logistic映射可按倍周期分岔和Hopf分岔走向混沌;在倍周期分岔过程中,系统在参数空间和相空间中都表现出自相似性和尺度变换下的不变性．对二维Logistic映射的吸引盆及其Mandelbrot-Julia集(简称M-J集)的研究表明：吸引盆中周期和非周期区域之间的边界是分形的,这意味着无法预测相平面上点运动的归宿;M-J集的结构由控制参数决定,且它们的边界是分形的．     

复合Logistic映射中的逆分岔与分形

利用分岔图，揭示出复合Logistic映射可按倍周期分岔走向混沌，且混沌区中存在混沌危机及逆分岔现象．同时，分析了复合Logistic映射临界点的轨道，给出了复合Logistic映射Mandelbrot-Julia集(简称M-J集)的定义，推广了Welstead和Cromer所提出的周期点查找技术，并利用该技术，构造出一系列复合Logistic映射的M-J集．在此基础上，研究了M-J集的对称性；探索了M集周期区域分布的拓扑不变性；通过定性地建立M集上J集的整体刻画，发现M集包含了J集构造的大量信息．     

On the signal degradation induced by the interpolation and the sampling rate reduction in aeroacoustics hybrid methods

Taking place within more extensive work that focuses on hybrid methods in aeroacoustics, the present study is devoted to the data transfer operations that are to be performed between two stages of a hybrid calculation. More precisely, the article focuses on two typical operations that usually accompany such data transfer, which are (i) the sampling rate reduction and (ii) the interpolation of the unsteady perturbations to be transmitted from one stage to the other. First part of the paper analyzes the two main issues of such operations, which are the spuriousing and the aliasing phenomena. For doing so, the usual notions of the interpolation theory are revisited before they are synthesized within an original approach. The here proposed formalism allows to understand better both the spuriousing and the aliasing phenomena, as well as to accurately predict the impact of the latter on the data to be transmitted in terms of signal degradation. Second part of the paper provides an illustration and a validation of these theoretical developments via a direct application to a typical aerodynamic noise problem (aeroacoustic emission by a 2D cylinder cross flow). There, it is further shown how the here proposed formalism can help in improving aeroacoustics hybrid calculations by predicting (and thus possibly minimizing) the bias to be induced on the acoustic extrapolation stage because of the aliasing and/or spuriousing effects inherited from the sampling rate reduction and/or interpolation of CFD data—which is likely to occur in any hybrid scenario. Copyright © 2012 John Wiley & Sons, Ltd.     

On the generalized Cauchy function and new Conjecture on its exterior singularities

This article studies on Cauchy’s function f (z) and its integral, (2πi)J[ f (z)] ≡■C f (t)dt/(t z) taken along a closed simple contour C, in regard to their comprehensive properties over the entire z = x + iy plane consisted of the simply connected open domain D + bounded by C and the open domain D outside C. (1) With f (z) assumed to be C n (n < ∞-times continuously differentiable) z ∈ D + and in a neighborhood of C, f (z) and its derivatives f (n) (z) are proved uniformly continuous in the closed domain D + = [D + + C]. (2) Cauchy’s integral formulas and their derivatives z ∈ D + (or z ∈ D ) are proved to converge uniformly in D + (or in D = [D +C]), respectively, thereby rendering the integral formulas valid over the entire z-plane. (3) The same claims (as for f (z) and J[ f (z)]) are shown extended to hold for the complement function F(z), defined to be C n z ∈ D and about C. (4) The uniform convergence theorems for f (z) and F(z) shown for arbitrary contour C are adapted to find special domains in the upper or lower half z-planes and those inside and outside the unit circle |z| = 1 such that the four general- ized Hilbert-type integral transforms are proved. (5) Further, the singularity distribution of f (z) in D is elucidated by considering the direct problem exemplified with several typ- ical singularities prescribed in D . (6) A comparative study is made between generalized integral formulas and Plemelj’s formulas on their differing basic properties. (7) Physical sig- nificances of these formulas are illustrated with applicationsto nonlinear airfoil theory. (8) Finally, an unsolved inverse problem to determine all the singularities of Cauchy function f (z) in domain D , based on the continuous numerical value of f (z) z ∈ D + = [D + + C], is presented for resolution as a conjecture.     

The hölder continuity of generalized solutions of a class quasilinear parabolic equations

Let A cud B satisfy the Structural conditions (2), the local Hölder continuityinterior to Q=G×(0, T) is proved for the generalized solutions of quasilinearparabolic equations as follows: u₂ - divA(x, t,u,∇u) + B(x, t, u, ∇u)=0     

The fractal structure of attractor for the generalized Kuramoto-Sivashinsky equations

I.IntroductionSincethereexistspectralbarriersandspectralgapconditions,theexistenceofaninertialmanifoldformanynonlineardissipativeevolutionequationsisstillamystery.Recently,Edenetal[5]havediscoveredthatfornonlinearsemigroup,definedbynonlineardissipativeevolutionequationsincludingZDNavier-Stokesequations,thereexistsatinliefractaldimensionalinertialsetwhichmayberepresentedbyaunionoffractillsetsandattractor,ifitisLipschitzcontinuousandissqueezingonacompacti,ositiveinvariantset.Ontileotherhand,S…     

Chebyshev approximation of the second kind of modified Bessel function of order zero

The second kind of modified Bessel function of order zero is the solutions of many problems in engineering. Modified Bessel equation is transformed by exponential transformation and expanded by J. P. Boyd‘ s rational Chebyshev basis.     

Numerical study of radial basis function interpolation for data transfer across discontinuous mesh interfaces

Allowing discontinuous or non‐matching mesh spacing across zonal interfaces within a computational domain offers many advantages, particularly in terms of easing the mesh generation process, reduction of required mesh densities, and relative motion between mesh zones. This paper presents a numerical study of a universal method for interpolating solution data across such interfaces. The method utilises radial basis functions (RBFs) for n‐dimensional volume interpolation, and treats the available solution data points simply as arbitrary clouds of points, eliminating all connectivity requirements and making it applicable to a wide range of computational problems. Properties of the developed meshless interface interpolation are investigated using analytic functions, and three issues are considered: the achievable order of spatial accuracy of the RBF interpolation alone and comparison with a variable order polynomial; the effect of a combined RBF and polynomial interpolation; and the ability of the method to recover frequency content. RBF interpolation alone is shown to achieve fourth‐order to sixth‐order spatial accuracy in one and two dimensions, and in three dimensions, using a small number of data points, third‐order and above is achievable even for a 3 : 1 discontinuous cell spacing ratio, that is a 27 : 1 volume ratio, across the interface. Hence, it is inefficient to include polynomial terms, since improving on the RBF spatial accuracy results in a significant increase in the system size and deterioration in conditioning. It is also shown that only five points per wavelength are required to capture both frequency and amplitude content of periodic solutions to less than 0.01% error.Copyright © 2013 John Wiley & Sons, Ltd.     

Accurate computation and interpolation technique of finite analytic coefficients

The finite analytic method (FA) developed in the last decade is an effective numerical method for solving fluid flow problems. However, because of the limitation in the present computer, large round-off errors are found in calculating FA coefficients when Reynolds number is large. This paper investigates the cause of this difficulty and presents a special programming technique in making an accurate computation of FA coefficients. Then a fundamental function known as "P_e" is tabulated by the accurate computation. In practical application the interpolation technique is employed so that the FA coefficients can be obtained reliably and quickly.     

The fractal research and predicating on the times series of sunspot relative number

In this paper, with the theory of nonlinear dynamic systems, It is analyzed that the dynamic behavior and the predictabilityfor the monthly mean variationsof the sunspot relative number recorded from January 1891 to December 1996. Inthe progress, the fractal dimension (D=3. 3±0.2) for the variation process wascomputed. This helped us to determine the embedded dimension [2×D+1]=7.By computing the Lyapunov index (λ₁=0.863), it was indicated that the variationprocess is a chaotic system. The Kolmogorov entropy (K=0.0260) was also computed, which provides, theoretically, the predicable time scale. And at the end, according to the result of the analysis above, an experimental predication is made,whose data was a part cut from the sample data.     

基于分形理论的煤层气非达西流动分析

在考虑了煤层的分形特征和启动压力梯度影响的基础上,建立了无限大煤层中气体低速非达西流动的数学模型,并求得了量纲为一的井底压力的Laplace空间解析解,并根据数值求解结果绘制了典型的井底压力动态曲线。研究结果表明:在定产量生产的情况下,分形维数和量纲为一的启动压力梯度对早期井底压力动态无显著影响;在生产的中后期,由于受二者的影响,压力导数曲线上的径向流水平直线段消失;量纲为一的井储系数的影响主要表现在续流阶段,而吸附系数则主要影响煤基质向裂隙扩散的过渡阶段。     

Similarities between the quasi‐bubble and the generalized wave continuity equation solutions to the shallow water equations

Two common strategies for solving the shallow water equations in the finite element community are the generalized wave continuity equation (GWCE) reformulation and the quasi‐bubble velocity approximation. The GWCE approach has been widely analysed in the literature. In this work, the quasi‐bubble equations are analysed and comparisons are made between the quasi‐bubble approximation of the primitive form of the shallow water equations and a linear finite element approximation of the GWCE reformulation of the shallow water equations. The discrete condensed quasi‐bubble continuity equation is shown to be identical to a discrete wave equation for a specific GWCE weighting parameter value. The discrete momentum equations are slightly different due to the bubble function. In addition, the dispersion relationships are shown to be almost identical and numerical experiments confirm that the two schemes compute almost identical results. Analysis of the quasi‐bubble formulation suggests a relationship that may guide selection of the optimal GWCE weighting parameter. Copyright © 2004 John Wiley & Sons, Ltd.     

基于单元相交的混合网格精确守恒插值方法

基于网格切割思想,发展了二维/三维混合网格条件下的单元相交算法,可精确计算任意两个多边形/多面体的交集。在此基础上,实现了基于单元相交(CIB/DC)的精确守恒插值算法。二维和三维验证算例表明,该方法能够保证插值过程中计算域内物理量的严格守恒,且具有比常规二阶插值更高的精度。