Pseudo-spectral method solution for wave equation PDE: d^2p/dt^2=​(c^2)*[d^2​p/dx^2 + d^2p/dy^2 + d^2p/dz^2]

20 views (last 30 days)
Hussein Muhammed
Hussein Muhammed on 9 May 2022
I'm trying to solving above wave equation via Pseudo-spectral method instead of finite-difference scheme, for the sake of solving very complicated type of wave equation called: wave equation in Riemannian coordinate system ( can seniors please help me y sugesstions and intial MATLAB codes. An initial MATLAB code for solving the wave equation in Cartesian coordinate system can be found here (

Answers (1)

Bjorn Gustavsson
Bjorn Gustavsson on 9 May 2022
This is just a linear wave-equation with constant speed of the waves. The general solution is just:
Harmonic expansion of waves...
So you just have to determine the complex-valued for all that you need to fit your initial and boundary-conditions. After that you'll have the solution as a number of propagating plane-wave-modes that describe the p-variation. From what you've given us so far it is difficult to give more specific advice.
  1 Comment
Hussein Muhammed
Hussein Muhammed on 9 May 2022
the key step in solving the acoustic wave equation in Riemannian coordinates is that the spatial axes of the laplacian operator have to approximated by Fourier Pseudo-spectral methods, all I need is a working code example of that method for any similer PDE then I can start my code.

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!