板问题多重网格法收敛性的低模估计

１．引言 近年来,多重网格法已成为行之有效的偏微分方程数值解法．对板问题有限元离散系统的多重网格法,也有不少的研究工作,如［４］,［５］,［１０］,［１３－１７］．在［４］,［１４－１７］中,作者讨论了Ｃ１协调元离散板问题的多重网格法,并在能量模（即 Ｈ２模）意义下获得了最优的收敛率．在［５］,［１０］中,作者讨论了非协调元离散问题的多重网格法,并在能量模意义下获得了最优的收敛率,同时在能量模意义下证明了套迭代多重网格法一阶收敛．但对板问题多重网格法的低模估计,即 Ｈ１模估计,至今尚未见研究,本文…     

非对称不定问题的非协调多重网格法

１．引言设　是Ｒｄ（ｄ＝２;３）中的有界多角形区域,α是它的边界．考虑下列模型问题此处ｆ∈ＥＬ２（Ω）,系数ＡＥ（Ｃ１（Ω））ｄ×ｄ满足下列一致椭园条件此处α０是正常数．此外假设Ｂ∈（Ｃ１（Ω）ｄ和ｃ∈Ｃ０（Ω）（［１４］）．（１．１）式的变分形式是：找ｕ∈Ｈ０（Ω）使得最近,非对称不定问题的非协调多重网格法吸引了众多的研究,详见问,［７］;［１０］．考虑非协调元多重网格的一个重要原因是混合元和非协调元之间存在着紧密的联系（详见【几问,问）．设ＦＩ是ｆｉ拟一致的Ｈ角形或矩形剖分,是由连接Ｆ'－'（…     

对流-扩散问题的特征──块中心差分法

１．引言１９８２年,Ｄｏｕｇｌａｓ和Ｒｕｓｓｅｌｌ［１］提出解对流一扩散问题的特征一差分方法,网格节点为均匀分布,求解区域为直线Ｒ．文中讨论了基于二次插值的特征一差分格式,但其近似解按离散Ｌ２模未达到最优阶误差估计．１９８８年Ｗｅｉｓｅｒ和Ｗｈｅｅｌｅｒ［２］提出解线性椭圆型和线性抛物型方程的块中心差分法,１９９１年王申林［３］讨论了解拟线性双曲型积分微分方程的块中心差分方法,其共同特点为近似解按离散的Ｌ２模达到最优阶误差估计,解的一阶导数的近似解达到超收敛误差估计．１９９３年由同顺［４］讨论了…     

一类非协调元的收敛性分析

１．引言最近,［１］在求解Ｓｔｏｋｅｓ方程时提出了一类非协调元．这类非协调元定义在矩形网格上,形式简单,自由度少,比熟知的多线性元还少．例如,对于三维问题,它们只有六个自由度,是型函数在各个表面中点处的函数值或者表面上的平均值．与一些熟知的非协调元不同之处是,由于这类非协调元简单的形式及其非常少的自由度,它们不包含任何协调的部分,但是,已有的结果表明它们具有良好的数值表现［１］,［２］中已有一些理论分析,用它们来求解晶体的微结构方程,也取得较好的结果．本文将详细分析这类非协调元的收敛性,所用的工…     

关于一类自由作业机器排序问题

文章研究文［１］中提出的加工时间依赖于机器的自由作业排序问题。Ｍ．Ｄｏｒｏｒ在［１］中提出了一个算法（算法３．４）。最近，Ａ．Ｊ．Ｖａｋｈａｒｉａ、Ｂ．Ｃａｔａｙ［２］及项思明、唐国春［３］均指出Ｍ．Ｄｏｒｏｒ的算法不是最优的。项思明和唐国春提出对这类问题在机器连续加工情形下的一种求解方法，即将排序问题化成指派问题。本文对这种解法作了简化，并回答文［３］中提出的几个问题。     

求解一类非单调线性互补问题的路径跟踪法及其计算复杂性

