Phylogenetic analyses of Saussurea sect. Pseudoeriocoryne (Asteraceae : Cardueae) based on chloroplast DNA trnL-F sequences

The genus Saussurea is distributed mainly in the temperate and subarctic regions of Eurasia and consists of about 300 species classified into six subgenera and 20 sections. Sect. Pseudoeriocoryne in the subgenus Eriocoryne comprises four species, and is delimited mainly by acaulescence and an inflorescence with congested capitula surrounded by a rosette of leaves. All of these species are endemic to the and Qinghai-Tibet Plateau. Sequences from the chloroplast DNA trnL-F region were obtained for the four species in this section and 26 other species from four subgenera of Saussurea to resolve phylogenetic relationships among these species and to determine whether the shared characters that define sect. Pseudoeriocoryne are synapomorphic or were acquired by convergent evolution. The resulting phylogenies indicated that Saussurea sect. Pseudoeriocoryne as traditionally defined does not constitute a monophyletic group and that each of its species belongs to separate clades. Furthermore, none of these species showed a close relationship with the other species of subgenus Eriocoryne. Our results further indicated that none of the investigated subgenera are monophyletic, and that species from different subgenera clustered together. All these conclusions are provisional and their confirmation would require stronger phylogenetic support. Two possible explanations are suggested for low sequence divergence, poor resolution of internal clades and clustering of species with the rather distinct morphology of Saussurea detected in the present study. The first is rapid radiation and diversification triggered by fast habitat fragmentation due to the recent lifting of the Qinghai-Tibet Plateau and the Quaternary climate oscillations. This could have led to rapid morphological divergence while sequences diverged very little, and also caused the convergent acquisition of similar characteristics in unrelated lineages due to similar selection pressures. The second possible explanation is that both introgressive hybridization and reticulate evolution might have caused the transferring of cpDNA sequences between morphologically dissimilar species, thus leading to homogenization of sequences between lineages. (C) 2004 Elsevier Ltd. All rights reserved.

转子系统油膜振荡的小波包分解与频带能量比例特征分析

在介绍小波包分解原理的基础上 ,对试验测得的单盘单跨转子系统的油膜振荡非平稳信号用小波包分解方法进行了研究。采用db44小波基函数进行 4层小波包分解。给出了各频带内分解信号的特点及频带能量比例 ,其中第 3频带是该转子系统在 960 0r/min时产生油膜振荡的特征频带。得到的试验数据及其分析结果对转子系统油膜振荡研究和旋转机械状态监测等具有重要意义。

一种基于运动状态反馈的机器人接触力控制

针对机器人从自由空间运动控制过渡到约束空间力控制的过程中，存在冲撞、振荡甚至不稳定等问题，提出了利用加速度反馈为力控制提供阻尼，克服单纯依靠速度反馈的局限，进而稳定力控制系统的方法．这种方法不增加系统的复杂性，易于实现．本文对控制系统结构及控制策略作了分析．在一个直接驱动关节上的实验结果证明了本文分析的正确性，与只采用速度反馈的接触力控制相比，显示了这种方法的有效性。

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

中新世秦安黄土一古土壤序列记录的古气候变化周期

Analysis of periodic oscillations of climate is very important in understanding the behavior of the climate system. Milankovitch hypothesis, which holds that the glacial-interglacial climatic cycles during the Quaternary were primarily driven by variations in orbital parameters, has been supported by substantial geological evidence. Continuous long-term and high-resolution records are crucial to detect how variations of Earth's orbital parameters affected climate before the Quaternary when the boundary conditions were significantly different. Qinan loess formed in the Miocene is nearly continuous aeolian deposit in northern China. Previous study has established a constrained chronology, which provides a basis to examine long-term climatic variations. One of important issues to untangle the mechanisms behind major climate changes is the investigation of climate cycles recorded in Qinan loess. In this paper, two climatic proxies, magnetic susceptibility and redness, are analyzed for QA-I section to evaluate climate cycles using Maximum entropy spectral analysis and Blackman-Tuckey method. Main conclusions are presented as following: Results exhibit significant peaks at periods of 100 ka, 64 ka, 41 ka, 30 ka and 23 ka, but also 1000 ka, 600 ka and 400 ka. These peaks correspond to the dominant periods of the Earth's orbit parameters, which indicates that the formation of the aeolian sediment in northern China might be primarily driven by variations in orbital parameters. Fluctuations with different cycles respectively dominated in different periods. Major shifts in the dominant cycles occurred at 20.3, 19.0, 17.9, 15.2, 12.5 and 11.3 Myr ago. The transition that happened at 17.9 Myr ago was synchronous with the uplift of the Tibetan Plateau, while others at 15.2, 12.5 and 11.3 Myr ago were in good agreement with the timing of the development of Antarctic ice sheet. Therefore we inferred that these shifts might be related to changes in global ice volume and/or the Tibetan uplift. 3. The strong period of 100 ka is observed between 17.9 and 15.2, and 12.5 and 11.3 Myr ago. Ice sheet-climate models that have been used to explain the cause of the 100 ka period since the middle Pleistocene couldn't be responsible for driving the 100 ka climate cycle in the Miocene in Northern China because of the different boundary of climatic conditions between the Quaternary and Miocene. Further investigation is needed to understand how this cycle became dominant in Qinan loess records during these two time segments.

