SMART: a program for simplified matrix algebra on remote terminals.

Thesis (M.S.)--University of Illinois.

Central elements of the elliptic Zn monodromy matrix algebra at roots of unity

The central elements of the algebra of monodromy matrices associated with the Z(n) R-matrix are studied. When the crossing parameter w takes a special rational value w = n/N, where N and n are positive coprime integers, the center is substantially larger than that in the generic case for which the quantum determinant provides the center. In the trigonometric limit, the situation corresponds to the quantum group at roots of unity. This is a higher rank generalization of the recent results by Belavin and Jimbo. (c) 2004 Elsevier B.V. All rights reserved.

Learning the kernel matrix with semidefinite programming

Kernel-based learning algorithms work by embedding the data into a Euclidean space, and then searching for linear relations among the embedded data points. The embedding is performed implicitly, by specifying the inner products between each pair of points in the embedding space. This information is contained in the so-called kernel matrix, a symmetric and positive semidefinite matrix that encodes the relative positions of all points. Specifying this matrix amounts to specifying the geometry of the embedding space and inducing a notion of similarity in the input space - classical model selection problems in machine learning. In this paper we show how the kernel matrix can be learned from data via semidefinite programming (SDP) techniques. When applied to a kernel matrix associated with both training and test data this gives a powerful transductive algorithm -using the labeled part of the data one can learn an embedding also for the unlabeled part. The similarity between test points is inferred from training points and their labels. Importantly, these learning problems are convex, so we obtain a method for learning both the model class and the function without local minima. Furthermore, this approach leads directly to a convex method for learning the 2-norm soft margin parameter in support vector machines, solving an important open problem.

A simple demonstration of the existence of the jordan canonical form for any square matrix

All the demonstrations known to this author of the existence of the Jordan Canonical Form are somewhat complex - usually invoking the use of new spaces, and what not. These demonstrations are usually too difficult for an average Mathematics student to understand how he or she can obtain the Jordan Canonical Form for any square matrix. The method here proposed not only demonstrates the existence of such forms but, additionally, shows how to find them in a step by step manner. I do not claim that the following demonstration is in any way “elegant” (by the standards of elegance in fashion nowadays among mathematicians) but merely simple (undergraduate students taking a fist course in Matrix Algebra would understand how it works).

Robust control to parametric uncertainties in smart structures using linear matrix inequalities

The study of algorithms for active vibrations control in flexible structures became an area of enormous interest, mainly due to the countless demands of an optimal performance of mechanical systems as aircraft and aerospace structures. Smart structures, formed by a structure base, coupled with piezoelectric actuators and sensor are capable to guarantee the conditions demanded through the application of several types of controllers. This article shows some steps that should be followed in the design of a smart structure. It is discussed: the optimal placement of actuators, the model reduction and the controller design through techniques involving linear matrix inequalities (LMI). It is considered as constraints in LMI: the decay rate, voltage input limitation in the actuators and bounded output peak (output energy). Two controllers robust to parametric variation are designed: the first one considers the actuator in non-optimal location and the second one the actuator is put in an optimal placement. The performance are compared and discussed. The simulations to illustrate the methodology are made with a cantilever beam with bonded piezoelectric actuators.

Eigenvalue analyses of two parallel lines using a single real transformation matrix

For a typical non-symmetrical system with two parallel three phase transmission lines, modal transformation is applied using some examples of single real transformation matrices. These examples are applied searching an adequate single real transformation matrix to two parallel three phase transmission line systems. The analyses are started with the eigenvector and eigenvalue studies, using Clarke's transformation or linear combinations of Clarke's elements. The Z C and parameters are analyzed for the case that presents the smallest errors between the exact eigenvalues and the single real transformation matrix application results. The single real transformation determined for this case is based on Clarke's matrix and its main characteristic is the use of a unique homopolar reference. So, the homopolar mode becomes a connector mode between the two three-phase circuits of the analyzed system. ©2005 IEEE.

Non-transposed three-phase line analyses with a single real transformation matrix

