Newton-Cotes公式的渐近展开及其应用

给出了(n+1)-点Newton-Cotes公式的渐近展开,作为特例,得到了梯形求积公式和Simpson求积公式及其复化形式的渐近展开,分析了与Euler-Maclaurin展开的区别与联系,指出从复化Simpson求积公式出发,可得到数值积分的高阶Romberg算法.

基于两步正则化Gauss-Newton迭代算法的ECT图像重建

电容层析成像(ECT)技术求解图像重建问题属于非线性问题,并且存在严重的不适定性。为提高图像重建精度,提出了一种基于两步正则化Gauss-Newton迭代算法的ECT图像重建方法。针对标准正则化Gauss-Newton迭代算法在图像重建中存在的不收敛问题,引入了两步迭代方法;改进了正则化矩阵,提高了解估计的精确度;考虑到Gauss-Newton算法对迭代初值的依赖性,加入了同伦算法。最后,进行仿真和静态实验,并与线性反投影(LBP)算法、Landweber算法、Tikhonov正则化算法进行对比。结果表明,该方法可有效提高图像重建精度。

自适应Newton-Thiele有理插值及应用

二元连分式插值是二元有理插值的重要组成部分;文章在前人研究的基础上,对Newton-Thiele有理插值构造过程进行改进。针对Newton-Thiele有理插值在插值过程出现逆差商不存在的情况,传统的解决方法是将相应的Thiele型插值连分式转换为Newton插值多项式,然而该处理方法会导致计算复杂度的增加。借鉴相关文献在一元有理插值上的选点方法,文章给出一种带终止条件的自适应贪婪选点算法,即在给定插值点中根据自适应条件筛选出局部点对函数进行构造,以提高Newton-Thiele有理插值函数构造过程的稳定性,提升运算效率。对非线性函数的插值结果表明:该算法的插值效果较好、误差较小;同时将该算法应用到图像修复中,并与其他相关算法的修复效果进行对比,进一步验证了该算法的有效性。

On the Polygonal Wear Evolution of Heavy-Haul Locomotive Wheels due to Wheel/Rail Flexibility and Its Mitigation Measures

Wheel polygonal wear can immensely worsen wheel/rail interactions and vibration performances of the train and track,and ultimately,lead to the shortening of service life of railway components.At present,wheel/rail medium-or high-frequency frictional interactions are perceived as an essential reason of the high-order polygonal wear of railway wheels,which are potentially resulted by the flexible deformations of the train/track system or other external excitations.In this work,the effect of wheel/rail flexibility on polygonal wear evolution of heavy-haul locomotive wheels is explored with aid of the long-term wheel polygonal wear evolution simulations,in which different flexible modeling of the heavy-haul wheel/rail coupled system is implemented.Further,the mitigation measures for the polygonal wear of heavy-haul locomotive wheels are discussed.The results point out that the evolution of polygonal wear of heavy-haul locomotive wheels can be veritably simulated with consideration of the flexible effect of both wheelset and rails.Execution of mixed-line operation of heavy-haul trains and application of multicut wheel re-profiling can effectively reduce the development of wheel polygonal wear.This research can provide a deep-going understanding of polygonal wear evolution mechanism of heavy-haul locomotive wheels and its mitigation measures.

基于Polygon-Nurbs混合建模的汽车造型数字化设计

文章阐述了汽车造型数字化设计中Polygon建模与Nurbs建模的特点与适用情境。提出了一种基于Polygon-Nurbs混合建模的汽车造型数字化设计方法。该建模方法的核心流程是先采用多边形确定形体,再使用Nurbs建模构建稳定的曲面。基于该混合建模思路,进行了汽车造型设计实践,证明了该建模方法能够发挥两种建模方式的优点,快速准确构建形体,并能生成高质量的曲面模型。

Special Relativity’s “Newtonization” in Complex “Para-Space”: The Two Theories Equivalence Question

