Two new algorithms are presented in the paper for the plus/minus factorization of a scalar discrete-time polynomial. The first method is based on the discrete Fourier transform theory (DFT) and its relationship to the...
详细信息
Two new algorithms are presented in the paper for the plus/minus factorization of a scalar discrete-time polynomial. The first method is based on the discrete Fourier transform theory (DFT) and its relationship to the Z-transform. Involving DFT computational techniques and the famous fast Fourier transform routine brings high computational efficiency and reliability. The method is applied in the case study of H-2-optimal inverse dynamic filter to an audio equipment. The second numerical procedure originates in a symmetric spectral factorization routine, namely the Bauer's method of the 1950s. As a by-product, a recursive LU factorization procedure for Toeplitz matrices is devised that is of more general impact and can be of use in other areas of applied mathematics as well. Performance of the method is demonstrated by an l(1) optimal controller design example. Copyright (c) 2006 John Wiley & Sons, Ltd.
We present numerical algorithms calculating compositions of the left and right fractional integrals. We apply quadratic interpolation and obtain the fractional Simpson's rule. We estimate the local truncation erro...
详细信息
We present numerical algorithms calculating compositions of the left and right fractional integrals. We apply quadratic interpolation and obtain the fractional Simpson's rule. We estimate the local truncation error of the proposed approximations, calculate errors generated by presented algorithms, and determine the experimental rate of convergence. Finally, we show examples of numerical evaluation of these operators.
We present several numerical algorithms for six-degree-of-freedom rigid-body registration of line fiducial objects to their marks in cross-sectional planar images, such as those obtained in CT and MRI, given the corre...
详细信息
We present several numerical algorithms for six-degree-of-freedom rigid-body registration of line fiducial objects to their marks in cross-sectional planar images, such as those obtained in CT and MRI, given the correspondence between the marks and line fiducials. The area of immediate application is frame-based stereotactic procedures, such as radiosurgery and functional neurosurgery. The algorithms are also suitable to problems where the fiducial pattern moves inside the imager, as is the case in robot-assisted image-guided surgical applications. We demonstrate the numerical methods on clinical CT images and computer-generated data and compare their performance in terms of robustness to missing data, robustness to noise, and speed. The methods show two unique strengths: (1) They provide reliable registration of incomplete fiducial patterns when up to two-thirds of the total fiducials are missing from the image;and (2) they are applicable to an arbitrary combination of line fiducials without algorithmic modification. The average speed of the fastest algorithm is 0.3236 s for six fiducial lines in real CT data in a Matlab implementation. (C) 2002 American Association of Physicists in Medicine.
We consider three integral observations inverse problems (IP1, IP2, IP3) for reconstruction of a time-dependent convection coefficient and a source in magnetohydrodynamics (MHD) model. On the first stage, using the in...
详细信息
We consider three integral observations inverse problems (IP1, IP2, IP3) for reconstruction of a time-dependent convection coefficient and a source in magnetohydrodynamics (MHD) model. On the first stage, using the integral observations, we reduce the inverse problems to nonclassical direct (forward) ones. The equivalence between the inverse and direct problems is established. Then, the well-posedness of nonclassical problems is proved. Further, to overcome the difficulties arising from the nonlinear nonlocal parabolic operators, we construct a linarization algorithm in time after their difference space discretizatons. Next, on each iteration, to solve the corresponding linear systems of algebraic equations, we propose adequate fast elimination algorithms. Computational results of test examples data are discussed.
The discrete-time optimal projection equations, which constitute necessary conditions fir optimal reduced-order LQG compensation, are strengthened. For the class of minimal stabilizing compensators the strengthened di...
详细信息
The discrete-time optimal projection equations, which constitute necessary conditions fir optimal reduced-order LQG compensation, are strengthened. For the class of minimal stabilizing compensators the strengthened discrete-time optimal projection equations are proved to he equivalent to first-order necessary optimality conditions for optimal reduced-order LQG compensation. The conventional discrete-time optimal projection equations are pro red to be weaker. As a result solutions of the conventional discrete-time optimal projection equations may not correspond to optimal reduced-order compensators. Through numerical examples it is demonstrated that, in fact, many solutions exist that do not correspond to optimal reduced-order compensators. To compute optimal reduced-order compensators two licit, algorithms are proposed. One is a homotopy algorithm and one is based oil iteration of the strengthened discrete-time optimal projection equations. The latter algorithm is a generalization of the algorithm that solves the two Riccati equations of full-order LQG control through iteration and therefore is highly efficient. Using different initializations of the iterative algorithm it is demonstrated that the reduced-order optimal LQG compensation problem, in general, may possess multiple extrema. Through two computer experiments it is demonstrated that the homotopy algorithm often, but not always, finds the global minimum.
Difficulties associated with ultrahigh-performance computer systems that will appear in the near future and possible ways of their solution are discussed. Examples of simulating magnetogasdynamics problems are given.
Difficulties associated with ultrahigh-performance computer systems that will appear in the near future and possible ways of their solution are discussed. Examples of simulating magnetogasdynamics problems are given.
The complexity and specificity of stability analysis applied to interpolative-type control structures makes the study of such property a difficult, almost impossible task - at least in the analytical manner. Therefore...
详细信息
ISBN:
(纸本)9781509017355
The complexity and specificity of stability analysis applied to interpolative-type control structures makes the study of such property a difficult, almost impossible task - at least in the analytical manner. Therefore, a better solution could be represented by the numerical approach. In this context, this paper starts with developing some methods and techniques with applicability for analysis of the interpolative-type controllers, based on Lyapunov method perspective. These methodological aspects are gathered together into a specific procedural algorithm. Based on this algorithm, a set of MATLAB-SIMULINK programs able to offer, in a flexible and user-interactive way, a possible solution to Lyapunov-stability analysis for a class of interpolative-type control systems with linear or non-linear processes of 2nd and 3rd order are developed. The solution based on the implemented software packages was finally validated through some practical examples.
The discrete-time optimal projection equations for systems with white parameters are strengthened. For the class of minimal ms (mean square) stabilizing compensators the strengthened discrete-time optimal projection e...
详细信息
ISBN:
(纸本)085296708X
The discrete-time optimal projection equations for systems with white parameters are strengthened. For the class of minimal ms (mean square) stabilizing compensators the strengthened discrete-time optimal projection equations are proved to be equivalent to first-order necessary optimality conditions for optimal reduced-order dynamic compensation of systems with white parameters. The conventional discrete-time optimal projection equations are proved to be weaker. As a result solutions of the conventional discrete-time optimal projection equations may not correspond to optimal reduced-order compensators. To compute optimal reduced-order compensators two numerical algorithms are proposed. One is a homotopy algorithm and one is based on iteration of the strengthened discrete-time optimal projection equations. The latter algorithm is a generalization of the algorithm that solves the full-order problem, which in turn is a generalization of the algorithm that solves the two Riccati equations of full-order LQG control through iteration. Therefore the efficiency of these three types of algorithms is comparable. It is demonstrated that, despite the strengthening of the optimal projection equations, the optimal reduced-order compensation problem, in general, may posses multiple extrema.
The Kantorovich analysis (Argyros in Convergence and applications of Newton-type iterations, Springer, New York, 2008;Argyros and Hi lout in Efficient methods for solving equations and variational inequalities, Polime...
详细信息
The Kantorovich analysis (Argyros in Convergence and applications of Newton-type iterations, Springer, New York, 2008;Argyros and Hi lout in Efficient methods for solving equations and variational inequalities, Polimetrica Publisher, Milano, 2009;Kantorovich and Akilov in Functional analysis, Pergamon Press, Oxford, 1982), and recurrent relation's approach (Gutierrez et al. in J Comput Appl Math 115:181-192,2000) are the most popular ways for generating sufficient conditions for the convergence of numerical algorithms to a solution of a nonlinear equations as well as providing the corresponding error estimates on the distances involved. We introduce the new approach of recurrent functions to show that a finer convergence analysis can be provided under the same hypotheses, and computational cost. numerical examples are provided where our results apply, hut not earlier ones.
This paper describes TANOR, an automated framework for designing hardware accelerators for numerical computation on reconfigurable platforms. Applications utilizing numerical algorithms on large-size data sets require...
详细信息
This paper describes TANOR, an automated framework for designing hardware accelerators for numerical computation on reconfigurable platforms. Applications utilizing numerical algorithms on large-size data sets require high-throughput computation platforms. The focus is on N-body interaction problems which have a wide range of applications spanning from astrophysics to molecular dynamics. The TANOR design flow starts with a MATLAB description of a particular interaction function, its parameters, and certain architectural constraints specified through a graphical user interface. Subsequently, TANOR automatically generates a configuration bitstream for a target FPGA along with associated drivers and control software necessary to direct the application from a host PC. Architectural exploration is facilitated through support for fully custom fixed-point and floating-point representations in addition to standard number representations such as single-precision floating point. Moreover, TANOR enables joint exploration of algorithmic and architectural variations in realizing efficient hardware accelerators. TANOR's capabilities have been demonstrated for three different N-body interaction applications: the calculation of gravitational potential in astrophysics, the diffusion or convolution with Gaussian kernel common in image processing applications, and the force calculation with vector-valued kernel function in molecular dynamics simulation. Experimental results show that TANOR-generated hardware accelerators achieve lower resource utilization without compromising numerical accuracy, in comparison to other existing custom accelerators.
暂无评论