Xcas

Un article de Wikipédia, l'encyclopédie libre.
Sauter à la navigation Sauter à la recherche
Xcas
Description de l'image Xcaslogo.png.
Description de l'image xcas.png.
Informations
Développé par Bernard Parisse
Dernière version 1.5.0 (décembre 2018)
Dépôt sourceforge.net/p/xcas/code/HEAD/treeVoir et modifier les données sur Wikidata
Écrit en C++Voir et modifier les données sur Wikidata
Système d'exploitation Linux, Microsoft Windows et macOSVoir et modifier les données sur Wikidata
Environnement Microsoft Windows, Linux, Mac OS X, FreeBSD
Type Calcul formel
Licence Licence publique générale GNUVoir et modifier les données sur Wikidata
Site web site officiel

Xcas (se prononce /ikskas/) est un logiciel de calcul formel libre pour Windows, Mac OS X, Linux et Unix[1].

Xcas est une interface de Giac, une bibliothèque C++ de calcul formel libre (licence GPL). Giac dispose d'un mode de compatibilité avec les logiciels Maple et Matlab[2] et MuPAD et Yacas et Qcas et WordMat et CPMP-Tools et ExpressionsinBar et les calculatrices TI-89, TI-92, Voyage 200 et TI-Nspire[3]. On peut donc utiliser Giac/Xcas aussi bien comme un logiciel gratuit compatible avec Maple, pour développer des algorithmes de calcul formel, ou l'utiliser dans un autre logiciel…

Xcas est intégré au complément CmathOOoCAS qui permet d'effectuer du calcul formel dans le tableur Calc et dans le traitement de textes Writer de la suite bureautique OpenOffice.org. Il est développé à l'université Joseph-Fourier de Grenoble[4].

Capacités[modifier | modifier le code]

Voici un bref aperçu de ce que Xcas est capable de faire :

  • calcul formel ;
  • géométrie dans le plan[5] ;
  • géométrie dans l'espace[6] ;
  • tableur ;
  • statistiques ;
  • graphiques ;
  • programmation.

Quelques exemples[modifier | modifier le code]

Simulation de la chute d'un objet[modifier | modifier le code]

On utilise les fonctions préprogrammées suivantes :

  • int() qui renvoie la primitive d'une fonction.
  • unapply() qui permet d'évaluer une expression selon un paramètre.
  • solve() qui permet de résoudre une équation.
  • plot() qui permet de créer et afficher un graphique.
  • string() qui permet de transformer un objet (exemple : nombre entier) en une chaine de caractères.
  • diff() Différenciation de la fonction
  • desolve() résoudre l'équation différentielle
  • factor() Polynôme d'affacturage
  • nPr() calculer les permutations
  • nCr() calculer des combinaisons
  • sqrt() racine carrée
  • cross([1,2,3],[4,3,2]) = [-5,10,-5] calculer le produit croisé de deux vecteurs
  • mean([3,4,2]) = 3 calculer la moyenne
  • stddev([3,4,2]) = sqrt(2/3) calculer l'écart type
  • variance([3,4,2]) = 2/3 calculer la variance
  • det([1,2],[3,4]) = -2 calculer le déterminant d'une matrice
  • extrema(-2*cos(x)-cos(x)^2,x) = [0],[pi] calculer les extrema locaux d'une fonction
  • line(x=1) = la ligne verticale x=1 tracez une ligne verticale dans le système de coordonnées

Trouver plus de commandes ici : http://www-fourier.ujf-grenoble.fr/~parisse/giac/cascmd_en.pdf

Premier modèle : Une chute sans frottement[modifier | modifier le code]

