What is the method of interpolation in the function pdeinterpolant?

6 visualizaciones (últimos 30 días)
Alex Degeling
Alex Degeling el 19 de Sept. de 2024
Respondida: Abhas el 30 de Sept. de 2024
What is the numerical method being used for interpolation in the function "pdeinterpolant"? I couldn't find any info on the actual algorithm being used (e.g. using the command doc pdeinterpolant). In particular I'd like to know whether I can expect C0 or C1 continuity in the interpolant when used on 2D quadratic triangulations.

Respuestas (1)

Abhas
Abhas el 30 de Sept. de 2024
Hi Alex,
The "pdeInterpolant" function in MATLAB is used to interpolate solutions of partial differential equations (PDEs) at arbitrary points within the geometry. However, the specific numerical method used for interpolation in pdeInterpolant is not explicitly documented in the MATLAB documentation.
From general principles of finite element methods (FEM) and interpolation, we can infer some details:
  • Continuity: For 2D quadratic triangulations, you can typically expect "C0 continuity". This means the interpolant is continuous, but its first derivative may not be continuous across element boundaries. While "C1 continuity" requires both the function and its first derivative to be continuous across element boundaries. Achieving "C1 continuity" typically involves using more complex elements or additional constraints, which is uncommon for standard quadratic elements in finite element methods.
  • Interpolation Method: The interpolation is likely based on shape functions used in "FEM". For quadratic elements, these shape functions is usually done using second-degree polynomials, which provide a higher degree of accuracy compared to linear elements.
You may refer to the below MathWorks documentation links if you need more precise control over the interpolation method or continuity:
  1. https://www.mathworks.com/help/pde/ug/pde.stationaryresults.interpolatesolution.html
  2. https://www.mathworks.com/help/pde/ug/pde.stationaryresults.evaluategradient.html

Categorías

Más información sobre Interpolation en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by