Laborka PPS 4

zadanie = 3;

if zadanie == 1
    % uno
    N = 1024; Fs = 100; n = (0:N - 1)'/Fs;
    x = sin(2 * pi * 13 * n) + 2 * sin(2 * pi * 37 * n) + 0.001 * randn(N, 1);
    X = fft(x);
    plotspec(X, Fs, 1);
elseif zadanie == 2
    % dos
    m = (0:1/Fs:10)';
    y = sin(2 * pi * 23 * m);
    Y = fft(y);
    plotspec(Y, Fs, 1);
elseif zadanie == 3
    % tres
    subplot(1, 1, 1);
    x = rand(50, 1); %plotspec(fft(x));
    y = x - mean(x);
    X = fft(x); Y = fft(y);
    plotspec([X, Y]);
elseif zadanie == 4
    % quattro
    N = 200; Fs = 100; n = (0:N - 1)'/Fs;
    x = sin(2 * pi * 20 * n); X = fft(x, 4096);
    %plotspec(X, Fs);
    N = 4095; w = [ones(N, 1); zeros(4096 - N, 1)]; W = fft(w);
    %plotspec([X, W], Fs);
    plotspec(W, Fs);
elseif zadanie == 5
    N = 2048; Fs = 100; n = (0:N - 1)'/Fs;
    x = sin(2 * pi * 27.1 * n) + 0.25 * sin(2 * pi * 27 * n);
    %w = gausswin(N, 1);
    w = hamming(N);
    y = x .* w;
    figure(1); plot(n, x, n, w, n, y);
    X = fft(x); Y = fft(y); W = fft(w);
    figure(2); plotspec([X, W, Y], Fs);
end