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方程 数值研究     

改善光纤倒像器的对比度传递特性

光纤倒像器是一种特殊类型的光纤板，可将传递图像直接倒转180°并应用于像增强器。由于光纤倒像器的特殊结构，造成光纤随着距扭转轴心的距离增加，数值孔径和光通量逐渐下降。相对于光纤板，光纤倒像器在透过率和对比度传递特性上均有一定程度的下降。通过改进玻璃系统，在一定程度上降低芯皮玻璃间离子扩散和相互渗透的程度，可提高光纤倒像器的实际数值孔径；通过调节EMA吸收量，既能满足像增强器对光纤倒像器荧光屏透过率的要求，又能保证最大程度地吸收光纤中逸出的杂散光，提高对比度。选择合理的芯皮比，不仅可弥补因扭转拉伸造成的皮层厚度减薄，还可进一步增加皮层厚度，抑制光从光纤中逸出。通过以上改进，可改善光纤倒像器的对比度传递特性。     

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.     

不完全信息群体决策专家权重的集结

本文对于属性权重信息和属性效用信息都不完全的群体多属性决策问题，通过构造属性值区间和运用系统聚类分析法，对群体决策中的专家进行分类，并确定每位专家的权重．     

Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

聚丙烯孔洞铁电驻极体膜的电极化及其电荷动态特性

根据栅控恒压电晕充电组合反极性电晕补偿充电法的实验结果计算出铁电驻极体的极化强度.结果说明，伴随着薄膜内孔洞气体的Paschen击穿，该铁电体的极化强度随栅压增加而显著上升.利用上述充电方法和热刺激放电(TSD)谱的分析讨论了这类空间电荷型宏观电偶极子，及与其补偿的空间电荷热退极化的电荷动态特性；阐明了这两类俘获电荷的能阱分布，即构成宏观电偶极子的位于孔洞上下介质层内的等值异号空间电荷分别被俘获在深、浅两种能值陷阱内，而位于薄膜表面层的注入空间电荷则被俘获在中等能值陷阱中. 关键词： 反极性电晕补偿充电法 铁电驻极体 充电电流 热刺激放电     

一类相空间中的准几率分布函数系

定义了一类相空间中的准几率分布函数系，这个准几率分布函数系直接建立在具有更加广泛意义的量子相空间Schr?dinger方程解的基础之上，其中定义_α=αp-i?q和_α=(1-α)q+i?p.发现了两个有趣的关系.(1)建立的量子相空间Schr?dinger方程的解实际上是对函数φ(λ)exp［i(1-α)qp］做窗口Fourier变换.(2)这个窗口函数g(λ)起着选择窗口形式的作用，而且不同的窗口对应着不同的分布函数.当g(λ)是一个代表Gauss窗的Gauss函数的时候，准几率分布函数就是一个类似于Husimi的分布函数f^HL_α(q,p)；当g(λ)是一个表示椭圆的复函数时，准几率分布函数就是一个椭圆分布函数f^E_α(q,p)；再在g(λ)为复函数的基础上附加α=0，就可得到标准序分布函数f^S(q,p)、反标准序分布函数f^AS(q,p)和Wigner分布函数f^W(q,p)，此时g(λ)表示高度为1/12π?而长度为λ的矩形窗. 关键词： 窗口Fourier变换 相空间 Wigner分布函数