Aller au contenu

Chaka

Membres
  • Compteur de contenus

    87
  • Inscription

  • Dernière visite

À propos de Chaka

  • Rang
    Posteur
  • Date de naissance 22/12/2000

Informations

  • Classe
    Autre
  • Sexe
    Garçon
  • Pays/Ville
    Chartres/Limoges

Visiteurs récents du profil

Le bloc de visiteurs récents est désactivé et il n’est pas visible pour les autres utilisateurs.

  1. Bonsoir, Je me permets de vous contacter car je peine à faire ces exercices de révision car je n'ai pas pu assister aux cours. Pour le 1), je trouve une matrice qui vaut $ A^{2}=\begin{pmatrix} 1 & 0 & 12\\ 0 & 1 & 0\\ -24 & -6 & 1 \end{pmatrix} $ ce qui ne correspond à rien de connu... Pour le 2), je vérifie que l'application est bien un projecteur, qu'il est de direction (1,-2) (élément invariant)mais je ne trouve pas la base. Merci par avance pour votre aide. Chaka exos revisions.pdf
  2. Bonjour, Je me permets de poster un message ici car je suis un peu perdu dans l'exercice 2 de ce TP (cf. Pièce-Jointe). Pas de difficulté sur l'exercice 1 mais je ne comprends pas trop ce que je dois faire dans l'exercice 2. Si vous avez des méthodes de résolutions / réponses... je suis preneur ! Bonne fin de journée tp2.pdf
  3. Chaka

    Méthode de Heun

    Oui tout à fait. En tout cas, je te remercie car tu es le seul à me répondre régulièrement qui plus est rapidement ! Bonne soirée.
  4. Chaka

    Méthode de Heun

    Bonsoir, Merci pour ton retour. On est censé trouvé que les méthodes d'ordre 1 (Euler) sont les moins précises, celle d'ordre 2 (Heun) un peu plus et celle d'ordre 4 le plus précis (RK4) Je vais chercher sur le web. Bonne soirée
  5. Chaka

    Méthode de Heun

    Bonjour, \(\left\{ \begin{array}{ll} y'(t)=y(t)\\ y(0)=1 \end{array} \right.\) Cependant, ma méthode de Heun me fais du n'importe quoi (du moins son approximation devrait être bien plus proche de la solution théorique). Voyez-vous mon erreur dans le script : def heun(fh,t0,tf,y0,N): y=[y0] t, pas=np.linspace(t0,tf,N+1), (tf-t0)/N #Initialisation k1=[fh(t[0],y[0])] k2=[fh(t[1],y[0]+pas*k1[0])] for i in range(N): k1.append(fh(t[i],y[i])) k2.append(fh(t[i+1],y[i]+pas*k1[i])) y.append(y[i]+pas*(k1[i]+k2[i])/2) retu
  6. Bonsoir! Merci pour ton retour ! En effet j'étais déjà tombé sur ce lien, mais en essayant de le mettre à ma sauce cela n'avait pas fonctionné. Autre question, j'essaye de comparer les approximations de différents algorithmes sur un problème de Cauchy simple : \(\left\{ \begin{array}{ll} y'(t)=y(t)\\ y(0)=1 \end{array} \right.\). Cependant, ma méthode de Heun me fais du n'importe quoi (du moins son approximation devrait être bien plus proche de la solution théorique). Voyez-vous mon erreur dans le script : def heun(fh,t0,tf,y0,N): y=[y0] t, pas=np.linsp
  7. Bonjour, je rencontre un problème sur la fonction fsolve. En effet python me renvoie cette erreur : TypeError: fsolve() missing 1 required positional argument: 'x0' (le problème viendrait donc du x0=0.01) J'ai lu sur le net que cet problème vennait du lambda composé avec le fsolve et qu'il vallait peut être mieux définir la fonction à côté, mais cette solution m'enquiquine un peu. Je vous poste ci-joint mon script si quelqu'un aurait une/des solution.s ... def eulerimplicite(fei,t0,tf,y0,N): y=[y0] t , pas=np.linspace(t0,tf,N+1) , (tf-t0)/N for i in range(N):
  8. On te demande une expression en fonction de x , tu cherches donc R sous la forme R(x)=ax²+bx+c, de même pour C(x) et B(x) (comme l'a dit anylor B(x)=R(x)-C(x)) Là tu nous sors , je cite, " B(x)=397 298 milliers" ok mais cela veut dire que tu as évalué ta fonction en un certain x=??? Peux-tu nous préciser?
  9. Bonjour, Pour le premier graphe : tout m'a l'air bon ! Pour la seconde question, le professeur te demande une fonction (une expression avec un x dedans). Rappel: la recette est l'argent que tu récoltes (donc sans compter le coût de production). Tu connais le prix d'un article et le nombre d'article (les x qui varient). Tu peux donc trouver l'argent que tu as récolté en vendant les x articles. Pour la troisième question, Rappel : le bénéfice se calcule comme ceci : recette - coup de production de l'ensemble des articles. Quatrième question : application directe de la formule donnée en Q
  10. Bonjour, Je viens de trouver mon erreur! dans la ligne z.append(...), il faut remplacer les k1,k2,... par m1,m2,... (j'ai voulu faire un copié collé pour gagner du temps mais finalement cela m'en a fait perdre) Bonne journée ! Chaka
  11. Je viens de finir mon script cependant je rencontre un problème par rapport aux résultats que j'obtiens : cf pièces jointes et scripts. La première est la figure que j'obtiens avec RK4 (convergence vers l'origine), la seconde avec la méthode de Newton (aspect chaotique). (avec les mêmes constantes et les mêmes conditions initiales) Trouvez-vous une explication à ces différences ? Si non, problèmes dans les scripts ? : ## Euler NS-Lorenz def lorenz(x, y, z, s=10, r=28, b=2.667): """ Entrée: x, y, z: point dans R3 auquel nous nous interessons s, r, b: par
  12. Merci beaucoup julesx !! Je viens de lire le dernier post de ce forum et cela à l'air de répondre à ma question. Entre temps j'avais un peu avancé et étais en bonne voie d'après ce que je viens de lire. Je t'en remercie ! Je posterai le code final pour ceux que cela intéresse Bonne soirée
  13. Je confirme que ce que tu as dit est mal rédigé. Grand bien te fasse de rester borné. Bref je ne veux pas tomber dans un débat stérile avec toi car il ne nous apporterai rien à tout les deux si ce n'est que des ennuis. Bien cordialement. Bonne continuation
  14. Lorenz avec la méthode d'Euler : def lorenz(x, y, z, s=10, r=28, b=2.667): """ Entrée: x, y, z: point dans R3 auquel nous nous interessons s, r, b: paramètres de l'attracteur de Lorenz Sortie: dérivées partiels aux points x, y, z liées aux équations de Lorenz """ x_dot = s*(y - x) y_dot = r*x - y - x*z z_dot = x*y - b*z return x_dot, y_dot, z_dot pas = 0.01 num_iter = 10000 # +1 pour prendre les CI x_sol = np.empty(num_iter + 1) y_sol = np.empty(num_iter + 1) z_sol = np.empty(n
  15. Bonsoir. Je te parlais du message que j'ai précédemment cité : "[...] Poser t^(n-1) = u --> (n-1).t^(n-2) dt = du et poser t * exp((-t^2)/25) dt = dv ---> v = -(25/2) *exp((-t^2)/25) [...] u = t^(n-1) et v' = t * exp((-t^2)/25)". Ce qui est écrit ici c'est un changement de variable. Ici tu as juste posé une variable u= une expression en fonction de t et non u_n(t)= une expression de t (qui serait donc bien une fonction de t). Sur ce ... Bonne soirée
×
×
  • Créer...
spam filtering