% forslag til løsning av Ukesoppgaver 4 oppgave 6
N=1024;
x = 2*pi*(0:(N-1))/N;
f = x;
f_ft = ifft(f);
loglog (0:(N-1),abs(f_ft),'b')
hold on
g = abs (x - pi);
g_ft = ifft(g);
loglog (0:(N-1),abs(g_ft),'r')
h = cos (x);
h_ft = ifft(h);
loglog (0:(N-1),abs(h_ft),'g')
xfeil = linspace(0,2*pi,N);
hfeil = cos (xfeil);
hfeil_ft = ifft(hfeil);
loglog ((N-1)/N*(0:(N-1)),abs(hfeil_ft),'c')
loglog ([1, N/2], [1, N/2].^(-1), '--k')
loglog ([1, N/2], [1, N/2].^(-2), '--y')
% Vi ser at:
% f_ft har k^(-1), parallell med svart stiplet linje
% g_ft har k^(-2), "toppene" parallell med gul stiplet linje
% h_ft har k^(-uendelig), verdier under 10e-15 er uendelig små
% Alt som vi forventet.
% For feilaktig x-akse ser vi at h_ft blir k^(-2), også som forventet, dog
% parallellforskjøvet litt ned.
hold off
