staka

% Definicija signala x[n]
n = 0:29;
x = (n + 1) .* (0.5).^n .* (n >= 6);

% Inicijalizacija promenljive za Furijeovu transformaciju
N = length(x);
X = zeros(1, N);

% Izracunavanje Furijeove transformacije
for k = 0:N-1
    X(k+1) = sum(x .* exp(-1j * 2 * pi * k * n / N));
end

% Prikaz rezultata
subplot(2, 1, 1);
stem(n, x);
title('Originalni signal x[n]');
xlabel('n');
ylabel('x[n]');

subplot(2, 1, 2);
stem(0:N-1, abs(X));
title('Diskretna Furijeova transformacija X[k]');
xlabel('k');
ylabel('|X[k]|');

% Prikaz faznog ugla
figure;
stem(0:N-1, angle(X));
title('Fazni ugao diskretne Furijeove transformacije');
xlabel('k');
ylabel('Phase angle (radians)');