Research Article Green s Theorem for Sign Data

Similar documents
The Order Relation and Trace Inequalities for. Hermitian Operators

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Research Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations

Case Study of Markov Chains Ray-Knight Compactification

Sharp integral inequalities involving high-order partial derivatives. Journal Of Inequalities And Applications, 2008, v. 2008, article no.

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Research Article Global Sufficient Optimality Conditions for a Special Cubic Minimization Problem

U-Pb Geochronology Practical: Background

Numerical Heat and Mass Transfer

Research Article Relative Smooth Topological Spaces

Research Article An Extension of Stolarsky Means to the Multivariable Case

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Existence of Two Conjugate Classes of A 5 within S 6. by Use of Character Table of S 6

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

Statistical inference for generalized Pareto distribution based on progressive Type-II censored data with random removals

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Error Probability for M Signals

One-sided finite-difference approximations suitable for use with Richardson extrapolation

A Solution of Porous Media Equation

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

Pulse Coded Modulation

More metrics on cartesian products

For now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.

11 Tail Inequalities Markov s Inequality. Lecture 11: Tail Inequalities [Fa 13]

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

On the spectral norm of r-circulant matrices with the Pell and Pell-Lucas numbers

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

A Robust Method for Calculating the Correlation Coefficient

Lecture 12: Discrete Laplacian

Composite Hypotheses testing

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Quality Assessment of Restored Satellite Data. Based on Signal to Noise Ratio

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

An Application of Fuzzy Hypotheses Testing in Radar Detection

Statistics Chapter 4

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Finite Element Modelling of truss/cable structures

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning

A random variable is a function which associates a real number to each element of the sample space

PHYS 705: Classical Mechanics. Calculus of Variations II

Mathematical Preparations

The internal structure of natural numbers and one method for the definition of large prime numbers

Research Article A Combination Method of Mixed Multiscale Finite-Element and Laplace Transform for Flow in a Dual-Permeability System

Multivariate Ratio Estimator of the Population Total under Stratified Random Sampling

