This function inverse a Vandermonde Matrix B.
Matrix B is a n-by-n matrix, its (i,j) entry is i^(j-1),
where i,j = 1,2,...,n
for example, n = 4
1 1 1 1
1 2 4 8
1 3 9 27
1 4 16 64
This routine uses a Stirling polynomial(the first kind) coefficients
For fast operation, a C Stirling coefficient function has
been posted with name: mStirling.c. The C-version of this
inverse function is also available upon request.
Steven Huang (2021). A Vandermonde matrix inversion (https://www.mathworks.com/matlabcentral/fileexchange/8048-a-vandermonde-matrix-inversion), MATLAB Central File Exchange. Retrieved .
The function does not invert the Vandermonde Matrix B that is given as an input, but returns the inverse a particular Vandermonde matrix of the same size, to wit the one of which the second column is (1:n).' .
Seems to work well for low dimension case. But gets worse quickly as n grows.
You can use the code to invert an Vandermonde matrix with n=20, and multiply it back and compare the result with the identity matrix. The error is huge.
First sentence should read: "This function inverts..." instead of "This function inverse..."
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!