RECONSTRUCTION OF NUMERICAL DERIVATIVES FROM SCATTERED NOISY DATA. 1. Introduction
|
|
- Justina Russell
- 5 years ago
- Views:
Transcription
1 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 regularization algorithm for the reconstruction of numerical derivatives from two-dimensional scattered noisy data. An error estimation that deduces a good regularization parameter is given. Numerical results show that the proposed method is efficient and stable.. Introduction The problem to reconstruct numerical derivatives from noisy observational data arises from several research disciplines, for instances, image processing [5], Volterra integral equation [4, 6], analyzing photoelectric response data [4] and identification [8]. This numerical differentiation problem is well known to be ill-posed in the sense that a small noise in measurement data can induce a huge error in the approximated derivatives. There are some computational methods proposed recently for one dimensional case [, 7, 5, 6, 7] but so far only a very few results on higher dimensional cases have been reported []. These papers fall into three categories: difference methods [, 7, 5, 7], interpolation methods [6] and regularization methods [, 9, 9]. The finite difference method is simple and effective for precise data. For noisy data, the chosen stepsize, as a regularization parameter, must depend on the level of noise. This restricts the choice of specified data points used in the computation. Lagrangian numerical differentiation technique is optimally stable [6] but suffers from the limited choice on the interpolation nodes in order to minimize the error. To solve effectively the ill-posed problems, regularization methods play an important role [8]. Most of the regularization procedures for Key words and phrases. Numerical derivatives, Thin plate spline, Tikhonov regularization. The work described in this paper was partially supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. CityU 78/P). The third author was also supported by the NSF of China (7).
2 T. WEI, Y. C. HON, AND Y. B. WANG numerical differentiation [, 9] make use of the variational approach to solve a regularized functional. Once the regularization parameter is chosen, the approximation to the unknown function and its derivatives are then obtained. Unfortunately, the determination of optimal parameter is generally a nontrivial task. Recently, Wang et. al [9] gave an effective method for treating the numerical differentiation by Tikhonov regularization in one dimension. With a simple strategy to choose the regularization parameter, the numerical results are computed much faster than using the discrepancy principle in [9]. Based on their work we extend the result to two-dimensional case and give an error estimation from which a good regularization parameter can be devised. Following the idea given in [, ], we use the thin plate spline approximation to obtain the numerical derivatives for noisy scattered data in two- dimension. The minimizer of the objective functional can be expressed in terms of a linear combination of thin plate spline functions and a linear polynomial (see [,,, ] for details). Furthermore, we apply the idea given in [9] and the results obtained in [] to give an error estimation for computing the numerical derivatives from which a priori strategy for choosing a good regularization parameter can be derived. Numerical results for various cases show that this selection method for the regularization parameter is effective and stable.. Formulation of problem and solution Let y = y(x) be a function in R. Suppose that the noisy scattered data ỹ i are collected at points x i R, i =,,, n. The problem is to find a function such that its partial derivatives approximate the exact partial derivatives of the function y(x). Let {x, x,, x n } be a finite set of distinct points in R which are not collinear and Ω be its convex hull. Suppose T h be a triangular mesh of Ω with nodes x, x,, x n and denote h = max T T h diam(t ). Let s i be the total area of the triangles that connect with node i and s be the area of the domain Ω. Suppose that the given noisy samples ỹ i of the values y(x i ) satisfy ỹ i y(x i ) σ i δ, i =,,, n,
3 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA where δ is a positive constant denoting the maximum noise level. Let I e = { i σ i = } be an index set of points with exact measurement data and I a = {,,, n} \ I e be an index set of the noisy data. We aim at finding an approximate function f (x) of y(x) so that f y L (Ω) converges to zero as δ and h tend to zero. Denote x = (ξ, η) and define a functional ( ) () R(f, g) = ( f g ξ )( ξ ) + ( f ξ η )( g ξ η ) + f ( η )( g η ) dξdη and R () I(f) = R(f, f). Define X = {f f C (R ), D f L (R ), Df < } where Df and D f denote the first and the second derivatives of f respectively; and Y = {f f X, f(x j ) = y j, j I e }. Analogous to the approach in [, ], the problem is changed to seek a solution f α in space Y that minimizes the following functional () Φ(f) = ( ) s i f(xi ) ỹ i + αi(f), s σ i i I a where ασ i are regularization parameters in the sense of Tikhonov [8]. In this paper, we consider the following two problems: () How to find f α Y such that Φ(f α ) Φ(f) for all f Y? () Assume that such f α has been obtained. How to choose a good parameter α so that f α (x) converges to the exact function y(x) at a convergence order as high as possible? The first problem has been solved in [, ] by using the reproducing kernel theory in Hilbert space. In this section we give a simpler proof by using mathematical analysis. For the second problem, we give a new convergence order estimate in the next section.
4 4 T. WEI, Y. C. HON, AND Y. B. WANG Theorem.. Suppose f α is a linear combination of thin plate splines and a linear polynomial: n (4) f α (x) = c j x x j log x x j + d + d ξ + d η, j= where x = (ξ, η), x x j is Euclidean distance in R and {c j } n, {d j } are constants satisfying the following linear system of equations (5) f α (x i ) + 4πα s s i σ i c i = ỹ i, i I a, (6) f α (x i ) = y(x i ), i I e, n n (7) c j = c j x j =. j= j= Then f α is the unique minimizer of the functional (), i.e. Φ(f α ) = min f Y Φ(f) Φ(y). We first need the following Lemmas. Lemma.. Let g be any function of the form (4) whose coefficients satisfy equation (7). Then, for sufficiently large x, g satisfies (8) g(x) c x and g(x) c, where c and c are positive constants. Furthermore, for sufficiently large x, the moduli of all second and third derivatives of g are bounded above by c x and c 4 x, respectively, where c and c 4 are also positive constants. Proof. See Powell s paper []. Lemma.. Let f be any function of the form (4) whose coefficients satisfy (7). Then for any g X, we have n (9) R(f, g) = 8π c j g(x j ). Proof. From Lemma., we know that the second derivative of f is square integrable over the part of R which is outside the circle {x : x M }, where M is a sufficiently large constant. Since f is infinitely differentiable on any region of R that excludes the points x j, j =,,, n, its second derivatives are square integrable on those regions which are closed and bounded. Since the second derivative j=
5 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA 5 of f is O(log x x i ) near the point x i, it is square integrable on any bounded region including point x i. It is then easy to see that f C (R ), and hence f X and R(f, g) is well defined. Since f is a function of the form (4), f = 8π n j= c jδ(x x j ) in the sense of distribution. That is, for every ϕ(x) C (R ), we have () R(f, ϕ) = ( f, ϕ) L = 8π n c j ϕ(x j ). From the results on Sobolev space [], we know that H (R ) C (R ) and j= C (R ) is dense in H (R ). Since the second derivative of f is also square integrable, equation () holds for any function in H (R ). Let ψ C (R ) satisfy ψ(x) = if x and ψ(x) = if x. Suppose that D α ψ M for any multi-index α whose degree α and M is a constant. For any g X, let g m = g(x)ψ m (x) where ψ m = ψ( x m ), m is a positive number, then g m H (R ). Thus we have () R(f, g m ) = 8π n c j g m (x j ), for m =,,. j= It is now clear that for sufficiently large m, n j= c jg m (x j ) = n j= c jg(x j ). In the following, we show that R(f, g m ) converges to R(f, g) as m tends to. Since () R(f, g m ) ( ) = ( f g R ξ )( ξ ) + ( f ξ η )( g ξ η ) + f ( η )( g η ) ( + ( f ψ m ξ )( ξ ) + ( f ψ m ξ η )( ξ η ) + f ψ m ( η )( η ) + R R ( ( f ξ )( g +( f g )( η η )( ψ m η ) ξ )( ψ m ξ ) + ( f ξ η )( g ξ ) dξdη. ψ m η + g ψ m η ξ ) ψ m dξdη ) gdξdη It suffices to show that the first part of integrations in the right-hand side of () converges to R(f, g) as m tends to and the second and third parts both converge to. Note that for any multi-index β, β, () D β ψ m = ( m ) β D β ψ( x m ) M( m ) β.
6 6 T. WEI, Y. C. HON, AND Y. B. WANG Each integrand in the first part of integrations in the right-hand side of () is bounded by M D f(x)d g(x) and hence is integrable in R. Using the dominant convergence theorem, we know that this first part of integrations converges to R(f, g) as m. (4) For sufficiently large m, from Lemma. and (), we have R f(x) ξ g(x) ξ ψ m (x) ξ dx CM m m x m x dx, as m. All the other terms in the third part of integration of the right-hand side of () have the same convergence. The bound of g in X implies that g(x) is at most a constant multiple of x when x is large. Similarly, we have (5) R f(x) ψ m ξ ξ g dx CM m m x m x x dx, as m. All the remained terms in the second integration of the right-hand side of () have the limit zero. Therefore, by letting m tend to in equation (), we then obtain (6) R(f, g) = 8π n c j g(x j ). Proof of Theorem.: For any f Y and f α given by (4)-(7), denote function j= g = f f α. From Lemma., we have n (7) R(f α, g) = 8π c i g(x i ) i= = 8π i I a c i g(x i ) = s i α s i I a (f α (x i ) ỹ i )g(x i ) σi. Note that Φ(f) Φ(f α ) = ( ) s i f(xi ) f α (x i ) + α (I(f) I(f α )) s σ i i I a + ( ) s i g(xi )(f α (x i ) ỹ i ) s σ i I a i = I + αi,
7 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA 7 and I = I(f) I(f α ) = I(f f α ) + R(f α, g) = I(f f α ) ( (fα (x i ) ỹ i )g(x i ) α i I a s i s σ i ). Thus, we have (8) Φ(f) Φ(f α ) = ( ) s i f(xi ) f α (x i ) + αi(f f α ), s σ i i I a which indicates that f α is a minimizer of Φ. Furthermore, if we assume that there exists another function f Y such that Φ(f ) = Φ(f α ), then substituting f into (8) and note that (6), we have I(f f α ) =, f (x i ) = f α (x i ), i =,,, n. Hence f f α is a linear polynomial. Since the points x, x,, x n are not collinear, we have f = f α, i.e. the minimizer of functional () is only given by (4) and (5)-(7). Theorem.4. There is a unique solution for the linear system of equations (5),(6) and (7). Proof. Let A = ( x i x j log x i x j ) n n and D = diag(β,, β n ), y = (y,, y n ) T in which β i and y i are given by follows 4παss i σi β i =, i I a, ỹ i, i I a,, y i =, i I e, y i, i I e. Denote the transpose of matrix P to be P T = ξ ξ ξ n η η η n, then the linear system of equations (5), (6) and (7) can be rewritten as (9) A + D P c = y, P T d
8 8 T. WEI, Y. C. HON, AND Y. B. WANG where c = (c, c,, c n ) T and d = (d, d,, d k ) T are vectors of coefficients in (4). Suppose (a T, b T ) = (a, a,, a n, b,, b k ) is any solution for the homogeneous equations of (9), i.e. (A + D)a + P b = and P T a =. This implies a T (A + D)a =. In the following we show that a T Aa. Since a T Da, we have a T Da = and hence a i = for all i I a. Denote g(x) = n i= a j x x j log x x j. From P T a = and Lemma., we have () I(g) = n a i g(x i ) = i= n a i a j x i x j log( x i x j ) = a T Aa. i,j= Combining these results gives () A P a = A + D P P T b P T a b =. When points x, x,, x n are not collinear, the (n + ) (n + ) matrix A P P T is nonsingular (see Powell []). Thus (a T, b T ) = and the matrix A + D P P T is nonsingular. The proof is completed. Note that a special case I a =, that is, formula (4) and (6), (7) give a thin plate splines interpolation on the data {(x i, y i ), i =,,, n}. Then from Theorem., we know the interpolation function minimizes the functional I(f) in space Y. In fact, we have Corollary.5. Suppose that I a = and y X and s is a thin plate spine interpolation function given by (4), (6) and (7). Then we have () I(s y) + I(s) = I(y) Proof. It is easy to complete the proof by using the fact R(s, s y) = that comes from Lemma..
9 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA 9. Error estimation In this section, we give an error estimation on the proposed method. The result will use the following lemma given by Powell s paper []. Lemma.. Let {x, x,, x m } be any finite set of distinct points in R that are not collinear and {c, c,, c m } be any real coefficients that satisfy the conditions () m m c i = c i x i =. i= i= Furthermore, let g be any function from R to R that has square integrable second derivatives. Then the linear functional m (4) L(g) = c i g(x i ) has the property m m (5) L(g) (8π) I(g) c i c j φ( x i x j ) where φ(r) = r log(r). i= i= j= /, Let e = f α y be the error function between the thin plate spline approximation f α and the exact function y. The error estimates on e L (Ω) and e L (Ω) are given in the following theorems. Hereafter e L (Ω) denotes the L (Ω) norm of function e(x). Theorem.. For any y X and f α given in Theorem., we have (6) f α y L (Ω) (6π) / (s log ) / ( α + I(y) ) / h +6s / ( + αi(y)) / δ. Proof. For every x Ω, suppose that x is in a triangle ijk with vertices x i, x j, x k. Then x = λ x i + λ x j + λ x k, λ + λ + λ = and λ, λ, λ. The conditions of Lemma. are satisfied if we let m = 4, c = λ, c = λ, c = λ, c 4 =, x = x i, x = x j, x = x k, x 4 = x and let g be the error function e = f α y. We have (7) 4 L(e) = c i e(x i ) i= = λ e(x i ) + λ e(x j ) + λ e(x k ) e(x),
10 T. WEI, Y. C. HON, AND Y. B. WANG and hence, (8) e(x) L(e) + λ e(x i ) + λ e(x j ) + λ e(x k ) [ (8π) I(e)C(λ) ] / + [ e (x i ) + e (x j ) + e (x k ) ] /, where (see Powell s paper []) (9) 4 4 C(λ) = λ i λ j φ( x i x j ), i= j= (h /) log. Hence, for every x ijk, we have () e (x) (4π) I(e)C(λ) + 6 [ e (x i ) + e (x j ) + e (x k ) ]. Integrating the above inequality on the triangle ijk, we get () ijk e dx (4π) I(e)C(λ)vol( ijk ) + 6 [ e (x i ) + e (x j ) + e (x k ) ] vol( ijk ). Summing over all triangles of Ω then gives n () e dx (4π) I(e)C(λ)s + 6 s i e (x i ). Furthermore, () n i= Ω s i s e (x i ) s i δ s i I a δ ( + Φ(y)), i= (f α (x i ) ỹ i ) + (ỹ i y i ), σ i δ ( + αi(y)). From (4) αi(f α ) Φ(f α ) Φ(y) + αi(y), we then have (5) I(f α ) α + I(y) and (6) I(e) (I(f α ) + I(y)) (α + I(y))
11 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA Substituting equations (), (6) and (9) into () yields (7) e dx (6π) s log ( α + I(y) ) h Ω +6sδ ( + αi(y)). Hence, we have (8) e L (Ω) (6π) / (s log ) / ( α + I(y) ) / h +6s / ( + αi(y)) / δ. Lemma.. Let Ω be a bounded domain in R satisfying the cone property. There exists a constant K = K(Ω) such that for < ε and for any u H (Ω), (9) u L (Ω) Kε(I(u)) / + Kε u L (Ω). Proof. see Adams []. Theorem.4. Suppose that f α is a thin plate spline approximation given by Theorem.. If y X, then for sufficiently small δ and h we have ( (4) f α y L (Ω) Ks/4 D h / + D δ /), where D and D are constants depending on parameter α and I(y). Proof. Denote e(x) = f α (x) y(x). For sufficiently small δ and h, from Theorem. we can get e L (Ω) Theorem., we have (4) e L (Ω) K e / where D =. By putting ε = e / L (Ω) L (Ω) [ ] (I(e)) / + [ K (6π) /4 (s log ) /4 ( α + I(y) ) /4 h / in Lemma. and using + 6 / s /4 ( + αi(y)) /4 δ /] [ ] / α / + (I(y)) / + = Ks /4 D h / + Ks /4 D δ /, [ (6π) /4 (log ) /4 ( α + I(y) ) ] [ ] /4 / α / + (I(y)) / +
12 T. WEI, Y. C. HON, AND Y. B. WANG and D = [ ] [ ] 6 / ( + αi(y)) /4 / α / + I / (y)) +. Remark.5. If the value of α = is chosen as one in the paper [9], then D and D are constants depending only on I(y). This implies that this choice of the parameter α is admissible in our proposed method. Remark.6. We can rewrite the right hand side of (4) as f α y L (Ω) Ks /4 [M α /4 + M + M α /4 + M 4 α / + M 5 α /4 ], where M = 6 / δ(i(y)) /4 + 6 / δ(i(y)) /4, M = /4 π /4 (log()) /4 h(i(y)) /4 +6 / 5/4 δ(i(y)) / + 6 / /4 δ +π /4 /4 (log()) /4 h(i(y)) /4, M = /4 π /4 /4 (log()) /4 h(i(y)) / +π /4 6 /4 (log()) /4 h + / δ(i(y)) /4, M 4 = /4 6 / δ + π /4 /4 / (log()) /4 h(i(y)) /4, M 5 = π /4 /4 / (log()) /4 h. By choosing a suitable parameter α > such that [M α /4 + M + M α /4 + M 4 α /4 + M 5 α + M 6 α 5/4 ] = min, the corresponding solution f α provides a more accurate approximation to the solution. In numerical computations, we can replace I(y) by I( s) in which s is the interpolation function of form (4) under the requirement of (7) on the noisy data {(x i, ỹ i )} n i=. If the domain Ω is a square and all measurement data are given on points {x j } n j= which generate a square partition on Ω, then we can modify the functional () to be (4) Φ(f) = ( ) f(xi ) ỹ i + αi(f). n i I a σ i
13 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA The minimizer f α of functional (4) is determined by (4), (6), (7) and following substitution of (5): (4) f α (x i ) + 8πnασi c i = ỹ i, i I a. Noting that the formula (9) under this special case will be changed to the following form as given in []: 5 5 (44) C(λ) = λ i λ j φ( x i x j ), i= j= 4 (log )h, where h is the side length of each small square. By the similar proof, we can then give an error estimation for the function and its derivatives in the following theorem. Theorem.7. If Ω is a square [a, b] [a, b] and all points {x j } n j= generate a square mesh with size h h, h = (b a)/( n ), then we have the following estimations (45) f α y L (Ω) (8π) / ( log ) / (b a) ( α + I(y) ) / h +6(b a) ( + αi(y)) / δ, and (46) f α y) L (Ω) K(b a) / [ / α / + (I(y)) / + ] [ (8π) /4 ( log ) /4 ( α + I(y) ) /4 h / + 4 ( + αi(y)) /4 δ /], where f α is given by (4), (6), (7) and (4). Remark.8. By choosing a suitable parameter α > such that the right hand side of (46) except constant K attains its minimum, the corresponding solution f α is a more accurate approximation. 4. Numerical results To verify the effectiveness of the proposed algorithm, we given in this section several tests under various cases. We take a smooth function y = sin(πξ)sin(πη)e (ξ +η ) for illustration.
14 4 T. WEI, Y. C. HON, AND Y. B. WANG To estimate the computational error of the proposed approximation, we choose N test points {t i } N i= over domain Ω and then compute the Root Mean Square Errors by the following formula (47) E(f α ) = N and (48) E( f α ) = N where is the -norm in R. N (f α (t i ) y(t i )), i= N f α (t i ) y(t i ), i= Case : Ω is a square domain [, ] [, ] and the points are uniformly distributed on Ω and generate a partition of h h square grids. The noisy data are generated by adding random numbers uniformly distributed over [ δ, δ]. In our computation we take h =. and δ =.5. The numerical results are shown in Figures -. Figure (a) displays the surface of exact function y and Figure (b) is the plot of its approximation f α while α =.54 given by Remark.8. In this case E(f α ) =.5, E( f α ) =.8. In Figure, we display five curves of function y(ξ, η and f α (ξ, η) for η =,,,, respectively. The solid line represents the exact function and the dotted line represents its numerical approximation. We can observe from the figures that the approximation f α matches the exact solution y well which indicates that the proposed priori choice of the regularization parameter works fine. Case : we partition the square domain [, ] [, ] into small square grids like Case and further divide each square grid into two triangles by using their diagonal. The length of the longest triangle side is h =.. The noisy data are generated by adding random numbers uniformly distributed over [ δ, δ]. In this case we take δ =.5. In Figure, we plot the curves for y(ξ, η) and f α (ξ, η) while η =,,,, respectively. In Figure (a), we choose α = and the root mean square errors are E(f ) =.86, E( f ) =.45. In Figure (b), we use the parameter α =.46 given in Remark.6 and the root mean square errors are then E(f α ) =., E( f α ) =.75. The numerical results in this case show that even for large noise, the proposed regularized parameter α gives a more accurate approximation than the constant parameter α =.
15 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA 5 Parameter α= y.5 f α.5.5 η ξ η ξ (a) (b) Figure. The exact gradient function and its approximation. (a) y (b) f α for α =.54. α=.54, η=,,,,.5 y(ξ,η) and f α (ξ,η) ξ Figure. The curves of functions y(ξ, η) and f α (ξ, η) for η =,,,, respectively. Case : we choose 665 scattered points on [, ] [, ] as shown in Figure 4. The length of longest triangle side is given by h =.58. The noisy data are generated by adding random numbers uniformly distributed over [.5,.5]. The computational results are displayed in Figure 5 in which curves have the same sense with Figure 4. In Figure 5(a), we take α = and RMSEs are E(f ) =.8, E( f ) =.96. In Figure 5(b), α =.4 and E(f α ) =.6, E( f α ) =.. All the numerical results show that for small noise f α match y well for two different parameters α. But the choice strategy given in Remark.6 works more accurate. Our proposed method is much more simple compared with other methods that use the discrepancy principle in [9].
16 6 T. WEI, Y. C. HON, AND Y. B. WANG α=, η=,,,, α=.466, η=,,,,.5.5 y(ξ,η) and f α (ξ,η).5 y(ξ,η) and f α (ξ,η) ξ ξ (a) (b) Figure. The curves of function y(ξ, η) and its approximation f α (ξ, η) for η =,,,,. (a) α = (b) α =.46. Figure 4. A triangular subdivision of the scattered points. Modifying the number n of specified points in Case and other parameters are fixed as δ =., ỹ i = y i +δsin(/πξ i )sin(/πη i ), i =,,, n. The root mean square errors for the function y and its gradient y with respect to n are shown in Table. The upper value is E(f α ) and the lower one is E( f α ) in every sheet of tables. Numerical results show that even for dense specified points, the proposed algorithm can given an accurate approximation when the finite difference method doesn t work. Table gives the root mean square errors for the function y and its gradient y under different noise levels δ for Case. In the computation, we take n =, ỹ i = y i + δsin(/πξ i )sin(/πη i ), i =,,, n. Numerical results shown in
17 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA 7 α=, η=,,,, α=.44, η=,,,,.5.5 y(ξ,η) and f α (ξ,η).5 y(ξ,η) and f α (ξ,η) ξ ξ (a) (b) Figure 5. The plots of the function y and its approximation f α for η =,,,,. (a) α = (b) α =.4. n = n = n = α = α = α α Table. Root mean square errors with respect to the numbers of total specified points. δ =. δ =. δ =.5 δ =. δ =.5 δ =. α =.69e α = α α Table. Root mean square errors with respect to noise levels δ. Table indicate that the proposed algorithm is stable for the choice of regularization parameters and effective for the data with these noise levels.
18 8 T. WEI, Y. C. HON, AND Y. B. WANG 5. Conclusion In this paper, we provide a thin plate spline approximation to a classical illposed problem numerical derivatives of two dimensional scattered data. A simple algorithm for the choice of a suitable regularization parameter is proposed. Proof on the convergence rate is given for the case when the second order derivative of the exact solution is square integrable. Numerical results for various cases show that the proposed method is effective and stable. References [] R. A. Adams, Sobolev spaces, Pure and Applied Mathematics, Vol. 65, Academic Press, New York-London, 975. [] R. S. Anderssen and M. Hegland, For numerical differentiation, dimensionality can be a blessing!, Mathematics of Computation, 68(999), -4. [] J. Cullum, Numerical differentiation and regularization, SIAM Journal on Numerical Analysis, 8(97), no., [4] J. Cheng, Y.C. Hon and Y.B. Wang, A numerical method for the discontinuous solutions of Abel intergral equations.(preprint) [5] S. R. Deans, Radon transform and its applications, A Wiley-Interscience Publication, 98. [6] R. Gorenflo and S. Vessella, Abel integral equations, Analysis and applications, Lecture Notes in Mathematics, 46, Springer-Verlag, Berlin, 99. [7] C. W. Groetsch, Diferentiation of approximately specified functions, The American Mathematical Monthly, 98(99), [8] M. Hanke and O. Scherzer, Error analysis of an equation error method for the identification of the diffusion coefficient in a quasi-linear parabolic differential equation, SIAM Journal on Applied Mathematics, 59(999), no., -7. [9] M. Hanke and O. Scherzer, Inverse problems light: numerical differentiation, The American Mathematical Monthly, 8(), no. 6, 5-5. [] F. J. Hickernell and Y.C. Hon, Radial basis function approximation of the surface wind filed from scattered data, International Journal of Applied Science and Computations, 4(998), no., -47. [] F. J. Hickernell and Y. C. Hon, Radial basis function approximations as smoothing splines, Applied Mathematics and Computation, (999), -4. [] W. Light, Advances in numerical analysis, Vol.II, Wavelets, Subdivision Algorithms, and Radial Basis Fuctions, Oxford Science Publication, 99. [] M. J. D. Powell, The uniform convergence of thin plate spline interpolation in two dimensions, Numerische mathematik, 68(994), 7-8.
19 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA 9 [4] C. K. Pallaghy and U. Luttge, Light-induced and H ion fluxes and bioelectric phenomena in mesophyll cells of atriplex spongiosa, Zeit. fuer Pflanz., 6(97), [5] R. Qu, A new approach to numerical differentiation and integration, Mathematical and Computer Modelling, 4(996), no., [6] T. J. Rivlin, Optimally stable Lagrangian numerical differentiation, SIAM Journal on Numerical Analysis, (975), no. 5, [7] A. G. Ramm and A. B. Smirnova, On stable numerical differentiation, Mathematics of Computation, 7 (), -5. [8] A. N. Tikhonov and V. Y. Arsenin, Solutions of ill-posed problems, Winston and Sons, Washington 977. [9] Y. B. Wang, X. Z. Jia and J. Cheng, A numerical differentiation method and its application to reconstruction of discontinuity, Inverse Problem, 8(), [] G. Wahba and J. Wendelberger, Some new mathematical methods for variational objective analysis using splines and cross validation, Monthly Weather Review, 8(98), -4. [] G.Wahba, Spline Models for Observational Data, SIAM, Philadelphia, 99.
20 T. WEI, Y. C. HON, AND Y. B. WANG.5 y.5.5 η ξ Figure (a)
21 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA Parameter α= f α.5 η ξ Figure (b)
22 T. WEI, Y. C. HON, AND Y. B. WANG α=.54, η=,,,,.5 y(ξ,η) and f α (ξ,η) ξ Figure
23 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA α=, η=,,,,.5 y(ξ,η) and f α (ξ,η) ξ Figure (a)
24 4 T. WEI, Y. C. HON, AND Y. B. WANG α=.466, η=,,,,.5 y(ξ,η) and f α (ξ,η) ξ Figure (b)
25 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA 5 Figure 4
26 6 T. WEI, Y. C. HON, AND Y. B. WANG α=, η=,,,,.5 y(ξ,η) and f α (ξ,η) ξ Figure 5(a)
27 NUMERICAL DERIVATIVES FOR SCATTERED NOISY DATA 7 α=.44, η=,,,,.5 y(ξ,η) and f α (ξ,η) ξ Figure 5(b)
28 8 T. WEI, Y. C. HON, AND Y. B. WANG Department of Mathematics, City University of Hong Kong and department of Mathematics, Lanzhou University, China address: Department of Mathematics, City University, Hong Kong, China address: Department of Mathematics, Fudan University, China address:
Two-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
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 informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationEstimation of transmission eigenvalues and the index of refraction from Cauchy data
Estimation of transmission eigenvalues and the index of refraction from Cauchy data Jiguang Sun Abstract Recently the transmission eigenvalue problem has come to play an important role and received a lot
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 informationAn l 1 Regularized Method for Numerical Differentiation Using Empirical Eigenfunctions
Journal of Matheatical Research with Applications Jul., 207, Vol. 37, No. 4, pp. 496 504 DOI:0.3770/j.issn:2095-265.207.04.0 Http://jre.dlut.edu.cn An l Regularized Method for Nuerical Differentiation
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationMath 660-Lecture 15: Finite element spaces (I)
Math 660-Lecture 15: Finite element spaces (I) (Chapter 3, 4.2, 4.3) Before we introduce the concrete spaces, let s first of all introduce the following important lemma. Theorem 1. Let V h consists of
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 informationREGULARIZATION PARAMETER SELECTION IN DISCRETE ILL POSED PROBLEMS THE USE OF THE U CURVE
Int. J. Appl. Math. Comput. Sci., 007, Vol. 17, No., 157 164 DOI: 10.478/v10006-007-0014-3 REGULARIZATION PARAMETER SELECTION IN DISCRETE ILL POSED PROBLEMS THE USE OF THE U CURVE DOROTA KRAWCZYK-STAŃDO,
More informationENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT
ENERGY METHODS IN IMAGE PROCESSING WITH EDGE ENHANCEMENT PRASHANT ATHAVALE Abstract. Digital images are can be realized as L 2 (R 2 objects. Noise is introduced in a digital image due to various reasons.
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationA CHARACTERIZATION OF STRICT LOCAL MINIMIZERS OF ORDER ONE FOR STATIC MINMAX PROBLEMS IN THE PARAMETRIC CONSTRAINT CASE
Journal of Applied Analysis Vol. 6, No. 1 (2000), pp. 139 148 A CHARACTERIZATION OF STRICT LOCAL MINIMIZERS OF ORDER ONE FOR STATIC MINMAX PROBLEMS IN THE PARAMETRIC CONSTRAINT CASE A. W. A. TAHA Received
More informationRKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets Class 22, 2004 Tomaso Poggio and Sayan Mukherjee
RKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets 9.520 Class 22, 2004 Tomaso Poggio and Sayan Mukherjee About this class Goal To introduce an alternate perspective of RKHS via integral operators
More information2 FRED J. HICKERNELL the sample mean of the y (i) : (2) ^ 1 N The mean square error of this estimate may be written as a sum of two parts, a bias term
GOODNESS OF FIT STATISTICS, DISCREPANCIES AND ROBUST DESIGNS FRED J. HICKERNELL Abstract. The Cramer{Von Mises goodness-of-t statistic, also known as the L 2 -star discrepancy, is the optimality criterion
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
More information(d). Why does this imply that there is no bounded extension operator E : W 1,1 (U) W 1,1 (R n )? Proof. 2 k 1. a k 1 a k
Exercise For k 0,,... let k be the rectangle in the plane (2 k, 0 + ((0, (0, and for k, 2,... let [3, 2 k ] (0, ε k. Thus is a passage connecting the room k to the room k. Let ( k0 k (. We assume ε k
More informationAN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION
J. KSIAM Vol.19, No.4, 409 416, 2015 http://dx.doi.org/10.12941/jksiam.2015.19.409 AN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION MORAN KIM 1 AND CHOHONG MIN
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
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 informationLecture 2: Tikhonov-Regularization
Lecture 2: Tikhonov-Regularization Bastian von Harrach harrach@math.uni-stuttgart.de Chair of Optimization and Inverse Problems, University of Stuttgart, Germany Advanced Instructional School on Theoretical
More informationANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS INVOLVING A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL
Electronic Journal of Differential Equations, Vol. 217 (217), No. 289, pp. 1 6. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationMorozov s discrepancy principle for Tikhonov-type functionals with non-linear operators
Morozov s discrepancy principle for Tikhonov-type functionals with non-linear operators Stephan W Anzengruber 1 and Ronny Ramlau 1,2 1 Johann Radon Institute for Computational and Applied Mathematics,
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 informationOn the quasi- interpolation
On the quasi- interpolation WU Zongmin Shanghai Key Laboratory for Contemporary Applied Mathematics School of Mathematical Sciences, Fudan University zmwu@fudan.edu.cn http://homepage.fudan.edu.cn/ ZongminWU/
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More informationConvergence rates in l 1 -regularization when the basis is not smooth enough
Convergence rates in l 1 -regularization when the basis is not smooth enough Jens Flemming, Markus Hegland November 29, 2013 Abstract Sparsity promoting regularization is an important technique for signal
More informationA Family of Preconditioned Iteratively Regularized Methods For Nonlinear Minimization
A Family of Preconditioned Iteratively Regularized Methods For Nonlinear Minimization Alexandra Smirnova Rosemary A Renaut March 27, 2008 Abstract The preconditioned iteratively regularized Gauss-Newton
More informationThe impact of a curious type of smoothness conditions on convergence rates in l 1 -regularization
The impact of a curious type of smoothness conditions on convergence rates in l 1 -regularization Radu Ioan Boț and Bernd Hofmann March 1, 2013 Abstract Tikhonov-type regularization of linear and nonlinear
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationRadial Basis Functions I
Radial Basis Functions I Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo November 14, 2008 Today Reformulation of natural cubic spline interpolation Scattered
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
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 informationCURVATURE ESTIMATES FOR WEINGARTEN HYPERSURFACES IN RIEMANNIAN MANIFOLDS
CURVATURE ESTIMATES FOR WEINGARTEN HYPERSURFACES IN RIEMANNIAN MANIFOLDS CLAUS GERHARDT Abstract. We prove curvature estimates for general curvature functions. As an application we show the existence of
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationLecture 10 Polynomial interpolation
Lecture 10 Polynomial interpolation Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationApplied Numerical Analysis Quiz #2
Applied Numerical Analysis Quiz #2 Modules 3 and 4 Name: Student number: DO NOT OPEN UNTIL ASKED Instructions: Make sure you have a machine-readable answer form. Write your name and student number in the
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationReproducing Kernel Hilbert Spaces Class 03, 15 February 2006 Andrea Caponnetto
Reproducing Kernel Hilbert Spaces 9.520 Class 03, 15 February 2006 Andrea Caponnetto About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 12, 2007 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationDetermination of a source term and boundary heat flux in an inverse heat equation
ISSN 1746-7659, England, UK Journal of Information and Computing Science Vol. 8, No., 013, pp. 3-114 Determination of a source term and boundary heat flux in an inverse heat equation A.M. Shahrezaee 1
More informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More informationReproducing Kernel Hilbert Spaces
9.520: Statistical Learning Theory and Applications February 10th, 2010 Reproducing Kernel Hilbert Spaces Lecturer: Lorenzo Rosasco Scribe: Greg Durrett 1 Introduction In the previous two lectures, we
More informationComplexity and regularization issues in kernel-based learning
Complexity and regularization issues in kernel-based learning Marcello Sanguineti Department of Communications, Computer, and System Sciences (DIST) University of Genoa - Via Opera Pia 13, 16145 Genova,
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationAnalysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient
Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Yifeng Xu 1 Jun Zou 2 Abstract Based on a new a posteriori error estimator, an adaptive finite element method is proposed
More informationLeast Squares Approximation
Chapter 6 Least Squares Approximation As we saw in Chapter 5 we can interpret radial basis function interpolation as a constrained optimization problem. We now take this point of view again, but start
More informationAn eigenvalue method using multiple frequency data for inverse scattering problems
An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationStochastic homogenization 1
Stochastic homogenization 1 Tuomo Kuusi University of Oulu August 13-17, 2018 Jyväskylä Summer School 1 Course material: S. Armstrong & T. Kuusi & J.-C. Mourrat : Quantitative stochastic homogenization
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 informationSzemerédi s Lemma for the Analyst
Szemerédi s Lemma for the Analyst László Lovász and Balázs Szegedy Microsoft Research April 25 Microsoft Research Technical Report # MSR-TR-25-9 Abstract Szemerédi s Regularity Lemma is a fundamental tool
More informationA Posteriori Estimates for Cost Functionals of Optimal Control Problems
A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany
More informationFinding discontinuities of piecewise-smooth functions
Finding discontinuities of piecewise-smooth functions A.G. Ramm Mathematics Department, Kansas State University, Manhattan, KS 66506-2602, USA ramm@math.ksu.edu Abstract Formulas for stable differentiation
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationKernels A Machine Learning Overview
Kernels A Machine Learning Overview S.V.N. Vishy Vishwanathan vishy@axiom.anu.edu.au National ICT of Australia and Australian National University Thanks to Alex Smola, Stéphane Canu, Mike Jordan and Peter
More informationAn interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes
An interpolation operator for H 1 functions on general quadrilateral and hexahedral meshes with hanging nodes Vincent Heuveline Friedhelm Schieweck Abstract We propose a Scott-Zhang type interpolation
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationDefinition and basic properties of heat kernels I, An introduction
Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationb i (µ, x, s) ei ϕ(x) µ s (dx) ds (2) i=1
NONLINEAR EVOLTION EQATIONS FOR MEASRES ON INFINITE DIMENSIONAL SPACES V.I. Bogachev 1, G. Da Prato 2, M. Röckner 3, S.V. Shaposhnikov 1 The goal of this work is to prove the existence of a solution to
More informationS. Lototsky and B.L. Rozovskii Center for Applied Mathematical Sciences University of Southern California, Los Angeles, CA
RECURSIVE MULTIPLE WIENER INTEGRAL EXPANSION FOR NONLINEAR FILTERING OF DIFFUSION PROCESSES Published in: J. A. Goldstein, N. E. Gretsky, and J. J. Uhl (editors), Stochastic Processes and Functional Analysis,
More informationOn the number of ways of writing t as a product of factorials
On the number of ways of writing t as a product of factorials Daniel M. Kane December 3, 005 Abstract Let N 0 denote the set of non-negative integers. In this paper we prove that lim sup n, m N 0 : n!m!
More informationA uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance tomography
c de Gruyter 28 J. Inv. Ill-Posed Problems 16 28), 381 396 DOI 1.1515 / JIIP.28.21 A uniqueness result and image reconstruction of the orthotropic conductivity in magnetic resonance electrical impedance
More informationSpectra of Multiplication Operators as a Numerical Tool. B. Vioreanu and V. Rokhlin Technical Report YALEU/DCS/TR-1443 March 3, 2011
We introduce a numerical procedure for the construction of interpolation and quadrature formulae on bounded convex regions in the plane. The construction is based on the behavior of spectra of certain
More informationarxiv: v1 [math.na] 1 May 2013
arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based
More informationA Power Penalty Method for a Bounded Nonlinear Complementarity Problem
A Power Penalty Method for a Bounded Nonlinear Complementarity Problem Song Wang 1 and Xiaoqi Yang 2 1 Department of Mathematics & Statistics Curtin University, GPO Box U1987, Perth WA 6845, Australia
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationInitial Temperature Reconstruction for a Nonlinear Heat Equation: Application to Radiative Heat Transfer.
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 Initial Temperature Reconstruction for a Nonlinear Heat Equation:
More informationADAPTIVE MULTISCALE MOMENT METHOD FOR SOLVING TWO-DIMENSIONAL FREDHOLM INTEGRAL EQUATION OF THE FIRST KIND
Progress In Electromagnetics Research, PIER 21, 173 201, 1999 ADAPTIVE MULTISCALE MOMENT METHOD FOR SOLVING TWO-DIMENSIONAL FREDHOLM INTEGRAL EQUATION OF THE FIRST KIND C. Su Dept. of Applied Math. Northwestern
More informationD. Shepard, Shepard functions, late 1960s (application, surface modelling)
Chapter 1 Introduction 1.1 History and Outline Originally, the motivation for the basic meshfree approximation methods (radial basis functions and moving least squares methods) came from applications in
More informationReproducing Kernel Hilbert Spaces
Reproducing Kernel Hilbert Spaces Lorenzo Rosasco 9.520 Class 03 February 9, 2011 About this class Goal To introduce a particularly useful family of hypothesis spaces called Reproducing Kernel Hilbert
More informationNedelec elements for computational electromagnetics
Nedelec elements for computational electromagnetics Per Jacobsson, June 5, 2007 Maxwell s equations E = jωµh () H = jωεe + J (2) (εe) = ρ (3) (µh) = 0 (4) J = jωρ (5) Only three of the equations are needed,
More informationLecture 2: Linear Algebra Review
EE 227A: Convex Optimization and Applications January 19 Lecture 2: Linear Algebra Review Lecturer: Mert Pilanci Reading assignment: Appendix C of BV. Sections 2-6 of the web textbook 1 2.1 Vectors 2.1.1
More informationA rational interpolation scheme with super-polynomial rate of convergence
Center for Turbulence Research Annual Research Briefs 2008 31 A rational interpolation scheme with super-polynomial rate of convergence By Q. Wang, P. Moin AND G. Iaccarino 1. Motivation and objectives
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationDue Giorni di Algebra Lineare Numerica (2GALN) Febbraio 2016, Como. Iterative regularization in variable exponent Lebesgue spaces
Due Giorni di Algebra Lineare Numerica (2GALN) 16 17 Febbraio 2016, Como Iterative regularization in variable exponent Lebesgue spaces Claudio Estatico 1 Joint work with: Brigida Bonino 1, Fabio Di Benedetto
More informationElements of Positive Definite Kernel and Reproducing Kernel Hilbert Space
Elements of Positive Definite Kernel and Reproducing Kernel Hilbert Space Statistical Inference with Reproducing Kernel Hilbert Space Kenji Fukumizu Institute of Statistical Mathematics, ROIS Department
More informationMultivariate Interpolation with Increasingly Flat Radial Basis Functions of Finite Smoothness
Multivariate Interpolation with Increasingly Flat Radial Basis Functions of Finite Smoothness Guohui Song John Riddle Gregory E. Fasshauer Fred J. Hickernell Abstract In this paper, we consider multivariate
More informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationFitting multidimensional data using gradient penalties and combination techniques
1 Fitting multidimensional data using gradient penalties and combination techniques Jochen Garcke and Markus Hegland Mathematical Sciences Institute Australian National University Canberra ACT 0200 jochen.garcke,markus.hegland@anu.edu.au
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationTikhonov Regularization for Weighted Total Least Squares Problems
Tikhonov Regularization for Weighted Total Least Squares Problems Yimin Wei Naimin Zhang Michael K. Ng Wei Xu Abstract In this paper, we study and analyze the regularized weighted total least squares (RWTLS)
More informationTHEOREMS, ETC., FOR MATH 515
THEOREMS, ETC., FOR MATH 515 Proposition 1 (=comment on page 17). If A is an algebra, then any finite union or finite intersection of sets in A is also in A. Proposition 2 (=Proposition 1.1). For every
More informationDIEUDONNE AGBOR AND JAN BOMAN
ON THE MODULUS OF CONTINUITY OF MAPPINGS BETWEEN EUCLIDEAN SPACES DIEUDONNE AGBOR AND JAN BOMAN Abstract Let f be a function from R p to R q and let Λ be a finite set of pairs (θ, η) R p R q. Assume that
More informationJacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation
Wei et al. SpringerPlus (06) 5:70 DOI 0.86/s40064-06-3358-z RESEARCH Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation Open Access
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
More information