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 :

un petit programme pour faire un calcul de commande par retour d'état


Sujet :

MATLAB

  1. #1
    Membre habitué
    Profil pro
    Inscrit en
    Août 2009
    Messages
    542
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Août 2009
    Messages : 542
    Points : 187
    Points
    187
    Par défaut un petit programme pour faire un calcul de commande par retour d'état
    Bonjour,

    Je n'arrive pas à faire fonctionner mon programme MATLAB qui est pourtant un petit programme tout simple (quelques calculs sur des matrices). Chaque ligne prise indépendamment fonctionne pourtant dans le mode Command.

    Je pense que j'ai du oublié d'entrée de jeu sur mon programme une ligne de commande de base à mettre.

    Voici mon code :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    %----------------------------------------------------------------------------- 
    %  Détermination K et L dans une commande à retour d'état
    %-----------------------------------------------------------------------------
    %
    % Définition des caractéristiques de la représentation d'état du système
     
    clear all;
    F=[1 0.1 0 0;0 1 0 0;1.9378 0 1.9378 1;-1 0 -1 0];
    g=[0.005;0.1;0;0];
     
    % ETAPE 1 : Définition du polynome caractéristique en boucle fermée. Il est donné à
    % partir d'un cahier des charges souhaitée (cf poly).
     
    alphat=[-0.2338,-0.76622,-0.2432,0.4038];
     
    % ETAPE 2 : Calcul du polynome caractéristique en boucle ouverte. Il est calculé
    % à partir de la représentation d'état
     
    pcbo=poly(F);
    at=pcbo(2:5);
    a1=pcbo(2);
    a2=pcbo(3);
    a3=pcbo(4);
     
    % ETAPE 3 : Détermination de la matrice de Toplitz
     
    Atop=[1,0,0,0;a1,1,0,0;a2,a1,1,0;a3,a2,a1,1];
    invAtop=inv(Atop);
     
    % ETAPE 4 : Détermination de gouvernabilité
     
    gv=ctrb(F,g);
    invgv=inv(gv);
     
    % ETAPE 5 : Calcul de K par la formule de Bass-Gura
     
    K=(at-alphat)*invAtop*invgv;
    Voici le message d'erreur que je récupère:

    Error: Missing operator, comma, or semicolon.

    ??? calcul du retour d'état K
    Je vous remercie d'avance pour votre aide

    Ben

  2. #2
    Expert éminent sénior
    Avatar de Caro-Line
    Profil pro
    Inscrit en
    Mars 2007
    Messages
    9 458
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Mars 2007
    Messages : 9 458
    Points : 14 828
    Points
    14 828
    Par défaut
    Le message d'erreur n'est pas complet.
    Peux-tu copier/coller exactement et en entier celui-ci ?
    (il doit normalement t'indiquer la ligne incriminée)

    Pour information ce message doit aussi apparaitre dans l'éditeur lorsque tu passes sur la ligne là où il y a des soulignements rouge ou orange) comme pour les fautes d'orthographe dans Word.

  3. #3
    Membre émérite
    Homme Profil pro
    Inscrit en
    Mai 2008
    Messages
    2 040
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations forums :
    Inscription : Mai 2008
    Messages : 2 040
    Points : 2 841
    Points
    2 841
    Par défaut
    Bonjour.
    Je n'ai pas d'erreur :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    invAtop =
        1.0000         0         0         0
        3.9378    1.0000    0.0000    0.0000
        9.6307    3.9378    1.0000         0
       18.7245    9.6307    3.9378    1.0000
    invgv =
     -174.6051   18.7303  238.3498  312.6640
      274.6051  -13.7303 -625.9098 -863.6778
     -125.3949    6.2697  536.7702  789.3635
       25.3949   -1.2697 -149.2102 -238.3498
    K =
      -66.8989  -16.3326   74.1359   49.3055

  4. #4
    Membre habitué
    Profil pro
    Inscrit en
    Août 2009
    Messages
    542
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Août 2009
    Messages : 542
    Points : 187
    Points
    187
    Par défaut
    Pour le message d'erreur, c'est bien çà.

    Sinon, vu qu'une personne a reussi à faire tourner le programme, c'est que surement je ne dois pas lancer le programme de la bonne façon.

    Pour le lancer, je fais tout simplement "run".

    Y a-t-il des packages supplémentaires à lancer quand on veut travailler avec des matrices ???

    Merci d'avance pour votre aide

    Ben

  5. #5
    Expert éminent sénior
    Avatar de Caro-Line
    Profil pro
    Inscrit en
    Mars 2007
    Messages
    9 458
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Mars 2007
    Messages : 9 458
    Points : 14 828
    Points
    14 828
    Par défaut
    Le truc qu'on ne voit que rarement....
    Ton fichier ne s'appelerait pas "calcul du retour d'état K" par le plus grand des hasards ?

    Si c'est ça trouve un nom de fichier qui soit accepter par MATLAB : alphanumériques + underscores avec au moins 1 lettre en 1ère position et pas d'accents (ben oui MATLAB c'est américain, ils ne connaissent pas les accents donc ne les gèrent pas).

  6. #6
    Membre habitué
    Avatar de mr_langelot
    Profil pro
    Inscrit en
    Août 2003
    Messages
    113
    Détails du profil
    Informations personnelles :
    Localisation : France, Paris (Île de France)

    Informations forums :
    Inscription : Août 2003
    Messages : 113
    Points : 150
    Points
    150
    Par défaut
    Bonjour,

    ça marche chez moi (matlab 6.1) à partir d'un .m...

    Essaie de mettre du code en commentaire jusqu'à voir ou ça pose problème...

  7. #7
    Expert éminent sénior
    Avatar de Caro-Line
    Profil pro
    Inscrit en
    Mars 2007
    Messages
    9 458
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Mars 2007
    Messages : 9 458
    Points : 14 828
    Points
    14 828
    Par défaut
    Oui le code fonctionne chez tout le monde, toutes versions, toutes plates-formes, c'est juste le nom de son fichier (information non fournie) qui n'est pas accepté.

  8. #8
    Membre habitué
    Profil pro
    Inscrit en
    Août 2009
    Messages
    542
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Août 2009
    Messages : 542
    Points : 187
    Points
    187
    Par défaut
    Merci à vous pour votre aide
    J'ai enlevé les accents et les espaces dans le nom de mon fichier et çà fonctionne!
    Ben

Discussions similaires

  1. Réponses: 1
    Dernier message: 07/09/2006, 10h03
  2. Réponses: 46
    Dernier message: 04/08/2006, 00h18
  3. Petit probléme pour faire un message d'erreur
    Par mygwel dans le forum AWT/Swing
    Réponses: 5
    Dernier message: 16/05/2006, 17h32
  4. utiliser données texte pour faire des calculs
    Par sarah67 dans le forum Access
    Réponses: 20
    Dernier message: 06/02/2006, 14h09
  5. aide petit programme pour débutant
    Par kartp0rqx dans le forum C
    Réponses: 16
    Dernier message: 14/10/2005, 19h31

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