C++

recherche algo balistique (bonus : jeux inside ! )

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : recherche algo balistique (bonus : jeux inside ! )

nico168

salut,

j'ai les constantes gravité, vent, masse de l'objet.
le vent est toujours parallele au sol, positif ou
negatif.

mon but ultime est d'avoir une fonction du type :

je lui transmet l'adresse d'une structure en entrée avec x, y, angle, velocité
ainsi qu'un timecode T.

la fonction me retourne un pointeur sur une structure contenant x,y,angle, velocité qu'aurra l'objet au timecode T.

a force de googliser j'ai reussi a trouver un algo qui pourrais me convenir, mais il ne contient pas la notion du vent.
je n'ai pas les competence en physique/mathematique pour modifier son algo.

si qq'un a vu un truc comme ca trainer qqpart (et GPL

merci

Message édité par nico168 le 05-01-2003 à 14:17:36

Publicité

sombresonge

Il faut resoudre l'equation différentiel du mouvement et ensuite tu implemente ca dans le langage que tu veux.

En gros tu resous

F1+F2+...+Fi = m.a = m.dV/dT = m.d²OM/dT²

[F1..Fi] etant l'expresion vectoriel des forces que subissent ton objet plus ya de force que tu prendra en compte plus la solution sera complexe.

m c la masse
V la vitesse
OM ton vecteur position cad le resultat que tu cherche
T le temps

Pour un truc simple tu peux prendre en compte comme force la pasenteur et la resistance de l'aire ca te donne une equation du type:

g - k.V/m = d²V/dT²
<=>
k.V/m + d²V/dT² = -g

ou k est le coeficient de frotement.

c une equation différentielle du second ordre avec second membre ya plus qu'à resoudre et t'aura t'as fonction de OM(T) cad f(t)->(x,y,z)

Taz

bisounours-codeur

mauvaise cat :pfff:

Kyle_Katarn

Pour le vent en fonction de la vitesse (en fait du nombre de Reynolds) ça sera une force en V ou en V²

nico168

houla....

je capte rien a se que je vous me raconter...
comme je le disais plus haut, je suis pas trop calé en math. je ne sais meme pas se qu'est une equation differentielles... alors pour se qui est de la resoudre...

soit l'algo suivant qui fait a eu pres se que je veux:

Code :

double g=-8, pi=3.141593, c=0.001; // define simulation constants -8 0.002
double fnf(double p, double q) {
return -c*p*sqrt(p*p+q*q);
}
double fng(double p, double q) {
return -c*q*sqrt(p*p+q*q)-g;
}
/////////////////////////////////////////////////////////////////////////////////////////////////////////
int ComputeBalistic( BalisticType *balistic )
{
#ifdef DEBUG_FONCNAME
fprintf(stderr,"ComputeBalistic\n" );
#endif
double t=0; // start time
double x=balistic->Start.X; // initial position
double y=balistic->Start.Y;
double v=balistic->Start.Velocity; // initial velocity
double a=balistic->Start.Angle; // fireing angle
double h=balistic->StepSize; // step size
double f1,g1,f2,g2,f3,g3,f4,g4,dp,dq,p,q;
int n;
BalisticDoubleRowType tmpRow;
BalisticIntRowType lastRow;
char truc[250];
lastRow.X=0;
lastRow.Y=0;
a=pi*a/180; // convert angle to rad
p=v*cos(a); // calc vert & horis. speeds.
q=v*sin(a);
for(n=1; n<=500000; n++)
{
f1=fnf(p,q);
g1=fng(p,q);
f2=fnf(p+h*f1/2, q+h*g1/2);
g2=fng(p+h*f1/2, q+h*g1/2);
f3=fnf(p+h*f2/2, q+h*g2/2);
g3=fng(p+h*f2/2, q+h*g2/2);
f4=fnf(p+h*f3, q+h*g3);
g4=fng(p+h*f3, q+h*g3);
dp=(f1+2*(f2+f3)+f4)/6;
dq=(g1+2*(g2+g3)+g4)/6;
x+=p*h+0.5*dp*h*h; // three-term-Taylor approximation
y+=q*h+0.5*dq*h*h;
p+=dp*h; q+=dq*h; t+=h;
if((lastRow.X!=(int)x)||(lastRow.Y!=(int)y)) //if it is drawn in another pixel
{
//fprintf(stderr,"added------------------" );
tmpRow.Timecode=t;
tmpRow.X=x;
tmpRow.Y=y;
tmpRow.Velocity=sqrt(p*p+q*q);
tmpRow.Angle=180*atan(q/p)/pi;
AddRowInBalistic(balistic, tmpRow);
lastRow.X=(int)x;
lastRow.Y=(int)y;
if( (IsFieldEmpty((int)x,(int)y)==0) ||((int)y>YReferenceSize)) //if new pixel is not empty
{
balistic->Impact.Timecode=h;
balistic->Impact.X=x;
balistic->Impact.Y=y;
balistic->Impact.Velocity=tmpRow.Velocity;
balistic->Impact.Angle=tmpRow.Angle;
break;
}
}
}
balistic->IsComputed=1;
return 1;
}