岱海岩芯沉积的粒度组成与湖区全新世降水变化历史

Sediment cores DH99a and DH99b recovered in the central part of Daihai Lake in north-central China were analysed at 2- to 4-crn intervals for grain-size distribution. Grain-size distributions of the lake sediments are inferred to be a proxy for past changes in East Asian monsoon precipitation, such that greater silt-size percentage and higher median grain size reflect increased monsoonal precipitation rates. The grain-size record of Daihai Lake sediments spanning the last ca 11,000 yr indicates that the monsoonal precipitation in the lake region can be divided into three stages: the Early, Middle and Late Holocene. During the Early Holocene before ca 7900 cal yr BP, the median grain size (Md) and the silt-fraction content were relatively low and constant, suggesting relatively low precipitation over the lake region. The Middle Holocene between ca 7900 and 3100 cal yr BP was marked by intensified and highly variable monsoonal precipitation, as indicated by high and variable Md values and silt contents of the lake sediments. During this period, average precipitation rate gradually increased from ca 7900 to 6900 cal yr BP, displayed intense oscillations between ca 6900 and 4400 cal yr BP, and exhibited a decreasing trend while fluctuating from ca 4400 to 3100 cal yr BP Although generally high during the Middle Holocene, both the Md and the silt content assumed distinctly low values at the short intervals of ca 6500-6400, 6000-5900, 5700-5600, 4400-4200 cal yr BP, implying that monsoonal precipitation might have been significantly reduced during these intervals. During the Late Holocene since ca 3100 cal yr BP, grain-size values suggest that precipitation decreased. However, during the Late Holocene, relatively higher Md values and silt contents occurring between ca 1700 to 1000 cal yr BP may denote an intensification of hydrological cycles in the lake area. Changes in the East Asian monsoonal precipitation were not only directly linked with the changing seasonality of solar insolation resulting from progressive changes in the Earth's orbital parameters, but also may have been closely related to variations in the temperature and size of the Western Pacific Warm Pool, in the intensity of the El Nino-Southern Oscillation, and in the path and strength of the North Equatorial Current in the western Pacific.

Descriptive Simulation: Combining Symbolic and Numerical Methods in the Analysis of Chemical Reaction Mechanisms

The Kineticist's Workbench is a computer program currently under development whose purpose is to help chemists understand, analyze, and simplify complex chemical reaction mechanisms. This paper discusses one module of the program that numerically simulates mechanisms and constructs qualitative descriptions of the simulation results. These descriptions are given in terms that are meaningful to the working chemist (e.g., steady states, stable oscillations, and so on); and the descriptions (as well as the data structures used to construct them) are accessible as input to other programs.

Residual Vibration Reduction in Computer Controlled Machines

Control of machines that exhibit flexibility becomes important when designers attempt to push the state of the art with faster, lighter machines. Three steps are necessary for the control of a flexible planet. First, a good model of the plant must exist. Second, a good controller must be designed. Third, inputs to the controller must be constructed using knowledge of the system dynamic response. There is a great deal of literature pertaining to modeling and control but little dealing with the shaping of system inputs. Chapter 2 examines two input shaping techniques based on frequency domain analysis. The first involves the use of the first deriviate of a gaussian exponential as a driving function template. The second, acasual filtering, involves removal of energy from the driving functions at the resonant frequencies of the system. Chapter 3 presents a linear programming technique for generating vibration-reducing driving functions for systems. Chapter 4 extends the results of the previous chapter by developing a direct solution to the new class of driving functions. A detailed analysis of the new technique is presented from five different perspectives and several extensions are presented. Chapter 5 verifies the theories of the previous two chapters with hardware experiments. Because the new technique resembles common signal filtering, chapter 6 compares the new approach to eleven standard filters. The new technique will be shown to result in less residual vibrations, have better robustness to system parameter uncertainty, and require less computation than other currently used shaping techniques.

Theory of Handwriting

