Zaszumianie sygnału

%Dane wejsciowe:
T=2;
N=1000;                     %czas trwania sygnalu
krok_t=T/N;                 %okres probkowania
t=0:krok_t:T-krok_t;        %wektor chwil czasow probkowania
a=2;                        %dekrement tlumienia
e=exp(-a*t);
A=17;
f=10;
x_s=A*sin(2*pi*f*t);
%x_s=square(2*pi*f*t);
%zaszumianie sygnalu szumem gaussowskim
y=x_s.*e;
y=y+1*randn(1,length(y));
subplot(3,1,1);
plot(t,e);
subplot(3,1,2);
plot(t,x_s);
subplot(3,1,3);
plot(t,y);