Discussion :

[Jragk]

Bonjour,

Je vous explique mon problème : j'ai besoin de générer une équation d'un socle pour faire tourner des simulations d'écoulement de glace dessus. Ce socle, je l'ai visualisé en théorie, je l'ai dessiné, j'ai récupéré les coordonnées de plusieurs points, et j'ai obtenu une équation polynomiale de ma surface via sftool. Pour info, voilà la forme grossière de ma surface :

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
clear
 
x=[0;20000;40000;60000;80000;100000;120000;140000;160000;180000;200000];
y=[0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 15000 20000 25000 30000];
 
     %Coefficients (with 95% confidence bounds):
       p00 =       9.754  %(-27.73, 47.24)
       p10 =    -0.01003  %(-0.01034, -0.009723)
       p01 =    -0.03569  %(-0.06244, -0.008948)
       p11 =   8.702e-08  %(-8.308e-08, 2.571e-07)
       p02 =   1.153e-05  %(5.714e-06, 1.735e-05)
       p12 =   6.269e-11  %(3.609e-11, 8.929e-11)
       p03 =  -1.167e-09  %(-1.661e-09, -6.728e-10)
       p13 =  -4.066e-15  %(-5.483e-15, -2.649e-15)
       p04 =   4.639e-14  %(2.826e-14, 6.451e-14)
       p14 =   6.894e-20  %(4.511e-20, 9.277e-20)
       p05 =  -6.334e-19  %(-8.738e-19, -3.93e-19)
 
       %Linear model Poly15:
       for j=1:length(y)
           for i=1:length(x)
                z(i,j) = p00 + p10*x(i) + p01*y(j) + p11*x(i)*y(j) + p02*y(j)^2 + p12*x(i)*y(j)^2 + p03*y(j)^3 + p13*x(i)*y(j)^3 + p04*y(j)^4 + p14*x(i)*y(j)^4 + p05*y(j)^5;
           end
       end
 
       surf(y,x,z)
       xlabel('x');
       ylabel('y');
       zlabel('z');

Le soucis, c'est que la pente située sur la coordonnées x=0 (pente la plus forte) est sensé être un axe de symétrie. Aussi j'aimerai savoir s'il est possible d'imposer lors du fit une contrainte sur la tangente au point y=2.10^(5) et x=0, par exemple, pour forcer le socle à être bien à plat afin de satisfaire la symétrie.

Je suppose que ça doit être possible via cet outil, mais je n'ai pas trouvé jusque là.

Merci !

[Jragk]

Ahah, je savais que je vous avais collé

Bon, pour info, j'ai finalement opté pour une autre fonction pour mon problème, une fonction "tanh".

Mais ça n'ôte pas ma question, qui s'applique aussi à la fonction cftool.

Voilà, si jamais quelqu'un passant par ici connait la réponse, je suis toujours intéressé