Handwriting production is viewed as a constrained modulation of an underlying oscillatory process. Coupled oscillations in horizontal and vertical directions produce letter forms, and when superimposed on a rightward constant velocity horizontal sweep result in spatially separated letters. Modulation of the vertical oscillation is responsible for control of letter height, either through altering the frequency or altering the acceleration amplitude. Modulation of the horizontal oscillation is responsible for control of corner shape through altering phase or amplitude. The vertical velocity zero crossing in the velocity space diagram is important from the standpoint of control. Changing the horizontal velocity value at this zero crossing controls corner shape, and such changes can be effected through modifying the horizontal oscillation amplitude and phase. Changing the slope at this zero crossing controls writing slant; this slope depends on the horizontal and vertical velocity zero amplitudes and on the relative phase difference. Letter height modulation is also best applied at the vertical velocity zero crossing to preserve an even baseline. The corner shape and slant constraints completely determine the amplitude and phase relations between the two oscillations. Under these constraints interletter separation is not an independent parameter. This theory applies generally to a number of acceleration oscillation patterns such as sinusoidal, rectangular and trapezoidal oscillations. The oscillation theory also provides an explanation for how handwriting might degenerate with speed. An implementation of the theory in the context of the spring muscle model is developed. Here sinusoidal oscillations arise from a purely mechanical sources; orthogonal antagonistic spring pairs generate particular cycloids depending on the initial conditions. Modulating between cycloids can be achieved by changing the spring zero settings at the appropriate times. Frequency can be modulated either by shifting between coactivation and alternating activation of the antagonistic springs or by presuming variable spring constant springs. An acceleration and position measuring apparatus was developed for measurements of human handwriting. Measurements of human writing are consistent with the oscillation theory. It is shown that the minimum energy movement for the spring muscle is bang-coast-bang. For certain parameter values a singular arc solution can be shown to be minimizing. Experimental measurements however indicate that handwriting is not a minimum energy movement.

Periodic Motion Detection and Estimation via Space-Time Sampling

A novel technique to detect and localize periodic movements in video is presented. The distinctive feature of the technique is that it requires neither feature tracking nor object segmentation. Intensity patterns along linear sample paths in space-time are used in estimation of period of object motion in a given sequence of frames. Sample paths are obtained by connecting (in space-time) sample points from regions of high motion magnitude in the first and last frames. Oscillations in intensity values are induced at time instants when an object intersects the sample path. The locations of peaks in intensity are determined by parameters of both cyclic object motion and orientation of the sample path with respect to object motion. The information about peaks is used in a least squares framework to obtain an initial estimate of these parameters. The estimate is further refined using the full intensity profile. The best estimate for the period of cyclic object motion is obtained by looking for consensus among estimates from many sample paths. The proposed technique is evaluated with synthetic videos where ground-truth is known, and with American Sign Language videos where the goal is to detect periodic hand motions.

On the Scalability-Performance Tradeoffs in MPLS and IP Routing

MPLS (Multi-Protocol Label Switching) has recently emerged to facilitate the engineering of network traffic. This can be achieved by directing packet flows over paths that satisfy multiple requirements. MPLS has been regarded as an enhancement to traditional IP routing, which has the following problems: (1) all packets with the same IP destination address have to follow the same path through the network; and (2) paths have often been computed based on static and single link metrics. These problems may cause traffic concentration, and thus degradation in quality of service. In this paper, we investigate by simulations a range of routing solutions and examine the tradeoff between scalability and performance. At one extreme, IP packet routing using dynamic link metrics provides a stateless solution but may lead to routing oscillations. At the other extreme, we consider a recently proposed Profile-based Routing (PBR), which uses knowledge of potential ingress-egress pairs as well as the traffic profile among them. Minimum Interference Routing (MIRA) is another recently proposed MPLS-based scheme, which only exploits knowledge of potential ingress-egress pairs but not their traffic profile. MIRA and the more conventional widest-shortest path (WSP) routing represent alternative MPLS-based approaches on the spectrum of routing solutions. We compare these solutions in terms of utility, bandwidth acceptance ratio as well as their scalability (routing state and computational overhead) and load balancing capability. While the simplest of the per-flow algorithms we consider, the performance of WSP is close to dynamic per-packet routing, without the potential instabilities of dynamic routing.

An Online Distributed Algorithm for Inferring Policy Routing Configurations