Complex model, say C3, of “para-space” as alternative to the real M4 Minkowski space-time for both relativistic and classical mechanics was shortly introduced as reference to our previous works on that subject. The actual aim, however, is an additional analysis of the physical and para-physical phenomena’ behavior as we formally transport observable mechanical phenomena [motion] to non-real interior of the complex domain. As it turns out, such procedure, when properly set, corresponds to transition from relativistic to more classic (or, possibly, just classic) kind of the motion. This procedure, we call the “Newtonization of relativistic physical quantities and phenomena”, first of all, includes the mechanical motion’s characteristics in the C3. The algebraic structure of vector spaces was imposed and analyzed on both: the set of all relativistic velocities and on the set of the corresponding to them “Galilean” velocities. The key point of the analysis is realization that, as a matter of fact, the relativistic theory and the classical are equivalent at least as for the kinematics. This conclusion follows the fact that the two defined structures of topological vector spaces i.e., the structure imposed on sets of all relativistic velocities and the structure on set of all “Galilean” velocities, are both diffeomorphic in their topological parts and are isomorphic as the vector spaces. As for the relativistic theory, the two approaches: the hyperbolic (“classical” SR) with its four-vector formalism and Euclidean, where SR is modeled by the complex para-space C3, were analyzed and compared.

Simulation of Steel Reinforcement on the Nonlinear Behaviour of Slender Glulam Beam Columns by Using the Newton-Raphson Method

The current theory in NF EN 1995-1-1/NA of Eurocode 5, which is based on maximum deflection, has been investigated on softwoods. Therefore, this theory is not adapted for slender glulam beam columns made of tropical hardwood species from the Congo Basin. This maximum deflection is caused by a set of loads applied to the structure. However, Eurocode 5 doesn’t provide how to predict this deflection in case of long-term load for such structures. This can be done by studying load-displacement (P-Δ) behaviour of these structures while taking into account second order effects. To reach this goal, a nonlinear analysis has been performed on a three-dimensional beam column embedded on both ends. Since conducting experimental investigations on large span structural products is time-consuming and expensive especially in developing countries, a numerical model has been implemented using the Newton-Raphson method to predict load-displacement (P-Δ) curve on a slender glulam beam column made of tropical hardwood species. On one hand, the beam has been analyzed without wood connection. On the other hand, the beam has been analyzed with a bolted wood connection and a slotted-in steel plate. The load cases considered include self-weight and a uniformly applied long-term load. Combinations of serviceability limit states (SLS) and ultimate limit states (ULS) have also been considered, among other factors. A finite-element software RFEM 5 has been used to implement the model. The results showed that the use of steel can reduce displacement by 20.96%. Additionally, compared to the maximum deflection provided by Eurocode 5 for softwoods, hardwoods can exhibit an increasing rate of 85.63%. By harnessing the plastic resistance of steel, the bending resistance of wood can be increased by 32.94%.

基于Newton插值与超松弛技术的铁路客运量预测研究

为提高铁路客运量预测精度,提出Newton插值法对客运量原始数据进行预处理以解决因节假日或重大事件造成的数据异常问题。另外,引入超松弛技术(OR)对铁路客运量预测结果进行修正,提出非线性递减权重改进粒子群算法以优化松弛因子。最后,将Newton插值法、超松弛技术与GM(1,1)和BP神经网络预测相结合,提出铁路客运量Newton-GM-BP-OR组合预测方法,并以北京市铁路客运量预测为例验证预测方法的有效性。研究结果表明,基于Newton插值法处理异常客运量数据的预测效果较基于原始数据序列更好,改进的粒子群算法在求解松弛因子过程中体现出更好的寻优能力和收敛速度,且超松弛技术对GM(1,1)和BP神经网络预测结果的修正也使得Newton-GM-BP-OR组合预测方法具有更高的预测精度。

基于Newton插值的具有前向安全性的可验证多秘密共享方案

