We propose a method for designing accurate interpolation formulas on the real axis for the purpose of function approximation in weighted Hardy spaces. In particular, we consider the Hardy space of functions that are analytic in a strip region around the real axis, being characterized by a weight function $w$ that determines the decay rate of its elements in the neighborhood of infinity. Such a space is considered as a set of functions that are transformed by variable transformations that realize a certain decay rate at infinity. Popular examples of such transformations are given by the single exponential (SE) and double exponential (DE) transformations for the SE-Sinc and DE-Sinc formulas, which are very accurate owing to the accuracy of sinc interpolation in the weighted Hardy spaces with single and double exponential weights $w$, respectively. However, it is not guaranteed that the sinc formulas are optimal in weighted Hardy spaces, although Sugihara has demonstrated that they are near optimal. An explicit form for an optimal approximation formula has only been given in weighted Hardy spaces with SE weights of a certain type. In general cases, explicit forms for optimal formulas have not been provided so far. We adopt a potential theoretic approach to obtain almost optimal formulas in weighted Hardy spaces in the case of general weight functions $w$. We formulate the problem of designing an optimal formula in each space as an optimization problem written in terms of a Green potential with an external field. By solving the optimization problem numerically, we obtain an almost optimal formula in each space. Furthermore, some numerical results demonstrate the validity of this method. In particular, for the case of a DE weight, the formula designed by our method outperforms the DE-Sinc formula.
Read full abstract