Applied Mathematics Commons^™

Noisy Signal Recovery Via Iterative Reweighted L1-Minimization, Deanna Needell

CMC Faculty Publications and Research

Compressed sensing has shown that it is possible to reconstruct sparse high dimensional signals from few linear measurements. In many cases, the solution can be obtained by solving an L1-minimization problem, and this method is accurate even in the presence of noise. Recent a modified version of this method, reweighted L1-minimization, has been suggested. Although no provable results have yet been attained, empirical studies have suggested the reweighted version outperforms the standard method. Here we analyze the reweighted L1-minimization method in the noisy case, and provide provable results showing an improvement in the error bound over the standard bounds.

Fractional Calculus: Definitions And Applications, Joseph M. Kimeu

Masters Theses & Specialist Projects

No abstract provided.

Determining The Orbit Locations Of Turkish Airborne Early Warning And Control Aircraft Over The Turkish Air Space, Nebi Sarikaya

Theses and Dissertations

The technology improvement affects the military needs of individual countries. The new doctrine of defense for many countries emphasizes detecting threats as far away as you can from your homeland. Today, the military uses both ground RADAR and Airborne Early Warning and Control (AEW&C) Aircraft. AEW&C aircraft has become vital to detect low altitude threats that a ground RADAR cannot detect because of obstacles on the earth. Turkey has ordered four AEW&C aircraft for her air defense system because of the lack of complete coverage by ground RADAR. This research provides optimal orbit locations that can be updated according to …

Efficient Evaluation Of Ranking Procedures When The Number Of Units Is Large With Application To Snp Identification, Thomas A. Louis, Ingo Ruczinski

Johns Hopkins University, Dept. of Biostatistics Working Papers

Simulation-based assessment is a popular and frequently necessary approach to evaluation of statistical procedures. Sometimes overlooked is the ability to take advantage of underlying mathematical relations and we focus on this aspect. We show how to take advantage of large-sample theory when conducting a simulation using the analysis of genomic data as a motivating example. The approach uses convergence results to provide an approximation to smaller-sample results, results that are available only by simulation. We consider evaluating and comparing a variety of ranking-based methods for identifying the most highly associated SNPs in a genome-wide association study, derive integral equation representations …

Thermocapillary Effects In Driven Dewetting And Self-Assembly Of Pulsed Laser-Irradiated Metallic Films, Agegnehu Atena, Mikhail Khenner

Mathematics Faculty Publications

In this paper the lubrication-type dynamical model is developed of a molten, pulsed laser-irradiated metallic film. The heat transfer problem that incorporates the absorbed heat from a single beam or interfering beams is solved analytically. Using this temperature field, we derive the 3D long-wave evolution PDE for the film height. To get insights into dynamics of dewetting, we study the 2D version of the evolution equation by means of a linear stability analysis and by numerical simulations. The stabilizing and destabilizing effects of various system parameters, such as the peak laser beam intensity, the film optical thickness, the Biot and …

Thermocapillary Effects In Driven Dewetting And Self-Assembly Of Pulsed Laser-Irradiated Metallic Films, Agegnehu Atena, Mikhail Khenner

Mathematics Faculty Publications

In this paper the lubrication-type dynamical model is developed of a molten, pulsed laser-irradiated metallic film. The heat transfer problem that incorporates the absorbed heat from a single beam or interfering beams is solved analytically. Using this temperature field, we derive the 3D long-wave evolution PDE for the film height. To get insights into dynamics of dewetting, we study the 2D version of the evolution equation by means of a linear stability analysis and by numerical simulations. The stabilizing and destabilizing effects of various system parameters, such as the peak laser beam intensity, the film optical thickness, the Biot and …

The Beauty Of Mathematics And The Mathematics Of Beauty: Continued Fractions And The Golden Ratio, Jessica Tush

Inquiry: The University of Arkansas Undergraduate Research Journal

