具有生产性与福利性公共支出的随机内生增长模型

文章建立一个随机内生增长模型来阐明主要政策参数对经济增长与社会福利的影响.若对生产函数、效用函数、偏好及随机干扰作一些特殊的假设,我们证明了主要政策参数的均衡值能被模型参数唯一决定.进一步我们还得到了期望增长率与储蓄的清晰解.文章的最后,我们证明了政府支出直接影响个体决策者的决策:即提高经济增长率将减少福利;反之,增加福利将减少增长率.     

How to assess and report the performance of a stochastic algorithm on a benchmark problem: mean or best result on a number of runs?

Some authors claim that reporting the best result obtained by a stochastic algorithm in a number of runs is more meaningful than reporting some central statistic. In this short note, we analyze and refute the main argument brought in favor of this statement.     

An optimal design of collateralized mortgage obligation with PAC-companion structure using dynamic cash reserve

This paper presents a model for optimally designing a collateralized mortgage obligation (CMO) with a planned amortization class (PAC)-companion structure using dynamic cash reserve. In this structure, the mortgage pool’s cash flow is allocated by rule to the two bond classes such that PAC bondholders receive substantial prepayment protection, that protection being provided by the companion bondholders. The structure we propose provides greater protection to the PAC bondholders than current structures during periods of rising interest rates when this class of bondholders faces greater extension risk. We do so by allowing a portion of the cash flow from the collateral to be reserved to meet the PAC’s scheduled cash flow in subsequent periods. The greater protection is provided by the companion bondholders exposure to interest loss. To tackle this problem, we transform the problem of designing the optimal PAC-companion structure into a standard stochastic linear programming problem which can be solved efficiently. Moreover, we present an extended model by considering the quality of the companion bond and by relaxing the PAC bondholder shortfall constraint. Based on numerical experiments through Monte Carlo simulation, we show the utility of the proposed model.     

关于三维各向同性谐振子径向矩阵元计算的讨论

根据广义Laguerre多项式的数学性质，导出了较为简单的三维各向同性谐振子径向矩阵元的普遍公式,并在这基础上计算了一些重要特殊情形的径向矩阵元: 矢径r整数次幂的平均值，电偶极跃迁矩阵元和电四极跃迁矩阵元. 关键词： 三维各向同性谐振子 径向矩阵元 广义Laguerre多项式 偶极跃迁 四极跃迁     

A parallel inexact newton method for stochastic programs with recourse

A parallel inexact Newton method with a line search is proposed for two-stage quadratic stochastic programs with recourse. A lattice rule is used for the numerical evaluation of multi-dimensional integrals, and a parallel iterative method is used to solve the quadratic programming subproblems. Although the objective only has a locally Lipschitz gradient, global convergence and local superlinear convergence of the method are established. Furthermore, the method provides an error estimate which does not require much extra computation. The performance of the method is illustrated on a CM5 parallel computer.This work was supported by the Australian Research Council and the numerical experiments were done on the Sydney Regional Centre for Parallel Computing CM5.     

Model-dependent nature of blowtorch effect on the escape and equilibration rates

We analyze the relaxation behavior of a bistable system when the background temperature profile is inhomogeneous due to the presence of a localized hot region (blowtorch) on one side of the potential barrier. Since the diffusion equation for inhomogeneous medium is model-dependent, we consider two physical models to study the kinetics of such system. Using a conventional stochastic method, we obtain the escape and equilibration rates of the system for the two physical models. For both models, we find that the hot region enhances the escape rate from the well where it is placed while it retards the escape rate from the other well. However, the value of the escape rate from the well where the hot region is placed differs for the two models while that of the escape rate from the other well is identical for both. This work, for the first time, gives a detailed report of the similarities and differences of the escape rates and, hence, exposes the common and distinct features of the two known physical models in determining the way the bistable system relaxes. Received 25 September 2001     

高功率振荡器锁相的粒子模拟

用粒子(PIC,Particle-in-Cell)模拟方法研究了当前颇感兴趣的高功率振荡器的锁相问题。给出了模拟计算结果。讨论了锁相(Phase Locking)、Priming和InjectionSeeding等概念的区别。     

Classical limit of relativistic quantal system with attractive Coulomb interaction

Using the appropriate harmonic oscillator states and reasonable approximations, we construct coherent wavepackets corresponding to the solutions of the Klein-Gordon equation for the attractive potentialV(r)=−k/r, k>0, in two and three space dimensions. We deduce the corresponding classical limit in two dimension by requiring that the expectation value 〈r〉 of the radial variable is large. In the case of three dimensions, besides the condition of large 〈r〉, we make the uncertainty Δr=[〈r ²〉 − 〈r〉²]^1/2 a minimum with respect to certain parameter of the wavepacket. We then investigate the trajectory traversed by the wavepacket in the classical limit. We find that the classical limit of this relativistic quantal problem gives, in the leading order, the same expression for the rate of motion of the perihelion as that given by the solution of the corresponding special relativistic classical dynamical problem. We also briefly discuss some of the subtle aspects of the classical limit of the relativistic quantal system, in general.     

New efficient numerical procedures for solving stochastic variational problems with a priori maximum pointwise error estimates

In a previous paper we gave a new formulation and derived the Euler equations and other necessary conditions to solve strong, pathwise, stochastic variational problems with trajectories driven by Brownian motion. Thus, unlike current methods which minimize the control over deterministic functionals (the expected value), we find the control which gives the critical point solution of random functionals of a Brownian path and then, if we choose, find the expected value.This increase in information is balanced by the fact that our methods are anticipative while current methods are not. However, our methods are more directly connected to the theory and meaningful examples of deterministic variational theory and provide better means of solution for free and constrained problems. In addition, examples indicate that there are methods to obtain nonanticipative solutions from our equations although the anticipative optimal cost function has smaller expected value.In this paper we give new, efficient numerical methods to find the solution of these problems in the quadratic case. Of interest is that our numerical solution has a maximal, a priori, pointwise error of O(h3/2) where h is the node size. We believe our results are unique for any theory of stochastic control and that our methods of proof involve new and sophisticated ideas for strong solutions which extend previous deterministic results by the first author where the error was O(h²).We note that, although our solutions are given in terms of stochastic differential equations, we are not using the now standard numerical methods for stochastic differential equations. Instead we find an approximation to the critical point solution of the variational problem using relations derived from setting to zero the directional derivative of the cost functional in the direction of simple test functions.Our results are even more significant than they first appear because we can reformulate stochastic control problems or constrained calculus of variations problems in the unconstrained, stochastic calculus of variations formulation of this paper. This will allow us to find efficient and accurate numerical solutions for general constrained, stochastic optimization problems. This is not yet being done, even in the deterministic case, except by the first author.     

A new formula for the saddle-to-scission time

We propose a new formula for the saddle-to-scission time that is more general that the one based on Kramers' approach. Its validity and applicability is then studied in detail. Such a formula is useful for the evaluation of the fission time of very heavy nuclei.