International Journal of Engineering Research and Modern Education (IJERME) Impact Factor: 7.018, ISSN (Online): (

Research Article On the Open Problem Related to Rank Equalities for the Sum of Finitely Many Idempotent Matrices and Its Applications

), it produces a response (output function g (x)

Formal solvers of the RT equation

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

A Note on Bound for Jensen-Shannon Divergence by Jeffreys

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Electronic Quantum Monte Carlo Calculations of Energies and Atomic Forces for Diatomic and Polyatomic Molecules

NUMERICAL DIFFERENTIATION

Chapter 13: Multiple Regression

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

Quantum Particle Motion in Physical Space

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

From Biot-Savart Law to Divergence of B (1)

Integrals and Invariants of Euler-Lagrange Equations

CHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD

A Hybrid Variational Iteration Method for Blasius Equation

T E C O L O T E R E S E A R C H, I N C.

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

On an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1

Voting Games with Positive Weights and. Dummy Players: Facts and Theory

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

On the correction of the h-index for career length

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

Differentiating Gaussian Processes

An Improved multiple fractal algorithm

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

A new Approach for Solving Linear Ordinary Differential Equations

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

Economics 130. Lecture 4 Simple Linear Regression Continued

ECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics

EEE 241: Linear Systems

Introduction to Random Variables

Lecture 4 Hypothesis Testing

Chapter Newton s Method

Feature Selection: Part 1

Difference Equations

arxiv: v1 [math.ho] 18 May 2008

Lossy Compression. Compromise accuracy of reconstruction for increased compression.

International Journal of Algebra, Vol. 8, 2014, no. 5, HIKARI Ltd,

4DVAR, according to the name, is a four-dimensional variational method.

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

The Synchronous 8th-Order Differential Attack on 12 Rounds of the Block Cipher HyRAL

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models

Limited Dependent Variables

Solution Thermodynamics

P R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering /

Transcription:

Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of Lousana at Lafayette, Lafayette, LA 70504-4210, USA Correspondence should be addressed to Lous M. Houston, houston@lousana.edu Receved 14 March 2012; Accepted 19 Aprl 2012 Academc Edtors: A. Bellouqud and M.-H. Hsu Copyrght q 2012 Lous M. Houston. Ths s an open access artcle dstrbuted under the Creatve Commons Attrbuton Lcense, whch permts unrestrcted use, dstrbuton, and reproducton n any medum, provded the orgnal work s properly cted. Sgn data are the sgns of sgnal added to nose. It s well known that a constant sgnal can be recovered from sgn data. In ths paper, we show that an ntegral over varant sgnal can be recovered from an ntegral over sgn data based on the varant sgnal. We refer to ths as a generalzed sgn data average. We use ths result to derve a Green s theorem for sgn data. Green s theorem s mportant to varous sesmc processng methods, ncludng sesmc mgraton. Results n ths paper generalze reported results for 2.5D data volumes n whch Green s theorem apples to sgn data based only on tradtonal sgn data recovery. 1. Introducton In certan cases, data, consstng of coherent sgnal and random nose, s summed repeatedly n order to enhance the sgnal and reduce the nose. It has been found that when the sgnal-tonose rato s between 0.1 and 1, f only the sgns of the data are retaned pror to summaton, then the sgnal can be recovered 1. We refer to ths type of data as sgn data. Sgn data has advantages n sgnfcantly reducng the amount of space needed to record the data. The reducton n space can amount to a rato of 20 to 1 bts of data storage 1. In the sesmc ndustry, t has been dscovered that most of the data processng methods used on regular data are also effectve on sgn data 2. Many of these data processng methods are dependent on Green s theorem. In partcular, Krchhoff mgraton s dependent on Green s theorem. Mgraton s a sesmc process that corrects data coordnates whch are dstorted by attrbutes of the sesmc experment. In accordance wth Green s theorem, we show that an ntegral of an operaton on sgn data over a volume s equal to an ntegral of an operaton of the orgnal data over a surface. Ths result generalzes the work of Houston and Rchard 3 n whch, based only on tradtonal sgn data recovery, t was shown that Green s theorem s satsfed when the sgn data encompasses a 2.5D volume. 2.5D data s data whch s three-dmensonal but only has

2 ISRN Appled Mathematcs 12 10 8 sgn (sgnal + unose) sgnal + unose 6 4 2 sgn(sgnal + unose) sgnal + unose 0 sgnal 2 40 30 20 10 0 10 20 30 40 x Fgure 1: A synthetc sgnal sn c x/3 e x/15 2 and the computer-generated unform random nose used to examne sgn-bt ampltude recovery. Ths test has nose of unt magntude a 1 and a sgnal-to-nose rato of one. Shown from bottom to top s the sgnal, the sgnal plus nose, the sgn of the sgnal plus nose, the average over 200 teratons of sgnal plus nose, and the average over 200 teratons of the sgn of sgnal plus nose. varatons along two dmensons. In ths paper, we fnd that Green s theorem s satsfed for sgn data when the data volume s largely symmetrc. For cases n whch the data volume s of arbtrary shape we have derved a varant of Green s theorem. 2. A Generalzed Sgn Data Average It was shown by O Bren et al. 1 that there s a recovery of ampltude from sgn data n the presence of unform random nose. If the random nose X has ampltude a, the sgnal s desgnated by the functon f k, and the number of teratons s M, then we can wrte f k a sgn ( ) f k X, M 2.1 where sgn ( y ) 1, y > 0 0, y 0 1, y < 0. 2.2 Fgure 1 llustrates the data recovery process ncludng sgn data.

ISRN Appled Mathematcs 3 Multply both sdes of 2.1 by the functon g k : g k f k a g k sgn ( ) f k X. M 2.3 Now sum both sdes of 2.3 over the k ndex: g k f k a g k sgn ( ) f k X. M k k 2.4 If we allow k, then 2.4 becomes g f a g sgn ( ) f X M 2.5 or g f a g sgn ( f X ). 2.6 It s clear that f f f and g g, then 2.6 becomes f a sgn ( ) f X, M 2.7 whch s essentally 2.1. A contnuous verson of 2.6 mght be wrtten as g v f v dv a g v sgn ( f v v ) dv. 2.8 An argument for the consstency of 2.8 s as follows. Let f v f. Then we have f a g v sgn ( f v ) dv g v dv. 2.9 Now ntegrate over all values of v: f a g v sgn( f v ) dv g v dv. 2.10 Let g v be a unform probablty densty, ρ v.

4 ISRN Appled Mathematcs That mples g v dv 1, 2.11 and 2.10 becomes f a ρ v sgn ( f v ) dv, 2.12 whch was shown by Houston et al. 4. 3. A Generalzed Average for Sgn Data Dervatves Consder the followng nth order forward fnte dfference 5 : n ( ) Δ n h f v 1 n f v n h, 0 3.1 wth n n!/! n!. If we make the varable, v dscrete by choosng a small real nterval, q and wrtng f f ( q ), 3.2 where s an nteger ndex, then 3.1 becomes Δ n h f n ( ) 1 n f ( q n h ). 0 3.3 Because the fnte dfference s a lnear operator, we can use 2.6 to derve g Δ n h f a g Δ n h sgn( f X ). 3.4 Equaton 3.4 takes nto account the fact that g ( q ) f ( q lh ) a g ( q ) sgn ( f ( q lh ) X ), 3.5 where l s an arbtrary ndex.

ISRN Appled Mathematcs 5 Snce d n lm dvn h 0 Δ n h h n, 3.6 3.4 suggests the followng equaton: g v dn f v dv a g v dn dvn dv n sgn( f v v ) dv. 3.7 Let d n /dv n d n /du n and f v f u. Then from 3.7 we have d n a g v d n /du n sgn ( f u v ) dv f u. 3.8 dun g v dv Once agan, ntegrate over all v,letg v ρ v,and 3.8 becomes d n f u a ρ v sgn ( f u v ) dv, dun 3.9 whch was shown by Houston et al. 4. 4. The Applcaton to Green s Theorem Equaton 3.7 mples the specal case: g v d2 f v dv a g v d2 dv2 dv sgn( f v v ) dv. 2 4.1 Employng three varables n 4.1 yelds g v 1,v 2,v 3 2 a g v 1,v 2,v 3 2 f v 1,v 2,v 3 dv 1 dv 2 dv 3 sgn ( f v 1,v 2,v 3 v ) dv1 dv 2 dv 3, 4.2 whch can be smplfed f g g v 1,v 2,v 3, f f v 1,v 2,v 3,anddV dv 1 dv 2 dv 3 : g 2 fdv a g 2 sgn ( f v ) dv. 4.3

6 ISRN Appled Mathematcs Summaton over the varables yelds g 2 fdv a g 2 sgn ( f v ) dv. 4.4 Consequently, 4.4 can be wrtten as g 2 fdv a g 2 sgn ( f v ) dv. 4.5 If we make the mappng g v 2 g, v 2 4.6 then 2.8 becomes f 2 gdv a sgn ( ) 2 f v gdv. 4.7 Summaton over the varables yelds f 2 gdv a sgn ( ) 2 f v gdv, 4.8 whch can be wrtten as f 2 gdv a sgn ( ) 2 f v gdv. 4.9 Subtractng 4.9 from 4.5 yelds ( ( g 2 f f g) 2 dv a g 2 sgn ( ) f v sgn ( ) ) 2 f v gdv. 4.10 Green s theorem s ( S U 1 U 2 n U 2 ) U ) 1 ds (U 1 2 U 2 U 2 2 U 1 dv. n V Consequently, we can wrte a varant of Green s theorem for sgn data as S ( g f n f g ) ( ds a g n 2 sgn ( ) f v sgn ( ) ) 2 f v g dv. 4.11 4.12

ISRN Appled Mathematcs 7 If the spatal ntervals are unform, that s, v 1 a 1,b 1, v 2 a 2,b 2, v 3 a 3,b 3 and a 1,b 1 a 2,b 2 a 3,b 3, then 4.12 becomes S ( g f n f g ) ( ds a g 2 sgn ( ) ( ) ) f v sgn f v 2 g dv. n 4.13 In the case of 4.13, we see that Green s theorem s satsfed by replacng the functon f wth sgn data and dvdng the surface ntegral by the nose ampltude. Therefore, all processng of real data based on Green s theorem wll be effectve for sgn data wth an assocated varance. We also note that the effectveness of Green s theorem on sgn data s enhanced by operatng over a symmetrc volume. When the volume s not symmetrc, we can apply the varant of Green s theorem gven by 4.12. 5. The 2.5D Case Now let us consder the specal case for whch the functons f and g have only twodmensonal varaton. That s, f f v 1,v 2 g g v 1,v 2. 5.1 Wthout usng a generalzed sgn data average, ths mples ( a g 2 sgn ( ) ( ) ) ( ) f v sgn f v 2 g dv b 3 a 3 g 2 f f 2 g dv 1 dv 2 5.2 or S ( g f n f g ) ( ) ds b 3 a 3 g 2 f f 2 g dv 1 dv 2. n 5.3 Equaton 5.3 s Green s theorem for sgn data when the data encompasses a 2.5D volume and s consstent wth results reported n Houston and Rchard 3. 6. Computatonal Tests Equaton 2.6 derves from a weghted average of the functon y sgn f X. We can thus wrte the expectaton value of y as E Y g y g. 6.1

8 ISRN Appled Mathematcs We fnd that E Y g f a g. 6.2 The varance, Var Y, can be wrtten as Var Y E Y E Y 2. 6.3 Ths can be reduced to a smpler form by usng the fact that ( E Y 2) 1. 6.4 Thus, we have Var Y 1 E Y 2. 6.5 We can demonstrate 2.6, g f a g sgn ( f X ), 6.6 computatonally. Let f sn ( q ), g cos ( q ). 6.7 Consequently, we want to demonstrate that cos ( q ) sn ( q ) a cos ( q ) sgn ( sn ( q ) X ). 6.8 We can compute the average percentage error β as ( cos( q ) sn ( q ) a cos( q ) sgn ( sn ( q ) ) ) 2 X β 100 ( cos( q ) sn ( q )). 6.9 2

ISRN Appled Mathematcs 9 Table 1: Twenty trals whch compare the values of the ndcated expressons for M teratons per tral. The trals are dvded nto ten trals for two dfferent nose ampltudes. Recall that M.LetM 10000 and q 0.0001. M 1 cos q sn q a M 1 cos q sgn sn q X a 1 3540.594404 3529.425407 1 3569.469208 2 3460.5978 3 3540.394914 4 3542.955134 5 3644.993387 6 3441.978126 7 3510.495952 8 3544.098984 9 3645.64803 10 β 0.56447%, Var Y 0.822949 a 2 3540.594404 3578.540737 1 3468.32837 2 3063.26568 3 3629.793227 4 3562.692845 5 3605.376332 6 3284.005063 7 3375.67384 8 3356.467161 9 3500.608956 10 β 1.731154%, Var Y 0.955737 We should see a correlaton between β and Var Y. The results of ths comparson are shown n Table 1. 7. Conclusons Usng the results of sgn data sgnal recovery leads to a dervaton of a generalzed sgn data average. Extendng these results to ncorporate dervatves leads to a varant of Green s theorem for sgn data. We fnd that Green s theorem drectly apples to sgn data when the data volume s symmetrc and the surface ntegral s dvded by the nose ampltude. A specfc applcaton of ths result s that Green s theorem apples to sgn data when the data volume s 2.5D and a generalzed sgn data average s not requred. Acknowledgment Dscussons wth Maxwell Lueckenhoff are apprecated.

10 ISRN Appled Mathematcs References 1 J. T. O Bren, W. P. Kamp, and G. M. Hoover, Sgn-bt ampltude recovery wth applcatons to sesmc data, Geophyscs, vol. 47, no. 11, pp. 1527 1539, 1982. 2 N. A. Anstey, Sesmc Prospectng Instruments, Gebruder Borntraeger, Berln, Germany, 2nd edton, 1981. 3 L. M. Houston and B. A. Rchard, The Helmholtz-Krchoff 2.5D ntegral theorem for sgn-bt data, Geophyscs and Engneerng, vol. 1, no. 1, pp. 84 87, 2004. 4 L. M. Houston, G. A. Glass, and A. D. Dymnkov, Sgn data dervatve recovery, ISRN Appled Mathematcs, vol. 2012, Artcle ID 630702, 7 pages, 2012. 5 W. G. Kelley and A. C. Peterson, Dfference Equatons: an Introducton wth Applcatons, Academc Press, San Dego, Calf, USA, 1st edton, 1991.

Advances n Operatons Research Advances n Decson Scences Appled Mathematcs Algebra Probablty and Statstcs The Scentfc World Journal Internatonal Dfferental Equatons Submt your manuscrpts at Internatonal Advances n Combnatorcs Mathematcal Physcs Complex Analyss Internatonal Mathematcs and Mathematcal Scences Mathematcal Problems n Engneerng Mathematcs Dscrete Mathematcs Dscrete Dynamcs n Nature and Socety Functon Spaces Abstract and Appled Analyss Internatonal Stochastc Analyss Optmzaton