USH1945H1  Sustaining chaos by using basin boundary saddles  Google Patents
Sustaining chaos by using basin boundary saddles Download PDFInfo
 Publication number
 USH1945H1 USH1945H1 US09204341 US20434198A USH1945H1 US H1945 H1 USH1945 H1 US H1945H1 US 09204341 US09204341 US 09204341 US 20434198 A US20434198 A US 20434198A US H1945 H1 USH1945 H1 US H1945H1
 Authority
 US
 Grant status
 Grant
 Patent type
 Prior art keywords
 system
 unstable
 state
 region
 parameter
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Abandoned
Links
Images
Classifications

 G—PHYSICS
 G05—CONTROLLING; REGULATING
 G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
 G05B21/00—Systems involving sampling of the variable controlled
 G05B21/02—Systems involving sampling of the variable controlled electric
Abstract
Description
This application has the priority of provisional application Ser. No. 60/067,050 filed Dec. 4, 1997.
In general, the present invention relates to a method and apparatus for sustaining chaos in a system. In particular, the present invention relates to a method and apparatus for preserving chaos in nonchaotic parameter regions by using the natural dynamics of the subject system to redirect flow towards the chaotic region along unstable manifolds of basin boundary saddles, by utilizing small, infrequent parameter perturbations.
In many nonlinear systems, such as chaotic vibrations in structures or lasers, there exist regions in which much of the energy is present within a small range of frequencies. Typically, this occurs when the system is operating at or near resonance. Systems, such as aluminum wings or combustion engines, which operate at resonance, may fail due to repetitive stress caused by driving such a system near resonance.
Although numerous areas in science are now known to exhibit chaos as a natural occurrence, many situations would benefit from the inducement of chaos. In biology, the disappearance of chaos may signal pathological phenomena. In mechanics, chaos could be induced in order to prevent resonance, such as with the aluminum wings or combustion engines noted above. For example, in a system of coupled pendulums, one can excite chaotic motion of several modes to spread the energy over a wide frequency range. In optics, material damage is caused by lasers having a peak intensity at a given temporal frequency, so chaos is desirable since it has broadband spectra. It has also been suggested that chaos occur for normal machine tool cutting, making chaos preservation a desired control for deeper than normal cutting.
A conventional method maintains chaos in a regime where only chaotic transients exist, based on accurate analytical knowledge of the dynamical system, and requiring not only a priori phase space knowledge of escape regions from chaos, but also of preiterates of these regions. Such conventional methods maintain chaos using time series which require monitoring and adjusting the system prior to entering an escape region of the attractor by several iterates, an overly burdensome process.
Moreover, for the flow considered herein, preimages of sets in the escape regions cover much of the chaotic transient region, so monitoring all preimages, as done conventionally, would be tedious.
It is therefore an object of the present invention to provide a method and apparatus for controlling and optimizing the sustain of chaotic transients merely by examining behavior near a saddle.
According to a preferred embodiment of the present invention, such a method for sustaining chaos includes monitoring output parameters of a system, defining a saddle region containing a basin saddle for the system, and redirecting the output parameters toward a chaotic region when the output parameters enter the saddle region, by applying changes to input parameters of the system based on information describing the chaotic region.
According to another aspect of the present invention, a method for maintaining a stable state of a system includes monitoring an adjustable parameter of a system, obtaining information describing an unstable state of the system, and applying parameter perturbations to the system based on the obtained information describing the unstable state. The information describing the unstable state of the system may be the location of the unstable state. Obtaining the information describing the unstable state of the system may include finding a branch connecting a current, stable state of the system to the unstable state, and estimating the location of the unstable state based on the connecting branch. Applying parameter perturbations to the system based on the obtained information describing the unstable state may include implementing a target procedure to direct the system from the stable state toward the unstable state. Implementing a targeting procedure may include varying the adjustable parameter. The method may further include using the location of the unstable state as a reference for control of the system by the application of parameter perturbations. The method may also include determining the parameter perturbations based on the location of the unstable state. Preferably, the stable state is observed, and the unstable state is unobservable.
According to a further aspect of the invention, an apparatus for sustaining chaos includes means for monitoring output parameters of a system, means for defining a saddle region containing a basin saddle for the system, and means for redirecting the output parameters toward a chaotic region when the output parameters enter the saddle region. The latter means includes means for applying changes to input parameters of the system based on received information describing the chaotic region.
According to yet another aspect of the invention, an apparatus for maintaining a stable state of a system includes means for monitoring an adjustable parameter of a system, means for obtaining information describing an unstable state of the system, and means for applying perturbations to the adjustable parameter based on the obtained information describing the unstable state. The information describing the unstable state of the system may be the location of the unstable state. The means for obtaining the information describing the unstable state of the system may include means for finding a branch connecting a current, stable state of the system to the unstable state, and means for estimating the location of the unstable state based on the connecting branch. The means for applying perturbations to the adjustable parameter based on the obtained information describing the unstable state may include means for implementing a target procedure to direct the system from the stable state toward the unstable state. The apparatus may also include means for using the location of the unstable state as a reference for control of the system by the application of parameter perturbations. The apparatus may further include means for determining the parameter perturbations based on the location of the unstable state. The means for implementing a targeting procedure may include means for varying the adjustable parameter. Preferably, the stable state is observed, and the unstable state is unobservable.
These and other features and advantages of the invention will be better understood upon consideration of the following detailed description of the preferred embodiments, as illustrated by the accompanying drawings.
FIG. 1(a) shows stable and unstable manifolds associated with basin boundary periodtwo saddle (S_{a}, S_{b}) at δ=1.88.
FIG. 1(b) shows time series for log (u) at δ=1.88 illustrating a chaotic transient landing on a periodfour attractor (every other iterate is shown).
FIG. 2 shows a probability distribution of escape times as a function of distance from the saddle S_{a }for points inside the chaotic region to the right of S_{a }near the unstable manifold.
FIG. 3 shows a sustained chaotic time series for log (u) and parameter fluctuations used to sustain chaos.
FIG. 4 shows a time series for log (u) and corresponding parameter fluctuations used to sustain chaos based on monitoring region E in FIG. 1a.
FIG. 5 shows a system having stable and unstable states, and an adjustable parameter P.
FIG. 6 shows a connection between the stable and unstable states of the system shown in FIG. 5.
FIG. 7 shows a generalized ndimensional system to which the method of the present invention may be applied.
FIG. 8 is a schematic of an embodiment according to the invention.
This invention provides a framework in which control may be used to spread the energy to a broadband by using small, infrequent perturbations to change the system to one in which chaos prevails. The technique may be applied to any system which was chaotic over a range of values, but is now resonant at a single frequency or over a few frequencies.
The technique of the present invention makes use of quantities which are measurable experimentally and thus can be applied to real data taken directly by time series. In particular, the present invention sustains chaos by using the dynamics of a governing saddle.
The situation addressed by the present invention occurs when there are chaotic transients in the presence of another nonchaotic attractor. The method of the invention maintains chaos by using the natural dynamics of unstable states lying on the basin boundary separating a periodic attractor from chaotic transients, which are referred to herein as basin saddles. Technically, there is only one attractor since chaos is a transient. However, there is still a stable manifold which separates the chaotic transient from the periodic attractor. According to the present invention, once system flow gets in a neighborhood of a basin saddle, small perturbations of an accessible system parameter may be used to redirect the flow towards the chaotic transient region. This may be done by a targeting technique which uses the linearization of the flow about the saddle. A probability distribution of escape times may be used to optimize chaos preservation by targeting regions having long chaotic transients, thus minimizing the number of parameter fluctuations.
In the context of the present invention, “escape region” means a region such that if an iterate falls within it, the subsequent iterates are rapidly drawn to an attractor. The same region is also commonly called a “loss region” in relevant technical fields. Instead of preventing escape to an attractor in advance, the approach of the present invention is to let the system enter a region containing a basin saddle and then redirect the flow back into the chaotic region, by using a targeting technique which uses the natural dynamics around the saddle. This makes the parameter changes very infrequent and, in fact minimal.
A basic advantage of the present invention is that it uses a saddle that is generic to many types of nonlinear systems to sustain chaotic behavior, even though the chaotic behavior itself is unstable. Another advantage is that a targeting algorithm can be made to target points which have very long chaotic transients, thus minimizing interventions by using adjustable parameters. Furthermore, by using a saddle as a reference point, frequency broadening may take place in the presence of parametric drift due to noise.
There are two basic premises that make the inventive method especially simple. First, the system has a saddle having a onedimensional unstable manifold, which crosses a stable manifold of arbitrary dimension, causing a boundary crisis to occur, which results in resonant periodic behavior as the only type of behavior. Secondly, the system has an accessible parameter to control the dynamics near the saddle. Thus, according to the present invention, the saddle may be used to control and optimize the sustain of chaotic transients.
The present invention is applicable to all systems in which it is advantageous to maintain chaos. Generally, consider the following ndimensional system
where y is an ndimensional vector, F is an ndimensional vector field, δ is an adjustable parameter, and t is a scalar representing time. The targeting method for sustaining chaos can be applied to the general system described above, but it is most easily exemplified as applied to a twodimensional system modeling a driven CO_{2 }laser, where
where Ω denotes the angular frequency of the laser output and φ its phase, and u and z denote (scaled) intensity and population inversion, respectively. The control parameter δ represents the amplitude of the drive, and ε_{1}, ε_{2}, and φ are fixed. At δ=δ_{c}=1.84, a crisis occurs between a periodtwo saddle and a chaotic attractor. Just before δ reaches δ_{c}, there exist two attractors, one periodic and one chaotic. The stable manifold of a periodtwo saddle lying on the basin boundary splits the local phase space about the saddle in the following manner: orbits having initial conditions lying to the right of and near one saddle converge to the chaotic attractor. Orbits starting to the left converge to the periodtwo attractor. For δ slightly past δ_{c}, a horseshoe is created from the right; that is, the unstable manifold to the right of the stable manifold crosses the stable manifold (see FIG. 1(a) near S_{a}). Almost all points in the region near the saddle now converge to a periodfour orbit which has perioddoubled off the periodtwo branch. FIG. 1(a) shows the manifold of phase space of the defining periodtwo basin saddle, labeled by S_{a }and S_{b}. The crosses represent the periodfour attractor. The chaotic transients lie near the unstable manifold. FIG. 1(b) shows a chaotic transient settling into a periodfour attractor after the iterate enters the escape region to the left of the upper saddle (every other iterate is shown).
For the inventive method of maintaining and optimizing the length of chaotic transients, the regions near the basin saddles which contain points of long chaotic transients are first identified. The flow is redirected towards these regions once it crosses the basin boundary of the attractor.
In FIG. 2, a distribution of escape times for trajectories starting near S_{a}, close to and to the right of the unstable manifold of S_{a }(that is, the side with chaotic transients) is shown. These points lie on a line segment, L_{a }near S_{a }as shown in FIG. 1. The horizontal axis represents equally spaced intervals on L_{a}. For the points in such an interval, the fraction of points having escape times in time intervals represented on the vertical axis are calculated. In other words, given a time interval on the vertical axis and a space interval on the horizontal axis, the intensity of the grey scale of the distribution indicates the fraction of points in that space interval having escape times in the chosen time interval. This data shows that there is a wide distribution of points near the saddle which escape very quickly to the attractor, as well as points with very long escape times. The fraction of points having a given mean escape time has an exponential distribution. Moreover, since preimages of points in the escape regions essentially cover the chaotic transient region, applying perturbations to prevent the flow from entering the escape regions, which would seem to be a logical approach to sustaining chaos, may actually kick the dynamics into faster escaping regions. Thus, the approach of the present invention, that is, to allow the system to enter the basin saddle region and redirect system flow back into the chaotic region, is more advantageous, although it may appear counterintuitive.
The relations set forth at (1) above may be written in generic form as: φ′=F(t,φ). A poincare map of this flow is obtained by sampling the system with the frequency of the drive, giving:
where δ is the parameter adjusted to maintain chaos. The Jacobian of the map T, calculated at the unstable orbit is given by:
An algorithm used according to the invention perturbs δ once an iterate of T^{2 }enters a neighborhood of the periodtwo basin boundary basin saddle (S_{a}, S_{b}), where S_{a}=T^{2}(S_{a}), and S_{b}=T(S_{a}) T^{2 }is linearized about both S_{a }and S_{b}. Based on the linearization, a formula may be derived for the parameter perturbation that can be used to regenerate chaos. The same procedure must be applied for both S_{a }and S_{b }(which contain escape regions), even when the map considered is T^{2}.
In order to make practical use of basin saddle information, the region of interest is restricted to a local region, D_{loc}, of phasespace near a saddle, for example S_{a }shown in FIG. 1. The region to the left of S_{a }contained in D_{loc }is called the basin escape region. A similar description holds for a local region near S_{b}.
The stable manifold of the periodtwo saddles separates the region with transient chaos from the periodfour attractor. Iterates are sent back to the opposite side of the stable manifold once they enter the basin escape region, preferably close to the unstable manifold of the saddle where the natural dynamics will send the iterates further into the chaotic region. By restricting the action of the dynamic to D_{loc}, the perturbations acting on a linearized system may be examined.
To derive the formula for the parameter perturbation let ξ_{n}=x_{n}−S_{a }and consider the map ξ_{n+1}=P(ξ_{n},δ)=T^{2}(x_{n},δ)−S_{a}(δ), which has fixed point ξ_{F}=0. This map is linearized about its fixed point to obtain:
where A has been expressed in terms of left and right eigenvectors, and
for some {overscore (δ)} close to the operational parameter δ_{0}.
A desired region is targeted near the saddle on the side with chaos, based on the distribution shown in FIG. 2, where it is known that there is a high probability of obtaining long chaotic transients. According to the premise by which (3) was derived, it is required that ξ_{tar}=x_{tar}−S_{a}, where x_{tar}εL_{a}∩D_{loc }is chosen as a target point. The value x_{tar }is chosen such that it is contained in a region where long chaotic transients emerge before the flow escapes again. Equation (3) becomes:
Multiplying through in (4) by f_{s}:
Applying small changes in parameters given by (5) may lead to large changes in preimages of escape regions, but the saddle location varies only slightly.
As an experimental example, the sustained chaotic time series for log (u) is shown in FIG. 3 along with the corresponding parameter perturbations. In this case, only five perturbations to δ were performed to maintain 500 iterates of T^{2}, that is, chaos was sustained for 1000 iterates of T. The noise level was 1% and the parameter perturbations occuring in this case were below 0.6. Certain regions were targeted near the unstable manifold which generate long chaotic transients. This is efficient since four of the five perturbations necessary to maintain chaos were larger than the mean chaotic transient time, which was found to be approximately 27 iterates.
One of the main assumptions made using the method of the invention is that just before entering the basin escape region, the iterates must pass near either saddle S_{a }or S_{b}. However, iterates contained in a region not local to the saddles, but nevertheless in an escape region, can be directed toward D_{loc}, which contains a basin saddle. A region (region E as shown in FIG. 1(a)) is chosen inside the chaotic transients and which lies in the basin of the periodfour attracting orbit.
The idea is to redirect the iterates from such a region by targeting a certain region or point inside the chaotic regime by using, for example, the algorithm introduced above. For example, it has been found that the region near the upper saddle to the right (region F as shown in FIG. 1(a)) is accessible from region E, in only one iterate by small amplitude perturbations of the parameter. The reason is that the image under the flow of a point in E lies inside F. If this were not the case, several intermediate points would have to be targeted before reaching region F. The phasespace sustained chaotic time series is shown in FIG. 4. The changes in the parameter used for targeting from region E to region F are one order of magnitude smaller than the other parameter changes and are not seen clearly in FIG. 4. Several other escape regions such as region E are mapped by the flow into one of the regions near the periodtwo saddle, so it is enough to apply parameter perturbations near the saddle.
Referring to FIG. 5, a further example is shown in which a system has a parameter P which is adjustable, and an observed behavior. In the crosshatched region R, the observed behavior is very complex, containing many frequencies and/or modes. Examples of such systems include combustion engines and fluidstructure interactions in resonance. After the parameter P passes a critical value of P_{c}, the behavior that is observed is periodic, indicated by the solid line PS. It is in this parameter region that the system is operating, and it is desirable to remain in this region. However, the goal is to get information about the unstable state, indicated by the dashed line U. Information about the unstable state may be used to apply parameter perturbations near the saddle. Thus, in addition to the behavior observed in FIG. 5, an unstable object exists, which is not observable because it is unstable. The unstable state U is located as follows.
FIG. 6 shows that the periodic state PS is connected to the unstable state U. By knowing the form of the connecting branch B, which is dictated by theory, one may estimate the location of the unstable state U to use as a reference for control. One then may implement a targeting procedure by varying the parameter P to get from point A of the observed periodic behavior PS to point B of the unobservable unstable behavior U. Once in the neighborhood of B, control may be implemented to stabilize that point, and discern the necessary local information for application to the system to avoid resonance.
The present invention has been described by way of example and in terms of preferred embodiments. However, it is to be understood that the present invention is not strictly limited to the disclosed embodiments. To the contrary, various modifications, as well as similar arrangements, are included within the spirit and scope of the present invention. The scope of the appended claims, therefore, should be accorded the broadest possible interpretation so as to encompass all such modifications and similar arrangements.
For example, the description may be generalized to include systems having any number of dimensions. Such an ndimensional system is described with reference to FIG. 7.
Crises typically occur when a saddle point has its unstable (W^{u}) and stable (W^{s}) manifold intersect, as shown in FIG. 7. The success of the control technique used in the system depends on being able to model the control system as a map. Specifically, the delay between the time when the system is monitored to when control is applied must be very small compared to the overall period. If the technique is attempted on fast systems relative to the control loop delay, such as electronics and optics systems, then the delay must be taken into account. Furthermore, the method of this example is restricted to one in which the parameter change is on for a whole drive period, which forces the system to retain the history of the parameter changes to target in higher dimensions. However, allowing the parameter delay as well as its duration to be adjustable allows one to accomodate more than one unstable dimension.
To see this, consider the steady state saddle of the ndimensional system shown in FIG. 7:
where p is a parameter, and (z,p)=(0,0) is the saddle point having unstable and stable dimensions greater than 1. Linearizing about the saddle yields
where A and B are the usual Jacobian approximations using delay embedding. Let T_{d }and T_{c }denote the delay of actuation and parameter duration, respectively. For a given vector x_{0 }near the saddle, the map, M, is defined by
Parameters are now chosen so that a given point, x_{T}, is targeted to lie in the unstable manifold which generates the chaotic transient. Assuming A=S·Λ·S^{−1}, where S is composed of right eigenvectors, e_{i}, and S^{−1 }of left eigenvectors, f_{i}, i=1 . . . n, Λ diagonal, up to three stable eigendirections, and n−3 unstable directions are allowed. To target x_{T}, M(x_{o},T_{c},T_{d},p)=S·k is solved for (T_{c},T_{d},p,k), where k=(0,0,0,k_{4 }. . . k_{n}). The elements k_{i }are used to target x_{T }on the unstable manifold. Since x_{T }is in the unstable manifold,
Using the fact that
it is easy to show that
FIG. 8 shows a setup for practicing the invention, in which the system under control is carbon dioxide laser 12. Laser diode 14 pumps laser 12, and is itself driven by adjustable drive 24. The magnitude of drive 24 corresponds to δ is the above equations. The output of laser 12 is transduced by photodetector 16, and digitized by member 18. The digitized output of laser 12 is fed to embedder/recorder 20 which may embed the output by any conventional technique, or record the data for later processing. Output from digitizer 18 also goes to comparator 22, which determines the difference between the output of laser 12 and a fixed reference 26, this difference controls the magnitude of offset driver 28, which adjusts the magnitude of input 24 of laser diode 14, thereby adjusting the intensity of the pumping of laser 12. Although the schematic of FIG. 8 shows as the system under control a carbon dioxide laser, one skilled in the art will recognize that within the scope of the invention one could as well employ in its place any other potentially chaotic system, so long as the system has in its phase space a periodic attractor separated by a boundary from a region of chaotic transients. As one skilled in the art will also recognize, there are innumerable such systems, a nonexhaustive list of which includes machine the aircraft wings and machine cutting tools mentioned above, metal fatigue, mechanical vibratory systems, and biological systems, most notably brain function.
Claims (18)
Priority Applications (2)
Application Number  Priority Date  Filing Date  Title 

