Finding discontinuities of piecewise-smooth functions
|
|
- Angela Hall
- 6 years ago
- Views:
Transcription
1 Finding discontinuities of piecewise-smooth functions A.G. Ramm Mathematics Department, Kansas State University, Manhattan, KS , USA Abstract Formulas for stable differentiation of piecewise-smooth functions are given. The data are noisy values of these functions. The locations of discontinuity points the sizes of the jumps across these points are not assumed known, but found stably from the noisy data. 1 Introduction Let f be a piecewise-c 2 ([0, 1]) function, 0 < x 1 < x 2 < < x J, 1 j J, are discontinuity points of f. We do not assume their locations x j their number J known a priori. We assume that the limits f(x j ± 0) exist, sup f (m) (x) M m, m = 0, 1, 2. (1.1) x x j,1 j J Assume that f δ is given, f f δ := sup x xj,1 j J f(x) f δ (x) δ, where f δ L (0, 1) are the noisy data. The problem is: given {f δ, δ}, where δ (0, δ 0 ) δ 0 > 0 is a small number, estimate stably f, find the locations of discontinuity points x j of f their number J, estimate the jumps p j := f(x j + 0) f(x j 0) of f across x j, 1 j J. A stable estimate R δ f δ of f is an estimate satisfying the relation lim δ 0 R δ f δ f = 0. There is a large literature on stable differentiation of noisy smooth functions (e.g., see references in [3]), but the problem stated above was not solved for piecewise-smooth functions by the method given below. A statistical estimation of the location of discontinuity Math subject classification: 65D35; 65D05 key words: inequalities, stable differentiation, noisy data, discontinuities, jumps, signal processing, edge detection.
2 points from noisy discrete data is given in [1]. In [5], [7], [2], various approaches to finding discontinuities of functions from the measured values of these functions are developed. The following formula was proposed originally (in 1968, see [4], [3]) for stable estimation of f (x), assuming f C 2 ([0, 1]), M 2 0, given noisy data f δ : R δ f δ := f ( ) 1 δ(x + h(δ)) f δ (x h(δ)) 2δ 2, h(δ) :=, h(δ) x 1 h(δ), (1.2) 2h(δ) M 2 R δ f δ f 2M 2 δ := ε(δ), (1.3) where the norm in (1.3) is L (0, 1) norm. Numerical efficiency stability of the stable differentiation method proposed in [4] has been demonstrated in [6]. Moreover, (cf [3]), inf T sup T f δ f ε(δ), (1.4) f K(M 2,δ) where T : L (0, 1) L (0, 1) runs through the set of all bounded operators, K(M 2, δ) := {f : f M 2, f f δ δ}. Therefore estimate (1.2) is the best possible estimate of f, given noisy data f δ, assuming f K(M 2, δ). In [3] this result was generalized to the case f K(M a, δ), f (a) M a, 1 < a 2, f (x) f (x ) x x a 1 where f (a) := f + f +sup x,x, 1 < a 2, f (a) is the fractional-order derivative of f. The aim of this paper is to extend the above results to the case of piecewise-smooth functions. In Section 2 the results are formulated, proofs are given. In Section 3 the case of continuous piecewise-smooth functions is treated. 2 Formulation of the result Theorem 1. Formula (1.2) gives stable estimate of f on the set S δ := [h(δ), 1 h(δ)] \ J j=1 (x j h(δ), x j + h(δ)), (1.3) holds with the norm taken on the set S δ. Assuming M 2 > 0 computing the quantities f j := f δ(jh+h) f δ (jh h), where h := h(δ) := ( ) 1 2δ 2 M 2, 1 j < [ 1 ], for sufficiently small δ, one finds the location of h ] is the integer discontinuity points of f with accuracy 2h, their number J. Here [ 1 h smaller than 1 closest to 1. The discontinuity points of f are located on the intervals h h (jh h, jh + h) such that f j 1 for sufficiently small δ, where ε(δ) is defined in (1.3). The size p j of the jump of f across the discontinuity point x j is estimated by the formula p j f δ (jh + h) f δ (jh h), the error of this estimate is O( δ). Let us assume that min j p j := p h(δ), where means much greater than. Then x j is located on the j-th interval [jh h, jh + h], h := h(δ), such that f j := f δ (jh + h) f δ (jh h) 2h 1, (2.1) 2 2h
3 so that x j is localized with the accuracy 2h(δ). More precisely, f j f(jh+h) f(jh h) δ, 2h h δ = 0.5ε(δ), where ε(δ) is defined in (1.3). One has h f(jh+h) f(jh h) p j f(jh+h) f(x j +0) f(jh h) f(x j 0) p j 2M 1 h. Thus, f j p j 2h M 1 0.5ε(δ) = c 1 p j δ c 2 1, where c 1 := M 2 2 2, c 2 := M ε(δ). The jump p j is estimated by the formula: the error estimate of this formula can be given: p j [f δ (jh + h) f δ (jh h)], (2.2) p j [f δ (jh + h) f δ (jh h)] 2δ + 2M 1 h = 2δ + 2M 1 2δ M 2 = O( δ). (2.3) Thus, the error of the calculation of p j by the formula p j f δ (jh + h) f δ (jh h) is O(δ 1 2 ) as δ 0. Proof of Theorem 1. If x S δ, then, using Taylor s formula, one gets: (R δ f δ )(x) f (x) δ h + M 2h 2. (2.4) Here we assume that M 2 > 0 the interval (x h(δ), x + h(δ)) S δ, i.e., this interval does not contain discontinuity points of f. If, for all sufficiently small h, not necessarily for h = h(δ), inequality (2.4) fails, i.e., if (R δ f δ )(x) f (x) > δ + M 2h for all sufficiently h 2 small h > 0, then the interval (x h, x + h) contains a point x j S δ, i.e., a point of discontinuity of f or f. This observation can be used for locating the position of an isolated discontinuity point x j of f with any desired accuracy provided that the size p j of the jump of f across x j is greater than kδ, where k > 2 is a constant, p j > kδ, that h can be taken as small as desirable. Indeed, if x j (x h, x + h), then we have p j 2δ 2hM 1 f δ (x + h) f δ (x h) p j + 2δ + 2hM 1. The above estimate follows from the relation f δ (x+h) f δ (x h) = f(x+h) f(x j +0)+p j +f(x j 0) f(x h)±2δ = p j ±(2δ+2hM 1 ). Here p ± b, where b > 0, denotes a quantity such that p b p ± b p + b. Thus, if h is sufficiently small p j > kδ, where k > 2, then the inequality (k 2)δ 2hM 1 f δ (x + h) f δ (x h) 3
4 can be checked, therefore the inclusion x j (x h, x + h) can be checked. Since h > 0 is arbitrarily small in this argument, it follows that the location of the discontinuity point x j of f is established with arbitrary accuracy. A discussion of the case when a discontinuity point x j belongs to the interval (x h(δ), x + h(δ)) will be given below. Minimizing the right-h side of (2.4) with respect to h yields formula (1.2) for the minimizer h = h(δ) defined in (1.2), estimate (1.3) for the minimum of the right-h side of (2.4). If p h(δ), (2.1) holds, then the discontinuity points are located with the accuracy 2h(δ), as we prove now. Consider the case when a discontinuity point x j of f belongs to the interval (jh h, jh + h), where h = h(δ). Then estimate (2.2) can be obtained as follows. For jh h x j jh + h, one has f(x j + 0) f(x j 0) f δ (jh + h) + f δ (jh h) 2δ+ + f(x j + 0) f(jh + h) + f(x j 0) f(jh h) 2δ + 2hM 1, h = h(δ). This yields formulas (2.2) (2.3). Computing the quantities f j for 1 j < [ 1 ], h finding the intervals on which (2.1) holds for sufficiently small δ, one finds the location of discontinuity points of f with accuracy 2h, the number J of these points. For a small fixed δ > 0 the above method allows one to recover the discontinuity points of f at which f j p j δ M 2h h 1 1. This is the inequality (2.1). If h = h(δ), then δ = 0.5ε(δ) = O( δ), 2hf h j p j = O( δ) as δ 0 provided that M 2 > 0. Theorem 1 is proved. Remark 1: Similar results can be derived if f (a) L (S δ ) := f (a) Sδ M a, 1 < a 2. [ ] 1 In this case h = h(δ) = c a δ 1 2 a a, where c a = M a(a 1), R δ f δ is defined in (1.2), the error of the estimate is: R δ f δ f Sδ ( ) a 1 am 1 2 a a δ a 1 a. a 1 The proof is similar to that given in Section 3. It is proved in [3] that for C a -functions given with noise it is possible to construct stable differentiation formulas if a > 1 it is impossible to construct such formulas if a 1. The obtained formulas are useful in applications. One can also use L p -norm on S δ in the estimate f (a) Sδ M a (cf. [3]). Remark 2: The case when M 2 = 0 requires a special discussion. In this case the last term on the right-h side of formula (2.4) vanishes the minimization with respect to h becomes void: it requires that h be as large as possible, but one cannot take h arbitrarily large because estimate (2.4) is valid only on the interval (x h, x + h) which does not contain discontinuity points of f, these points are unknown. If M 2 = 0, then f is a piecewise-linear function. The discontinuity points of a piecewise-linear function 4
5 can be found if the sizes p j of the jumps of f across these points satisfy the inequality p j >> 2δ +2M 1 h for some choice of h. For instance, if h = δ M 1, then 2δ +2M 1 h = 4δ. So, if p j >> 4δ, then the location of discontinuity points of f can be found in the case when M 2 = 0. These points are located on the intervals for which f δ (jh+h) f δ (jh h) >> 4δ, where h = δ M 1. The size p j of the jump of f across a discontinuity point x j can be estimated by formula (2.2) with h = δ M 1, one assumes that x j (jh h, jh + h) is the only discontinuity point on this interval. The error of the formula (2.2) is estimated as in the proof of Theorem 1. This error is not more than 2δ + 2M 1 h = 4δ for the above choice of h = δ M 1. One can estimate the derivative of f at the point of smoothness of f assuming M 2 = 0 provided that this derivative is not too small. If M 2 = 0, then f = a j x + b j on every interval j between the discontinuity points x j, where a j b j are some constants. If (jh h, jh+h) j, f j := f δ(jh+h) f δ (jh h), then f 2h j a j δ tδ. Choose h = h M 1, where t > 0 is a parameter, M 1 = max j a j. Then the relative error of the approximate formula a j f j for the derivative f = a j on j equals to f j a j a j M 1 t a j. Thus, if, e.g., a j M 1 t = 20, then the relative error of the above approximate formula is not 2 more than Continuous piecewise-smooth functions Suppose now that ξ (mh h, mh + h), where m > 0 is an integer, ξ is a point at which f is continuous but f (ξ) does not exist. Thus, the jump of f across ξ is zero, but ξ is not a point of smoothness of f. How does one locate the point ξ? The algorithm we propose consists of the following. We assume that M 2 > 0 on S δ. Calculate the numbers f j := f δ(jh+h) f δ (jh h) f 2h j+1 f j, j = 1, 2,..., h = h(δ) = 2δ M 2. Inequality (1.3) implies f j ε(δ) f (jh) f j + ε(δ), where ε(δ) is defined in (1.3). Therefore, if f j > ε(δ), then sign f j = sign f (jh). One has: J 2δ h f j+1 f j J + 2δ where δ h = 0.5ε(δ) J := h, f(jh + 2h) f(jh) f(jh + h) + f(jh h). 2h Using Taylor s formula, one derives the estimate: 0.5[J 1 ε(δ)] J 0.5[J 1 + ε(δ)], (3.1) where J 1 := f (jh + h) f (jh). 5
6 If the interval (jh h, jh + 2h) belongs to S δ, then In this case J ε(δ), so J 1 = f (jh + h) f (jh) M 2 h = ε(δ). f j+1 f j 2ε(δ) if (jh h, jh + 2h) S δ. (3.2) Conclusion: If f j+1 f j > 2ε(δ), then the interval (jh h, jh + 2h) does not belong to S δ, that is, there is a point ξ (jh h, jh + 2h) at which the function f is not twice continuously differentiable with f M 2. Since we assume that either at a point ξ the function is twice differentiable, or at this point f does not exist, it follows that if f j+1 f j > 2ε(δ), then there is a point ξ (jh h, jh + 2h) at which f does not exist. If f j f j+1 < 0, (3.3) min( f j+1, f j ) > ε(δ), (3.4) then (3.3) implies f (jh)f (jh + h) < 0, so the interval (jh, jh + h) contains a critical point ξ of f, or a point ξ at which f does not exist. To determine which one of these two cases holds, let us use the right inequality (3.1). If ξ is a critical point of f ξ (jh, jh + h) S δ, then J 1 ε(δ), in this case the right inequality (3.1) yields f j+1 f j 2ε(δ). (3.5) Conclusion: If (3.3)-(3.5) hold, then ξ is a critical point. If (3.3) (3.4) hold f j+1 f j > 2ε(δ) then ξ is a point of discontinuity of f. If ξ is a point of discontinuity of f, we would like to estimate the jump Using Taylor s formula one gets P := f (ξ + 0) f (ξ 0). f j+1 f j = P 2 ± 2.5ε(δ). (3.6) The expression A = B ± b, b > 0, means that B b A B + b. Therefore, We have proved the following theorem: P = 2(f j+1 f j ) ± 5ε(δ). (3.7) Theorem 2. If ξ (jh h, jh + 2h) is a point of continuity of f f j+1 f j > 2ε(δ), then ξ is a point of discontinuity of f. If (3.3) (3.4) hold, f j+1 f j 2ε(δ), then ξ is a critical point of f. If (3.3) (3.4) hold f j+1 f j > 2ε(δ), then ξ (jh, jh + h) is a point of discontinuity of f. The jump P of f across ξ is estimated by formula (3.7). 6
7 4 Finding nonsmoothness points of piecewise-linear functions Assume that f is a piecewise-linear function on the interval [0, 1] 0 < x 1 <...x J < 1 is its nonsmoothness points, i.e, the discontinuity points of f or these of f. Assume that f δ is known at a grid mh, m = 0, 1, 2,..., M, h = 1, f M δ,m = f δ (mh), f(mh) f δ,m δ m, f m = f(mh). If mh is a discontinuity point, mh = x j, then we define its value as f(x j 0) or f(x j + 0), depending on which of these two numbers satisfy the inequalty f(mh) f δ,m δ. The problem is: Given f δ,m m, estimate the location of the discontinuity points x j, their number J, find out which of these points are points of discontinuity of f which are points of discontinuity of f but points of continuity of f, estimate the sizes of the jumps p j = f(x j + 0) f(x j 0) the sizes of the jumps q j = f (x j + 0) f (x j 0) at the continuity points of f which are discontinuity points of f. Let us solve this problem. Consider the quantities where G m := f δ,m+1 2f δ,m + f δ,m 1 2h 2 := g m + w m, g m := f m+1 2f m + f m 1 2h 2, w m := f δ,m+1 f m+1 2(f δ,m f m ) + f δ,m 1 f m 2h 2. We have Therefore, if min j x j+1 x j > 2h w m 4δ 2h 2 = 2δ h 2, g m = 0 if x j (mh h, mh + h) j. G m > 2δ h 2, (4.1) then the interval (mh h, mh+h) must contain a discontinuity point of f. This condition is sufficient for the interval (mh h, mh+h) to contain a discontinuity point of f, but not a necessary one: if the condition min j x j+1 x j > 2h is dropped, because then it may happen that the interval (mh h, mh + h) contains more than one discontinuity points without changing g m or G m, so that one cannot detect these points by the above method. We have proved the following result. Theorem 3. Condition (4.1) is a sufficient condition for the interval (mh h, mh+h) to contain a nonsmoothness point of f. If one knows a priori that x j+1 x j > 2h then condition (4.1) is a necessary sufficient condition for the interval (mh h, mh + h) to contain exactly one point of nonsmoothness of f. 7
8 Let us estimate the size of the jump p j. Let us assume that (4.1) holds, x j+1 x j > 2h x j (mh h, mh). The case when x j (mh, mh + h) is treated similarly. Let f(x) = a j x + b j when mh < x < x j, f(x) = a j+1 x + b j+1 when x j < x < (m + 1)h, where a j, b j are constants. One has Thus g m = (a j+1 a j )(mh h) (b j+1 b j ) 2h 2, p j = (a j+1 a j )x j + b j+1 b j. g m = (a j+1 a j )x j (b j+1 b j ) (a j+1 a j )(mh h x j ) 2h 2 = p j 2h ± a j+1 a j x j (mh h), 2 2h 2 where the symbol a ± b means a b a ± b a + b. The quantity a j+1 a j = q j, x j (mh h) h if mh h < x j < mh. Thus, G m = p ( j 2h ± qj h 2 2h + 2δ ), 2 h 2 provided that p j > 0. If then If p j = 0 then x j = b j b j+1 a j+1 a j, G m = p ( j 1 ± q ) jh + 4δ, 2h 2 p j q j h + 4δ p j << 1 p j > 0, p j 2h 2 G m. G m = q j(x j mh + h) 2h 2 because x j > mh h by the assumption. Thus, q j 2h G m, q j > 0. ± 2δ h 2, If min j x j+1 x j > 2h, then the number of the nonsmoothness points of f can be determined as the number of intervals on which (4.1) holds. 8
9 References [1] A. I. Katsevich, A. G. Ramm, Nonparametric estimation of the singularities of a signal from noisy measurements, Proc. AMS, 120, N8, (1994), [2] A. I. Katsevich, A. G. Ramm, Multidimensional algorithm for finding discontinuities of functions from noisy data. Math. Comp. Modelling, 18, N1, (1993), [3] A. G. Ramm, Inverse problems, Springer, New York, [4] A. G. Ramm, On numerical differentiation, Mathematics, Izvestija VUZOV, 11, (1968), [5] A. G. Ramm, A. I. Katsevich, The Radon Transform Local Tomography, CRC Press, Boca Raton, [6] A. G. Ramm, A. Smirnova, On stable numerical differentiation, Math. of Comput., 70, (2001), [7] A. G. Ramm, A.Zaslavsky, Reconstructing singularities of a function given its Radon transform, Math. Comp.Modelling, 18, N1, (1993),
444/,/,/,A.G.Ramm, On a new notion of regularizer, J.Phys A, 36, (2003),
444/,/,/,A.G.Ramm, On a new notion of regularizer, J.Phys A, 36, (2003), 2191-2195 1 On a new notion of regularizer A.G. Ramm LMA/CNRS, 31 Chemin Joseph Aiguier, Marseille 13402, France and Mathematics
More informationIn the book: Ten Mathematical Essays on Approximation in Analysis and Topology, Elsevier, Boston, 2005, pp (Eds J.Perrera, J.Lopez-Gomez, F.
In the book: Ten Mathematical Essays on Approximation in Analysis and Topology, Elsevier, Boston, 2005, pp.245-262 (Eds J.Perrera, J.Lopez-Gomez, F. Ruiz del Portal). 1 An essay on some problems of approximation
More informationComm. Nonlin. Sci. and Numer. Simul., 12, (2007),
Comm. Nonlin. Sci. and Numer. Simul., 12, (2007), 1390-1394. 1 A Schrödinger singular perturbation problem A.G. Ramm Mathematics Department, Kansas State University, Manhattan, KS 66506-2602, USA ramm@math.ksu.edu
More informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationA collocation method for solving some integral equations in distributions
A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA sapto@math.ksu.edu A G Ramm
More informationInverse scattering problem with underdetermined data.
Math. Methods in Natur. Phenom. (MMNP), 9, N5, (2014), 244-253. Inverse scattering problem with underdetermined data. A. G. Ramm Mathematics epartment, Kansas State University, Manhattan, KS 66506-2602,
More informationA numerical algorithm for solving 3D inverse scattering problem with non-over-determined data
A numerical algorithm for solving 3D inverse scattering problem with non-over-determined data Alexander Ramm, Cong Van Department of Mathematics, Kansas State University, Manhattan, KS 66506, USA ramm@math.ksu.edu;
More informationHow large is the class of operator equations solvable by a DSM Newton-type method?
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. How large
More informationA method for creating materials with a desired refraction coefficient
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. A method
More informationDynamical Systems Method for Solving Ill-conditioned Linear Algebraic Systems
Dynamical Systems Method for Solving Ill-conditioned Linear Algebraic Systems Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 6656-6, USA sapto@math.ksu.edu A G Ramm
More informationQ. Zou and A.G. Ramm. Department of Mathematics, Kansas State University, Manhattan, KS 66506, USA.
Computers and Math. with Applic., 21 (1991), 75-80 NUMERICAL SOLUTION OF SOME INVERSE SCATTERING PROBLEMS OF GEOPHYSICS Q. Zou and A.G. Ramm Department of Mathematics, Kansas State University, Manhattan,
More informationCubic Splines. Antony Jameson. Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305
Cubic Splines Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305 1 References on splines 1. J. H. Ahlberg, E. N. Nilson, J. H. Walsh. Theory of
More information2.4 The Precise Definition of a Limit
2.4 The Precise Definition of a Limit Reminders/Remarks: x 4 < 3 means that the distance between x and 4 is less than 3. In other words, x lies strictly between 1 and 7. So, x a < δ means that the distance
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationApplication of local operators for numerical reconstruction of the singular support of a vector field by its known ray transforms
Application of local operators for numerical reconstruction of the singular support of a vector field by its known ray transforms E Yu Derevtsov 1, V V Pickalov 2 and T Schuster 3 1 Sobolev Institute of
More informationDynamical Systems Method for Solving Operator Equations
Dynamical Systems Method for Solving Operator Equations Alexander G. Ramm Department of Mathematics Kansas State University Manhattan, KS 6652 email: ramm@math.ksu.edu URL: http://www.math.ksu.edu/ ramm
More informationAnn. Polon. Math., 95, N1,(2009),
Ann. Polon. Math., 95, N1,(29), 77-93. Email: nguyenhs@math.ksu.edu Corresponding author. Email: ramm@math.ksu.edu 1 Dynamical systems method for solving linear finite-rank operator equations N. S. Hoang
More informationA Lower Bound Theorem. Lin Hu.
American J. of Mathematics and Sciences Vol. 3, No -1,(January 014) Copyright Mind Reader Publications ISSN No: 50-310 A Lower Bound Theorem Department of Applied Mathematics, Beijing University of Technology,
More informationProof. We indicate by α, β (finite or not) the end-points of I and call
C.6 Continuous functions Pag. 111 Proof of Corollary 4.25 Corollary 4.25 Let f be continuous on the interval I and suppose it admits non-zero its (finite or infinite) that are different in sign for x tending
More informationPhys.Let A. 360, N1, (2006),
Phys.Let A. 360, N1, (006), -5. 1 Completeness of the set of scattering amplitudes A.G. Ramm Mathematics epartment, Kansas State University, Manhattan, KS 66506-60, USA ramm@math.ksu.edu Abstract Let f
More informationProjection Theorem 1
Projection Theorem 1 Cauchy-Schwarz Inequality Lemma. (Cauchy-Schwarz Inequality) For all x, y in an inner product space, [ xy, ] x y. Equality holds if and only if x y or y θ. Proof. If y θ, the inequality
More informationProblem 3. Give an example of a sequence of continuous functions on a compact domain converging pointwise but not uniformly to a continuous function
Problem 3. Give an example of a sequence of continuous functions on a compact domain converging pointwise but not uniformly to a continuous function Solution. If we does not need the pointwise limit of
More informationON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS
MATHEMATICS OF OPERATIONS RESEARCH Vol. 28, No. 4, November 2003, pp. 677 692 Printed in U.S.A. ON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS ALEXANDER SHAPIRO We discuss in this paper a class of nonsmooth
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationMath 5520 Homework 2 Solutions
Math 552 Homework 2 Solutions March, 26. Consider the function fx) = 2x ) 3 if x, 3x ) 2 if < x 2. Determine for which k there holds f H k, 2). Find D α f for α k. Solution. We show that k = 2. The formulas
More informationComments on the letter of P.Sabatier
Comments on the letter of P.Sabatier ALEXANDER G. RAMM Mathematics Department, Kansas State University, Manhattan, KS 66506-2602, USA ramm@math.ksu.edu Aug. 22, 2003. Comments on the letter of P.Sabatier,
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationAN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT
AN INVERSE PROBLEM FOR THE WAVE EQUATION WITH A TIME DEPENDENT COEFFICIENT Rakesh Department of Mathematics University of Delaware Newark, DE 19716 A.G.Ramm Department of Mathematics Kansas State University
More informationA Note on the Differential Equations with Distributional Coefficients
MATEMATIKA, 24, Jilid 2, Bil. 2, hlm. 115 124 c Jabatan Matematik, UTM. A Note on the Differential Equations with Distributional Coefficients Adem Kilicman Department of Mathematics, Institute for Mathematical
More informationDoes negative refraction make a perfect lens?
Does negative refraction make a perfect lens? arxiv:0809.080v1 [cond-mat.mtrl-sci] 1 Sep 008 A. G. Ramm Mathematics Department, Kansas State University, Manhattan, KS 66506-60, USA email: ramm@math.ksu.edu
More information0.1 Pointwise Convergence
2 General Notation 0.1 Pointwise Convergence Let {f k } k N be a sequence of functions on a set X, either complex-valued or extended real-valued. We say that f k converges pointwise to a function f if
More information1-Bit Matrix Completion
1-Bit Matrix Completion Mark A. Davenport School of Electrical and Computer Engineering Georgia Institute of Technology Yaniv Plan Mary Wootters Ewout van den Berg Matrix Completion d When is it possible
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More information1-Bit Matrix Completion
1-Bit Matrix Completion Mark A. Davenport School of Electrical and Computer Engineering Georgia Institute of Technology Yaniv Plan Mary Wootters Ewout van den Berg Matrix Completion d When is it possible
More information1-Bit Matrix Completion
1-Bit Matrix Completion Mark A. Davenport School of Electrical and Computer Engineering Georgia Institute of Technology Yaniv Plan Mary Wootters Ewout van den Berg Matrix Completion d When is it possible
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationDynamical systems method (DSM) for selfadjoint operators
Dynamical systems method (DSM) for selfadjoint operators A.G. Ramm Mathematics Department, Kansas State University, Manhattan, KS 6656-262, USA ramm@math.ksu.edu http://www.math.ksu.edu/ ramm Abstract
More informationElectromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient
168 Int. J. Computing Science and Mathematics, Vol. Electromagnetic wave scattering by many small bodies and creating materials with a desired refraction coefficient A.G. Ramm Mathematics epartment, Kansas
More informationSection 8.7. Taylor and MacLaurin Series. (1) Definitions, (2) Common Maclaurin Series, (3) Taylor Polynomials, (4) Applications.
Section 8.7 Taylor and MacLaurin Series (1) Definitions, (2) Common Maclaurin Series, (3) Taylor Polynomials, (4) Applications. MATH 126 (Section 8.7) Taylor and MacLaurin Series The University of Kansas
More informationIt follows from the above inequalities that for c C 1
3 Spaces L p 1. In this part we fix a measure space (, A, µ) (we may assume that µ is complete), and consider the A -measurable functions on it. 2. For f L 1 (, µ) set f 1 = f L 1 = f L 1 (,µ) = f dµ.
More informationHeat Transfer in a Medium in Which Many Small Particles Are Embedded
Math. Model. Nat. Phenom. Vol. 8, No., 23, pp. 93 99 DOI:.5/mmnp/2384 Heat Transfer in a Medium in Which Many Small Particles Are Embedded A. G. Ramm Department of Mathematics Kansas State University,
More informationIowa State University. Instructor: Alex Roitershtein Summer Homework #1. Solutions
Math 501 Iowa State University Introduction to Real Analysis Department of Mathematics Instructor: Alex Roitershtein Summer 015 EXERCISES FROM CHAPTER 1 Homework #1 Solutions The following version of the
More informationNumerical differentiation by means of Legendre polynomials in the presence of square summable noise
www.oeaw.ac.at Numerical differentiation by means of Legendre polynomials in the presence of square summable noise S. Lu, V. Naumova, S. Pereverzyev RICAM-Report 2012-15 www.ricam.oeaw.ac.at Numerical
More informationThe Laplace Transform
The Laplace Transform Laplace Transform Philippe B. Laval KSU Today Philippe B. Laval (KSU) Definition of the Laplace Transform Today 1 / 16 Outline General idea behind the Laplace transform and other
More informationAN ITERATIVE METHOD FOR COMPUTING ZEROS OF OPERATORS SATISFYING AUTONOMOUS DIFFERENTIAL EQUATIONS
Electronic Journal: Southwest Journal of Pure Applied Mathematics Internet: http://rattler.cameron.edu/swjpam.html ISBN 1083-0464 Issue 1 July 2004, pp. 48 53 Submitted: August 21, 2003. Published: July
More informationAPPROXIMATE ISOMETRIES ON FINITE-DIMENSIONAL NORMED SPACES
APPROXIMATE ISOMETRIES ON FINITE-DIMENSIONAL NORMED SPACES S. J. DILWORTH Abstract. Every ε-isometry u between real normed spaces of the same finite dimension which maps the origin to the origin may by
More informationComplex geometrical optics solutions for Lipschitz conductivities
Rev. Mat. Iberoamericana 19 (2003), 57 72 Complex geometrical optics solutions for Lipschitz conductivities Lassi Päivärinta, Alexander Panchenko and Gunther Uhlmann Abstract We prove the existence of
More informationA nonlinear singular perturbation problem
A nonlinear singular perturbation problem arxiv:math-ph/0405001v1 3 May 004 Let A.G. Ramm Mathematics epartment, Kansas State University, Manhattan, KS 66506-60, USA ramm@math.ksu.edu Abstract F(u ε )+ε(u
More informationFunctional Limits and Continuity
Chapter 4 Functional Limits and Continuity 4.1 Discussion: Examples of Dirichlet and Thomae Although it is common practice in calculus courses to discuss continuity before differentiation, historically
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematics ON A HYBRID FAMILY OF SUMMATION INTEGRAL TYPE OPERATORS VIJAY GUPTA AND ESRA ERKUŞ School of Applied Sciences Netaji Subhas Institute of Technology
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More informationInfinite Limits. By Tuesday J. Johnson
Infinite Limits By Tuesday J. Johnson Suggested Review Topics Algebra skills reviews suggested: Evaluating functions Graphing functions Working with inequalities Working with absolute values Trigonometric
More informationGeneral Theory of Large Deviations
Chapter 30 General Theory of Large Deviations A family of random variables follows the large deviations principle if the probability of the variables falling into bad sets, representing large deviations
More informationProblems for Chapter 3.
Problems for Chapter 3. Let A denote a nonempty set of reals. The complement of A, denoted by A, or A C is the set of all points not in A. We say that belongs to the interior of A, Int A, if there eists
More informationREGULARITY THROUGH APPROXIMATION FOR SCALAR CONSERVATION LAWS
SIAM J. MATH. ANAL. c 1988 Society for Industrial and Applied Mathematics Vol. 19, No. 4, pp. 1 XX, July 1988 003 REGULARITY THROUGH APPROXIMATION FOR SCALAR CONSERVATION LAWS BRADLEY J. LUCIER Abstract.
More informationSubdifferential representation of convex functions: refinements and applications
Subdifferential representation of convex functions: refinements and applications Joël Benoist & Aris Daniilidis Abstract Every lower semicontinuous convex function can be represented through its subdifferential
More informationThe Hilbert Transform and Fine Continuity
Irish Math. Soc. Bulletin 58 (2006), 8 9 8 The Hilbert Transform and Fine Continuity J. B. TWOMEY Abstract. It is shown that the Hilbert transform of a function having bounded variation in a finite interval
More informationProblem Set 5: Solutions Math 201A: Fall 2016
Problem Set 5: s Math 21A: Fall 216 Problem 1. Define f : [1, ) [1, ) by f(x) = x + 1/x. Show that f(x) f(y) < x y for all x, y [1, ) with x y, but f has no fixed point. Why doesn t this example contradict
More informationCESARO OPERATORS ON THE HARDY SPACES OF THE HALF-PLANE
CESARO OPERATORS ON THE HARDY SPACES OF THE HALF-PLANE ATHANASIOS G. ARVANITIDIS AND ARISTOMENIS G. SISKAKIS Abstract. In this article we study the Cesàro operator C(f)() = d, and its companion operator
More informationyou expect to encounter difficulties when trying to solve A x = b? 4. A composite quadrature rule has error associated with it in the following form
Qualifying exam for numerical analysis (Spring 2017) Show your work for full credit. If you are unable to solve some part, attempt the subsequent parts. 1. Consider the following finite difference: f (0)
More informationOptimal stopping time formulation of adaptive image filtering
Optimal stopping time formulation of adaptive image filtering I. Capuzzo Dolcetta, R. Ferretti 19.04.2000 Abstract This paper presents an approach to image filtering based on an optimal stopping time problem
More informationNormal splines in reconstruction of multi-dimensional dependencies
Normal splines in reconstruction of multi-dimensional dependencies V. K. GORBUNOV, Mathematical Economics Department, Ulyanovsk State University, Ulyanovsk, RUSSIA, K. S. MAKEDONSKY, Mathematical Economics
More informationCreating materials with a desired refraction coefficient: numerical experiments
Creating materials with a desired refraction coefficient: numerical experiments Sapto W. Indratno and Alexander G. Ramm Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA
More informationF (x) = P [X x[. DF1 F is nondecreasing. DF2 F is right-continuous
7: /4/ TOPIC Distribution functions their inverses This section develops properties of probability distribution functions their inverses Two main topics are the so-called probability integral transformation
More informationSYMMETRIC INTEGRALS DO NOT HAVE THE MARCINKIEWICZ PROPERTY
RESEARCH Real Analysis Exchange Vol. 21(2), 1995 96, pp. 510 520 V. A. Skvortsov, Department of Mathematics, Moscow State University, Moscow 119899, Russia B. S. Thomson, Department of Mathematics, Simon
More informationTHE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
Asian-European Journal of Mathematics Vol. 3, No. 1 (2010) 57 105 c World Scientific Publishing Company THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. Hoang
More informationNONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian UNDER NONHOMOGENEOUS NEUMANN BOUNDARY CONDITION
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 210, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian
More informationC.7. Numerical series. Pag. 147 Proof of the converging criteria for series. Theorem 5.29 (Comparison test) Let a k and b k be positive-term series
C.7 Numerical series Pag. 147 Proof of the converging criteria for series Theorem 5.29 (Comparison test) Let and be positive-term series such that 0, for any k 0. i) If the series converges, then also
More informationOptimal global rates of convergence for interpolation problems with random design
Optimal global rates of convergence for interpolation problems with random design Michael Kohler 1 and Adam Krzyżak 2, 1 Fachbereich Mathematik, Technische Universität Darmstadt, Schlossgartenstr. 7, 64289
More informationCollocation and iterated collocation methods for a class of weakly singular Volterra integral equations
Journal of Computational and Applied Mathematics 229 (29) 363 372 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationFunctional Analysis (2006) Homework assignment 2
Functional Analysis (26) Homework assignment 2 All students should solve the following problems: 1. Define T : C[, 1] C[, 1] by (T x)(t) = t x(s) ds. Prove that this is a bounded linear operator, and compute
More informationLocal strong convexity and local Lipschitz continuity of the gradient of convex functions
Local strong convexity and local Lipschitz continuity of the gradient of convex functions R. Goebel and R.T. Rockafellar May 23, 2007 Abstract. Given a pair of convex conjugate functions f and f, we investigate
More informationOn Discontinuous Differential Equations
On Discontinuous Differential Equations Alberto Bressan and Wen Shen S.I.S.S.A., Via Beirut 4, Trieste 3414 Italy. Department of Informatics, University of Oslo, P.O. Box 18 Blindern, N-316 Oslo, Norway.
More informationWEAK ASYMPTOTIC SOLUTION FOR A NON-STRICTLY HYPERBOLIC SYSTEM OF CONSERVATION LAWS-II
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 94, pp. 1 14. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu WEAK ASYMPTOTIC
More informationMETRIC HEIGHTS ON AN ABELIAN GROUP
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 44, Number 6, 2014 METRIC HEIGHTS ON AN ABELIAN GROUP CHARLES L. SAMUELS ABSTRACT. Suppose mα) denotes the Mahler measure of the non-zero algebraic number α.
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationThe Split Hierarchical Monotone Variational Inclusions Problems and Fixed Point Problems for Nonexpansive Semigroup
International Mathematical Forum, Vol. 11, 2016, no. 8, 395-408 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2016.6220 The Split Hierarchical Monotone Variational Inclusions Problems and
More informationALGEBRAIC GEOMETRY HOMEWORK 3
ALGEBRAIC GEOMETRY HOMEWORK 3 (1) Consider the curve Y 2 = X 2 (X + 1). (a) Sketch the curve. (b) Determine the singular point P on C. (c) For all lines through P, determine the intersection multiplicity
More informationIf you make reference to this version of the manuscript, use the following information:
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. Inverse
More informationDivision of the Humanities and Social Sciences. Supergradients. KC Border Fall 2001 v ::15.45
Division of the Humanities and Social Sciences Supergradients KC Border Fall 2001 1 The supergradient of a concave function There is a useful way to characterize the concavity of differentiable functions.
More information1 Functions of normal random variables
Tel Aviv University, 200 Gaussian measures : results formulated Functions of normal random variables It is of course impossible to even think the word Gaussian without immediately mentioning the most important
More informationA Simple Proof of the Fredholm Alternative and a Characterization of the Fredholm Operators
thus a n+1 = (2n + 1)a n /2(n + 1). We know that a 0 = π, and the remaining part follows by induction. Thus g(x, y) dx dy = 1 2 tanh 2n v cosh v dv Equations (4) and (5) give the desired result. Remarks.
More informationCHAPTER 3 Further properties of splines and B-splines
CHAPTER 3 Further properties of splines and B-splines In Chapter 2 we established some of the most elementary properties of B-splines. In this chapter our focus is on the question What kind of functions
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationThis article was published in an Elsevier journal. The attached copy is furnished to the author for non-commercial research and education use, including for instruction at the author s institution, sharing
More informationHow to Distinguish True Dependence from Varying Independence?
University of Texas at El Paso DigitalCommons@UTEP Departmental Technical Reports (CS) Department of Computer Science 8-1-2013 How to Distinguish True Dependence from Varying Independence? Marketa Krmelova
More informationSolutions to Tutorial 1 (Week 2)
THE UNIVERSITY OF SYDNEY SCHOOL OF MATHEMATICS AND STATISTICS Solutions to Tutorial 1 (Week 2 MATH3969: Measure Theory and Fourier Analysis (Advanced Semester 2, 2017 Web Page: http://sydney.edu.au/science/maths/u/ug/sm/math3969/
More informationBounded variation and Helly s selection theorem
Bounded variation and Helly s selection theorem Alexander P. Kreuzer ENS Lyon Universität der Bundeswehr München, October 10th, 2013 Outline 1 Functions of bounded variation Representation 2 Helly s selection
More informationat time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))
Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time
More informationAlgorithms for Nonsmooth Optimization
Algorithms for Nonsmooth Optimization Frank E. Curtis, Lehigh University presented at Center for Optimization and Statistical Learning, Northwestern University 2 March 2018 Algorithms for Nonsmooth Optimization
More informationA PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
IJMMS 25:6 2001) 397 409 PII. S0161171201002290 http://ijmms.hindawi.com Hindawi Publishing Corp. A PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More information212a1214Daniell s integration theory.
212a1214 Daniell s integration theory. October 30, 2014 Daniell s idea was to take the axiomatic properties of the integral as the starting point and develop integration for broader and broader classes
More informationA LITTLE REAL ANALYSIS AND TOPOLOGY
A LITTLE REAL ANALYSIS AND TOPOLOGY 1. NOTATION Before we begin some notational definitions are useful. (1) Z = {, 3, 2, 1, 0, 1, 2, 3, }is the set of integers. (2) Q = { a b : aεz, bεz {0}} is the set
More informationSOME PROPERTIES ON THE CLOSED SUBSETS IN BANACH SPACES
ARCHIVUM MATHEMATICUM (BRNO) Tomus 42 (2006), 167 174 SOME PROPERTIES ON THE CLOSED SUBSETS IN BANACH SPACES ABDELHAKIM MAADEN AND ABDELKADER STOUTI Abstract. It is shown that under natural assumptions,
More informationRECONSTRUCTION OF NUMERICAL DERIVATIVES FROM SCATTERED NOISY DATA. 1. Introduction
RECONSTRUCTION OF NUMERICAL DERIVATIVES FROM SCATTERED NOISY DATA T. WEI, Y. C. HON, AND Y. B. WANG Abstract. Based on the thin plate spline approximation theory, we propose in this paper an efficient
More informationSome Background Math Notes on Limsups, Sets, and Convexity
EE599 STOCHASTIC NETWORK OPTIMIZATION, MICHAEL J. NEELY, FALL 2008 1 Some Background Math Notes on Limsups, Sets, and Convexity I. LIMITS Let f(t) be a real valued function of time. Suppose f(t) converges
More informationComputational method for acoustic wave focusing
Int. J. Computing Science and Mathematics Vol. 1 No. 1 7 1 Computational method for acoustic wave focusing A.G. Ramm* Mathematics epartment Kansas State University Manhattan KS 6656-6 USA E-mail: ramm@math.ksu.edu
More information