jeux : lui rajouter la notion du vent

Message édité par nico168 le 05-01-2003 à 14:52:25

verdoux

And I'm still waiting

Faut voir si c'est réellement pertinent avec le modèle que tu utilises.
Déjà la vitesse initiale d'un obus c'est typiquement 500 m/s pour un obus de canon de 155mm et plus de 1000m/s pour un obus antichar de 120mm.

Alors un petit vent de 50km/h (=13.9m/s) va pas vraiment modifier grand chose.

Et si tu recherches la précision (disons 50m pour une cible à 20km), il faudra aussi tenir compte des caractéristiques (pression, température) des couches atmosphériques traversées car l'obus va monter assez haut.

Message édité par verdoux le 05-01-2003 à 14:35:22

nico168

je ne cherche pas vraiment une precision.
tout ca c'est pour un jeux.
tu connais le jeux worms ?
tu deplace un petit personnage, tu choisis une arme, tu lui donne un angle de tir, une force, et tu tire.
tout ca se passe en 2d.
il peut y avoir du vent, par exemple si tu tire a 45° aux bazooka et qu'il y a un fort vent contraire, il peut te revenir dessus.

bref..j'utilise l'algo si dessus, tout marche bien...mais y'a pas de vent.

moi je verrais bien le vent comme :
il n'agit que sur le position X (parallele au sol)
X=X+(timecode actuel-timecode de depart)*vent.
il serais donc exprimé en pixel par unité de timecode ?
mais je ne sais pas si c'est correct et comment l'ajouter dans l'algo

Message édité par nico168 le 05-01-2003 à 14:57:37

LeGreg

Citation :

moi je verrais bien le vent comme :
il n'agit que sur le position X (parallele au sol)
X=X+(timecode actuel-timecode de depart)*vent.
il serais donc exprimé en pixel par unité de timecode ?
mais je ne sais pas si c'est correct et comment l'ajouter dans l'algo

non ce n'est pas correct.
Dans ton equation, tu as simplement vitesse de ton objet = vent
ce qui est correct quand il n'est soumis a aucune autre force
or ton objet a une vitesse propre est soumis a la gravite et il avance contre ou dans le vent.

j'ecrirais quelque chose comme:

x = x + vx * dt;
y = y + vy * dt;
vx = vx + f * (vvx-vx) * dt;
vy = vy - (f * vy + g) * dt;