基于Newton插值,提出一个新的可验证的具有前向安全性的(t,n)门限秘密共享方案。方案中,利用Newton插值法构造多项式进行秘密的分发和恢复秘密;利用椭圆曲线上的双线性对理论,验证分发者的诚实性和分发过程的有效性;方案具有前向安全性,假设敌手窃取前一时间段的秘密也无法获取任何有效信息,有效地降低了秘密泄露带来的风险。

GENERIC NEWTON POLYGON OF THE L-FUNCTION OF n VARIABLES OF THE LAURENT POLYNOMIAL I

The Hodge bound for the Newton polygon of L-functions of T-adic exponential sums associated to a Laurent polynomial is established.We improve the lower bound and study the properties of this new bound.We also study when this new bound is reached with large p arbitrarily,and hence the generic Newton polygon is determined.

非精确Newton方法中线性迭代收敛判据研究

在计算流体力学中,采用隐式时间推进方法时通常需要采用Newton类迭代方法求解大型非线性离散系统。每步非线性迭代需求解由非线性系统Jacobian矩阵组成的大型线性方程组,其中线性方程组求解误差会对非线性系统的收敛性产生显著影响,然而对存在Jacobian矩阵误差情况下的线性迭代收敛判据缺乏深入的研究。本文针对上述问题,首先给出了存在Jacobian矩阵误差和线性迭代误差情况下Newton迭代式的形式,并通过数值测试验证了Jacobian矩阵误差对迭代产生较大影响的可能性;其次对常见的两种不同类型的线性迭代收敛判据进行了数值测试,重点研究了存在Jacobian矩阵误差情况下容易产生的过度求解问题;最后,结合上述两类判据的特点发展了一种新的线性迭代收敛判据,结果表明:新提出的迭代收敛判据能够有效缓解过度求解问题,从而提高计算效率。

Drag Coefficient of a Non-Convex Polygonal Plate during Free Fall

Waterside creatures or aquatic organisms use a fin or web to generate a thrust force. These fins or webs have a non-convex section, referred to as a non-convex shape. We investigate the drag force acting on a non-convex plate during unsteady motion. We perform the experiment in a water tank during free fall. We fabricate the non-convex plate by cutting isosceles triangles from the side of a convex hexagonal plate. The base angle of the triangle is between 0° to 45°. The base angle is 0 indicates the convex hexagonal thin plate. We estimate the drag coefficient with the force balance acting on the model based on the image analysis technique. The results indicate that increasing the base angle by more than 30° increased the drag coefficient. The drag coefficient during unsteady motion changed with the growth of the vortex behind the model. The vortex has small vortices in the shear layer, which is related to the Kelvin-Helmholtz instabilities.

基于非精确Newton迭代正则化的EIT图像重构算法

图像重构在电阻抗成像技术的实际应用中起着至关重要的作用。针对EIT图像重构问题存在的不适定性和病态性,提出了一种快速的非精确Newton迭代正则化方法,该方法包括内外两层迭代:外层迭代为非精确Newton方法且以偏差原则作为停止准则;内层迭代采用迭代正则化格式并结合选取适当的步长来加速,如Landweber迭代和同伦摄动迭代等。结果表明,与传统的迭代正则化方法相比,非精确Newton迭代正则化方法大大改进了成像效率,可达到较好的图像重构质量和鲁棒性。

Differences of polygonal faults related to upper Miocene channels:a case study from the Beijiao sag of Qiongdongnan basin,South China Sea

