Mercurial > repos > public > sbplib
comparison four.m @ 707:0de70ec8bf60 feature/quantumTriangles
merge with feature/optim
| author | Ylva Rydin <ylva.rydin@telia.com> |
|---|---|
| date | Fri, 10 Nov 2017 14:22:56 +0100 |
| parents | f6ada1719420 |
| children |
comparison
equal
deleted
inserted
replaced
| 696:7c16b5af8d98 | 707:0de70ec8bf60 |
|---|---|
| 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 |