%% 07/03/2012
% Confronto tutte le prove
% prove a +-1g +-0.5g e +-0.375g
% R=500,1000,1500,2000,3000,4000 ohm
% R=1,2,3,4,5,6,7,8,9,9.86 kohm
% R=10,20,30,40,50,60,70,80,90,100 kohm
% m=0,2.4,7.8,15.6 grammi

close all
clear all
clc


% vettore masse
m=[0 2.4 7.8 15.6]; %g
f=[120 75 50 40]; %Hz frequenze nominali


%lista  file di input e rispettive funzioni del datasheet
%funzioni per 1g e 0.5g
lista1=struct('funzione',{'-0.0480*x.^2+0.5338*x';'-0.0446*x.^2+0.2897*x';...
                         '-0.0353*x.^2+0.5861*x';'-0.0353*x.^2+0.2687*x';...
                         '-0.0412*x.^2+1.0892*x';'-0.0334*x.^2+0.5256*x';...
                         '-0.0413*x.^2+1.2357*x';'-0.0408*x.^2+0.6919*x'});

%funzioni per 0.5g e 0.375g                   
lista2=struct('funzione',{'-0.0446*x.^2+0.2897*x';'-0.0404*x.^2+0.2119*x';...
                         '-0.0353*x.^2+0.2687*x';'-0.0405*x.^2+0.2188*x';...
                         '-0.0334*x.^2+0.5256*x';'-0.0308*x.^2+0.3787*x';...
                         '-0.0407*x.^2+0.6902*x';'-0.0375*x.^2+0.5077*x'});

 
figure('position',[100,100,800,600])

x=[0:.1:30];

dV1=load('tensioni1.dat');
dV2=load('tensioni2.dat');
dV3=load('tensioni3.dat');
P1=load('potenze1.dat');
P2=load('potenze2.dat');
P3=load('potenze3.dat');

%% plot curve datasheet e risultati sperimentali

col=['brbrbrbrbrbrbrbr'];
for i=1:length(lista1)
    
    funz=inline(lista1(i).funzione);
    figure(1)    
    subplot(2,2,(i+mod(i,2))/2)
    hold on
    plot(x,funz(x),col(i),'linewidth',2)
    

    plot(dV1(i,:),P1(i,:),[col(i) 'v'])
    plot(dV2(i,:),P2(i,:),[col(i) 'o'])

end



col=['rgrgrgrgrgrgrgrg'];
for i=1:length(lista2)
    
    funz=inline(lista2(i).funzione);
    figure(1)    
    subplot(2,2,(i+mod(i,2))/2)
    hold on
    plot(x,funz(x),col(i),'linewidth',2)
    
        plot(dV3(i,:),P3(i,:),[col(i) 'd'])
    
    ylim([0,inf])
    grid on
    
    titolo=sprintf('m= %.1f g, f_0~ %.1f Hz',m((i+mod(i,2))/2),f((i+mod(i,2))/2));
    title(titolo);
    
    xlabel('V [V_R_M_S] ')
    ylabel('P [mW]')
    
    
end

print -painters -dpdf confronto_prove.pdf