This project begins with a look at the history of simple continued fractions and how we have arrived where we are today. We then move through a study of simple continued fractions, beginning first with rational numbers and moving to irrational numbers. Continuing further in the pursuit of joining mathematics and art, we define the specific continued fraction that gives rise to the Fibonacci sequence and the Golden Ratio~ (phi, pronounced 'Jai"). These two notions form a direct link to art and the properties that we hope to examine. I have taken an analytic approach to showing that the Golden …

An H-Adaptive Finite-Element Technique For Constructing 3d Wind Fields, Darrell Pepper, Xiuling Wang

Mechanical Engineering Faculty Research

An h-adaptive, mass-consistent finite-element model (FEM) has been developed for constructing 3D wind fields over irregular terrain utilizing sparse meteorological tower data. The element size in the computational domain is dynamically controlled by an a posteriori error estimator based on the L2 norm. In the h-adaptive FEM algorithm, large element sizes are typically associated with smooth flow regions and small errors; small element sizes are attributed to fast-changing flow regions and large errors. The adaptive procedure employed in this model uses mesh refinement–unrefinement to satisfy error criteria. Results are presented for wind fields using sparse data obtained from two regions …

Polynomial Extension Operators. Part Ii, Leszek Demkowicz, Jay Gopalakrishnan, Joachim Schöberl

Mathematics and Statistics Faculty Publications and Presentations

Consider the tangential trace of a vector polynomial on the surface of a tetrahedron. We construct an extension operator that extends such a trace function into a polynomial on the tetrahedron. This operator can be continuously extended to the trace space of H(curl ). Furthermore, it satisfies a commutativity property with an extension operator we constructed in Part I of this series. Such extensions are a fundamental ingredient of high order finite element analysis.

A Covert Encryption Method For Applications In Electronic Data Interchange, Jonathan Blackledge, Dmitry Dubovitskiy

Articles

A principal weakness of all encryption systems is that the output data can be ‘seen’ to be encrypted. In other words, encrypted data provides a ‘flag’ on the potential value of the information that has been encrypted. In this paper, we provide a new approach to ‘hiding’ encrypted data in a digital image.

In conventional (symmetric) encryption, the plaintext is usually represented as a binary stream and encrypted using an XOR type operation with a binary cipher. The algorithm used is ideally designed to: (i) generate a maximum entropy cipher so that there is no bias with regard to any …

Singular Superposition/Boundary Element Method For Reconstruction Of Multi-Dimensional Heat Flux Distributions With Application To Film Cooling Holes, Mahmood Silieti, Eduardo Divo, Alain J. Kassab

Publications

A hybrid singularity superposition/boundary element-based inverse problem method for the reconstruction of multi-dimensional heat flux distributions is developed. Cauchy conditions are imposed at exposed surfaces that are readily reached for measurements while convective boundary conditions are unknown at surfaces that are not amenable to measurements such as the walls of the cooling holes. The purpose of the inverse analysis is to determine the heat flux distribution along cooling hole surfaces. This is accomplished in an iterative process by distributing a set of singularities (sinks) inside the physical boundaries of the cooling hole (usually along cooling hole centerline) with a given …

Greedy Signal Recovery Review, Deanna Needell, Joel A. Tropp, Roman Vershynin

CMC Faculty Publications and Research

The two major approaches to sparse recovery are L1-minimization and greedy methods. Recently, Needell and Vershynin developed Regularized Orthogonal Matching Pursuit (ROMP) that has bridged the gap between these two approaches. ROMP is the first stable greedy algorithm providing uniform guarantees.

Even more recently, Needell and Tropp developed the stable greedy algorithm Compressive Sampling Matching Pursuit (CoSaMP). CoSaMP provides uniform guarantees and improves upon the stability bounds and RIC requirements of ROMP. CoSaMP offers rigorous bounds on computational cost and storage. In many cases, the running time is just O(NlogN), where N is the ambient dimension of the signal. This …

Electrothermal Imaging In One And Two Dimensions, Michael Janus, David Kibbling

Mathematical Sciences Technical Reports (MSTR)

Developing methods for the nondestructive testing of materials is an important area of research for industry. Situations often arise in which the integrity of an object is questioned, but testing it is very difficult. For example, a support bar may be embedded in a larger structure so that testing the bar’s integrity directly would require the impractical task of breaking down the larger structure. Instead, the ends of the bar might be accessible without dismantling the enclosing structure. The goal of nondestructive testing is to use methods that require taking measurements at the ends of the bar alone to give …

Determining The Shape Of A Resistor Grid, Esther Chiew, Vincent Selhorst-Jones

Mathematical Sciences Technical Reports (MSTR)

Impedance imaging has received a lot of attention in the past two decades, as a means for non-destructively imaging the interior of a conductive object. One injects a known electrical current pattern into an object at the exterior boundary, then measures the induced potential (voltage) on some portion of the boundary. The goal is to recover information about the interior conductivity of the object, which (we hope) influences the voltages we measure. Of course one can also use multiple input currents and measured voltages. A variation on this problem is that of "boundary identification," in which some portion of the …

Transient Analysis Of Heat And Mass Transfer By Natural Convection In Power-Law Fluid Past A Vertical Plate Immersed In A Porous Medium (Numerical Study), Nasser S. Elgazery

Applications and Applied Mathematics: An International Journal (AAM)

This paper attempted a numerical examination of the problem of unsteady free convection with heat and mass transfer from an isothermal vertical flat plate to a non-Newtonian fluid saturated porous medium. The flow in the porous medium was described via the Darcy-Brinkman-Forchheimer model. The simultaneous development of the problem of boundary layers was obtained numerically by finite difference method. Boundary layer and Boussinesq approximations had been incorporated. Numerical calculations were carried out for the various parameters entering into the problem. Velocity, temperature and concentration profiles were shown graphically and the physical quantities of the problem were given in tables. It …

Homotopy Perturbation Method And Padé Approximants For Solving Flierl-Petviashivili Equation, Syed T. Mohynd-Din, Muhammad A. Noor

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we present a reliable combination of homotopy perturbation method and Padé approximants to investigate the Flierl-Petviashivili (FP) equation. The approach introduces a new transformation necessary for the conversion of the Flierl-Petviashivili equation to a first order initial value problem and a reliable framework designed to overcome the difficulty of the singular point at x = 0. The proposed homotopy perturbation method is applied to the reformulated first order initial value problem which leads the solution in terms of transformed variable. The desired series solution is obtained by making use of the inverse transformation. The suggested algorithm may …

On A-Ary Subdivision For Curve Design Ii. 3-Point And 5-Point Interpolatory Schemes, Jian-Ao Lian

Applications and Applied Mathematics: An International Journal (AAM)

The a-ary 3-point and 5-point interpolatery subdivision schemes for curve design are introduced for arbitrary odd integer a greater than or equal to 3. These new schemes further extend the family of the classical 4- and 6-point interpolatory schemes.

A Reliable Approach For Higher-Order Integro-Differential Equations, Muhammad A. Noor, Syed T. Mohyud-Din

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we apply the variational iteration method (VIM) for solving higher-order integro differential equations by converting the problems into system of integral equations. The proposed technique is applied to the re-formulated system of integro-differential equations. Numerical results show the accuracy and efficiency of the suggested algorithm. The fact that the VIM solves nonlinear problems without calculating Adomian’s polynomials is a clear advantage of this technique over the decomposition method.

Chebyshev Collocation Method For The Effect Of Variable Thermal Conductivity On Micropolar Fluid Flow Over Vertical Cylinder With Variable Surface Temperature, Nasser S. Elgazery, Nader Y. Abd Elazem

Applications and Applied Mathematics: An International Journal (AAM)

An analysis is performed to study the role of a variable thermal conductivity on unsteady free convection in a micro-polar fluid past a semi-infinite vertical cylinder with variable surface temperature in the presence of magnetic filed and radiation. The surface temperature is measured to vary as a power of the axial coordinate measured from the leading edge of the cylinder. The governing non-linear partial differential equations are transformed into a linear algebraic system utilizing Chebyshev collocation method in spatial and Crank-Nicolson method in time. Numerical results for the velocity, angular velocity and temperature profiles as well as for the local …

Solving Higher Dimensional Initial Boundary Value Problems By Variational Iteration Decomposition Method, Muhammad A. Noor, Syed T. Mohyud-Din

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we apply a relatively new technique which is called the variational iteration decomposition method (VIDM) by combining the traditional variational iteration and the decomposition methods for solving higher dimensional initial boundary value problems. The proposed method is an elegant combination of variational iteration and the decomposition methods. The analytical results of the problems have been obtained in terms of convergent series with easily computable components. The method is quite efficient and is practically well suited for use in these problems. Several examples are given to verify the accuracy and efficiency of the proposed technique.

Model-Based Clustering Of Methylation Array Data: A Recursive-Partitioning Algorithm For High-Dimensional Data Arising As A Mixture Of Beta Distributions, E. Andres Houseman, Brock C. Christensen, Ru-Fang Yeh, Carmen J. Marsit, Margaret R. Karagas, Margaret Wrensch, Heather H. Nelson, Joseph Wiemels, Shichun Zheng, John K. Wiencke, Karl T. Kelsey

Harvard University Biostatistics Working Paper Series

No abstract provided.

Uniform Uncertainty Principle And Signal Recovery Via Regularized Orthogonal Matching Pursuit, Deanna Needell, Roman Vershynin

CMC Faculty Publications and Research

This paper seeks to bridge the two major algorithmic approaches to sparse signal recovery from an incomplete set of linear measurements—L1-minimization methods and iterative methods (Matching Pursuits). We find a simple regularized version of Orthogonal Matching Pursuit (ROMP) which has advantages of both approaches: the speed and transparency of OMP and the strong uniform guarantees of L1-minimization. Our algorithm, ROMP, reconstructs a sparse signal in a number of iterations linear in the sparsity, and the reconstruction is exact provided the linear measurements satisfy the uniform uncertainty principle.

On A-Ary Subdivision For Curve Design: I. 4-Point And 6-Point Interpolatory Schemes, Jian-Ao Lian

Applications and Applied Mathematics: An International Journal (AAM)

The classical binary 4-point and 6-point interpolatery subdivision schemes are generalized to a-ary setting for any integer a greater than or equal to 3. These new a-ary subdivision schemes for curve design are derived easily from their corresponding two-scale scaling functions, a notion from the context of wavelets.

Variational Iteration Method For Solving Initial And Boundary Value Problems Of Bratu-Type, Muhammad A. Noor, Syed T. Mohyud-Din

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we present a reliable framework to solve the initial and boundary value problems of Bratu-type which are widely applicable in fuel ignition of the combustion theory and heat transfer. The algorithm rests mainly on a relatively new technique, the variational iteration method. Several examples are given to confirm the efficiency and the accuracy of the proposed algorithm.

Signed Decomposition Of Fully Fuzzy Linear Systems, Tofigh Allahviranloo, Nasser Mikaeilvand, Narsis A. Kiani, Rasol M. Shabestari

Applications and Applied Mathematics: An International Journal (AAM)

System of linear equations is applied for solving many problems in various areas of applied sciences. Fuzzy methods constitute an important mathematical and computational tool for modeling real-world systems with uncertainties of parameters. In this paper, we discuss about fully fuzzy linear systems in the form AX = b (FFLS). A novel method for finding the non-zero fuzzy solutions of these systems is proposed. We suppose that all elements of coefficient matrix A are positive and we employ parametric form linear system. Finally, Numerical examples are presented to illustrate this approach and its results are compared with other methods.

Methods Of Assessing And Ranking Probable Sources Of Error, Nataniel Greene

Publications and Research

A classical method for ranking n potential events as sources of error is Bayes' theorem. However, a ranking based on Bayes' theorem lacks a fundamental symmetry: the ranking in terms of blame for error will not be the reverse of the ranking in terms of credit for lack of error. While this is not a flaw in Bayes' theorem, it does lead one to inquire whether there are related methods which have such symmetry. Related methods explored here include the logical version of Bayes' theorem based on probabilities of conditionals, probabilities of biconditionals, and ratios or differences of credit to …

Caching And Distributing Statistical Analyses In R, Roger D. Peng

Johns Hopkins University, Dept. of Biostatistics Working Papers

We present the cacher package for \proglang{R}, which provides tools for caching statistical analyses and for distributing these analyses to others in an efficient manner. The cacher package takes objects created by evaluating R expressions and stores them in key-value databases. These databases of cached objects can subsequently be assembled into packages for distribution over the web. The cacher package also provides tools to help readers examine the data and code in a statistical analysis and reproduce, modify, or improve upon the results. In addition, readers can easily conduct alternate analyses of the data. We describe the design and implementation …

Study Of Multibody Dynamics Using Matlab, Kalyani Soman

Mechanical & Aerospace Engineering Theses & Dissertations

Research in computational dynamics has gained importance in the past few years because of the demand for analysis and simulation of various multi-body systems in the growing mechanical and aerospace industry. These multi-body systems are made of individual bodies that are interconnected via mechanical joints. Mathematically, these joints that connect the bodies can be described as constraint equations imposed upon the motions of the involved free bodies. This process will result in an equation of motion expressed in the form of a differential-algebraic equation. It is the main objective of this thesis to develop a solution algorithm that can take …

A Wavelet-Based Method For Overcoming The Gibbs Phenomenon, 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. Here we describe a numerical procedure for overcoming the Gibbs phenomenon called the inverse wavelet reconstruction method. 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.

An Overview Of Conditionals And Biconditionals In Probability, Nataniel Greene

Publications and Research

Conditional and biconditional statements are a standard part of symbolic logic but they have only recently begun to be explored in probability for applications in artificial intelligence. Here we give a brief overview of the major theorems involved and illustrate them using two standard model problems from conditional probability.