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 :

Programme pour variogramme


Sujet :

MATLAB

  1. #1
    Futur Membre du Club
    Programme pour variogramme
    Bonjour,

    J'ai trouvé sur internet un programme sensé pouvoir réaliser un variogramme cependant quand je le lance dans l'interface matlab j'ai une erreur que je n'arrive pas à corriger...

    ??? Error: File: Variogram.m Line: 33 Column: 31
    Unbalanced or unexpected parenthesis or bracket.
    Voici le programme en question (j'ai souligné en rouge l'erreur):

    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
    % Computes 2D variogram {directional or omnidirectional}
    % Copyright alghalandis.com | Younes Fadakar Alghalandis
    function [h,g,c] = variogram(x,y,v,nL,dir,td)
    n = numel(x);
    t = zeros(n*(n-1)/2,1);
    r = t;
    k = 1;
    for o=1:n-1                          % computes distances and angles
       for p=o+1:n
          dx = x(p)-x(o);
          dy = y(p)-y(o);
          t(k) = atan2(dy,dx);
          r(k) = sqrt(dx^2+dy^2);
          k = k+1;
       end
    end
    L = max(r)/nL;                       % initialize lag length
    g = zeros(1,nL);                     % initialize gammas with zeros
    c = g;                               %    and count of pairs
    b = 1;
    ta = deg2rad(td);                    % degree to radian
    for a=deg2rad(dir)
       for s=1:nL
          h = s*L;
          th = L/2;
          k = 1;
          q = 0;
          g(b,s) = 0;
          c(b,s) = 0;
          for o=1:n-1                    % computes variogram for different h's
             for p=o+1:n
                if (abs(r(k)-h)<=th) && (abs(t(k)-a)<=ta)
                   g(b,s) = g(b,s)+(v(p)-v(o))^2;
                   q = q+1;
                end
                k = k+1;
             end
          end
          g(b,s) = g(b,s)/(2*q);
          c(b,s) = q;
       end
       b = b+1;
    end
    h = (1:nL)*L;                        % returns lags as h vector


    Merci d'avance !

    Ludovic

  2. #2
    Rédacteur/Modérateur

    Remplace simplement &lt; par <



    C'est un problème de génération du HTML de son site internet (voir les entités de caractère HTML pour plus d'information)
    Mes compétences :
    • conception mécanique 3D (Autodesk Fusion 360)
    • développement informatique (MATLAB, Python, C, VBA)
    • impression 3D (Ultimaker)
    • programmation de microcontrôleur (Microchip PIC et 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)

  3. #3
    Futur Membre du Club
    Merci, le programme tourne bien maintenant...mais le problème c'est que le variogramme n'est pas du tout bon. J'ai juste en sortie un point en plein milieu du graphe...et vu que je ne suis pas l'auteur du programme je vais avoir les plus grandes difficultés à trouver d'où vient l'erreur.

    Si vous arrivez à la trouver ça m'aiderais vraiment beaucoup !

    Merci d'avance

  4. #4
    Futur Membre du Club
    Je viens de résoudre l'erreur, très bête en plus, j'ai juste défini mes données X, Y et Z comme étant les 3 premières lignes et non les 3 premières colonnes -_-.

  5. #5
    Nouveau Candidat au Club
    erreur
    Citation Envoyé par LudOL56 Voir le message
    Je viens de résoudre l'erreur, très bête en plus, j'ai juste défini mes données X, Y et Z comme étant les 3 premières lignes et non les 3 premières colonnes -_-.
    bonsoir j'ai vu ton programme et rectifier l'erreur comme tu l'a demandé mais je n'arrive toujours pas à implementer; je sais pas peut etre je m'y prend mal
    besoin d'aide s'il vous plait. ça m'affiche ceci:

    >> [x,y,v]=variogram(x,y,v,nL,dir,td)
    Undefined function or variable 'nL'.