Mercurial > repos > public > sbplib
comparison four.m @ 820:501750fbbfdb
Merge with feature/grids
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Fri, 07 Sep 2018 14:40:58 +0200 |
| parents | f6ada1719420 |
| children |
comparison
equal
deleted
inserted
replaced
| 819:fdf0ef9150f4 | 820:501750fbbfdb |
|---|---|
| 1 % four returns the fourier transform u_hat of the function u and the frequencies w | |
| 2 function [w, u_hat] = four(x, u) | |
| 3 u_hat = fft(u); | |
| 4 | |
| 5 N = length(x); | |
| 6 L = x(end) - x(1); | |
| 7 | |
| 8 k = shift_k(0:N-1); | |
| 9 | |
| 10 u_hat = fftshift(u_hat); | |
| 11 | |
| 12 dw = 2*pi/L; | |
| 13 w = dw*k; | |
| 14 end | |
| 15 | |
| 16 function k_shifted = shift_k(k) | |
| 17 N = length(k); | |
| 18 k_shifted = [-floor(N/2):-1, 0, 1:ceil(N/2)-1]; | |
| 19 end |