% forslag til løsning av Ukesoppgaver 5 oppgave 1
a=[0 1 0 -1];
N = length(a);
x = 2*pi*(0:(N-1))/N;
plot (x,a,'*')
hold on
aft = ifft(a); % Denne deler på N, så vi trenger ikke tenke mer på det!

% riktig derivertformel
daft = aft.*[0:(N/2), (-N/2+1):-1]*(-i);
da = fft(daft);
plot (x,da,'o');

% interpolasjon
NN = 100;
aaft = [aft(1:(N/2)), zeros(1,NN-4), aft((N/2+1):N)];
aa = fft(aaft);
xx = 2*pi*(0:(NN-1))/NN;
plot (xx,aa,'r-')

% riktig derivertformel
daaft = aaft.*[0:NN/2, (-NN/2+1):-1]*(-i);
daa=fft(daaft);
plot (xx,daa,'b-')

% feil derivertformel, prøv selv!
% feildaaft = aaft.*[0:(N-1)]*(-i);
% feildaa=fft(feildaaft);
% plot (xx,feildaa,':')
hold off
