7
$\begingroup$

Let $\{x_1,\dots,x_n\}$ be pairwise distinct complex numbers and $l_1+l_2+\dots+l_n=N$. The $N\times N$ confluent Vandermonde matrix is defined as $$V= \begin{bmatrix} v_{1,0}&v_{2,0}&\dots&v_{n,0}\\\\ v_{1,1}&v_{2,1}&\dots&v_{n,1}\\\\ \vdots\\\\ v_{1,N-1}&v_{2,N-1}&\dots&v_{n,N-1} \end{bmatrix}$$ where $v_{j,k}=\begin{bmatrix}x_j^k,&kx_j^{k-1},&\dots&k(k-1)\times\dots\times (k-l_j+1) x_j^{k-l_j+1}\end{bmatrix}$. Let $\|\cdot\|$ denote the row sum matrix norm. In some applications (e.g. interpolation, signal processing) one would like to estimate the quantity $\|V^{-1}\|$.

Gautschi [1] has shown that for $l_1=\dots=l_n=2$ one has $$ \|V^{-1}\| \leq \max_{1\leq \lambda\leq n} \beta_{\lambda} \prod_{\nu=1,\nu\neq\lambda}^n \biggl(\frac{1+|x_{\lambda}|}{|x_{\nu}-x_{\lambda}|}\biggr)^2$$ where $\beta_{\lambda}=\max\biggl(1+|x_{\lambda}|,1+2(1+|x_{\lambda}|)\sum_{\nu\neq\lambda}{1\over |x_\nu-x_\lambda|}\biggr)$.

I am interested in a somewhat cruder estimates, as follows: if $|x_j|\leq 1$ and $|x_i-x_j|\geq \delta$, then for the above case we have $$\|V^{-1}\| \leq C n 2^N\delta^{-N+1}\qquad (*)$$ for some absolute constant $C$.

Is it true that something like $(*)$ holds for the general configuration $\{l_1,\dots,l_n\}$?

EDIT: using [2], this seems to boil down to the following. Let $$ h_j(x)=\prod_{i \neq j}(x-x_i)^{-l_i}.$$ For $t=0,1,\dots,l_j$ evaluate $h_j^{(t)}(x_j).$

[1] W.Gautschi, "On Inverses of Vandermonde and Confluent Vandermonde matrices II", Numerische Mathematik 5, 425-430, 1963.

[2] R.Schapelle, "The Inverse of the Confluent Vandermonde Matrix", IEEE Trans. on Automatic Control, October 1972, pp.724-725.

$\endgroup$
1
  • $\begingroup$ Also posted on MSE. $\endgroup$
    – Did
    Dec 1, 2012 at 12:09

1 Answer 1

2
$\begingroup$

For $j=1,\dots,n$ and $k=0,1,\dots,l_j-1$ denote by $u_{j,k}$ the row with index $l_1+\dots +l_{j-1}+k$ of the matrix $V^{-1}$. By using a generalization of the Hermite interpolation formula (see [3]), in [2] it is shown that the elements of $u_{j,k}$ are the coefficients of the polynomial $$ {1\over k!} \sum_{t=0}^{l_j-1-k} {1\over t!} h_j^{(t)}(x_j) (x-x_j)^{k+t} \prod_{i\neq j} (x-x_i)^{l_i} $$

Now thanks to the answer to this MSE question, one has $$ |h_j^{(t)}(x_j)|\leq N(N+1)\cdots (N+t-1)\delta^{-N-t}. $$

The sum of absolute values of the coefficients of the polynomials $(x-x_j)^{k+t} \prod_{i\neq j} (x-x_i)^{l_i}$ is at most (see [4, Lemma]) $$ (1+|x_j|)^{k+t} \prod_{i\neq j}(1+|x_i|)^{l_i} \leq 2^{N-(l_j-k-t)}. $$

So now $$ \|u_{j,k}\| \leq {1\over k!}\sum_{t=0}^{l_j-1-k} {1\over t!} {N(N+1)\cdots (N+t-1) \over {\delta^{N+t}}}2^{N-l_j+k+t}\\\\ \leq \biggl({2\over \delta}\biggr)^N {1\over {2^{l_j-k}k!}}\sum_{t=0}^{l_j-1-k} {l_j-1-k \choose t} {N(N+1)\cdots(N+t-1)\over (l_j-k-t)\cdots(l_j-k-2)(l_j-k-1)} \biggl({2\over \delta}\biggr)^t\\\\ \leq \biggl({2\over \delta}\biggr)^N {1\over {2^{l_j-k}k!}} \biggl(1+{2N\over \delta}\biggr)^{l_j-1-k}\\\\ =\biggl({2\over \delta}\biggr)^N {2\over k!} \biggl({1\over 2}+{N\over\delta}\biggr)^{l_j-1-k}. $$

[3] A.Spitzbart, "A Generalization of Hermite's Interpolation Formula", The American Mathematical Monthly, Vol.67 No.1, p.42-46, 1960.

[4] W.Gautschi, "On Inverses of Vandermonde and Confluent Vandermonde matrices", Numerische Mathematik 4, p.117-123, 1962.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.