非线性等式与不等式问题的信赖域算法

５１．引言 本文讨论如下非线性等式与不等式问题的求解问题：其中,吐。）是＊”＋Ｂ的连续可微函数． 非线性问题（１．１）广泛出现于优化问题、互补问题及变分不等式等应用数学领域中,因而对其算法的研究也是计算数学的一个重要部分．许多作者考察了该问题的求解,如Ｄｅｎｎｉｓ和Ｓｃｈｎａｂｅ［６］, Ｐｓｈｅｎｉｃｈｌｌｙｉ［９］, Ｒｏｂｉｎｓｏｎ［１０１等人讨论了该问题的牛顿法求解,梯度法求解,Ｂｕｒｋｅ和 Ｈａｎ［‘］讨论了高斯一牛顿法求解． Ｄｅｎｎｉｓ在［５］中首次采用了信赖域算法的求解,［５］中通过引…     

基于图的任意域内点集的Delaunay三角剖分算法

本文提出了一种基于图的二维任意域内点集的Ｄｅｌａｕｎａｙ三角剖分算法．该算法首先求出任意域内点集的约束最小生成树，然后逐次加入一边构造三角形网格，最后通过局部优化变换，得到二维任意域内点集的Ｄｅｌａｕｎａｙ三角剖分．本文还给出了该算法在有限元网格自动生成过程中的应用.     

平面及空间区域渐变无结构网格的自动生成

１．前 言 网格生成是许多数值计算首先要解决的问题．规则区域上的均匀网格,比较容易生成．但许多工程实际问题,求解区域边界形状极不规则,且由于物理参数的剧烈变化及解的性态复杂,对求解区域的网格疏密变化有某些特殊的要求．如何在非规则区域上自动生成符合使用者特殊要求的网格,是科学计算过程中人们所关心的问题． 本文基于 Ｄｅｌａｕｎａｙ三角划分原理,在二维（三维）区域给出一种自动生成渐变无结构三角形（四面体）网格的方法．对事先给定的指定结点集合和对应的间隔值集合,算法将首先自动生成全部边界结点和内部结点,然…     

平面弹性方程外问题的非重叠型区域分解算法

１．引言 区域分解算法是八十年代兴起的偏微分方程求解新技术．基于有限元法的区域分解算法对求解有界区域问题行之有效［２,４,９］．边界元方法则是处理无界区域问题的强有力的工具［１,１０,１７］,有限元与边界元耦合法得到广泛应用 ［３,５,７］．近年又发展了基于自然边界归化的区域分解算法,特别适用于无界区域问题［８,１１,１２］．迄今这方面的文章主要是针对二维Ｐｏｉｓｓｏｎ方程及双调和方程的［１３－１６］． 本文讨论平面弹性方程的Ｄｉｒｉｃｈｌｅｔ外边值问题其中Ω是充分光滑闭曲线Г０之外的无界区域,ｕ…     

2－DRoesser模型的静态干扰解耦

本文讨论了２－ＤＲｏｅｓｓｅｒ模型［１］（ＲＭ）的静态干扰解耦问题［２］（简称为２－ＤＤＤＰ），即寻求２－Ｄ状态反馈使相应的闭环系统具有抗干扰的能力，得到了问题有解的充分条件和计算相应反馈阵的算法．     

一种常用的二维任意域的Delaunay三角剖分算法的健壮性补充

由于对任意给定的平面点集通过Ｄｅｌａｕｎａｙ三角剖分进行处理可得到具有整体最优性的三角形网络,因而该方法得到了广泛的重视。但研究发现,常用的二维任意域Ｄｅｌａｕｎａｙ三角剖分算法＾［１,２］是有缺陷的,它在构成Ｄｅｌａｕｎａｙ三角形候选点的选择过程中,可以使候选点出现“位置违约”的错误,即在候造节点链表中,虽然可出现依据算法的判据有条件成为Ｄｅｌａｕｎａｙ三角形的构成点,但采用该点构成Ｄｅｌａｕｎ     

一个m维整数瓶颈运输问题及其算法

５１．引言 整数瓶颈问题 （ＩＢＰ的研究最早可追朔到 １９５９年Ｇｒｏｓｓ关于瓶颈指派问题 （ＢＡＰ的研究［’］以及 Ｂａｒｓｏｗ关于时间运输问题（ＴＴＰ）的研究［’］．对 ＢＡＰ相继又出现了改进算法［３］、图算法［４］以及将其化为一般指派问题的算法［’］． １９６９年以来， Ｈａｍｍｅｒ等人又研究了 ＴＴＰ［６－８１．近十几年来，对 ＩＢＰ仍有一些新成果面世 ［９－１５１．但无论国内或者国外，所有成果其目标函数均局限在１，２维的简单情形．本文提出并研究了目标函数为。维的整数瓶颈运输问题， ５２．问题的提出…     

