关于AOR迭代法的研究

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

Field Theory Methods in Classical Dynamics

A Dirac picture perturbation theory is developed for the time evolution operator in classical dynamics in the spirit of the Schwinger–Feynman–Dyson perturbation expansion and detailed rules are derived for computations. Complexification formalisms are given for the time evolution operator suitable for phase space analyses, and then extended to a two-dimensional setting for a study of the geometrical Berry phase as an example. Finally a direct integration of Hamilton's equations is shown to lead naturally to a path integral expression, as a resolution of the identity, as applied to arbitrary functions of generalized coordinates and momenta.     

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

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

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

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

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

克尔介质中超连续谱的产生与自聚焦的抑制

 从（3+1）维非线性薛定谔方程出发，理论上分析了超短脉冲频谱展宽与自聚焦的影响因素。分析得出：通过改变泵浦光的功率和光束口径，可以实现光谱的极大展宽并避免自聚焦成丝。数值模拟了小口径强泵浦光束在BK7玻璃中的传输过程并进行了实验验证。模拟结果显示在超连续谱产生的同时小尺度调制被完全抑制。实验结果表明：降低泵浦光功率，使光束不会因为全光束自聚焦而发生塌陷，同时还能控制除自聚焦外的其它非线性效应，进而改善近场光束质量。由于自相位调制是超短超强脉冲产生超连续谱的重要机制之一，需要维持传输过程中的泵浦光功率，由此最佳的入射光功率应选在全光束自聚焦功率阈值附近。     

β-MnO2高压相的从头计算模拟

 基于密度泛函理论（DFT），用完全势能线性缀加平面波方法（FPLAPW）计算模拟了MnO₂的同质多相变体（金红石型结构、CaCl₂型结构、黄铁矿型结构），通过计算，预测了β-MnO₂在5 GPa时从金红石型结构转变为CaCl₂型结构，在20 GPa时进一步转变为黄铁矿型结构。另外，还总结和比较了几种金红石型结构二氧化物的压致相变特点，得出第Ⅳ主族元素的金红石型结构氧化物有很好的规律性，发生的相变序列基本一致；并随着金属阳离子半径的增大，发生相变的压力值也相应地递减，各氧化物的体积模量值相应地减小。