f=2;
T=1/f;
t=[0:0.01:5];
y=5+6.*sqrt(3).*exp(-2.*t).*cos(3.*t);
AA=[2,1,5,6,7];
A=mean(AA);
omega=2*pi*f;
z=A.*sin(omega.*t);
f=y+z;
szum=randn(size(t));
g=f+szum;

figure;
grid on;
subplot(5,1,1);
plot(t,y);
title('Sygnal y(t)');
ylabel('czas');
xlabel('amplituda');

subplot(5,1,2);
plot(t,z);
title('Sygnal z(t)');
ylabel('czas');
xlabel('amplituda');

subplot(5,1,3);
plot(t,f)
title('Sygnal f(t)');
ylabel('czas');
xlabel('amplituda');

subplot(5,1,4);
plot(t,g);
title('Sygnal g(t)');
ylabel('czas');
xlabel('amplituda');

subplot(5,1,5);
stem(t,z);
title('Sygnal sprobkowany z(t)');
ylabel('czas');
xlabel('amplituda');

