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 :

Boucle avec matrices dans un script de fonction


Sujet :

MATLAB

  1. #1
    Candidat au Club
    Homme Profil pro
    No
    Inscrit en
    juin 2015
    Messages
    4
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 33
    Localisation : France, Yvelines (Île de France)

    Informations professionnelles :
    Activité : No

    Informations forums :
    Inscription : juin 2015
    Messages : 4
    Points : 3
    Points
    3
    Par défaut Boucle avec matrices dans un script de fonction
    Bonjour,

    Je suis nouvelle sur ce forum. Jusqu'à présent je trouvais la solution à mes problèmes en parcourant les sujets des autres participants mais cette fois-ci, ce n'est pas le cas !

    Mon problème est le suivant. Je dois résoudre un système d'équations non linéaires et je dois le faire plusieurs fois de suite.
    J'ai d'abord créé des matrices contenant les variables utilisées dans la fonction

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    nb=35;
    for i=1:nb
    eval(['A' num2str(i) '=[m(i,1),n(i,1),o(i,1),p(i,1); q(i,1),r(i,1),s(i,1),t(i,1); u(i,1),v(i,1),w(i,1),x(i,1); y(i,1),z(i,1),k(i,1),h(i,1)]']);
    end
    où m, n, ..., h sont des variables (nb,1) provenant d'un fichier de données.

    Je veux faire appel à ces matrices pour résoudre le système à chaque date (1 à 35). Le problème est que je ne peux faire appel aux matrices avec un indice (ça ne fonctionne pas).
    Ma fonction :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
     
    function F=fun(x)
     
    B = [x(1), x(2), x(3), x(4) ; x(5), x(6), x(7), x(8) ; x(9), x(10), x(11), x(12) ; x(13), x(14), x(15), x(16)];
    C = [x(1), x(2), x(3), x(4) ; x(5), x(6), x(7), x(8) ; x(9), x(10), x(11), x(12) ; x(13), x(14), x(15), x(16)];
    D = [x(1), x(2), x(3), x(4) ; x(5), x(6), x(7), x(8) ; x(9), x(10), x(11), x(12) ; x(13), x(14), x(15), x(16)];
     
    F=(B*C*D)- A(indice i de 1 à 35); 
     
    end
    En solution je veux des matrices (4,4) (comme les matrices A(indice i)) pour tous les i allant de 1 à 35.

    Cette fonction fonctionne (je ne demande pas de boucle) :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
     
    function F=fun(x)
     
    global A1;
     
    B = [x(1), x(2), x(3), x(4) ; x(5), x(6), x(7), x(8) ; x(9), x(10), x(11), x(12) ; x(13), x(14), x(15), x(16)];
    C = [x(1), x(2), x(3), x(4) ; x(5), x(6), x(7), x(8) ; x(9), x(10), x(11), x(12) ; x(13), x(14), x(15), x(16)];
    D = [x(1), x(2), x(3), x(4) ; x(5), x(6), x(7), x(8) ; x(9), x(10), x(11), x(12) ; x(13), x(14), x(15), x(16)];
     
    F=(B*C*D)- A1; 
     
    end
    quand je la résouds comme ceci :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
     
    nb=35;
    for i=1:nb
    eval(['OQT' num2str(i) '=[PP(i,1),PT(i,1),PU(i,1),PI(i,1); TP(i,1),TT(i,1),TU(i,1),TI(i,1); UP(i,1),UT(i,1),UU(i,1),UI(i,1); IP(i,1),IT(i,1),IU(i,1),II(i,1)]']);
    end
     
    global A1;
     
    x0 = [25, 0.1, 0.2, 0.3 ; 1.6, 19.2, 2, 1.5 ; 1.5, 2.2, 15, 4 ; 0.15, 0.15, 0.55, 25];
     
     
    % Make a starting guess at the solution
    options=optimset('Display','iter');   % Option to display output
    [x,fval,EXITFLAG] = fsolve(@fun,x0,options);  % Call solver
    Mais impossible de le faire pour les 35 dates en même temps.

    Merci pour votre aide !!!

  2. #2
    Modérateur

    Homme Profil pro
    Ingénieur R&D
    Inscrit en
    août 2014
    Messages
    1 295
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 30
    Localisation : France, Nord (Nord Pas de Calais)

    Informations professionnelles :
    Activité : Ingénieur R&D
    Secteur : High Tech - Électronique et micro-électronique

    Informations forums :
    Inscription : août 2014
    Messages : 1 295
    Points : 2 385
    Points
    2 385
    Par défaut
    Bonjour,

    Je pense que tu peux déjà aller consulter la FAQ sur la création dynamique des variables nommées A1, A2, A3... AN.
    Si vous cherchez des réponses sur ce forum il faudra avant tout expliquer clairement votre problème et exposer la démarche que vous avez entreprise pour le résoudre. Fournissez une base de travail et de réflexion à vos interlocuteurs!

  3. #3
    Candidat au Club
    Homme Profil pro
    No
    Inscrit en
    juin 2015
    Messages
    4
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 33
    Localisation : France, Yvelines (Île de France)

    Informations professionnelles :
    Activité : No

    Informations forums :
    Inscription : juin 2015
    Messages : 4
    Points : 3
    Points
    3
    Par défaut
    Merci pour votre réponse si rapide.

    J'ai consulté cette rubrique mais j'ai l'impression que cela ne résoud pas mon problème.
    On y trouve que l'indexage peut se faire avec des {} mais ça ne fonctionne pas dans ma fonction...

    Je vais quand même retourner consulter cette rubrique !

Discussions similaires

  1. [Bash] Problème avec locate dans un script
    Par kopros2 dans le forum Administration système
    Réponses: 2
    Dernier message: 07/12/2011, 08h06
  2. Boucle avec String dans JApplet.
    Par Antoniom dans le forum Langage
    Réponses: 4
    Dernier message: 12/11/2010, 17h32
  3. Réponses: 2
    Dernier message: 04/12/2009, 17h27
  4. Intégration de la fusion d'images avec "imagecopymerge" dans un script
    Par foxprox dans le forum EDI, CMS, Outils, Scripts et API
    Réponses: 4
    Dernier message: 24/03/2009, 16h52
  5. Problème dans boucle avec matrice et vecteur
    Par lilyla dans le forum MATLAB
    Réponses: 9
    Dernier message: 26/11/2007, 20h45

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