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

[invitedb5b75cd]

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);
-----

[invitedb5b75cd]

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 !

[invitedb5b75cd]

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 !