IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

MATLAB Discussion :

programmer une fonction


Sujet :

MATLAB

  1. #1
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut programmer une fonction
    Bonjour,

    comment puisse je écrire le programme suivant sur Matlab:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    for i=1:n
        w(i)=w(i-1)*i+w_0;
    end
    avec
    en fait j'ai la formule suivante:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    k=T/n;                    % le pas
    w(i)=w(i-1)*i+w_0;
    C_m(i)=J/k*(i*w(i)+w_0+(k*a/J)*(w(i))^2
    et j'essaie d'avoir C_m j'ai testé un programme sur matlab, mais ça marche pas, il m'affiche une erreur avec le w(i)=w(i-1)*i+w_0;

    cordialement.

  2. #2
    Modérateur
    Avatar de le fab
    Homme Profil pro
    Développeur informatique
    Inscrit en
    Mars 2005
    Messages
    1 881
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 48
    Localisation : France, Isère (Rhône Alpes)

    Informations professionnelles :
    Activité : Développeur informatique
    Secteur : Industrie

    Informations forums :
    Inscription : Mars 2005
    Messages : 1 881
    Points : 3 429
    Points
    3 429
    Par défaut
    salut

    Citation Envoyé par Rita94 Voir le message
    mais ça marche pas, il m'affiche une erreur avec le w(i)=w(i-1)*i+w_0;
    quelle erreur ?

    probablement que l'indice 0 n'existe pas (pour i = 1 )

  3. #3
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut
    il m'affiche cet erreur:
    Error in test_a_vitesse_variable (line 29)
        w(i)=w(i-1)*i+w_0;

  4. #4
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Recopie le message d'erreur complet
    Voir la FAQ : Comment exploiter les messages d'erreurs ?
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  5. #5
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut
    voici la fonction que j'ai fait:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    function [time w]=vitesse (w_0,T,n)
    k=T/n;                            % n subdivisions
    w=zeros(1,n);
    t=0 ;                             % temps initial
    Time=zeros(1,n);
    for i=1:n
        Time(i)=t+(i-1)*k;
        w(i)=w(i-1)*i+w_0;
        t=t+k;
    end
    end
    et le message d'erreur en testant le programme avec la fonction de w est:

    Subscript indices must either be real positive integers or logicals.
    
    Error in vitesse (line 8)
        w(i)=w(i-1)*i+w_0;
    
    Error in test_a_vitesse_variable (line 28)
    [Time w]=vitesse(w_0,T,n);

  6. #6
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  7. #7
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut
    désolé, mais j'ai pas comprit.

  8. #8
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Fabien t'avait pourtant mis sur la piste :

    Citation Envoyé par le fab Voir le message
    probablement que l'indice 0 n'existe pas (pour i = 1 )
    Les indices des matrices commencent à 1 avec MATLAB et non pas à 0 comme avec d'autres langage (C, Python…)

    Tu ne peux pas utiliser w(0) dans ton code.
    Or c'est ce que tu fais pour i = 1
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  9. #9
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut
    J'ai essayé de changer les variables, mais ça marche pas.
    en fait, je suis débutant en MATLAB et du coup de sais pas trop comment manipuler les erreurs.
    cdt.

  10. #10
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    w_0 est équivalent à w(0), c'est ça ?
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  11. #11
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut
    oui

  12. #12
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Tu peux donc considérer la première itération de la boucle comme un cas particuliers :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    for i = 1:n
        Time(i) = t+(i-1)*k;
        t = t+k;
        if i==1
            w(i) = w_0;
        else
            w(i) = w(i-1)*i+w_0;
        end    
    end
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  13. #13
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut
    j'ai ce problème qui vient d'une équation mathématique:

    w(i+1)=w(i)*(i+1)+w(0)
    C_m(i)=J*[i*w(i)+k*a/J*(w(i))^2+w(0)]/k

    avec w(0) constante donnée (notant b) condition initiale
    et J,a sont des constantes donnée aussi
    et le k c'est le pas (=T/n)
    et on prend pour c_m(0)=v ( une valeur donnée)

    comment je peux le résoudre via matlab, afin d'avoir les variation de w et C_m en fonction du temps, et C_m en fonction de w.

    le programme que j'ai fait ne marche pas, si vous pouvez m'aider.

  14. #14
    Nouveau membre du Club
    Homme Profil pro
    Technicien optique intrumentale et optronique
    Inscrit en
    Avril 2015
    Messages
    35
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Paris (Île de France)

    Informations professionnelles :
    Activité : Technicien optique intrumentale et optronique
    Secteur : Industrie

    Informations forums :
    Inscription : Avril 2015
    Messages : 35
    Points : 35
    Points
    35
    Par défaut
    Si la solution donnée par Jerome ne marche pas, tu peux essayer de transformer ton code ainsi :

    Pour n subdivisions (ou pas) tu as probablement n+1 valeurs en comptant celle de départ.

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    function [time w]=vitesse (w_0,T,n)
    k=T/n;                            % n subdivisions
    w=zeros(1,n+1);
    t=0 ;                             % temps initial
    Time=zeros(1,n+1);
    for i=1:n+1
        Time(i+1)=t+(i)*k;
        w(i+1)=w(i)*i+w_0;
        t=t+k;
    end
    end
    Je n'ai pas transformé le i multiplicateur en (i+1) à la ligne du calcul de w mais tu pourrais en avoir besoin! Essaye avec et sans.
    Je ne pense pas que tu doives rajouter un +1 à la ligne k=T/(n+1) mais essaye aussi si le code marche mais que tu as des erreurs de résultat avec les deux options précédentes.

    EDIT: Tu auras peut-être une erreur de type dépassement d'indice genre "Index exceed matrix dimension", si c'est le cas réessaye une boucle for i=1:n.
    Et il faudra de toute façon que tu initialise ton vecteur w en tapant avant ou après la boucle for w(1)=w_0.
    Si tout ça ne marche pas essaye aussi de reprendre ton code comme il est et de faire une boucle for i=2:n, mais il se peut que tu doive aussi remettre tes matrices zeros et ta boucle à n+1.

  15. #15
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut
    je vous remercie infiniment, les deus programmes marchent (mais ils donnent pas les mêmes résultats).

    maintenant, si vous permettez, comment je peux avoir le Cm qui est en relation avec w tel que:

    C_m(i)=J/k*(i*w(i)+k*a(w(i))^2+w_0) ??

    est ce que je dois définir une autre fonction dont j'appelle la fonction w?

    cordialement.

  16. #16
    Nouveau membre du Club
    Homme Profil pro
    Technicien optique intrumentale et optronique
    Inscrit en
    Avril 2015
    Messages
    35
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Paris (Île de France)

    Informations professionnelles :
    Activité : Technicien optique intrumentale et optronique
    Secteur : Industrie

    Informations forums :
    Inscription : Avril 2015
    Messages : 35
    Points : 35
    Points
    35
    Par défaut
    D'après moi tu n'as pas besoin de refaire une fonction mais tu peux si tu veux diviser ton programme. Comme dit plus haut, il faut que tu rajoutes un point devant tes opérateurs * / et ^ lorsqu'ils sont utilisés pour mettre en relation deux matrices (ou vecteurs) (l'opérateur ^ mettant en relation une matrice avec elle même) mais attention il est nécessaire que les deux fassent la même taille. Cela fait traité les deux entités valeur correspondante par valeur correspondante.

    Il faudrait que tu tapes :

    C_m(i)=J/k*(i.*w(i)+k*a(w(i)).^2+w_0)

  17. #17
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut Il faudrait que tu tapes :
    Re-bonjour,

    Gaspar-koening, J'ai essayé d'appliquer ce que vous m'avez expliqué. Alors j'ai créé la fonction suivante:
    code de la fonction:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    function [Time Cm]=couple (T,n,w_0,J)
    Cm=zeros(1,n+1);
    t=0 ;                             % temps initial
    Time=zeros(1,n+1);
    for i=1:n  
        Time(i) = t+(i-1)*T/n;
        t = t+i*T/n;
        C_m(i+1)=J/T/n*((i+1).*w(i+1)+T/n*a(w(i+1)).^2+w_0);
    end
    end
    et dans le programme global, j'ai appelé cette fonction:
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    clear all
    close all
    clc
     
    w_0=24000;
    T=0.3*3600;
    n=10000000;
    J=120;
    [Time w]=vitesse(w_0,T,n);
    [Time Cm]=couple(T,n,w_0,J);
    plot(Time,w)
    hold on
    plot(Time,Cm)
    et il m 'affiche les erreurs suivantes:

    Not enough input arguments.
    
    Error in w (line 3)
    w=12000*ones(1,n+1);
    
    Error in couple (line 8)
        C_m(i+1)=J/T/n*((i+1).*w(i+1)+T/n*a(w(i+1)).^2+w_0);
    
    Error in test_vitesse (line 10)
    [Time Cm]=couple(T,n,w_0,J);
    je pense il y a un problème en argument au niveau de la fonction w, j'ai définie w comme suit
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    function [Time w]=w (w_0,T,n)
     
    w=12000*ones(1,n+1);
    t=0 ;                             % temps initial
    Time=zeros(1,n+1);
    for i=1:n+1
        Time(i+1)=t+i*T/n;
        w(i+1)=w(i)*(i+1)*T/n+w_0;    
        t=t+T/n;
    end
    end
    je sais pas comment résoudre c'est problème, que ce que vous pensez?

  18. #18
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    Regarde bien la ligne suivante :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    function [Time w]=w (w_0,T,n)
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

  19. #19
    Futur Membre du Club
    Femme Profil pro
    Étudiant
    Inscrit en
    Septembre 2017
    Messages
    18
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 29
    Localisation : France, Seine et Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Septembre 2017
    Messages : 18
    Points : 7
    Points
    7
    Par défaut
    la vitesse w a comme paramètre d'entré (w_0,T,n), mais je sais pas comment dire à MATLAB que Cm(i+1) dépend de w(i+1) (w_0,T,n)?

  20. #20
    Rédacteur/Modérateur

    Avatar de Jerome Briot
    Homme Profil pro
    Freelance mécatronique - Conseil, conception et formation
    Inscrit en
    Novembre 2006
    Messages
    20 302
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Garonne (Midi Pyrénées)

    Informations professionnelles :
    Activité : Freelance mécatronique - Conseil, conception et formation

    Informations forums :
    Inscription : Novembre 2006
    Messages : 20 302
    Points : 53 160
    Points
    53 160
    Par défaut
    J'ajoute un peu de couleur :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    function [Time w]= w (w_0,T,n)
    Ingénieur indépendant en mécatronique - Conseil, conception et formation
    • Conception mécanique (Autodesk Fusion 360)
    • Impression 3D (Ultimaker)
    • Développement informatique (Python, MATLAB, C)
    • Programmation de microcontrôleur (Microchip PIC, ESP32, Raspberry Pi, Arduino…)

    « J'étais le meilleur ami que le vieux Jim avait au monde. Il fallait choisir. J'ai réfléchi un moment, puis je me suis dit : "Tant pis ! J'irai en enfer" » (Saint Huck)

+ Répondre à la discussion
Cette discussion est résolue.
Page 1 sur 2 12 DernièreDernière

Discussions similaires

  1. programmer une fonction par VB
    Par bahja_83 dans le forum Windows
    Réponses: 1
    Dernier message: 22/04/2008, 00h41
  2. Programmer une fonction réciproque
    Par bichou dans le forum MATLAB
    Réponses: 4
    Dernier message: 20/06/2007, 16h34
  3. programmer une fonction
    Par Maria1505 dans le forum Visual C++
    Réponses: 3
    Dernier message: 30/04/2007, 20h00
  4. Réponses: 10
    Dernier message: 11/01/2007, 22h45
  5. Programmer une fonction joindre_fichier
    Par leCcsympas dans le forum Réseau
    Réponses: 5
    Dernier message: 03/12/2006, 20h51

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo