Discussion :

[chouicha]

Bonjour à tous
je cherche à élaborer un modèle mathématique y=f(x)
C'est pour ce faire j'ai utilisé la fonction polyfit et polyval (MATLAB) le problème c'est que je n'arrive pas à comprendre le résultat obtenu.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
61
62
63
x =                           
  294.1624
  294.6455
  295.2028
  295.7994
  296.4475
  297.1439
  297.9043
  298.7370
  299.6505
  300.6718
  301.8202
  303.1392
  304.6991
  306.4747
  308.7776
y =
 
    0.0100
    0.0150
    0.0200
    0.0250
    0.0300
    0.0350
    0.0400
    0.0450
    0.0500
    0.0550
    0.0600
    0.0650
    0.0700
    0.0750
    0.0800
>> polyfit(x,y,length(x-1))
Warning: Polynomial is not unique; degree >= number of data points.
> In polyfit at 72
 
ans =
 
  1.0e-026 *
 
  Columns 1 through 15
 
    0.0000   -0.0000    0.0012   -0.1224         0         0         0         0         0         0         0         0         0         0         0
 
  Column 16
 
         0
 
Warning: Polynomial is not unique; degree >= number of data points.
> In polyfit at 72
 
ans =
 
  1.0e+011 *
 
  Columns 1 through 15
 
         0         0         0         0         0         0         0    6.1440   -2.1373    0.3125   -0.0249    0.0012   -0.0000    0.0000   -0.0000
 
  Column 16
 
    0.0000

[Jerome Briot]

La documentation de la fonction polyfit est pourtant claire :

p = polyfit(x,y,n) finds the coefficients of a polynomial p(x) of degree n

Par exemple, pour un polynôme de degré 2 :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
x = [294.1624
    294.6455
    295.2028
    295.7994
    296.4475
    297.1439
    297.9043
    298.7370
    299.6505
    300.6718
    301.8202
    303.1392
    304.6991
    306.4747
    308.7776];
 
y = [0.0100
    0.0150
    0.0200
    0.0250
    0.0300
    0.0350
    0.0400
    0.0450
    0.0500
    0.0550
    0.0600
    0.0650
    0.0700
    0.0750
    0.0800];
 
p = polyfit(x, y, 2);
 
xx = linspace(min(x), max(x), 100);
yy = polyval(p, xx);
 
figure
plot(x, y, 'r*', xx, yy, 'b-')

[chouicha]

Merci