Apprendre Maple

[ronan.dupuy]

Voila mon problème :

Je cherche a modéliser sur maple ce que percoit un radar doppler et par la suite le champ de vecteurs du vent par recoupement des données acquises par plusieurs radars...

Le problème c'est que je n'arrive pas a afficher un champ de vecteurs (fonction fieldplot dans la librairie plot) défini par une procédure...

Du coup je me rabats sur un affichage de sequences de procédures afin d'afficher les vecteurs 1 par 1, ce qui est moins elegant mais qui revient a peu pres au meme, mais la encore j'ai des problèmes : on m'annonce que des lignes ne sont pas définies...

Si quelqu'un pouvait jeter un coup d'oeil a mon programme et trouver l'erreur (moi je n'y arrive pas, et pourtant je tourne, retourne et reretourne le problème dans tous les sens...), cela me serait d'un grand secours !!!!

Merci a vous...

> restart:
> with(geometry):
> with (plots):
Warning, the name changecoords has been redefined
définition des coordonnées des vecteurs (fonctions prises au hasard)
> Vx:=-3*x;
> Vy:=sqrt (y**2)+18;
calcul du champ percu par le radar 1
> a1:=arctan(Vy/Vx):
> b1:=arctan(y/x):
> c1:=sqrt((a1-b1)**2):
> Nx:=sqrt(Vx**2+Vy**2)*cos(c1):
> Fx:=-Nx*(x/sqrt(x**2+y**2)):
> Fy:=-Nx*(y/sqrt(x**2+y**2)):
calcul du champ percu par le radar 2
> x2:=x+7:
> y2:=y+7:
> Vx2:=-3*(x+7):
> Vy2:=sqrt ((y+7)**2)+18:
> a2:=arctan(Vy2/Vx2):
> b2:=arctan((y2)/(x2)):
> c2:=sqrt((a2-b2)**2):
> Nx2:=sqrt(Vx2**2+Vy2**2)*cos(c2):
> Fx2:=-Nx2*((x2)/sqrt(x2**2+y2**2)):
> Fy2:=-Nx2*(y2/sqrt(x2**2+y2**2)):
procédure pour retrouver le vrai champ
> recoup:= proc(x,y)
> local Vx3,Vy3,ligne;
> point(O,[x,y]):point(F1,[x+Vx,y+Vy]):point(F2,[x+Vx2,y+Vy2]):
> line(d1,[O,F1]):line(d2,[O,F2]):
> PerpendicularLine(p1,O,d1):PerpendicularLine(p2,O,d2):
> intersection(i,p1,p2):
> Vx3:=coordinates(i)[1]:
> Vy3:=coordinates(i)[2]:
> ligne:=plot([[x,y][x+2,y+3]]):
> end proc:
>
> recoup(2,3);
geometry/checkline: One of the following conditions must be satisfied 3*x <> 0 -18-(y^2)^(1/2) <> 0
Error, (in geometry/checkline) not enough information: the line is not defined
> display([seq(seq(recoup(i,j),j=1..10),i=1..10)]);
geometry/checkline: One of the following conditions must be satisfied 3*x <> 0 -18-(y^2)^(1/2) <> 0
Error, (in geometry/checkline) not enough information: the line is not defined

[ronan.dupuy]

un grand merci !!!
votre conseil est effectivement le "dépannage" que je cherchais deseperement...