We present an online distributed algorithm, the Causation Logging Algorithm (CLA), in which Autonomous Systems (ASes) in the Internet individually report route oscillations/flaps they experience to a central Internet Routing Registry (IRR). The IRR aggregates these reports and may observe what we call causation chains where each node on the chain caused a route flap at the next node along the chain. A chain may also have a causation cycle. The type of an observed causation chain/cycle allows the IRR to infer the underlying policy routing configuration (i.e., the system of economic relationships and constraints on route/path preferences). Our algorithm is based on a formal policy routing model that captures the propagation dynamics of route flaps under arbitrary changes in topology or path preferences. We derive invariant properties of causation chains/cycles for ASes which conform to economic relationships based on the popular Gao-Rexford model. The Gao-Rexford model is known to be safe in the sense that the system always converges to a stable set of paths under static conditions. Our CLA algorithm recovers the type/property of an observed causation chain of an underlying system and determines whether it conforms to the safe economic Gao-Rexford model. Causes for nonconformity can be diagnosed by comparing the properties of the causation chains with those predicted from different variants of the Gao-Rexford model.

Consciousness CLEARS the Mind

A full understanding of consciouness requires that we identify the brain processes from which conscious experiences emerge. What are these processes, and what is their utility in supporting successful adaptive behaviors? Adaptive Resonance Theory (ART) predicted a functional link between processes of Consciousness, Learning, Expectation, Attention, Resonance, and Synchrony (CLEARS), includes the prediction that "all conscious states are resonant states." This connection clarifies how brain dynamics enable a behaving individual to autonomously adapt in real time to a rapidly changing world. The present article reviews theoretical considerations that predicted these functional links, how they work, and some of the rapidly growing body of behavioral and brain data that have provided support for these predictions. The article also summarizes ART models that predict functional roles for identified cells in laminar thalamocortical circuits, including the six layered neocortical circuits and their interactions with specific primary and higher-order specific thalamic nuclei and nonspecific nuclei. These prediction include explanations of how slow perceptual learning can occur more frequently in superficial cortical layers. ART traces these properties to the existence of intracortical feedback loops, and to reset mechanisms whereby thalamocortical mismatches use circuits such as the one from specific thalamic nuclei to nonspecific thalamic nuclei and then to layer 4 of neocortical areas via layers 1-to-5-to-6-to-4.

Spikes, Synchrony, and Attentive Learning by Laminar Thalamocortical Circuits

This article develops the Synchronous Matching Adaptive Resonance Theory (SMART) neural model to explain how the brain may coordinate multiple levels of thalamocortical and corticocortical processing to rapidly learn, and stably remember, important information about a changing world. The model clarifies how bottom-up and top-down processes work together to realize this goal, notably how processes of learning, expectation, attention, resonance, and synchrony are coordinated. The model hereby clarifies, for the first time, how the following levels of brain organization coexist to realize cognitive processing properties that regulate fast learning and stable memory of brain representations: single cell properties, such as spiking dynamics, spike-timing-dependent plasticity (STDP), and acetylcholine modulation; detailed laminar thalamic and cortical circuit designs and their interactions; aggregate cell recordings, such as current-source densities and local field potentials; and single cell and large-scale inter-areal oscillations in the gamma and beta frequency domains. In particular, the model predicts how laminar circuits of multiple cortical areas interact with primary and higher-order specific thalamic nuclei and nonspecific thalamic nuclei to carry out attentive visual learning and information processing. The model simulates how synchronization of neuronal spiking occurs within and across brain regions, and triggers STDP. Matches between bottom-up adaptively filtered input patterns and learned top-down expectations cause gamma oscillations that support attention, resonance, and learning. Mismatches inhibit learning while causing beta oscillations during reset and hypothesis testing operations that are initiated in the deeper cortical layers. The generality of learned recognition codes is controlled by a vigilance process mediated by acetylcholine.

Spiking Dynamics during Perceptual Grouping in the Laminar Circuits of Visual Cortex

Grouping of collinear boundary contours is a fundamental process during visual perception. Illusory contour completion vividly illustrates how stable perceptual boundaries interpolate between pairs of contour inducers, but do not extrapolate from a single inducer. Neural models have simulated how perceptual grouping occurs in laminar visual cortical circuits. These models predicted the existence of grouping cells that obey a bipole property whereby grouping can occur inwardly between pairs or greater numbers of similarly oriented and co-axial inducers, but not outwardly from individual inducers. These models have not, however, incorporated spiking dynamics. Perceptual grouping is a challenge for spiking cells because its properties of collinear facilitation and analog sensitivity to inducer configurations occur despite irregularities in spike timing across all the interacting cells. Other models have demonstrated spiking dynamics in laminar neocortical circuits, but not how perceptual grouping occurs. The current model begins to unify these two modeling streams by implementing a laminar cortical network of spiking cells whose intracellular temporal dynamics interact with recurrent intercellular spiking interactions to quantitatively simulate data from neurophysiological experiments about perceptual grouping, the structure of non-classical visual receptive fields, and gamma oscillations.