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.     

混沌实验数据处理及仿真

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

非傍轴平顶高斯光束M2因子两种定义的比较研究

基于功率密度的二阶矩方法,推导出了非傍轴平顶高斯(FG)光束束宽和远场发散角的解析表达式·研究表明,当w0/λ→0时,远场发散角趋于渐近值θmax=63.435°,与阶数无关·使用非傍轴高斯光束代替傍轴高斯光束作为理想光束,研究了非傍轴FG光束的M2因子,并与传统定义的M2因子作了比较·在非傍轴范畴,非傍轴FG光束的M2因子不仅与阶数N有关,而且与w0/λ有关·按照定义,当w0/λ→0时,非傍轴FG光束的M2因子不等于0,对阶数N=1,2,3时,M2因子分别趋于0.913,0.882和0.886·当N→∞时,M2因子取最小值M2min=0.816·     

Finding the response matrix to the external action from a subspace for a discrete linear stochastic dynamical system

For a discrete linear stochastic dynamical system, computation of the response matrix to the external action from a subspace using given observational data is examined. An algorithm is proposed and substantiated that makes it possible to improve the numerical accuracy and to reduce the amount of observational data compared to the general case where an arbitrary external action is allowed. As an illustration, a discrete system arising in the analysis of a linear stochastic dynamical continuous-time system is considered more thoroughly. Some numerical results are presented.     

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)谱的分析讨论了这类空间电荷型宏观电偶极子，及与其补偿的空间电荷热退极化的电荷动态特性；阐明了这两类俘获电荷的能阱分布，即构成宏观电偶极子的位于孔洞上下介质层内的等值异号空间电荷分别被俘获在深、浅两种能值陷阱内，而位于薄膜表面层的注入空间电荷则被俘获在中等能值陷阱中. 关键词： 反极性电晕补偿充电法 铁电驻极体 充电电流 热刺激放电     

提高微晶硅薄膜太阳电池效率的研究

采用甚高频等离子体增强化学气相沉积技术制备了系列微晶硅薄膜太阳电池，指出了气体总流量和背反射电极的类型对电池性能参数的影响.电池的I-V测试结果表明：随反应气体总流量的增加，对应电池的短路电流密度、开路电压和填充因子都有很大程度的提高，结果使得电池的光电转换效率得以提高.另外，ZnO/Ag/Al背反射电极能明显提高电池的短路电流密度，进而也提高了电池的光电转换效率.对气体总流量和背反射电极类型影响电池效率的原因进行了分析. 关键词： 微晶硅薄膜太阳电池 气体流量 ZnO/Ag/Al背反射电极     

非傍轴平顶高斯光束M2因子两种定义的比较研究

基于功率密度的二阶矩方法，推导出了非傍轴平顶高斯(FG)光束束宽和远场发散角的解析表达式.研究表明，当w0/λ→0时，远场发散角趋于渐近值θmax=63.435°，与阶数无关.使用非傍轴高斯光束代替傍轴高斯光束作为理想光束，研究了非傍轴FG光束的M2因子，并与传统定义的M2因子作了比较.在非傍轴范畴，非傍轴FG光束的M2因子不仅与阶数N有关，而且与w0/λ有关.按照定义，当w0/λ→0时，非傍轴FG光束的M2因子不等于0，对阶数N=1, 2, 3时，M2因子分别趋于0.913,0.882和0.886.当N→∞时，M2因子取最小值M2min=0.816.