１．引言及记号 线性互补问题的一般形式是;求（ｘ,ｓ） 　　　　　　　 使其中 众所周知,当Ω＋非空时,单调线性互补问题可在多项式时间内求解,而且人们已经设计出了多种求解单调线性互补问题的有效的内点算法（见［１］和［７］）．然而,对于求解非单调线性互补问题的内点算法的研究可以说才刚刚开始．文［２］讨论了当Ｍ为Ｐ矩阵时问题（１）的中心路径的存在唯一性;文［３］给出了设计求解一类非单调线性互补问题的内点算法的一般框架;文［４］给出了求解一类非单调线性互补问题的一种势能函数约减法并讨论了其算法的计算复杂…     

含弥散核废料污染问题有限元方法的收敛性

１引言核废料污染的物理过程包括四个方面：（１）流体流动,（２）热量迁移,（３）主要污染元素的相混溶驱动,（４）微量元素的相混溶驱动．其数学模型为耦合拟线性抛物型方程组的初边值问题［２．５］．在前期数值分析工作中,［３,４］研究了不可压缩问题的有限元混合元等方法,［２,７,８］研究了可压缩问题的有限元混合元等方法,其收敛性分析限于仅考虑分子扩散而忽略弥散的情形．［９］研究了可压缩问题,考虑分子扩散和弥散,对原始数学模型进行了转化和变形,给出了有限元格式和Ｌ２（Ｊ;Ｈ１（））收敛性．本文进一步研究包…     

HILBERT空间中散逸动力系统一般线性方法的散逸稳定性

１．引言 １９９４年,Ｓｔｕａｒｔ与 Ｈｕｍｐｈｒｉｅｓ［４,５］首先考察了用 Ｒｕｎｇｅ－Ｋｕｔｔａ方法求解 Ｒｍ中的散逸动力系统（２．１）－（２．２）时数值解是否继承真解具有的散逸稳定性,并表明代数稳定且不可约的 Ｒｕｎｇｅ－Ｋｕｔｔａ方法是散逸稳定的且有一有界吸引集．１９９６年,本文作者［１］把这一工作推广到了两类特殊的一般线性方法．１９９７年,Ｈｉｌｌ在［３］中证明了Ａ－稳定是单支方法散逸稳定的充要条件,在［２］中又把文［４,５］的工作推广到了 Ｈｉｌｂｅｒｔ空间中的散逸动力系统（２．１）－（…     

一个第二类变分不等式的有限元逼近

本短文讨论下述第二类变分不等式（见 ［２, ４］）的有限元逼近及其误差分析：其中是平面凸多边形区域的的边界, 且而 ． 诸如热量控制问题,流体通过半可透性壁的扩散问题以及简化库仑摩擦接触问题的正则化方法等均可归为上述变分不等式（１）（见［２,３］）．在文［２］中给出了上述变分不等式的有限元逼近格式,作出了收敛性分析及误差估计．本文的目的是进一步用数值积分简化上述有限元逼近格式并改进原有的估计误差． 设Ｔｈ是的拟一致三角形部分,Ｖｈ是对应的线性元空间,且使得ｖｈ＝０在上．［２］中用数值积分代替其中 Ｍｉ…     

退化椭圆问题的最小二乘混合元逼近

１．引言我们将考虑具有退化系数的椭圆问题其中 Ω为 ＩＲ２中的一个凸多边形区域,定义为这里的ｇ＞０为分片线性连续函数．从物理背景来看,问题（１．１）来源于轴对称共振结构中的电磁场研究．Ｍａｒｉｎｉ,Ｐｉｅｔｒａ（１９９５）［４］研究了问题（１．１）的混合有限元逼近,并得到了最优误差估计． 此文,我们采用一种新的混合元,即最小二乘混合元方法［５］;对退化问题（１．１）进行逼近,利用插值投影证明了近似解具有最优阶精度的收敛性．比较起经典混合元方法来,最小二乘混合元方法有两个优越性：有限元空间不必满足Ｌ…     

