PPS zaszumienie sygnalu sinusoidalnego

%skrypt sin_szum
%zaszumianie sygnalu sinusoidalnego

clf
f=5;
t=0:0.001:0.999;
x=2*sin(2*pi*f*t);
szum=randn(1,length(x));
y=x+szum;

subplot(3,1,1);
plot(t,x);
xlabel('czas [s]');
ylabel('sygnal 1');
subplot(3,1,2);
plot(t,szum);
xlabel('czas [s]');
ylabel('szum bialy');
subplot(3,1,3);
plot(t,y);
xlabel('czas [s]');
ylabel('sygnal + szum');