Répondre à la discussion
Affichage des résultats 1 à 3 sur 3

Profilage géographique pour trouver un assassin ,densité de probabilité avec Maple



  1. #1
    seb234

    Smile Profilage géographique pour trouver un assassin ,densité de probabilité avec Maple


    ------

    Bonjour,

    Le but de mon TIPE est de trouver les coordonnées de l'habitat d'un criminel à partir des coordonnées des lieux de meurtres qu'il a commis. Pour cela, j'utilise une formule de densité de probabilité (probabilité pour une surface). Cette formule peut être assimilée à une fonction à deux variables. Cependant, elle admet des conditions. POur ce faire, je crée la fonction à partir d'une procédure avec des boucles if. Pour faire afficher les zones de probabilité, j'utilise contourplot. Cependant, j'avais comme message d'erreur contourplot ne pouvait pas faire d'évaluation en booléen (à cause des if). Je rajoute alors la ligne suivante pour transformer ma procédure en fonction
    Prob:=(x,y)->Prob(x,y); et la commande whattype me confirme que Prob est bien une fonction. Je réessaie contourplot et là je reçois le message d'erreur suivant Error, (in plot/iplot2d/levelcurve) could not evaluate expression
    Et c'est là que je n'arrive pas à comprendre et que j'ai besoin d'aide
    Alors que Prob est une fonction, il n'arrive pas à la tracer !
    Pour vous aider, voici l'algorithme
    c:=9: Nombre de meurtres
    > x[1]:=2.15:y[1]:=-0.46: Coordonnées des lieux de meurtres
    > x[2]:=2.83:y[2]:=-1.16:
    > x[3]:=2.41:y[3]:=-2.53:
    > x[4]:=2.96:y[4]:=-3.03:
    > x[5]:=2.81:y[5]:=-3.35:
    > x[6]:=2.32:y[6]:=-3.28:
    > x[7]:=1.82:y[7]:=-3.52:
    > x[8]:=1.82:y[8]:=-3.55:
    > x[9]:=2.46:y[9]:=-7.13:
    > k:=(1/c): Constante
    > p:=4.13: q:= -7.18: Taille de la carte

    > B:=0.2: Constante (zone tampon)
    > f:=0.2: g:=0.4: Constantes
    Procédure pour créer la fonction
    > Prob:=proc(x,y)
    > local m;
    > for m from 1 to c do
    > if (abs(x-x[m])+abs(y-y[m])=0)
    > then 0;
    > else
    > if (abs(x-x[m])+abs(y-y[m])<=B)
    > then k*sum((B^(g-f))/((2*B-abs(x-x[n])-abs(y-y[n]))^g),n=1..c);
    > else k*sum(1/(abs(x-x[n])+abs(y-y[n]))^f,n=1..c);
    > fi;
    > fi;
    > od;
    > end;
    >
    > Prob:=(x,y)->Prob(x,y);
    > with(plots):
    > contourplot(eval(Prob(x,y)),x= 0..p,y=q..0,filled=true);
    > plot3d(eval(Prob(x,y)),x=0..p, y=q..0);

    -----

  2. #2
    seb234

    Re : Profilage géographique pour trouver un assassin ,densité de probabilité avec Maple

    Bonjour,

    Au niveau de la procédure, j'ai rectifié quelques trucs car il manquait quelquechose
    Voici ma nouvelle procédure
    Prob:=proc()
    > local m;
    > global x,y,p,q,B,c,f,g,k,Rossmo,K;
    > x[1]:=7.11:y[1]:=-0.30:
    > x[2]:=3.05:y[2]:=-7.07:
    > x[3]:=3.22:y[3]:=-5.27:
    > x[4]:=2.93:y[4]:=-5.05:
    > x[5]:=2.36:y[5]:=-5.82:
    > x[6]:=2.37:y[6]:=-5.74:
    > x[7]:=1.77:y[7]:=-5.54:
    > x[8]:=1.63:y[8]:=-5.39:
    > x[9]:=0.14:y[9]:=-5.98:
    > p:=7.23: q:= -7.12:
    > B:=0.6:
    > f:=0.02: g:=0.03:
    > c:=9:
    > for m from 1 to c do
    > if ((abs(s-x[m])+abs(t-y[m]))<=B and (abs(s-x[m])+abs(t-y[m]))>0)
    > then K[m]:=0: B[m]:=B:
    > elif ((abs(s-x[m])+abs(t-y[m]))>B and (abs(s-x[m])+abs(t-y[m]))<2*B)
    > then K[m]:=1: B[m]:=B:
    > elif ((abs(s-x[m])+abs(t-y[m]))>2*B)
    > then K[m]:=1: B[m]:=B:
    > else K[m]:=0: B[m]:=0:
    > fi;
    > od;
    > k:=(1/c);
    > Rossmo:=(s,t)->k*sum((K[n]/(abs(s-x[n])+abs(t-y[n]))^f)+((1-K[n])*(B[n]^(g-f)))/((2*B[n]-abs(s-x[n])-abs(t-y[n]))^g),n=1..c); with(plots):
    > contourplot(Rossmo(s,t),s=0..p ,t=q..0,filled=true);
    > plot3d(Rossmo(s,t),s=0..p,t=q. .0);
    > end;


    Cependant quand je fais Prob(9,7.23,-7.12,3,0.2,0.4);
    , il me répond toujours le même refrain Error, (in Prob) cannot evaluate boolean
    En même temps, pour mieux comprendre, je vous donne le lien où vous pourrez trouver la formule de densité de probabilité que l'on utilise dans le profilage géographique


    Ceux qui peuvent m'aider à comprendre sont les bienvenus !

  3. #3
    seb234

    Re : Profilage géographique pour trouver un assassin ,densité de probabilité avec Maple

    Bonjour,

    Au niveau de la procédure, j'ai rectifié quelques trucs car il manquait quelquechose
    Voici ma nouvelle procédure
    Prob:=proc()
    > local m;
    > global x,y,p,q,B,c,f,g,k,Rossmo,K;
    > x[1]:=7.11:y[1]:=-0.30:
    > x[2]:=3.05:y[2]:=-7.07:
    > x[3]:=3.22:y[3]:=-5.27:
    > x[4]:=2.93:y[4]:=-5.05:
    > x[5]:=2.36:y[5]:=-5.82:
    > x[6]:=2.37:y[6]:=-5.74:
    > x[7]:=1.77:y[7]:=-5.54:
    > x[8]:=1.63:y[8]:=-5.39:
    > x[9]:=0.14:y[9]:=-5.98:
    > p:=7.23: q:= -7.12:
    > B:=0.6:
    > f:=0.02: g:=0.03:
    > c:=9:
    > for m from 1 to c do
    > if ((abs(s-x[m])+abs(t-y[m]))<=B and (abs(s-x[m])+abs(t-y[m]))>0)
    > then K[m]:=0: B[m]:=B:
    > elif ((abs(s-x[m])+abs(t-y[m]))>B and (abs(s-x[m])+abs(t-y[m]))<2*B)
    > then K[m]:=1: B[m]:=B:
    > elif ((abs(s-x[m])+abs(t-y[m]))>2*B)
    > then K[m]:=1: B[m]:=B:
    > else K[m]:=0: B[m]:=0:
    > fi;
    > od;
    > k:=(1/c);
    > Rossmo:=(s,t)->k*sum((K[n]/(abs(s-x[n])+abs(t-y[n]))^f)+((1-K[n])*(B[n]^(g-f)))/((2*B[n]-abs(s-x[n])-abs(t-y[n]))^g),n=1..c); with(plots):
    > contourplot(Rossmo(s,t),s=0..p ,t=q..0,filled=true);
    > plot3d(Rossmo(s,t),s=0..p,t=q. .0);
    > end;


    Cependant quand je fais Prob(9,7.23,-7.12,3,0.2,0.4);
    , il me répond toujours le même refrain Error, (in Prob) cannot evaluate boolean
    En même temps, pour mieux comprendre, je vous donne le lien où vous pourrez trouver la formule de densité de probabilité que l'on utilise dans le profilage géographique


    Ceux qui peuvent m'aider à comprendre sont les bienvenus !

Sur le même thème :

Discussions similaires

  1. Densité probabilité gaussienne
    Par hammer313 dans le forum Mathématiques du supérieur
    Réponses: 0
    Dernier message: 19/01/2008, 12h20
  2. Fonction densité de probabilité
    Par x1999 dans le forum Mathématiques du supérieur
    Réponses: 2
    Dernier message: 17/12/2007, 11h43
  3. probabilité de présence et densité radiale
    Par Jabba le Hutt dans le forum Physique
    Réponses: 4
    Dernier message: 04/10/2007, 20h55
  4. Densité radiale de probabilité de présence
    Par epsilone88 dans le forum Chimie
    Réponses: 6
    Dernier message: 06/11/2006, 02h36
  5. densité de probabilité gaussienne
    Par Seirios dans le forum Mathématiques du collège et du lycée
    Réponses: 13
    Dernier message: 25/08/2006, 19h37