exercice octave problème
Répondre à la discussion
Affichage des résultats 1 à 3 sur 3

exercice octave problème



  1. #1
    invite4d86fd80

    exercice octave problème


    ------

    salut

    voila je dois modifier un programme octave

    mais le probléme est que je ne comprends pas et je n'y arrive pas , j'ai testé mais rien ne marche

    voila les question qui me pose probleme

    3. Modifier le programme energy.m pour qu’il imprime également l’énergie totale de la masse m
    pour les valeurs initiales de la position données par la question 1.

    4. Modifier le programme energy.m pour qu’il utilise la méthode d’Euler à la place de celle
    d’Euler-Cromer. Refaire la question 1 avec ce nouveau programme et comparer les résultats
    avec la méthode d’Euler-Cromer. Conclusion.

    5. Modifier le programme energy.m pour qu’il présente la vitesse et la position en fonction du
    temps. Exécuter votre programme pour les valeurs de la question 1.

    6. Ecrire un programme qui permet de calculer la période d’oscillation de la particule. Représenter
    graphiquement cette période pour les valeurs initiales de la question 1. Utiliser un graphique
    semi-logarithmique.

    voila le programme en question

    Code:
    % This matlab program computes the kinetic and potential energies
    % of a particle under the action of a conservative force
    clear all;
    clf;
    % Initialisation of alpha and beta
    alpha =1.0;
    beta=0.1;
    % Input the initial position
    x=input('Enter the initial value of the displacement in meters:');
    % set up the initial velocity, mass and time step
    v=0; % Initial velocity in meters/s.
    mass = 1.0 % Mass of particle
    step = 0.02 %time step in seconds
    Kenergy(1) =0.5 * mass * v^2; % kinitic energy
    Uenergy(1) = -0.5 *alpha * x^2 + beta * x^4; % potential energy
    Tenergy(1) = Kenergy(1) + Uenergy(1); % total energy
    Nstep = 601
    % start the calculation
    % start preparing the graphics
    xlabel('Time (s)'); % x-axis label
    ylabel('Enrgy (J)'); % y-axis label
    title('Kinitic energy (-) and potential energy (+) and total (*)')
    hold on;
    for istep = 1: Nstep % make the number of steps necessary
    Kenergy (istep) =0.5 * mass * v^2; % kinitic energy
    Uenergy (istep)= -0.5 *alpha * x^2 + beta * x^4; % potential energy
    t (istep) = (istep -1) * step; % time
    Force = alpha * x - 4 * beta * x^3; % compute the force on the particle
    accel = Force / mass; % compute the accel. on the particle
    % compute the velocity and position using Euler-Cromer algo.
    v = v + accel * step; x = x + v * step;
    endfor
    Energylimit1 =max(Kenergy) + 0.1
    Energylimit2 =min(Uenergy) - 0.1
    %Energylimit =int32 (Energylimit+0.5)
    axis([0, Nstep * step, Energylimit2(1), Energylimit1(1)]); % Set axis limits
    plot(t,Kenergy,'r-'); % plot the Kinetic energy
    plot(t,Uenergy,'b+'); % plot the potential energy
    merci d'avance pour toute aide

    -----
    Dernière modification par yoda1234 ; 29/11/2009 à 04h53.

  2. #2
    inviteba120af5

    Re : exercice octave probléme

    bonjour j'ai un programme à faire merci de m'aider

    Soit un tableau T de taille N où chaque cellule ci avec 1 ≤ i ≤ N contient un nombre. Trier T par ordre croissant consiste à réordonner les valeurs des cellules de manière à ce que :
    √ 1≤i ≤ N et 1 ≤j≤ N : ci ≤ cj

    Pour passer d’un tableau non trié à un tableau trié par ordre croissant, il existe plusieurs méthodes de tri. L’objectif de ce projet est de programmer deux de ces méthodes pour pouvoir ensuite comparer leur efficacité.
    Exercice 3 : Fonction de Tri par insertion
    Le tri par insertion consiste à classer les deux premiers éléments du tableau. Une fois que les deux premiers sont ordonnés, on prend l’élément qui suit et on le classe à son tour dans ce qui a déjà été classé. Pour chaque élément i du tableau, on sait que les éléments de 1 à i-1 sont déjà classés. On va chercher la position j parmi les i-1 premiers éléments de manière à ce que T[i] < T[j]. On insère alors T[i] à la position j.
    Soit le tableau de départ suivant
    On commence par classer les 2 premiers éléments
    Puis on prend l’élément à la position i=3 (T[3] = 7), on cherche la position j dans le sous tableau déjà trié ([5 ; 8]) : j=2, x=T[3]=7, on décale tous les éléments d’une case entre j et i-1 soit ici entre 2 et 2. Donc seul 8 est décalé d’une case. Enfin on insère x à la position j.
    On prend ensuite l’élément à la position i=4 (T[4] = 1), on cherche la position j dans le sous tableau déjà trié ([5 ; 7 ; 8]) : j=1, x=T[4]=1, on décale tous les éléments d’une case entre j et i-1 soit ici entre 1 et 3. Donc 5, 7 et 8 sont décalés d’une case. Enfin on insère x à la position j.
    Enfin, on prend l’élément à la position i=5 (T[5] = 2), on cherche la position j dans le sous tableau déjà trié ([1 ; 5 ; 7 ; 8]) : j=2, x=T[5] = 2, on décale tous les éléments d’une case entre j et i-1 soit ici entre 2 et 4. Donc 5, 7 et 8 sont décalés d’une case. Enfin on insère x à la position j.
    Programmez une fonction Octave qui réalise le tri par insertion d’un tableau.

  3. #3
    whoami

    Re : exercice octave probléme

    Bonjour,

    Comme nous ne sommes pas sur www.OnFaitTesDevoirs.com, montre-nous ce que tu as fait, explique-nous ce qui te bloque.

    Si tu n'as pas compris les explications de l'algorithme incluses dans l'énoncé, la bonne démarche est d'en chercher d'autres, et pas de vouloir que quelqu'un fasse ton travail.

Discussions similaires

  1. Matlab/Octave
    Par invite3dfb4bc4 dans le forum Mathématiques du supérieur
    Réponses: 2
    Dernier message: 20/10/2015, 11h26
  2. Tutoriel Octave.
    Par antoineg dans le forum Mathématiques du supérieur
    Réponses: 0
    Dernier message: 14/03/2010, 17h33
  3. Octave
    Par inviteec8c9286 dans le forum Physique
    Réponses: 1
    Dernier message: 29/04/2009, 20h56
  4. [OCTAVE] Problème 'quad' et argument
    Par Youry dans le forum Logiciel - Software - Open Source
    Réponses: 1
    Dernier message: 14/01/2009, 17h00
  5. Octave
    Par julien_4230 dans le forum Logiciel - Software - Open Source
    Réponses: 2
    Dernier message: 26/03/2007, 20h23