Formulation of linear problems and solution by a universal machine

Using the predicate language for ordered fields a class of problems referred to aslinear problems is defined. This class contains, for example, all systems of linear equations and inequalities, all linear programming problems, all integer programming problems with bounded variables, all linear complementarity problems, the testing of whether sets that are defined by linear inequalities are semilattices, all satisfiability problems in sentenial logic, the rank-computation of matrices, the computation of row-reduced echelon forms of matrices, and all quadratic programming problems with bounded variables. A single, one, algorithm, to which we refer as theUniversal Linear Machine, is described. It solves any instance of any linear problem. The Universal Linear Machine runs in two phases. Given a linear problem, in the first phase a Compiler running on a Turing Machine generates alinear algorithm for the problem. Then, given an instance of the linear problem, in the second phase the linear algorithm solves the particular instance of the linear problem. The linear algorithm is finite, deterministic, loopless and executes only the five ordered field operations — additions, multiplications, subtractions, divisions and comparisons. Conversely, we show that for each linear algorithm there is a linear problem which the linear algorithm solves uniquely. Finally, it is shown that with a linear algorithm for a linear problem, one can solve certain parametric instances of the linear problem.Research was supported in part by the National Science Foundation Grant DMS 92-07409, by the Department of Energy Grant DE-FG03-87-ER-25028, by the United States—Israel Binational Science Foundation Grant 90-00434 and by ONR Grant N00014-92-J1142.Corresponding author.     

A double-stage piecewise recursive filter for real-time applications

This paper examines the basic problem of estimating the state of a system described by a set of dynamical equations. The state estimation is performed by means of optimal filtering techniques in which the state observation is given by a set of nonlinear equations.The basic issue encountered in most engineering problems is the availability of more observable information than the processor can effectively process in real-time, while performing full optimal filtering on all available observable data. To alleviate this deficiency, an algorithm is developed in which the rate of information processing is kept at a high value, while the system statistics are evaluated at a much slower rate.The method utilizes a combination of fast and slow filtering loops, in which the observable data are processed at a high rate during the fast loop, while the system error covariance, gain computations, and all other system statistics are processed at a lower rate in the slow loop. Methods are provided to resolve any incompatibility in the system statistics resulting from this fast-loop/slow-loop processing combination.Applications describing problems pertaining to aircraft navigation are presented. Specifically, applications to aircraft navigation through a satellite network are studied. Appropriate simulation results from the above studies are shown.This work was performed under contract with the Department of the Air Force, Space and Missile Organization (SAMSO), Los Angeles, California, Contract No. F04-701-75-C-0180.The authors are indebted to Dr. C. Johnson, GPS System Engineer, Mr. G. Consolver, GPS Software Manager, and Mr. W. Riley, Communications/Navigation Department, Texas Instruments for their continuous encouragement in the course of this work. Also, they are indebted to Major M. Birnbaum, SAMSO, for his constructive criticism on the conceptual design of this work, and to Mr. A. Bierman, Aerospace Corporation, for kindly providing simulation data of aircraft trajectories. Finally, the authors would like to express their gratitude to Dr. N. Carlson, Intermetrics, for valuable consultation during the progress of this work.Portions of this paper were presented by the senior author at the GSP Conference, Plan 76, IEEE Position, Location, and Navigation Symposium, San Diego, California, 1976.     

提高维吾尔文搜索质量的一些关键技术研究