Deep-water coarse-grained channels are embedded within a polygonal fault tier,and the polygonal faults(PFs)present non-polygonal geometries rather than classic polygonal geometry in plan view.However,PFs present differences when they encounter deep-water(coarse-grained vs.fine-grained)channels with different lithology,which has not been further studied to date.Three-dimensional(3D)seismic data and a drilling well from Beijiao sag of Qiongdongnanbasin,South China Sea were utilized to document the plan view and cross-sectional properties of the PFs and their differences and genetic mechanism were investigated.Results show that,first,PFs can be divided morphologically into channel-segmenting PFs and channel-bounding PFs in plan view.The former virtually cuts or segments the axes of channels in highand low-amplitudes,and the latter nearly parallels the boundaries of the channels.Both are approximately perpendicular to each other.Secondly,channel-bounding PFs that related to low-amplitude channels are much longer than those of high-amplitude ones;channel-segmenting PFs related to low-amplitude channels are slightly longer than the counterparts related to high-amplitude channels.Lastly,the magnitudes(e.g.,heights)of the PFs are proportional to the scales(e.g.,widths and heights)of low-amplitude channels,whereas the magnitudes of the PFs are inversely proportional to the scales of high amplitude channels.Coarse-grained(high amplitude)channels act as a mechanical barrier to the propagation of PFs,whereas fine-grained(low-amplitude)channels are beneficial to the propagation and nucleation of PFs.Additionally,the genetic mechanism of PFs is discussed and reckoned as combined geneses of gravitational spreading and overpressure hydrofracture.The differences of the PFs can be used to reasonably differentiate coarse-grained channels from fine-grained channels.This study provides new insights into understanding the different geometries of the PFs related to coarse-grained and fine-grained channels and their genetic mechanism.

Radio-Anatomical Study of Anterior Variants of the Polygon of Willis

The anterior system of the polygon of Willis plays an anastomotic role between the carotid territories and communicates with the posterior system. It ensures adequate cerebral blood flow in the event of internal occlusion. Its ability to redistribute blood flow depends on its morphology, the presence and size of the vessels that constitute it and their variants, knowledge of which is essential. The morphology of the carotid system is not constant and a number of variants are well recognized with greater hemodynamic importance. In this context, we undertook this work which aims to seek the anatomical variants of the anterior part of the PW;to determine the diameters of the vessels of the latter and to correlate these results with epidemiological data. We carried out a retrospective, descriptive and analytical study over a period of 4 months at Fann University Hospital, targeting patients who had undergone brain MRI with a 3D TOF sequence whatever the indication. Patients with lesions of cerebrovascular pathologies such as ischemia, hemorrhage, aneurysm, arteriovenous malformations or presenting a limited analysis examination due to the presence of kinetic or other artifact, were not included. Demographic data and anatomical variants were studied as well as the measurements of any continuous arterial segment greater than 0.8 mm in diameter. Those less than 0.8 mm in diameter were considered hypoplastic. The C3 segment of the internal carotid arteries, the A1 segment of the anterior cerebral arteries, the anterior communicating artery, constituted our main measurement sites. The sex ratio was 1.06, the mean age of the patients was 43 years ± 17.9 with extremes of 11 and 85 years. The anterior configuration of the polygon was complete in 57.5% representing 50% of men and 65.5% of women, in 56% of subjects under 40 years old and 58% in subjects over 40 years old. We noted a predominance of type a in 47% of patients followed by type g found in 35% of patients. The morphology of the anterior system of the Willis polygon is a function of its variants, some of which have greater hemodynamic importance than others. Its knowledge is essential for the management of cerebrovascular diseases.

Differences of Polygonal Faults with Irregularly Polygonal Geometries: A Case Study from the Changchang Sag of Qiongdongnan Basin, Northern South China Sea

Polygonal faults(PFs)generally have a classic polygonal geometry in map view.However,under the influence of tectonic faults,diapirs,channels,and slopes,the classic polygonal geometry of PFs is not preserved,demonstrating differences(different characters)in map-view 3D seismic data covering an area of 334km^(2) of the Changchang(CC)sag,are used to document the mapview and cross-sectional characteristics of PFs.These data also help investigate the irregularly polygonal geometries of PFs due to the presence of influence factors,such as transtensional faults,submarine fans,channels,diapirs/gas chimneys,and the basal slope within the lower-middle Miocene strata.Results show that various irregularly polygonal geometries of PFs can be classified into enechelon and arcuate PFs,channel-segmenting and-bounding PFs,radial PFs,and rectangular PFs in map-view.En-echelon and arcuate PFs are induced by transtensional faults and exhibit a unique‘flower’structure in NE-and SE-trending cross-sections in the NW area of the study area.This finding is documented for the first time.Channel-segmenting PFs occur in the(northwest)low-amplitude muddy channel and are inhibited in the(southeast)high-amplitude sandy channel in the SW area.Radial PFs are radially aligned around a gas chimney/diapir containing some high-amplitude anomalies(HAAs)in the middle area.The presence of intrusive sandstones with HAAs along the periphery of the diapirs restricts the occurrence of PFs.Two high-amplitude submarine fans act as a mechanical barrier to the propagation of PFs.Meanwhile,the(moderate)slope in the NE area induces rectangular PFs.Additionally,the geneses of the PFs in the current study are comprehensively discussed.This study adds to our understanding of the differences between PFs with irregularly polygonal geometries.