In transmission line transient analyses, a single real transformation matrix can obtain exact modes when the analyzed line is transposed. For non-transposed lines, the results are not exact. In this paper, non-symmetrical and non transposed three-phase line samples are analyzed with a single real transformation matrix application (Clarke's matrix). Some interesting characteristics of this matrix application are: single, real, frequency independent, line parameter independent, identical for voltage and current determination. With Clarke's matrix use, mathematical simplifications are obtained and the developed model can be applied directly in programs based on time domain. This model works without convolution procedures to deal with phase-mode transformation. In EMTP programs, Clarke's matrix can be represented by ideal transformers and the frequency dependent line parameters can be represented by modified-circuits. With these representations, the electrical values at any line point can be accessed for phase domain or mode domain using the Clarke matrix or its inverse matrix. For symmetrical and non-transposed lines, the model originates quite small errors. In addition, the application of the proposed model to the non-symmetrical and non-transposed three phase transmission lines is investigated. ©2005 IEEE.

Active flutter suppression in a 2-D airfoil using linear matrix inequalities techniques

Flutter is an in-flight vibration of flexible structures caused by energy in the airstream absorbed by the lifting surface. This aeroelastic phenomenon is a problem of considerable interest in the aeronautic industry, because flutter is a potentially destructive instability resulting from an interaction between aerodynamic, inertial, and elastic forces. To overcome this effect, it is possible to use passive or active methodologies, but passive control adds mass to the structure and it is, therefore, undesirable. Thus, in this paper, the goal is to use linear matrix inequalities (LMIs) techniques to design an active state-feedback control to suppress flutter. Due to unmeasurable aerodynamic-lag states, one needs to use a dynamic observer. So, LMIs also were applied to design a state-estimator. The simulated model, consists of a classical flat plate in a two-dimensional flow. Two regulators were designed, the first one is a non-robust design for parametric variation and the second one is a robust control design, both designed by using LMIs. The parametric uncertainties are modeled through polytopic uncertainties. The paper concludes with numerical simulations for each controller. The open-loop and closed-loop responses are also compared and the results show the flutter suppression. The perfomance for both controllers are compared and discussed. Copyright © 2006 by ABCM.

Analysis of numeric conditioning of Hessian matrix of Lagrangian via singular values

This paper presents an analyze of numeric conditioning of the Hessian matrix of Lagrangian of modified barrier function Lagrangian method (MBFL) and primal-dual logarithmic barrier method (PDLB), which are obtained in the process of solution of an optimal power flow problem (OPF). This analyze is done by a comparative study through the singular values decomposition (SVD) of those matrixes. In the MBLF method the inequality constraints are treated by the modified barrier and PDLB methods. The inequality constraints are transformed into equalities by introducing positive auxiliary variables and are perturbed by the barrier parameter. The first-order necessary conditions of the Lagrangian function are solved by Newton's method. The perturbation of the auxiliary variables results in an expansion of the feasible set of the original problem, allowing the limits of the inequality constraints to be reached. The electric systems IEEE 14, 162 and 300 buses were used in the comparative analysis. ©2007 IEEE.

Modal transformation obtained from Clarke's matrix - Asymmetrical three-phase case

Some constant matrices can be used as phase-mode transformation matrices for transposed three-phase transmission lines. Clarke's matrix is one of these options. Its application as a phase-mode transformation matrix for untransposed three-phase transmission lines has been analyzed through error and frequency scan comparisons. Based on an actual untransposed asymmetrical three-phase transmission line example, a correction procedure is applied searching for better results from the Clarke's matrix applicaton as a phase-mode transformation matrix. The error analyses are carried out using Clarke's matrix and the new transformation matrices obtained from the correction procedure. Applying Clarke's matrix, the relative errors of the eigenvalue matrix elements can be considered negligible and the relative values of the off-diagonal elements are significant. If the the corrected transformation matrices are used, the relative values of the off-diagonal elements are decreased. Based on the results of these analyses, the homopolar mode is more sensitive to the frequency influence than the two other modes related to three-phase lines. © 2007 IEEE.

An enriched radial point interpolation method (e-RPIM) for analysis of crack tip fields

In this paper, an enriched radial point interpolation method (e-RPIM) is developed the for the determination of crack tip fields. In e-RPIM, the conventional RBF interpolation is novelly augmented by the suitable trigonometric basis functions to reflect the properties of stresses for the crack tip fields. The performance of the enriched RBF meshfree shape functions is firstly investigated to fit different surfaces. The surface fitting results have proven that, comparing with the conventional RBF shape function, the enriched RBF shape function has: (1) a similar accuracy to fit a polynomial surface; (2) a much better accuracy to fit a trigonometric surface; and (3) a similar interpolation stability without increase of the condition number of the RBF interpolation matrix. Therefore, it has proven that the enriched RBF shape function will not only possess all advantages of the conventional RBF shape function, but also can accurately reflect the properties of stresses for the crack tip fields. The system of equations for the crack analysis is then derived based on the enriched RBF meshfree shape function and the meshfree weak-form. Several problems of linear fracture mechanics are simulated using this newlydeveloped e-RPIM method. It has demonstrated that the present e-RPIM is very accurate and stable, and it has a good potential to develop a practical simulation tool for fracture mechanics problems.

Control of an underactuated-slender-hull unmanned underwater vehicle using Port-Hamiltonian theory

This paper presents a control design for tracking of attitude and speed of an underactuated slender-hull unmanned underwater vehicle (UUV). The control design is based on Port-Hamiltonian theory. The target dynamics (desired dynamic response) is shaped with particular attention to the target mass matrix so that the influence of the unactuated dynamics on the controlled system is suppressed. This results in achievable dynamics independent of uncontrolled states. Throughout the design, insight of the physical phenomena involved is used to propose the desired target dynamics. The performance of the design is demonstrated through simulation with a high-fidelity model.

Stationary random response of multidegree-of-freedom systems

An approximate approach is presented for determining the stationary random response of a general multidegree-of-freedom nonlinear system under stationary Gaussian excitation. This approach relies on defining an equivalent linear system for the nonlinear system. Two particular systems which possess exact solutions have been solved by this approach, and it is concluded that this approach can generate reasonable solutions even for systems with fairly large nonlinearities. The approximate approach has also been applied to two examples for which no exact or approximate solutions were previously available.

Also presented is a matrix algebra approach for determining the stationary random response of a general multidegree-of-freedom linear system. Its derivation involves only matrix algebra and some properties of the instantaneous correlation matricies of a stationary process. It is therefore very direct and straightforward. The application of this matrix algebra approach is in general simpler than that of commonly used approaches.

A probabilistic index for completely positive maps and an application

Gohm, Rolf, (2003) 'A probabilistic index for completely positive maps and an application', Journal of Operator Theory 54(2) pp.339-361 RAE2008