Listato jitter2.m

function [ z ] = jitter2( n )

%JITTER2 Effetto del jitter di campionamento su un'onda a dente di sega

 

%Inizializzazione

Fs=10000;

D=0.05;

x=[0:1/Fs:1-1/Fs];

N=length(x);

s=sawtooth(2*pi*3*x);

 

%Effetto del jitter di campionamento

for i=1:n

jti=sqrt(1/1e4)*randn;

sJ(i,:)=sawtooth(2*pi*3*(x+D+jti));

plot(x,s);

axis([0 1 -1.1 1.1]);

hold on;

plot(x,sJ(i,:),'m');

end

 

title('Jitter di campionamento');

hold off;

pause;

 

%Media sulle acquisizioni

for i=1:N

sM(i)=sum(sJ(:,i))/n;

end

 

plot(x,s);

axis([0 1 -1.1 1.1]);

hold on;

plot(x,sM,'m');

title('Media sulle acquisizioni');

hold off;

 

pause;

close;