Numerical Analysis and Computation Commons^™

Generalized Least-Squares Regressions I: Efficient Derivations, Nataniel Greene

Publications and Research

Ordinary least-squares regression suffers from a fundamental lack of symmetry: the regression line of y given x and the regression line of x given y are not inverses of each other. Alternative symmetric regression methods have been developed to address this concern, notably: orthogonal regression and geometric mean regression. This paper presents in detail a variety of least squares regression methods which may not have been known or fully explicated. The derivation of each method is made efficient through the use of Ehrenberg's formula for the ordinary least-squares error and through the extraction of a weight function g(b) which characterizes …

Pressure Poisson Method For The Incompressible Navier-Stokes Equations Using Galerkin Finite Elements, John Cornthwaite

College of Graduate Studies: Theses & Dissertations

In this thesis we examine the Navier-Stokes equations (NSE) with the continuity equation replaced by a pressure Poisson equation (PPE). Appropriate boundary conditions are developed for the PPE, which allow for a fully decoupled numerical scheme to recover the pressure. The variational form of the NSE with PPE is derived and used in the Galerkin Finite Element discretization. The Galerkin finite element method is then used to solve the NSE with PPE. Moderate accuracy is shown.

Full Newton Step Interior Point Method For Linear Complementarity Problem Over Symmetric Cones, Andrii Berdnikov

College of Graduate Studies: Theses & Dissertations

In this thesis, we present a new Feasible Interior-Point Method (IPM) for Linear Complementarity Problem (LPC) over Symmetric Cones. The advantage of this method lies in that it uses full Newton-steps, thus, avoiding the calculation of the step size at each iteration. By suitable choice of parameters we prove the global convergence of iterates which always stay in the the central path neighborhood. A global convergence of the method is proved and an upper bound for the number of iterations necessary to find ε-approximate solution of the problem is presented.

Nabla Fractional Calculus And Its Application In Analyzing Tumor Growth Of Cancer, Fang Wu

Masters Theses & Specialist Projects

This thesis consists of six chapters. In the first chapter, we review some basic definitions and concepts of fractional calculus. Then we introduce fractional difference equations involving the Riemann-Liouville operator of real number order between zero and one. In the second chapter, we apply the Brouwer fixed point and Contraction Mapping Theorems to prove that there exists a solution for up to the first order nabla fractional difference equation with an initial condition. In chapter three, we define a lower and an upper solution for up to the first order nabla fractional difference equation with an initial condition. Under certain …

A Constructive Proof Of Fundamental Theory For Fuzzy Variable Linear Programming Problems, A. Ebrahimnejad

Applications and Applied Mathematics: An International Journal (AAM)

Two existing methods for solving fuzzy variable linear programming problems based on ranking functions are the fuzzy primal simplex method proposed by Mahdavi-Amiri et al. (2009) and the fuzzy dual simplex method proposed by Mahdavi-Amiri and Nasseri (2007). In this paper, we prove that in the absence of degeneracy these fuzzy methods stop in a finite number of iterations. Moreover, we generalize the fundamental theorem of linear programming in a crisp environment to a fuzzy one. Finally, we illustrate our proof using a numerical example.

The Octonions And The Exceptional Lie Algebra G2, Ian M. Anderson

Research Vignettes

The octonions O are an 8-dimensional non-commutative, non-associative normed real algebra. The set of all derivations of O form a real Lie algebra. It is remarkable fact, first proved by E. Cartan in 1908, that the the derivation algebra of O is the compact form of the exceptional Lie algebra G2. In this worksheet we shall verify this result of Cartan and also show that the derivation algebra of the split octonions is the split real form of G2.

PDF and Maple worksheets can be downloaded from the links below.

Modeling And Mathematical Analysis Of Plant Models In Ecology, Eric A. Eager

Department of Mathematics: Dissertations, Theses, and Student Research

Population dynamics tries to explain in a simple mechanistic way the variations of the size and structure of biological populations. In this dissertation we use mathematical modeling and analysis to study the various aspects of the dynamics of plant populations and their seed banks.

In Chapter 2 we investigate the impact of structural model uncertainty by considering different nonlinear recruitment functions in an integral projection model for Cirsium canescens. We show that, while having identical equilibrium populations, these two models can elicit drastically different transient dynamics. We then derive a formula for the sensitivity of the equilibrium population to …

Degree Constrained Triangulation, Roshan Gyawali

UNLV Theses, Dissertations, Professional Papers, and Capstones

