hanekawa

post

[NP] Fourier Analysis with Plancherel-Parseval Identity

NP

How does Fourier transform work? Fourier transform is just reforming the shape of function. In some cases, it acts as a smoother, e.g. exp(-abs(x)). In the sense of wave propagation, Fourier transform is the superposition of waves with the amplitude equal to our plugin function. But what really make it become so useful, especially in the numerical analysis of PDE? A very important theoretical foundation is Plancherel Theorem, a.k.a Parseval Theorem, which shows that the norm of original function is equal to the norm of its Fourier transformation. With that relation, we can apply Fourier transform on stability analysis: instead of analyzing the original function, we focus on its Fourier transformation. As long as that one satisfies the stable condition, the original one would automatically gain the same result. A question here is, why it's easier? The most well-known case is that using the transformation to turn a differentiation operation (original func space) into a multiplication operation (Fourier trans space). If we used that on the hyperbolic equations, or more specifically, on one way wave equation, we then turn a multivariable PDE into a singlevariable ODE, which we can solve directly. Another case is Von Neumann analysis, where we put Fourier tranform into the analysis of stability of finite difference schemes. We then derive that after Fourier transform, the computed FDM hyperbolic solution can be inductively represented by the scaling of initial data. That scaling, depends on the time step, is called amplification factor, which includes all the info we need for the stability analysis. Hence, at the beginning of Von Neumann analysis, we wrote the Fourier transformed solution into inductive formula, tracking back to the initial data with the amplification factor. And then we just work on that factor.