1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55
| [d]=xlsread('Projet1.xlsx','Sheet1');
mv1 = d(:,2);
mv2 = d(:,3);
vol = d(:,1);
lignesatrou = any(isnan([mv1,mv2]),2);
vol(lignesatrou) = [];
mv1(lignesatrou) = [];
mv2(lignesatrou) = [];
pas_de_trou = ~isnan(mv1);
mv1 = interp1(vol(pas_de_trou),mv1(pas_de_trou), vol);
pas_de_trou = ~isnan(mv2);
mv2 = interp1(vol(pas_de_trou),mv2(pas_de_trou), vol);
outliers =([mv1 mv2] - mean([mv1;mv2]) >= 3*std([mv1;mv2]));
mv1(outliers(:,1))= NaN;
mv2(outliers(:,2))= NaN;
span = 0;
lissage = conv([mv1;mv2],ones(span,1)/span, 'same');
moyenne = mean([mv1 mv2]')';
ecart_type = std([mv1 mv2]')';
intervallevraisemblance = ecart_type/sqrt(2)*t(0.05,1);
errorbar(vol,moyenne,intervallevraisemblance);
title('Dosage halogénure [KI et KCl]')
xlabel('Volume [mL]')
ylabel('Potentiel [mV]')
pv1 = polyfit(vol(10:11), mv1(10:11),1);
ph1 = polyfit(vol(15:22), mv1(15:22),1);
pv2 = polyfit(vol(10:11), mv2(10:11),1);
ph2 = polyfit(vol(15:22), mv2(15:22),1);
x1 = (ph1(2)-pv1(2))/(pv1(1)-ph1(1));
x2 = (ph2(2)-pv2(2))/(pv2(1)-ph2(1));
moyennePtEqui = mean(x1:x2);
d1_exp1 = gradient(mv1, vol);
d2_exp1 = gradient(d1_exp1, vol);
tmp = [ vol d1_exp1 d2_exp1 ];
d1_exp2 = gradient(mv2, vol);
d2_exp2 = gradient(d1_exp1, vol);
tmp = [ vol d1_exp2 d2_exp2 ]; |
Partager