Triangulation of simple polygons or sets of points in two dimensions is a widely investigated problem in computational geometry. Some researchers have considered variations of triangulation problems that include minimum weight triangulation, de-launay triangulation and triangulation refinement. In this thesis we consider a constrained version of the triangulation problem that asks for triangulating a given domain (polygon or point sites) so that the resulting triangulation has an increased number of even degree vertices. This problem is called Degree Constrained Triangulation (DCT). We propose four algorithms to solve DCT problems. We also present experimental results based on the implementation of the …

Applying Differential Transform Method To Nonlinear Partial Differential Equations: A Modified Approach, Marwan T. Alquran

Applications and Applied Mathematics: An International Journal (AAM)

This paper proposes another use of the Differential transform method (DTM) in obtaining approximate solutions to nonlinear partial differential equations (PDEs). The idea here is that a PDE can be converted to an ordinary differential equation (ODE) upon using a wave variable, then applying the DTM to the resulting ODE. Three equations, namely, Benjamin-Bona-Mahony (BBM), Cahn-Hilliard equation and Gardner equation are considered in this study. The proposed method reduces the size of the numerical computations and use less rules than the usual DTM method used for multi-dimensional PDEs. The results show that this new approach gives very accurate solutions.

Geometric Programming Subject To System Of Fuzzy Relation Inequalities, Elyas Shivanian, Mahdi Keshtkar, Esmaile Khorram

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, an optimization model with geometric objective function is presented. Geometric programming is widely used; many objective functions in optimization problems can be analyzed by geometric programming. We often encounter these in resource allocation and structure optimization and technology management, etc. On the other hand, fuzzy relation equalities and inequalities are also used in many areas. We here present a geometric programming model with a monomial objective function subject to the fuzzy relation inequality constraints with maxproduct composition. Simplification operations have been given to accelerate the resolution of the problem by removing the components having no effect on …

A New Four Point Circular-Invariant Corner-Cutting Subdivision For Curve Design, Jian-Ao Lian

Applications and Applied Mathematics: An International Journal (AAM)

A 4-point nonlinear corner-cutting subdivision scheme is established. It is induced from a special C-shaped biarc circular spline structure. The scheme is circular-invariant and can be effectively applied to 2-dimensional (2D) data sets that are locally convex. The scheme is also extended adaptively to non-convex data. Explicit examples are demonstrated.

Numerical Methods For Fluid-Structure Interaction - A Review, Gene Hou, Jin Wang, Anita Layton

Mechanical & Aerospace Engineering Faculty Publications

The interactions between incompressible fluid flows and immersed structures are nonlinear multi-physics phenomena that have applications to a wide range of scientific and engineering disciplines. In this article, we review representative numerical-methods based on conforming and non-conforming meshes that are currently available for computing fluid-structure interaction problems, with an emphasis on some of the recent developments in the field. A goal is to categorize the selected methods and assess their accuracy and efficiency. We discuss challenges faced by researchers in this field, and we emphasize the importance of interdisciplinary effort for advancing the study in fluid-structure interactions

Robust 𝒍𝟏 And 𝒍∞ Solutions Of Linear Inequalities, Maziar Salahi

Applications and Applied Mathematics: An International Journal (AAM)

Infeasible linear inequalities appear in many disciplines. In this paper we investigate the 𝑙₁ and 𝑙_∞ solutions of such systems in the presence of uncertainties in the problem data. We give equivalent linear programming formulations for the robust problems. Finally, several illustrative numerical examples using the cvx software package are solved showing the importance of the robust model in the presence of uncertainties in the problem data.

Improved Algorithms For Ear-Clipping Triangulation, Bartosz Kajak

UNLV Theses, Dissertations, Professional Papers, and Capstones

We consider the problem of improving ear-slicing algorithm for triangulating a simple polygon. We propose two variations of ear-slicing technique for generating “good-quality” triangulation. The first approach is based on searching for the best triangle along the boundary. The second approach considers polygon partitioning on a pre-process before applying the ear-slicing. Experimental investigation reveals that both approaches yield better quality triangulation than the standard ear-slicing method.

Reliability Analysis Of A Series And Parallel Network Using Triangular Intuitionistic Fuzzy Sets, D. Pandey, S. K. Tyagi, Vinesh Kumar

Applications and Applied Mathematics: An International Journal (AAM)

This paper describes a novel approach, based on intuitionistic fuzzy set theory for reliability analysis of series and parallel network. The triangular intuitionistic fuzzy sets are used to represent the failure possibility of each basic (terminal) event to get more comprehensive results for the failure possibility of the top event. The proposed technique is demonstrated on a web server LOG data used to illustrate HTTP (Hyper Text Transfer Protocol) failure

