On numerical ranges and roots

Existence of the fractional powers is established in Banach algebra setting, in terms of the numerical ranges of elements involved. The behavior of the spectra and (for Hermitian ∗-algebras satisfying some additional hypotheses) the ∗-numerical range under taking these powers also is investigated.     

Constraint incorporation in optimization

Numerical methods for solving constrained optimization problems need to incorporate the constraints in a manner that satisfies essentially competing interests; the incorporation needs to be simple enough that the solution method is tractable, yet complex enough to ensure the validity of the ultimate solution. We introduce a framework for constraint incorporation that identifies a minimal acceptable level of complexity and defines two basic types of constraint incorporation which (with combinations) cover nearly all popular numerical methods for constrained optimization, including trust region methods, penalty methods, barrier methods, penalty-multiplier methods, and sequential quadratic programming methods. The broad application of our framework relies on addition and chain rules for constraint incorporation which we develop here.     

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.     

Numerical analysis of time-dependent Boussinesq models

In this paper we analyse numerical models for time-dependent Boussinesq equations. These equations arise when so-called Boussinesq terms are introduced into the shallow water equations. We use the Boussinesq terms proposed by Katapodes and Dingemans. These terms generalize the constant depth terms given by Broer. The shallow water equations are discretized by using fourth-order finite difference formulae for the space derivatives and a fourth-order explicit time integrator. The effect on the stability and accuracy of various discrete Boussinesq terms is investigated. Numerical experiments are presented in the case of a fourth-order Runge-Kutta time integrator.     

Simplified analytic solutions and a novel fast algorithm for Yb-doped double-clad fiber lasers

In this paper, continuous wave Yb³⁺-doped double-clad fiber lasers (DCFLs) with linear-cavity are investigated theoretically and numerically using the rate equations. Under the steady state conditions, the simplified analytic solutions of Yb³⁺-doped DCFLs under considering the scattering loss are deduced in the strongly pump condition. Compared with the known analytic solutions in published literatures, our analytic solutions are more accurate, especially, at higher reflectivity of output mirror. In addition, a fast and stable algorithm based on the Newton-Raphson method is proposed to simulate numerically Yb³⁺-doped DCFLs. The results by simplified analytic solutions are in good agreement with those by the numerical simulation. Moreover, we have performed the optimization of an Yb³⁺-doped DCFL using the simplified analytic solutions and the numerical simulations, respectively.     

Preparation and characterization of ZnO nanorods from NaOH solutions with assisted electrical field

ZnO naorods on ZnO-coated seed substrates were fabricated by solution chemical method from Zn(NO₃)₂/NaOH under assisted electrical field. The working mechanism of electrical field was analyzed and the factors affecting the rod growth such as potential, precursor concentration and growth temperature were elucidated. The structural and optical properties are characterized by SEM, TEM, XRD, HRTEM and UV-vis. The results indicated that the nanorods have wurtzite structure without electrical field and are primarily of zincite structure under electrical field; when the electrical field is 1.1-1.3 V, not only the elevation of ion diffusion and adsorption lower the crystallite/solution interfacial energy and then the crystal nucleation barrier by increasing charge intensity, but also the production of H⁺ through oxidation of OH⁻ increases properly the degree of solution supersaturation near the substrate, and thus lowers the activation energy. Both the two processes do favor to rod growth. With increasing precursor concentration in this system, the average diameter and length of ZnO nanorods increase, leading to decreasing of optical transmittance. The maximum rod growth rate at given concentration of Zn²⁺ occurs at a specific temperature.     

Fuzzy矩阵方程的摄动问题

讨论了Fuzzy矩阵A的同解简化矩阵A^(2)，指出陈贻源论文《解Fuzzy关系方程》中定量3的错误。研究Fuzzy矩阵方程的摄动问题，解决了汤服成（2000）提出的未解决问题。     

On the Numerical Range of Operator Polynomials

An upper bound for the number of components of the numerical range of matrix polynomials is established. We also establish a necessary condition and a sufficient condition for the connectedness of the numerical range of a quadratic selfadjoint matrix polynomial. The boundary of the numerical range of linear matrix polynomials is also considered.     

A note on stability of the Douglas splitting method

In this note some stability results are derived for the Douglas splitting method. The relevance of the theoretical results is tested for an advection-reaction equation.

     

灰色Verhulst模型的样条插值函数的残差修正

本用样条函数对灰色Verhulst模型的残差序列进行插值拟合，然后作用于二阶线性微分方程，并以此修正原模型，得到一种新的预测模型的数值解，提高了拟合的精确度。