任意平面域渐变三角形网格的自动划分

结合前沿生成法和Ｄｅｌａｕｎａｙ三角化方法的优点,利用节点间距函数来控制区域内网格尺寸变化,并优先处理前沿上的最长边,尽可能在局部生成边长逐渐减小的Ｄｅｌａｕｎａｙ三角形,最终实现区域内网格的疏密过渡。     

一种基于图的平面点集Delaunay三角剖分算法

本文提出了一种基于图的平面点集Ｄｅｌａｕｎａｙ三角剖分算法。该算法首先求出平面点集的欧几里得最小生成树，然后逐次加入一边构造三角形网格，最后按最小内角最大的三角化准则，通过局部变换，得到平面点集的Ｄｅｌａｕｎａｙ三角剖分。本文同时阐述了它的对偶图；平面点集的Ｖｏｒｏｎｏｉ图的概念和性质。     

DIB位图影像方件的生成

本文在论文［１］、［２］的基础，论述利用扫描获得的影像参数和影像数据生成ＷＩＮ－ＤＯＷＳ ＤＩＢ格式影像文件的方法。     

一种基于逐点插入Delaunay三角剖分生成Voronoi图的算法

采用改进的逐点插入算法生成Voronoi图。该算法在逐点插入的过程中生成凸壳,进而生成Delaunay三角剖分。在生成Voronoi图的实现过程中,通过遍历三角形的边顶点快速识别相关的三角形组,进而生成Voronoi图。试验结果表明,该算法能实现,成功生成Voronoi图。     

Efficient minimum spanning tree construction without Delaunay triangulation

Given n points in a plane, a minimum spanning tree is a set of edges which connects all the points and has a minimum total length. A naive approach enumerates edges on all pairs of points and takes at least Ω(n2) time. More efficient approaches find a minimum spanning tree only among edges in the Delaunay triangulation of the points. However, Delaunay triangulation is not well defined in rectilinear distance. In this paper, we first establish a framework for minimum spanning tree construction which is based on a general concept of spanning graphs. A spanning graph is a natural definition and not necessarily a Delaunay triangulation. Based on this framework, we then design an O(nlogn) sweep-line algorithm to construct a rectilinear minimum spanning tree without using Delaunay triangulation.     

Automatic mesh generation for complex three-dimensional regions using a constrained Delaunay triangulation

This paper describes a method for generating tetrahedral meshes. The algorithm, based on the Delaunay triangulation, can treat objects of essentially arbitrary complexity. In order to preserve the surface triangulation of solid objects, it is necessary to override the Delaunay property and redefine the triangulation when points are introduced that are close to solid boundaries. Details of this constrained Delaunay algorithm are presented and an efficient implementation of the triangulation method is described. Techniques for controlling the distribution of mesh points and tetrahedron quality are also discussed.     

Ed带权点集的Regular三角化的构造算法

