t2=(0:0.1:40);
f2=(0:0.025:10);
y4=square(t2);
widmo_sqr = fft(y4);
widmo_sqr_amp = abs(widmo_sqr);
widmo_sqr_faz = angle(widmo_sqr);
figure;
plot(t2, y4);
title('Prostokąt');
xlabel('sec');
ylabel('amp');
figure;
plot(f2, widmo_sqr_amp);
title('Widmo sygnału prostokątnego');
xlabel('freq');
ylabel('amp');