High-Order Central WENO Schemes for 1D Hamilton-Jacobi Equations
|
|
- Theodore Poole
- 5 years ago
- Views:
Transcription
1 High-Order Central WENO Schemes for D Hamilton-Jacobi Equations Steve Bryson and Doron Levy Program in Scientific Computing/Computational Mathematics, Stanford University and the NASA Advanced Supercomputing Division, NASA Ames Research Center, Moffett Field, CA ; bryson@nas.nasa.gov Department of Mathematics, Stanford University, Stanford, CA ; dlevy@math.stanford.edu Introduction We consider high-order central approximations for solutions of one-dimensional Hamilton-Jacobi HJ) equations of the form t φx, t)+h φ x,x)=0, x R, ) subject to the initial data φx, t=0) = φ 0 x). Solutions for ) with smooth initial data typically remain continuous but develop discontinuous derivatives in finite time. Such solutions are not unique; the physically relevant solution is known as the viscosity solution see[,3,4,5,8,5]andthereferences therein). Various numerical methods were proposed in order to approximate the solutions of ). Examples for such methods are the high-order Godunov-type schemes that were introduced in [0, ], and were based on an Essentially Non-Oscillatory ENO) reconstruction step [7] that was evolved in time with a first-order monotone flux. The least dissipative monotone flux, the Godunov flux, requires solving Riemann problems at cell interfaces. A fifth-order Weighted ENO WENO) scheme, based on [0, 8], was introduced by Jiang and Peng [9]. Recently, Lin and Tadmor introduced in [6, 7] central schemes for approximating solutions of the HJ equation. These schemes are based on the Nessyahu-Tadmor scheme for approximating solutions of hyperbolic conservation laws [9]. Unlike upwind schemes, central schemes do not require Riemann solvers, which makes them attractive for solving systems of equations and for multi-dimensional problems. A second-order semi-discrete version of these schemes was introduced by Kurganov and Tadmor in []. While less dissipative, the semi-discrete scheme requires the estimation of the local speed of propagation, which is computationally intensive in particular in multi-dimensional problems. In a later work [], the numerical viscosity was further reduced by computing more precise information about local speed of propagation. To address the problem of schemes that are too computationally intensive, we introduced in [] efficient first- and second-order central schemes for approximating the solutions of multi-dimensional versions of ).
2 Steve Bryson and Doron Levy Unlike the previous attempts, our schemes in [] scale well with increasing dimension. In this paper we derive fully-discrete Central WENO CWENO) schemes for approximating solutions of ), which combine our previous works [, 3, 4]. We introduce third- and fifth-order accurate schemes, which are the first central schemes for the HJ equations of order higher than two. The core ingredient in the derivation of our schemes is a high-order CWENO reconstructions in space. Acknowledgment: We would like to thank Volker Elling for helpful discussions. CWENO Schemes for HJ Equations We are interested in approximating solutions of ) subject to the initial data φx, t =0) = φ 0 x). For simplicity we assume a uniform grid grid in space and time with mesh spacings, h := x and t. We denote the grid points by x i = i x, t n = n t, and the fixed mesh ratio by λ = t/ x. Let ϕ n i denote the approximate value of φ x i,t n ), and ϕ x ) n i denote the approximate value of the derivative φ x x i,t n ). We define ϕ n i := ϕn i+ ϕn i, ϕ n i := ϕn i ϕn i and 0 ϕ n i := ϕn i+ ϕn i. We assume that the approximate solution at time t n, ϕ n i is given. In order to approximate the solution at the next time step t n+, ϕ n+ i, we start by reconstructing a continuous piecewise-polynomial from the data, ϕ n i,and sample it at the half-integer points, {x i+/ }, in order to obtain the pointvalues of the interpolant at these points ϕ n i+/ as well as the derivative, ϕ i+/. We then evolve ϕn from time t n to time t n+ according to ), i+ t n+ )) ϕ x i+,tn+) = ϕ x i+,tn) H ϕ x x i+,t dt. ) t n This evolution is done at the half-integer grid points where the reconstruction is smooth as long as the CFL condition λ H ϕ x ) / is satisfied). Finally, in order to return to the original grid, we project ϕ n+ i+/ back onto the integer grid points {x i } toendupwithϕ n+ i. Since the evolution step ) is done at points where the solution is smooth, we can approximate the time integral at the RHS of ) using a sufficiently accurate quadrature rule. For example, for a third- and fourth-order method, this integral can be replaced by a Simpson s quadrature, t n+ t n )) H ϕ x x i+,t dt t [ H ϕ x x 6 i+,tn)) 3) )) +4H ϕ x x i+,tn+ +H ϕ x x i+,tn+))].
3 CWENO Schemes for HJ Equations 3 The intermediate values of the derivative in time, ϕ x xi+/,t n+/),and ϕ x xi+/,t n+), which are required in the quadrature 3), can be predicted using a Taylor expansion or with a Runge-Kutta RK) method. For details we refer the reader to [3, 9] and the references therein. The remaining ingredient is the piecewise-polynomial reconstruction in space. A careful study of the above procedure reveals that there are actually three different quantities that should be recovered in every time step. First, given ϕ i at time t n we need to reconstruct the point-values at the half-integer grid points, ϕ i+/, at the same time t n. This is the first term on the RHS of ). The second term on the RHS of ) requires evaluating the Hamiltonian H at the derivative ϕ i+/. Hence, the second quantity we should recover is ϕ i+/ from ϕ i. Finally, the predictor step that provides the values at the quadrature nodes in 3), require us to estimate ϕ i+/ from ϕ i+/ at every step of the RK method. In the next two sections we will focus on the reconstruction of these three quantities, first for a third-order method and then for a fifth-order method. The projection from ϕ n+ i+/ onto the original grid points to get ϕn+ i is accomplished using the same reconstruction used to approximate ϕ n i+/ from ϕ n i.. A Third-Order Scheme Following the above procedure, a third-order scheme can be generated by combining a third-order accurate ODE solver in time with a sufficiently high-order reconstruction in space. Here we present fourth-order CWENO reconstructions of the point values of ϕ i+/ and its derivative ϕ i+/. The reconstruction of ϕ i+/ from ϕ i. In order to obtain a fourth-order reconstruction of ϕ i+/ we will write a convex combination of two quadratic polynomials, ϕ [] constructed on a stencil which is left-biased with respect to x i+/, and the right-biased ϕ [] +, ϕ [] x) =ϕ i + ) ϕ i x xi )+ h ϕ [] + x) =ϕ i + h h ) ϕ i x xi )x x i+ )+O h 3), ) ϕ i x xi )+ h ) ϕ i x xi )x x i+ )+O h 3). An evaluation of these approximations at {x i+ } reads ) ϕ [] = ) 8 ϕ i +6ϕ i +3ϕ i+ ), ϕ [] + = 8 3ϕ i+6ϕ i+ ϕ i+ ). x i+ A straightforward computation shows that x i+ ϕ[] x i+ )+ ϕ[] + x i+ )=ϕ i+ + O h 4).
4 4 Steve Bryson and Doron Levy The fourth-order WENO estimate of ϕ i+/ is therefore given by the convex combination ) ) ) ϕ [4] w x i+ = w ϕ [] i+ x i+ + w + ϕ [] i+ + x i+, where the weights satisfy w i+/ + w+ i+/ =,w± i+/ 0, i. Insmooth regions we would like to satisfy wi w i + to attain an O ) h 4) error, while when the stencil {x i,x i,x i+,x i+ } supporting ϕ w x i+ contains a discontinuity, the weight of the more oscillatory polynomial should vanish. Following [0, 8], we meet these requirements by setting w k i+ = αk i+ l αl i+, α k i+ = c k ) p 4) ɛ + S k i+ where k, l {+, } k and l will range over a larger space of symbols when we use more interpolants). The constants c ± =/ and are independent of the grid-point. We choose ɛ as 0 6 to prevents the denominator in 4) from vanishing, and set p = see [0]). The smoothness measures S i ± should be large when ϕ is nearly singular. Following the standard practice with WENOtype schemes [0], we take S ± i to be the sum of the L -norms of the first and second derivatives on the stencil supporting ϕ [] ±. If we approximate the first derivative at x i+/ by h + ϕ i+/, the second derivative by h ϕ i+/, and define the smoothness measure S i+ [r, s] =h s j=r h + ϕ i+j+ ) + h s j=r+ then for the fourth-order interpolation of ϕ w x i+ S i+/ [, 0] and S + i+/ = S i+/ [0, ]. ) h + ϕ i+j+, 5) ) we have S i+/ = The reconstruction of ϕ i+/ from ϕ i. To obtain a fourth-order estimate of the derivative ϕ x i+/ ) from ϕx i ), we start from the cubic interpolants x) =ϕ i + ) ϕ i x xi )+ h h ) ϕ i x xi )x x i+ ) + 6h 3 ) ϕ i x xi )x x i+ )x x i )+O h 4), + x) =ϕ i + ) ϕ i x xi )+ h h ) ϕ i x xi )x x i+ ) + 6h 3 ) ϕ i x xi )x x i+ )x x i+ )+O h 4). Differentiating ± at x i+
5 ϕ [3],i+ ϕ [3] +,i+ CWENO Schemes for HJ Equations 5 = 4h ϕ i 3ϕ i ϕ i +3ϕ i+ ), = 4h 3ϕ i +ϕ i+ +3ϕ i+ ϕ i+3 ). Again, ϕ [3] +,i+ ϕ [3] = ϕ +,i+ i+ + O h 4), ) and a fourth-order WENO reconstruction of ϕ x i+ is ϕ [4] i+/ = w ϕ [3] + w + ϕ [3] i+,i+ i+ +,i+ where the weights are of the form 4) with c ± = / and S i+/ = S i+/ [, 0] and S + i+/ = S i+/ [0, ]. The reconstruction of ϕ i+/ from ϕ i+/. Repeating the above procedure, this time with three quadratic interpolants ϕ [] x) =ϕ i+ + ) ) ϕ h i+ x x i+ + ) ) ) h ϕ i+ x x i+ x x i+ 3 + O h 3), ϕ [] 0 x) =ϕ i+ + ) ) 0 ϕ h i+ x x i+ + ) ) ) h ϕ i+ x x i x x i+ 3 + O h 3), ϕ [] + x) =ϕ i+ + ) ) ϕ h i+ x x i+ + ) ) ) h ϕ i+ x x i+ x x i+ 3 + O h 3), results with 6 ϕ [] +,i+ 3 ϕ [] + 0,i+ 6 ϕ [] = ϕ +,i+ i+ + O h 4), where ϕ [],i+ ϕ [] +,i+ = h ϕ i 3 4ϕ i +3ϕ i+ ), ϕ [] 0,i+ = h 3ϕ i+ +4ϕ i+ 3 ϕ i+ 5 ). The fourth-order WENO estimate of ϕ i+/ is = h ϕ i+ 3 ϕ i ), ϕ [4] i+/ = w ϕ [] + w 0 i+,i+ i+ ϕ [] + w + ϕ [] 0,i+ i+ +,i+ where the weights w are of the form 4) with c = c + =/6,c 0 =/3, and the oscillatory indicators S i+/ = S i+/ [, ], S i+/ = S i+/ [, 0], and S + i+/ = S i+/ [0, ].
6 6 Steve Bryson and Doron Levy. A Fifth-Order Scheme Once again, similarly to the third-order scheme, we need to reconstruct the point-values of ϕ and ϕ. We start with the reconstruction of ϕ i+/ and ϕ i+/ from ϕ i. We write sixth-order interpolants as a convex combination of cubic interpolants, x) and + x) introduced above and 0 x) =ϕ i + ) ϕ i x xi )+ h h ) ϕ i x xi )x x i+ ) + 6h 3 ) ϕ i x xi )x x i+ )x x i+ )+O h 4). In this case where In a similar way, where 3 6 ϕ[3] + 5,i+ 8 ϕ[3] + 3 0,i+ 6 ϕ[3] +,i+,i+ 0,i+ +,i+ = ϕ i+ + O h 6), = 6 ϕ i 5ϕ i +5ϕ i +5ϕ i+ ), = 6 ϕ i +9ϕ i +9ϕ i+ ϕ i+ ), = 6 5ϕ i +5ϕ i+ 5ϕ i+ + ϕ i+3 ) ϕ [3] + 49,i+ 40 ϕ [3] 9 0,i+ 80 ϕ [3] = ϕ +,i+ i+/ + O h 6), ϕ [3],i+ ϕ [3] 0,i+ ϕ [3] +,i+ = 4h ϕ i 3ϕ i ϕ i +3ϕ i+ ), = 4h ϕ i 7ϕ i +7ϕ i+ ϕ i+ ), = 4h 3ϕ i +ϕ i+ +3ϕ i+ ϕ i+3 ). The sixth-order WENO estimates for ϕ i+/ and ϕ i+/ are ϕ [6] i+ ϕ [6] i+ = w + w 0 i+,i+ i+ + w +, 0,i+ i+ +,i+ = w ϕ [3] + w 0 i+,i+ i+ ϕ [3] + w + ϕ [3], 0,i+ i+ +,i+ where the weights for ϕ are given by 4), with c = c + =3/6,c 0 =5/8 and the oscillatory indicators are S i+/ = S i+/ [, 0], Si+/ 0 = S i+/ [, ] and S + i+/ = S i+/ [0, ]. The negative weights for ϕ require special treatment see [] for details). Following [] we split the positive and negative
7 CWENO Schemes for HJ Equations 7 weights in the following way: first, we set γ = γ+ =9/40, γ0 =49/40 and γ+ = γ+ + =9/80, γ0 + =49/0. Then, For k, l {, 0, +}, setσ ± = k γk ± so that similarly to 4), α k ±,i+ = γ k ± σ ± ɛ + S k i+ ) p and w k i+ α k +,i+ = σ + l αl +,i+ α k,i+ σ. l αl,i+ Because i+/ and ϕ [3] i+/ are defined on the same stencils, they use the same smoothness measures S i+/. All that is left is the reconstruction of ϕ i+/ from ϕ i+/. In this case a sixth-order approximation to ϕ i+/ requires a weighted sum of four cubic interpolants. This reconstruction is similar to the previous ones. We skip the details and summarize the result: where ϕ [6] = w ϕ [3] + w 0 ϕ [3] + w 0+ ϕ [3] + w + ϕ [3], i+ i+,i+ i+ 0,i+ i+ 0+,i+ i+ +,i+ ϕ [3],i+ ϕ [3] 0,i+ ϕ [3] 0+,i+ ϕ [3] +,i+ = 6h ϕ i 5 +9ϕ i 3 8ϕ i +ϕ i+ ), = 6h ϕ i 3 6ϕ i +3ϕ i+ +ϕ i+ 3 ), = 6h ϕ i 3ϕ i+ +6ϕ i+ 3 ϕ i+ 5 ), = 6h ϕ i+ +8ϕ i+ 3 9ϕ i+ 5 +ϕ i+ 7 ). Here, c = c + =/0,c 0 = c 0 +=9/0, S i+/ = S i+/ [ 3, ], S 0 i+/ = S i+/ [, 0], S 0+ i+/ = S i+/ [, ] and S + i+/ = S i+/ [0, ]. 3 Numerical Examples In all our numerical simulations, the ODE solvers we use are the non-linear fourth-order Strong-Stability Preserving Runge-Kutta SSP-RK) methods of [6]. We start by testing the accuracy of our new CWENO methods when approximating the solution of the linear advection equation, ϕ t + ϕ x =0. The initial data is taken as ϕ x, 0) = sin 4 πx), the mesh ratio λ =0.9 and the time T = 4. The results obtained with the fifth-order method of. are shown in Table.
8 8 Steve Bryson and Doron Levy Table. Error and convergence rate for linear advection with initial condition ϕ x, 0) = sin 4 πx) N L error L order Next, we test the CWENO methods with two nonlinear Hamiltonians: a convex Hamiltonian ϕ t + ϕ x +) = 0 and a non-convex Hamiltonian ϕ t cos ϕ x + ) = 0. The interval is [0, ], the boundary conditions are periodic and the initial conditions for both Hamiltonians are taken as ϕ x, 0) = cos πx). The exact solution to both problems is smooth until t /π, after which a singularity forms. A second singularity forms in the non-convex H example at t.9/π. The results of the accuracy test with the fifth-order method are shown in Table, and the solution at time T =.5/π is plotted in Figure. Following [9] the errors in Table after the formation of the singularity are computed at a distance of 0. away from any singularities. Table. L Error and convergence rate estimates for convex and non-convex Hamiltonians. top: T =0.5/π, bottom: T =.5/π. λ =0.3 N convex convex non-convex non-convex L error L order L error L order N convex convex non-convex non-convex L error L order L error L order
9 CWENO Schemes for HJ Equations approximation exact.5 approximation exact x x com- Fig.. left: Convex Hamiltonian right: non-convex Hamiltonian at T =.5 π pared with the exact solution, N = 00. References. Barles G., Solution de viscosité des équations de Hamilton-Jacobi, Springer- Verlag, Berlin, Bryson S., Levy D., Central Schemes for Multi-Dimensional Hamilton-Jacobi Equations, NASA Technical Report NAS 0 0, 00, submitted. 3. Crandall M.G., Evans L.C., Lions P.-L., Some properties of viscosity solutions of Hamilton-Jacobi equations, Trans. Amer. Math. Soc., 8, 984), pp Crandall M.G., Ishii H., Lions P.-L., User s guide to viscosity solutions of second order partial differential equations, Bull. Amer. Math. Soc., 7, 99), pp Crandall M.G., Lions P.-L., Viscosity solutions of Hamilton-Jacobi equations, Trans. Amer. Math. Soc., ), pp Gottlieb S., Shu C.-W., Tadmor E., Strong stability-preserving high order time discretization methods, SIAMReview, 43, 00), pp Harten A., Engquist B., Osher S., Chakravarthy S., Uniformly High Order Accurate Essentially Non-oscillatory Schemes III, JCP, 7, 987), pp Kruzkov S.N., The Cauchy problem in the large for nonlinear equations and for certain quasilinear systems of the first order with several variables, Soviet Math. Dokl., 5, 964), pp Jiang G.-S., Peng D., Weighted ENO schemes for Hamilton-Jacobi equations, SIAM J. Sci. Comp.,, 000), pp Jiang G.-S., Shu C.-W., Efficient Implementation of Weighted ENO Schemes, JCP, 6, 996), pp Kurganov A., Noelle S., Petrova G., Semi-discrete central-upwind schemes for hyperbolic conservation laws and Hamilton-Jacobi equations, SIAMJ.Sci. Comp., to appear.. Kurganov A., Tadmor E., New high-resolution semi-discrete central schemes for Hamilton-Jacobi equations, JCP, 60, 000), pp Levy D., Puppo G., Russo G., Central WENO schemes for hyperbolic systems of conservation laws, Math. Model. and Numer. Anal., 33, no ), pp Levy D., Puppo G., Russo G., Compact central WENO schemes for multidimensional conservation laws, SIAMJ. Sci. Comp.,, 000), pp Lions P.L., Generalized solutions of Hamilton-Jacobi equations, Pitman, London, 98.
10 0 Steve Bryson and Doron Levy 6. Lin C.-T., Tadmor E., L -stability and error estimates for approximate Hamilton-Jacobi solutions, Numer. Math., 87, 00), pp Lin C.-T., Tadmor E., High-resolution non-oscillatory central schemes for approximate Hamilton-Jacobi equations, SIAMJ. Sci. Comp.,, no. 6, 000), pp Liu X.-D., Osher S., Chan T., Weighted Essentially Non-oscillatory Schemes, JCP, 5, 994), pp Nessyahu H., Tadmor E., Non-oscillatory central differencing for hyperbolic conservation laws, JCP, 87, no. 990), pp Osher S., Sethian J., Fronts propagating with curvature dependent speed: algorithms based on Hamilton-Jacobi formulations, JCP, 79, 988), pp Osher S., Shu C.-W., High-order essentially nonoscillatory schemes for Hamilton-Jacobi equations, SIAMJ. Numer. Anal., 8, 99), pp Shi J., Hu C., Shu C.-W., A technique of treating negative weights in WENO schemes, JCP, to appear.
Fast sweeping methods and applications to traveltime tomography
Fast sweeping methods and applications to traveltime tomography Jianliang Qian Wichita State University and TRIP, Rice University TRIP Annual Meeting January 26, 2007 1 Outline Eikonal equations. Fast
More informationNumerical Methods for Optimal Control Problems. Part II: Local Single-Pass Methods for Stationary HJ Equations
Numerical Methods for Optimal Control Problems. Part II: Local Single-Pass Methods for Stationary HJ Equations Ph.D. course in OPTIMAL CONTROL Emiliano Cristiani (IAC CNR) e.cristiani@iac.cnr.it (thanks
More informationA Toolbox of Hamilton-Jacobi Solvers for Analysis of Nondeterministic Continuous and Hybrid Systems
A Toolbox of Hamilton-Jacobi Solvers for Analysis of Nondeterministic Continuous and Hybrid Systems Ian Mitchell Department of Computer Science University of British Columbia Jeremy Templeton Department
More informationA Local Ordered Upwind Method for Hamilton-Jacobi and Isaacs Equations
A Local Ordered Upwind Method for Hamilton-Jacobi and Isaacs Equations S. Cacace E. Cristiani M. Falcone Dipartimento di Matematica, SAPIENZA - Università di Roma, Rome, Italy (e-mail: cacace@mat.uniroma1.it).
More informationPostprocessing of nonuniform MRI
Postprocessing of nonuniform MRI Wolfgang Stefan, Anne Gelb and Rosemary Renaut Arizona State University Oct 11, 2007 Stefan, Gelb, Renaut (ASU) Postprocessing October 2007 1 / 24 Outline 1 Introduction
More informationA second-order fast marching eikonal solver a
A second-order fast marching eikonal solver a a Published in SEP Report, 100, 287-292 (1999) James Rickett and Sergey Fomel 1 INTRODUCTION The fast marching method (Sethian, 1996) is widely used for solving
More informationAn adaptive finite-difference method for traveltimes and amplitudes
GEOPHYSICS, VOL. 67, NO. (JANUARY-FEBRUARY 2002); P. 67 76, 6 FIGS., 2 TABLES. 0.90/.45472 An adaptive finite-difference method for traveltimes and amplitudes Jianliang Qian and William W. Symes ABSTRACT
More informationAn efficient discontinuous Galerkin method on triangular meshes for a. pedestrian flow model. Abstract
An efficient discontinuous Galerkin method on triangular meshes for a pedestrian flow model Yinhua Xia 1,S.C.Wong 2, Mengping Zhang 3,Chi-WangShu 4 and William H.K. Lam 5 Abstract In this paper, we develop
More informationParaxial Eikonal Solvers for Anisotropic Quasi-P Travel Times
Journal of Computational Physics 73, 256 278 (200) doi:0.006/jcph.200.6875, available online at http://www.idealibrary.com on Paraxial Eikonal Solvers for Anisotropic Quasi-P Travel Times Jianliang Qian
More informationRICE UNIVERSITY. 3 D First Arrival Traveltimes and Amplitudes via. Eikonal and Transport Finite Dierence Solvers. Maissa A.
RICE UNIVERSITY 3 D First Arrival Traveltimes and Amplitudes via Eikonal and Transport Finite Dierence Solvers by Maissa A. Abd El-Mageed A Thesis Submitted in Partial Fulfillment of the Requirements for
More informationTHE EIKONAL EQUATION: NUMERICAL EFFICIENCY VS. ALGORITHMIC COMPLEXITY ON QUADRILATERAL GRIDS. 1. Introduction. The Eikonal equation, defined by (1)
Proceedings of ALGORITMY 2005 pp. 22 31 THE EIKONAL EQUATION: NUMERICAL EFFICIENCY VS. ALGORITHMIC COMPLEXITY ON QUADRILATERAL GRIDS SHU-REN HYSING AND STEFAN TUREK Abstract. This paper presents a study
More information, SIAM GS 13 Conference, Padova, Italy
2013-06-18, SIAM GS 13 Conference, Padova, Italy A Mixed Order Scheme for the Shallow Water Equations on the GPU André R. Brodtkorb, Ph.D., Research Scientist, SINTEF ICT, Department of Applied Mathematics,
More informationAn Adjoint State Method for Numerical Approximation of Continuous Traffic Congestion Equilibria
An Adjoint State Method for Numerical Approximation of Continuous Traffic Congestion Equilibria Songting Luo Shingyu Leung Jianliang Qian Abstract The equilibrium metric for minimizing a continuous congested
More informationSeismology and Seismic Imaging
Seismology and Seismic Imaging 5. Ray tracing in practice N. Rawlinson Research School of Earth Sciences, ANU Seismology lecture course p.1/24 Introduction Although 1-D whole Earth models are an acceptable
More informationFast-marching eikonal solver in the tetragonal coordinates
Stanford Exploration Project, Report 97, July 8, 1998, pages 241 251 Fast-marching eikonal solver in the tetragonal coordinates Yalei Sun and Sergey Fomel 1 keywords: fast-marching, Fermat s principle,
More informationFast-marching eikonal solver in the tetragonal coordinates
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 499?? Fast-marching eikonal solver in the tetragonal coordinates Yalei Sun and Sergey Fomel 1 ABSTRACT Accurate and efficient traveltime
More informationBeamforming in Interference Networks for Uniform Linear Arrays
Beamforming in Interference Networks for Uniform Linear Arrays Rami Mochaourab and Eduard Jorswieck Communications Theory, Communications Laboratory Dresden University of Technology, Dresden, Germany e-mail:
More informationA Fast Marching Method for Hamilton-Jacobi Equations Modeling Monotone Front Propagations
A Fast Marching Method for Hamilton-Jacobi Equations Modeling Monotone Front Propagations Emiliano Cristiani November 15, 2008 Abstract In this paper we present a generalization of the Fast Marching method
More informationPOWER FLOW SOLUTION METHODS FOR ILL- CONDITIONED SYSTEMS
104 POWER FLOW SOLUTION METHODS FOR ILL- CONDITIONED SYSTEMS 5.1 INTRODUCTION: In the previous chapter power flow solution for well conditioned power systems using Newton-Raphson method is presented. The
More informationA Primer on Image Segmentation. Jonas Actor
A Primer on Image Segmentation It s all PDE s anyways Jonas Actor Rice University 21 February 2018 Jonas Actor Segmentation 21 February 2018 1 Table of Contents 1 Motivation 2 Simple Methods 3 Edge Methods
More informationWavelet Transform. From C. Valens article, A Really Friendly Guide to Wavelets, 1999
Wavelet Transform From C. Valens article, A Really Friendly Guide to Wavelets, 1999 Fourier theory: a signal can be expressed as the sum of a series of sines and cosines. The big disadvantage of a Fourier
More informationMath 148 Exam III Practice Problems
Math 48 Exam III Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture, all homework problems, all lab
More informationAn Exact Algorithm for Calculating Blocking Probabilities in Multicast Networks
An Exact Algorithm for Calculating Blocking Probabilities in Multicast Networks Eeva Nyberg, Jorma Virtamo, and Samuli Aalto Laboratory of Telecommunications Technology Helsinki University of Technology
More informationResource Allocation Challenges in Future Wireless Networks
Resource Allocation Challenges in Future Wireless Networks Mohamad Assaad Dept of Telecommunications, Supelec - France Mar. 2014 Outline 1 General Introduction 2 Fully Decentralized Allocation 3 Future
More informationAn Adjoint State Method for Numerical Approximation of Continuous Traffic Congestion Equilibria
3 4 5 6 7 8 9 3 4 5 6 7 8 9 Commun. Comput. Phys. doi: 8/cicp..3a Vol. x, No. x, pp. -9 xxx An Adjoint State Method for Numerical Approximation of Continuous Traffic Congestion Equilibria Songting Luo,,
More informationDynamic Programming in Real Life: A Two-Person Dice Game
Mathematical Methods in Operations Research 2005 Special issue in honor of Arie Hordijk Dynamic Programming in Real Life: A Two-Person Dice Game Henk Tijms 1, Jan van der Wal 2 1 Department of Econometrics,
More informationAn Adjoint State Method For Three-dimensional Transmission Traveltime Tomography Using First-Arrivals
An Adjoint State Method For Three-dimensional Transmission Traveltime Tomography Using First-Arrivals Shingyu Leung Jianliang Qian January 3, 6 Abstract Traditional transmission travel-time tomography
More informationLIMIT CYCLES FROM A CUBIC REVERSIBLE SYSTEM VIA THE THIRD-ORDER AVERAGING METHOD
Electronic Journal of Differential Equations, Vol. 5 5, No., pp. 7. ISSN: 7-669. UL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu LIMIT CYCLES FOM A CUBIC EVESIBLE
More informationSoliton-effect compression and dispersive radiation
1 March 000 Optics Communications 175 000 469 475 www.elsevier.comrlocateroptcom Soliton-effect compression and dispersive radiation Noel F. Smyth ) Department of Mathematics and Statistics, UniÕersity
More informationDigital Processing of Continuous-Time Signals
Chapter 4 Digital Processing of Continuous-Time Signals 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 4-1-1 Digital Processing of Continuous-Time Signals Digital
More informationDigital Processing of
Chapter 4 Digital Processing of Continuous-Time Signals 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 4-1-1 Digital Processing of Continuous-Time Signals Digital
More informationWavelet Transform. From C. Valens article, A Really Friendly Guide to Wavelets, 1999
Wavelet Transform From C. Valens article, A Really Friendly Guide to Wavelets, 1999 Fourier theory: a signal can be expressed as the sum of a, possibly infinite, series of sines and cosines. This sum is
More informationTHE idea that a signal which has energy in only a limited
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 45, NO. 5, JULY 1999 1555 Minimum Rate Sampling and Reconstruction of Signals with Arbitrary Frequency Support Cormac Herley, Member, IEEE, and Ping Wah Wong,
More information2.1 Partial Derivatives
.1 Partial Derivatives.1.1 Functions of several variables Up until now, we have only met functions of single variables. From now on we will meet functions such as z = f(x, y) and w = f(x, y, z), which
More informationA STUDY OF EULERIAN NUMBERS FOR PERMUTATIONS IN THE ALTERNATING GROUP
INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 6 (2006), #A31 A STUDY OF EULERIAN NUMBERS FOR PERMUTATIONS IN THE ALTERNATING GROUP Shinji Tanimoto Department of Mathematics, Kochi Joshi University
More informationElemental Image Generation Method with the Correction of Mismatch Error by Sub-pixel Sampling between Lens and Pixel in Integral Imaging
Journal of the Optical Society of Korea Vol. 16, No. 1, March 2012, pp. 29-35 DOI: http://dx.doi.org/10.3807/josk.2012.16.1.029 Elemental Image Generation Method with the Correction of Mismatch Error by
More informationSpectral Feature of Sampling Errors for Directional Samples on Gridded Wave Field
Spectral Feature of Sampling Errors for Directional Samples on Gridded Wave Field Ming Luo, Igor G. Zurbenko Department of Epidemiology and Biostatistics State University of New York at Albany Rensselaer,
More informationNonuniform multi level crossing for signal reconstruction
6 Nonuniform multi level crossing for signal reconstruction 6.1 Introduction In recent years, there has been considerable interest in level crossing algorithms for sampling continuous time signals. Driven
More informationReal- Time Computer Vision and Robotics Using Analog VLSI Circuits
750 Koch, Bair, Harris, Horiuchi, Hsu and Luo Real- Time Computer Vision and Robotics Using Analog VLSI Circuits Christof Koch Wyeth Bair John. Harris Timothy Horiuchi Andrew Hsu Jin Luo Computation and
More informationRouting in Massively Dense Static Sensor Networks
Routing in Massively Dense Static Sensor Networks Eitan ALTMAN, Pierre BERNHARD, Alonso SILVA* July 15, 2008 Altman, Bernhard, Silva* Routing in Massively Dense Static Sensor Networks 1/27 Table of Contents
More informationEFFECTS OF PHASE AND AMPLITUDE ERRORS ON QAM SYSTEMS WITH ERROR- CONTROL CODING AND SOFT DECISION DECODING
Clemson University TigerPrints All Theses Theses 8-2009 EFFECTS OF PHASE AND AMPLITUDE ERRORS ON QAM SYSTEMS WITH ERROR- CONTROL CODING AND SOFT DECISION DECODING Jason Ellis Clemson University, jellis@clemson.edu
More informationTime-average constraints in stochastic Model Predictive Control
Time-average constraints in stochastic Model Predictive Control James Fleming Mark Cannon ACC, May 2017 James Fleming, Mark Cannon Time-average constraints in stochastic MPC ACC, May 2017 1 / 24 Outline
More informationEnumeration of Two Particular Sets of Minimal Permutations
3 47 6 3 Journal of Integer Sequences, Vol. 8 (05), Article 5.0. Enumeration of Two Particular Sets of Minimal Permutations Stefano Bilotta, Elisabetta Grazzini, and Elisa Pergola Dipartimento di Matematica
More informationApplications of Monte Carlo Methods in Charged Particles Optics
Sydney 13-17 February 2012 p. 1/3 Applications of Monte Carlo Methods in Charged Particles Optics Alla Shymanska alla.shymanska@aut.ac.nz School of Computing and Mathematical Sciences Auckland University
More informationConvolution Pyramids. Zeev Farbman, Raanan Fattal and Dani Lischinski SIGGRAPH Asia Conference (2011) Julian Steil. Prof. Dr.
Zeev Farbman, Raanan Fattal and Dani Lischinski SIGGRAPH Asia Conference (2011) presented by: Julian Steil supervisor: Prof. Dr. Joachim Weickert Fig. 1.1: Gradient integration example Seminar - Milestones
More informationChapter-2 SAMPLING PROCESS
Chapter-2 SAMPLING PROCESS SAMPLING: A message signal may originate from a digital or analog source. If the message signal is analog in nature, then it has to be converted into digital form before it can
More informationTRANSFORMS / WAVELETS
RANSFORMS / WAVELES ransform Analysis Signal processing using a transform analysis for calculations is a technique used to simplify or accelerate problem solution. For example, instead of dividing two
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:0.038/nature727 Table of Contents S. Power and Phase Management in the Nanophotonic Phased Array 3 S.2 Nanoantenna Design 6 S.3 Synthesis of Large-Scale Nanophotonic Phased
More informationEikonal equations on the Sierpinski gasket 1. Fabio Camilli SBAI-"Sapienza" Università di Roma
Eikonal equations on the Sierpinski gasket 1 Fabio Camilli SBAI-"Sapienza" Università di Roma 1 F. CAMILLI, R.CAPITANELLI, C. MARCHI, Eikonal equations on the Sierpinski gasket, arxiv:1404.3692, 2014 Fabio
More informationENO morphological wavelet and its application in signal processing
Available online www.jocpr.com Journal of Chemical and Pharmaceutical Research, 2014, 6(6):1339-1346 Research Article ISSN : 0975-7384 CODEN(USA) : JCPRC5 ENO morphological wavelet and its application
More informationROBUST SUPERDIRECTIVE BEAMFORMER WITH OPTIMAL REGULARIZATION
ROBUST SUPERDIRECTIVE BEAMFORMER WITH OPTIMAL REGULARIZATION Aviva Atkins, Yuval Ben-Hur, Israel Cohen Department of Electrical Engineering Technion - Israel Institute of Technology Technion City, Haifa
More informationReview guide for midterm 2 in Math 233 March 30, 2009
Review guide for midterm 2 in Math 2 March, 29 Midterm 2 covers material that begins approximately with the definition of partial derivatives in Chapter 4. and ends approximately with methods for calculating
More informationOn Waveform Design for MIMO Radar with Matrix Completion
On Waveform Design for MIMO Radar with Matrix Completion Shunqiao Sun and Athina P. Petropulu ECE Department, Rutgers, The State University of New Jersey, Piscataway, NJ, 08854 Email: {shunq.sun, athinap}@rutgers.edu
More informationA FAST EIKONAL EQUATION SOLVER USING THE SCHRÖDINGER WAVE EQUATION
A FAST EIKONAL EQUATION SOLVER USING THE SCHRÖDINGER WAVE EQUATION KARTHIK S. GURUMOORTHY AND ANAND RANGARAJAN Abstract. We use a Schrödinger wave equation formalism to solve the eikonal equation. We show
More informationOn Optimum Communication Cost for Joint Compression and Dispersive Information Routing
2010 IEEE Information Theory Workshop - ITW 2010 Dublin On Optimum Communication Cost for Joint Compression and Dispersive Information Routing Kumar Viswanatha, Emrah Akyol and Kenneth Rose Department
More informationThis exam contains 9 problems. CHECK THAT YOU HAVE A COMPLETE EXAM.
Math 126 Final Examination Winter 2012 Your Name Your Signature Student ID # Quiz Section Professor s Name TA s Name This exam contains 9 problems. CHECK THAT YOU HAVE A COMPLETE EXAM. This exam is closed
More informationHigher-order schemes for 3D first-arrival traveltimes and amplitudes
GEOPHYSICS, VOL. 77, NO. (MARCH-APRIL ); P. T47 T56, FIGS..9/GEO-363. Higher-order schemes for 3D first-arrival traveltimes and amplitudes Songting Luo, Jianliang Qian, and Hongkai Zhao ABSTRACT In the
More informationPerformance Analysis of a 1-bit Feedback Beamforming Algorithm
Performance Analysis of a 1-bit Feedback Beamforming Algorithm Sherman Ng Mark Johnson Electrical Engineering and Computer Sciences University of California at Berkeley Technical Report No. UCB/EECS-2009-161
More informationMINIMIZING SELECTIVE AVAILABILITY ERROR ON TOPEX GPS MEASUREMENTS. S. C. Wu*, W. I. Bertiger and J. T. Wu
MINIMIZING SELECTIVE AVAILABILITY ERROR ON TOPEX GPS MEASUREMENTS S. C. Wu*, W. I. Bertiger and J. T. Wu Jet Propulsion Laboratory California Institute of Technology Pasadena, California 9119 Abstract*
More informationMaximum Likelihood Detection of Low Rate Repeat Codes in Frequency Hopped Systems
MP130218 MITRE Product Sponsor: AF MOIE Dept. No.: E53A Contract No.:FA8721-13-C-0001 Project No.: 03137700-BA The views, opinions and/or findings contained in this report are those of The MITRE Corporation
More informationA moment-preserving approach for depth from defocus
A moment-preserving approach for depth from defocus D. M. Tsai and C. T. Lin Machine Vision Lab. Department of Industrial Engineering and Management Yuan-Ze University, Chung-Li, Taiwan, R.O.C. E-mail:
More informationFast Computation for Secure Communication with Full-Duplex Radio
Fast Computation for Secure Communication with Full-Duplex Radio Lei Chen, Qiping Zhu, Yingbo Hua Dept of Electrical and Computer Eng, University of California, Riverside, CA 9252, USA Email: lchen@ucredu,
More informationAdaptive Rate Transmission for Spectrum Sharing System with Quantized Channel State Information
Adaptive Rate Transmission for Spectrum Sharing System with Quantized Channel State Information Mohamed Abdallah, Ahmed Salem, Mohamed-Slim Alouini, Khalid A. Qaraqe Electrical and Computer Engineering,
More information1813. Two-way collinear interaction of longitudinal waves in an elastic medium with quadratic nonlinearity
83. Two-way collinear interaction of longitudinal waves in an elastic medium with quadratic nonlinearity Zhenghao Sun, Fucai Li 2, Hongguang Li 3 State Key Laboratory of Mechanical System and Vibration,
More informationT he Parrondo s paradox describes the counterintuitive situation where combining two individually-losing
OPEN SUBJECT AREAS: APPLIED MATHEMATICS COMPUTATIONAL SCIENCE Received 6 August 013 Accepted 11 February 014 Published 8 February 014 Correspondence and requests for materials should be addressed to J.-J.S.
More informationGauss and AGM. Burton Rosenberg. January 30, 2004
Gauss and AGM Burton Rosenberg January 3, 24 Introduction derivation of equation. what has it to do w/ the lemniscate agm properties of I elliptic integrals The Elliptic Integral of the First Kind Define
More informationThe fast marching method in Spherical coordinates: SEG/EAGE salt-dome model
Stanford Exploration Project, Report 97, July 8, 1998, pages 251 264 The fast marching method in Spherical coordinates: SEG/EAGE salt-dome model Tariq Alkhalifah 1 keywords: traveltimes, finite difference
More informationWavelet-based image compression
Institut Mines-Telecom Wavelet-based image compression Marco Cagnazzo Multimedia Compression Outline Introduction Discrete wavelet transform and multiresolution analysis Filter banks and DWT Multiresolution
More informationINFLUENCE OF VORTEX STRUCTURES ON PRESSURE AND ULTRASOUND IN VORTEX FLOW-METERS
INFLUENCE OF VORTEX STRUCTURES ON PRESSURE AND ULTRASOUND IN VORTEX FLOW-METERS V. Hans*, H. Windorfer*, S. Perpeet** *Institute of Measurement and Control **Institute of Turbomachinery University of Essen,
More informationCounting Snakes, Differentiating the Tangent Function, and Investigating the Bernoulli-Euler Triangle by Harold Reiter
Counting Snakes, Differentiating the Tangent Function, and Investigating the Bernoulli-Euler Triangle by Harold Reiter In this paper we will examine three apparently unrelated mathematical objects One
More informationTopic 7f Time Domain FDM
Course Instructor Dr. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: rcrumpf@utep.edu Topic 7f Time Domain FDM EE 4386/5301 Computational Methods in EE Topic 7f Time Domain FDM 1 Outline
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 12, December-2013 ISSN
International Journal of Scientific & Engineering Research, Volume, Issue, December- ISSN 9-558 9 Application of Error s by Generalized Neuron Model under Electric Short Term Forecasting Chandragiri Radha
More informationPHYSICS-BASED THRESHOLD VOLTAGE MODELING WITH REVERSE SHORT CHANNEL EFFECT
Journal of Modeling and Simulation of Microsystems, Vol. 2, No. 1, Pages 51-56, 1999. PHYSICS-BASED THRESHOLD VOLTAGE MODELING WITH REVERSE SHORT CHANNEL EFFECT K-Y Lim, X. Zhou, and Y. Wang School of
More informationON PARALLEL ALGORITHMS FOR SOLVING THE DIRECT AND INVERSE PROBLEMS OF IONOSPHERIC SOUNDING
MATHEMATICA MONTISNIGRI Vol XXXII (2015) 23-30 Dedicated to the 80th anniversary of professor V. I. Gavrilov Dedicated to the 80th anniversary of professor V. I. Gavrilov ON PARALLEL ALGORITHMS FOR SOLVING
More informationNH 67, Karur Trichy Highways, Puliyur C.F, Karur District DEPARTMENT OF INFORMATION TECHNOLOGY DIGITAL SIGNAL PROCESSING UNIT 3
NH 67, Karur Trichy Highways, Puliyur C.F, 639 114 Karur District DEPARTMENT OF INFORMATION TECHNOLOGY DIGITAL SIGNAL PROCESSING UNIT 3 IIR FILTER DESIGN Structure of IIR System design of Discrete time
More informationPermutations with short monotone subsequences
Permutations with short monotone subsequences Dan Romik Abstract We consider permutations of 1, 2,..., n 2 whose longest monotone subsequence is of length n and are therefore extremal for the Erdős-Szekeres
More informationOrthogonal vs Non-Orthogonal Multiple Access with Finite Input Alphabet and Finite Bandwidth
Orthogonal vs Non-Orthogonal Multiple Access with Finite Input Alphabet and Finite Bandwidth J. Harshan Dept. of ECE, Indian Institute of Science Bangalore 56, India Email:harshan@ece.iisc.ernet.in B.
More informationInfinite Impulse Response Filters
6 Infinite Impulse Response Filters Ren Zhou In this chapter we introduce the analysis and design of infinite impulse response (IIR) digital filters that have the potential of sharp rolloffs (Tompkins
More informationBlind Blur Estimation Using Low Rank Approximation of Cepstrum
Blind Blur Estimation Using Low Rank Approximation of Cepstrum Adeel A. Bhutta and Hassan Foroosh School of Electrical Engineering and Computer Science, University of Central Florida, 4 Central Florida
More informationThree-Mirror Anastigmat Telescope with an Unvignetted Flat Focal Plane
Three-Mirror Anastigmat Telescope with an Unvignetted Flat Focal Plane arxiv:astro-ph/0504514v1 23 Apr 2005 Kyoji Nariai Department of Physics, Meisei University, Hino, Tokyo 191-8506 nariai.kyoji@gakushikai.jp
More informationA Fast-Marching Approach to Cardiac Electrophysiology Simulation for XMR Interventional Imaging
A Fast-Marching Approach to Cardiac Electrophysiology Simulation for XMR Interventional Imaging M. Sermesant 1, Y. Coudière 2,V.Moreau-Villéger 3,K.S.Rhode 1, D.L.G. Hill 4,, and R.S. Razavi 1 1 Division
More informationA Signal Space Theory of Interferences Cancellation Systems
A Signal Space Theory of Interferences Cancellation Systems Osamu Ichiyoshi Human Network for Better 21 Century E-mail: osamu-ichiyoshi@muf.biglobe.ne.jp Abstract Interferences among signals from different
More informationA Weighted Least Squares Algorithm for Passive Localization in Multipath Scenarios
A Weighted Least Squares Algorithm for Passive Localization in Multipath Scenarios Noha El Gemayel, Holger Jäkel, Friedrich K. Jondral Karlsruhe Institute of Technology, Germany, {noha.gemayel,holger.jaekel,friedrich.jondral}@kit.edu
More informationSolutions to the problems from Written assignment 2 Math 222 Winter 2015
Solutions to the problems from Written assignment 2 Math 222 Winter 2015 1. Determine if the following limits exist, and if a limit exists, find its value. x2 y (a) The limit of f(x, y) = x 4 as (x, y)
More informationREFLECTION AND TRANSMISSION OF LAMB WAVES AT DISCONTINUITY IN PLATE Z. Liu NDT Systems & Services AG, Stutensee, Germany
REFLECTION AND TRANSMISSION OF LAMB WAVES AT DISCONTINUITY IN PLATE Z. Liu NDT Systems & Services AG, Stutensee, Germany Abstract: Lamb waves can be used for testing thin plate and pipe because they provide
More information(i) Understanding the basic concepts of signal modeling, correlation, maximum likelihood estimation, least squares and iterative numerical methods
Tools and Applications Chapter Intended Learning Outcomes: (i) Understanding the basic concepts of signal modeling, correlation, maximum likelihood estimation, least squares and iterative numerical methods
More informationGroundwave Propagation, Part One
Groundwave Propagation, Part One 1 Planar Earth groundwave 2 Planar Earth groundwave example 3 Planar Earth elevated antenna effects Levis, Johnson, Teixeira (ESL/OSU) Radiowave Propagation August 17,
More informationMinimising Time-Stepping Errors in Numerical Models of the Atmosphere and Ocean
Minimising Time-Stepping Errors in Numerical Models of the Atmosphere and Ocean University of Reading School of Mathematics, Meteorology and Physics Robert J. Smith August 2010 This dissertation is submitted
More informationFast Placement Optimization of Power Supply Pads
Fast Placement Optimization of Power Supply Pads Yu Zhong Martin D. F. Wong Dept. of Electrical and Computer Engineering Dept. of Electrical and Computer Engineering Univ. of Illinois at Urbana-Champaign
More informationCitation for published version (APA): Nutma, T. A. (2010). Kac-Moody Symmetries and Gauged Supergravity Groningen: s.n.
University of Groningen Kac-Moody Symmetries and Gauged Supergravity Nutma, Teake IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please
More informationOptimization Techniques for Alphabet-Constrained Signal Design
Optimization Techniques for Alphabet-Constrained Signal Design Mojtaba Soltanalian Department of Electrical Engineering California Institute of Technology Stanford EE- ISL Mar. 2015 Optimization Techniques
More informationContinuity of the Norm of a Composition Operator
Integr. equ. oper. theory 45 (003) 35 358 0378-60X/03035-8 $.50+0.0/0 c 003 Birkhäuser Verlag Basel/Switzerl Integral Equations Operator Theory Continuity of the Norm of a Composition Operator David B.
More informationThe Casey angle. A Different Angle on Perspective
A Different Angle on Perspective Marc Frantz Marc Frantz (mfrantz@indiana.edu) majored in painting at the Herron School of Art, where he received his.f.a. in 1975. After a thirteen-year career as a painter
More informationA New Hybrid Multitoning Based on the Direct Binary Search
IMECS 28 19-21 March 28 Hong Kong A New Hybrid Multitoning Based on the Direct Binary Search Xia Zhuge Yuki Hirano and Koji Nakano Abstract Halftoning is an important task to convert a gray scale image
More informationSignal Analysis Using The Solitary Chirplet
Signal Analysis Using The Solitary Chirplet Sai Venkatesh Balasubramanian Sree Sai Vidhya Mandhir, Mallasandra, Bengaluru-560109, Karnataka, India saivenkateshbalasubramanian@gmail.com Abstract: In the
More informationTwo-Dimensional Wavelets with Complementary Filter Banks
Tendências em Matemática Aplicada e Computacional, 1, No. 1 (2000), 1-8. Sociedade Brasileira de Matemática Aplicada e Computacional. Two-Dimensional Wavelets with Complementary Filter Banks M.G. ALMEIDA
More informationQ(173)Q(177)Q(188)Q(193)Q(203)
MATH 313: SOLUTIONS HW3 Problem 1 (a) 30941 We use the Miller-Rabin test to check if it prime. We know that the smallest number which is a strong pseudoprime both base 2 and base 3 is 1373653; hence, if
More informationOn uniquely k-determined permutations
On uniquely k-determined permutations Sergey Avgustinovich and Sergey Kitaev 16th March 2007 Abstract Motivated by a new point of view to study occurrences of consecutive patterns in permutations, we introduce
More informationOPTIMAL EXCITATION FREQUENCY FOR DELAMINATION IDENTIFICATION OF LAMINATED BEAMS USING A 0 LAMB MODE
OPTIMAL EXCITATION FREQUENCY FOR DELAMINATION IDENTIFICATION OF LAMINATED BEAMS USING A 0 LAMB MODE N. Hu 1 *, H. Fukunaga 2, Y. Liu 3 and Y. Koshin 2 1 Department of Mechanical Engineering, Chiba University,
More informationName: ID: Section: Math 233 Exam 2. Page 1. This exam has 17 questions:
Page Name: ID: Section: This exam has 7 questions: 5 multiple choice questions worth 5 points each. 2 hand graded questions worth 25 points total. Important: No graphing calculators! Any non scientific
More information