Listato sweep.m

function [ z ] = sweep(  )

%SWEEP Effetto spettrale dell'errore di velocità di sweep

 

%Inizializzazione

Fs=10000;

D=0.05;

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

N=length(x);

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

 

%Spostamento rigido dello spettro a causa

%dell'errore di velocità di sweep

avi=avifile('spettro.avi');

avi.Compression='None';

avi.FPS=10;

 

cs=fft(s);

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

 

for i=0:1/100:3

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

csE=fft(sE);

stem(f(1:16),abs(cs(1:16))/N,'-b.');

title('Spostamento rigido dello spettro');

xlabel('Frequenza [Hz]');

axis([0 15 0 .6]);

axis manual;

hold on;

stem(f(1:16),abs(csE(1:16))/N,'-m.');

hold off;

F=getframe;

avi=addframe(avi,F);

end

hold off;

avi=close(avi);

 

pause;

close;