Exact Optimal Solution Of Fuzzy Critical Path Problems, Amit Kumar, Parmpreet Kumar

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, a fuzzy critical path problem is chosen to show that the results, obtained by using the existing method [Liu, S.T.: Fuzzy activity times in critical path and project crashing problems. Cybernetics and Systems 34 (2), 161-172 (2003)], could be improved to reflect, more appropriate real life situations. To obtain more accurate results of fuzzy critical path problems, a new method that modifies the existing one is proposed here. To demonstrate the advantages of the proposed method it is used to solve a specific fuzzy critical path problem.

Shooting Neural Networks Algorithm For Solving Boundary Value Problems In Odes, Kais I. Ibraheem, Bashir M. Khalaf

Applications and Applied Mathematics: An International Journal (AAM)

The objective of this paper is to use Neural Networks for solving boundary value problems (BVPs) in Ordinary Differential Equations (ODEs). The Neural networks use the principle of Back propagation. Five examples are considered to show effectiveness of using the shooting techniques and neural network for solving the BVPs in ODEs. The convergence properties of the technique, which depend on the convergence of the integration technique and accuracy of the interpolation technique are considered.

A Nonnegative Analysis Of Politics, Tim Chartier, Charles D. Wessell

Math Faculty Publications

The article investigates how linear algebra can recover mathematical information from the electronic messages using the Enron Email Sets in Pennsylvania. It states that a term-by-email matrix has been created to cluster algorithm, which allows one to mine through data and discover meaning. Moreover, nonnegative matrix factorization (NMF) enables one to interpret the resulting factorization in terms of the original problem.

Quantitative Stability Of Linear Infinite Inequality Systems Under Block Perturbations With Applications To Convex Systems, M J. Cánovas, M A. Lopez, Boris S. Mordukhovich, J Parra

Mathematics Research Reports

The original motivation for this paper was to provide an efficient quantitative analysis of convex infinite (or semi-infinite) inequality systems whose decision variables run over general infinite-dimensional (resp. finite-dimensional) Banach spaces and that are indexed by an arbitrary fixed set J. Parameter perturbations on the right-hand side of the inequalities are required to be merely bounded, and thus the natural parameter space is loo(J). Our basic strategy consists of linearizing the parameterized convex system via splitting convex inequalities into linear ones by using the Fenchel-Legendre conjugate. This approach yields that arbitrary bounded right-hand side perturbations of the convex system turn …

Complete Characterizations Of Local Weak Sharp Minima With Applications To Semi-Infinite Optimization And Complementarity, Boris S. Mordukhovich, Naihua Xiu, Jinchuan Zhou

Mathematics Research Reports

In this paper we identify a favorable class of nonsmooth functions for which local weak sharp minima can be completely characterized in terms of normal cones and subdifferentials, or tangent cones and subderivatives, or their mixture in finite-dimensional spaces. The results obtained not only significantly extend previous ones in the literature, but also allow us to provide new types of criteria for local weak sharpness. Applications of the developed theory are given to semi-infinite programming and to semi-infinite complementarity problems.

A New Method For Fuzzy Critical Path Analysis In Project Networks With A New Representation Of Triangular Fuzzy Numbers, Amit Kumar, Parmpreet Kaur

Applications and Applied Mathematics: An International Journal (AAM)

The method for finding fuzzy optimal solution of fully fuzzy critical path (FFCP) problems i.e., critical path problems in which all the parameters are represented by fuzzy numbers, is at best scant; possibly non-existent. In this paper, a method is proposed to find the fuzzy optimal solution of FFCP problems, together with a new representation of triangular fuzzy numbers. This paper will show the advantages of using, the proposed representation over the existing representations of triangular fuzzy numbers and will present with great clarity the proposed method and illustrate its application to FFCP problems occurring in real life situations.

An Analytical Technique For Solving Nonlinear Heat Transfer Equations, Hossein Aminikhah, Milad Hemmatnezhad

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, an analytic technique, namely the New Homotopy Perturbation Method (NHPM) is applied for solving the nonlinear differential equations arising in the field of heat transfer. In this method, the solution is considered as an infinite series expansion where converges rapidly to the exact solution. The nonlinear convective–radioactive cooling equation and nonlinear equation of conduction heat transfer with the variable physical properties are chosen as illustrative examples and the exact solutions have been found for each case.

Duality In Fuzzy Linear Programming With Symmetric Trapezoidal Numbers, S. H. Nasseri, E. Ebrahimnejad, S. Mizuno

Applications and Applied Mathematics: An International Journal (AAM)

Linear programming problems with trapezoidal fuzzy numbers have recently attracted much interest. Various methods have been developed for solving these types of problems. Here, following the work of Ganesan and Veeramani and using the recent approach of Mahdavi-Amiri and Nasseri, we introduce the dual of the linear programming problem with symmetric trapezoidal fuzzy numbers and establish some duality results. The results will be useful for post optimality analysis.

Optimal Correction Of Infeasible System In Linear Equality Via Genetic Algorithm, S. Ketabchi, H. Moosaei, S. Fallahi

Applications and Applied Mathematics: An International Journal (AAM)

This work is focused on the optimal correction of infeasible system of linear equality. In this paper, for correcting this system, we will make the changes just in the coefficient matrix by using l 􀬶 norm and show that solving this problem is equivalent to solving a fractional quadratic problem. To solve this problem, we use the genetic algorithm. Some examples are provided to illustrate the efficiency and validity of the proposed method.

Approximate Analytical Solutions For Fractional Space- And Time- Partial Differential Equations Using Homotopy Analysis Method, Subir, Das, R. Kumar, P. K. Gupta, Hossein Jafari

Applications and Applied Mathematics: An International Journal (AAM)

This article presents the approximate analytical solutions of first order linear partial differential equations (PDEs) with fractional time- and space- derivatives. With the aid of initial values, the explicit solutions of the equations are solved making use of reliable algorithm like homotopy analysis method (HAM). The speed of convergence of the method is based on a rapidly convergent series with easily computable components. The fractional derivatives are described in Caputo sense. Numerical results show that the HAM is easy to implement and accurate when applied to space- time- fractional PDEs.

Generalized Newton's Method Based On Graphical Derivatives, T Hoheisel, C Kanzow, Boris S. Mordukhovich, Hung M. Phan

Mathematics Research Reports

This paper concerns developing a numerical method of the Newton type to solve systems of nonlinear equations described by nonsmooth continuous functions. We propose and justify a new generalized Newton algorithm based on graphical derivatives, which have never been used to derive a Newton-type method for solving nonsmooth equations. Based on advanced techniques of variational analysis and generalized differentiation, we establish the well-posedness of the algorithm, its local superlinear convergence, and its global convergence of the Kantorovich type. Our convergence results hold with no semismoothness assumption, which is illustrated by examples. The algorithm and main results obtained in the paper …

Wavelet Transform Of Fractional Integrals For Integrable Boehmians, Deshna Loonker, P. K. Banerji, S. L. Kalla

Applications and Applied Mathematics: An International Journal (AAM)

The present paper deals with the wavelet transform of fractional integral operator (the Riemann- Liouville operators) on Boehmian spaces. By virtue of the existing relation between the wavelet transform and the Fourier transform, we obtained integrable Boehmians defined on the Boehmian space for the wavelet transform of fractional integrals.

Convergence Of The Sinc Method Applied To Volterra Integral Equations, M. Zarebnia, J. Rashidinia

Applications and Applied Mathematics: An International Journal (AAM)

A collocation procedure is developed for the linear and nonlinear Volterra integral equations, using the globally defined Sinc and auxiliary basis functions. We analytically show the exponential convergence of the Sinc collocation method for approximate solution of Volterra integral equations. Numerical examples are included to confirm applicability and justify rapid convergence of our method.

An Algorithm To Generate Two-Dimensional Drawings Of Conway Algebraic Knots, Jen-Fu Tung

Masters Theses & Specialist Projects

The problem of finding an efficient algorithm to create a two-dimensional embedding of a knot diagram is not an easy one. Typically, knots with a large number of crossings will not nicely generate two-dimensional drawings. This thesis presents an efficient algorithm to generate a knot and to create a nice two-dimensional embedding of the knot. For the purpose of this thesis a drawing is “nice” if the number of tangles in the diagram consisting of half-twists is minimal. More specifically, the algorithm generates prime, alternating Conway algebraic knots in O(n) time where n is the number of crossings …

Multiresolution Inverse Wavelet Reconstruction From A Fourier Partial Sum, Nataniel Greene

Publications and Research

The Gibbs phenomenon refers to the lack of uniform convergence which occurs in many orthogonal basis approximations to piecewise smooth functions. This lack of uniform convergence manifests itself in spurious oscillations near the points of discontinuity and a low order of convergence away from the discontinuities.In previous work [11,12] we described a numerical procedure for overcoming the Gibbs phenomenon called the Inverse Wavelet Reconstruction method (IWR). The method takes the Fourier coefficients of an oscillatory partial sum and uses them to construct the wavelet coefficients of a non-oscillatory wavelet series. However, we only described the method standard wavelet series and …