%
% corso COSTRUZIONI OPTOELETTRONICHE
% data 04/11/2004
%
% 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 = k * f_vibr * amp_vibr_picco
%


lambda = 0.6328; % [um]
k = 18.75; % [um^-1]



freq = [			10			20		50		100	115	130	145	160	175	200 	500	1000]; % [Hz]
V_drive = [		32			32		33		30		36		48		52		41		31		52		350	1400]; % [mW] picco-picco
f_max_fft = [	.6			1.14	2.898	7.66	12.3	20.6	22.65	15.9	10.56	12.9	23		55]; % [kHz] 
f_max_fft = f_max_fft * 1e3; % [Hz]

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

figure(1)
clf
set(1,'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(2)
clf
set(2,'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