维吾尔语和土耳其语、乌兹别克语、哈萨克语等同属于阿尔泰语系突厥语分支,是一种粘着型语言,具有丰富的构词合构形附加成分.主要分布于中国新疆维吾尔自治区,以及中西亚地区的维吾尔族使用.目前使用基于阿拉伯文字母的维吾尔文.当前用google,bing等搜索引擎搜索维吾尔文时搜索结果不包括该词语的构词合构形附加成分等其它形式.严重影响了搜索结果的准确性.通过词干提取和词缀附加生成词语组合的方法以解决提高维吾尔文搜索质量的问题.此问题解决,将在大幅度提高维吾尔文搜索的准确性,并针对提高其它突厥语族语言(土耳其语、乌兹别克语、哈萨克语等)的搜索质量、数据挖掘、数据分析、数据安全、自然语言理解等研究领域有一定的参考价值.     

绿色产品定价与市场入侵研究:基于市场权力结构的视角

为研究绿色产品定价与市场入侵问题,针对由两个制造商和一个零售商组成的系统,构建了制造商主导、零售商主导以及供需双方均势等三种权力结构模型,分析了权力结构、消费者绿色偏好等因素对绿色产品定价和市场入侵的影响。研究发现:绿色产品入侵市场存在可行条件,并且绿色产品市场入侵会导致普通产品制造商的利润损失;市场入侵对零售商有利,说明零售商会有动机引入绿色产品;绿色产品入侵市场和市场权力结构均会影响各个企业最优定价策略,其中拥有主导权的企业会选择高价策略;权力结构会影响供需双方的利润分配,以及导致系统利润的损失。     

金融风险在股票市场的传染效应及联动行为分析

近年来金融危机频发并表现出了易传染性,引起了众多学者的高度关注。以动态条件相关模型研究美欧股市与中、日、韩股市间的时变相关性,并结合内生多重结构突变模型划分危机传染阶段,选用溢出指数模型分析股市间的风险溢出特性;随后,定义股市间相互影响的联动模式并构建不同传染阶段的加权有向网络图分析股市间的联动行为。研究表明:美欧股市对中日韩股市有明显的传染效应,被传染的速度和持续时间均不相同;金融传染和风险溢出展现出一定的不一致性,危机期间日股的风险溢出效应强于美股;传染效应在联动网络中表现为联动模式的高聚类性和高联动性,相比欧债危机,次贷危机时期股市间展现出更强的联动行为;日股与美欧股市在两次危机中均表现出最强的联动性,其所受影响也最大。     

Personae at play. ‘Men of Mathematics’ in commentary

Mathematical writers, above all, Euclid, tend to present their theorems as decontextualized, abstract propositions, which has become the standard modus of textual presentation in theoretical mathematics. Mathematical commentators, however, provide their readers with personal names and historical facts in order to elucidate problems, provide contexts of discovery, or construct doxographies, among other things. Modern readers have used such information for the construction of histories of science. When we look at these passages, however, we see that personal names and information about mathematicians of the past can serve quite a range of different objectives, such as the strategic self-positioning of the commentator vis-à-vis the past or present of mathematics, the education of the reader, mathematical or moral, the construction of the history of the field, etc. Not only does the commentator present a persona of himself to the reader, he can also turn colleagues and predecessors into personae. This paper attempts to elucidate the practice, by offering four examples of such plays of and with personae, in Pappus, Eutocius, al-Nayrīzī, and Proclus.     

On the mechanical and thermomechanical properties of low-density polyethylene/ethylene-α-octene copolymer blends

Blends of low-density polyethylene (LDPE) and ethylene-octene copolymer (EOC) were obtained. The effect of EOC content and absorbed radiation dose on the mechanical and thermomechanical properties of LDPE/EOC blends are investigated. Particular attention is given to a tensile stress-strain analysis and the “form-memory” effect of the blends. With growing LDPE content, the elastic modulus, the yield stress, and the thermorelaxation and residual stresses of the blends increase, but the ultimate elongation at break decreases, which is caused by the higher crystallinity of polyethylene. As a result of radiation-induced cross-linking, the elastic modulus, the yield stress (at a 1% strain), the ultimate yield strength, and the thermorelaxation and residual stresses increase, while the ultimate elongation at break and the melt flow-behavior index decrease, which is confirmed by the growing gel fraction in the blend. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 44, No. 2, pp. 279–286, March–April, 2008.     