US6705097 true  19971204  19971204  
US09204341 USH1945H1 (en)  19971204  19981204  Sustaining chaos by using basin boundary saddles 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

US09204341 USH1945H1 (en)  19971204  19981204  Sustaining chaos by using basin boundary saddles 
Publications (1)
Publication Number  Publication Date 

USH1945H1 true USH1945H1 (en)  20010306 
Family
ID=22073399
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

US09204341 Abandoned USH1945H1 (en)  19971204  19981204  Sustaining chaos by using basin boundary saddles 
Country Status (2)
Country  Link 

US (1)  USH1945H1 (en) 
WO (1)  WO1999027765A3 (en) 
Cited By (1)
Publication number  Priority date  Publication date  Assignee  Title 

US6763271B2 (en) *  20010412  20040713  The United States Of America As Represented By The Secretary Of The Navy  Tracking sustained chaos 
Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

US5402334A (en)  19920511  19950328  The United States Of America As Represented By The Secretary Of The Navy  Method and apparatus for pseudoperiodic drive 
US5442510A (en) *  19930623  19950815  The United States Of America As Represented By The Secretary Of The Navy  Control system for tracking nonlinear systems 
US5447520A (en)  19920819  19950905  The Regents Of The University Of California  Real time stabilizing system for pulsating activity 
US5486280A (en)  19941020  19960123  Martin Marietta Energy Systems, Inc.  Process for applying control variables having fractal structures 
US5706192A (en)  19951116  19980106  The United States Of America As Represented By The Secretary Of The Navy  Antiphase switching in arrays of globally coupled oscillators 
US5729388A (en)  19950619  19980317  Massachusetts Institute Of Technology  System employing dissipative pseudorandum dynamics for communications and measurement 
Family Cites Families (1)
Publication number  Priority date  Publication date  Assignee  Title 