avec f la force du vent rapportee a la masse
de ton objet (un objet dans une matiere legere subira
fortement l'influence du vent, contrairement a un objet
dans une matiere lourde). Ca depend aussi de l'aerodynamisme
et de l'angle d'attaque de l'objet dans l'air (qui influera egalement sur son moment mais on laisse tomber ca pour l'instant)

g est la gravite, c'est une constante si tu parcours de petites distances et elle vaut 9.81 m/(s*s)

dt est le temps qui s'ecoule entre deux calculs successifs de position, plus il est petit plus le calcul est precis. Le probleme avec cette methode c'est que pour arriver a l'instant T, il faut calculer une floppee de positions successives avant ce point.

Si tu ne veux pas faire ca ET si la vitesse du vent est une constante. Tu as la formule qui te donne la position directement en fonction de la date et qui decoule directement de l'equation differentielle mais.. le resultat est particulierement lourd
et fait intervenir des exponentielles si tu te contentes d'une force du vent dependant lineairement de la vitesse..
Je posterais la formule si besoin est..

LeGreg

---------------
voxel terrain render engine | animation mentor

nico168

merci pour vos reponses,

Citation :

non ce n'est pas correct.
Dans ton equation, tu as simplement vitesse de ton objet = vent

je me suis mal exprimé...dans mon equation...je supposais que X etait deja la position calculé pour le timecode actuel, je n'ai donc plus qu'a lui ajout l'action du vent. je vais donc reformuler :
X1= position X au timecode demandé sans l'action du vent
X2= position X au timecode demandé avec l'action du vent
Vent= exprimé en unité de distance/unité de timecode
sachant que X1 est calculé par l'algo plus haut, je pensais que :
X2=X1+(timecode actuel-timecode de depart)*vent

j'ai d'ailleur trouver ca sur le net :

Citation :

To take into account wind drag, add to the horizontal distance moved a small extra distance term proportional to the wind speed at each time step. Not quite physically realistic, but good enough I expect.

Message édité par nico168 le 07-01-2003 à 12:21:14

BettaSplendens

Tout cul tendu mérite son dû

nico168 a écrit :

tu connais le jeux worms ?

:love: :love: :love:

---------------
Tout cul tendu mérite son dû

Publicité

nico168

bon..j'ai mis une version windows pre-alpha de mon projets ici

http://c500.dyndns.org/test.zip (147ko)

les touches gauche et droite pour deplacer
les touches haut et bas pour modifier l'angle
les touches page up page down pour changer la velocité
la touche tab pour previsualiser le tir
la touche espace pour tirer (marche pas encore )

comme vous pourrez le constater...j'ai des petits pb avec le calcul des rebonds mais je desespere pas.

si ca peut vous permettre de mieux comprendre le pb.

je suis ouvert a tous conseils/remarques.

Message édité par nico168 le 07-01-2003 à 16:33:30

nico168

bon...je vois que mon projets a attiré des foules...

(p'ti up déguisé)

LeGreg

Citation :

ah ben si tu t'en fiches que ce soit physiquement correct
tu fais ce que tu veux effectivement :sweat:

Enfin le type qui a ecrit ca il ne s'est vraiment pas fatigue..

LeGreg

---------------
voxel terrain render engine | animation mentor

sombresonge

C klr que Worms c pas le summon du réalisme Moi qd j'ai vue que tu demandais une equation de ballistique je pensais que ct pour faire un trucs sérieux [:zaib3k]

BettaSplendens

Tout cul tendu mérite son dû

nico168 a écrit :

oups... les touches haut bas gauche droite ne font rien chez moi... sinon, ça commence pas mal...

---------------
Tout cul tendu mérite son dû

nico168

ouai mais bon, d'un autre coté, si les touches haut/bas/gauche/droite ne marchaient pas, tu n'as pas du voir grand choses

je comprends pas pourquoi, ca marche pas chez toi. 2 potes l'ont testé et ils n'ont pas ce pb.

HappyHarry

BeTtASpLeNdEnS a écrit :

oups... les touches haut bas gauche droite ne font rien chez moi... sinon, ça commence pas mal...

nico168

HappyHarry a écrit :

tu veux dire que toi aussi les touches de directions ne marche pas ?

HappyHarry

nico168 a écrit :

tu veux dire que toi aussi les touches de directions ne marche pas ?

bah oué

BettaSplendens

Tout cul tendu mérite son dû

et perso, j'ai essayé les touches standard.. et aussi en passant par le pavé numérique (activé comme déséactivé)

bah.. j'sais pas trop d'où ça vient ton pb là faudrait essayer de définir d'autres touches pour voir peut etre si c plus un code page pas pareil qui pourrait faire ça ou pas...

---------------
Tout cul tendu mérite son dû

FORUM HardWare.fr

Programmation

C++

recherche algo balistique (bonus : jeux inside ! )

Sujets relatifs
[PHP] Recherche programmeur	recherche script pour bare de defilment ...
Besoin d'aide pour algo procédures/tableau	deux requetes en une! [SQL inside]
A la recherche d'un cobol-mode pour Emacs	[HTML] please, Aidez moi a centrer mes calques... (débutant inside)
Moteur de recherche efficace site intranet	recherche votre avis sur 2 pitites questions projet info
ident failed for user machin?? d'où ça vient (Postgres inside)[résolu]	algo de langage evenementiel HELP
Plus de sujets relatifs à : recherche algo balistique (bonus : jeux inside ! )

Page générée en 0.093 secondes