Representation of the solution to a model problem in semiconductor physics

We study a mixed type problem for the Poisson equation arising in the modeling of charge transport in semiconductor devices [V. Romano, 2D simulation of a silicon MESFET with a non-parabolic hydrodynamical model based on the maximum entropy principle, J. Comput. Phys. 176 (2002) 70-92; A.M. Blokhin, R.S. Bushmanov, A.S. Rudometova, V. Romano, Linear asymptotic stability of the equilibrium state for the 2D MEP hydrodynamical model of charge transport in semiconductors, Nonlinear Anal. 65 (2006) 1018-1038]. Unlike well-studied elliptic boundary-value problems in domains with smooth boundaries (see, for example, [O.A. Ladyzhenskaya, N.N. Uralceva, Linear and Quasilinear Elliptic Equations, Nauka, Moscow, 1973; D. Gilbarg, N.S. Trudinger, Elliptic Partial Differential Equations of Second Order, Springer-Verlag, Berlin, 1983]), our problem has two significant features: firstly, the boundary is not a smooth curve and, secondly, the type of boundary conditions is mixed (the Dirichlet condition is satisfied on the one part of the boundary whereas the Neumann condition on the other part). The well-posedness of the problem in Hölder and Sobolev spaces is proved. The representation of the solution to the problem is obtained in an explicit form.     

树脂基复合材料粘接的优化设计

本文在试验数据和理论分析的基础上，提出树脂基复合材料单搭接的优化设计的原则，胶粘剂的选择、被粘物的设计，粘接长度的选择、胶层厚度的选择等．认为通过优化设计可以提高粘接强度，减轻产品的重量，提高其质量．     

Wave parameter identification problem for ocean test structure data,part 2, discrete formulation

This paper deals with the solution of the wave parameter identification problem for ocean test structure data. A discrete formulation is assumed. An ocean test structure is considered, and wave elevation and velocities are assumed to be measured with a number of sensors. Within the frame of linear wave theory, a Fourier series model is chosen for the wave elevation and velocities. Then, the following problem is posed: Find the amplitudes of the various wave components of specified frequency and direction, so that the assumed model of wave elevation and velocities provides the best fit to the measured data. Here, the term best fit is employed in the least-square sense over a given time interval.At each time instant, the wave representation involves four indexes (frequency, direction, instrument, time); hence, four-dimensional arrays are required. This formal difficulty can be avoided by switching to an alternative representation involving only two indexes (frequency-direction, instrument-time); hence, standard vector-matrix notation can be used. Within this frame, optimality conditions are derived for the amplitudes of the assumed wave model.A characteristic of the wave parameter identification problem is that the condition number of the system matrix can be large. Therefore, the numerical solution is not an easy task and special procedures must be employed. Specifically, Gaussian elimination is avoided and advantageous use is made of the Householder transformation, in the light of the least-square nature of the problem and the discretized approach to the problem.Numerical results are presented. The effect of various system parameters (number of frequencies, number of directions, sampling time, number of sensors, and location of sensors) is investigated in connection with global or strong accuracy, local or weak accuracy, integral accuracy, and condition number of the system matrix.From the numerical experiments, it appears that the wave parameter identification problem has a unique solution if the number of directions is smaller than or equal to the number of sensors; it has an infinite number of solutions otherwise. In the case where a unique solution exists, the condition number of the system matrix increases as the size of the system increases, and this has a detrimental effect on the accuracy. However, the accuracy can be improved by proper selection of the sampling time and by proper choice of the number and location of the sensors.Generally speaking, the computations done for the discrete case exhibit better accuracy than the computations done for the continuous case (Ref. 5). This improved accuracy is a direct consequence of having used advantageously the Householder transformation and is obtained at the expense of increased memory requirements and increased CPU time.This work was supported by Exxon Production Research Company, Houston, Texas. This paper is based partly on Refs. 1–4.