Generate a Vandermonde matrix of the Laguerre polynomial with float array of points in Python


To generate a pseudo Vandermonde matrix of the Laguerre polynomial, use the laguerre.lagvander() in Python Numpy. The method returns the pseudo-Vandermonde matrix. The shape of the returned matrix is x.shape + (deg + 1,), where The last index is the degree of the corresponding Laguerre polynomial. The dtype will be the same as the converted x.

The parameter, x returns an Array of points. The dtype is converted to float64 or complex128 depending on whether any of the elements are complex. If x is scalar it is converted to a 1-D array. The parameter, deg is the degree of the resulting matrix.

Steps

At first, import the required library −

import numpy as np
from numpy.polynomial import laguerre as L

Create an array −

x = np.array([0, 3.5, -1.4, 2.5])

Display the array −

print("Our Array...\n",c)

Check the Dimensions −

print("\nDimensions of our Array...\n",c.ndim)

Get the Datatype −

print("\nDatatype of our Array object...\n",c.dtype)

Get the Shape −

print("\nShape of our Array object...\n",c.shape)

To generate a pseudo Vandermonde matrix of the Laguerre polynomial, use the laguerre.lagvander() in Python −

print("\nResult...\n",L.lagvander(x, 2))

Example

import numpy as np
from numpy.polynomial import laguerre as L

# Create an array
x = np.array([0, 3.5, -1.4, 2.5])

# Display the array
print("Our Array...\n",x)

# Check the Dimensions
print("\nDimensions of our Array...\n",x.ndim)

# Get the Datatype
print("\nDatatype of our Array object...\n",x.dtype)

# Get the Shape
print("\nShape of our Array object...\n",x.shape)

# To generate a pseudo Vandermonde matrix of the Laguerre polynomial, use the laguerre.lagvander() in Python Numpy
print("\nResult...\n",L.lagvander(x, 2))

Output

Our Array...
   [ 0. 3.5 -1.4 2.5]

Dimensions of our Array...
1

Datatype of our Array object...
float64

Shape of our Array object...
(4,)

Result...
   [[ 1.  1.    1. ]
   [ 1.  -2.5   0.125]
   [ 1.   2.4   4.78 ]
   [ 1.  -1.5  -0.875]]

Updated on: 04-Mar-2022

79 Views

Kickstart Your Career

Get certified by completing the course

Get Started
Advertisements