US5815413A (en) *  19970508  19980929  Lockheed Martin Energy Research Corporation  Integrated method for chaotic time series analysis 
Patent Citations (7)
Publication number  Priority date  Publication date  Assignee  Title 

US5402334A (en)  19920511  19950328  The United States Of America As Represented By The Secretary Of The Navy  Method and apparatus for pseudoperiodic drive 
US5447520A (en)  19920819  19950905  The Regents Of The University Of California  Real time stabilizing system for pulsating activity 
US5522863A (en)  19920819  19960604  The United States Of America As Represented By The Secretary Of The Navy  Pulsating behavior monitoring and modification system for neural networks 
US5442510A (en) *  19930623  19950815  The United States Of America As Represented By The Secretary Of The Navy  Control system for tracking nonlinear systems 
US5486280A (en)  19941020  19960123  Martin Marietta Energy Systems, Inc.  Process for applying control variables having fractal structures 
US5729388A (en)  19950619  19980317  Massachusetts Institute Of Technology  System employing dissipative pseudorandum dynamics for communications and measurement 
US5706192A (en)  19951116  19980106  The United States Of America As Represented By The Secretary Of The Navy  Antiphase switching in arrays of globally coupled oscillators 
NonPatent Citations (1)
Title 

Ira B. Schwartz and Ioana Triandaf, "Sustaining Chaos by Using Basin Boundary Saddles," Physical Review Letters, vol. 77, No. 23, pp 47404743, Dec. 2, 1996. * 
Cited By (1)
Publication number  Priority date  Publication date  Assignee  Title 

