Discussion :

[DJEcalcul]

Bonjour a tous,

Quelqu'un pourrait-il me renseigner sur comment trace une "tranche" de gaussienne. C'est a dire j'aimerai avoir une equation qui decrive un gradient centre en un point C, en gros c'est une gaussienne mais une tranche. J'aimerai avoir l'image de gauche, quel en est donc l'equation ?



Pour le moment j'utilise cette equation classique d'une gaussienne, mais qui me donne pas ce que je veux.

Code :

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

1
2
3
4
5
6
7
x = -3:1:3;
y = -3:1:3;
xc=0;
yc=0;
exponent = 10.*(2.*(x.^2) + 5.*(y.^2) );
val       = (exp(-exponent));
imagesc(val)

Merci a vous.

[OlivANoyaux]

Hum... Et en essayant : f(x,y) = exp(-(x²+y²)) ça ne donne pas ce que tu souhaites ? Un autre moyen de trouver des fonctions qui conviennent (tu n'auras pas d'expression analytique mais ça se fait numériquement) est de résoudre l'équation de Poisson : - alpha * Laplacien(u) = 0, avec u = 0 au bord. En jouant sur la valeur alpha tu auras une solution qui "diffuse" plus ou moins.

[DJEcalcul]

Hum... Et en essayant : f(x,y) = exp(-(x²+y²)) ça ne donne pas ce que tu souhaites ? Un autre moyen de trouver des fonctions qui conviennent (tu n'auras pas d'expression analytique mais ça se fait numériquement) est de résoudre l'équation de Poisson : - alpha * Laplacien(u) = 0, avec u = 0 au bord. En jouant sur la valeur alpha tu auras une solution qui "diffuse" plus ou moins.

Merci a toi oliv,

En fait en utilisant "imagesc" avec la fonction classique gaussienne, qu tu cites plus haut, ca donne ce que je veux.

Merci a toi.

Voici le code:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
clear all
x=[-2:0.1:2];
y=[-2:0.1:2];
[X,Y]=meshgrid(x,y);
 
std=1;
 
model=exp( (-X.^2-Y.^2)/std^2 );
imagesc(model)
title('Model')
colorbar