distorsione

function [ z ] = distorsione( n )

%DISTORSIONE Spettro di una sinusoide distorta

 

%Inizializzazione

Fs=10000;

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

N=length(x);

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

 

%Distorsione

disty=distn(y,n);

 

%FFT

cDy=fft(disty);

 

%Vettore frequenza

f=(0:N-1)*Fs/N;

 

%Visualizzazione risultati

plot(x,y);

title('Segnale di ingresso');

pause;

 

plot(x,x);

hold on;

plot(x,distn(x,n),'r');

title('Distorsione non lineare');

hold off;

pause;

 

stem(f(1:30),abs(cDy(1:30))/N,'-b.');

title('Spettro del segnale distorto');

 

pause;

close;