引导混沌运动到周期运动的自适应控制策略

提出一种自适应控制策略，对控制参数作线性反馈将非线性动力系统由混沌运动引导到指定的周期运动．所解决的关键问题是将反馈控制强度的确定转化为扩维相空间中的极点配置问题．给出了将该策略用于控制Ｌｏｇｉｓｔｉｃ映射和受迫Ｄｕｆｉｎｇ振子的仿真．     

Low rank update of singular values

The notion of a low rank update arises in many important applications. This paper deals with the inverse problem of updating a rectangular matrix by additive low rank matrices so as to reposition the associated singular values. The setting is analogous to the classical pole assignment problem where eigenvalues of a square matrix are relocated. Precise and easy-to-check necessary and sufficient conditions under which the problem is solvable are completely characterized, generalizing some traditional Weyl inequalities for singular values. The constructive proof makes it possible to compute such a solution numerically. A pseudo algorithm is outlined.

     

二元对角向量值有理插值的算法

首先提出了二元对角向量值有理插值问题,它包括主对角和副对角两种向量值有理插值,并分别给出了主对角线和副对角线上向量值有理插值的两种算法,即直接求系数bi,j的算法和基于Samelson广义逆所定义的特殊初等变换的矩阵算法.然后构造了在预给极点情况下求主对角线和副对角线上向量值有理插值的矩阵算法.最后给出多个数值例子说明上述算法的有效性.     

The index of triangular operator matrices

For any triangular operator matrix acting in a direct sum of complex Banach spaces, the order of a pole of the resolvent (i.e. the index) is determined as a function of the coefficients in the Laurent series for all the (resolvents of the) operators on the diagonal and of the operators below the diagonal. This result is then applied to the case of certain nonnegative operators in Banach lattices. We show how simply these results imply the Rothblum Index Theorem (1975) for nonnegative matrices. Finally, examples for calculating the index are presented.

     

Applying numerical continuation to the parameter dependence of solutions of the Schrödinger equation

In molecular reactions at the microscopic level, the appearance of resonances has an important influence on the reactivity. It is important to predict when a bound state transitions into a resonance and how these transitions depend on various system parameters such as internuclear distances. The dynamics of such systems are described by the time-independent Schrödinger equation and the resonances are modeled by poles of the S-matrix.Using numerical continuation methods and bifurcation theory, techniques which find their roots in the study of dynamical systems, we are able to develop efficient and robust methods to study the transitions of bound states into resonances. By applying Keller’s Pseudo-Arclength continuation, we can minimize the numerical complexity of our algorithm. As continuation methods generally assume smooth and well-behaving functions and the S-matrix is neither, special care has been taken to ensure accurate results.We have successfully applied our approach in a number of model problems involving the radial Schrödinger equation.     

一类极点的留数计算

本文探讨了一类比较特殊的极点类型的留数计算并从理论上给出证明.     

Piezomagnetic and piezoelectric poling effects on mode I and II crack initiation behavior of magnetoelectroelastic materials

The ferrite and ferroelectric phase of magnetoelectroelastic (MEE) material can be selected and processed to control the macroscopic behavior of electron devices using continuum mechanics models. Once macro- and/or microdefects appear, the highly intensified magnetic and electric energy localization could alter the response significantly to change the design performance. Alignment of poling directions of piezomagnetic and piezoelectric materials can add to the complexity of the MEE material behavior to which this study will be concerned with.Appropriate balance of distortional and dilatational energy density is no longer obvious when a material possesses anisotropy and/or nonhomogeneity. An excess of the former could result in unwanted geometric change while the latter may lead to unexpected fracture initiation. Such information can be evaluated quantitatively from the stationary values of the energy density function dW/dV. The maxima and minima have been known to coincide, respectively, with possible locations of permanent shape change and crack initiation regardless of material and loading type. The direction of poling with respect to a line crack and the material microstructure described by the constitutive coefficients will be specified explicitly with reference to the applied magnetic field, electric field and mechanical stress, both normal and shear. The crack initiation load and direction could be predicted by finding the direction for which the volume change is the largest. In contrast to intuition, change in poling directions can influence the cracking behavior of MEE dramatically. This will be demonstrated by the numerical results for the BaTiO₃–CoFe₂O₄ composite having different volume fractions where BaTiO₃ and CoFe₂O₄ are, respectively, the inclusion and matrix.To be emphasized is that mode I and II crack behavior will not have the same definition as that in classical fracture mechanics where load and crack extension symmetry would coincide. A striking result is found for a mode II crack. By keeping the magnetic poling fixed, a reversal of electric poling changed the crack initiation angle from θ₀=+80° to θ₀=−80° using the line extending ahead of the crack as the reference. This effect is also sensitive to the distance from the crack tip. Displayed and discussed are results for r/a=10⁻⁴ and 10⁻¹. Because the theory of magnetoelectroelasticity used in the analysis is based on the assumption of equilibrium where the influence of material microstructure is homogenized, the local space and temporal effects must be interpreted accordingly. Among them are the maximum values of (dW/dV)_max and (dW/dV)_min which refer to as possible sites of yielding and fracture. Since time and size are homogenized, it is implicitly understood that there is more time for yielding as compared to fracture being a more sudden process. This renders a higher dW/dV in contrast to that for fracture. Put it differently, a lower dW/dV with a shorter time for release could be more detrimental.     

MATLAB辅助数字信号处理课程教学的分析研究

数字信号处理是一门理论性和实践性都很强的课程。以离散时间系统的极零分析为例,内容抽象、不易理解。如将MATLAB软件的滤波器设计工具引入课堂教学过程中,可对教学效果差、学生对知识点理解困难等问题有所改善,进一步提高教学质量。     

Pole movement in electronic and optoelectronic oscillators

An RLC circuit with poles on the left half of the complex frequency plane is capable of executing transient oscillations. During this period, energy conversion from potential to kinetic and from kinetic to potential continuously goes on, until the stored energy is lost in dissipation through the resistance. On the other hand, in an electronic or opto-electronic oscillator with an embedded RLC circuit, the poles are forcibly placed on the right-half plane (RHP) and as far as practicable away from the imaginary axis in order to help the growth of oscillation as quickly as possible. And ultimately, it is imagined that, like the case of an ideal linear harmonic oscillator, the poles are frozen on the imaginary axis so that the oscillation neither grows nor decays. The authors feel that this act of holding the poles right on the imaginary axis is a theoretical conjecture in a soft or hard self-excited oscillator. In this article, a detailed discussion on pole movement in an electronic and opto-electronic oscillator is carried out from the basic concept. A new analytical method for estimating the time-dependent part of the pole is introduced here.