clc
clear all;
T = 2;                                          %okres
U = 1;                                          %amplituda
N = 100;                                        %liczba próbek
f= 1/T;                                         %częstotliwość
w = 2* pi*f;                                    %pulsacja
t = -T/2:T/N:T/2-T/N;                           %przedział czasu
for k=1:N
    u = abs((2*U*t)/T);                         %funkcja u1(t)
end
t1 = [t, t+T, t+2*T];
u1 = [u, u, u];
subplot(2,2,1);                                 %rozmieszczenie wykresu
plot(t1, u1);                                   %stworzenie wykresu funkcji
title("Sygnał u1 w postaci funkcji czasu");     %tytuł wykresu
xlabel("Czas [s]");                             %podpis osi X
ylabel("Amplituda [V]");                        %podpis osi Y
grid on;                                        %włączenie siatki
a=fft(u1);
m=abs(a);
n=0:N/2;
f_a=n*f/(3*N/2);
subplot(2,2,2);
stem(f_a(1:N/2),m(1:N/2),'.');
u2 = U/2-4*U/pi*pi*(cos(w*t)+cos(3*w*t)/3*3+cos(5*w*t)/5*5+cos(7*w*t)/7*7+cos(9*w*t)/9*9+cos(11*w*t)/11*11);
u2 = [u2, u2, u2];
subplot(2,2,3);                                 %rozmieszczenie wykresu
plot(t1, u2);