File Exchange

image thumbnail

PMPack - Parameterized Matrix Package

version 1.0 (686 KB) by

Polynomial spectral methods for solving parameterized matrix equations.

9 Downloads

Updated

View License

Contains a suite for approximating the solution of a parameterized matrix equation using either a residual minimizing spectral Galerkin method or pseudospectral method. Both methods employ a basis of orthogonal polynomials -- multivariate polynomials are constructed at products of univariate polynomials.

Demos are provided, although the demo for solving the elliptic PDE with a Karhunen-Loeve expansion for the log of the coefficients requires the MATLAB PDE Toolbox.

Also contains utilities for working with the orthogonal polynomials and associated Gaussian quadrature rules.

Comments and Ratings (3)

Ok, from the code I understood that the computed polynomials (Legendre, Chebyshev, Jacobi etc.) are in fact orthonormal and not only orthogonal w.r.t. a certain weight.

So were the legendre polynomials constructed using the uniform weight of 0.5 and chebyshev with 1/pi*(1-x^2)^-0.5? How can we extract the resulting polynomial coefficients?

The Legendre polynomials should assume fixed values at the points x = -1 and x = 1. However, when I plot them using the evaluate_expansion function I see that this is not the case. Is this a bug?

João

João (view profile)

MATLAB Release
MATLAB 7.9 (R2009b)
Acknowledgements

Inspired by: Random Field Simulation

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Win prizes and improve your MATLAB skills

Play today