Listato fase.m

function [ z ] = fase( n )

%FASE Effetto del rumore di fase su una sinusoide

 

%Inizializzazione

Fs=1000;

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

N=length(x);

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

 

%Visualizzazione del rumore di fase

nti=sqrt(1/1e4)*randn(1,N);

sE=sin(2*pi*3*(x+nti));

 

plot(x,s);

axis([0 1 -1.1 1.1]);

hold on;

plot(x,sE,'m');

title('Rumore di fase');

hold off;

 

pause;

 

%Inizializzazione 2

Fs=10000;

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

N=length(x);

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

 

%Aggiunta del rumore di fase

for i=1:n

nti=sqrt(1/1e4)*randn(1,N);

sR(i,:)=sin(2*pi*3*(x+nti));

end

 

%Media sulle acquisizioni

for i=1:N

sM(i)=sum(sR(:,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;