A Distributed Newton Method for Processing Signals Defined on the Large-Scale Networks

In the graph signal processing(GSP)framework,distributed algorithms are highly desirable in processing signals defined on large-scale networks.However,in most existing distributed algorithms,all nodes homogeneously perform the local computation,which calls for heavy computational and communication costs.Moreover,in many real-world networks,such as those with straggling nodes,the homogeneous manner may result in serious delay or even failure.To this end,we propose active network decomposition algorithms to select non-straggling nodes(normal nodes)that perform the main computation and communication across the network.To accommodate the decomposition in different kinds of networks,two different approaches are developed,one is centralized decomposition that leverages the adjacency of the network and the other is distributed decomposition that employs the indicator message transmission between neighboring nodes,which constitutes the main contribution of this paper.By incorporating the active decomposition scheme,a distributed Newton method is employed to solve the least squares problem in GSP,where the Hessian inverse is approximately evaluated by patching a series of inverses of local Hessian matrices each of which is governed by one normal node.The proposed algorithm inherits the fast convergence of the second-order algorithms while maintains low computational and communication cost.Numerical examples demonstrate the effectiveness of the proposed algorithm.

求解大规模混合互补问题的Newton-GMRES方法

针对大规模混合互补问题,借助Chen-Harker-Kanzow-Smale(CHKS)光滑逼近函数,将其转化为光滑方程组。与非精确牛顿法相结合,得到了一种求解混合互补问题的Newton-GMRES算法,证明了该算法的全局收敛性。数值实验结果表明,该算法可行有效。

Not Relying on the Newton Gravitational Constant Gives More Accurate Gravitational Predictions

The Newton gravitational constant is considered a cornerstone of modern gravity theory. Newton did not invent or use the gravity constant;it was invented in 1873, about the same time as it became standard to use the kilogram mass definition. We will claim that G is just a term needed to correct the incomplete kilogram definition so to be able to make gravity predictions. But there is another way;namely, to directly use a more complete mass definition, something that in recent years has been introduced as collision-time and a corresponding energy called collision-length. The collision-length is quantum gravitational energy. We will clearly demonstrate that by working with mass and energy based on these new concepts, rather than kilogram and the gravitational constant, one can significantly reduce the uncertainty in most gravity predictions.

基于多分辨率-多边形单元建模策略的多材料结构动刚度拓扑优化方法

利用多边形有限单元的高精度求解优势,融合多分辨率拓扑优化方法,实现粗糙位移网格条件下的高分辨率构型设计,由此提出多材料结构动刚度问题的拓扑优化方法。将多边形单元(位移场求解单元)劈分为精细的小单元,构造设计变量与密度变量的重叠网格,形成多分辨率-多边形单元的优化建模策略;以平均动柔度最小化为目标和多材料的体积占比为约束,建立多材料结构的动力学拓扑优化模型,通过HHT-α方法求解结构动响应,采用伴随变量法推导目标函数和约束的灵敏度表达式,利用基于敏度分离技术的ZPR设计变量更新方案构建多区域体积约束问题的优化迭代格式;通过典型数值算例分析优化方法的可行性和动态载荷作用时间对优化结果的影响机制。