%
% corso COSTRUZIONI OPTOELETTRONICHE
% data 04/11/2005
%
% misura della risposta in frequenza di un altoparlante tramite un interferometro di Michelson
%
% ** metodo della misura della massima frequenza nello spettro FFT del segnale interferometrico **
%
% formula:  f_max_fft = C * f_vibr * amp_vibr_picco-picco
%


lambda = 0.6328; % [um]
C = 10; % [um^-1]



freq = [			10			20		50		75		100	125	150	175	200 	500	1000	2000]; % [Hz]
V_drive = [		114		115	60		24		11.5	10		40		27		29		280	830	1700]; % [mW] picco-picco
f_max_fft = [	5.03		9.08	13		8.57	5.44	5.31	22.4	11.77	10.56	30.5	61.4	30]; % [kHz] 
f_max_fft = f_max_fft * 1e3; % [Hz]

spost = f_max_fft ./ freq / C; % [um] picco-picco
resp = spost ./ V_drive; % [um/mV]

figure(3)
clf
set(3,'color',[1 1 1])
plot(freq,resp,'-b','linewidth',2)
hold on
plot(freq,resp,'ob','MarkerSize',8,'LineWidth',2)
grid on
xlabel('Frequency  [Hz]')
ylabel('Loudspeaker Responsivity  [\mum/mV]')
sizefig(14) % call alla funzione sizefig.m DA SCARICARE ED INSERIRE NELLA DIRECTORY CORRENTE
zoom on

figure(4)
clf
set(4,'color',[1 1 1])
loglog(freq,resp,'-b','linewidth',2)
hold on
loglog(freq,resp,'ob','MarkerSize',8,'LineWidth',2)
grid on
xlabel('Frequency  [Hz]')
ylabel('Loudspeaker Responsivity  [\mum/mV]')
sizefig(14) % call alla funzione sizefig.m DA SCARICARE ED INSERIRE NELLA DIRECTORY CORRENTE
axis([10 1000 0.001 1])
zoom on