US6763271B2 (en) *  20010412  20040713  The United States Of America As Represented By The Secretary Of The Navy  Tracking sustained chaos 
Also Published As
Publication number  Publication date  Type 

WO1999027765A2 (en)  19990610  application 
WO1999027765A3 (en)  19990715  application 
Similar Documents
Publication  Publication Date  Title 

Zhang et al.  Two improved algorithms for reliability analysis  
Jou et al.  Microstructural evolution in inhomogeneous elastic media  
Zhang et al.  Dynamic response sensitivity of inelastic structures  
US6578006B1 (en)  Project work management method and system  
Mori et al.  Timedependent system reliability analysis by adaptive importance sampling  
Felber  Development of a hybrid bridge evaluation system  
Hastings  Complex interactions between dispersal and dynamics: lessons from coupled logistic equations  
US6224636B1 (en)  Speech recognition using nonparametric speech models  
Parker et al.  Bend theory of river meanders. Part 2. Nonlinear deformation of finiteamplitude bends  
Harman et al.  Improving evolutionary testing by flag removal  
US6675127B2 (en)  Computerized systems and methods for managing project issues and risks  
US20020196292A1 (en)  Graph data visualization and graphics preparation  
Ashpis et al.  The vibrating ribbon problem revisited  
Aha  Generalizing from case studies: A case study  
Zhu et al.  An optimal nonlinear feedback control strategy for randomly excited structural systems  
Mohanty et al.  A modified Ibrahim time domain algorithm for operational modal analysis including harmonic excitation  
Wagner et al.  Symmetric statespace method for a class of nonviscously damped systems  
US5966414A (en)  System and method for processing signals to determine their stochastic properties  
Zavodney et al.  The response of a singledegreeoffreedom system with quadratic and cubic nonlinearities to a principal parametric resonance  
Grigoryan et al.  Dispersionmanaged soliton dynamics  
Lardies et al.  Modal parameter estimation based on the wavelet transform of output data  
Wright  Why a pure primal Newton barrier step may be infeasible  
Karlsen et al.  Convergence of finite difference schemes for viscous and inviscid conservation laws with rough coefficients  
Yu  Simplest normal forms of Hopf and generalized Hopf bifurcations  
US6252369B1 (en)  Motor speed controller and gain setting method of the controller 