We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Effective and accurate high-degree spline interpolation is still a challenging task in today’s applications. Higher degree spline interpolation is not so commonly used, because it requires the knowledge of higher order derivatives at the nodes of a function on a given mesh.
In this article, our goal is to demonstrate the continuity of the piecewise polynomials and their derivatives at the connecting points, obtained with a method initially developed by Beaudoin (1998, 2003) and Beauchemin (2003). This new method, involving the discrete Fourier transform (DFT/FFT), leads to higher degree spline interpolation for equally spaced data on an interval $[0,T]$. To do this, we analyze the singularities that may occur when solving the system of equations that enables the construction of splines of any degree. We also note an important difference between the odd-degree splines and even-degree splines. These results prove that Beaudoin and Beauchemin’s method leads to spline interpolation of any degree and that this new method could eventually be used to improve the accuracy of spline interpolation in traditional problems.
The Chebyshev conjecture posits that Chebyshev subsets of a real Hilbert space $X$ are convex. Works by Asplund, Ficken and Klee have uncovered an equivalent formulation of the Chebyshev conjecture in terms of uniquely remotal subsets of $X$. In this tradition, we develop another equivalent formulation in terms of Chebyshev subsets of the unit sphere of $X\times \mathbb{R}$. We characterise such sets in terms of the image under stereographic projection. Such sets have superior structure to Chebyshev sets and uniquely remotal sets.
One of the approaches to the Riemann Hypothesis is the Nyman–Beurling criterion. Cotangent sums play a significant role in this criterion. Here we investigate the values of these cotangent sums for various shifts of the argument.
We investigate the Gibbs–Wilbraham phenomenon for generalized sampling series, and related interpolation series arising from cardinal functions. We prove the existence of the overshoot characteristic of the phenomenon for certain cardinal functions, and characterize the existence of an overshoot for sampling series.
We study the asymptotic behaviour of a class of small-noise diffusions driven by fractional Brownian motion, with random starting points. Different scalings allow for different asymptotic properties of the process (small-time and tail behaviours in particular). In order to do so, we extend some results on sample path large deviations for such diffusions. As an application, we show how these results characterise the small-time and tail estimates of the implied volatility for rough volatility models, recently proposed in mathematical finance.
Transform inversions, in which density and survival functions are computed from their associated moment generating function $\mathcal{M}$, have largely been based on methods which use values of $\mathcal{M}$ in its convergence region. Prominent among such methods are saddlepoint approximations and Fourier-series inversion methods, including the fast Fourier transform. In this paper we propose inversion methods which make use of values for $\mathcal{M}$ which lie outside of its convergence region and in its analytic continuation. We focus on the simplest and perhaps richest setting for applications in which $\mathcal{M}$ is either a meromorphic function in its analytic continuation, so that all of its singularities are poles, or else the singularities are isolated essential. Asymptotic expansions of finite- and infinite-orders are developed for density and survival functions using the poles of $\mathcal{M}$ in its analytic continuation. For finite-order expansions, the expansion error is a contour integral in the analytic continuation, which we approximate using the saddlepoint method based on following the path of steepest descent. Such saddlepoint error approximations accurately determine expansion errors and, thus, provide the means for determining the order of the expansion needed to achieve some preset accuracy. They also provide an additive correction term which increases accuracy of the expansion. Further accuracy is achieved by computing the expansion errors numerically using a contour path which ultimately tracks the steepest descent direction. Important applications include Wilks’ likelihood ratio test in MANOVA, compound distributions, and the Sparre Andersen and Cramér–Lundberg ruin models.
In the present paper, an inverse result of approximation, i.e. a saturation theorem for the sampling Kantorovich operators, is derived in the case of uniform approximation for uniformly continuous and bounded functions on the whole real line. In particular, we prove that the best possible order of approximation that can be achieved by the above sampling series is the order one, otherwise the function being approximated turns out to be a constant. The above result is proved by exploiting a suitable representation formula which relates the sampling Kantorovich series with the well-known generalized sampling operators introduced by Butzer. At the end, some other applications of such representation formulas are presented, together with a discussion concerning the kernels of the above operators for which such an inverse result occurs.
We establish the mapping properties of Fourier-type transforms on rearrangement-invariant quasi-Banach function spaces. In particular, we have the mapping properties of the Laplace transform, the Hankel transforms, the Kontorovich-Lebedev transform and some oscillatory integral operators. We achieve these mapping properties by using an interpolation functor that can explicitly generate a given rearrangement-invariant quasi-Banach function space via Lebesgue spaces.
Bilinear fractal interpolation surfaces were introduced by Ruan and Xu in 2015. In this paper, we present the formula for their box dimension under certain constraint conditions.
Let $\unicode[STIX]{x1D6FA}$ be a member of a certain class of convex ellipsoids of finite/infinite type in $\mathbb{C}^{2}$. In this paper, we prove that every holomorphic function in $L^{p}(\unicode[STIX]{x1D6FA})$ can be approximated by holomorphic functions on $\bar{\unicode[STIX]{x1D6FA}}$ in $L^{p}(\unicode[STIX]{x1D6FA})$-norm, for $1\leq p<\infty$. For the case $p=\infty$, the continuity up to the boundary is additionally required. The proof is based on $L^{p}$ bounds in the additive Cousin problem.
We define fractal interpolation on unbounded domains for a certain class of topological spaces and construct local fractal functions. In addition, we derive some properties of these local fractal functions, consider their tensor products, and give conditions for local fractal functions on unbounded domains to be elements of Bochner–Lebesgue spaces.
A new minimization principle for the Poisson equation using two variables – the solution and the gradient of the solution – is introduced. This principle allows us to use any conforming finite element spaces for both variables, where the finite element spaces do not need to satisfy the so-called inf–sup condition. A numerical example demonstrates the superiority of this approach.
In this paper we consider the algorithm for recovering sparse orthogonal polynomials using stochastic collocation via ℓq minimization. The main results include: 1) By using the norm inequality between ℓq and ℓ2 and the square root lifting inequality, we present several theoretical estimates regarding the recoverability for both sparse and non-sparse signals via ℓq minimization; 2) We then combine this method with the stochastic collocation to identify the coefficients of sparse orthogonal polynomial expansions, stemming from the field of uncertainty quantification. We obtain recoverability results for both sparse polynomial functions and general non-sparse functions. We also present various numerical experiments to show the performance of the ℓq algorithm. We first present some benchmark tests to demonstrate the ability of ℓq minimization to recover exactly sparse signals, and then consider three classical analytical functions to show the advantage of this method over the standard ℓ1 and reweighted ℓ1 minimization. All the numerical results indicate that the ℓq method performs better than standard ℓ1 and reweighted ℓ1 minimization.
In this paper we prove a short time asymptotic expansion of a hypoelliptic heat kernel on a Euclidean space and a compact manifold. We study the ‘cut locus’ case, namely, the case where energy-minimizing paths which join the two points under consideration form not a finite set, but a compact manifold. Under mild assumptions we obtain an asymptotic expansion of the heat kernel up to any order. Our approach is probabilistic and the heat kernel is regarded as the density of the law of a hypoelliptic diffusion process, which is realized as a unique solution of the corresponding stochastic differential equation. Our main tools are S. Watanabe’s distributional Malliavin calculus and T. Lyons’ rough path theory.
We address the construction and approximation for feed-forward neural networks (FNNs) with zonal functions on the unit sphere. The filtered de la Vallée-Poussin operator and the spherical quadrature formula are used to construct the spherical FNNs. In particular, the upper and lower bounds of approximation errors by the FNNs are estimated, where the best polynomial approximation of a spherical function is used as a measure of approximation error.
In this paper, a novel second-order two-scale (SOTS) computational method is developed for nonlinear dynamic thermo-mechanical problems of composites with cylindrical periodicity. The non-linearities of these multi-scale problems were caused by the temperature-dependent properties of the composites. Firstly, the formal SOTS solutions for these problems are constructed by the multiscale asymptotic analysis. Then we theoretically explain the importance of the SOTS solutions by the error analysis in the pointwise sense. In addition, a SOTS numerical algorithm is proposed in detail to effectively solve these problems. Finally, some numerical examples verify the feasibility and effectiveness of the SOTS numerical algorithm we proposed.
For a prescribed set of lacunary data with equally spaced knot sequence in the unit interval, we show the existence of a family of fractal splines satisfying for v = 0, 1, … ,N and suitable boundary conditions. To this end, the unique quintic spline introduced by A. Meir and A. Sharma [SIAM J. Numer. Anal. 10(3) 1973, pp. 433-442] is generalized by using fractal functions with variable scaling parameters. The presence of scaling parameters that add extra “degrees of freedom”, self-referentiality of the interpolant, and “fractality” of the third derivative of the interpolant are additional features in the fractal version, which may be advantageous in applications. If the lacunary data is generated from a function Φ satisfying certain smoothness condition, then for suitable choices of scaling factors, the corresponding fractal spline satisfies , as the number of partition points increases.
We transpose the parametric geometry of numbers, recently created by Schmidt and Summerer, to fields of rational functions in one variable and analyze, in that context, the problem of simultaneous approximation to exponential functions.
Let $X$ be a vector space and let $\unicode[STIX]{x1D711}:X\rightarrow \mathbb{R}\cup \{-\infty ,+\infty \}$ be an extended real-valued function. For every function $f:X\rightarrow \mathbb{R}\cup \{-\infty ,+\infty \}$, let us define the $\unicode[STIX]{x1D711}$-envelope of $f$ by
where $\begin{array}{@{}c@{}}-\\ \\ \\ \cdot \end{array}$ denotes the lower subtraction in $\mathbb{R}\cup \{-\infty ,+\infty \}$. The main purpose of this paper is to study in great detail the properties of the important generalized conjugation map $f\mapsto f^{\unicode[STIX]{x1D711}}$. When the function $\unicode[STIX]{x1D711}$ is closed and convex, $\unicode[STIX]{x1D711}$-envelopes can be expressed as Legendre–Fenchel conjugates. By particularizing with $\unicode[STIX]{x1D711}=(1/p\unicode[STIX]{x1D706})\Vert \cdot \Vert ^{p}$, for $\unicode[STIX]{x1D706}>0$ and $p\geqslant 1$, this allows us to derive new expressions of the Klee envelopes with index $\unicode[STIX]{x1D706}$ and power $p$. Links between $\unicode[STIX]{x1D711}$-envelopes and Legendre–Fenchel conjugates are also explored when $-\unicode[STIX]{x1D711}$ is closed and convex. The case of Moreau envelopes is examined as a particular case. In addition to the $\unicode[STIX]{x1D711}$-envelopes of functions, a parallel notion of envelope is introduced for subsets of $X$. Given subsets $\unicode[STIX]{x1D6EC}$, $C\subset X$, we define the $\unicode[STIX]{x1D6EC}$-envelope of $C$ as $C^{\unicode[STIX]{x1D6EC}}=\bigcap _{x\in C}(x+\unicode[STIX]{x1D6EC})$. Connections between the transform $C\mapsto C^{\unicode[STIX]{x1D6EC}}$ and the aforestated $\unicode[STIX]{x1D711}$-conjugation are investigated.
We introduce a multiple interval Chebyshev-Gauss-Lobatto spectral collocation method for the initial value problems of the nonlinear ordinary differential equations (ODES). This method is easy to implement and possesses the high order accuracy. In addition, it is very stable and suitable for long time calculations. We also obtain the hp-version bound on the numerical error of the multiple interval collocation method under H1-norm. Numerical experiments confirm the theoretical expectations.