关于AOR迭代法的研究

本文论证了严格对角占优矩阵之AOR法的误差估计式中的误差估计常数hγ,ω（0≤γ≤ω0）的最小值是h1,1.     

Finite Element Solution of Conical Diffraction Problems

This paper is devoted to the numerical study of diffraction by periodic structures of plane waves under oblique incidence. For this situation Maxwell's equations can be reduced to a system of two Helmholtz equations in R ² coupled via quasiperiodic transmission conditions on the piecewise smooth interfaces between different materials. The numerical analysis is based on a strongly elliptic variational formulation of the differential problem in a bounded periodic cell involving nonlocal boundary operators. We obtain existence and uniqueness results for discrete solutions and provide the corresponding error analysis.     

线性分式规划最优解集的求法

本文使用多面集的表示定理，导出了线性分式规划最优解集的结构，并给出确定全部最优解的计算步骤。     

Truncated Newton method for sparse unconstrained optimization using automatic differentiation

When solving large complex optimization problems, the user is faced with three major problems. These are (i) the cost in human time in obtaining accurate expressions for the derivatives involved; (ii) the need to store second derivative information; and (iii), of lessening importance, the time taken to solve the problem on the computer. For many problems, a significant part of the latter can be attributed to solving Newton-like equations. In the algorithm described, the equations are solved using a conjugate direction method that only needs the Hessian at the current point when it is multiplied by a trial vector. In this paper, we present a method that finds this product using automatic differentiation while only requiring vector storage. The method takes advantage of any sparsity in the Hessian matrix and computes exact derivatives. It avoids the complexity of symbolic differentiation, the inaccuracy of numerical differentiation, the labor of finding analytic derivatives, and the need for matrix store. When far from a minimum, an accurate solution to the Newton equations is not justified, so an approximate solution is obtained by using a version of Dembo and Steihaug's truncated Newton algorithm (Ref. 1).This paper was presented at the SIAM National Meeting, Boston, Massachusetts, 1986.     

金属银熔体快冷过程的计算机模拟

利用分子动力学模拟技术研究了由周期性边界条件控制的500个银(Ag)原子构成的金属Ag熔体快速冷凝过程.冷却速率为8×1013K/s.模拟在FS(Finnis-Sinclair)相互作用势的基础上,通过双体分布函数、键对分析技术、键取向序等多种方法,对液银快速冷凝过程的微观结构转变特性作了分析,给出了连续快速冷凝过程中液银原子间依靠相互作用力形成的独特的微观结构图像.并考察了冷却过程中体系能量和元胞体积随温度的变化.模拟结果表明在快速冷凝过程中液Ag没有形成bcc结构的倾向.     

混沌实验数据处理及仿真

大学物理实验中混沌实验大多采用观察现象的方法进行,本实验采用蔡氏电路(Chua s circuit)产生混沌行为.在观察不同初始值条件下出现的倍周期分岔、阵发混沌、奇异吸引子等相图及现象的基础上,通过对采集数据进行处理,对负电阻伏安特性进行分段线性拟合,用功率频谱法、计算机仿真方法(龙格-库塔数值积分法)对混沌现象进行描绘,将实验数据与非线性方程组的数值解相结合,呈现出混沌现象的本质.     

准晶体中八方晶系点群的对称性与矩阵表示

从理论上对准晶体中八方晶系各点群进行了研究.运用八方晶系各点群的极赤投影图,列出了各点群的所有对称操作;填出了固有点群822的群乘表.运用坐标变换和群论在自定义的八方坐标系中,推导出八方晶系点群所有对称操作的矩阵.这32个3×3矩阵的结构是相当简洁的,它们的矩阵元只有5种可能取值:0,±1,±2.其中2是反映八方晶系准晶体所具有的准周期性的特殊无理数.     

强流脉冲离子束辐照靶及其喷发的数值研究

利用拟合实测的TEMP Ⅱ型加速器磁绝缘二极管(MID)电压波形及其焦点附近束流密度曲线，建立了Gaussian分布模型，据此计算了与靶作用的离子的能量及数量，采用Monte Carlo(MC)方法计算了沉积在靶内的能量.并以此作为热源，与流体动力学(HD)模型相结合，对不同的靶状态采用相应的状态方程，模拟计算了靶内压力演化情况; 同时对烧蚀产生的等离子体采用理想气体状态方程, 结合HD方程组, 模拟计算了喷发过程中压力的空间演化过程. 关键词： 强流脉冲离子束 Gaussian模型 HD方程 数值研究     

飞秒孤子传输演化的数值研究

利用分步傅里叶法数值模拟了飞秒孤子在光纤中的传输演化过程。对光纤中单个孤子的传输及二阶孤子的自陡峭效应和自频移效应进行了分析,指出在一定的参数取值范围内,自频移效应对二阶fs孤子传输的影响要比自陡峭效应大,占主导地位,且对自陡峭效应有一定的抑制作用。     

High-order accurate equations describing vibrations of thin bars

A method for deriving one-dimensional wave propagation equations in thin inhomogeneous anisotropic bars based on the mathematical homogenization theory for periodic media is used to obtain equations governing the longitudinal and transverse vibrations of a homogeneous circular bar. The equations are derived up to O(ε⁸) terms and take into account variable body forces and surface loads. Here, ε is the ratio of the bar’s typical thickness to the typical wavelength.