Première méthode : fonction qui retourne un nombre (la vitesse à l'impact)[modifier | modifier le code]

On sait que la gravité terrestre nous fournit une accélération de 9,81 m s−2, donc il suffit de l'intégrer deux fois pour avoir la position.

À la ligne 6, on remplit L avec deux solutions : une négative et une positive. On ne prendra que la solution positive, qui se trouve être la deuxième (on utilisera l'indice 1, c'est à dire L[1] (L'indice 0 étant la première solution)).

Code utilisé :

1 vitesse_chute_1_1(h):={
2   local a,v,x,L,temps_chute;
3   a(t):=9.81;
4   v:=unapply(int(a(t),t),t);
5   x:=unapply(int(v(t),t),t);
6   L:=solve(x(t)=h,t);
7   temps_chute:=L[1]
8   retourne v(temps_chute);
9 }:;

On a donc ici l'appel de la fonction :

vitesse_chute_1_1(9)

Qui renvoie :

13.2883407542

Qui est bien sûr la vitesse (en m s−1) que l'objet atteint au niveau du sol après une chute sans frottement d'une hauteur de 9 m.

Deuxième méthode : fonction qui retourne une phrase et un graphique[modifier | modifier le code]

On garde la même structure que la fonction vitesse_chute_1_1, mais on ajoute une phrase et un graphique. On met un moins devant la vitesse et la position, car l'objet tombe vers le bas, qui est considéré négatif (l'altitude étant positive vers le haut).

Code utilisé :

 1 vitesse_chute_1_2(h):={
 2   local a,v,x,L,temps_chute;
 3   a(t):=9.81;
 4   v:=unapply(int(a(t),t),t);
 5   x:=unapply(int(v(t),t),t);
 6   L:=solve(x(t)=h,t);
 7   temps_chute:=L[1];
 8   title="Altitude et vitesse en fonction du temps";
 9   plot(-x(t),t,0,temps_chute,couleur=2+line_width_6);
10   plot(-v(t),t,0,temps_chute,couleur=1+line_width_6);
11   retourne "Chute de "+string(h)+" mètres : Vitesse au niveau du sol après "+string(temps_chute)+" secondes de chute : "+string(v(temps_chute))+" m.s^(-1) = "+string(v(temps_chute)*3.6)+" km.h^(-1)";
12 }:;

On a donc l'appel de la fonction :

vitesse_chute_1_2(9)

Qui renvoie :

Renvoi

Ainsi qu'un graphique :

Vitesse (rouge) et position (vert) d'un objet faisant une chute de 9 mètres sans frottement (dans une gravité terrestre).

Deuxième modèle : Une chute dans l'air avec frottement du type proportionnel (à la vitesse) : [modifier | modifier le code]

Cette fois ci, on prend en compte la force de frottement  : les fonctions donnant la vitesse et la position ont été calculées manuellement, grâce au principe fondamental de la dynamique.

On trouve les fonctions (du temps) suivantes :

    • Avec
    • Avec

NB : On prend en supposant que et que la vitesse terminale (vitesse maximale lors de la chute) de cette masses est de 200 km h−1.

De plus, on ne met pas de moins devant la vitesse et la position, car le signe a déjà été pris en compte dans la réalisation des formules.

Code utilisé :

 1 vitesse_chute_2(h):={
 2   local g,v,x,L,temps_chute;
 3   g:=9.81:;
 4   k:=12.4;
 5   m:=70:;
 6   v:=unapply((m*g)/k*(exp((-t)*(k/m))-1),t);
 7   x:=unapply(g*(m/k)^2*(1-e^((-t)*(k/m)))-(m*g*t)/k,t):;
 8   L:=solve(x(t)=-h,t)
 9   temps_chute:=L[0];
10   title="Chute de 70 kg de 9 mètre : altitude (vert) et vitesse (rouge) en fonction du temps : (épais = sans frottement ; fin = avec frottements)";
11   plot(v(t),t,0,temps_chute,couleur=1);
12   plot(x(t),t,0,temps_chute,couleur=2);
13   retourne 0;
14 }:;

On a donc l'appel de la fonction :

vitesse_chute_2(9)

Qui renvoie :

0

(Qui est juste présent pour vérifier que la fonction a été lue jusqu'au bout.)

Ainsi qu'un graphique :

Vitesse (rouge) et position (vert) d'un objet faisant une chute de 9 mètres avec frottement proportionnel à la vitesse (dans une gravité terrestre).

Superposition des deux modèles[modifier | modifier le code]

Vitesse (rouge) et position (vert) d'un objet de 70 kg faisant une chute de 9 mètres avec (trais fins) et sans (trais épais) frottement proportionnel à la vitesse, ceci dans une gravité terrestre.

On constate que les deux graphiques (avec et sans frottement) sont effectivement différents ; l'accélération n'est pas constante dans le cas du frottement : la vitesse (rouge) tend à se stabiliser (pour devenir horizontale).

Notes et références[modifier | modifier le code]

  1. « Xcas Calcul Formel Lycee | Intégral | Variable (Mathématiques) », sur Scribd (consulté le 20 octobre 2019)
  2. (en-US) « Systeme | » (consulté le 8 janvier 2020)
  3. « Xcas - Mathematical software - swMATH », sur swmath.org (consulté le 21 décembre 2019)
  4. « About: Xcas », sur dbpedia.org (consulté le 21 décembre 2019)
  5. « Download Limit Exceeded », sur citeseerx.ist.psu.edu (consulté le 21 décembre 2019)
  6. « Download Limit Exceeded », sur citeseerx.ist.psu.edu (consulté le 21 décembre 2019)

Liens externes[modifier | modifier le code]