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 56 57 58 59 60
|
close all
clear all
clc
format long
mat = [
0.000 0.000 0.000
0.100 -0.800 -0.250
0.200 -1.600 -0.500
0.300 -2.400 -0.750
0.400 -3.200 -1.000
0.500 -4.000 -1.250
0.600 -4.800 -1.500
0.700 -5.600 -1.750
0.800 -6.400 -2.000
0.900 -7.200 -2.250
1.000 -8.000 -2.500
1.100 -8.800 -2.750
1.200 -9.600 -3.000
1.300 -10.400 -3.250
1.400 -11.200 -3.500
1.500 -12.000 -3.750
1.600 -12.800 -4.000
1.700 -13.600 -4.250
1.800 -14.400 -4.500
1.900 -15.200 -4.750
2.000 -16.000 -5.000
2.100 -16.800 -5.250
2.200 -17.600 -5.500
2.300 -18.400 -5.750
2.400 -19.200 -6.000
2.500 -20.000 -6.250
];
abscisse = mat(:,1); [fin b] = size(abscisse);
m0 = mat(:,2);
m1 = mat(:,3);
plot(abscisse,m0,'r',abscisse,m1,'g')
grid on ;
f = polyfit(abscisse,m0,1)
r = polyval(f,abscisse);
plot(abscisse,m0,'o',abscisse,r,'-')
f2 = polyfit(abscisse,m1,1)
r2 = polyval(f2,abscisse);
plot(abscisse,m0,'go',abscisse,r,'g-',abscisse,m1,'ro',abscisse,r2,'r-')
syms x real ;
int1 = int((f(1,1)*x+f(1,2))*(f2(1,1)*x+f2(1,2)),abscisse(1,1),abscisse(fin,1)) * 1;
h=[int1]
syms x unreal ;
pause
csvwrite('delta.csv',h); |
Partager