Worst-case analysis and linear parameter-varying gain-scheduled control of aerospace systems
NASA Astrophysics Data System (ADS)
Shin, Jong-Yeob
In this thesis, two main subjects are discussed. The first is a worst-case performance analysis, the second is a linear parameter varying (LPV) synthesis using a blending approach. On the first subject, a linear fractional transformation (LFT) model of the linearized X-38 Crew Return Vehicle (CRV) has been developed to facilitate the analysis of its flight control system. The LFT model represents uncertainty in nine aerodynamic stability derivatives at a given flight condition. The X-38 LFT model, combined with a controller at specific flight conditions, is used to determine the aerodynamic coefficients within a predefined set that result in the worst-case performance and worst-case gain/phase margins of the closed-loop system. LPV and mu controllers are synthesized for the X-38 CRV lateral-directional axes over the candidate flight envelope and compared with the baseline gain-scheduled classical control design. Worst-case analysis of the LPV and mu controllers are compared with the baseline gain-scheduled classical control design. Analysis and time simulations show that the LPV controller achieves significant performance and robustness improvements when compared to a linear mu controller and the baseline gain-scheduled controller. On the second subject, a quasi-LPV model of the F-16 longitudinal axes was developed using three methods: Jacobian linearization, state transformation and function substitution. Time simulations of quasi-LPV models show that the quasi-LPV models developed using state transformation and function substitution accurately represent the nonlinear dynamics of the F-16 longitudinal axes. In designing an LPV controller for the F-16 longitudinal axes, the function substitution quasi-LPV models are used since these quasi-LPV models can represent the nonlinear dynamics at non-trim points. Two LPV controllers are synthesized for the F-16 longitudinal axes for two separated flight envelopes: low and high altitude regions. Blending these controllers
Gain Scheduling for the Orion Launch Abort Vehicle Controller
NASA Technical Reports Server (NTRS)
McNamara, Sara J.; Restrepo, Carolina I.; Madsen, Jennifer M.; Medina, Edgar A.; Proud, Ryan W.; Whitley, Ryan J.
2011-01-01
One of NASAs challenges for the Orion vehicle is the control system design for the Launch Abort Vehicle (LAV), which is required to abort safely at any time during the atmospheric ascent portion of ight. The focus of this paper is the gain design and scheduling process for a controller that covers the wide range of vehicle configurations and flight conditions experienced during the full envelope of potential abort trajectories from the pad to exo-atmospheric flight. Several factors are taken into account in the automation process for tuning the gains including the abort effectors, the environmental changes and the autopilot modes. Gain scheduling is accomplished using a linear quadratic regulator (LQR) approach for the decoupled, simplified linear model throughout the operational envelope in time, altitude and Mach number. The derived gains are then implemented into the full linear model for controller requirement validation. Finally, the gains are tested and evaluated in a non-linear simulation using the vehicles ight software to ensure performance requirements are met. An overview of the LAV controller design and a description of the linear plant models are presented. Examples of the most significant challenges with the automation of the gain tuning process are then discussed. In conclusion, the paper will consider the lessons learned through out the process, especially in regards to automation, and examine the usefulness of the gain scheduling tool and process developed as applicable to non-Orion vehicles.
Prompt critical control of the ACRR using a linear quadratic regulator design
NASA Astrophysics Data System (ADS)
Gilkey, Jeffrey C.
1993-01-01
This paper describes the application of ``Modern Control'' design techniques to the problem of nuclear reactor control. The control algorithm consists of generating a nominal trajectory within the control authority of the reactor rod drives, and then following this trajectory with a gain scheduled linear quadratic regulator (LQR). A controller based on this algorithm has generated power pulses up to 100 mW on Sandia's Annular Core Research Reactor (ACRR). Prompt critical control at 1.02 net reactivity and at start-up rates over 350 decades per minute (DPM) has also been demonstrated using this controller.
Communications circuit including a linear quadratic estimator
Ferguson, Dennis D.
2015-07-07
A circuit includes a linear quadratic estimator (LQE) configured to receive a plurality of measurements a signal. The LQE is configured to weight the measurements based on their respective uncertainties to produce weighted averages. The circuit further includes a controller coupled to the LQE and configured to selectively adjust at least one data link parameter associated with a communication channel in response to receiving the weighted averages.
Extended Decentralized Linear-Quadratic-Gaussian Control
NASA Technical Reports Server (NTRS)
Carpenter, J. Russell
2000-01-01
A straightforward extension of a solution to the decentralized linear-Quadratic-Gaussian problem is proposed that allows its use for commonly encountered classes of problems that are currently solved with the extended Kalman filter. This extension allows the system to be partitioned in such a way as to exclude the nonlinearities from the essential algebraic relationships that allow the estimation and control to be optimally decentralized.
LPV Controller Interpolation for Improved Gain-Scheduling Control Performance
NASA Technical Reports Server (NTRS)
Wu, Fen; Kim, SungWan
2002-01-01
In this paper, a new gain-scheduling control design approach is proposed by combining LPV (linear parameter-varying) control theory with interpolation techniques. The improvement of gain-scheduled controllers can be achieved from local synthesis of Lyapunov functions and continuous construction of a global Lyapunov function by interpolation. It has been shown that this combined LPV control design scheme is capable of improving closed-loop performance derived from local performance improvement. The gain of the LPV controller will also change continuously across parameter space. The advantages of the newly proposed LPV control is demonstrated through a detailed AMB controller design example.
A Note on the Linearly and Quadratically Weighted Kappa Coefficients.
Li, Pingke
2016-09-01
The linearly and quadratically weighted kappa coefficients are popular statistics in measuring inter-rater agreement on an ordinal scale. It has been recently demonstrated that the linearly weighted kappa is a weighted average of the kappa coefficients of the embedded 2 by 2 agreement matrices, while the quadratically weighted kappa is insensitive to the agreement matrices that are row or column reflection symmetric. A rank-one matrix decomposition approach to the weighting schemes is presented in this note such that these phenomena can be demonstrated in a concise manner. PMID:27246436
Design of Linear Quadratic Regulators and Kalman Filters
NASA Technical Reports Server (NTRS)
Lehtinen, B.; Geyser, L.
1986-01-01
AESOP solves problems associated with design of controls and state estimators for linear time-invariant systems. Systems considered are modeled in state-variable form by set of linear differential and algebraic equations with constant coefficients. Two key problems solved by AESOP are linear quadratic regulator (LQR) design problem and steady-state Kalman filter design problem. AESOP is interactive. User solves design problems and analyzes solutions in single interactive session. Both numerical and graphical information available to user during the session.
Finding Optimal Gains In Linear-Quadratic Control Problems
NASA Technical Reports Server (NTRS)
Milman, Mark H.; Scheid, Robert E., Jr.
1990-01-01
Analytical method based on Volterra factorization leads to new approximations for optimal control gains in finite-time linear-quadratic control problem of system having infinite number of dimensions. Circumvents need to analyze and solve Riccati equations and provides more transparent connection between dynamics of system and optimal gain.
Analysis of integral controls in linear quadratic regulator design
NASA Technical Reports Server (NTRS)
Slater, G. L.
1979-01-01
The application of linear optimal control to the design of systems with integral control action on specified outputs is considered. Using integral terms in a quadratic performance index, an asymptotic analysis is used to determine the effect of variable quadratic weights on the eigenvalues and eigenvectors of the closed loop system. It is shown that for small integral terms the placement of integrator poles and gain calculation can be effectively decoupled from placement of the primary system eigenvalues. This technique is applied to the design of integral controls for a STOL aircraft outer loop guidance system.
Reaction Wheel Control Design Using Linear Quadratic Controller
NASA Astrophysics Data System (ADS)
Nubli Muhamad, Nur; Susanto, Erwin; Syihabuddin, Budi; Prasetya Dwi Wibawa, Ig.
2016-01-01
This paper studies the design of active attitude control system of a nanosatellite in a single axis. In this paper, we consider dc motor based reaction wheel as an actuator, because of its pointing accuracy. However, the power consumption of the dc motor is often relatively large and needed to be optimized. Linear quadratic controller is supposed to have an ability to minimize power consumption and able to enhance the system performance. To show the advantage of this method, simulation result of attitude response, state trajectory, and trajectory of DC motor voltage are presented.
Linear-Quadratic-Gaussian Regulator Developed for a Magnetic Bearing
NASA Technical Reports Server (NTRS)
Choi, Benjamin B.
2002-01-01
Linear-Quadratic-Gaussian (LQG) control is a modern state-space technique for designing optimal dynamic regulators. It enables us to trade off regulation performance and control effort, and to take into account process and measurement noise. The Structural Mechanics and Dynamics Branch at the NASA Glenn Research Center has developed an LQG control for a fault-tolerant magnetic bearing suspension rig to optimize system performance and to reduce the sensor and processing noise. The LQG regulator consists of an optimal state-feedback gain and a Kalman state estimator. The first design step is to seek a state-feedback law that minimizes the cost function of regulation performance, which is measured by a quadratic performance criterion with user-specified weighting matrices, and to define the tradeoff between regulation performance and control effort. The next design step is to derive a state estimator using a Kalman filter because the optimal state feedback cannot be implemented without full state measurement. Since the Kalman filter is an optimal estimator when dealing with Gaussian white noise, it minimizes the asymptotic covariance of the estimation error.
Robust linear quadratic designs with respect to parameter uncertainty
NASA Technical Reports Server (NTRS)
Douglas, Joel; Athans, Michael
1992-01-01
The authors derive a linear quadratic regulator (LQR) which is robust to parametric uncertainty by using the overbounding method of I. R. Petersen and C. V. Hollot (1986). The resulting controller is determined from the solution of a single modified Riccati equation. It is shown that, when applied to a structural system, the controller gains add robustness by minimizing the potential energy of uncertain stiffness elements, and minimizing the rate of dissipation of energy through uncertain damping elements. A worst-case disturbance in the direction of the uncertainty is also considered. It is proved that performance robustness has been increased with the robust LQR when compared to a mismatched LQR design where the controller is designed on the nominal system, but applied to the actual uncertain system.
Sensitivity Analysis of Parameters in Linear-Quadratic Radiobiologic Modeling
Fowler, Jack F.
2009-04-01
Purpose: Radiobiologic modeling is increasingly used to estimate the effects of altered treatment plans, especially for dose escalation. The present article shows how much the linear-quadratic (LQ) (calculated biologically equivalent dose [BED] varies when individual parameters of the LQ formula are varied by {+-}20% and by 1%. Methods: Equivalent total doses (EQD2 = normalized total doses (NTD) in 2-Gy fractions for tumor control, acute mucosal reactions, and late complications were calculated using the linear- quadratic formula with overall time: BED = nd (1 + d/ [{alpha}/{beta}]) - log{sub e}2 (T - Tk) / {alpha}Tp, where BED is BED = total dose x relative effectiveness (RE = nd (1 + d/ [{alpha}/{beta}]). Each of the five biologic parameters in turn was altered by {+-}10%, and the altered EQD2s tabulated; the difference was finally divided by 20. EQD2 or NTD is obtained by dividing BED by the RE for 2-Gy fractions, using the appropriate {alpha}/{beta} ratio. Results: Variations in tumor and acute mucosal EQD ranged from 0.1% to 0.45% per 1% change in each parameter for conventional schedules, the largest variation being caused by overall time. Variations in 'late' EQD were 0.4% to 0.6% per 1% change in the only biologic parameter, the {alpha}/{beta} ratio. For stereotactic body radiotherapy schedules, variations were larger, up to 0.6 to 0.9 for tumor and 1.6% to 1.9% for late, per 1% change in parameter. Conclusions: Robustness occurs similar to that of equivalent uniform dose (EUD), for the same reasons. Total dose, dose per fraction, and dose-rate cause their major effects, as well known.
AESOP- INTERACTIVE DESIGN OF LINEAR QUADRATIC REGULATORS AND KALMAN FILTERS
NASA Technical Reports Server (NTRS)
Lehtinen, B.
1994-01-01
AESOP was developed to solve a number of problems associated with the design of controls and state estimators for linear time-invariant systems. The systems considered are modeled in state-variable form by a set of linear differential and algebraic equations with constant coefficients. Two key problems solved by AESOP are the linear quadratic regulator (LQR) design problem and the steady-state Kalman filter design problem. AESOP is designed to be used in an interactive manner. The user can solve design problems and analyze the solutions in a single interactive session. Both numerical and graphical information are available to the user during the session. The AESOP program is structured around a list of predefined functions. Each function performs a single computation associated with control, estimation, or system response determination. AESOP contains over sixty functions and permits the easy inclusion of user defined functions. The user accesses these functions either by inputting a list of desired functions in the order they are to be performed, or by specifying a single function to be performed. The latter case is used when the choice of function and function order depends on the results of previous functions. The available AESOP functions are divided into several general areas including: 1) program control, 2) matrix input and revision, 3) matrix formation, 4) open-loop system analysis, 5) frequency response, 6) transient response, 7) transient function zeros, 8) LQR and Kalman filter design, 9) eigenvalues and eigenvectors, 10) covariances, and 11) user-defined functions. The most important functions are those that design linear quadratic regulators and Kalman filters. The user interacts with AESOP when using these functions by inputting design weighting parameters and by viewing displays of designed system response. Support functions obtain system transient and frequency responses, transfer functions, and covariance matrices. AESOP can also provide the user
A Novel Approach to Noise-Filtering Based on a Gain-Scheduling Neural Network Architecture
NASA Technical Reports Server (NTRS)
Troudet, T.; Merrill, W.
1994-01-01
A gain-scheduling neural network architecture is proposed to enhance the noise-filtering efficiency of feedforward neural networks, in terms of both nominal performance and robustness. The synergistic benefits of the proposed architecture are demonstrated and discussed in the context of the noise-filtering of signals that are typically encountered in aerospace control systems. The synthesis of such a gain-scheduled neurofiltering provides the robustness of linear filtering, while preserving the nominal performance advantage of conventional nonlinear neurofiltering. Quantitative performance and robustness evaluations are provided for the signal processing of pitch rate responses to typical pilot command inputs for a modern fighter aircraft model.
Linear versus quadratic portfolio optimization model with transaction cost
NASA Astrophysics Data System (ADS)
Razak, Norhidayah Bt Ab; Kamil, Karmila Hanim; Elias, Siti Masitah
2014-06-01
Optimization model is introduced to become one of the decision making tools in investment. Hence, it is always a big challenge for investors to select the best model that could fulfill their goal in investment with respect to risk and return. In this paper we aims to discuss and compare the portfolio allocation and performance generated by quadratic and linear portfolio optimization models namely of Markowitz and Maximin model respectively. The application of these models has been proven to be significant and popular among others. However transaction cost has been debated as one of the important aspects that should be considered for portfolio reallocation as portfolio return could be significantly reduced when transaction cost is taken into consideration. Therefore, recognizing the importance to consider transaction cost value when calculating portfolio' return, we formulate this paper by using data from Shariah compliant securities listed in Bursa Malaysia. It is expected that, results from this paper will effectively justify the advantage of one model to another and shed some lights in quest to find the best decision making tools in investment for individual investors.
Linear quadratic optimal controller for cable-driven parallel robots
NASA Astrophysics Data System (ADS)
Abdolshah, Saeed; Shojaei Barjuei, Erfan
2015-12-01
In recent years, various cable-driven parallel robots have been investigated for their advantages, such as low structural weight, high acceleration, and large work-space, over serial and conventional parallel systems. However, the use of cables lowers the stiffness of these robots, which in turn may decrease motion accuracy. A linear quadratic (LQ) optimal controller can provide all the states of a system for the feedback, such as position and velocity. Thus, the application of such an optimal controller in cable-driven parallel robots can result in more efficient and accurate motion compared to the performance of classical controllers such as the proportional- integral-derivative controller. This paper presents an approach to apply the LQ optimal controller on cable-driven parallel robots. To employ the optimal control theory, the static and dynamic modeling of a 3-DOF planar cable-driven parallel robot (Feriba-3) is developed. The synthesis of the LQ optimal control is described, and the significant experimental results are presented and discussed.
Linear state feedback, quadratic weights, and closed loop eigenstructures. M.S. Thesis
NASA Technical Reports Server (NTRS)
Thompson, P. M.
1979-01-01
Results are given on the relationships between closed loop eigenstructures, state feedback gain matrices of the linear state feedback problem, and quadratic weights of the linear quadratic regulator. Equations are derived for the angles of general multivariable root loci and linear quadratic optimal root loci, including angles of departure and approach. The generalized eigenvalue problem is used for the first time to compute angles of approach. Equations are also derived to find the sensitivity of closed loop eigenvalues and the directional derivatives of closed loop eigenvectors (with respect to a scalar multiplying the feedback gain matrix or the quadratic control weight). An equivalence class of quadratic weights that produce the same asymptotic eigenstructure is defined, sufficient conditions to be in it are given, a canonical element is defined, and an algorithm to find it is given. The behavior of the optimal root locus in the nonasymptotic region is shown to be different for quadratic weights with the same asymptotic properties.
A decentralized linear quadratic control design method for flexible structures
NASA Technical Reports Server (NTRS)
Su, Tzu-Jeng; Craig, Roy R., Jr.
1990-01-01
A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass
Singular linear-quadratic control problem for systems with linear delay
Sesekin, A. N.
2013-12-18
A singular linear-quadratic optimization problem on the trajectories of non-autonomous linear differential equations with linear delay is considered. The peculiarity of this problem is the fact that this problem has no solution in the class of integrable controls. To ensure the existence of solutions is required to expand the class of controls including controls with impulse components. Dynamical systems with linear delay are used to describe the motion of pantograph from the current collector with electric traction, biology, etc. It should be noted that for practical problems fact singularity criterion of quality is quite commonly occurring, and therefore the study of these problems is surely important. For the problem under discussion optimal programming control contained impulse components at the initial and final moments of time is constructed under certain assumptions on the functional and the right side of the control system.
Self-consistent linearization of non-linear BEM formulations with quadratic convergence
NASA Astrophysics Data System (ADS)
Fernandes, G. R.; de Souza Neto, E. A.
2013-11-01
In this work, a general technique to obtain the self-consistent linearization of non-linear formulations of the boundary element method (BEM) is presented. In the incremental-iterative procedure required to solve the non-linear problem the convergence is quadratic, being the solution obtained from the consistent tangent operator. This technique is applied to non-linear BEM formulations for plates where two independent problems are discussed: the plate bending and the stretching problem. For both problems an equilibrium equation is written in terms of strains and internal forces and then the consistent tangent operator is derived by applying the Newton-Raphson’s scheme. The Von Mises criterion is adopted to govern the elasto-plastic material behaviour checked at points along the plate thickness, although the presented formulations can be used with any non-linear model. Numerical examples are presented showing the accuracy of the results as well as the high convergence rate of the iterative procedure.
Closed-loop structural stability for linear-quadratic optimal systems
NASA Technical Reports Server (NTRS)
Wong, P. K.; Athans, M.
1975-01-01
This paper contains an explicit parameterization of a subclass of linear constant gain feedback maps that never destabilize an originally open-loop stable system. These results can then be used to obtain several new structural stability results for multi-input linear-quadratic feedback optimal designs.
Propagator for the time-dependent charged oscillator via linear and quadratic invariants
Abdalla, M. Sebawe Choi, Jeong-Ryeol
2007-12-15
The problem of a charged particle in the presence of a variable magnetic field is considered. Using the linear and the quadratic invariants as a tool, the wave functions in Fock state as well as in coherent state are obtained. The corresponding propagators which propagate the wave functions in the space-time are derived. Using numerical computations we have managed to draw some plots for the real, imaginary, and absolute values of the propagators. This has been used to analyze the properties of the propagators associated with both of the linear and the quadratic invariants. It has been shown that there is no essential difference between the behavior of the absolute value of the propagators in both of the linear and the quadratic invariants.
Nonadiabatic effects in ultracold molecules via anomalous linear and quadratic Zeeman shifts.
McGuyer, B H; Osborn, C B; McDonald, M; Reinaudi, G; Skomorowski, W; Moszynski, R; Zelevinsky, T
2013-12-13
Anomalously large linear and quadratic Zeeman shifts are measured for weakly bound ultracold 88Sr2 molecules near the intercombination-line asymptote. Nonadiabatic Coriolis coupling and the nature of long-range molecular potentials explain how this effect arises and scales roughly cubically with the size of the molecule. The linear shifts yield nonadiabatic mixing angles of the molecular states. The quadratic shifts are sensitive to nearby opposite f-parity states and exhibit fourth-order corrections, providing a stringent test of a state-of-the-art ab initio model. PMID:24483652
Nonadiabatic Effects in Ultracold Molecules via Anomalous Linear and Quadratic Zeeman Shifts
NASA Astrophysics Data System (ADS)
McGuyer, B. H.; Osborn, C. B.; McDonald, M.; Reinaudi, G.; Skomorowski, W.; Moszynski, R.; Zelevinsky, T.
2013-12-01
Anomalously large linear and quadratic Zeeman shifts are measured for weakly bound ultracold $^{88}$Sr$_2$ molecules near the intercombination-line asymptote. Nonadiabatic Coriolis coupling and the nature of long-range molecular potentials explain how this effect arises and scales roughly cubically with the size of the molecule. The linear shifts yield nonadiabatic mixing angles of the molecular states. The quadratic shifts are sensitive to nearby opposite $f$-parity states and exhibit fourth-order corrections, providing a stringent test of a state-of-the-art \\textit{ab initio} model.
NASA Technical Reports Server (NTRS)
Fleming, P.
1985-01-01
A design technique is proposed for linear regulators in which a feedback controller of fixed structure is chosen to minimize an integral quadratic objective function subject to the satisfaction of integral quadratic constraint functions. Application of a non-linear programming algorithm to this mathematically tractable formulation results in an efficient and useful computer-aided design tool. Particular attention is paid to computational efficiency and various recommendations are made. Two design examples illustrate the flexibility of the approach and highlight the special insight afforded to the designer.
A path-following interior-point algorithm for linear and quadratic problems
Wright, S.J.
1993-12-01
We describe an algorithm for the monotone linear complementarity problem that converges for many positive, not necessarily feasible, starting point and exhibits polynomial complexity if some additional assumptions are made on the starting point. If the problem has a strictly complementary solution, the method converges subquadratically. We show that the algorithm and its convergence extend readily to the mixed monotone linear complementarity problem and, hence, to all the usual formulations of the linear programming and convex quadratic programming problems.
Linear quadratic servo control of a reusable rocket engine
NASA Technical Reports Server (NTRS)
Musgrave, Jeffrey L.
1991-01-01
A design method for a servo compensator is developed in the frequency domain using singular values. The method is applied to a reusable rocket engine. An intelligent control system for reusable rocket engines was proposed which includes a diagnostic system, a control system, and an intelligent coordinator which determines engine control strategies based on the identified failure modes. The method provides a means of generating various linear multivariable controllers capable of meeting performance and robustness specifications and accommodating failure modes identified by the diagnostic system. Command following with set point control is necessary for engine operation. A Kalman filter reconstructs the state while loop transfer recovery recovers the required degree of robustness while maintaining satisfactory rejection of sensor noise from the command error. The approach is applied to the design of a controller for a rocket engine satisfying performance constraints in the frequency domain. Simulation results demonstrate the performance of the linear design on a nonlinear engine model over all power levels during mainstage operation.
Graphical Description of Johnson-Neyman Outcomes for Linear and Quadratic Regression Surfaces.
ERIC Educational Resources Information Center
Schafer, William D.; Wang, Yuh-Yin
A modification of the usual graphical representation of heterogeneous regressions is described that can aid in interpreting significant regions for linear or quadratic surfaces. The standard Johnson-Neyman graph is a bivariate plot with the criterion variable on the ordinate and the predictor variable on the abscissa. Regression surfaces are drawn…
NASA Astrophysics Data System (ADS)
Benhalouche, Fatima Zohra; Karoui, Moussa Sofiane; Deville, Yannick; Ouamri, Abdelaziz
2015-10-01
In this paper, a new Spectral-Unmixing-based approach, using Nonnegative Matrix Factorization (NMF), is proposed to locally multi-sharpen hyperspectral data by integrating a Digital Surface Model (DSM) obtained from LIDAR data. In this new approach, the nature of the local mixing model is detected by using the local variance of the object elevations. The hyper/multispectral images are explored using small zones. In each zone, the variance of the object elevations is calculated from the DSM data in this zone. This variance is compared to a threshold value and the adequate linear/linearquadratic spectral unmixing technique is used in the considered zone to independently unmix hyperspectral and multispectral data, using an adequate linear/linear-quadratic NMF-based approach. The obtained spectral and spatial information thus respectively extracted from the hyper/multispectral images are then recombined in the considered zone, according to the selected mixing model. Experiments based on synthetic hyper/multispectral data are carried out to evaluate the performance of the proposed multi-sharpening approach and literature linear/linear-quadratic approaches used on the whole hyper/multispectral data. In these experiments, real DSM data are used to generate synthetic data containing linear and linear-quadratic mixed pixel zones. The DSM data are also used for locally detecting the nature of the mixing model in the proposed approach. Globally, the proposed approach yields good spatial and spectral fidelities for the multi-sharpened data and significantly outperforms the used literature methods.
Gain-scheduled controller synthesis for a nonlinear PDE
NASA Astrophysics Data System (ADS)
Mahdi Hashemi, Seyed; Werner, Herbert
2012-01-01
Linear parameter-varying (LPV) modelling and control of a nonlinear partial differential equation (PDE) is considered in this article. The one-dimensional viscous Burgers' equation is discretised using a finite difference scheme; the boundary conditions are taken as control inputs and the velocities at two grid points are assumed to be measurable. A nonlinear high-order state space model is generated and proper orthogonal decomposition is used for model order reduction. After assessing the accuracy of the reduced model, a low-order functional observer is designed to estimate the reduced states which are linear combinations of the velocities at all grid points. A discrete-time quasi-LPV model that is affine in scheduling parameters is derived based on the reduced model. A polytopic LPV controller is synthesised based on a generalised plant containing the LPV model and the functional observer. More generally, the proposed method can be used to design an LPV controller for a quasi-LPV system with non-measurable scheduling parameters. Simulation results demonstrate the high tracking performance and disturbance and measurement noise rejection capabilities of the designed LPV controller compared with a linear quadratic Gaussian (LQG) controller based on a linearised model.
Gain-Scheduled Complementary Filter Design for a MEMS Based Attitude and Heading Reference System
Yoo, Tae Suk; Hong, Sung Kyung; Yoon, Hyok Min; Park, Sungsu
2011-01-01
This paper describes a robust and simple algorithm for an attitude and heading reference system (AHRS) based on low-cost MEMS inertial and magnetic sensors. The proposed approach relies on a gain-scheduled complementary filter, augmented by an acceleration-based switching architecture to yield robust performance, even when the vehicle is subject to strong accelerations. Experimental results are provided for a road captive test during which the vehicle dynamics are in high-acceleration mode and the performance of the proposed filter is evaluated against the output from a conventional linear complementary filter. PMID:22163824
Clifford group, stabilizer states, and linear and quadratic operations over GF(2)
Dehaene, Jeroen; Moor, Bart de
2003-10-01
We describe stabilizer states and Clifford group operations using linear operations and quadratic forms over binary vector spaces. We show how the n-qubit Clifford group is isomorphic to a group with an operation that is defined in terms of a (2n+1)x(2n+1) binary matrix product and binary quadratic forms. As an application we give two schemes to efficiently decompose Clifford group operations into one- and two-qubit operations. We also show how the coefficients of stabilizer states and Clifford group operations in a standard basis expansion can be described by binary quadratic forms. Our results are useful for quantum error correction, entanglement distillation, and possibly quantum computing.
Model Reduction by Balanced Truncation of Linear Systems with a Quadratic Output
NASA Astrophysics Data System (ADS)
Van Beeumen, Roel; Meerbergen, Karl
2010-09-01
Balanced truncation is a widely used and appreciated projection-based model reduction technique for linear systems. This technique has the following two important properties: approximations by balanced truncation preserve the stability and the H∞-norm (the maximum of the frequency response) of the error system is bounded above by twice the sum of the neglected singular values. This paper tries to extend the framework of linear balanced truncation to systems with a quadratic output. For such systems, the controllability Gramian remains the same. The observability Gramian is computed from a linear system with multiple outputs that is derived from the quadratic output of the original system. We give a numerical example for a large-scale system arising from structural analysis.
Linear state feedback, quadratic weights, and closed loop eigenstructures. M.S. Thesis. Final Report
NASA Technical Reports Server (NTRS)
Thompson, P. M.
1980-01-01
Equations are derived for the angles of general multivariable root loci and linear quadratic optimal root loci, including angles of departure and approach. The generalized eigenvalue problem is used to compute angles of approach. Equations are also derived to find the sensitivity of closed loop eigenvalue and the directional derivatives of closed loop eigenvectors. An equivalence class of quadratic weights that produce the same asymptotic eigenstructure is defined, a canonical element is defined, and an algorithm to find it is given. The behavior of the optimal root locus in the nonasymptotic region is shown to be different for quadratic weights with the same asymptotic properties. An algorithm is presented that can be used to select a feedback gain matrix for the linear state feedback problem which produces a specified asymptotic eigenstructure. Another algorithm is given to compute the asymptotic eigenstructure properties inherent in a given set of quadratic weights. Finally, it is shown that optimal root loci for nongeneric problems can be approximated by generic ones in the nonasymptotic region.
NASA Technical Reports Server (NTRS)
Gibson, J. S.; Rosen, I. G.
1986-01-01
An abstract approximation framework is developed for the finite and infinite time horizon discrete-time linear-quadratic regulator problem for systems whose state dynamics are described by a linear semigroup of operators on an infinite dimensional Hilbert space. The schemes included the framework yield finite dimensional approximations to the linear state feedback gains which determine the optimal control law. Convergence arguments are given. Examples involving hereditary and parabolic systems and the vibration of a flexible beam are considered. Spline-based finite element schemes for these classes of problems, together with numerical results, are presented and discussed.
NASA Technical Reports Server (NTRS)
Banks, H. T.; Silcox, R. J.; Keeling, S. L.; Wang, C.
1989-01-01
A unified treatment of the linear quadratic tracking (LQT) problem, in which a control system's dynamics are modeled by a linear evolution equation with a nonhomogeneous component that is linearly dependent on the control function u, is presented; the treatment proceeds from the theoretical formulation to a numerical approximation framework. Attention is given to two categories of LQT problems in an infinite time interval: the finite energy and the finite average energy. The behavior of the optimal solution for finite time-interval problems as the length of the interval tends to infinity is discussed. Also presented are the formulations and properties of LQT problems in a finite time interval.
Application of gain scheduling to the control of batch bioreactors
NASA Technical Reports Server (NTRS)
Cardello, Ralph; San, Ka-Yiu
1987-01-01
The implementation of control algorithms to batch bioreactors is often complicated by the inherent variations in process dynamics during the course of fermentation. Such a wide operating range may render the performance of fixed gain PID controllers unsatisfactory. In this work, a detailed study on the control of batch fermentation is performed. Furthermore, a simple batch controller design is proposed which incorporates the concept of gain-scheduling, a subclass of adaptive control, with oxygen uptake rate as an auxiliary variable. The control of oxygen tension in the biorector is used as a vehicle to convey the proposed idea, analysis and results. Simulation experiments indicate significant improvement in controller performance can be achieved by the proposed approach even in the presence of measurement noise.
NASA Technical Reports Server (NTRS)
Ito, K.; Teglas, R.
1984-01-01
The numerical scheme based on the Legendre-tau approximation is proposed to approximate the feedback solution to the linear quadratic optimal control problem for hereditary differential systems. The convergence property is established using Trotter ideas. The method yields very good approximations at low orders and provides an approximation technique for computing closed-loop eigenvalues of the feedback system. A comparison with existing methods (based on averaging and spline approximations) is made.
NASA Technical Reports Server (NTRS)
Ito, Kazufumi; Teglas, Russell
1987-01-01
The numerical scheme based on the Legendre-tau approximation is proposed to approximate the feedback solution to the linear quadratic optimal control problem for hereditary differential systems. The convergence property is established using Trotter ideas. The method yields very good approximations at low orders and provides an approximation technique for computing closed-loop eigenvalues of the feedback system. A comparison with existing methods (based on averaging and spline approximations) is made.
A linear-quadratic-Gaussian control problem with innovations-feedthrough solution
NASA Technical Reports Server (NTRS)
Platzman, L. K.; Johnson, T. L.
1976-01-01
The structure of the separation-theorem solution to the standard linear-quadratic-Gaussian (LQG) control problem does not involve direct output feedback as a consequence of the form of the performance index. It is shown that the performance index may be generalized in a natural fashion so that the optimal control law involves output feedback or, equivalently, innovations feedthrough (IF). Applications where this formulation may be advantageous are indicated through an examination of properties of the IF control law.
ORACLS - A linear-quadratic-Gaussian computer-aided design package
NASA Technical Reports Server (NTRS)
Armstrong, E. S.
1982-01-01
ORACLS, an acronym denoting Optimal Regular Algorithms for the Control of Linear Systems, is a collection of FORTRAN coded subroutines dedicated to the formulation and solution of the Linear-Quadratic-Gaussian (LQG) design problem modeled in both continuous and discrete form. The ORACLS system is under continuous development at the NASA Langley Research Center, Hampton, Virginia, and is widely used by universities and industry within the U.S.A. The current (operational) ORACLS version as well as new software under development is described.
A linear quadratic tracker for Control Moment Gyro based attitude control of the Space Station
NASA Technical Reports Server (NTRS)
Kaidy, J. T.
1986-01-01
The paper discusses a design for an attitude control system for the Space Station which produces fast response, with minimal overshoot and cross-coupling with the use of Control Moment Gyros (CMG). The rigid body equations of motion are linearized and discretized and a Linear Quadratic Regulator (LQR) design and analysis study is performed. The resulting design is then modified such that integral and differential terms are added to the state equations to enhance response characteristics. Methods for reduction of computation time through channelization are discussed as well as the reduction of initial torque requirements.
Anisotropic models with two fluids in linear and quadratic forms of f( T) gravitational theories
NASA Astrophysics Data System (ADS)
Nashed, Gamal G. L.
2015-06-01
Recent astronomical observations show that the universe may be anisotropic on large scales. The Union2 SnIa data hint that the universe has a preferred direction. If such a cosmological privileged axis indeed exists, one has to consider an anisotropic expanding universe, instead of the isotropic cosmological model. In this study, we apply the field equations of quadratic form of the modified teleparallel gravitational theories, f( T)= T+ ɛT 2, to anisotropic model. We assume two fluid components, the matter components have two equation of states (EoS). We study different equation of states for the linear case and show that there is no recombination era between the two fluids. For the quadratic one, we assume two equations of state corresponding to dark matter. In this model we obtain an inflation model and show that the values of the parameter, in the early universe, ɛ are depend on the sign of the cosmological constant.
Quadratic Herman-Wallis factors in the fundamental bands of linear molecules
NASA Astrophysics Data System (ADS)
Watson, James K. G.
1987-10-01
General theoretical formulas are derived for the coefficients in the terms M˜12 and M˜13 of the effective molecular dipole moment operator, and applied to the parallel and perpendicular fundamentals of linear molecules. The Herman-Wallis factors for P- and R-branch lines are F PR = [1 + A 1m + A 2PRm 2] 2, m = δ J( J' + J″ + 1)/2 and for Q-branch lines F Q = [1 + A 2QJ ( J + 1)] 2 The quadratic coefficients A2PR and A2Q depend on up to cubic potential derivatives and quadratic dipole derivatives. Calculated A2PR and A2Q values for the fundamentals of CO 2 do not agree well with recent measurements of Johns, and possible reasons for the discrepancies are discussed.
Steering of Frequency Standards by the Use of Linear Quadratic Gaussian Control Theory
NASA Technical Reports Server (NTRS)
Koppang, Paul; Leland, Robert
1996-01-01
Linear quadratic Gaussian control is a technique that uses Kalman filtering to estimate a state vector used for input into a control calculation. A control correction is calculated by minimizing a quadratic cost function that is dependent on both the state vector and the control amount. Different penalties, chosen by the designer, are assessed by the controller as the state vector and control amount vary from given optimal values. With this feature controllers can be designed to force the phase and frequency differences between two standards to zero either more or less aggressively depending on the application. Data will be used to show how using different parameters in the cost function analysis affects the steering and the stability of the frequency standards.
Schneider, Uwe
2009-04-15
A simple mechanistic model for predicting cancer induction after fractionated radiotherapy is developed. The model is based upon the linear-quadratic model. The inductions of carcinomas and sarcomas are modeled separately. The linear-quadratic model of cell kill is applied to normal tissues which are unintentionally irradiated during a cancer treatment with radiotherapy. Tumor induction is modeled such that each transformation process results in a tumor cell. The microscopic transformation parameter was chosen such that in the limit of low dose and acute exposure, the parameters of the linear-no-threshold model for tumor induction were approached. The differential equations describing carcinoma and sarcoma inductions can be solved analytically. Cancer induction in this model is a function of treatment dose, the cell kill parameters ({alpha},{beta}), the tumor induction variable ({mu}), and the repopulation parameter ({xi}). Carcinoma induction shows a bell shaped behavior as long as cell repopulation is small. Assuming large cell repopulation rates, a plateaulike function is approached. In contrast, sarcoma induction is negligible for low doses and increases with increasing dose up to a constant value. The proposed model describes carcinoma and sarcoma inductions after fractionated radiotherapy as an analytical function of four parameters. In the limit of low dose and for an instant irradiation it reproduces the results of the linear-no-threshold model. The obtained dose-response curves for cancer induction can be implemented with other models such as the organ-equivalent dose model to predict second cancers after radiotherapy.
Gantry cranes gain scheduling feedback control with friction compensation
NASA Astrophysics Data System (ADS)
Omar, Hanafy M.; Nayfeh, Ali H.
2005-03-01
We designed a controller based on gain-scheduling feedback to move a load on a gantry crane from point to point within one oscillation cycle and without inducing large swings. The settling time of the system is taken to be equal to the period of oscillation of the load. This criterion enables calculation of the controller feedback gains for varying load weight and cable length. Numerical simulations show that the controller is effective in reducing load oscillations and transferring the load in a reasonable time compared with that of optimal control. To experimentally validate the theory, we had to compensate for friction. To this end, we estimated the friction, and then applied an opposite control action to cancel it. To estimate the friction force, we assumed a mathematical model, and then we estimated the model coefficients using an off-line identification technique, such as the method of least squares (LS). First, the process of identification is applied to a theoretical model of a DC motor with known friction coefficients. From this example, some guidelines and rules are deduced for the choice of the LS parameters. Then, the friction coefficients of the gantry crane model are estimated and validated.
Quadratic blind linear unmixing: A graphical user interface for tissue characterization.
Gutierrez-Navarro, O; Campos-Delgado, D U; Arce-Santana, E R; Jo, Javier A
2016-02-01
Spectral unmixing is the process of breaking down data from a sample into its basic components and their abundances. Previous work has been focused on blind unmixing of multi-spectral fluorescence lifetime imaging microscopy (m-FLIM) datasets under a linear mixture model and quadratic approximations. This method provides a fast linear decomposition and can work without a limitation in the maximum number of components or end-members. Hence this work presents an interactive software which implements our blind end-member and abundance extraction (BEAE) and quadratic blind linear unmixing (QBLU) algorithms in Matlab. The options and capabilities of our proposed software are described in detail. When the number of components is known, our software can estimate the constitutive end-members and their abundances. When no prior knowledge is available, the software can provide a completely blind solution to estimate the number of components, the end-members and their abundances. The characterization of three case studies validates the performance of the new software: ex-vivo human coronary arteries, human breast cancer cell samples, and in-vivo hamster oral mucosa. The software is freely available in a hosted webpage by one of the developing institutions, and allows the user a quick, easy-to-use and efficient tool for multi/hyper-spectral data decomposition. PMID:26589467
Quadratic Blind Linear Unmixing: A Graphical User Interface for Tissue Characterization
Gutierrez-Navarro, O.; Campos-Delgado, D.U.; Arce-Santana, E. R.; Jo, Javier A.
2016-01-01
Spectral unmixing is the process of breaking down data from a sample into its basic components and their abundances. Previous work has been focused on blind unmixing of multi-spectral fluorescence lifetime imaging microscopy (m-FLIM) datasets under a linear mixture model and quadratic approximations. This method provides a fast linear decomposition and can work without a limitation in the maximum number of components or end-members. Hence this work presents an interactive software which implements our blind end-member and abundance extraction (BEAE) and quadratic blind linear unmixing (QBLU) algorithms in Matlab. The options and capabilities of our proposed software are described in detail. When the number of components is known, our software can estimate the constitutive end-members and their abundances. When no prior knowledge is available, the software can provide a completely blind solution to estimate the number of components, the end-members and their abundances. The characterization of three case studies validates the performance of the new software: ex-vivo human coronary arteries, human breast cancer cell samples, and in-vivo hamster oral mucosa. The software is freely available in a hosted webpage by one of the developing institutions, and allows the user a quick, easy-to-use and efficient tool for multi/hyper-spectral data decomposition. PMID:26589467
NASA Technical Reports Server (NTRS)
Becus, G. A.; Lui, C. Y.; Venkayya, V. B.; Tischler, V. A.
1987-01-01
A method for simultaneous structural and control design of large flexible space structures (LFSS) to reduce vibration generated by disturbances is presented. Desired natural frequencies and damping ratios for the closed loop system are achieved by using a combination of linear quadratic regulator (LQR) synthesis and numerical optimization techniques. The state and control weighing matrices (Q and R) are expressed in terms of structural parameters such as mass and stiffness. The design parameters are selected by numerical optimization so as to minimize the weight of the structure and to achieve the desired closed-loop eigenvalues. An illustrative example of the design of a two bar truss is presented.
A new linear quadratic optimal controller for the 34-meter high efficiency antenna position loop
NASA Technical Reports Server (NTRS)
Nickerson, J. A.
1987-01-01
The design of a new position loop controller for the 34-meter High Efficiency Deep Space antennas using linear quadratic (LQ) optimal control techniques is discussed. The LQ optimal control theory is reviewed, and model development and verification are discussed. Families of optimal gain vectors are generated by varying weight parameters. Performance specifications were used to select a final gain vector. Estimator dynamics were selected and the corresponding gain vectors were computed. Final estimator selection was based on position, commanded rate, and estimator error responses.
NASA Technical Reports Server (NTRS)
Gibson, J. S.; Rosen, I. G.
1986-01-01
An abstract approximation theory and computational methods are developed for the determination of optimal linear-quadratic feedback control, observers and compensators for infinite dimensional discrete-time systems. Particular attention is paid to systems whose open-loop dynamics are described by semigroups of operators on Hilbert spaces. The approach taken is based on the finite dimensional approximation of the infinite dimensional operator Riccati equations which characterize the optimal feedback control and observer gains. Theoretical convergence results are presented and discussed. Numerical results for an example involving a heat equation with boundary control are presented and used to demonstrate the feasibility of the method.
A user oriented microcomputer facility for designing linear quadratic Gaussian feedback compensators
NASA Technical Reports Server (NTRS)
Houpt, P. K.; Wahid, J.; Johnson, T. L.; Ward, S. A.
1979-01-01
The paper describes a laboratory design facility for digital microprocessor implementation of Linear-Quadratic-Gaussian feedback compensators. Outputs from user interactive programs for solving infinite time horizon LQ regulator and Kalman filter problems are conditioned for implementation on a laboratory microcomputer system. The software consists of two parts: (1) an off-line high-level program for solving the LQ Ricatti equations and generating associated feedback and filter gains, and (2) a cross compiler/macro assembler which generates object code for the target microprocessor system. Application to the control of a two dimensional inverted pendulum and expanding the design/prototyping system to other target machine architectures are discussed.
ORACLS: A system for linear-quadratic-Gaussian control law design
NASA Technical Reports Server (NTRS)
Armstrong, E. S.
1978-01-01
A modern control theory design package (ORACLS) for constructing controllers and optimal filters for systems modeled by linear time-invariant differential or difference equations is described. Numerical linear-algebra procedures are used to implement the linear-quadratic-Gaussian (LQG) methodology of modern control theory. Algorithms are included for computing eigensystems of real matrices, the relative stability of a matrix, factored forms for nonnegative definite matrices, the solutions and least squares approximations to the solutions of certain linear matrix algebraic equations, the controllability properties of a linear time-invariant system, and the steady state covariance matrix of an open-loop stable system forced by white noise. Subroutines are provided for solving both the continuous and discrete optimal linear regulator problems with noise free measurements and the sampled-data optimal linear regulator problem. For measurement noise, duality theory and the optimal regulator algorithms are used to solve the continuous and discrete Kalman-Bucy filter problems. Subroutines are also included which give control laws causing the output of a system to track the output of a prescribed model.
Gain-scheduling multivariable LPV control of an irrigation canal system.
Bolea, Yolanda; Puig, Vicenç
2016-07-01
The purpose of this paper is to present a multivariable linear parameter varying (LPV) controller with a gain scheduling Smith Predictor (SP) scheme applicable to open-flow canal systems. This LPV controller based on SP is designed taking into account the uncertainty in the estimation of delay and the variation of plant parameters according to the operating point. This new methodology can be applied to a class of delay systems that can be represented by a set of models that can be factorized into a rational multivariable model in series with left/right diagonal (multiple) delays, such as, the case of irrigation canals. A multiple pool canal system is used to test and validate the proposed control approach. PMID:27184416
Robust Gain-Scheduled Fault Tolerant Control for a Transport Aircraft
NASA Technical Reports Server (NTRS)
Shin, Jong-Yeob; Gregory, Irene
2007-01-01
This paper presents an application of robust gain-scheduled control concepts using a linear parameter-varying (LPV) control synthesis method to design fault tolerant controllers for a civil transport aircraft. To apply the robust LPV control synthesis method, the nonlinear dynamics must be represented by an LPV model, which is developed using the function substitution method over the entire flight envelope. The developed LPV model associated with the aerodynamic coefficient uncertainties represents nonlinear dynamics including those outside the equilibrium manifold. Passive and active fault tolerant controllers (FTC) are designed for the longitudinal dynamics of the Boeing 747-100/200 aircraft in the presence of elevator failure. Both FTC laws are evaluated in the full nonlinear aircraft simulation in the presence of the elevator fault and the results are compared to show pros and cons of each control law.
McKenna, Frederick W.; Ahmad, Salahuddin
2011-01-01
The linear quadratic is the standard model for calculating isoeffects in the range of conventional dose per fraction. However, the use of hypofractionation and stereotactic body radiation therapy can call for isoeffect calculations for large doses per fraction. The purpose of this work is to investigate the linear quadratic at large doses per fraction. The linear quadratic is compared to models that incorporate effects such as dose protraction, whose purpose is to extend the useful range of the linear quadratic to larger doses. The linear quadratic and extended linear quadratic models are fit to 4 data sets. The model-predicted isoeffects for these data sets are calculated. It is found that the linear quadratic and extended linear quadratic predict different isoeffect curves for certain data sets. However, for these data sets, by appropriate selection of a α/β ratio, the linear quadratic can well approximate the extended linear quadratic models. In particular, it is found that a α/β ratio of 0.5 well approximates the extended linear quadratic isoeffect curve for 2 prostate cell lines for conventional and moderate doses per fraction. PMID:21731226
Approximation theory for LQG (Linear-Quadratic-Gaussian) optimal control of flexible structures
NASA Technical Reports Server (NTRS)
Gibson, J. S.; Adamian, A.
1988-01-01
An approximation theory is presented for the LQG (Linear-Quadratic-Gaussian) optimal control problem for flexible structures whose distributed models have bounded input and output operators. The main purpose of the theory is to guide the design of finite dimensional compensators that approximate closely the optimal compensator. The optimal LQG problem separates into an optimal linear-quadratic regulator problem and an optimal state estimation problem. The solution of the former problem lies in the solution to an infinite dimensional Riccati operator equation. The approximation scheme approximates the infinite dimensional LQG problem with a sequence of finite dimensional LQG problems defined for a sequence of finite dimensional, usually finite element or modal, approximations of the distributed model of the structure. Two Riccati matrix equations determine the solution to each approximating problem. The finite dimensional equations for numerical approximation are developed, including formulas for converting matrix control and estimator gains to their functional representation to allow comparison of gains based on different orders of approximation. Convergence of the approximating control and estimator gains and of the corresponding finite dimensional compensators is studied. Also, convergence and stability of the closed-loop systems produced with the finite dimensional compensators are discussed. The convergence theory is based on the convergence of the solutions of the finite dimensional Riccati equations to the solutions of the infinite dimensional Riccati equations. A numerical example with a flexible beam, a rotating rigid body, and a lumped mass is given.
NASA Technical Reports Server (NTRS)
Frost, Susan A.; Bodson, Marc; Acosta, Diana M.
2009-01-01
The Next Generation (NextGen) transport aircraft configurations being investigated as part of the NASA Aeronautics Subsonic Fixed Wing Project have more control surfaces, or control effectors, than existing transport aircraft configurations. Conventional flight control is achieved through two symmetric elevators, two antisymmetric ailerons, and a rudder. The five effectors, reduced to three command variables, produce moments along the three main axes of the aircraft and enable the pilot to control the attitude and flight path of the aircraft. The NextGen aircraft will have additional redundant control effectors to control the three moments, creating a situation where the aircraft is over-actuated and where a simple relationship does not exist anymore between the required effector deflections and the desired moments. NextGen flight controllers will incorporate control allocation algorithms to determine the optimal effector commands and attain the desired moments, taking into account the effector limits. Approaches to solving the problem using linear programming and quadratic programming algorithms have been proposed and tested. It is of great interest to understand their relative advantages and disadvantages and how design parameters may affect their properties. In this paper, we investigate the sensitivity of the effector commands with respect to the desired moments and show on some examples that the solutions provided using the l2 norm of quadratic programming are less sensitive than those using the l1 norm of linear programming.
Modelling Ocean Dissipation in Icy Satellites: A Comparison of Linear and Quadratic Friction
NASA Astrophysics Data System (ADS)
Hay, H.; Matsuyama, I.
2015-12-01
Although subsurface oceans are confirmed in Europa, Ganymede, Callisto, and strongly suspected in Enceladus and Titan, the exact mechanism required to heat and maintain these liquid reservoirs over Solar System history remains a mystery. Radiogenic heating can supply enough energy for large satellites whereas tidal dissipation provides the best explanation for the presence of oceans in small icy satellites. The amount of thermal energy actually contributed to the interiors of these icy satellites through oceanic tidal dissipation is largely unquantified. Presented here is a numerical model that builds upon previous work for quantifying tidally dissipated energy in the subsurface oceans of the icy satellites. Recent semi-analytical models (Tyler, 2008 and Matsuyama, 2014) have solved the Laplace Tidal Equations to estimate the time averaged energy flux over an orbital period in icy satellite oceans, neglecting the presence of a solid icy shell. These models are only able to consider linear Rayleigh friction. The numerical model presented here is compared to one of these semi-analytical models, finding excellent agreement between velocity and displacement solutions for all three terms to the tidal potential. Time averaged energy flux is within 2-6% of the analytical values. Quadratic (bottom) friction is then incorporated into the model, replacing linear friction. This approach is commonly applied to terrestrial ocean dissipation studies where dissipation scales nonlinearly with velocity. A suite of simulations are also run for the quadratic friction case which are then compared to and analysed against recent scaling laws developed by Chen and Nimmo (2013).
NASA Technical Reports Server (NTRS)
Turso, James A.; Litt, Jonathan S.
2004-01-01
A method for accommodating engine deterioration via a scheduled Linear Parameter Varying Quadratic Lyapunov Function (LPVQLF)-Based controller is presented. The LPVQLF design methodology provides a means for developing unconditionally stable, robust control of Linear Parameter Varying (LPV) systems. The controller is scheduled on the Engine Deterioration Index, a function of estimated parameters that relate to engine health, and is computed using a multilayer feedforward neural network. Acceptable thrust response and tight control of exhaust gas temperature (EGT) is accomplished by adjusting the performance weights on these parameters for different levels of engine degradation. Nonlinear simulations demonstrate that the controller achieves specified performance objectives while being robust to engine deterioration as well as engine-to-engine variations.
Regions of attraction and ultimate boundedness for linear quadratic regulators with nonlinearities
NASA Technical Reports Server (NTRS)
Joshi, S. M.
1984-01-01
The closed-loop stability of multivariable linear time-invariant systems controlled by optimal linear quadratic (LQ) regulators is investigated for the case when the feedback loops have nonlinearities N(sigma) that violate the standard stability condition, sigma N(sigma) or = 0.5 sigma(2). The violations of the condition are assumed to occur either (1) for values of sigma away from the origin (sigma = 0) or (2) for values of sigma in a neighborhood of the origin. It is proved that there exists a region of attraction for case (1) and a region of ultimate boundedness for case (2), and estimates are obtained for these regions. The results provide methods for selecting the performance function parameters to design LQ regulators with better tolerance to nonlinearities. The results are demonstrated by application to the problem of attitude and vibration control of a large, flexible space antenna in the presence of actuator nonlinearities.
NASA Astrophysics Data System (ADS)
Zhu, Huai-nian; Zhang, Cheng-ke; Bin, Ning
2014-05-01
In this paper, we deal with the problem of stochastic Nash differential games of Markov jump linear systems governed by Itô-type equation. Combining the stochastic stabilizability with the stochastic systems, a necessary and sufficient condition for the existence of the Nash strategy is presented by means of a set of cross-coupled stochastic algebraic Riccati equations. Moreover, the stochastic H2/H∞ control for stochastic Markov jump linear systems is discussed as an immediate application and an illustrative example is presented.
Xia, Youshen; Feng, Gang; Wang, Jun
2004-09-01
This paper presents a recurrent neural network for solving strict convex quadratic programming problems and related linear piecewise equations. Compared with the existing neural networks for quadratic program, the proposed neural network has a one-layer structure with a low model complexity. Moreover, the proposed neural network is shown to have a finite-time convergence and exponential convergence. Illustrative examples further show the good performance of the proposed neural network in real-time applications. PMID:15312842
NASA Technical Reports Server (NTRS)
Young, J. W.; Hamer, H. A.; Johnson, K. G.
1984-01-01
A decoupled-control analysis was performed for a large flexible space antenna. Control involved commanding changes in the rigid-body modes or nulling disturbances in the flexible modes. The study provides parametric-type data which could be useful in the final design of a large space antenna control system. Results are presented to illustrate the effect on control requirements of (1) the number of modes controlled; (2) the number, type, and location of control actuators; and (3) variations in the closed-loop dynamics of the control system. Comparisons are given between the decoupled-control results and those obtained by using a linear quadratic regulator approach. Time history responses are presented to illustrate the effects of the control procedures.
NASA Technical Reports Server (NTRS)
Milman, Mark H.
1987-01-01
The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary systems. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.
A user oriented microcomputer facility for designing linear quadratic Gaussian feedback compensators
NASA Technical Reports Server (NTRS)
Houpt, P. K.; Wahid, J.; Johnson, T. L.; Ward, S. A.
1978-01-01
A laboratory design facility for digital microprocessor implementation of linear-quadratic-Gaussian feedback compensators is described. Outputs from user interactive programs for solving infinite time horizon LQ regulator and Kalman filter problems were conditioned for implementation on the laboratory microcomputer system. The software consisted of two parts: an offline high-level program for solving the LQ Ricatti equations and generating associated feedback and filter gains and a cross compiler/macro assembler which generates object code for the target microprocessor system. A PDP 11/70 with a UNIX operating system was used for all high level program and data management, and the target microprocessor system is an Intel MDS (8080-based processor). Application to the control of a two dimensional inverted pendulum is presented and issues in expanding the design/prototyping system to other target machine architectures are discussed.
Frequency locking of an optical cavity using linear-quadratic Gaussian integral control
NASA Astrophysics Data System (ADS)
Sayed Hassen, S. Z.; Heurs, M.; Huntington, E. H.; Petersen, I. R.; James, M. R.
2009-09-01
We show that a systematic modern control technique such as linear-quadratic Gaussian (LQG) control can be applied to a problem in experimental quantum optics which has previously been addressed using traditional approaches to controller design. An LQG controller which includes integral action is synthesized to stabilize the frequency of the cavity to the laser frequency and to reject low frequency noise. The controller is successfully implemented in the laboratory using a dSpace digital signal processing board. One important advantage of the LQG technique is that it can be extended in a straightforward way to control systems with multiple measurements and multiple feedback loops. This work is expected to pave the way for extremely stable lasers with fluctuations approaching the quantum noise limit and which could be potentially used in a wide range of applications.
NASA Astrophysics Data System (ADS)
Aras, Mohd Shahrieel Mohd; Abdullah, Shahrum Shah; Kamarudin, Muhammad Nizam; Rahman, Ahmad Fadzli Nizam Abdul; Azis, Fadilah Abd; Jaafar, Hazriq Izzuan
2015-05-01
This paper describes the effectiveness of observer-based output feedback for Unmanned Underwater Vehicle (UUV) with Linear Quadratic Regulation (LQR) performance. Tuning of observer parameters is crucial for tracking purpose. Prior to tuning facility, the ranges of observer and LQR parameters are obtained via system output cum error. The validation of this technique using unmanned underwater vehicles called Remotely Operated Vehicle (ROV) modelling helps to improve steady state performance of system response. The ROV modeling is focused for depth control using ROV 1 developed by the Underwater Technology Research Group (UTeRG). The results are showing that this technique improves steady state performances in term of overshoot and settling time of the system response.
Mohammed, W; Meier, J; Galle, M; Qian, L; Aitchison, J S; Smith, P W E
2007-11-15
We measured linear and quadratic dispersion on millimeter-length fibers, waveguides, and nanowires based on common-path spectral interferometry. We obtained the linear dispersion parameter, beta', with a relative precision of 1.45 x 10(-4), and extracted the quadratic dispersion parameter, beta'', from the Taylor expansion of beta' x beta'' values show a discrepancy of < 1% when compared with simulation as well as with measurement results obtained by a conventional Michelson interferometer. Using this method, we experimentally confirmed the sign inversion of the group velocity dispersion of AlGaAs nanowires for what is believed to be the first time. PMID:18026291
NASA Technical Reports Server (NTRS)
Gibson, J. S.; Rosen, I. G.
1985-01-01
In the optimal linear quadratic regulator problem for finite dimensional systems, the method known as an alpha-shift can be used to produce a closed-loop system whose spectrum lies to the left of some specified vertical line; that is, a closed-loop system with a prescribed degree of stability. This paper treats the extension of the alpha-shift to hereditary systems. As infinite dimensions, the shift can be accomplished by adding alpha times the identity to the open-loop semigroup generator and then solving an optimal regulator problem. However, this approach does not work with a new approximation scheme for hereditary control problems recently developed by Kappel and Salamon. Since this scheme is among the best to date for the numerical solution of the linear regulator problem for hereditary systems, an alternative method for shifting the closed-loop spectrum is needed. An alpha-shift technique that can be used with the Kappel-Salamon approximation scheme is developed. Both the continuous-time and discrete-time problems are considered. A numerical example which demonstrates the feasibility of the method is included.
NASA Technical Reports Server (NTRS)
Gibson, J. S.; Rosen, I. G.
1987-01-01
In the optimal linear quadratic regulator problem for finite dimensional systems, the method known as an alpha-shift can be used to produce a closed-loop system whose spectrum lies to the left of some specified vertical line; that is, a closed-loop system with a prescribed degree of stability. This paper treats the extension of the alpha-shift to hereditary systems. As infinite dimensions, the shift can be accomplished by adding alpha times the identity to the open-loop semigroup generator and then solving an optimal regulator problem. However, this approach does not work with a new approximation scheme for hereditary control problems recently developed by Kappel and Salamon. Since this scheme is among the best to date for the numerical solution of the linear regulator problem for hereditary systems, an alternative method for shifting the closed-loop spectrum is needed. An alpha-shift technique that can be used with the Kappel-Salamon approximation scheme is developed. Both the continuous-time and discrete-time problems are considered. A numerical example which demonstrates the feasibility of the method is included.
NASA Technical Reports Server (NTRS)
Lehtinen, B.; Geyser, L. C.
1984-01-01
AESOP is a computer program for use in designing feedback controls and state estimators for linear multivariable systems. AESOP is meant to be used in an interactive manner. Each design task that the program performs is assigned a "function" number. The user accesses these functions either (1) by inputting a list of desired function numbers or (2) by inputting a single function number. In the latter case the choice of the function will in general depend on the results obtained by the previously executed function. The most important of the AESOP functions are those that design,linear quadratic regulators and Kalman filters. The user interacts with the program when using these design functions by inputting design weighting parameters and by viewing graphic displays of designed system responses. Supporting functions are provided that obtain system transient and frequency responses, transfer functions, and covariance matrices. The program can also compute open-loop system information such as stability (eigenvalues), eigenvectors, controllability, and observability. The program is written in ANSI-66 FORTRAN for use on an IBM 3033 using TSS 370. Descriptions of all subroutines and results of two test cases are included in the appendixes.
Linear quadratic Gaussian and feedforward controllers for the DSS-13 antenna
NASA Technical Reports Server (NTRS)
Gawronski, W. K.; Racho, C. S.; Mellstrom, J. A.
1994-01-01
The controller development and the tracking performance evaluation for the DSS-13 antenna are presented. A trajectory preprocessor, linear quadratic Gaussian (LQG) controller, feedforward controller, and their combination were designed, built, analyzed, and tested. The antenna exhibits nonlinear behavior when the input to the antenna and/or the derivative of this input exceeds the imposed limits; for slewing and acquisition commands, these limits are typically violated. A trajectory preprocessor was designed to ensure that the antenna behaves linearly, just to prevent nonlinear limit cycling. The estimator model for the LQG controller was identified from the data obtained from the field test. Based on an LQG balanced representation, a reduced-order LQG controller was obtained. The feedforward controller and the combination of the LQG and feedforward controller were also investigated. The performance of the controllers was evaluated with the tracking errors (due to following a trajectory) and the disturbance errors (due to the disturbances acting on the antenna). The LQG controller has good disturbance rejection properties and satisfactory tracking errors. The feedforward controller has small tracking errors but poor disturbance rejection properties. The combined LQG and feedforward controller exhibits small tracking errors as well as good disturbance rejection properties. However, the cost for this performance is the complexity of the controller.
Neural network-based nonlinear model predictive control vs. linear quadratic gaussian control
Cho, C.; Vance, R.; Mardi, N.; Qian, Z.; Prisbrey, K.
1997-01-01
One problem with the application of neural networks to the multivariable control of mineral and extractive processes is determining whether and how to use them. The objective of this investigation was to compare neural network control to more conventional strategies and to determine if there are any advantages in using neural network control in terms of set-point tracking, rise time, settling time, disturbance rejection and other criteria. The procedure involved developing neural network controllers using both historical plant data and simulation models. Various control patterns were tried, including both inverse and direct neural network plant models. These were compared to state space controllers that are, by nature, linear. For grinding and leaching circuits, a nonlinear neural network-based model predictive control strategy was superior to a state space-based linear quadratic gaussian controller. The investigation pointed out the importance of incorporating state space into neural networks by making them recurrent, i.e., feeding certain output state variables into input nodes in the neural network. It was concluded that neural network controllers can have better disturbance rejection, set-point tracking, rise time, settling time and lower set-point overshoot, and it was also concluded that neural network controllers can be more reliable and easy to implement in complex, multivariable plants.
Gain-scheduled H∞ control via parameter-dependent Lyapunov functions
NASA Astrophysics Data System (ADS)
Chumalee, Sunan; Whidborne, James F.
2015-01-01
Synthesising a gain-scheduled output feedback H∞ controller via parameter-dependent Lyapunov functions for linear parameter-varying (LPV) plant models involves solving an infinite number of linear matrix inequalities (LMIs). In practice, for affine LPV models, a finite number of LMIs can be achieved using convexifying techniques. This paper proposes an alternative approach to achieve a finite number of LMIs. By simple manipulations on the bounded real lemma inequality, a symmetric matrix polytope inequality can be formed. Hence, the LMIs need only to be evaluated at all vertices of such a symmetric matrix polytope. In addition, a construction technique of the intermediate controller variables is also proposed as an affine matrix-valued function in the polytopic coordinates of the scheduled parameters. Computational results on a numerical example using the approach were compared with those from a multi-convexity approach in order to demonstrate the impacts of the approach on parameter-dependent Lyapunov-based stability and performance analysis. Furthermore, numerical simulation results show the effectiveness of these proposed techniques.
Robust gain-scheduling for smart-structures in parallel robots
NASA Astrophysics Data System (ADS)
Algermissen, Stephan; Rose, Michael; Keimer, Ralf; Sinapius, Michael
2009-03-01
In the past years parallel robots demonstrated their capability in applications with high-dynamic trajectories. Smart-structures offer the potential to further increase the productivity of parallel robots by reducing disturbing vibrations caused by high dynamic loads effectively. To investigate parallel robots and their applications, including suitable control concepts for smart-structures, the Collaborative Research Center 562 was founded by the German Research Council (DFG). The latest prototype within this research center is called Triglide. It is a four degree of freedom (DOF) robot with three translational and one rotational DOF. It realizes an acceleration of 10 g* at the effector. In the structure of the robot six active rods and a tri-axial accelerometer are integrated to control effector vibrations in three translational DOF. The main challenge of this control application is the position dependent vibration behavior. A single robust controller is not able to gain satisfying performance within the entire workspace. Therefore a strategy for describing the vibration behavior by linearization at several operating points is developed. Behavior in-between is approximated by a linear approach. On a trajectory robust controllers in all operating points are smoothly switched by robust gain-scheduling. The scheduling parameters are fast varying and though a suitable stability proof is defined, based on Small-Gain approach. Several transformations enhance the results from Small-Gain Theorem and reduce the usual conservatism. Experimental data is used to show the improvements made.
ERIC Educational Resources Information Center
Molenje, Levi
2012-01-01
This study was designed to explore secondary mathematics teachers' beliefs about graphing calculators, their practices with the graphing calculators when teaching linear and quadratic functions, and the relationship between the teachers' beliefs and their practices. The study was conducted in two phases. In the first phase, 81 teachers…
A nanodosimetry-based linear-quadratic model of cell survival for mixed-LET radiations
NASA Astrophysics Data System (ADS)
Wang, C.-K. Chris; Zhang, Xin
2006-12-01
A new nanodosimetry-based linear-quadratic (LQ) model of cell survival for mixed-LET radiations has been developed. The new model employs three physical quantities and three biological quantities. The three physical quantities are related to energy depositions at two nanometre scales, 5 nm and 25 nm. The three biological quantities are related to the lesion production and interaction probabilities and the lesion repair rate. The coefficients α and β of the LQ formula (αD + βD2) are explicitly expressed in terms of the three physical quantities and the three biological quantities. The new model is shown to be consistent with the previously published cell survival curves of V-79 cells. The advantage of this new model is that it can be conveniently adopted to estimate the iso-effect for radiotherapies that involve ionizing radiation of mixed LET. An example is given to estimate the cell survival fractions for a high-dose-rate mixed neutron and gamma-ray field from a 252Cf source.
Spacecraft Formation Flying Maneuvers Using Linear Quadratic Regulation With No Radial Axis Inputs
NASA Technical Reports Server (NTRS)
Starin, Scott R.; Yedavalli, R. K.; Sparks, Andrew G.; Bauer, Frank H. (Technical Monitor)
2001-01-01
Regarding multiple spacecraft formation flying, the observation has been made that control thrust need only be applied coplanar to the local horizon to achieve complete controllability of a two-satellite (leader-follower) formation. A formulation of orbital dynamics using the state of one satellite relative to another is used. Without the need for thrust along the radial (zenith-nadir) axis of the relative reference frame, propulsion system simplifications and weight reduction may be accomplished. This work focuses on the validation of this control system on its own merits, and in comparison to a related system which does provide thrust along the radial axis of the relative frame. Maneuver simulations are performed using commercial ODE solvers to propagate the Keplerian dynamics of a controlled satellite relative to an uncontrolled leader. These short maneuver simulations demonstrate the capacity of the controller to perform changes from one formation geometry to another. Control algorithm performance is evaluated based on measures such as the fuel required to complete a maneuver and the maximum acceleration required by the controller. Based on this evaluation, the exclusion of the radial axis of control still allows enough control authority to use Linear Quadratic Regulator (LQR) techniques to design a gain matrix of adequate performance over finite maneuvers. Additional simulations are conducted including perturbations and using no radial control inputs. A major conclusion presented is that control inputs along the three axes have significantly different relationships to the governing orbital dynamics that may be exploited using LQR.
NASA Astrophysics Data System (ADS)
Chun, Tae Yoon; Lee, Jae Young; Park, Jin Bae; Choi, Yoon Ho
2016-03-01
In this paper, we analyse the convergence and stability properties of generalised policy iteration (GPI) applied to discrete-time linear quadratic regulation problems. GPI is one kind of the generalised adaptive dynamic programming methods used for solving optimal control problems, and is composed of policy evaluation and policy improvement steps. To analyse the convergence and stability of GPI, the dynamic programming (DP) operator is defined. Then, GPI and its equivalent formulas are presented based on the notation of DP operator. The convergence of the approximate value function to the exact one in policy evaluation is proven based on the equivalent formulas. Furthermore, the positive semi-definiteness, stability, and the monotone convergence (PI-mode and VI-mode convergence) of GPI are presented under certain conditions on the initial value function. The online least square method is also presented for the implementation of GPI. Finally, some numerical simulations are carried out to verify the effectiveness of GPI as well as to further investigate the convergence and stability properties.
Linear quadratic modeling of increased late normal-tissue effects in special clinical situations
Jones, Bleddyn . E-mail: b.jones.1@bham.ac.uk; Dale, Roger G.; Gaya, Andrew M.
2006-03-01
Purpose: To extend linear quadratic theory to allow changes in normal-tissue radiation tolerance after exposure to cytotoxic chemotherapy, after surgery, and in elderly patients. Methods: Examples of these situations are analyzed by use of the biologic effective dose (BED) concept. Changes in tolerance can be allowed for by: estimation of either the contribution of the additional factor as an equivalent BED or the equivalent dose in 2-Gy fractions or by the degree of radiosensitization by a mean dose-modifying factor (x). Results: The estimated x value is 1.063 (95% confidence limits for the mean, 1.056 to 1.070) for subcutaneous fibrosis after cyclophosphamide, methotrexate, and fluorouracil (CMF) chemotherapy and radiotherapy in breast cancer. The point estimate of x is 1.18 for the additional risk of gastrointestinal late-radiation effects after abdominal surgery in lymphoma patients (or 10.62 Gy at 2 Gy per fraction). For shoulder fibrosis in patients older than 60 years after breast and nodal irradiation, x is estimated to be 1.033 (95% confidence limits for the mean, 1.028 to 1.0385). The equivalent BED values were CMF chemotherapy (6.48 Gy{sub 3}), surgery (17.73 Gy{sub 3}), and age (3.61 Gy{sub 3}). Conclusions: The LQ model can, in principle, be extended to quantify reduced normal-tissue tolerance in special clinical situations.
NASA Technical Reports Server (NTRS)
Gawronski, W.
2004-01-01
Wind gusts are the main disturbances that depreciate tracking precision of microwave antennas and radiotelescopes. The linear-quadratic-Gaussian (LQG) controllers - as compared with the proportional-and-integral (PI) controllers significantly improve the tracking precision in wind disturbances. However, their properties have not been satisfactorily understood; consequently, their tuning is a trial-and-error process. A control engineer has two tools to tune an LQG controller: the choice of coordinate system of the controller model and the selection of weights of the LQG performance index. This article analyzes properties of an open- and closed-loop antenna. It shows that the proper choice of coordinates of the open-loop model simplifies the shaping of the closed-loop performance. The closed-loop properties are influenced by the LQG weights. The article shows the impact of the weights on the antenna closed-loop bandwidth, disturbance rejection properties, and antenna acceleration. The bandwidth and the disturbance rejection characterize the antenna performance, while the acceleration represents the performance limit set by the antenna hardware (motors). The article presents the controller tuning procedure, based on the coordinate selection and the weight properties. The procedure rationally shapes the closed-loop performance, as an alternative to the trial-and-error approach.
Linear Quadratic Tracking Design for a Generic Transport Aircraft with Structural Load Constraints
NASA Technical Reports Server (NTRS)
Burken, John J.; Frost, Susan A.; Taylor, Brian R.
2011-01-01
When designing control laws for systems with constraints added to the tracking performance, control allocation methods can be utilized. Control allocations methods are used when there are more command inputs than controlled variables. Constraints that require allocators are such task as; surface saturation limits, structural load limits, drag reduction constraints or actuator failures. Most transport aircraft have many actuated surfaces compared to the three controlled variables (such as angle of attack, roll rate & angle of side slip). To distribute the control effort among the redundant set of actuators a fixed mixer approach can be utilized or online control allocation techniques. The benefit of an online allocator is that constraints can be considered in the design whereas the fixed mixer cannot. However, an online control allocator mixer has a disadvantage of not guaranteeing a surface schedule, which can then produce ill defined loads on the aircraft. The load uncertainty and complexity has prevented some controller designs from using advanced allocation techniques. This paper considers actuator redundancy management for a class of over actuated systems with real-time structural load limits using linear quadratic tracking applied to the generic transport model. A roll maneuver example of an artificial load limit constraint is shown and compared to the same no load limitation maneuver.
A nanodosimetry-based linear-quadratic model of cell survival for mixed-LET radiations.
Chris Wang, C-K; Zhang, Xin
2006-12-01
A new nanodosimetry-based linear-quadratic (LQ) model of cell survival for mixed-LET radiations has been developed. The new model employs three physical quantities and three biological quantities. The three physical quantities are related to energy depositions at two nanometre scales, 5 nm and 25 nm. The three biological quantities are related to the lesion production and interaction probabilities and the lesion repair rate. The coefficients alpha and beta of the LQ formula (alpha D + beta D(2)) are explicitly expressed in terms of the three physical quantities and the three biological quantities. The new model is shown to be consistent with the previously published cell survival curves of V-79 cells. The advantage of this new model is that it can be conveniently adopted to estimate the iso-effect for radiotherapies that involve ionizing radiation of mixed LET. An example is given to estimate the cell survival fractions for a high-dose-rate mixed neutron and gamma-ray field from a (252)Cf source. PMID:17110772
Spacecraft Formation Flying Maneuvers Using Linear-Quadratic Regulation with No Radial Axis Inputs
NASA Technical Reports Server (NTRS)
Starin, Scott R.; Yedavalli, R. K.; Sparks, Andrew G.; Bauer, Frank H. (Technical Monitor)
2001-01-01
Regarding multiple spacecraft formation flying, the observation has been made that control thrust need only be applied coplanar to the local horizon to achieve complete controllability of a two-satellite (leader-follower) formation. A formulation of orbital dynamics using the state of one satellite relative to another is used. Without the need for thrust along the radial (zenith-nadir) axis of the relative reference frame ' propulsion system simplifications and weight reduction may be accomplished. Several linear-quadratic regulators (LQR) are explored and compared based on performance measures likely to be important to many missions, but not directly optimized in the LQR designs. Maneuver simulations are performed using commercial ODE solvers to propagate the Keplerian dynamics of a controlled satellite relative to an uncontrolled leader. These short maneuver simulations demonstrate the capacity of the controller to perform changes from one formation geometry to another. This work focusses on formations in which the controlled satellite has a relative trajectory which projects onto the local horizon of the uncontrolled satellite as a circle. This formation has potential uses for distributed remote sensing systems.
Fitting the linear quadratic model to detailed data sets for different dose ranges
NASA Astrophysics Data System (ADS)
Garcia, L. M.; Leblanc, J.; Wilkins, D.; Raaphorst, G. P.
2006-06-01
Survival curve behaviour and degree of correspondence between the linear-quadratic (LQ) model and experimental data in an extensive dose range for high dose rates were analysed. Detailed clonogenic assays with irradiation given in 0.5 Gy increments and a total dose range varying from 10.5 to 16 Gy were performed. The cell lines investigated were: CHOAA8 (Chinese hamster fibroblast cells), U373MG (human glioblastoma cells), CP3 and DU145 (human prostate carcinoma cell lines). The analyses were based on χ2-statistics and Monte Carlo simulation of the experiments. A decline of LQ fit quality at very low doses (<2 Gy) is observed. This result can be explained by the hypersensitive effect observed in CHOAA8, U373MG and DU145 data and an adaptive-type response in the CP3 cell line. A clear improvement of the fit is discerned by removing the low dose data points. The fit worsening at high doses also shows that LQ cannot explain this region. This shows that the LQ model fits better the middle dose region of the survival curve. The analysis conducted in our study reveals a dose dependency of the LQ fit in different cell lines.
Irrigation Control in the Presence of Salinity: Extended Linear Quadratic Approach
NASA Astrophysics Data System (ADS)
Bras, Rafael L.; Seo, Dong-Jun
1987-07-01
An intraseasonal irrigation scheduling problem is dealt with via extended linear quadratic (ELQ) control. The ELQ control is well-suited for constrained multidimensional problems and provides openloop feedback control rules over the control horizon. A conceptual model is developed to describe the dynamics of water allocation and salt movement in the root zone of a crop. Moisture stress and osmotic stress are combined to obtain the integrated inhibitory effect of salinity on transpiration. For the intraseasonal model to be effective against perennial salt accumulation in the root zone, it should be able to yield control laws which will lead to favorable root zone conditions at the end of an irrigation season, thus avoiding any significant leaching prior to the next growing season. This long-term aspect of salinity control is handled via probabilistic state constraints which impose desired salinity and moisture levels with desired confidence level. The ELQ control is employed in a case study of expected net benefit maximization over an irrigation season of corn in Fort Morgan, Colorado. The results, in general, correspond well with expected irrigation schedules under different conditions and provide valuable information on both short- and long-term aspects of irrigation control under saline conditions. The ELQ control, being an analytic iterative solution scheme with theoretically guaranteed fast convergence, has a distinct computational advantage over state-of-the-art procedures.
NASA Astrophysics Data System (ADS)
Xu, Jiuping; Li, Jun
2002-09-01
In this paper a class of stochastic multiple-objective programming problems with one quadratic, several linear objective functions and linear constraints has been introduced. The former model is transformed into a deterministic multiple-objective nonlinear programming model by means of the introduction of random variables' expectation. The reference direction approach is used to deal with linear objectives and results in a linear parametric optimization formula with a single linear objective function. This objective function is combined with the quadratic function using the weighted sums. The quadratic problem is transformed into a linear (parametric) complementary problem, the basic formula for the proposed approach. The sufficient and necessary conditions for (properly, weakly) efficient solutions and some construction characteristics of (weakly) efficient solution sets are obtained. An interactive algorithm is proposed based on reference direction and weighted sums. Varying the parameter vector on the right-hand side of the model, the DM can freely search the efficient frontier with the model. An extended portfolio selection model is formed when liquidity is considered as another objective to be optimized besides expectation and risk. The interactive approach is illustrated with a practical example.
NASA Astrophysics Data System (ADS)
Kheiri, R.
2016-09-01
As an undergraduate exercise, in an article (2012 Am. J. Phys. 80 780–14), quantum and classical uncertainties for dimensionless variables of position and momentum were evaluated in three potentials: infinite well, bouncing ball, and harmonic oscillator. While original quantum uncertainty products depend on {{\\hslash }} and the number of states (n), a dimensionless approach makes the comparison between quantum uncertainty and classical dispersion possible by excluding {{\\hslash }}. But the question is whether the uncertainty still remains dependent on quantum number n. In the above-mentioned article, there lies this contrast; on the one hand, the dimensionless quantum uncertainty of the potential box approaches classical dispersion only in the limit of large quantum numbers (n\\to ∞ )—consistent with the correspondence principle. On the other hand, similar evaluations for bouncing ball and harmonic oscillator potentials are equal to their classical counterparts independent of n. This equality may hide the quantum feature of low energy levels. In the current study, we change the potential intervals in order to make them symmetric for the linear potential and non-symmetric for the quadratic potential. As a result, it is shown in this paper that the dimensionless quantum uncertainty of these potentials in the new potential intervals is expressed in terms of quantum number n. In other words, the uncertainty requires the correspondence principle in order to approach the classical limit. Therefore, it can be concluded that the dimensionless analysis, as a useful pedagogical method, does not take away the quantum feature of the n-dependence of quantum uncertainty in general. Moreover, our numerical calculations include the higher powers of the position for the potentials.
LPV gain-scheduled control of SCR aftertreatment systems
NASA Astrophysics Data System (ADS)
Meisami-Azad, Mona; Mohammadpour, Javad; Grigoriadis, Karolos M.; Harold, Michael P.; Franchek, Matthew A.
2012-01-01
Hydrocarbons, carbon monoxide and some of other polluting emissions produced by diesel engines are usually lower than those produced by gasoline engines. While great strides have been made in the exhaust aftertreatment of vehicular pollutants, the elimination of nitrogen oxide (NO x ) from diesel vehicles is still a challenge. The primary reason is that diesel combustion is a fuel-lean process, and hence there is significant unreacted oxygen in the exhaust. Selective catalytic reduction (SCR) is a well-developed technology for power plants and has been recently employed for reducing NO x emissions from automotive sources and in particular, heavy-duty diesel engines. In this article, we develop a linear parameter-varying (LPV) feedforward/feedback control design method for the SCR aftertreatment system to decrease NO x emissions while keeping ammonia slippage to a desired low level downstream the catalyst. The performance of the closed-loop system obtained from the interconnection of the SCR system and the output feedback LPV control strategy is then compared with other control design methods including sliding mode, and observer-based static state-feedback parameter-varying control. To reduce the computational complexity involved in the control design process, the number of LPV parameters in the developed quasi-LPV (qLPV) model is reduced by applying the principal component analysis technique. An LPV feedback/feedforward controller is then designed for the qLPV model with reduced number of scheduling parameters. The designed full-order controller is further simplified to a first-order transfer function with a parameter-varying gain and pole. Finally, simulation results using both a low-order model and a high-fidelity and high-order model of SCR reactions in GT-POWER interfaced with MATLAB/SIMULINK illustrate the high NO x conversion efficiency of the closed-loop SCR system using the proposed parameter-varying control law.
Robust gain-scheduling energy-to-peak control of vehicle lateral dynamics stabilisation
NASA Astrophysics Data System (ADS)
Zhang, Hui; Zhang, Xinjie; Wang, Junmin
2014-03-01
In this paper, we investigate the vehicle lateral dynamics stabilisation problem to enhance vehicle handling by considering time-varying longitudinal velocity. The longitudinal velocity is described by a polytope with finite vertices and a novel technique is proposed to reduce the number of vertices. Since the tyre dynamics is nonlinear, the cornering stiffness is represented via the norm-bounded uncertainty. Concerning the time-varying velocity and the nonlinear tyre model, a linear parameter-varying vehicle model is obtained. As the velocity and the states are measurable, a gain-scheduling state-feedback controller is introduced. In the lateral control, the sideslip angle is required to be as small as possible and the yaw rate is constrained to a certain level. Thus, the control objective is to minimise the sideslip angle while the yaw rate is under a prescribed level or constrain both the sideslip angle and the yaw rate to prescribed levels. To consider the transient response of the closed-loop system, the ?-stability is also employed in the energy-to-peak control. The optimal controller can be obtained by solving a set of linear matrix inequalities. A nonlinear vehicle model is utilised to illustrate the design procedure and the effectiveness of the proposed design method. Finally, simulations and comparisons are carried out to show the significant advantage of the designed controller. Compared to the open-loop system, the closed-loop system with the designed controller can achieve much smaller sideslip angle and the yaw rate is closer to the desired yaw rate from a reference model. Therefore, the vehicle safety and the handling are both improved in our simulation cases.
NASA Technical Reports Server (NTRS)
Athans, M.; Kapasouris, P.; Kappos, E.; Spang, H. A., III
1984-01-01
The design of a multivariable feedback control system for the Pratt and Whitney F-100 turbofan jet engine is a challenging task for control engineers. This paper employs a linearized model of the F-100 engine to demonstrate the use of the newly developed Linear Quadratic Gaussian/Loop Transfer Recovery (LQG/LTR) design methodology, which adopts an integrated frequency-domain and time-domain approach to multivariable feedback control synthesis so as to meet stability-robustness, command-following, and disturbance-rejection specifications.
NASA Astrophysics Data System (ADS)
Savić, Stevan; Milovanović, Boško; Lužanin, Zorana; Lazić, Lazar; Dolinaj, Dragan
2015-08-01
This paper presents an analysis of the homogenised mean maximum ( T max) and minimum ( T min) temperatures. The data used in the analysis were collected at eight stations in the Autonomous Province of Vojvodina (Serbia) during the 1949-2008 period. The trends obtained from the slopes of the regression lines using the least square method show 0.9 °C/60 years for T max and 1.1 °C/60 years for T min; the non-parametric Mann-Kendall test was used to determine the statistically significant increasing trends of these two extreme parameters. In this paper, we analyse the influence of the Vangengeim-Girs classification of atmospheric circulation on the T max and T min trends in the Autonomous Province of Vojvodina (Serbia) using linear and quadratic models based on the least square method. Linear stepwise regression and the forward method reveal the highest dependence of T max and T min when the W or E circulation types are included in the model. Non-linear models show a greater contribution of T max and T min at W, E and C circulation types, respectively. The correction of the variance contribution of quadratic models ranges from approximately 16 to 44 % for T max and 32 to 38 % for T min.
Haley, Benjamin M.; Paunesku, Tatjana; Grdina, David J.; Woloschak, Gayle E.
2015-01-01
Introduction The US government regulates allowable radiation exposures relying, in large part, on the seventh report from the committee to estimate the Biological Effect of Ionizing Radiation (BEIR VII), which estimated that most contemporary exposures- protracted or low-dose, carry 1.5 fold less risk of carcinogenesis and mortality per Gy than acute exposures of atomic bomb survivors. This correction is known as the dose and dose rate effectiveness factor for the life span study of atomic bomb survivors (DDREFLSS). It was calculated by applying a linear-quadratic dose response model to data from Japanese atomic bomb survivors and a limited number of animal studies. Methods and Results We argue that the linear-quadratic model does not provide appropriate support to estimate the risk of contemporary exposures. In this work, we re-estimated DDREFLSS using 15 animal studies that were not included in BEIR VII’s original analysis. Acute exposure data led to a DDREFLSS estimate from 0.9 to 3.0. By contrast, data that included both acute and protracted exposures led to a DDREFLSS estimate from 4.8 to infinity. These two estimates are significantly different, violating the assumptions of the linear-quadratic model, which predicts that DDREFLSS values calculated in either way should be the same. Conclusions Therefore, we propose that future estimates of the risk of protracted exposures should be based on direct comparisons of data from acute and protracted exposures, rather than from extrapolations from a linear-quadratic model. The risk of low dose exposures may be extrapolated from these protracted estimates, though we encourage ongoing debate as to whether this is the most valid approach. We also encourage efforts to enlarge the datasets used to estimate the risk of protracted exposures by including both human and animal data, carcinogenesis outcomes, a wider range of exposures, and by making more radiobiology data publicly accessible. We believe that these steps will
NASA Technical Reports Server (NTRS)
Gangsaas, D.; Ly, U.; Norman, D. C.
1981-01-01
A modified linear quadratic Gaussian (LQG) synthesis procedure has been used to design low-order robust multiloop controllers for a flexible airplane. The introduction of properly constructed fictitious Gauss-Markov processes in the control loops allowed meeting classical frequency-domain stability criteria using the direct synthesis procedures of modern time-domain control theory. Model reduction was used to simplify the control laws to the point where they could be easily implemented on onboard flight computers. These control laws provided excellent gust load and flutter mode control with good stability margins and compared very favorably to other control laws synthesized by the classical root-locus technique.
CAD of control systems: Application of nonlinear programming to a linear quadratic formulation
NASA Technical Reports Server (NTRS)
Fleming, P.
1983-01-01
The familiar suboptimal regulator design approach is recast as a constrained optimization problem and incorporated in a Computer Aided Design (CAD) package where both design objective and constraints are quadratic cost functions. This formulation permits the separate consideration of, for example, model following errors, sensitivity measures and control energy as objectives to be minimized or limits to be observed. Efficient techniques for computing the interrelated cost functions and their gradients are utilized in conjunction with a nonlinear programming algorithm. The effectiveness of the approach and the degree of insight into the problem which it affords is illustrated in a helicopter regulation design example.
Gain scheduling adaptive control strategies for HVDC systems to accommodate large disturbances
Reeve, J.; Sultan, M. )
1994-02-01
Techniques have been developed to permit the response of the controls for dc transmission systems to adapt to large system changes. A gain scheduling approach tunes the control as an on-line function of the effective short-circuit ratio and contingency indicators. The method has been tested by digital simulation, based on EMTP, of a back-to-back dc system. It has been found to be robust and control performance has been enhanced.
Guerrero, Mariana; Carlone, Marco
2010-08-15
Purpose: In recent years, several models were proposed that modify the standard linear-quadratic (LQ) model to make the predicted survival curve linear at high doses. Most of these models are purely phenomenological and can only be applied in the particular case of acute doses per fraction. The authors consider a mechanistic formulation of a linear-quadratic-linear (LQL) model in the case of split-dose experiments and exponentially decaying sources. This model provides a comprehensive description of radiation response for arbitrary dose rate and fractionation with only one additional parameter. Methods: The authors use a compartmental formulation of the LQL model from the literature. They analytically solve the model's differential equations for the case of a split-dose experiment and for an exponentially decaying source. They compare the solutions of the survival fraction with the standard LQ equations and with the lethal-potentially lethal (LPL) model. Results: In the case of the split-dose experiment, the LQL model predicts a recovery ratio as a function of dose per fraction that deviates from the square law of the standard LQ. The survival fraction as a function of time between fractions follows a similar exponential law as the LQ but adds a multiplicative factor to the LQ parameter {beta}. The LQL solution for the split-dose experiment is very close to the LPL prediction. For the decaying source, the differences between the LQL and the LQ solutions are negligible when the half-life of the source is much larger than the characteristic repair time, which is the clinically relevant case. Conclusions: The compartmental formulation of the LQL model can be used for arbitrary dose rates and provides a comprehensive description of dose response. When the survival fraction for acute doses is linear for high dose, a deviation of the square law formula of the recovery ratio for split doses is also predicted.
NASA Astrophysics Data System (ADS)
Yubai, Kazuhiro; Okuhara, Kazunori; Hirai, Junji
Gain-scheduling control is one of effective methods for plants whose dynamics changes significantly according to its operating point. A frozen parameter method is known to be a practical gain-scheduling controller synthesis, which interpolates the controllers designed at the prespecified (frozen) operating points according to the current operation point. Hyde et al. proposed a gain-scheduling control that H∞ loop shaping procedure is adopted as a controller synthesis at each operating point. H∞ loop shaping procedure is based on loop shaping of an open loop characteristic by frequency weights and is known to be effective for plants with bad condition number. However, weight selection satisfying control specifications is hard job for a designer. This paper describes the design of a suboptimal weight and a controller by means of algorithm that maximizes the robust stability margin and shapes the open loop characteristic into the desired shape at each operating point. Moreover, we formulate a weight optimization problem as a generalized eigenvalue minimization problem, which reduces the designer's burden of weight selection. Finally, we realize robust and high performance control system by scheduling both weights and controllers. The effectiveness of the proposed control system is verified in terms of the achieved robust stability margin and experimental time responses of a rotary inverted pendulum which involves strong nonlinear dynamics.
Mizuta, Masahiro; Takao, Seishin; Date, Hiroyuki; Kishimoto, Naoki; Sutherland, Kenneth L.; Onimaru, Rikiya; Shirato, Hiroki
2012-11-01
Purpose: Hypofractionated irradiation is often used in precise radiotherapy instead of conventional multifractionated irradiation. We propose a novel mathematical method for selecting a hypofractionated or multifractionated irradiation regimen based on physical dose distribution adding to biologic consideration. Methods and Materials: The linear-quadratic model was used for the radiation effects on tumor and normal tissues, especially organs at risk (OARs). On the basis of the assumption that the OAR receives a fraction of the dose intended for the tumor, the minimization problem for the damage effect on the OAR was treated under the constraint that the radiation effect on the tumor is fixed. Results: For an N-time fractionated irradiation regimen, the constraint of tumor lethality was described by an N-dimensional hypersphere. The total dose of the fractionated irradiations was considered for minimizing the damage effect on the OAR under the hypersphere condition. It was found that the advantage of hypofractionated or multifractionated irradiation therapies depends on the magnitude of the ratio of {alpha}/{beta} parameters for the OAR and tumor in the linear-quadratic model and the ratio of the dose for the OAR and tumor. Conclusions: Our mathematical method shows that multifractionated irradiation with a constant dose is better if the ratio of {alpha}/{beta} for the OAR and tumor is less than the ratio of the dose for the OAR and tumor, whereas hypofractionated irradiation is better otherwise.
NASA Technical Reports Server (NTRS)
Rosen, I. G.; Wang, C.
1990-01-01
The convergence of solutions to the discrete or sampled time linear quadratic regulator problem and associated Riccati equation for infinite dimensional systems to the solutions to the corresponding continuous time problem and equation, as the length of the sampling interval (the sampling rate) tends toward zero (infinity) is established. Both the finite and infinite time horizon problems are studied. In the finite time horizon case, strong continuity of the operators which define the control system and performance index together with a stability and consistency condition on the sampling scheme are required. For the infinite time horizon problem, in addition, the sampled systems must be stabilizable and detectable, uniformly with respect to the sampling rate. Classes of systems for which this condition can be verified are discussed. Results of numerical studies involving the control of a heat/diffusion equation, a hereditary of delay system, and a flexible beam are presented and discussed.
Haley, Benjamin M.; Paunesku, Tatjana; Grdina, David J.; Woloschak, Gayle E.; Aravindan, Natarajan
2015-12-09
The US government regulates allowable radiation exposures relying, in large part, on the seventh report from the committee to estimate the Biological Effect of Ionizing Radiation (BEIR VII), which estimated that most contemporary exposures- protracted or low-dose, carry 1.5 fold less risk of carcinogenesis and mortality per Gy than acute exposures of atomic bomb survivors. This correction is known as the dose and dose rate effectiveness factor for the life span study of atomic bomb survivors (DDREFLSS). As a result, it was calculated by applying a linear-quadratic dose response model to data from Japanese atomic bomb survivors and a limitedmore » number of animal studies.« less
Haley, Benjamin M.; Paunesku, Tatjana; Grdina, David J.; Woloschak, Gayle E.; Aravindan, Natarajan
2015-12-09
The US government regulates allowable radiation exposures relying, in large part, on the seventh report from the committee to estimate the Biological Effect of Ionizing Radiation (BEIR VII), which estimated that most contemporary exposures- protracted or low-dose, carry 1.5 fold less risk of carcinogenesis and mortality per Gy than acute exposures of atomic bomb survivors. This correction is known as the dose and dose rate effectiveness factor for the life span study of atomic bomb survivors (DDREF_{LSS}). As a result, it was calculated by applying a linear-quadratic dose response model to data from Japanese atomic bomb survivors and a limited number of animal studies.
NASA Technical Reports Server (NTRS)
Rosen, I. G.; Wang, C.
1992-01-01
The convergence of solutions to the discrete- or sampled-time linear quadratic regulator problem and associated Riccati equation for infinite-dimensional systems to the solutions to the corresponding continuous time problem and equation, as the length of the sampling interval (the sampling rate) tends toward zero(infinity) is established. Both the finite-and infinite-time horizon problems are studied. In the finite-time horizon case, strong continuity of the operators that define the control system and performance index, together with a stability and consistency condition on the sampling scheme are required. For the infinite-time horizon problem, in addition, the sampled systems must be stabilizable and detectable, uniformly with respect to the sampling rate. Classes of systems for which this condition can be verified are discussed. Results of numerical studies involving the control of a heat/diffusion equation, a hereditary or delay system, and a flexible beam are presented and discussed.
Using Linear versus Quadratic Rules in Predictive and Descriptive Discriminant Analysis.
ERIC Educational Resources Information Center
McGee, Jennifer
Both predictive discriminant analysis (PDA) and descriptive discriminant analysis (DDA) require a decision to pool group covariance matrices, or alternatively, to retain separate group covariance matrices when the group covariance matrices are too dissimilar to pool. Pooling the group covariance matrices involves the so-called "linear" rule,…
Lefkoff, L.J.; Gorelick, S.M.
1987-01-01
A FORTRAN-77 computer program code that helps solve a variety of aquifer management problems involving the control of groundwater hydraulics. It is intended for use with any standard mathematical programming package that uses Mathematical Programming System input format. The computer program creates the input files to be used by the optimization program. These files contain all the hydrologic information and management objectives needed to solve the management problem. Used in conjunction with a mathematical programming code, the computer program identifies the pumping or recharge strategy that achieves a user 's management objective while maintaining groundwater hydraulic conditions within desired limits. The objective may be linear or quadratic, and may involve the minimization of pumping and recharge rates or of variable pumping costs. The problem may contain constraints on groundwater heads, gradients, and velocities for a complex, transient hydrologic system. Linear superposition of solutions to the transient, two-dimensional groundwater flow equation is used by the computer program in conjunction with the response matrix optimization method. A unit stress is applied at each decision well and transient responses at all control locations are computed using a modified version of the U.S. Geological Survey two dimensional aquifer simulation model. The program also computes discounted cost coefficients for the objective function and accounts for transient aquifer conditions. (Author 's abstract)
Ayvaz, Muzaffer; Demiralp, Metin
2011-09-14
In this study, the optimal control equations for one dimensional quantum harmonic oscillator under the quadratic control operators together with linear dipole polarizability effects are constructed in the sense of Heisenberg equation of motion. A numerical technique based on the approximation to the non-commuting quantum mechanical operators from the fluctuation free expectation value dynamics perspective in the classical limit is also proposed for the solution of optimal control equations which are ODEs with accompanying boundary conditions. The dipole interaction of the system is considered to be linear, and the observable whose expectation value will be suppressed during the control process is considered to be quadratic in terms of position operator x. The objective term operator is also assumed to be quadratic.
Kohli, Nidhi; Sullivan, Amanda L; Sadeh, Shanna; Zopluoglu, Cengiz
2015-04-01
Effective instructional planning and intervening rely heavily on accurate understanding of students' growth, but relatively few researchers have examined mathematics achievement trajectories, particularly for students with special needs. We applied linear, quadratic, and piecewise linear mixed-effects models to identify the best-fitting model for mathematics development over elementary and middle school and to ascertain differences in growth trajectories of children with learning disabilities relative to their typically developing peers. The analytic sample of 2150 students was drawn from the Early Childhood Longitudinal Study - Kindergarten Cohort, a nationally representative sample of United States children who entered kindergarten in 1998. We first modeled students' mathematics growth via multiple mixed-effects models to determine the best fitting model of 9-year growth and then compared the trajectories of students with and without learning disabilities. Results indicate that the piecewise linear mixed-effects model captured best the functional form of students' mathematics trajectories. In addition, there were substantial achievement gaps between students with learning disabilities and students with no disabilities, and their trajectories differed such that students without disabilities progressed at a higher rate than their peers who had learning disabilities. The results underscore the need for further research to understand how to appropriately model students' mathematics trajectories and the need for attention to mathematics achievement gaps in policy. PMID:25746821
Post-Stall Aerodynamic Modeling and Gain-Scheduled Control Design
NASA Technical Reports Server (NTRS)
Wu, Fen; Gopalarathnam, Ashok; Kim, Sungwan
2005-01-01
A multidisciplinary research e.ort that combines aerodynamic modeling and gain-scheduled control design for aircraft flight at post-stall conditions is described. The aerodynamic modeling uses a decambering approach for rapid prediction of post-stall aerodynamic characteristics of multiple-wing con.gurations using known section data. The approach is successful in bringing to light multiple solutions at post-stall angles of attack right during the iteration process. The predictions agree fairly well with experimental results from wind tunnel tests. The control research was focused on actuator saturation and .ight transition between low and high angles of attack regions for near- and post-stall aircraft using advanced LPV control techniques. The new control approaches maintain adequate control capability to handle high angle of attack aircraft control with stability and performance guarantee.
NASA Astrophysics Data System (ADS)
Na, Xiaoxiang; Cole, David J.
2013-02-01
This paper is concerned with the modelling of strategic interactions between the human driver and the vehicle active front steering (AFS) controller in a path-following task where the two controllers hold different target paths. The work is aimed at extending the use of mathematical models in representing driver steering behaviour in complicated driving situations. Two game theoretic approaches, namely linear quadratic game and non-cooperative model predictive control (non-cooperative MPC), are used for developing the driver-AFS interactive steering control model. For each approach, the open-loop Nash steering control solution is derived; the influences of the path-following weights, preview and control horizons, driver time delay and arm neuromuscular system (NMS) dynamics are investigated, and the CPU time consumed is recorded. It is found that the two approaches give identical time histories as well as control gains, while the non-cooperative MPC method uses much less CPU time. Specifically, it is observed that the introduction of weight on the integral of vehicle lateral displacement error helps to eliminate the steady-state path-following error; the increase in preview horizon and NMS natural frequency and the decline in time delay and NMS damping ratio improve the path-following accuracy.
NASA Astrophysics Data System (ADS)
Bromley, Regina; Oliver, Lyn; Davey, Ross; Harvie, Rozelle; Baldock, Clive
2009-01-01
In this study we present two prediction methods, mean dose and summed dose, for predicting the number of A549 cells that will survive after modulated x-ray irradiation. The prediction methods incorporate the dose profile from the modulated x-ray fluence map applied across the cell sample and the linear quadratic (LQ) model. We investigated the clonogenic survival of A549 cells when irradiated using two different modulated x-ray fluence maps. Differences between the measured and predicted surviving fraction were observed for modulated x-ray irradiation. When the x-ray fluence map produced a steep dose gradient across the sample, fewer cells survived in the unirradiated region than expected. When the x-ray fluence map produced a less steep dose gradient across the sample, more cells survived in the unirradiated region than expected. Regardless of the steepness of the dose gradient, more cells survived in the irradiated region than expected for the reference dose range of 1-10 Gy. The change in the cell survival for the unirradiated regions of the two different dose gradients may be an important factor to consider when predicting the number of cells that will survive at the edge of modulated x-ray fields. This investigation provides an improved method of predicting cell survival for modulated x-ray radiation treatment. It highlights the limitations of the LQ model, particularly in its ability to describe the biological response of cells irradiated under these conditions.
NASA Astrophysics Data System (ADS)
McAneney, H.; O'Rourke, S. F. C.
2007-02-01
The standard linear-quadratic survival model for radiotherapy is used to investigate different schedules of radiation treatment planning to study how these may be affected by different tumour repopulation kinetics between treatments. The laws for tumour cell repopulation include the logistic and Gompertz models and this extends the work of Wheldon et al (1977 Br. J. Radiol. 50 681), which was concerned with the case of exponential re-growth between treatments. Here we also consider the restricted exponential model. This has been successfully used by Panetta and Adam (1995 Math. Comput. Modelling 22 67) in the case of chemotherapy treatment planning.Treatment schedules investigated include standard fractionation of daily treatments, weekday treatments, accelerated fractionation, optimized uniform schedules and variation of the dosage and α/β ratio, where α and β are radiobiological parameters for the tumour tissue concerned. Parameters for these treatment strategies are extracted from the literature on advanced head and neck cancer, prostate cancer, as well as radiosensitive parameters. Standardized treatment protocols are also considered. Calculations based on the present analysis indicate that even with growth laws scaled to mimic initial growth, such that growth mechanisms are comparable, variation in survival fraction to orders of magnitude emerged. Calculations show that the logistic and exponential models yield similar results in tumour eradication. By comparison the Gompertz model calculations indicate that tumours described by this law result in a significantly poorer prognosis for tumour eradication than either the exponential or logistic models. The present study also shows that the faster the tumour growth rate and the higher the repair capacity of the cell line, the greater the variation in outcome of the survival fraction. Gaps in treatment, planned or unplanned, also accentuate the differences of the survival fraction given alternative growth
Habibullah, H. Pota, H. R. Petersen, I. R.
2014-03-15
This paper demonstrates a high-speed spiral imaging technique for an atomic force microscope (AFM). As an alternative to traditional raster scanning, an approach of gradient pulsing using a spiral line is implemented and spirals are generated by applying single-frequency cosine and sine waves of slowly varying amplitudes to the X and Y-axes of the AFM’s piezoelectric tube scanner (PTS). Due to these single-frequency sinusoidal input signals, the scanning process can be faster than that of conventional raster scanning. A linear quadratic Gaussian controller is designed to track the reference sinusoid and a vibration compensator is combined to damp the resonant mode of the PTS. An internal model of the reference sinusoidal signal is included in the plant model and an integrator for the system error is introduced in the proposed control scheme. As a result, the phase error between the input and output sinusoids from the X and Y-PTSs is reduced. The spirals produced have particularly narrow-band frequency measures which change slowly over time, thereby making it possible for the scanner to achieve improved tracking and continuous high-speed scanning rather than being restricted to the back and forth motion of raster scanning. As part of the post-processing of the experimental data, a fifth-order Butterworth filter is used to filter noises in the signals emanating from the position sensors and a Gaussian image filter is used to filter the images. A comparison of images scanned using the proposed controller (spiral) and the AFM PI controller (raster) shows improvement in the scanning rate using the proposed method.
Control of microbial fuel cell voltage using a gain scheduling control strategy
NASA Astrophysics Data System (ADS)
Boghani, Hitesh C.; Michie, Iain; Dinsdale, Richard M.; Guwy, Alan J.; Premier, Giuliano C.
2016-08-01
Recent microbial fuel cell (MFC) research frequently addresses matters associated with scale and deployability. Modularisation is often needed to reduce ohmic losses with increasing volume. Series/parallel is then often an obvious strategy to enhance power quality during operation, to make best use of generated electricity. Hence, voltage reversal resulting from power and voltage mismatch between cells become virtually unavoidable. Control MFC voltages could be used to stabilise MFC stacks. Here, nonlinear MFCs are controlled using simple gain scheduled Proportional + Integral actions. Parsimonious control may be necessary for implementation in MFC arrays, so minimising costs. Controller parameterisation used several linearised models over the dynamic operating range of the MFCs. Controller gains were then scheduled according to the operating conditions. A digital potentiometer was used to actuate the control, varying the current sourced from the MFC. The results show that the controller was able to control MFC voltages, rejecting the disturbances. It was shown that the controller was transferable between MFCs with different power performances. This study demonstrates that the control of MFCs can be achieved with relatively simple digital approaches, plausibly implementable using low cost microcontrollers, and likely to be useful in the effective deployment of MFCs in large scale arrays.
Feedback shape control for deployable mesh reflectors using gain scheduling method
NASA Astrophysics Data System (ADS)
Xie, Yangmin; Shi, Hang; Alleyne, Andrew; Yang, Bingen
2016-04-01
This paper presents a theoretical study on the dynamic shape control problem of deployable mesh reflectors (DMRs) via feedback approaches. The reflector structure is simplified from a nonlinear model to be quasi-static with respect to temperature variations but dynamic with respect to mechanical vibrations. The orbital cycle is segmented into multiple temperature zones, and an H∞ robust state feedback controller is designed for each zone to guarantee the local stability of the system under the model uncertainty caused by thermal effects and to reject external force disturbances. At the same time, gain scheduling control method is adopted to compensate thermal distortions and to ensure smooth transition response when switching among the local robust controllers. A DMR model is considered in the case study to show the effectiveness of the control approach. The structural vibrations caused by external force disturbances can be sufficiently suppressed in a much shorter time. The closed loop response of the DMR structure shows that much higher surface accuracy is obtained during the orbiting mission compared to the open-loop configuration, and transient focal length and transient de-focus of the reflector are well controlled within the satisfactory bounds, demonstrating the numerical feasibility of the proposed method to solve the dynamic shape control problem of DMRs.
ERIC Educational Resources Information Center
Fay, Temple H.
2012-01-01
Quadratic friction involves a discontinuous damping term in equations of motion in order that the frictional force always opposes the direction of the motion. Perhaps for this reason this topic is usually omitted from beginning texts in differential equations and physics. However, quadratic damping is more realistic than viscous damping in many…
NASA Astrophysics Data System (ADS)
Hao, Lili; Wang, Qiang; Hou, Chunfeng
2016-05-01
The temperature effects on the intensity profile, self-deflection process, and stability of (1 + 1)-dimensional steady-state bright solitons resulting from both the linear and quadratic electro-optic effects are comprehensively analyzed. Moreover, three physical factors, i.e. diffusion effect, dark irradiance, and the dielectric constant, have been investigated through the theoretical analysis to determine which one dominates the temperature dependence of intensity profile and self-deflection of bright solitons. It is also found that the incident beam evolves into stable bright solitons in the vicinity of initial temperature, but oscillates or even collapses when the crystal temperature deviates significantly from the initial temperature.
Liu, Shichao; Liu, Xiaoping P; El Saddik, Abdulmotaleb
2014-03-01
In this paper, we investigate the modeling and distributed control problems for the load frequency control (LFC) in a smart grid. In contrast with existing works, we consider more practical and real scenarios, where the communication topology of the smart grid changes because of either link failures or packet losses. These topology changes are modeled as a time-varying communication topology matrix. By using this matrix, a new closed-loop power system model is proposed to integrate the communication topology changes into the dynamics of a physical power system. The globally asymptotical stability of this closed-loop power system is analyzed. A distributed gain scheduling LFC strategy is proposed to compensate for the potential degradation of dynamic performance (mean square errors of state vectors) of the power system under communication topology changes. In comparison to conventional centralized control approaches, the proposed method can improve the robustness of the smart grid to the variation of the communication network as well as to reduce computation load. Simulation results show that the proposed distributed gain scheduling approach is capable to improve the robustness of the smart grid to communication topology changes. PMID:24200162
Vehicle rollover avoidance by application of gain-scheduled LQR controllers using state observers
NASA Astrophysics Data System (ADS)
Dal Poggetto, Vinicius F.; Serpa, Alberto L.
2016-02-01
Many researches have been conducted in the area of control applied to vehicle dynamics, aiming at reducing the possibility of the occurrence of the type of accident known as rollover. In this research, based on a common nonlinear model and its linearisation, a method for properly selecting matrices for solving the Riccati equation considering different speeds was proposed. The method showed in which ways speed really influences the choice of controller gains. By developing the dynamic equations for the yaw- and roll-coupled motions and modelling of controllers and state observers, it is possible to compare the efficacy of this control strategy using both linear and nonlinear simulations using Matlab. Significant results were obtained regarding the reduction of the rollover coefficient for a double-lane change manoeuvre at different speeds, thus indicating advantages of using this controller in practical cases.
NASA Technical Reports Server (NTRS)
Hamer, H. A.; Johnson, K. G.
1986-01-01
An analysis was performed to determine the effects of model error on the control of a large flexible space antenna. Control was achieved by employing two three-axis control-moment gyros (CMG's) located on the antenna column. State variables were estimated by including an observer in the control loop that used attitude and attitude-rate sensors on the column. Errors were assumed to exist in the individual model parameters: modal frequency, modal damping, mode slope (control-influence coefficients), and moment of inertia. Their effects on control-system performance were analyzed either for (1) nulling initial disturbances in the rigid-body modes, or (2) nulling initial disturbances in the first three flexible modes. The study includes the effects on stability, time to null, and control requirements (defined as maximum torque and total momentum), as well as on the accuracy of obtaining initial estimates of the disturbances. The effects on the transients of the undisturbed modes are also included. The results, which are compared for decoupled and linear quadratic regulator (LQR) control procedures, are shown in tabular form, parametric plots, and as sample time histories of modal-amplitude and control responses. Results of the analysis showed that the effects of model errors on the control-system performance were generally comparable for both control procedures. The effect of mode-slope error was the most serious of all model errors.
Meigooni, Ali S.; Dou Kai; Meigooni, Navid J.; Gnaster, Michael; Awan, Shahid; Dini, Sharifeh; Johnson, Ellis L.
2006-09-15
Grid radiation therapy with megavoltage x-ray beam has been proven to be an effective technique for management of large, bulky malignant tumors. The clinical advantage of GRID therapy, combined with conventional radiation therapy, has been demonstrated using a prototype GRID block [Mohiuddin, Curtis, Grizos, and Komarnicky, Cancer 66, 114-118 (1990)]. Recently, a new GRID block design with improved dosimetric properties has become commercially available from Radiation Product Design, Inc. (Albertive, MN). This GRID collimator consists of an array of focused apertures in a cerrobend block arranged in a hexagonal pattern having a circular cross-section with a diameter and center-to-center spacing of 14.3 and 21.1 mm, respectively, in the plane of isocenter. In this project, dosimetric characteristics of the newly redesigned GRID block have been investigated for a Varian 21EX linear accelerator (Varian Associates, Palo Alto, CA). These determinations were performed using radiographic films, thermoluminescent dosimeters in Solid Water trade mark sign phantom materials, and an ionization chamber in water. The output factor, percentage depth dose, beam profiles, and isodose distributions of the GRID radiation as a function of field size and beam energy have been measured using both 6 and 18 MV x-ray beams. In addition, the therapeutic advantage obtained from this treatment modality with the new GRID block design for a high, single fraction of dose has been calculated using the linear quadratic model with {alpha}/{beta} ratios for typical tumor and normal cells. These biological characteristics of the new GRID block design will also be presented.
Meigooni, Ali S; Dou, Kai; Meigooni, Navid J; Gnaster, Michael; Awan, Shahid; Dini, Sharifeh; Johnson, Ellis L
2006-09-01
Grid radiation therapy with megavoltage x-ray beam has been proven to be an effective technique for management of large, bulky malignant tumors. The clinical advantage of GRID therapy, combined with conventional radiation therapy, has been demonstrated using a prototype GRID block [Mohiuddin, Curtis, Grizos, and Komarnicky, Cancer 66, 114-118 (1990)]. Recently, a new GRID block design with improved dosimetric properties has become commercially available from Radiation Product Design, Inc. (Albertive, MN). This GRID collimator consists of an array of focused apertures in a cerrobend block arranged in a hexagonal pattern having a circular cross-section with a diameter and center-to-center spacing of 14.3 and 21.1 mm, respectively, in the plane of isocenter. In this project, dosimetric characteristics of the newly redesigned GRID block have been investigated for a Varian 21EX linear accelerator (Varian Associates, Palo Alto, CA). These determinations were performed using radiographic films, thermoluminescent dosimeters in Solid Water phantom materials, and an ionization chamber in water. The output factor, percentage depth dose, beam profiles, and isodose distributions of the GRID radiation as a function of field size and beam energy have been measured using both 6 and 18 MV x-ray beams. In addition, the therapeutic advantage obtained from this treatment modality with the new GRID block design for a high, single fraction of dose has been calculated using the linear quadratic model with alpha/beta ratios for typical tumor and normal cells. These biological characteristics of the new GRID block design will also be presented. PMID:17022209
A model-based gain scheduling approach for controlling the common-rail system for GDI engines
NASA Astrophysics Data System (ADS)
di Gaeta, Alessandro; Montanaro, Umberto; Fiengo, Giovanni; Palladino, Angelo; Giglio, Veniero
2012-04-01
The progressive reduction in vehicle emission requirements have forced the automotive industry to invest in research for developing alternative and more efficient control strategies. All control features and resources are permanently active in an electronic control unit (ECU), ensuring the best performance with respect to emissions, fuel economy, driveability and diagnostics, independently from engine working point. In this article, a considerable step forward has been achieved by the common-rail technology which has made possible to vary the injection pressure over the entire engine speed range. As a consequence, the injection of a fixed amount of fuel is more precise and multiple injections in a combustion cycle can be made. In this article, a novel gain scheduling pressure controller for gasoline direct injection (GDI) engine is designed to stabilise the mean fuel pressure into the rail and to track demanded pressure trajectories. By exploiting a simple control-oriented model describing the mean pressure dynamics in the rail, the control structure turns to be simple enough to be effectively implemented in commercial ECUs. Experimental results in a wide range of operating points confirm the effectiveness of the proposed control method to tame efficiently the mean value pressure dynamics of the plant showing a good accuracy and robustness with respect to unavoidable parameters uncertainties, unmodelled dynamics, and hidden coupling terms.
Nikzad, Safoora; Hashemi, Bijan; Hasan, Zuhair Saraf; Mozdarani, Hossein; Baradaran-Ghahfarokhi, Milad; Amini, Payam
2016-01-01
Purpose To investigate the effect of increasing the overall treatment time as well as delivering the compensating doses on the Balb/c breast adenocarcinoma (4T1) tumor. Materials and methods A total of 72 mice were divided into two aliquots (classes A and B) based on the initial size of their induced tumor. Each class was divided into a control and several treatment groups. Among the treatment groups, group 1 was continuously exposed to 2 Gy irradiation, and groups 2 and 3 received two subfractions of 1 Gy over the total treatment times of 30 and 60 min, respectively. To investigate the effect of compensating doses, calculated based on the developed linear quadratic model (LQ) model, the remaining two groups (groups 4 and 5) received two subfractions of 1.16 and 1.24 Gy over the total treatment times of 30 and 60 min, respectively. The growing curves, Tumor Growth Time (TGT), Tumor Growth Delay Time (TGDT) and the survival of the animals were studied. Results For class A (tumor size ≤ 30 mm(3)), the average tumor size in the irradiated groups 1-5 was considerably different compared to the control group as one unit (day) change in time, by amount of -160.8, -158.9, +39.4 and +44.0, respectively. While these amounts were +22.0, +17.9, -21.7 and -0.1 for class B (tumor size ≥ 400 mm(3)). For the class A of animals, the TGT and TGDT parameters were significantly lower (0 ≤ 0.05) for the groups 2 and 3, compared to group 1. There was no significant difference (p > 0.05) between groups 1, 4 and 5 in this class. There was no significant difference (p > 0.05) between all the treated groups in class B. Conclusions Increasing total treatment time affects the radiobiological efficiency of treatment especially in small-sized tumor. The compensating doses derived from the LQ model can be used to compensate the effects of prolonged treatment times at in vivo condition. PMID:26630280
ERIC Educational Resources Information Center
Withers, Christopher S.; Nadarajah, Saralees
2012-01-01
We show that there are exactly four quadratic polynomials, Q(x) = x [superscript 2] + ax + b, such that (x[superscript 2] + ax + b) (x[superscript 2] - ax + b) = (x[superscript 4] + ax[superscript 2] + b). For n = 1, 2, ..., these quadratic polynomials can be written as the product of N = 2[superscript n] quadratic polynomials in x[superscript…
NASA Astrophysics Data System (ADS)
Sandhu, Amit
A sequential quadratic programming method is proposed for solving nonlinear optimal control problems subject to general path constraints including mixed state-control and state only constraints. The proposed algorithm further develops on the approach proposed in [1] with objective to eliminate the use of a high number of time intervals for arriving at an optimal solution. This is done by introducing an adaptive time discretization to allow formation of a desirable control profile without utilizing a lot of intervals. The use of fewer time intervals reduces the computation time considerably. This algorithm is further used in this thesis to solve a trajectory planning problem for higher elevation Mars landing.
Identification and robust control of linear parameter-varying systems
NASA Astrophysics Data System (ADS)
Lee, Lawton Hubert
This dissertation deals with linear parameter-varying (LPV) systems: linear dynamic systems that depend on time-varying parameters. These systems appear in gain scheduling problems, and much recent research has been devoted to their prospective usefulness for systematic gain scheduling. We primarily focus on robust control of uncertain LPV systems and identification of LPV systems that are modelable as linear-fractional transformations (LFTs). Using parameter-dependent quadratic Lyapunov functions, linear matrix inequalities (LMIs), and scaled small-gain arguments, we define notions of stability and induced-{cal L}sb2 performance for uncertain LPV systems whose parameters and rates of parameter variation satisfy given bounds. The performance criterion involves integral quadratic constraints and implies naturally parameter-dependent induced-{cal L}sb2 norm bounds. We formulate and solve an {cal H}sb{infty}-like control problem for an LPV plant with measurable parameters and an "Output/State Feedback" structure: the feedback outputs include some noiselessly measured states. Necessary and sufficient solvability conditions reduce to LMIs that can be solved approximately using finite-dimensional convex programming. Reduced-order LPV controllers are constructed from the LMI solutions. A D-K iteration-like procedure provides robustness to structured, time-varying, parametric uncertainty. The design method is applied to a motivating example: flight control for the F-16 VISTA throughout its subsonic flight envelope. Parameter-dependent weights and {cal H}sb{infty} design principles describe the performance objectives. Closed-loop responses exhibited by nonlinear simulations indicate satisfactory flying qualities. Identification of linear-fractional LPV systems is treated using maximum-likelihood parameter estimation. Computing the gradient and Hessian of a maximum-likelihood cost function reduces to simulating one LPV filter per identified parameter. We use nonlinear
Quadratic boundedness of uncertain nonlinear dynamic systems
NASA Astrophysics Data System (ADS)
Brockman, Mark Lawrence
Physical systems are often perturbed by unknown external disturbances or contain important system parameters which are difficult to model exactly. However, engineers are expected to design systems which perform well even in the presence of uncertainties. For example, an airplane designer can never know the precise direction or magnitude of wind gusts, or the exact mass distribution inside the aircraft, but passengers expect to arrive on time after a smooth ride. This thesis will first present the concept of quadratic boundedness of an uncertain nonlinear dynamic system, and then develop analysis techniques and control design methods for systems containing unknown disturbances and parameters. For a class of nonlinear systems, conditions for quadratic boundedness are given, and the relationship between quadratic boundedness and quadratic stability is explored. An important consequence of quadratic boundedness is the ability to calculate an upper bound on the system gain of an uncertain nonlinear system. For nominally linear systems, necessary and sufficient conditions for quadratic boundedness are given. The innovative use of linear matrix inequalities in an iterative algorithm provides a means to analyze the quadratic boundedness properties of systems containing parameter uncertainties. The analysis results establish a framework for the development of design methods which integrate performance specifications into the control design process for all the types of systems considered. Numerous examples illustrate the major results of the thesis.
Quadratic eigenvalue problems.
Walsh, Timothy Francis; Day, David Minot
2007-04-01
In this report we will describe some nonlinear eigenvalue problems that arise in the areas of solid mechanics, acoustics, and coupled structural acoustics. We will focus mostly on quadratic eigenvalue problems, which are a special case of nonlinear eigenvalue problems. Algorithms for solving the quadratic eigenvalue problem will be presented, along with some example calculations.
NASA Astrophysics Data System (ADS)
Withers, Christopher S.; Nadarajah, Saralees
2012-06-01
We show that there are exactly four quadratic polynomials, Q(x) = x 2 + ax + b, such that
NASA Astrophysics Data System (ADS)
Zhou, J. X.; Zhang, L.
2005-01-01
Incremental harmonic balance (IHB) formulations are derived for general multiple degrees of freedom (d.o.f.) non-linear autonomous systems. These formulations are developed for a concerned four-d.o.f. aircraft wheel shimmy system with combined Coulomb and velocity-squared damping. A multi-harmonic analysis is performed and amplitudes of limit cycles are predicted. Within a large range of parametric variations with respect to aircraft taxi velocity, the IHB method can, at a much cheaper cost, give results with high accuracy as compared with numerical results given by a parametric continuation method. In particular, the IHB method avoids the stiff problems emanating from numerical treatment of aircraft wheel shimmy system equations. The development is applicable to other vibration control systems that include commonly used dry friction devices or velocity-squared hydraulic dampers.
The Mystical "Quadratic Formula."
ERIC Educational Resources Information Center
March, Robert H.
1993-01-01
Uses projectile motion to explain the two roots found when using the quadratic formula. An example is provided for finding the time of flight for a projectile which has a negative root implying a negative time of flight. This negative time of flight also has a useful physical meaning. (MVL)
ERIC Educational Resources Information Center
Fay, Temple H.
2010-01-01
Through numerical investigations, we study examples of the forced quadratic spring equation [image omitted]. By performing trial-and-error numerical experiments, we demonstrate the existence of stability boundaries in the phase plane indicating initial conditions yielding bounded solutions, investigate the resonance boundary in the [omega]…
Quadratic elongation: A quantitative measure of distortion in coordination polyhedra
Robinson, Kelly F.; Gibbs, G.V.; Ribbe, P.H.
1971-01-01
Quadratic elongation and the variance of bond angles are linearly correlated for distorted octahedral and tetrahedral coordination complexes, both of which show variations in bond length and bond angle. The quadratic elonga tion is dimensionless, giving a quantitative measure of polyhedral distortion which is independent of the effective size of the polyhedron.
Galactic chemical evolution and nucleocosmochronology - Analytic quadratic models
NASA Technical Reports Server (NTRS)
Clayton, D. D.
1985-01-01
Quadratic models of the chemical evolution of the Galaxy for a star formation rate proportional to the square of the gas mass are studied. The search for analytic solutions to the gas mass and star mass for time-dependent rates of gaseous infall onto the disk is examined. The quadratic models are compared to models having linear star formation rates. The mass, metallicity, number of stars, and U-235/U-238 isotopic ratio for the models which are subjected to the same infall rate, the same initial disk mass, and the same final gas fraction are compared. The results of the comparison indicate that: (1) the average dwarf age is greater in the quadratic model, (2) the metallicity grows initially faster in the quadratic model, (3) the quadratic model has a smaller percentage of low-Z dwarfs, and (4) the U-235/U-238 isotopic ratio indicates a younger quadratic model.
Fourier analysis of quadratic phase interferograms
NASA Astrophysics Data System (ADS)
Muñoz-Maciel, Jesús; Mora-González, Miguel; Casillas-Rodríguez, Francisco J.; Peña-Lecona, Francisco G.
2015-06-01
A phase demodulation method from a single interferogram with a quadratic phase term is developed. The fringe pattern being analysed may contain circular, elliptic or astigmatic fringes. The Fourier transform of such interferograms is seen to be also a sine or a cosine of a second order polynomial in both the real and imaginary parts. In this work we take a discrete Fourier transform of the fringe patterns and then we take separate inverse discrete transforms of the real and imaginary parts of the frequency spectrum. This results in two new interferograms corresponding to the sine and cosine of the quadratic term of the phase modulated by the sine and cosine of the linear term. The linear term of these interferograms may be recovered with similar procedures of fringe analysis from open fringe interferograms. Once the linear term is retrieved the quadratic phase of the interferogram being analysed can also be calculated. The present approach is also being investigated for interferograms with nearly circularly symmetry given that the phase contains some tilt. The described procedure of Fourier analysis from quadratic phase interferograms of nearly symmetric interferograms could be used instead of complex and time consuming algorithms for phase recovery from fringe patterns with closed fringes. Finally, the method is tested in simulated and real data.
Single-photon quadratic optomechanics
Liao, Jie-Qiao; Nori, Franco
2014-01-01
We present exact analytical solutions to study the coherent interaction between a single photon and the mechanical motion of a membrane in quadratic optomechanics. We consider single-photon emission and scattering when the photon is initially inside the cavity and in the fields outside the cavity, respectively. Using our solutions, we calculate the single-photon emission and scattering spectra, and find relations between the spectral features and the system's inherent parameters, such as: the optomechanical coupling strength, the mechanical frequency, and the cavity-field decay rate. In particular, we clarify the conditions for the phonon sidebands to be visible. We also study the photon-phonon entanglement for the long-time emission and scattering states. The linear entropy is employed to characterize this entanglement by treating it as a bipartite one between a single mode of phonons and a single photon. PMID:25200128
AdS waves as exact solutions to quadratic gravity
Guellue, Ibrahim; Sisman, Tahsin Cagri; Tekin, Bayram; Guerses, Metin
2011-04-15
We give an exact solution of the quadratic gravity in D dimensions. The solution is a plane-fronted wave metric with a cosmological constant. This metric solves not only the full quadratic gravity field equations but also the linearized ones which include the linearized equations of the recently found critical gravity. A subset of the solutions change the asymptotic structure of the anti-de Sitter space due to their logarithmic behavior.
Optimal Approximation of Quadratic Interval Functions
NASA Technical Reports Server (NTRS)
Koshelev, Misha; Taillibert, Patrick
1997-01-01
Measurements are never absolutely accurate, as a result, after each measurement, we do not get the exact value of the measured quantity; at best, we get an interval of its possible values, For dynamically changing quantities x, the additional problem is that we cannot measure them continuously; we can only measure them at certain discrete moments of time t(sub 1), t(sub 2), ... If we know that the value x(t(sub j)) at a moment t(sub j) of the last measurement was in the interval [x-(t(sub j)), x + (t(sub j))], and if we know the upper bound D on the rate with which x changes, then, for any given moment of time t, we can conclude that x(t) belongs to the interval [x-(t(sub j)) - D (t - t(sub j)), x + (t(sub j)) + D (t - t(sub j))]. This interval changes linearly with time, an is, therefore, called a linear interval function. When we process these intervals, we get an expression that is quadratic and higher order w.r.t. time t, Such "quadratic" intervals are difficult to process and therefore, it is necessary to approximate them by linear ones. In this paper, we describe an algorithm that gives the optimal approximation of quadratic interval functions by linear ones.
NASA Astrophysics Data System (ADS)
Colin, M.; Di Menza, L.; Saut, J. C.
2016-03-01
In this paper, we investigate the properties of solitonic structures arising in quadratic media. First, we recall the derivation of systems governing the interaction process for waves propagating in such media and we check the local and global well-posedness of the corresponding Cauchy problem. Then, we look for stationary states in the context of normal or anomalous dispersion regimes, that lead us to either elliptic or non-elliptic systems and we address the problem of orbital stability. Finally, some numerical experiments are carried out in order to compute localized states for several regimes and to study dynamic stability as well as long-time asymptotics.
Quadratic algebras for three-dimensional superintegrable systems
Daskaloyannis, C. Tanoudis, Y.
2010-02-15
The three-dimensional superintegrable systems with quadratic integrals of motion have five functionally independent integrals, one among them is the Hamiltonian. Kalnins, Kress, and Miller have proved that in the case of nondegenerate potentials with quadratic integrals of motion there is a sixth quadratic integral, which is linearly independent of the other integrals. The existence of this sixth integral implies that the integrals of motion form a ternary parafermionic-like quadratic Poisson algebra with five generators. In this contribution we investigate the structure of this algebra. We show that in all the nondegenerate cases there is at least one subalgebra of three integrals having a Poisson quadratic algebra structure, which is similar to the two-dimensional case.
Quadratic spatial soliton interactions
NASA Astrophysics Data System (ADS)
Jankovic, Ladislav
Quadratic spatial soliton interactions were investigated in this Dissertation. The first part deals with characterizing the principal features of multi-soliton generation and soliton self-reflection. The second deals with two beam processes leading to soliton interactions and collisions. These subjects were investigated both theoretically and experimentally. The experiments were performed by using potassium niobate (KNBO 3) and periodically poled potassium titanyl phosphate (KTP) crystals. These particular crystals were desirable for these experiments because of their large nonlinear coefficients and, more importantly, because the experiments could be performed under non-critical-phase-matching (NCPM) conditions. The single soliton generation measurements, performed on KNBO3 by launching the fundamental component only, showed a broad angular acceptance bandwidth which was important for the soliton collisions performed later. Furthermore, at high input intensities multi-soliton generation was observed for the first time. The influence on the multi-soliton patterns generated of the input intensity and beam symmetry was investigated. The combined experimental and theoretical efforts indicated that spatial and temporal noise on the input laser beam induced multi-soliton patterns. Another research direction pursued was intensity dependent soliton routing by using of a specially engineered quadratically nonlinear interface within a periodically poled KTP sample. This was the first time demonstration of the self-reflection phenomenon in a system with a quadratic nonlinearity. The feature investigated is believed to have a great potential for soliton routing and manipulation by engineered structures. A detailed investigation was conducted on two soliton interaction and collision processes. Birth of an additional soliton resulting from a two soliton collision was observed and characterized for the special case of a non-planar geometry. A small amount of spiraling, up to 30
Quadratic soliton self-reflection at a quadratically nonlinear interface
NASA Astrophysics Data System (ADS)
Jankovic, Ladislav; Kim, Hongki; Stegeman, George; Carrasco, Silvia; Torner, Lluis; Katz, Mordechai
2003-11-01
The reflection of bulk quadratic solutions incident onto a quadratically nonlinear interface in periodically poled potassium titanyl phosphate was observed. The interface consisted of the boundary between two quasi-phase-matched regions displaced from each other by a half-period. At high intensities and small angles of incidence the soliton is reflected.
Degenerate nonlinear programming with a quadratic growth condition.
Anitescu, M.; Mathematics and Computer Science
2000-01-01
We show that the quadratic growth condition and the Mangasarian-Fromovitz constraint qualification (MFCQ) imply that local minima of nonlinear programs are isolated stationary points. As a result, when started sufficiently close to such points, an L1 exact penalty sequential quadratic programming algorithm will induce at least R-linear convergence of the iterates to such a local minimum. We construct an example of a degenerate nonlinear program with a unique local minimum satisfying the quadratic growth and the MFCQ but for which no positive semidefinite augmented Lagrangian exists. We present numerical results obtained using several nonlinear programming packages on this example and discuss its implications for some algorithms.
Students' Understanding of Quadratic Equations
ERIC Educational Resources Information Center
López, Jonathan; Robles, Izraim; Martínez-Planell, Rafael
2016-01-01
Action-Process-Object-Schema theory (APOS) was applied to study student understanding of quadratic equations in one variable. This required proposing a detailed conjecture (called a genetic decomposition) of mental constructions students may do to understand quadratic equations. The genetic decomposition which was proposed can contribute to help…
Quadratic Programming for Allocating Control Effort
NASA Technical Reports Server (NTRS)
Singh, Gurkirpal
2005-01-01
A computer program calculates an optimal allocation of control effort in a system that includes redundant control actuators. The program implements an iterative (but otherwise single-stage) algorithm of the quadratic-programming type. In general, in the quadratic-programming problem, one seeks the values of a set of variables that minimize a quadratic cost function, subject to a set of linear equality and inequality constraints. In this program, the cost function combines control effort (typically quantified in terms of energy or fuel consumed) and control residuals (differences between commanded and sensed values of variables to be controlled). In comparison with prior control-allocation software, this program offers approximately equal accuracy but much greater computational efficiency. In addition, this program offers flexibility, robustness to actuation failures, and a capability for selective enforcement of control requirements. The computational efficiency of this program makes it suitable for such complex, real-time applications as controlling redundant aircraft actuators or redundant spacecraft thrusters. The program is written in the C language for execution in a UNIX operating system.
Holographic entropy increases in quadratic curvature gravity
NASA Astrophysics Data System (ADS)
Bhattacharjee, Srijit; Sarkar, Sudipta; Wall, Aron C.
2015-09-01
Standard methods for calculating the black hole entropy beyond general relativity are ambiguous when the horizon is nonstationary. We fix these ambiguities in all quadratic curvature gravity theories, by demanding that the entropy be increasing at every time, for linear perturbations to a stationary black hole. Our result matches with the entropy formula found previously in holographic entanglement entropy calculations. We explicitly calculate the entropy increase for Vaidya-like solutions in Ricci-tensor gravity to show that (unlike the Wald entropy) the holographic entropy obeys a second law.
On the connection of the quadratic Lienard equation with an equation for the elliptic functions
NASA Astrophysics Data System (ADS)
Kudryashov, Nikolay A.; Sinelshchikov, Dmitry I.
2015-07-01
The quadratic Lienard equation is widely used in many applications. A connection between this equation and a linear second-order differential equation has been discussed. Here we show that the whole family of quadratic Lienard equations can be transformed into an equation for the elliptic functions. We demonstrate that this connection can be useful for finding explicit forms of general solutions of the quadratic Lienard equation. We provide several examples of application of our approach.
Students' understanding of quadratic equations
NASA Astrophysics Data System (ADS)
López, Jonathan; Robles, Izraim; Martínez-Planell, Rafael
2016-05-01
Action-Process-Object-Schema theory (APOS) was applied to study student understanding of quadratic equations in one variable. This required proposing a detailed conjecture (called a genetic decomposition) of mental constructions students may do to understand quadratic equations. The genetic decomposition which was proposed can contribute to help students achieve an understanding of quadratic equations with improved interrelation of ideas and more flexible application of solution methods. Semi-structured interviews with eight beginning undergraduate students explored which of the mental constructions conjectured in the genetic decomposition students could do, and which they had difficulty doing. Two of the mental constructions that form part of the genetic decomposition are highlighted and corresponding further data were obtained from the written work of 121 undergraduate science and engineering students taking a multivariable calculus course. The results suggest the importance of explicitly considering these two highlighted mental constructions.
A quadratic weight selection algorithm. [for optimal flight control
NASA Technical Reports Server (NTRS)
Broussard, J. R.
1981-01-01
A new numerical algorithm is presented which determines a positive semi-definite state weighting matrix in the linear-quadratic optimal control design problem. The algorithm chooses the weighting matrix by placing closed-loop eigenvalues and eigenvectors near desired locations using optimal feedback gains. A simplified flight control design example is used to illustrate the algorithms capabilities.
Quadratic invariants for discrete clusters of weakly interacting waves
NASA Astrophysics Data System (ADS)
Harper, Katie L.; Bustamante, Miguel D.; Nazarenko, Sergey V.
2013-06-01
We consider discrete clusters of quasi-resonant triads arising from a Hamiltonian three-wave equation. A cluster consists of N modes forming a total of M connected triads. We investigate the problem of constructing a functionally independent set of quadratic constants of motion. We show that this problem is equivalent to an underlying basic linear problem, consisting of finding the null space of a rectangular M × N matrix {A} with entries 1, -1 and 0. In particular, we prove that the number of independent quadratic invariants is equal to J ≡ N - M* ⩾ N - M, where M* is the number of linearly independent rows in {A}. Thus, the problem of finding all independent quadratic invariants is reduced to a linear algebra problem in the Hamiltonian case. We establish that the properties of the quadratic invariants (e.g., locality) are related to the topological properties of the clusters (e.g., types of linkage). To do so, we formulate an algorithm for decomposing large clusters into smaller ones and show how various invariants are related to certain parts of a cluster, including the basic structures leading to M* < M. We illustrate our findings by presenting examples from the Charney-Hasegawa-Mima wave model, and by showing a classification of small (up to three-triad) clusters.
A transient, quadratic nodal method for triangular-Z geometry
DeLorey, T.F.
1993-06-01
Many systematically-derived nodal methods have been developed for Cartesian geometry due to the extensive interest in Light Water Reactors. These methods typically model the transverse-integrated flux as either an analytic or low order polynomial function of position within the node. Recently, quadratic nodal methods have been developed for R-Z and hexagonal geometry. A static and transient quadratic nodal method is developed for triangular-Z geometry. This development is particularly challenging because the quadratic expansion in each node must be performed between the node faces and the triangular points. As a consequence, in the 2-D plane, the flux and current at the points of the triangles must be treated. Quadratic nodal equations are solved using a non-linear iteration scheme, which utilizes the corrected, mesh-centered finite difference equations, and forces these equations to match the quadratic equations by computing discontinuity factors during the solution. Transient nodal equations are solved using the improved quasi-static method, which has been shown to be a very efficient solution method for transient problems. Several static problems are used to compare the quadratic nodal method to the Coarse Mesh Finite Difference (CMFD) method. The quadratic method is shown to give more accurate node-averaged fluxes. However, it appears that the method has difficulty predicting node leakages near reactor boundaries and severe material interfaces. The consequence is that the eigenvalue may be poorly predicted for certain reactor configurations. The transient methods are tested using a simple analytic test problem, a heterogeneous heavy water reactor benchmark problem, and three thermal hydraulic test problems. Results indicate that the transient methods have been implemented correctly.
Quadratic expressions by means of `summing all the matchsticks'
NASA Astrophysics Data System (ADS)
Faaiz Gierdien, M.
2012-09-01
This note presents demonstrations of quadratic expressions that come about when particular problems are posed with respect to matchsticks that form regular triangles, squares, pentagons and so on. Usually when such 'matchstick' problems are used as ways to foster algebraic thinking, the expressions for the number of matchstick quantities are linear and not quadratic. It will be shown that a pedagogy of 'summing all the matchsticks' is central to the emergence of quadratic expressions. This pedagogy involves generational and transformational activities which are considered as some of the main activities of algebra. Key elements to these activities are processes such as recognizing and extending patterns, and specializing and generalizing particular functional relationships. Implications of these processes in terms of algebraic thinking are considered.
Orthogonality preserving infinite dimensional quadratic stochastic operators
Akın, Hasan; Mukhamedov, Farrukh
2015-09-18
In the present paper, we consider a notion of orthogonal preserving nonlinear operators. We introduce π-Volterra quadratic operators finite and infinite dimensional settings. It is proved that any orthogonal preserving quadratic operator on finite dimensional simplex is π-Volterra quadratic operator. In infinite dimensional setting, we describe all π-Volterra operators in terms orthogonal preserving operators.
Compact stars with quadratic equation of state
NASA Astrophysics Data System (ADS)
Ngubelanga, Sifiso A.; Maharaj, Sunil D.; Ray, Subharthi
2015-05-01
We provide new exact solutions to the Einstein-Maxwell system of equations for matter configurations with anisotropy and charge. The spacetime is static and spherically symmetric. A quadratic equation of state is utilised for the matter distribution. By specifying a particular form for one of the gravitational potentials and the electric field intensity we obtain new exact solutions in isotropic coordinates. In our general class of models, an earlier model with a linear equation of state is regained. For particular choices of parameters we regain the masses of the stars PSR J1614-2230, 4U 1608-52, PSR J1903+0327, EXO 1745-248 and SAX J1808.4-3658. A comprehensive physical analysis for the star PSR J1903+0327 reveals that our model is reasonable.
Feedback linearization application for LLRF control system
Kwon, S.; Regan, A.; Wang, Y.M.; Rohlev, T.
1999-06-01
The Low Energy Demonstration Accelerator (LEDA) being constructed at Los Alamos National Laboratory will serve as the prototype for the low energy section of Acceleration Production of Tritium (APT) accelerator. This paper addresses the problem of the LLRF control system for LEDA. The authors propose a control law which is based on exact feedback linearization coupled with gain scheduling which reduces the effect of the deterministic klystron cathode voltage ripple that is due to harmonics of the high voltage power supply and achieves tracking of desired set points. Also, they propose an estimator of the ripple and its time derivative and the estimates based feedback linearization controller.
Feedback linearization application for LLRF control system
Kwon, S.; Regan, A.; Wang, Y.M.; Rohlev, T.
1998-12-31
The Low Energy Demonstration Accelerator (LEDA) being constructed at Los Alamos National Laboratory will serve as the prototype for the low energy section of Acceleration Production of Tritium (APT) accelerator. This paper addresses the problem of the LLRF control system for LEDA. The authors propose a control law which is based on exact feedback linearization coupled with gain scheduling which reduces the effect of the deterministic klystron cathode voltage ripple that is due to harmonics of the high voltage power supply and achieves tracking of desired set points. Also, they propose an estimator of the ripple and its time derivative and the estimates based feedback linearization controller.
Coherent states for quadratic Hamiltonians
NASA Astrophysics Data System (ADS)
Contreras-Astorga, Alonso; Fernández C, David J.; Velázquez, Mercedes
2011-01-01
The coherent states for a set of quadratic Hamiltonians in the trap regime are constructed. A matrix technique which allows us to directly identify the creation and annihilation operators will be presented. Then, the coherent states as simultaneous eigenstates of the annihilation operators will be derived, and will be compared with those attained through the displacement operator method. The corresponding wavefunction will be found, and a general procedure for obtaining several mean values involving the canonical operators in these states will be described. The results will be illustrated through the asymmetric Penning trap.
Observers for Systems with Nonlinearities Satisfying an Incremental Quadratic Inequality
NASA Technical Reports Server (NTRS)
Acikmese, Ahmet Behcet; Corless, Martin
2004-01-01
We consider the problem of state estimation for nonlinear time-varying systems whose nonlinearities satisfy an incremental quadratic inequality. These observer results unifies earlier results in the literature; and extend it to some additional classes of nonlinearities. Observers are presented which guarantee that the state estimation error exponentially converges to zero. Observer design involves solving linear matrix inequalities for the observer gain matrices. Results are illustrated by application to a simple model of an underwater.
Geometrical and Graphical Solutions of Quadratic Equations.
ERIC Educational Resources Information Center
Hornsby, E. John, Jr.
1990-01-01
Presented are several geometrical and graphical methods of solving quadratic equations. Discussed are Greek origins, Carlyle's method, von Staudt's method, fixed graph methods and imaginary solutions. (CW)
Multivariable quadratic synthesis of an advanced turbofan engine controller
NASA Technical Reports Server (NTRS)
Dehoff, R. L.; Hall, W. E., Jr.
1978-01-01
A digital controller for an advanced turbofan engine utilizing multivariate feedback is described. The theoretical background of locally linearized control synthesis is reviewed briefly. The application of linear quadratic regulator techniques to the practical control problem is presented. The design procedure has been applied to the F100 turbofan engine, and details of the structure of this system are explained. Selected results from simulations of the engine and controller are utilized to illustrate the operation of the system. It is shown that the general multivariable design procedure will produce practical and implementable controllers for modern, high-performance turbine engines.
Binary Quadratic Forms: A Historical View
ERIC Educational Resources Information Center
Khosravani, Azar N.; Beintema, Mark B.
2006-01-01
We present an expository account of the development of the theory of binary quadratic forms. Beginning with the formulation and proof of the Two-Square Theorem, we show how the study of forms of the type x[squared] + ny[squared] led to the discovery of the Quadratic Reciprocity Law, and how this theorem, along with the concept of reduction relates…
An Unexpected Influence on a Quadratic
ERIC Educational Resources Information Center
Davis, Jon D.
2013-01-01
Using technology to explore the coefficients of a quadratic equation can lead to an unexpected result. This article describes an investigation that involves sliders and dynamically linked representations. It guides students to notice the effect that the parameter "a" has on the graphical representation of a quadratic function in the form…
Factorising a Quadratic Expression with Geometric Insights
ERIC Educational Resources Information Center
Joarder, Anwar H.
2015-01-01
An algorithm is presented for factorising a quadratic expression to facilitate instruction and learning. It appeals to elementary geometry which may provide better insights to some students or teachers. There have been many methods for factorising a quadratic expression described in school text books. However, students often seem to struggle with…
A Quadratic Closure for Compressible Turbulence
Futterman, J A
2008-09-16
We have investigated a one-point closure model for compressible turbulence based on third- and higher order cumulant discard for systems undergoing rapid deformation, such as might occur downstream of a shock or other discontinuity. In so doing, we find the lowest order contributions of turbulence to the mean flow, which lead to criteria for Adaptive Mesh Refinement. Rapid distortion theory (RDT) as originally applied by Herring closes the turbulence hierarchy of moment equations by discarding third order and higher cumulants. This is similar to the fourth-order cumulant discard hypothesis of Millionshchikov, except that the Millionshchikov hypothesis was taken to apply to incompressible homogeneous isotropic turbulence generally, whereas RDT is applied only to fluids undergoing a distortion that is 'rapid' in the sense that the interaction of the mean flow with the turbulence overwhelms the interaction of the turbulence with itself. It is also similar to Gaussian closure, in which both second and fourth-order cumulants are retained. Motivated by RDT, we develop a quadratic one-point closure for rapidly distorting compressible turbulence, without regard to homogeneity or isotropy, and make contact with two equation turbulence models, especially the K-{var_epsilon} and K-L models, and with linear instability growth. In the end, we arrive at criteria for Adaptive Mesh Refinement in Finite Volume simulations.
NASA Astrophysics Data System (ADS)
Louboutin, Stephane
1992-07-01
Starting from the analytic class number formula involving its L-function, we first give an expression for the class number of an imaginary quadratic field which, in the case of large discriminants, provides us with a much more powerful numerical technique than that of counting the number of reduced definite positive binary quadratic forms, as has been used by Buell in order to compute his class number tables. Then, using class field theory, we will construct a periodic character &chi , defined on the ring of integers of a field K that is a quadratic extension of a principal imaginary quadratic field k, such that the zeta function of K is the product of the zeta function of k and of the L-function L(s,χ) . We will then determine an integral representation of this L-function that enables us to calculate the class number of K numerically, as soon as its regulator is known. It will also provide us with an upper bound for these class numbers, showing that Hua's bound for the class numbers of imaginary and real quadratic fields is not the best that one could expect. We give statistical results concerning the class numbers of the first 50000 quadratic extensions of {Q}(i) with prime relative discriminant (and with K/Q a non-Galois quartic extension). Our analytic calculation improves the algebraic calculation used by Lakein in the same way as the analytic calculation of the class numbers of real quadratic fields made by Williams and Broere improved the algebraic calculation consisting in counting the number of cycles of reduced ideals. Finally, we give upper bounds for class numbers of K that is a quadratic extension of an imaginary quadratic field k which is no longer assumed to be of class number one.
NASA Astrophysics Data System (ADS)
Hansen, Mikkel Bo; Christiansen, Ove; Hättig, Christof
2009-10-01
Quadratic response functions are derived and implemented for a vibrational configuration interaction state. Combined electronic and vibrational quadratic response functions are derived using Born-Oppenheimer vibronic product wave functions. Computational tractable expressions are derived for determining the total quadratic response contribution as a sum of contributions involving both electronic and vibrational linear and quadratic response functions. In the general frequency-dependent case this includes a new and more troublesome type of electronic linear response function. Pilot calculations for the FH, H2O, CH2O, and pyrrole molecules demonstrate the importance of vibrational contributions for accurate comparison to experiment and that the vibrational contributions in some cases can be very large. The calculation of transition properties between vibrational states is combined with sum-over-states expressions for analysis purposes. On the basis of this some simple analysis methods are suggested. Also, a preliminary study of the effect of finite lifetimes on quadratic response functions is presented.
Seven Wonders of the Ancient and Modern Quadratic World.
ERIC Educational Resources Information Center
Taylor, Sharon E.; Mittag, Kathleen Cage
2001-01-01
Presents four methods for solving a quadratic equation using graphing calculator technology: (1) graphing with the CALC feature; (2) quadratic formula program; (3) table; and (4) solver. Includes a worksheet for a lab activity on factoring quadratic equations. (KHR)
Quadratic Stochastic Operators with Countable State Space
NASA Astrophysics Data System (ADS)
Ganikhodjaev, Nasir
2016-03-01
In this paper, we provide the classes of Poisson and Geometric quadratic stochastic operators with countable state space, study the dynamics of these operators and discuss their application to economics.
Schur Stability Regions for Complex Quadratic Polynomials
ERIC Educational Resources Information Center
Cheng, Sui Sun; Huang, Shao Yuan
2010-01-01
Given a quadratic polynomial with complex coefficients, necessary and sufficient conditions are found in terms of the coefficients such that all its roots have absolute values less than 1. (Contains 3 figures.)
Quantum integrability of quadratic Killing tensors
Duval, C.; Valent, G.
2005-05-01
Quantum integrability of classical integrable systems given by quadratic Killing tensors on curved configuration spaces is investigated. It is proven that, using a 'minimal' quantization scheme, quantum integrability is ensured for a large class of classic examples.
Measurement of quadratic electrogyration effect in castor oil
NASA Astrophysics Data System (ADS)
Izdebski, Marek; Ledzion, Rafał; Górski, Piotr
2015-07-01
This work presents a detailed analysis of electrogyration measurement in liquids with the usage of an optical polarimetric technique. Theoretical analysis of the optical response to an applied electric field is illustrated by experimental data for castor oil which exhibits natural optical activity, quadratic electro-optic effect and quadratic electrogyration effect. Moreover, the experimental data show that interaction of the oil with a pair of flat electrodes induces a significant dichroism and natural linear birefringence. The combination of these effects occurring at the same time complicates the procedure of measurements. It has been found that a single measurement is insufficient to separate the contribution of the electrogyration effect, but it is possible on the basis of several measurements performed with various orientations of the polarizer and the analyser. The obtained average values of the quadratic electrogyration coefficient β13 in castor oil at room temperature are from - 0.92 ×10-22 to - 1.44 ×10-22m2V-2 depending on the origin of the oil. Although this study is focused on measurements in castor oil, the presented analysis is much more general.
Weight of quadratic forms and graph states
NASA Astrophysics Data System (ADS)
Cosentino, Alessandro; Severini, Simone
2009-11-01
We prove a connection between Schmidt rank and weight of quadratic forms. This provides a new tool for the classification of graph states based on entanglement. Our main tool arises from a reformulation of previously known results concerning the weight of quadratic forms in terms of graph states properties. As a byproduct, we obtain a straightforward characterization of the weight of functions associated with pivot-minor of bipartite graphs.
Phase recovery based on quadratic programming
NASA Astrophysics Data System (ADS)
Zhang, Quan Bing; Ge, Xiao Juan; Cheng, Ya Dong; Ni, Na
2014-11-01
Most of the information of optical wavefront is encoded in the phase which includes more details of the object. Conventional optical measuring apparatus is relatively easy to record the intensity of light, but can not measure the phase of light directly. Thus it is important to recovery the phase from the intensity measurements of the object. In recent years, the methods based on quadratic programming such as PhaseLift and PhaseCut can recover the phase of general signal exactly for overdetermined system. To retrieve the phase of sparse signal, the Compressive Phase Retrieval (CPR) algorithm combines the l1-minimization in Compressive Sensing (CS) with low-rank matrix completion problem in PhaseLift, but the result is unsatisfied. This paper focus on the recovery of the phase of sparse signal and propose a new method called the Compressive Phase Cut Retrieval (CPCR) by combining the CPR algorithm with the PhaseCut algorithm. To ensure the sparsity of the recovered signal, we use CPR method to solve a semi-definite programming problem firstly. Then apply linear transformation to the recovered signal, and set the phase of the result as the initial value of the PhaseCut problem. We use TFOCS (a library of Matlab-files) to implement the proposed CPCR algorithm in order to improve the recovered results of the CPR algorithm. Experimental results show that the proposed method can improve the accuracy of the CPR algorithm, and overcome the shortcoming of the PhaseCut method that it can not recover the sparse signal effectively.
The Factorability of Quadratics: Motivation for More Techniques
ERIC Educational Resources Information Center
Bosse, Michael J.; Nandakumar, N. R.
2005-01-01
Typically, secondary and college algebra students attempt to utilize either completing the square or the quadratic formula as techniques to solve a quadratic equation only after frustration with factoring has arisen. While both completing the square and the quadratic formula are techniques which can determine solutions for all quadratic equations,…
Fast approximate quadratic programming for graph matching.
Vogelstein, Joshua T; Conroy, John M; Lyzinski, Vince; Podrazik, Louis J; Kratzer, Steven G; Harley, Eric T; Fishkind, Donniell E; Vogelstein, R Jacob; Priebe, Carey E
2015-01-01
Quadratic assignment problems arise in a wide variety of domains, spanning operations research, graph theory, computer vision, and neuroscience, to name a few. The graph matching problem is a special case of the quadratic assignment problem, and graph matching is increasingly important as graph-valued data is becoming more prominent. With the aim of efficiently and accurately matching the large graphs common in big data, we present our graph matching algorithm, the Fast Approximate Quadratic assignment algorithm. We empirically demonstrate that our algorithm is faster and achieves a lower objective value on over 80% of the QAPLIB benchmark library, compared with the previous state-of-the-art. Applying our algorithm to our motivating example, matching C. elegans connectomes (brain-graphs), we find that it efficiently achieves performance. PMID:25886624
Fast Approximate Quadratic Programming for Graph Matching
Vogelstein, Joshua T.; Conroy, John M.; Lyzinski, Vince; Podrazik, Louis J.; Kratzer, Steven G.; Harley, Eric T.; Fishkind, Donniell E.; Vogelstein, R. Jacob; Priebe, Carey E.
2015-01-01
Quadratic assignment problems arise in a wide variety of domains, spanning operations research, graph theory, computer vision, and neuroscience, to name a few. The graph matching problem is a special case of the quadratic assignment problem, and graph matching is increasingly important as graph-valued data is becoming more prominent. With the aim of efficiently and accurately matching the large graphs common in big data, we present our graph matching algorithm, the Fast Approximate Quadratic assignment algorithm. We empirically demonstrate that our algorithm is faster and achieves a lower objective value on over 80% of the QAPLIB benchmark library, compared with the previous state-of-the-art. Applying our algorithm to our motivating example, matching C. elegans connectomes (brain-graphs), we find that it efficiently achieves performance. PMID:25886624
Limit cycles near hyperbolas in quadratic systems
NASA Astrophysics Data System (ADS)
Artés, Joan C.; Dumortier, Freddy; Llibre, Jaume
In this paper we introduce the notion of infinity strip and strip of hyperbolas as organizing centers of limit cycles in polynomial differential systems on the plane. We study a strip of hyperbolas occurring in some quadratic systems. We deal with the cyclicity of the degenerate graphics DI2a from the programme, set up in [F. Dumortier, R. Roussarie, C. Rousseau, Hilbert's 16th problem for quadratic vector fields, J. Differential Equations 110 (1994) 86-133], to solve the finiteness part of Hilbert's 16th problem for quadratic systems. Techniques from geometric singular perturbation theory are combined with the use of the Bautin ideal. We also rely on the theory of Darboux integrability.
Quadratic forms of projective spaces over rings
NASA Astrophysics Data System (ADS)
Levchuk, V. M.; Starikova, O. A.
2006-06-01
In the passage from fields to rings of coefficients quadratic forms with invertible matrices lose their decisive role. It turns out that if all quadratic forms over a ring are diagonalizable, then in effect this is always a local principal ideal ring R with 2\\in R^*. The problem of the construction of a `normal' diagonal form of a quadratic form over a ring R faces obstacles in the case of indices \\vert R^*:R^{*2}\\vert greater than 1. In the case of index 2 this problem has a solution given in Theorem 2.1 for 1+R^{*2}\\subseteq R^{*2} (an extension of the law of inertia for real quadratic forms) and in Theorem 2.2 for 1+R^2 containing an invertible non-square. Under the same conditions on a ring R with nilpotent maximal ideal the number of classes of projectively congruent quadratic forms of the projective space associated with a free R-module of rank n is explicitly calculated (Proposition 3.2). Up to projectivities, the list of forms is presented for the projective plane over R and also (Theorem 3.3) over the local ring F\\lbrack\\lbrack x,y\\rbrack\\rbrack/\\langle x^{2},xy,y^{2}\\rangle with non-principal maximal ideal, where F=2F is a field with an invertible non-square in 1+F^{2} and \\vert F^{*}:F^{*2}\\vert=2. In the latter case the number of classes of non-diagonalizable quadratic forms of rank 0 depends on one's choice of the field F and is not even always finite; all the other forms make up 21 classes.
Quintessence with quadratic coupling to dark matter
Boehmer, Christian G.; Chan, Nyein; Caldera-Cabral, Gabriela; Lazkoz, Ruth; Maartens, Roy
2010-04-15
We introduce a new form of coupling between dark energy and dark matter that is quadratic in their energy densities. Then we investigate the background dynamics when dark energy is in the form of exponential quintessence. The three types of quadratic coupling all admit late-time accelerating critical points, but these are not scaling solutions. We also show that two types of coupling allow for a suitable matter era at early times and acceleration at late times, while the third type of coupling does not admit a suitable matter era.
Heredity in one-dimensional quadratic maps
NASA Astrophysics Data System (ADS)
Romera, M.; Pastor, G.; Alvarez, G.; Montoya, F.
1998-12-01
In an iterative process, as is the case of a one-dimensional quadratic map, heredity has never been mentioned. In this paper we show that the pattern of a superstable orbit of a one-dimensional quadratic map can be expressed as the sum of the gene of the chaotic band where the pattern is to be found, and the ancestral path that joins all its ancestors. The ancestral path holds all the needed genetic information to calculate the descendants of the pattern. The ancestral path and successive descendant generations of the pattern constitute the family tree of the pattern, which is important to study and understand the orbit's ordering.
On orthogonality preserving quadratic stochastic operators
Mukhamedov, Farrukh; Taha, Muhammad Hafizuddin Mohd
2015-05-15
A quadratic stochastic operator (in short QSO) is usually used to present the time evolution of differing species in biology. Some quadratic stochastic operators have been studied by Lotka and Volterra. In the present paper, we first give a simple characterization of Volterra QSO in terms of absolutely continuity of discrete measures. Further, we introduce a notion of orthogonal preserving QSO, and describe such kind of operators defined on two dimensional simplex. It turns out that orthogonal preserving QSOs are permutations of Volterra QSO. The associativity of genetic algebras generated by orthogonal preserving QSO is studied too.
Quadratic-Like Dynamics of Cubic Polynomials
NASA Astrophysics Data System (ADS)
Blokh, Alexander; Oversteegen, Lex; Ptacek, Ross; Timorin, Vladlen
2016-02-01
A small perturbation of a quadratic polynomial f with a non-repelling fixed point gives a polynomial g with an attracting fixed point and a Jordan curve Julia set, on which g acts like angle doubling. However, there are cubic polynomials with a non-repelling fixed point, for which no perturbation results into a polynomial with Jordan curve Julia set. Motivated by the study of the closure of the Cubic Principal Hyperbolic Domain, we describe such polynomials in terms of their quadratic-like restrictions.
Guises and disguises of quadratic divergences
Cherchiglia, A.L.; Vieira, A.R.; Hiller, Brigitte; Baêta Scarpelli, A.P.; Sampaio, Marcos
2014-12-15
In this contribution, we present a new perspective on the control of quadratic divergences in quantum field theory, in general, and in the Higgs naturalness problem, in particular. Our discussion is essentially based on an approach where UV divergences are parameterized, after being reduced to basic divergent integrals (BDI) in one internal momentum, as functions of a cutoff and a renormalization group scale λ. We illustrate our proposal with well-known examples, such as the gluon vacuum self energy of QCD and the Higgs decay in two photons within this approach. We also discuss frameworks in effective low-energy QCD models, where quadratic divergences are indeed fundamental.
Quadratic minima and modular forms II
NASA Astrophysics Data System (ADS)
Brent, Barry
We give upper bounds on the size of the gap between a non-zero constant term and the next non-zero Fourier coefficient of an entire level two modular form. We give upper bounds for the minimum positive integer represented by a level two even positive-definite quadratic form. These bounds extend partial results in part I.
Curious Consequences of a Miscopied Quadratic
ERIC Educational Resources Information Center
Poet, Jeffrey L.; Vestal, Donald L., Jr.
2005-01-01
The starting point of this article is a search for pairs of quadratic polynomials x[superscript 2] + bx plus or minus c with the property that they both factor over the integers. The search leads quickly to some number theory in the form of primitive Pythagorean triples, and this paper develops the connection between these two topics.
Integration of the Quadratic Function and Generalization
ERIC Educational Resources Information Center
Mitsuma, Kunio
2011-01-01
We will first recall useful formulas in integration that simplify the calculation of certain definite integrals with the quadratic function. A main formula relies only on the coefficients of the function. We will then explore a geometric proof of one of these formulas. Finally, we will extend the formulas to more general cases. (Contains 3…
User's guide for SOL/QPSOL: a Fortran package for quadratic programming
Gill, P.E.; Murray, W.; Saunders, M.A.; Wright, M.H.
1983-07-01
This report forms the user's guide for Version 3.1 of SOL/QPSOL, a set of Fortran subroutines designed to locate the minimum value of an arbitrary quadratic function subject to linear constraints and simple upper and lower bounds. If the quadratic function is convex, a global minimum is found; otherwise, a local minimum is found. The method used is most efficient when many constraints or bounds are active at the solution. QPSOL treats the Hessian and general constraints as dense matrices, and hence is not intended for large sparse problems. This document replaces the previous user's guide of June 1982.
Campoamor-Stursberg, R.
2008-05-15
By means of contractions of Lie algebras, we obtain new classes of indecomposable quasiclassical Lie algebras that satisfy the Yang-Baxter equations in its reformulation in terms of triple products. These algebras are shown to arise naturally from noncompact real simple algebras with nonsimple complexification, where we impose that a nondegenerate quadratic Casimir operator is preserved by the limiting process. We further consider the converse problem and obtain sufficient conditions on integrable cocycles of quasiclassical Lie algebras in order to preserve nondegenerate quadratic Casimir operators by the associated linear deformations.
NASA Astrophysics Data System (ADS)
Swaidan, Waleeda; Hussin, Amran
2015-10-01
Most direct methods solve finite time horizon optimal control problems with nonlinear programming solver. In this paper, we propose a numerical method for solving nonlinear optimal control problem with state and control inequality constraints. This method used quasilinearization technique and Haar wavelet operational matrix to convert the nonlinear optimal control problem into a quadratic programming problem. The linear inequality constraints for trajectories variables are converted to quadratic programming constraint by using Haar wavelet collocation method. The proposed method has been applied to solve Optimal Control of Multi-Item Inventory Model. The accuracy of the states, controls and cost can be improved by increasing the Haar wavelet resolution.
An Instability Index Theory for Quadratic Pencils and Applications
NASA Astrophysics Data System (ADS)
Bronski, Jared; Johnson, Mathew A.; Kapitula, Todd
2014-04-01
Primarily motivated by the stability analysis of nonlinear waves in second-order in time Hamiltonian systems, in this paper we develop an instability index theory for quadratic operator pencils acting on a Hilbert space. In an extension of the known theory for linear pencils, explicit connections are made between the number of eigenvalues of a given quadratic operator pencil with positive real parts to spectral information about the individual operators comprising the coefficients of the spectral parameter in the pencil. As an application, we apply the general theory developed here to yield spectral and nonlinear stability/instability results for abstract second-order in time wave equations. More specifically, we consider the problem of the existence and stability of spatially periodic waves for the "good" Boussinesq equation. In the analysis our instability index theory provides an explicit, and somewhat surprising, connection between the stability of a given periodic traveling wave solution of the "good" Boussinesq equation and the stability of the same periodic profile, but with different wavespeed, in the nonlinear dynamics of a related generalized Korteweg-de Vries equation.
Confidence set inference with a prior quadratic bound
NASA Technical Reports Server (NTRS)
Backus, George E.
1989-01-01
In the uniqueness part of a geophysical inverse problem, the observer wants to predict all likely values of P unknown numerical properties z=(z sub 1,...,z sub p) of the earth from measurement of D other numerical properties y (sup 0) = (y (sub 1) (sup 0), ..., y (sub D (sup 0)), using full or partial knowledge of the statistical distribution of the random errors in y (sup 0). The data space Y containing y(sup 0) is D-dimensional, so when the model space X is infinite-dimensional the linear uniqueness problem usually is insoluble without prior information about the correct earth model x. If that information is a quadratic bound on x, Bayesian inference (BI) and stochastic inversion (SI) inject spurious structure into x, implied by neither the data nor the quadratic bound. Confidence set inference (CSI) provides an alternative inversion technique free of this objection. Confidence set inference is illustrated in the problem of estimating the geomagnetic field B at the core-mantle boundary (CMB) from components of B measured on or above the earth's surface.
Factorization using the quadratic sieve algorithm
Davis, J.A.; Holdridge, D.B.
1983-12-01
Since the cryptosecurity of the RSA two key cryptoalgorithm is no greater than the difficulty of factoring the modulus (product of two secret primes), a code that implements the Quadratic Sieve factorization algorithm on the CRAY I computer has been developed at the Sandia National Laboratories to determine as sharply as possible the current state-of-the-art in factoring. Because all viable attacks on RSA thus far proposed are equivalent to factorization of the modulus, sharper bounds on the computational difficulty of factoring permit improved estimates for the size of RSA parameters needed for given levels of cryptosecurity. Analysis of the Quadratic Sieve indicates that it may be faster than any previously published general purpose algorithm for factoring large integers. The high speed of the CRAY I coupled with the capability of the CRAY to pipeline certain vectorized operations make this algorithm (and code) the front runner in current factoring techniques.
Factorization using the quadratic sieve algorithm
Davis, J.A.; Holdridge, D.B.
1983-01-01
Since the cryptosecurity of the RSA two key cryptoalgorithm is no greater than the difficulty of factoring the modulus (product of two secret primes), a code that implements the Quadratic Sieve factorization algorithm on the CRAY I computer has been developed at the Sandia National Laboratories to determine as sharply as possible the current state-of-the-art in factoring. Because all viable attacks on RSA thus far proposed are equivalent to factorization of the modulus, sharper bounds on the computational difficulty of factoring permit improved estimates for the size of RSA parameters needed for given levels of cryptosecurity. Analysis of the Quadratic Sieve indicates that it may be faster than any previously published general purpose algorithm for factoring large integers. The high speed of the CRAY I coupled with the capability of the CRAY to pipeline certain vectorized operations make this algorithm (and code) the front runner in current factoring techniques.
Characterization of a Quadratic Function in Rn
ERIC Educational Resources Information Center
Xu, Conway
2010-01-01
It is proved that a scalar-valued function "f"(x) defined in "n"-dimensional space must be quadratic, if the intersection of tangent planes at x[subscript 1] and x[subscript 2] always contains the midpoint of the line joining x[subscript 1] and x[subscript 2]. This is the converse of a result of Stenlund proved in this JOURNAL in 2001.
An alternative method on quadratic programming problems
NASA Astrophysics Data System (ADS)
Dasril, Y.; Mohd, I. B.; Mustaffa, I.; Aminuddin, MMM.
2015-05-01
In this paper we proposed an alternative approach to find the optimum solution of quadratic programming problems (QPP) in its original form without additional information such as slack variable, surplus variable or artificial variable as done in other favourite methods. This approached is based on the violated constraints by the unconstrained optimum. The optimal solution of QPP obtained by searching from initial point to another point alongside of feasible region.
THE EFFECTIVENESS OF QUADRATS FOR MEASURING VASCULAR PLANT DIVERSITY
Quadrats are widely used for measuring characteristics of vascular plant communities. It is well recognized that quadrat size affects measurements of frequency and cover. The ability of quadrats of varying sizes to adequately measure diversity has not been established. An exha...
Graphical Solution of the Monic Quadratic Equation with Complex Coefficients
ERIC Educational Resources Information Center
Laine, A. D.
2015-01-01
There are many geometrical approaches to the solution of the quadratic equation with real coefficients. In this article it is shown that the monic quadratic equation with complex coefficients can also be solved graphically, by the intersection of two hyperbolas; one hyperbola being derived from the real part of the quadratic equation and one from…
Developed Adomian method for quadratic Kaluza-Klein relativity
NASA Astrophysics Data System (ADS)
Azreg-Aïnou, Mustapha
2010-01-01
We develop and modify the Adomian decomposition method (ADecM) to work for a new type of nonlinear matrix differential equations (MDE's) which arise in general relativity (GR) and possibly in other applications. The approach consists in modifying both the ADecM linear operator with highest order derivative and ADecM polynomials. We specialize in the case of a 4 × 4 nonlinear MDE along with a scalar one describing stationary cylindrically symmetric metrics in quadratic five-dimensional GR, derive some of their properties using ADecM and construct the most general unique power series solutions. However, because of the constraint imposed on the MDE by the scalar one, the series solutions terminate in closed forms exhausting all possible solutions.
Cosmology for quadratic gravity in generalized Weyl geometry
NASA Astrophysics Data System (ADS)
Beltrán Jiménez, Jose; Heisenberg, Lavinia; Koivisto, Tomi S.
2016-04-01
A class of vector-tensor theories arises naturally in the framework of quadratic gravity in spacetimes with linear vector distortion. Requiring the absence of ghosts for the vector field imposes an interesting condition on the allowed connections with vector distortion: the resulting one-parameter family of connections generalises the usual Weyl geometry with polar torsion. The cosmology of this class of theories is studied, focusing on isotropic solutions wherein the vector field is dominated by the temporal component. De Sitter attractors are found and inhomogeneous perturbations around such backgrounds are analysed. In particular, further constraints on the models are imposed by excluding pathologies in the scalar, vector and tensor fluctuations. Various exact background solutions are presented, describing a constant and an evolving dark energy, a bounce and a self-tuning de Sitter phase. However, the latter two scenarios are not viable under a closer scrutiny.
Wind turbine power tracking using an improved multimodel quadratic approach.
Khezami, Nadhira; Benhadj Braiek, Naceur; Guillaud, Xavier
2010-07-01
In this paper, an improved multimodel optimal quadratic control structure for variable speed, pitch regulated wind turbines (operating at high wind speeds) is proposed in order to integrate high levels of wind power to actively provide a primary reserve for frequency control. On the basis of the nonlinear model of the studied plant, and taking into account the wind speed fluctuations, and the electrical power variation, a multimodel linear description is derived for the wind turbine, and is used for the synthesis of an optimal control law involving a state feedback, an integral action and an output reference model. This new control structure allows a rapid transition of the wind turbine generated power between different desired set values. This electrical power tracking is ensured with a high-performance behavior for all other state variables: turbine and generator rotational speeds and mechanical shaft torque; and smooth and adequate evolution of the control variables. PMID:20434153
Schwarz and multilevel methods for quadratic spline collocation
Christara, C.C.; Smith, B.
1994-12-31
Smooth spline collocation methods offer an alternative to Galerkin finite element methods, as well as to Hermite spline collocation methods, for the solution of linear elliptic Partial Differential Equations (PDEs). Recently, optimal order of convergence spline collocation methods have been developed for certain degree splines. Convergence proofs for smooth spline collocation methods are generally more difficult than for Galerkin finite elements or Hermite spline collocation, and they require stronger assumptions and more restrictions. However, numerical tests indicate that spline collocation methods are applicable to a wider class of problems, than the analysis requires, and are very competitive to finite element methods, with respect to efficiency. The authors will discuss Schwarz and multilevel methods for the solution of elliptic PDEs using quadratic spline collocation, and compare these with domain decomposition methods using substructuring. Numerical tests on a variety of parallel machines will also be presented. In addition, preliminary convergence analysis using Schwarz and/or maximum principle techniques will be presented.
Spectroscopy of one-dimensionally inhomogeneous media with quadratic nonlinearity
Golubkov, A A; Makarov, Vladimir A
2011-11-30
We present a brief review of the results of fifty years of development efforts in spectroscopy of one-dimensionally inhomogeneous media with quadratic nonlinearity. The recent original results obtained by the authors show the fundamental possibility of determining, from experimental data, the coordinate dependences of complex quadratic susceptibility tensor components of a onedimensionally inhomogeneous (along the z axis) medium with an arbitrary frequency dispersion, if the linear dielectric properties of the medium also vary along the z axis and are described by a diagonal tensor of the linear dielectric constant. It is assumed that the medium in question has the form of a plane-parallel plate, whose surfaces are perpendicular to the direction of the inhomogeneity. Using the example of several components of the tensors X{sup (2)}(z, {omega}{sub 1} {+-} {omega}{sub 2}; {omega}{sub 1}, {+-} {omega}{sub 2}), we describe two methods for finding their spatial profiles, which differ in the interaction geometry of plane monochromatic fundamental waves with frequencies {omega}{sub 1} and {omega}{sub 2}. The both methods are based on assessing the intensity of the waves propagating from the plate at the sum or difference frequency and require measurements over a range of angles of incidence of the fundamental waves. Such measurements include two series of additional estimates of the intensities of the waves generated under special conditions by using the test and additional reference plates, which eliminates the need for complicated phase measurements of the complex amplitudes of the waves at the sum (difference) frequency.
NASA Astrophysics Data System (ADS)
Wang, Xiao-Lu; Fan, Xiang-Yu; He, Yong-Ming; Nie, Ren-Shi; Huang, Quan-Hua
2013-08-01
Based on material balance and Darcy's law, the governing equation with the quadratic pressure gradient term was deduced. Then the nonlinear model for fluid flow in a multiple-zone composite reservoir including the quadratic gradient term was established and solved using a Laplace transform. A series of standard log-log type curves of 1-zone (homogeneous), 2-zone and 3-zone reservoirs were plotted and nonlinear flow characteristics were analysed. The type curves governed by the coefficient of the quadratic gradient term (β) gradually deviate from those of a linear model with time elapsing. Qualitative and quantitative analyses were implemented to compare the solutions of the linear and nonlinear models. The results showed that differences of pressure transients between the linear and nonlinear models increase with elapsed time and β. At the end, a successful application of the theoretical model data against the field data shows that the nonlinear model will be a good tool to evaluate formation parameters more accurately.
Li, Xiangrong; Zhao, Xupei; Duan, Xiabin; Wang, Xiaoliang
2015-01-01
It is generally acknowledged that the conjugate gradient (CG) method achieves global convergence—with at most a linear convergence rate—because CG formulas are generated by linear approximations of the objective functions. The quadratically convergent results are very limited. We introduce a new PRP method in which the restart strategy is also used. Moreover, the method we developed includes not only n-step quadratic convergence but also both the function value information and gradient value information. In this paper, we will show that the new PRP method (with either the Armijo line search or the Wolfe line search) is both linearly and quadratically convergent. The numerical experiments demonstrate that the new PRP algorithm is competitive with the normal CG method. PMID:26381742
Smits, Iris A M; Timmerman, Marieke E; Stegeman, Alwin
2016-05-01
Maximum likelihood estimation of the linear factor model for continuous items assumes normally distributed item scores. We consider deviations from normality by means of a skew-normally distributed factor model or a quadratic factor model. We show that the item distributions under a skew-normal factor are equivalent to those under a quadratic model up to third-order moments. The reverse only holds if the quadratic loadings are equal to each other and within certain bounds. We illustrate that observed data which follow any skew-normal factor model can be so well approximated with the quadratic factor model that the models are empirically indistinguishable, and that the reverse does not hold in general. The choice between the two models to account for deviations of normality is illustrated by an empirical example from clinical psychology. PMID:26566696
Homogeneous systems with quadratic integrals, Lie-Poisson quasibrackets, and Kovalevskaya's method
NASA Astrophysics Data System (ADS)
Bizyaev, I. A.; Kozlov, V. V.
2015-12-01
We consider differential equations with quadratic right-hand sides that admit two quadratic first integrals, one of which is a positive-definite quadratic form. We indicate conditions of general nature under which a linear change of variables reduces this system to a certain 'canonical' form. Under these conditions, the system turns out to be divergenceless and can be reduced to a Hamiltonian form, but the corresponding linear Lie-Poisson bracket does not always satisfy the Jacobi identity. In the three-dimensional case, the equations can be reduced to the classical equations of the Euler top, and in four-dimensional space, the system turns out to be superintegrable and coincides with the Euler-Poincaré equations on some Lie algebra. In the five-dimensional case we find a reducing multiplier after multiplying by which the Poisson bracket satisfies the Jacobi identity. In the general case for n>5 we prove the absence of a reducing multiplier. As an example we consider a system of Lotka-Volterra type with quadratic right-hand sides that was studied by Kovalevskaya from the viewpoint of conditions of uniqueness of its solutions as functions of complex time. Bibliography: 38 titles.
Wiedemann, H.
1981-11-01
Since no linear colliders have been built yet it is difficult to know at what energy the linear cost scaling of linear colliders drops below the quadratic scaling of storage rings. There is, however, no doubt that a linear collider facility for a center of mass energy above say 500 GeV is significantly cheaper than an equivalent storage ring. In order to make the linear collider principle feasible at very high energies a number of problems have to be solved. There are two kinds of problems: one which is related to the feasibility of the principle and the other kind of problems is associated with minimizing the cost of constructing and operating such a facility. This lecture series describes the problems and possible solutions. Since the real test of a principle requires the construction of a prototype I will in the last chapter describe the SLC project at the Stanford Linear Accelerator Center.
Inter-annual precipitation changes as quadratic signals in the GRACE time-variable gravity
NASA Astrophysics Data System (ADS)
Ogawa, R.; Chao, B. F.; Heki, K.
2009-04-01
The Gravity Recovery and Climate Experiment (GRACE) satellite mission has been producing scientific results on mass variations on inter-annual timescales, e.g. melting of ice sheet in Greenland and mountain glaciers in Alaska, Eastern Africa drought, water level increase in Caspian Sea, etc. In these discussions only linear trends and the seasonal components have been analyzed in the monthly GRACE time series, whereas little attention has been paid so far to the existence of the quadratic changes which signify the temporal accelerations. With over 6 years of GRACE data and revisiting the time-variable gravity field of various regions, we find that such acceleration/deceleration terms are quite often significantly different from zero. They include East Africa, near Obi River, Caspian Sea, Black Sea, Central Asia, and southern South America, whereof discussions of linear trends without specifying the epochs are inadequate. Here we investigate geophysical implication of these quadratic terms; in particular gravity changes in land areas reflect, to a large extent, soil moisture variations. Soil moisture is the time integration of water fluxes, i.e. precipitation, evapotranspiration and runoff. Here we consider that the linear trend in precipitation is responsible for the quadratic change in gravity, and examine trends of observed precipitation in various regions from CMAP (Climate Prediction Center Merged Analysis of Precipitation). Thus, in order to compare linear trend in CMAP and acceleration in GRACE, we calculate month-to-month difference of equivalent water depth at GRACE grid points, and modeled them with seasonal variations and linear trends. We found good agreement between their geographical distributions although amplitudes are smaller in GRACE, meaning the quadratic gravity changes in the GRACE data do reflect inter-annual changes of precipitation fairly faithfully.
Reduced order parameter estimation using quasilinearization and quadratic programming
NASA Astrophysics Data System (ADS)
Siade, Adam J.; Putti, Mario; Yeh, William W.-G.
2012-06-01
The ability of a particular model to accurately predict how a system responds to forcing is predicated on various model parameters that must be appropriately identified. There are many algorithms whose purpose is to solve this inverse problem, which is often computationally intensive. In this study, we propose a new algorithm that significantly reduces the computational burden associated with parameter identification. The algorithm is an extension of the quasilinearization approach where the governing system of differential equations is linearized with respect to the parameters. The resulting inverse problem therefore becomes a linear regression or quadratic programming problem (QP) for minimizing the sum of squared residuals; the solution becomes an update on the parameter set. This process of linearization and regression is repeated until convergence takes place. This algorithm has not received much attention, as the QPs can become quite large, often infeasible for real-world systems. To alleviate this drawback, proper orthogonal decomposition is applied to reduce the size of the linearized model, thereby reducing the computational burden of solving each QP. In fact, this study shows that the snapshots need only be calculated once at the very beginning of the algorithm, after which no further calculations of the reduced-model subspace are required. The proposed algorithm therefore only requires one linearized full-model run per parameter at the first iteration followed by a series of reduced-order QPs. The method is applied to a groundwater model with about 30,000 computation nodes where as many as 15 zones of hydraulic conductivity are estimated.
A quadratic analog-to-digital converter
NASA Technical Reports Server (NTRS)
Harrison, D. C.; Staples, M. H.
1980-01-01
An analog-to-digital converter with a square root transfer function has been developed for use with a pair of CCD imaging detectors in the White Light Coronagraph/X-ray XUV Telescope experiment to be flown as part of the Internal Solar Polar Mission. It is shown that in background-noise-limited instrumentation systems a quadratic analog-to-digital converter will allow a maximum dynamic range with a fixed number of data bits. Low power dissipation, moderately fast conversion time, and reliability are achieved in the proposed design using standard components and avoiding nonlinear elements.
Finite element simulation of articular contact mechanics with quadratic tetrahedral elements.
Maas, Steve A; Ellis, Benjamin J; Rawlins, David S; Weiss, Jeffrey A
2016-03-21
Although it is easier to generate finite element discretizations with tetrahedral elements, trilinear hexahedral (HEX8) elements are more often used in simulations of articular contact mechanics. This is due to numerical shortcomings of linear tetrahedral (TET4) elements, limited availability of quadratic tetrahedron elements in combination with effective contact algorithms, and the perceived increased computational expense of quadratic finite elements. In this study we implemented both ten-node (TET10) and fifteen-node (TET15) quadratic tetrahedral elements in FEBio (www.febio.org) and compared their accuracy, robustness in terms of convergence behavior and computational cost for simulations relevant to articular contact mechanics. Suitable volume integration and surface integration rules were determined by comparing the results of several benchmark contact problems. The results demonstrated that the surface integration rule used to evaluate the contact integrals for quadratic elements affected both convergence behavior and accuracy of predicted stresses. The computational expense and robustness of both quadratic tetrahedral formulations compared favorably to the HEX8 models. Of note, the TET15 element demonstrated superior convergence behavior and lower computational cost than both the TET10 and HEX8 elements for meshes with similar numbers of degrees of freedom in the contact problems that we examined. Finally, the excellent accuracy and relative efficiency of these quadratic tetrahedral elements was illustrated by comparing their predictions with those for a HEX8 mesh for simulation of articular contact in a fully validated model of the hip. These results demonstrate that TET10 and TET15 elements provide viable alternatives to HEX8 elements for simulation of articular contact mechanics. PMID:26900037
Confidence set inference with a prior quadratic bound
NASA Technical Reports Server (NTRS)
Backus, George E.
1988-01-01
In the uniqueness part of a geophysical inverse problem, the observer wants to predict all likely values of P unknown numerical properties z = (z sub 1,...,z sub p) of the earth from measurement of D other numerical properties y(0)=(y sub 1(0),...,y sub D(0)) knowledge of the statistical distribution of the random errors in y(0). The data space Y containing y(0) is D-dimensional, so when the model space X is infinite-dimensional the linear uniqueness problem usually is insoluble without prior information about the correct earth model x. If that information is a quadratic bound on x (e.g., energy or dissipation rate), Bayesian inference (BI) and stochastic inversion (SI) inject spurious structure into x, implied by neither the data nor the quadratic bound. Confidence set inference (CSI) provides an alternative inversion technique free of this objection. CSI is illustrated in the problem of estimating the geomagnetic field B at the core-mantle boundary (CMB) from components of B measured on or above the earth's surface. Neither the heat flow nor the energy bound is strong enough to permit estimation of B(r) at single points on the CMB, but the heat flow bound permits estimation of uniform averages of B(r) over discs on the CMB, and both bounds permit weighted disc-averages with continous weighting kernels. Both bounds also permit estimation of low-degree Gauss coefficients at the CMB. The heat flow bound resolves them up to degree 8 if the crustal field at satellite altitudes must be treated as a systematic error, but can resolve to degree 11 under the most favorable statistical treatment of the crust. These two limits produce circles of confusion on the CMB with diameters of 25 deg and 19 deg respectively.
Blind deconvolution estimation of fluorescence measurements through quadratic programming
NASA Astrophysics Data System (ADS)
Campos-Delgado, Daniel U.; Gutierrez-Navarro, Omar; Arce-Santana, Edgar R.; Skala, Melissa C.; Walsh, Alex J.; Jo, Javier A.
2015-07-01
Time-deconvolution of the instrument response from fluorescence lifetime imaging microscopy (FLIM) data is usually necessary for accurate fluorescence lifetime estimation. In many applications, however, the instrument response is not available. In such cases, a blind deconvolution approach is required. An iterative methodology is proposed to address the blind deconvolution problem departing from a dataset of FLIM measurements. A linear combination of a base conformed by Laguerre functions models the fluorescence impulse response of the sample at each spatial point in our formulation. Our blind deconvolution estimation (BDE) algorithm is formulated as a quadratic approximation problem, where the decision variables are the samples of the instrument response and the scaling coefficients of the basis functions. In the approximation cost function, there is a bilinear dependence on the decision variables. Hence, due to the nonlinear nature of the estimation process, an alternating least-squares scheme iteratively solves the approximation problem. Our proposal searches for the samples of the instrument response with a global perspective, and the scaling coefficients of the basis functions locally at each spatial point. First, the iterative methodology relies on a least-squares solution for the instrument response, and quadratic programming for the scaling coefficients applied just to a subset of the measured fluorescence decays to initially estimate the instrument response to speed up the convergence. After convergence, the final stage computes the fluorescence impulse response at all spatial points. A comprehensive validation stage considers synthetic and experimental FLIM datasets of ex vivo atherosclerotic plaques and human breast cancer cell samples that highlight the advantages of the proposed BDE algorithm under different noise and initial conditions in the iterative scheme and parameters of the proposal.
Quadratic spline collocation and parareal deferred correction method for parabolic PDEs
NASA Astrophysics Data System (ADS)
Liu, Jun; Wang, Yan; Li, Rongjian
2016-06-01
In this paper, we consider a linear parabolic PDE, and use optimal quadratic spline collocation (QSC) methods for the space discretization, proceed the parareal technique on the time domain. Meanwhile, deferred correction technique is used to improve the accuracy during the iterations. The error estimation is presented and the stability is analyzed. Numerical experiments, which is carried out on a parallel computer with 40 CPUs, are attached to exhibit the effectiveness of the hybrid algorithm.
Endicott, Julia S.; Joubert-Doriol, Loïc; Izmaylov, Artur F.
2014-07-21
We consider a fully quadratic vibronic model Hamiltonian for studying photoinduced electronic transitions through conical intersections. Using a second order perturbative approximation for diabatic couplings, we derive an analytical expression for the time evolution of electronic populations at a given temperature. This formalism extends upon a previously developed perturbative technique for a linear vibronic coupling Hamiltonian. The advantage of the quadratic model Hamiltonian is that it allows one to use separate quadratic representations for potential energy surfaces of different electronic states and a more flexible representation of interstate couplings. We explore features introduced by the quadratic Hamiltonian in a series of 2D models, and then apply our formalism to the 2,6-bis(methylene) adamantyl cation and its dimethyl derivative. The Hamiltonian parameters for the molecular systems have been obtained from electronic structure calculations followed by a diabatization procedure. The evolution of electronic populations in the molecular systems using the perturbative formalism shows a good agreement with that from variational quantum dynamics.
When does brain aging accelerate? Dangers of quadratic fits in cross-sectional studies.
Fjell, Anders M; Walhovd, Kristine B; Westlye, Lars T; Østby, Ylva; Tamnes, Christian K; Jernigan, Terry L; Gamst, Anthony; Dale, Anders M
2010-05-01
Many brain structures show a complex, non-linear pattern of maturation and age-related change. Often, quadratic models (beta(0) + beta(1)age + beta(2)age(2) + epsilon) are used to describe such relationships. Here, we demonstrate that the fitting of quadratic models is substantially affected by seemingly irrelevant factors, such as the age-range sampled. Hippocampal volume was measured in 434 healthy participants between 8 and 85 years of age, and quadratic models were fit to subsets of the sample with different age-ranges. It was found that as the bottom of the age-range increased, the age at which volumes appeared to peak was moved upwards and the estimated decline in the last part of the age-span became larger. Thus, whether children were included or not affected the estimated decline between 60 and 85 years. We conclude that caution should be exerted in inferring age-trajectories from global fit models, e.g. the quadratic model. A nonparametric local smoothing technique (the smoothing spline) was found to be more robust to the effects of different starting ages. The results were replicated in an independent sample of 309 participants. PMID:20109562
Quadratic relations in continuous and discrete Painlevé equations
NASA Astrophysics Data System (ADS)
Ramani, A.; Grammaticos, B.; Tamizhmani, T.
2000-04-01
The quadratic relations between the solutions of a Painlevé equation and that of a different one, or the same one with a different set of parameters, are investigated in the continuous and discrete cases. We show that the quadratic relations existing for the continuous PII , PIII , PV and PVI have analogues as well as consequences in the discrete case. Moreover, the discrete Painlevé equations have quadratic relations of their own without any reference to the continuous case.
Quadratic dynamical decoupling with nonuniform error suppression
Quiroz, Gregory; Lidar, Daniel A.
2011-10-15
We analyze numerically the performance of the near-optimal quadratic dynamical decoupling (QDD) single-qubit decoherence errors suppression method [J. West et al., Phys. Rev. Lett. 104, 130501 (2010)]. The QDD sequence is formed by nesting two optimal Uhrig dynamical decoupling sequences for two orthogonal axes, comprising N{sub 1} and N{sub 2} pulses, respectively. Varying these numbers, we study the decoherence suppression properties of QDD directly by isolating the errors associated with each system basis operator present in the system-bath interaction Hamiltonian. Each individual error scales with the lowest order of the Dyson series, therefore immediately yielding the order of decoherence suppression. We show that the error suppression properties of QDD are dependent upon the parities of N{sub 1} and N{sub 2}, and near-optimal performance is achieved for general single-qubit interactions when N{sub 1}=N{sub 2}.
Large-scale sequential quadratic programming algorithms
Eldersveld, S.K.
1992-09-01
The problem addressed is the general nonlinear programming problem: finding a local minimizer for a nonlinear function subject to a mixture of nonlinear equality and inequality constraints. The methods studied are in the class of sequential quadratic programming (SQP) algorithms, which have previously proved successful for problems of moderate size. Our goal is to devise an SQP algorithm that is applicable to large-scale optimization problems, using sparse data structures and storing less curvature information but maintaining the property of superlinear convergence. The main features are: 1. The use of a quasi-Newton approximation to the reduced Hessian of the Lagrangian function. Only an estimate of the reduced Hessian matrix is required by our algorithm. The impact of not having available the full Hessian approximation is studied and alternative estimates are constructed. 2. The use of a transformation matrix Q. This allows the QP gradient to be computed easily when only the reduced Hessian approximation is maintained. 3. The use of a reduced-gradient form of the basis for the null space of the working set. This choice of basis is more practical than an orthogonal null-space basis for large-scale problems. The continuity condition for this choice is proven. 4. The use of incomplete solutions of quadratic programming subproblems. Certain iterates generated by an active-set method for the QP subproblem are used in place of the QP minimizer to define the search direction for the nonlinear problem. An implementation of the new algorithm has been obtained by modifying the code MINOS. Results and comparisons with MINOS and NPSOL are given for the new algorithm on a set of 92 test problems.
Milyutin, A A
2002-04-30
An integral quadratic functional with constant coefficients on a half-axis is considered. A necessary and sufficient condition for its non-negativity at all square integrable pairs of functions related by a linear ODE is proposed, which is based on the Hamilton-Jacobi inequality. A connection between this condition and the well-known frequency criterion is established.
A gain-scheduled observer under transmissions without delivery acknowledgment
NASA Astrophysics Data System (ADS)
Dolz, Daniel; Peñarrocha, Ignacio; Sanchis, Roberto
2015-11-01
This paper addresses the estimation problem for discrete-time systems where both measurements and control commands are sent to a central station through a lossy network without delivery acknowledgment. The central unit implements the estimation and control algorithms. We propose a jump observer that uses the expected value of the unknown control input at the actuator to run an open loop estimation. Then, the absence of acknowledgment in the control input transmission is dealt with through the introduction of an unknown disturbance. The observer update is performed by means of jumping gains when there are available measurements. We employ an statistic of the control error (new disturbance), i.e., the difference between the control inputs at the plant and at the observer, to schedule the observer gains in real time. The observer is designed to minimize the H∞ norm from disturbances, measurement noises and control input errors, to estimation error. The infinite-dimensional design problem is turn into an optimization problem over polynomials using sum-of-squares decomposition techniques. Benefits of the proposal are shown in a simulation example.
Gain-Scheduled Fault Tolerance Control Under False Identification
NASA Technical Reports Server (NTRS)
Shin, Jong-Yeob; Belcastro, Christine (Technical Monitor)
2006-01-01
An active fault tolerant control (FTC) law is generally sensitive to false identification since the control gain is reconfigured for fault occurrence. In the conventional FTC law design procedure, dynamic variations due to false identification are not considered. In this paper, an FTC synthesis method is developed in order to consider possible variations of closed-loop dynamics under false identification into the control design procedure. An active FTC synthesis problem is formulated into an LMI optimization problem to minimize the upper bound of the induced-L2 norm which can represent the worst-case performance degradation due to false identification. The developed synthesis method is applied for control of the longitudinal motions of FASER (Free-flying Airplane for Subscale Experimental Research). The designed FTC law of the airplane is simulated for pitch angle command tracking under a false identification case.
Analysis of Students' Error in Learning of Quadratic Equations
ERIC Educational Resources Information Center
Zakaria, Effandi; Ibrahim; Maat, Siti Mistima
2010-01-01
The purpose of the study was to determine the students' error in learning quadratic equation. The samples were 30 form three students from a secondary school in Jambi, Indonesia. Diagnostic test was used as the instrument of this study that included three components: factorization, completing the square and quadratic formula. Diagnostic interview…
Tangent Lines without Derivatives for Quadratic and Cubic Equations
ERIC Educational Resources Information Center
Carroll, William J.
2009-01-01
In the quadratic equation, y = ax[superscript 2] + bx + c, the equation y = bx + c is identified as the equation of the line tangent to the parabola at its y-intercept. This is extended to give a convenient method of graphing tangent lines at any point on the graph of a quadratic or a cubic equation. (Contains 5 figures.)
Effects of Classroom Instruction on Students' Understanding of Quadratic Equations
ERIC Educational Resources Information Center
Vaiyavutjamai, Pongchawee; Clements, M. A.
2006-01-01
Two hundred and thirty-one students in six Grade 9 classes in two government secondary schools located near Chiang Mai, Thailand, attempted to solve the same 18 quadratic equations before and after participating in 11 lessons on quadratic equations. Data from the students' written responses to the equations, together with data in the form of…
Visualising the Roots of Quadratic Equations with Complex Coefficients
ERIC Educational Resources Information Center
Bardell, Nicholas S.
2014-01-01
This paper is a natural extension of the root visualisation techniques first presented by Bardell (2012) for quadratic equations with real coefficients. Consideration is now given to the familiar quadratic equation "y = ax[superscript 2] + bx + c" in which the coefficients "a," "b," "c" are generally…
Geometric quadratic stochastic operator on countable infinite set
Ganikhodjaev, Nasir; Hamzah, Nur Zatul Akmar
2015-02-03
In this paper we construct the family of Geometric quadratic stochastic operators defined on the countable sample space of nonnegative integers and investigate their trajectory behavior. Such operators can be reinterpreted in terms of of evolutionary operator of free population. We show that Geometric quadratic stochastic operators are regular transformations.
An Authentic Task That Models Quadratics
ERIC Educational Resources Information Center
Baron, Lorraine M.
2015-01-01
As students develop algebraic reasoning in grades 5 to 9, they learn to recognize patterns and understand expressions, equations, and variables. Linear functions are a focus in eighth-grade mathematics, and by algebra 1, students must make sense of functions that are not linear. This article describes how students worked through a classroom task…
Quadratic divergences and quantum gravitational contributions to gauge coupling constants
NASA Astrophysics Data System (ADS)
Toms, David J.
2011-10-01
The calculation of quadratic divergences in Einstein-Maxwell theory with a possible cosmological constant is considered. We describe a method of calculation, using the background-field method, that is sensitive to quadratic divergences, is respectful of gauge invariance, and is independent of gauge conditions. A standard renormalization group analysis is applied to the result where it is shown that the quadratic divergences do lead to asymptotic freedom as found in the original paper of Robinson and Wilczek. The role and nature of these quadratic divergences is critically evaluated in light of recent criticism. Within the context of the background-field method, it is shown that it is possible to define the charge in a physically motivated way in which the quadratic divergences do not play a role. This latter view is studied in more depth in a toy model described in an appendix.
Impact of a global quadratic potential on galactic rotation curves.
Mannheim, Philip D; O'Brien, James G
2011-03-25
We present a conformal gravity fit to the 20 largest of a sample of 110 spiral galaxies. We identify the presence of a universal quadratic potential V(κ)(r)=-κc²r²/2 with κ=9.54×10⁻⁵⁴ cm⁻² induced by cosmic inhomogeneities. When V(κ)(r) is taken in conjunction with both a universal linear potential V(γ₀)(r)=γ₀c²r/2 with γ₀=3.06×10⁻³⁰ cm⁻¹ generated by the homogeneous cosmic background and the contribution generated by the local luminous matter in galaxies, the theory then accounts for the rotation curve systematics observed in the entire 110 galaxies, without the need for any dark matter whatsoever. Our study suggests that using dark matter may be nothing more than an attempt to describe global effects in purely local galactic terms. With V(κ)(r) being negative, galaxies can only support bound orbits up to distances of order γ₀/κ=100kpc, with global physics imposing a limit on the size of galaxies. PMID:21517292
Quadratic Measurement and Conditional State Preparation in an Optomechanical System
NASA Astrophysics Data System (ADS)
Brawley, George; Vanner, Michael; Bowen, Warwick; Schmid, Silvan; Boisen, Anja
2014-03-01
An important requirement in the study of quantum systems is the ability to measure non-linear observables at the level of quantum fluctuations. Such measurements enable the conditional preparation of highly non-classical states. Nonlinear measurement, although achieved in a variety of quantum systems including microwave cavity modes and optical fields, remains an outstanding problem in both electromechanical and optomechanical systems. To the best of our knowledge, previous experimental efforts to achieve nonlinear measurement of mechanical motion have not yielded strong coupling, nor the observation of quadratic mechanical motion. Here using a new technique reliant on the intrinsic nonlinearity of the optomechanical interaction, we experimentally observe for the first time a position squared (x2) measurement of the room-temperature Brownian motion of a nanomechanical oscillator. We utilize this measurement to conditionally prepare non-Gaussian bimodal states, which are the high temperature classical analogue of quantum macroscopic superposition states, or cat states. In the future with the aid of cryogenics and state-of-the-art optical cavities, our approach will provide a viable method of generating quantum superposition states of mechanical oscillators. This research was funded by the ARC Center of Excellence for Engineered Quantum Systems.
Phase Transitions in the Quadratic Contact Process on Complex Networks
NASA Astrophysics Data System (ADS)
Varghese, Chris; Durrett, Rick
2013-03-01
The quadratic contact process (QCP) is a natural extension of the well studied linear contact process where a single infected (1) individual can infect a susceptible (0) neighbor and infected individuals are allowed to recover (1 --> 0). In the QCP, a combination of two 1's is required to effect a 0 --> 1 change. We extend the study of the QCP, which so far has been limited to lattices, to complex networks as a model for the change in a population via sexual reproduction and death. We define two versions of the QCP - vertex centered (VQCP) and edge centered (EQCP) with birth events 1 - 0 - 1 --> 1 - 1 - 1 and 1 - 1 - 0 --> 1 - 1 - 1 respectively, where ` -' represents an edge. We investigate the effects of network topology by considering the QCP on regular, Erdős-Rényi and power law random graphs. We perform mean field calculations as well as simulations to find the steady state fraction of occupied vertices as a function of the birth rate. We find that on the homogeneous graphs (regular and Erdős-Rényi) there is a discontinuous phase transition with a region of bistability, whereas on the heavy tailed power law graph, the transition is continuous. The critical birth rate is found to be positive in the former but zero in the latter.
The application of quadratic optimal cooperative control synthesis to a CH-47 helicopter
NASA Technical Reports Server (NTRS)
Townsend, Barbara K.
1987-01-01
A control-system design method, quadratic optimal cooperative control synthesis (CCS), is applied to the design of a stability and control augmentation system (SCAS). The CCS design method is different from other design methods in that it does not require detailed a priori design criteria, but instead relies on an explicit optimal pilot-model to create desired performance. The design method, which was developed previously for fixed-wing aircraft, is simplified and modified for application to a Boeing CH-47 helicopter. Two SCAS designs are developed using the CCS design methodology. The resulting CCS designs are then compared with designs obtained using classical/frequency-domain methods and linear quadratic regulator (LQR) theory in a piloted fixed-base simulation. Results indicate that the CCS method, with slight modifications, can be used to produce controller designs which compare favorably with the frequency-domain approach.
Magneto-optical conductivity of Weyl semimetals with quadratic term in momentum
NASA Astrophysics Data System (ADS)
Shao, J. M.; Yang, G. W.
2016-02-01
Weyl semimetal is a three-dimensional Dirac material whose low energy dispersion is linear in momentum. Adding a quadratic (Schrödinger) term to the Weyl node breaks the original particle-hole symmetry and also breaks the mirror symmetry between the positive and negative Landau levels in present of magnetic field. This asymmetry splits the absorption line of the longitudinal magneto-optical conductivity into a two peaks structure. It also results in an oscillation pattern in the absorption part of the Hall conductivity. The two split peaks in Reσxx (or the positive and negative oscillation in Imσxy) just correspond to the absorptions of left-handed (σ-) and right-handed (σ+) polarization light, respectively. The split in Reσxx and the displacement between the absorption of σ+ and σ- are decided by the magnitude of the quadratic term and the magnetic field.
Polychromatic solitons in a quadratic medium.
Towers, I N; Malomed, B A
2002-10-01
We introduce the simplest model to describe parametric interactions in a quadratically nonlinear optical medium with the fundamental harmonic containing two components with (slightly) different carrier frequencies [which is a direct analog of wavelength-division multiplexed models, well known in media with cubic nonlinearity]. The model takes a closed form with three different second-harmonic components, and it is formulated in the spatial domain. We demonstrate that the model supports both polychromatic solitons (PCSs), with all the components present in them, and two types of mutually orthogonal simple solitons, both types being stable in a broad parametric region. An essential peculiarity of PCS is that its power is much smaller than that of a simple (usual) soliton (taken at the same values of control parameters), which may be an advantage for experimental generation of PCSs. Collisions between the orthogonal simple solitons are simulated in detail, leading to the conclusion that the collisions are strongly inelastic, converting the simple solitons into polychromatic ones, and generating one or two additional PCSs. A collision velocity at which the inelastic effects are strongest is identified, and it is demonstrated that the collision may be used as a basis to design a simple all-optical XOR logic gate. PMID:12443362
Optimal channels for channelized quadratic estimators.
Kupinski, Meredith K; Clarkson, Eric
2016-06-01
We present a new method for computing optimized channels for estimation tasks that is feasible for high-dimensional image data. Maximum-likelihood (ML) parameter estimates are challenging to compute from high-dimensional likelihoods. The dimensionality reduction from M measurements to L channels is a critical advantage of channelized quadratic estimators (CQEs), since estimating likelihood moments from channelized data requires smaller sample sizes and inverting a smaller covariance matrix is easier. The channelized likelihood is then used to form ML estimates of the parameter(s). In this work we choose an imaging example in which the second-order statistics of the image data depend upon the parameter of interest: the correlation length. Correlation lengths are used to approximate background textures in many imaging applications, and in these cases an estimate of the correlation length is useful for pre-whitening. In a simulation study we compare the estimation performance, as measured by the root-mean-squared error (RMSE), of correlation length estimates from CQE and power spectral density (PSD) distribution fitting. To abide by the assumptions of the PSD method we simulate an ergodic, isotropic, stationary, and zero-mean random process. These assumptions are not part of the CQE formalism. The CQE method assumes a Gaussian channelized likelihood that can be a valid for non-Gaussian image data, since the channel outputs are formed from weighted sums of the image elements. We have shown that, for three or more channels, the RMSE of CQE estimates of correlation length is lower than conventional PSD estimates. We also show that computing CQE by using a standard nonlinear optimization method produces channels that yield RMSE within 2% of the analytic optimum. CQE estimates of anisotropic correlation length estimation are reported to demonstrate this technique on a two-parameter estimation problem. PMID:27409452
On a 'Mysterious' Case of a Quadratic Hamiltonian
NASA Astrophysics Data System (ADS)
Sakovich, Sergei
2006-07-01
We show that one of the five cases of a quadratic Hamiltonian, which were recently selected by Sokolov and Wolf who used the Kovalevskaya-Lyapunov test, fails to pass the Painlevé test for integrability.
Direct Orthogonal Distance to Quadratic Surfaces in 3D.
Lott, Gus K
2014-09-01
Discovering the orthogonal distance to a quadratic surface is a classic geometric task in vision, modeling, and robotics. I describe a simple, efficient, and stable direct solution for the orthogonal distance (foot-point) to an arbitrary quadratic surface from a general finite 3D point. The problem is expressed as the intersection of three quadratic surfaces, two of which are derived from the requirement of orthogonality of two non-coincident planes with the tangent plane to the quadric. A sixth order single-variable polynomial is directly generated in one coordinate of the surface point. The method detects intersection points at infinity and operates smoothly across all real quadratic surface classes. The method also geometrically detects continuums of orthogonal points (i.e., from the exact center of a sphere). I discuss algorithm performance, compare it to a state-of-the-art estimator, demonstrate the algorithm on synthetic data, and describe extension to arbitrary dimension. PMID:26352239
Use of non-quadratic yield surfaces in design of optimal deep-draw blank geometry
Logan, R.W.
1995-12-01
Planar anisotropy in the deep-drawing of sheet can lead to the formation of ears in cylindrical cups and to undesirable metal flow in the blankholder in the general case. For design analysis purposes in non-linear finite-element codes, this anisotropy is characterized by the use of an appropriate yield surface which is then implemented into codes such as DYNA3D . The quadratic Hill yield surface offers a relatively straightforward implementation and can be formulated to be invariant to the coordinate system. Non-quadratic yield surfaces can provide more realistic strength or strain increment ratios, but they may not provide invariance and thus demand certain approximations. Forms due to Hosford and Badat et al. have been shown to more accurately address the earning phenomenon. in this work, use is made of these non-quadratic yield surfaces in order to determine the optimal blank shape for cups and other shapes using ferrous and other metal blank materials with planar anisotropy. The analyses are compared to previous experimental studies on non-uniform blank motion due to anisotropy and asymmetric geometry.
Culpepper, Steven Andrew
2016-06-01
Standardized tests are frequently used for selection decisions, and the validation of test scores remains an important area of research. This paper builds upon prior literature about the effect of nonlinearity and heteroscedasticity on the accuracy of standard formulas for correcting correlations in restricted samples. Existing formulas for direct range restriction require three assumptions: (1) the criterion variable is missing at random; (2) a linear relationship between independent and dependent variables; and (3) constant error variance or homoscedasticity. The results in this paper demonstrate that the standard approach for correcting restricted correlations is severely biased in cases of extreme monotone quadratic nonlinearity and heteroscedasticity. This paper offers at least three significant contributions to the existing literature. First, a method from the econometrics literature is adapted to provide more accurate estimates of unrestricted correlations. Second, derivations establish bounds on the degree of bias attributed to quadratic functions under the assumption of a monotonic relationship between test scores and criterion measurements. New results are presented on the bias associated with using the standard range restriction correction formula, and the results show that the standard correction formula yields estimates of unrestricted correlations that deviate by as much as 0.2 for high to moderate selectivity. Third, Monte Carlo simulation results demonstrate that the new procedure for correcting restricted correlations provides more accurate estimates in the presence of quadratic and heteroscedastic test score and criterion relationships. PMID:25953477
Adiabatic femtosecond pulse compression and control by using quadratic cascading nonlinearity
NASA Astrophysics Data System (ADS)
Zeng, Xianglong; Ashihara, Satoshi; Shimura, Tsutomu; Kuroda, Kazuo
2008-01-01
We experimentally demonstrate that adiabatic compression of femtosecond pulse can be achieved by employing the management of quadratic cascading nonlinearity in quasi-phase-matching gratings. Cascading nonlinearity is not a simple analogy with third-order optical nonlinearity in term of the engineering properties of the magnitude and focusing (or defocusing) nonlinearity. Femtosecond pulse compression is investigated based on type-I (e: o + o) collinear QPM geometry of aperiodically poled MgO-doped LiNbO 3 (MgO: LN). Group-velocity-matching condition is chosen to generate quadratic femtosecond soliton consisting of fundamental (FF) and second harmonic (SH) pulses. Adiabatic-like compression process is observed in the length of 50 mm linearly chirped QPM. Cascading nonlinearity is local managed, instead of dispersion management used in fiber adiabatic soliton compression. Quadratic soliton including FF and SH pulses are obtained from the compression of 95 fs FF pulse in the initial experiments. Dependence on the phase mismatch and group velocity mismatch, cascading nonlinearity has a flexible property and presents a new challenge for exploring femtosecond pulse shaping and control. The demonstrated pulse compression and control based on cascading nonlinearity is useful for generation of shorter pulses with clean temporal profiles, efficient femtosecond second harmonic generation and group-velocity control.
Finite-element analysis of earing using non-quadratic yield surfaces
Logan, R.W.
1995-06-18
During deep draw cupping, the phenomenon known as earing may occur as the cup wall is formed, resulting in a periodic variation of cup wall height around the perimeter of the finished cup. This is generally due to planar anisotropy of flow in rolled sheet product. It is generally observed that the anisotropy parameter R will vary in the plane of the sheet when ears are observed in cupping, with a parameter {Delta}R describing the variation of R in the plane of the sheet. For many common textures in face-centered and body-centered materials, the ears form relative to the sheet rolling direction at 0{degrees} and 90{degrees} around the perimeter if {Delta}R>0, and at -45{degrees} and +45{degrees} if {Delta}R<0. There is extensive experimental evidence that ear height shows a linear correlation with {Delta}R/R, but attempts to duplicate this using the finite-element method are highly dependent on both the methodology and yield surface used. It was shown previously that using a coarse mesh and the quadratic Hill yield surface tends to greatly under predict earing. In this study, we have used two different finite-element codes developed at LLNL to examine the predicted earing using both quadratic Hill and alternative non-quadratic yield surfaces. These results are compared to experimental data and conclusions drawn about the most desirable closed-form yield surfaces to duplicate the observed earing phenomena.
Quadratic function approaching method for magnetotelluric soundingdata inversion
Liangjun, Yan; Wenbao, Hu; Zhang, Keni
2004-04-05
The quadratic function approaching method (QFAM) is introduced for magnetotelluric sounding (MT) data inversion. The method takes the advantage of that quadratic function has single extreme value, which avoids leading to an inversion solution for local minimum and ensures the solution for global minimization of an objective function. The method does not need calculation of sensitivity matrix and not require a strict initial earth model. Examples for synthetic data and field measurement data indicate that the proposed inversion method is effective.
Gravity waves from non-minimal quadratic inflation
Pallis, Constantinos; Shafi, Qaisar
2015-03-12
We discuss non-minimal quadratic inflation in supersymmetric (SUSY) and non-SUSY models which entails a linear coupling of the inflaton to gravity. Imposing a lower bound on the parameter c{sub R}, involved in the coupling between the inflaton and the Ricci scalar curvature, inflation can be attained even for subplanckian values of the inflaton while the corresponding effective theory respects the perturbative unitarity up to the Planck scale. Working in the non-SUSY context we also consider radiative corrections to the inflationary potential due to a possible coupling of the inflaton to bosons or fermions. We find ranges of the parameters, depending mildly on the renormalization scale, with adjustable values of the spectral index n{sub s}, tensor-to-scalar ratio r≃(2−4)⋅10{sup −3}, and an inflaton mass close to 3⋅10{sup 13} GeV. In the SUSY framework we employ two gauge singlet chiral superfields, a logarithmic Kähler potential including all the allowed terms up to fourth order in powers of the various fields, and determine uniquely the superpotential by applying a continuous R and a global U(1) symmetry. When the Kähler manifold exhibits a no-scale-type symmetry, the model predicts n{sub s}≃0.963 and r≃0.004. Beyond no-scale SUGRA, n{sub s} and r depend crucially on the coefficient involved in the fourth order term, which mixes the inflaton with the accompanying non-inflaton field in the Kähler potential, and the prefactor encountered in it. Increasing slightly the latter above (−3), an efficient enhancement of the resulting r can be achieved putting it in the observable range. The inflaton mass in the last case is confined in the range (5−9)⋅10{sup 13} GeV.
NASA Astrophysics Data System (ADS)
Hong, Xiaochun; Xie, Shaolong; Chen, Longwei
In this study, we determine the associated number of zeros for Abelian integrals in four classes of quadratic reversible centers of genus one. Based on the results of [Li et al., 2002b],, we prove that the upper bounds of the associated number of zeros for Abelian integrals with orbits formed by conics, cubics, quartics, and sextics, under polynomial perturbations of arbitrary degree n, depend linearly on n.
Design of a candidate flutter suppression control law for DAST ARW-2
NASA Technical Reports Server (NTRS)
Adams, W. M., Jr.; Tiffany, S. H.
1984-01-01
A control law is developed to suppress symmetric flutter for a mathematical model of an aeroelastic research vehicle. An implementable control law is attained by including modified LQC (Linear Quadratic Gaussian) design techniques, controller order reduction, and gain scheduling. An alternate (complementary) design approach is illustrated for one flight condition wherein nongradient-based constrained optimization techniques are applied to maximize controller robustness.