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

  1. #1
    Nouveau Candidat au Club
    Homme Profil pro
    Étudiant
    Inscrit en
    octobre 2018
    Messages
    2
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 25
    Localisation : France, Val de Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant
    Secteur : Enseignement

    Informations forums :
    Inscription : octobre 2018
    Messages : 2
    Points : 1
    Points
    1

    Par défaut Résolution d'un système d'équations non linéaires 4 équations à 4 inconnus

    Bonsoir tout le monde
    J'ai un problème en Matlab j'ai essayé de résoudre un système d'équations de 4 équations à 4 inconnus avec "Solve ", mais les résultats sont illogiques pouvez-vous m'aider s'il vous plaît c'est URGENT et voila mon code tel que M c'est une matrice des coordonnés des points n*3:

    Nom : 12.PNG
Affichages : 62
Taille : 15,4 Ko

    les résultats sont parfois en complexe!!

    Nom : 22.PNG
Affichages : 51
Taille : 36,8 Ko

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

    Informations forums :
    Inscription : mai 2008
    Messages : 1 910
    Points : 2 631
    Points
    2 631

    Par défaut

    Bonjour,

    Qu’as-tu pris pour X ?
    Utilise plutôt # pour mettre ton code.

  3. #3
    Nouveau Candidat au Club
    Homme Profil pro
    Étudiant
    Inscrit en
    octobre 2018
    Messages
    2
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 25
    Localisation : France, Val de Marne (Île de France)

    Informations professionnelles :
    Activité : Étudiant
    Secteur : Enseignement

    Informations forums :
    Inscription : octobre 2018
    Messages : 2
    Points : 1
    Points
    1

    Par défaut

    On faite je vais vous expliquer un peu le problème c'est de détecter un cylindre dans un nuage de points donné par la matrice X, et pour constituer un cylindre de 3 points arbitraires du nuage je suis obligé de résoudre ces équations non linéaires mais le problème que mon code prend beaucoup de temps!!

    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
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    % Detection d'un cylindre
    clear all;
    close all;
    % Charger les points
    [X1]= stlread('cylinder.stl');
    X=unique(X1,'rows');
    %Initialisation
    gamma = 0.01;
    L=[];
    Lmeilleur=[];
    for it=1:100
        i=1;
        % pick up 3 random points
        M= X(floor(rand(3,1)*(size(X,1)))+1,:);
        %Identifaction of plan=>three points (xi, yi, zi), which are not collinear, are necessary: a*xi + b*yi + c*zi +d
        V1 = transpose(M(1,:)- M(2,:));
        V2 = transpose(M(1,:)- M(3,:));
        N = cross(V1,V2);
        d = (N(1)*M(1,1)+N(2)*M(1,2)+N(3)*M(1,3));
        % the radius (r) and the center c = (x0, y0, z0) are obtained by solving the following equations.
        syms x0 y0 z0 r;
        eq=[((M(1,1)-x0)^2)+((M(1,2)-y0)^2)+((M(1,3)-z0)^2)==r^2;
            ((M(2,1)-x0)^2)+((M(2,2)-y0)^2)+((M(2,3)-z0)^2)==r^2;
            ((M(3,1)-x0)^2)+((M(3,2)-y0)^2)+((M(3,3)-z0)^2==r^2);
            (N(1)*x0)+(N(2)*y0)+(N(3)*z0)-d==0];
        sol=solve(eq,[x0,y0,z0,r],'IgnoreAnalyticConstraints',true,'PrincipalValue', true);
        x0 = vpa(sol.x0); y0 = vpa(sol.y0);z0 = vpa(sol.z0);r = vpa(sol.r);
        C=[x0 y0 z0];
        % the minimum distance (dminj ) between each point of D and the cylinder axis
        for ind= 1:length(X)
            P = transpose(C-X(ind,:));
            D = norm(cross(N,P))/norm(N);
            if D < r+gamma
                L(i,:)=X(ind,:);
                i=i+1;
            end
        end
        if length(L)>length(Lmeilleur)
            Lmeilleur=L;
        end
    end
    % Plot
    figure; hold on;
    plot3(X(:,1),X(:,2),X(:,3),'or');
    plot3(Mmeilleur(:,1),Mmeilleur(:,2),Mmeilleur(:,3),'og');
    xlabel('x');ylabel('y');zlabel('z');
    title('RANSAC results for cylinder estimation');
    legend('Outliers','Inliers');
    axis equal tight;

Discussions similaires

  1. Résolution d'un système d'équations non linéaires+
    Par Volimps15 dans le forum MATLAB
    Réponses: 0
    Dernier message: 04/11/2015, 01h36
  2. [Débutant] Résolution d'un système d'équations non-linéaires
    Par nicollivier dans le forum MATLAB
    Réponses: 9
    Dernier message: 09/12/2009, 17h32
  3. Réponses: 1
    Dernier message: 22/05/2008, 13h56
  4. Système équations non linéaires, résolution ?
    Par smercier2 dans le forum Mathématiques
    Réponses: 23
    Dernier message: 23/10/2007, 04h00
  5. Réponses: 1
    Dernier message: 14/02/2007, 11h12

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