该文提出了一种构造带权点集的Regular三角化的算法，此算法统一了构造点集的Delauany三角化的Bowyer/Waston算法，如果采用一种称为Delauany树的数据结构来辅助点的定位，则算法的效率为O(nlogn n^d[d/s]。与Edelsbrunner和Shah提出的局部变换算法相比，此算法从理论和实现两方面都要简单一些，算法可以应用到曲线曲面重构和分子建模等领域。     

Delaunay三角网高效构建及地形仿真应用

针对基于离散点的Delaunay三角网构建过程中待插入点的定位耗时问题,提出Delaunay三角网高效构建算法,并将其用于三维地形仿真应用中。对大量数据点进行分块排序预处理后,运用空间自相关理论使下一个待插入点总是紧邻新近插入点,融合最短路径定位算法和三角形面积法,结合三角形重心与点、有向线段的关系遍历三角形,减少遍历时间。在对三角网进行LOP局部优化时,采用Delaunay四叉树保存待调整的所有边的节点信息,提高遍历效率。实验结果证明,该算法构建的三维地表真实感较强,并且具有较低的时间复杂度。     

Delaunay三角网通用合并算子及分治算法的简化

Delaunay三角网在未来地学数值模拟中将发挥重要作用。分治算法是一种著名的经典构网算法,但其子网合并过程十分复杂,限制了其应用。提出使用通用算子的概念,并用从以往算法中独立出来的算子和3个新算子来简化分治算法的子网合并。扩展三角形算子用于构造每个新三角形并维护三角网的拓扑关系和边界链表。凹边界填充算子对边界链表用递归来自动完成凹边界的智能三角形填充。子网合并算子先用一个新三角形连接两个子三角网,再合并边界链表,调用凹边界填充算子填充子网间的缝隙区域。所有算子都基于有向边的数据结构和用链表管理的三角网外边界,借助链表操作,使算法的构建简洁而又高效。除分治法外,这些算子还被成功用于构建其他算法。由随机点集以及LiDAR点云的测试表明,所有算法的构网均准确无误且分治算法的执行效率较高。     

A Flexible Algorithm for Generating All the Spanning Trees in Undirected Graphs

In this paper we propose an algorithm for generating all the spanning trees in undirected graphs. The algorithm requires O (n+m+ τ n) time where the given graph has n vertices, m edges, and τ spanning trees. For outputting all the spanning trees explicitly, this time complexity is optimal. Our algorithm follows a special rooted tree structure on the skeleton graph of the spanning tree polytope. The rule by which the rooted tree structure is traversed is irrelevant to the time complexity. In this sense, our algorithm is flexible. If we employ the depth-first search rule, we can save the memory requirement to O (n+m). A breadth-first implementation requires as much as O (m+ τ n) space, but when a parallel computer is available, this might have an advantage. When a given graph is weighted, the best-first search rule provides a ranking algorithm for the minimum spanning tree problem. The ranking algorithm requires O (n+ m + τ n) time and O (m+ τ n) space when we have a minimum spanning tree. Received January 21, 1995; revised February 19, 1996.     

改进的自连接Delaunay三角网生成算法

本文提出了一个改进的自连接Ｄｅｌａｕｎａｙ三角网生成算法。在原算法的基础上引进了封闭点的概念,在三角网的生成过程中动态地剔除封闭点,从而大大加快了生成新三角形时对点的查找过程。其次,通过对边扩展过程的研究,发现对除了第一个三角形以外的其他三角形实际上只有两条可扩展的边,新算法对这一点也进行了改进。最后,给出了实验的结果数据。     

A near optimal isosurface extraction algorithm using the span space

Presents the “Near Optimal IsoSurface Extraction” (NOISE) algorithm for rapidly extracting isosurfaces from structured and unstructured grids. Using the span space, a new representation of the underlying domain, we develop an isosurface extraction algorithm with a worst case complexity of o(√n+k) for the search phase, where n is the size of the data set and k is the number of cells intersected by the isosurface. The memory requirement is kept at O(n) while the preprocessing step is O(n log n). We utilize the span space representation as a tool for comparing isosurface extraction methods on structured and unstructured grids. We also present a fast triangulation scheme for generating and displaying unstructured tetrahedral grids     

Generation of tetrahedral mesh of variable element size by sphere packing over an unbounded 3D domain

This paper describes an algorithm for the generation of tetrahedral mesh of specified element size over an unbounded three-dimensional domain. Starting from an arbitrary point in space (defined as the origin) and guided by the concept of advancing front, spheres of size compatible with the specified element size are packed tightly together one by one to form a cluster of spheres of different sizes. The compactness of the cluster of spheres is achieved by packing spheres at a site closest to the origin in a densest manner with tangent but no overlapping with as many other spheres as possible. In view of these criteria, a rotational mechanism between spheres is innovated, which allows the newly inserted sphere to follow the path by rotation between existing spheres until the lowest point is reached. The centres of the packing spheres provide ideal locations for Delaunay point insertion to form a triangulation of tetrahedral elements of size compatible with the specified value. Spheres of random size distribution or of size specified by a node spacing function are packed by the proposed algorithm. In all test examples, high-quality tetrahedral elements of size consistent with the specified node spacing are generated. The process is fast and robust, and the time complexity for mesh generation is expected to be almost linear; however, in the present implementation, a quasi-linear time relationship is observed as a search for the nearest node on the front is proposed in the sphere packing process. The memory requirement has been kept to a minimum as no additional data structure other than the adjacency relationship of the tetrahedral elements as required by the Delaunay triangulation is stored. The structure of the advancing frontal surface needs not be explicitly constructed nor updated as it is simply not required in the sphere packing process.