FR 3060814 A1 20180622 FR 1663067 A 20161221 1 PROCEDE DE DETERMINATION DES PARAMETRES PHYSIQUES D'UN SYSTEME, PROCEDE DE REGULATION D'UN SYSTEME A PARTIR DESDITS PARAMETRES ET DISPOSITIF METTANT EN OEUVRE LESDITS PROCEDES DOMAINE TECHNIQUE DE L'INVENTION Le domaine technique de l'invention est celui de la mesure des paramètres physiques d'un système. La présente invention concerne un procédé de détermination des paramètres physiques d'un système et en particulier un procédé de détermination des paramètres physiques d'un système à l'aide d'une pluralité de points d'exploration. ARRIERE-PLAN TECHNOLOGIQUE DE L'INVENTION Dans le cadre du développement des systèmes de gestion d'énergie évolutifs, il est souvent nécessaire d'optimiser la consommation en énergie d'installations, par exemple des bâtiments. Cette optimisation se fait en général sur la base d'un modèle physique non linéaire qu'il est nécessaire de calibrer avant de pouvoir être utilisé à des fins de gestion d'énergie. Afin de calibrer ces modèles, plusieurs méthodes de type méta-heuristique existent, comme les algorithmes de génétique ou bien encore les colonies de fourmis. Ces méthodes présentent cependant l'inconvénient d'être sensibles à la présence de minima locaux qui fausse la calibration du modèle. Il en résulte une moins bonne gestion de l'énergie, le modèle physique utilisé pour ladite gestion étant trop éloigné de la réalité de l'installation que l'on cherche à réguler. La qualité de la gestion du système d'énergie s'en trouve dégradée. Il existe donc un besoin d'une méthode permettant de calibrer le modèle physique sans être sensible à la présence de minima locaux dans la fonction d'optimisation. 3060814 2 RESUME DE L'INVENTION L'invention offre une solution aux problèmes évoqués précédemment, en permettant de déterminer les paramètres physiques d'un système par le biais d'un procédé qui est peu sensible aux minima locaux. Ceci est notamment rendu 5 possible par la mise en place, lors de la mise en oeuvre de l'invention, d'une pluralité de points d'exploration et d'une pluralité de groupes, chaque point d'exploration pouvant appartenir à un ou plusieurs groupes. Pour cela, un premier aspect de l'invention concerne un procédé de détermination des paramètres physiques d'un système en fonction d'un critère d'arrêt 10 comprenant : - une étape d'acquisition d'une première pluralité de grandeurs physiques à un premier instant k et d'une deuxième pluralité de grandeurs physique à un deuxième instant k+ 1, la valeur des grandeurs physiques de la deuxième pluralité à l'instant k+ 1 pouvant être déduite de la valeur des 15 grandeurs physiques de la première pluralité à l'instant k par l'intermédiaire d'une équation non linéaire dépendante desdits paramètres physiques ; - une étape de création d'une pluralité de points d'explorations ci, chaque point d'exploration étant associé à un premier vecteur vi, dit vecteur vitesse, un deuxième vecteur xi, dit vecteur position et à un scalaire pi, dit 20 éloignement, chaque coordonnée du vecteur position xi étant associé à un paramètre physique à déterminer et l'éloignement étant fonction du vecteur position, de la valeur de la première pluralité de grandeurs physiques à l'instant k, de la valeur de la deuxième pluralité de grandeurs physiques à l'instant k+ 1 et de l'équation non linéaire ; 25 - une étape de création d'une pluralité de groupes de points d'exploration, chaque groupe comprenant un nombre de points d'exploration supérieur à 1, un point d'exploration ci pouvant appartenir à plusieurs groupes ; - pour chaque point d'exploration ci, une étape d'initialisation aléatoire des coordonnées des vecteurs vitesse vi et position pi ; 3060814 3 - au moins une itération d'une étape de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci de la pluralité de points d'exploration comprenant, pour chaque itération t : ^ une sous-étape de détermination du point d'exploration parmi la pluralité 5 de point d'exploration dont l'éloignement à l'itération t est le moins élevé, dit point leader global ; ^ pour chaque groupe, une sous-étape de détermination du point d'exploration dont l'éloignement à l'itération t est le moins élevé, dit point d'exploration leader du groupe ; 10 ^ pour chaque point ci de la pluralité de point, une sous-étape de détermination des coordonnées des vecteurs dudit point à l'itération t + 1, lesdites coordonnées étant déterminées en fonction des coordonnées des vecteurs dudit point ci, du ou des points leaders du ou des groupes auxquels appartient ledit point ci et du point leader global à 15 l'itération t. Le procédé comprend en outre, lorsque le critère d'arrêt est rempli à l'issue de l'étape de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration, une étape de détermination des paramètres du système à partir des coordonnées du vecteur position du point leader global, le critère d'arrêt étant 20 vérifié à chaque itération de l'étape de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci. Ainsi, à l'issue du procédé, on obtient les paramètres physiques du modèle. De plus, lors de la mise en oeuvre du procédé, les points d'exploration sont non seulement influencés par le leader global, mais également par le ou les leaders de 25 leurs groupes, ce qui permet d'éviter une trop forte influence d'un leader global dont la position s'orienterait vers un minimum local. Cela rend le procédé selon l'invention peu sensible à la présence de minima locaux dans le modèle dudit système. On entend par « critère rempli » le fait que la condition fixée par ledit critère est bien satisfaite. Par exemple, si le critère concerne le fait qu'une grandeur donnée est au-dessus d'un seuil, alors le critère sera rempli lorsque ladite valeur sera au-dessus dudit seuil. On entend par « le critère d'arrêt est 3060814 4 vérifié » le fait de vérifier si la condition fixée par ledit critère est satisfaite ou non afin de déterminer si ledit critère est rempli ou non. Outre les caractéristiques qui viennent d'être évoquées dans le paragraphe précédent, le procédé selon un premier aspect de l'invention peut présenter une 5 ou plusieurs caractéristiques complémentaires parmi les suivantes, considérées individuellement ou selon toutes les combinaisons techniquement possibles. Avantageusement, le critère d'arrêt est rempli lorsqu'au moins une des conditions suivantes est vérifiée : - l'éloignement pi du point leader global est inférieur à une première valeur 10 seuil prédéterminée ; - le nombre d'itérations de l'étape de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci est supérieur à une deuxième valeur seuil prédéterminée ; - la variation de l'éloignement pi du point leader entre J itérations 15 consécutives de l'étape de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci est inférieure à une troisième valeur seuil prédéterminée, J étant un entier positif fixé en fonction de la précision souhaitée. Ainsi, il est possible de modifier le critère d'arrêt en fonction du système considéré 20 et/ou de la précision désirée. Avantageusement, l'éloignement pi du point ci à l'itération t est égal à la norme de la différence entre la valeur de la deuxième pluralité de grandeurs physiques mesurée à l'instant k + 1 et la valeur de la deuxième pluralité de grandeurs physiques calculée à partir de l'équation non linéaire et en fonction de la première 25 pluralité de grandeur physique mesurée à l'instant k, les paramètres de ladite équation étant égaux aux coordonnées du vecteur position xi du point ci à l'itération t. Avantageusement, la détermination des coordonnées des vecteurs d'un point à l'itération t + 1 comprend : 3060814 5 - pour chaque groupe auquel appartient ledit point, une étape de réévaluation de la vitesse à l'itération t, cette dernière étant donnée par : vi(t) = wxvi(t) + ai(t)xR(c1)x (49(0 - xi(t)) + R(c2)x(xi(t) - xi(t)) où w est un paramètre compris entre 0 et 1, R(c) est un paramètre aléatoire 5 à partir d'une distribution uniforme en fonction de c et c = c1 = c2 = 0.25, x19(t) est la position du point leader du groupe relatif au point d'exploration considéré à l'itération t, xi(t) est la position du point leader global à l'itération t, et bxi(t) est donnée par bxi(t) = xi(e) avec pi(e) = min(pi(t"))Vt" E [0, t], la réévaluation de la vitesse dans un groupe donné 10 tenant compte des réévaluations faites auparavant dans les autres groupes auxquels appartient ledit point ; - une étape de détermination de la vitesse à l'itération t+ 1, cette dernière étant donnée par vi(t + 1) = vi(t) + R(c3)x(bxi(t) - xi(t)) avec c3 = 0.5 ; - une étape de détermination de la position à l'itération t + 1, cette dernière 15 étant donnée par xi(t + 1) = xi(t) + vi(t + 1) ; - une étape de calcul de l'éloignement à partir de la position nouvellement déterminée. Avantageusement, le nombre de points d'exploration I est égal à 40 + VT\1, N étant le nombre de paramètres du système à déterminer. 20 Avantageusement, le nombre de groupes M est égal à -2i - V I étant le nombre de points d'exploration et N le nombre de paramètres du système à déterminer. Avantageusement, l'étape d'acquisition comprend l'acquisition d'une première pluralité de grandeurs physiques aux instants k à k + R et d'une deuxième pluralité de grandeurs physiques aux instants k + 1 à k+(R + 1), la valeur des grandeurs 25 physiques de la deuxième pluralité à l'instant k + (r + 1) pouvant être déduite de la valeur des grandeurs physiques de la première pluralité à l'instant k + r avec r E [0, R] par l'intermédiaire d'une équation non linéaire dépendante desdits paramètres physiques, et en ce que l'éloignement pi du point ci à l'itération t est égal à la somme des normes des différences entre la valeur de la deuxième 3060814 6 pluralité de grandeurs physiques mesurée à l'instant k + (r + 1) et la valeur de la deuxième pluralité de grandeurs physiques calculée à partir de l'équation non linéaire et en fonction de la première pluralité de grandeur physique mesurée à l'instant k + r pour tout r E [0, R], les paramètres de ladite équation étant égaux 5 aux coordonnées de la position xi du point ci à l'itération t. Ainsi, l'éloignement pi n'est plus évalué en fonction de l'accord entre le modèle et la mesure pour une seule mesure, mais pour une pluralité de mesures de la première pluralité grandeurs physiques et la deuxième pluralité de grandeurs physiques. 10 Avantageusement, un nombre de points d'exploration minimum Imin et maximum Imax ainsi qu'un nombre d'itérations maximal Tmax sont définis et, lorsque le nombre d'itérations T de l'étape de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration est égal au nombre itérations maximal Tmax, le procédé comprend une étape de suppression et de réévaluation 15 comportant : - une sous-étape de suppression du point d'exploration ayant l'éloignement le plus élevé ; - une sous-étape de réorganisation aléatoire des groupes ; - au moins une nouvelle itération de l'étape de détermination de nouvelles 20 valeurs des vecteurs de chaque point d'exploration ci de la pluralité de points d'exploration, une nouvelle étape de suppression et de réévaluation étant mise en oeuvre lorsque le nombre d'itérations T de l'étape de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration est égale au nombre d'itérations maximum Tmax ; 25 le procédé comprenant, lorsque le nombre de points d'exploration I est égal au nombre de points d'exploration minimum Imin, une étape d'ajout et de réévaluation comportant : - une sous-étape d'insertion d'un point d'exploration dont les coordonnées position correspondent à une région non encore explorée ; 30 - une sous-étape de réorganisation aléatoire des groupes ; 3060814 7 au moins une nouvelle itération de l'étape de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration ci de la pluralité de points d'exploration, une nouvelle étape d'ajout et de réévaluation étant mise en oeuvre lorsque le nombre d'itérations T de l'étape de détermination 5 de nouvelles valeurs des vecteurs de chaque point d'exploration est égale au nombre d'itérations maximum Tm' ; le procédé comprenant, lorsque le nombre de points d'exploration I est égal au nombre de points d'exploration maximum Im', une nouvelle étape de suppression et de réévaluation. 10 Ainsi, lorsque le procédé selon un premier aspect de l'invention ne permet pas de converger de manière satisfaisante, il est possible de modifier le nombre des points d'exploration en supprimant certains d'entre eux et/ou en insérant de nouveaux points d'exploration. Un deuxième aspect de l'invention concerne un procédé de régulation thermique 15 d'un bâtiment à partir d'un modèle dudit bâtiment, les paramètres dudit modèle étant obtenus à l'aide d'un procédé selon un premier aspect de l'invention. Un troisième aspect de l'invention concerne un dispositif comprenant des moyens de calculs et des moyens d'acquisition, lesdits moyens étant configurés pour mettre en oeuvre un procédé selon un premier aspect de l'invention ou un 20 deuxième aspect de l'invention. Un quatrième aspect de l'invention concerne un produit programme d'ordinateur comprenant des instructions qui conduisent le dispositif selon un troisième aspect de l'invention à exécuter les étapes du procédé selon un premier aspect de l'invention ou selon un second aspect de l'invention. 25 Un cinquième aspect de l'invention concerne un support lisible par ordinateur, sur lequel est enregistré le programme d'ordinateur selon un quatrième aspect de l'invention. L'invention et ses différentes applications seront mieux comprises à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. 3060814 8 BREVE DESCRIPTION DES FIGURES Les figures sont présentées à titre indicatif et nullement limitatif de l'invention. - Les figures 1A et 1B montrent un ordinogramme d'un premier mode de réalisation d'un procédé selon un premier aspect de l'invention. 5 - La figure 2 montre un ordinogramme d'un deuxième mode de réalisation d'un procédé selon un premier aspect de l'invention. DESCRIPTION DETAILLEE D'AU MOINS UN MODE DE REALISATION DE L'INVENTION Sauf précision contraire, un même élément apparaissant sur des figures 10 différentes présente une référence unique. Un premier mode de réalisation selon un premier aspect de l'invention illustré aux figures 1A et 1B concerne un procédé de détermination des paramètres physiques d'un système en fonction d'un critère d'arrêt. Pour cela, le procédé comprend une étape 101 d'acquisition d'une première pluralité de grandeurs physiques à un 15 premier instant k et d'une deuxième pluralité de grandeurs physiques à un deuxième instant k + 1, la valeur des grandeurs physiques de la deuxième pluralité à l'instant k + 1 pouvant être déduite de la valeur des grandeurs physiques de la première pluralité à l'instant k par l'intermédiaire d'une équation non linéaire dépendante desdits paramètres physiques. 20 Le procédé comprend en outre une étape 102 de création d'une pluralité de points d'exploration ci, chaque point d'exploration étant associé à un premier vecteur v;, dit vecteur vitesse, un deuxième vecteur xi, dit vecteur position et un scalaire pi, dit éloignement, chaque coordonnée du vecteur position xi étant associée à un paramètre physique à déterminer et l'éloignement pi étant fonction du vecteur 25 position xi, de la valeur de la première pluralité de grandeurs physiques à l'instant k, de la valeur de la deuxième pluralité de grandeurs physiques à l'instant k + 1 et de l'équation non linéaire. Ainsi, la connaissance de la position xi d'un point ci permet également de connaître l'éloignement pi dudit point. 3060814 9 Dans un mode de réalisation, l'éloignement pi est défini comme la norme de la différence entre la valeur de la deuxième pluralité de grandeurs physiques mesurée à l'instant k + 1 et la valeur de la deuxième pluralité de grandeurs physiques calculée à partir de l'équation non linéaire et en fonction de la première 5 pluralité de grandeurs physiques mesurée à l'instant k. Par exemple, en prenant F(k) la première pluralité de valeurs à l'instant k, G(k + 1) la deuxième pluralité de valeur à l'instant k + 1, ces deux valeurs étant liée par G(k + 1) = AxF(k), l'éloignement pourra être défini comme : pi = II G(k + 1) - NiF(k) II avec Ni = f(xi) 10 Dans cet exemple, les paramètres à déterminer sont donc les coefficients de la matrice A. De plus, chaque coefficient de la matrice A'i est donné par une coordonnée du vecteur xi. Autrement dit, plus l'éloignement est faible, plus les coefficients de la matrice A'i sont proches des coefficients recherchés de la matrice A. Or, les coefficients de la matrice A'i sont égaux aux coefficients du 15 vecteur position xi du point ci considéré. Autrement dit, plus l'éloignement pi est faible, plus les coordonnées du vecteur position xi du point ci considéré sont proches des paramètres physiques que l'on cherche à déterminer. Le procédé comprend également une étape 103 de création d'une pluralité de groupes de points d'exploration, chaque groupe comprenant un nombre de points 20 d'exploration supérieur à 1, un point d'exploration pouvant appartenir à plusieurs groupes. Puis, pour chaque point d'exploration ci, le procédé comprend une étape 104 d'initialisation aléatoire des coordonnées du vecteur vitesse vi et du vecteur position xi. Lorsque l'on dispose de valeurs approximatives des paramètres à déterminer, les valeurs aléatoires pourront être bornées en fonction 25 de ces valeurs approximatives et de l'incertitude associée à ces valeurs. Le procédé comprend ensuite au moins une itération d'une étape 105 de détermination de nouvelles valeurs des vecteurs vitesse vi et position xi et de l'éloignement pi de chaque point d'exploration ci de la pluralité de points d'exploration. Cependant, de manière générale, le procédé comprend une pluralité 3060814 10 d'itérations car il est rare qu'une itération suffise à déterminer les paramètres du système. Plus particulièrement, comme illustré à la figure 1B, cette étape de détermination comprend pour chaque itération t une sous-étape 1051 de détermination du point 5 d'exploration parmi la pluralité de points d'exploration dont l'éloignement à l'itération t est le moins élevé, dit point leader global. Cette sous-étape 1051 permet donc de connaître le point d'exploration ayant l'éloignement le plus faible ou, autrement dit, le point d'exploration dont les coordonnées du vecteur position se rapprochent le plus de la valeur des paramètres physiques que l'on cherche à 10 déterminer. L'étape de détermination comprend également, pour chaque itération t et pour chaque groupe, une sous-étape 1052 de détermination du point d'exploration dont l'éloignement à l'itération t est le moins élevé du groupe considéré, dit point d'exploration leader du groupe. 15 Enfin, elle comprend pour chaque itération t et pour chaque point ci de la pluralité de point, une sous-étape 1053 de détermination des coordonnées des vecteurs vitesse vi et position xi et de l'éloignement pi dudit point ci à l'itération t + 1, lesdites coordonnées et ledit éloignement étant déterminés en fonction des coordonnées des vecteurs dudit point ci, du ou des points leaders du ou des 20 groupes auxquels appartient ledit point et du point leader global à l'itération t. Ainsi, lors de ces itérations, les points d'exploration vont avoir tendance à adopter la vitesse des leaders de leurs groupes. De plus, les points d'exploration vont également avoir tendance à adopter la vitesse du leader global. Or, les leaders, qu'ils soient leader global ou leader de groupe, sont les points ci possédant 25 l'éloignement pi le moins élevé. Autrement dit, ils correspondent aux points d'exploration dont la valeur des coordonnées du vecteur position xi se rapprochent le plus de la valeur des paramètres physiques que l'on cherche à déterminer. L'ensemble des points d'exploration va donc avoir tendance à converger vers la position conduisant à l'éloignement pi le moins élevé et donc vers les valeurs des 3060814 11 paramètres physiques que l'on cherche à déterminer. Il est important de noter que, lors de cette convergence, la notion de groupe est importante. En effet, les points d'exploration sont non seulement influencés par le leader global, mais également par le ou les leaders de leurs groupes, ce qui permet d'éviter une trop 5 forte influence d'un leader global dont la position s'orienterait vers un minimum local. Cela rend le procédé selon l'invention peu sensible à la présence de minima locaux dans le modèle dudit système. On comprend donc que durant ces itérations, les points vont voir, au moins pour une partie d'entre eux, leur éloignement pi tendre vers une valeur faible ou nulle. 10 Enfin, le procédé comprend, lorsque le critère d'arrêt est rempli, une étape de détermination des paramètres du système à partir des coordonnées du vecteur position du point leader global, le critère d'arrêt étant vérifié à chaque itération de étape 105 de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci. On entend par « critère rempli » 15 le fait que la condition fixée par ledit critère est bien satisfaite. Par exemple, si le critère concerne le fait qu'une grandeur donnée est au-dessus d'un seuil, alors le critère sera rempli lorsque ladite valeur sera au-dessus dudit seuil. On entend par « le critère d'arrêt est vérifié » le fait de vérifier si la condition fixée par ledit critère est satisfaite ou non afin de déterminer si ledit critère est rempli ou non. 20 Dans un mode de réalisation le critère d'arrêt est rempli lorsque l'éloignement pi du point leader global est inférieur à une première valeur seuil prédéterminée. Cette valeur seuil peut être choisie en fonction de la précision souhaitée concernant les paramètres physiques du système à déterminer. Dans un mode de réalisation alternatif ou complémentaire, le critère d'arrêt est 25 rempli lorsque le nombre d'itérations de l'étape 105 de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci est supérieur à une deuxième valeur seuil prédéterminée. Dans un mode de réalisation alternatif ou complémentaire, le critère d'arrêt est rempli lorsque la variation de l'éloignement pi du point leader entre J itérations 30 consécutives de l'étape 105 de détermination de nouvelles valeurs des vecteurs et 3060814 12 de l'éloignement pi de chaque point d'exploration ci est inférieure à une troisième valeur seuil prédéterminée, J étant un entier positif fixé en fonction de la précision souhaitée. Il est important de noter ici que les critères précédents peuvent être combinés 5 entre eux. Par exemple, le critère sera rempli lorsqu'au moins une des conditions énumérées ci-dessus est satisfaite. On pourra également envisager de ne considérer le critère rempli que lorsque deux conditions au moins sont satisfaites. Dans un mode de réalisation, la sous-étape 1053 de détermination des coordonnées des vecteurs et de l'éloignement d'un point à l'itération t + 1 10 comprend pour chaque groupe, une étape de réévaluation de la vitesse à l'itération t, cette dernière étant donnée par : vi(t) = wxvi(t) + ai(t)xR(c1)x (xlg(t) - xi(t)) + R(c2)x(xi(t) - xi(t)) où w est un paramètre compris entre 0 et 1, R(c) est un paramètre aléatoire à partir d'une distribution en fonction de c et c = cl = c2 = 0.25, xig(t) est la position 15 du point leader du groupe relatif au point d'exploration considéré à l'itération t, xl(t) est la position du point leader global à l'itération t, et bxi(t) est donnée par : bxi(t) = xi(V) avec pi(e) = min(pi(t")) Vt" E [0, t]. Autrement dit, bxi(t) est la position pour laquelle l'éloignement a été le plus faible en considérant l'ensemble des itérations déjà effectuées ou encore la meilleure 20 position connue à cet instant pour le point considéré. De plus, la réévaluation dans un groupe donné tient compte des réévaluations faites auparavant dans les autres groupes auxquels appartient le point ci considéré. Ainsi, lorsqu'un point d'exploration appartient à plusieurs groupes, la valeur vi(t) du membre de droite tient compte des réévaluations effectuées en 25 considérant les groupes déjà évalués. Par exemple, si le point ce appartient à un premier groupe g1 et à un deuxième groupe g2, la vitesse du point sera d'abord évaluée en fonction du premier groupe avec : vie-eevaim = wxve (t) + ae (t) x R( ci) x (xlegl(t) - xe(t)) + R(c2)x(xl(t) - xe(t)) 3060814 13 Puis ensuite réévaluer en fonction du deuxième groupe avec : Ve reeval2 et) = wxvi-eeval ) OEe(t) xR(c1)x (xieg2(t) - xe(t)) + R(c2)x(xi(t) - xe(t)) Ainsi, un point appartenant à un nombre de groupe NG verra sa vitesse à 5 l'itération t réévalué NG fois. Une fois la vitesse vi réévaluée pour l'ensemble des points d'exploration, la sous-étape 1051 de détermination des coordonnées des vecteurs d'un point à l'itération t + 1 comprend, pour chaque point d'exploration, une étape de détermination de la vitesse vi à l'itération t + 1, cette dernière étant donnée par : 10 vi(t + 1) = Vfeeval(t) R(C3) X (bXi(t) - Xi (t)) où, comme expliqué précédemment, bxi(t) est la position pour laquelle l'éloignement a été le plus faible en considérant l'ensemble des itérations déjà effectuées ou encore la meilleure position connue à cet instant pour le point considéré. Il est ensuite possible d'en déduire le vecteur position xi de chaque 15 point ci lors d'une étape de détermination de la position xi à l'itération t + 1, cette dernière étant donnée par : xi (t + 1) = xi (t) + vi(t + 1) Enfin, une fois la position xi de chaque point d'exploration ci connue, la détermination des coordonnées des vecteurs d'un point à l'itération t + 1 20 comprend une étape de calcul de l'éloignement pi à partir de la position xi nouvellement déterminée. Ainsi, à l'issue de ces étapes, le vecteur vitesse vi, le vecteur position xi et l'éloignement pi à l'itération t + 1 ont été déterminés. Dans un mode de réalisation, le nombre de points d'exploration est fonction du 25 nombre de paramètres à déterminer. Par exemple, s'il y a N paramètres à déterminer alors le nombre de points d'exploration I sera donné par la relation I = 40 + V. Dans un mode de réalisation, le nombre de groupes est fonction du nombre N de 3060814 14 paramètres à déterminer. Par exemple, s'il y a N paramètre à déterminer, alors le nombre de de groupes M sera donné par la relation M = -21 + VT\1. Dans un mode de réalisation, l'étape 101 d'acquisition comprend l'acquisition d'une première pluralité de grandeurs physiques aux instants k à k + R et d'une 5 deuxième pluralité de grandeurs physiques aux instants k + 1 à k+(R + 1), la valeur des grandeurs physiques de la deuxième pluralité à l'instant k + (r + 1) pouvant être déduite de la valeur des grandeurs physiques de la première pluralité à l'instant k + r avec r E [0, R] par l'intermédiaire d'une équation non linéaire dépendante desdits paramètres physiques. 10 Dans ce mode de réalisation, l'éloignement pi du point ci à l'itération t est égal à la somme des normes des différences entre la valeur de la deuxième pluralité de grandeurs physiques mesurée à l'instant k + (r + 1) et la valeur de la deuxième pluralité de grandeurs physiques calculée à partir de l'équation non linéaire et en fonction de la première pluralité de grandeur physique mesurée à l'instant k + r 15 pour tout r E [0, R], les paramètres de ladite équation étant égaux aux coordonnées de la position xi du point ci à l'itération t. Par exemple, en prenant F(k) la première pluralité de valeurs à l'instant k, G(k + 1) la deuxième pluralité de valeurs à l'instant k + 1, ces deux valeurs étant liées par G(k + 1) = AxF(k), l'éloignement pi peut être défini comme : 20 pi = + r + 1) - A'iF(k + 011 avec Ni = «xi). r=1 De manière générale, toute définition mathématique d'une distance entre deux points quelconques peut être utilisée pour définir l'éloignement. Ainsi, l'éloignement pi n'est plus évalué en fonction de l'accord entre le modèle et la mesure pour une seule mesure, mais pour une pluralité de mesures de la 25 première pluralité grandeurs physiques et la deuxième pluralité de grandeurs physiques. On notera que si l'instant k est choisi comme l'instant 0, alors la formule précédente peut être mise sous la forme suivante : 3060814 15 pi = + 1) - NiF(r) II avec Ni = «xi). r Dans certaines circonstances, lorsque l'on constate que le procédé ne permet pas de converger de manière satisfaisante, il peut être avantageux de modifier le nombre des points d'exploration en supprimant certains d'entre eux et/ou en 5 insérant de nouveaux points d'exploration. Pour cela, dans un mode de réalisation illustré à la figure 2, un nombre de points d'exploration minimum Imin et maximum Imax ainsi qu'un nombre d'itérations maximal Tmax sont définis. De plus, lorsque le nombre d'itérations T de l'étape 105 de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration égal à Tmax, le procédé selon 10 un premier aspect de l'invention comprend une étape107 de suppression et de réévaluation. Cette étape 107 suppression et de réévaluation s'effectue en deux temps. Elle comprend tout d'abord une sous-étape 1071 de suppression du point d'exploration ayant l'éloignement pi le plus élevée. Elle comprend ensuite une sous-étape 1072 15 de réorganisation aléatoire des groupes. A l'issue de cette réorganisation, elle comprend au moins une nouvelle d'itérations de l'étape 105 de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration de la pluralité de points d'exploration. Lorsque le nombre d'itérations T de l'étape 105 de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration 20 égal à Tmax, le procédé selon un premier aspect de l'invention comprend une nouvelle sous-étape 1071 de suppression du point d'exploration ayant l'éloignement le plus élevé. Le procédé comprend également, lorsque le nombre de points d'exploration I est égal au nombre de points d'exploration minimal Imin, une étape 108 d'ajout et de 25 réévaluation. Cette étape 108 d'ajout et de réévaluation s'effectue en trois temps. Tout d'abord elle comprend une sous-étape 1081 d'insertion d'un point d'exploration dont les coordonnées position correspondent à une région non encore explorée. Elle comprend ensuite une sous-étape 1082 de réorganisation 3060814 16 aléatoire des groupes. Enfin, elle comprend au moins une nouvelle itération de l'étape 105 de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration de la pluralité de points d'exploration. Lorsque le nombre d'itérations T de l'étape 105 de détermination de nouvelles valeurs des vecteurs 5 de chaque point d'exploration égal à Tm', le procédé selon un premier aspect de l'invention comprend une nouvelle sous-étape 1071 de suppression du point d'exploration ayant l'éloignement le plus élevé. Le procédé comprend également, lorsque le nombre de points d'exploration I est égal au nombre de points d'exploration maximal Im', une nouvelle étape 106 de suppression et de 10 réévaluation. Le compteur d'itération (autrement dit le nombre d'itérations) de l'étape 105 de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration est remis à zéro lors de chaque étape sous-étape 1071 de suppression ou sous-étape 1081 d'insertion. Ainsi, le procédé alterne entre une étape 107 de suppression et de réévaluation et une étape 108 d'ajout et de 15 réévaluation. Comme précédemment le procédé comprend, lorsque le critère d'arrêt est rempli à l'issue de l'étape 105 de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration, une étape de détermination des paramètres du système à partir des coordonnées du vecteur position dudit point leader global, 20 l'éloignement étant comparé à la valeur seuil à chaque itération de étape 105 de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci. Dans un mode de réalisation, le système est un bâtiment et les paramètres que l'on cherche à déterminer sont relatifs au comportement en température dudit 25 bâtiment. Le procédé comprend alors une étape 101 d'acquisition de la température T en différents points du bâtiment ainsi que les vecteurs énergies U (le chauffage, l'ensoleillement etc..) qui affectent le bâtiment à un instant k ainsi que la température T en différents points du bâtiment à un instant k + 1. Ces différentes variables sont en outre liées entre elles par 30 l'équation non linéaire suivante : T (k + 1) = AxT (k) + B xU (k) 3060814 17 où A et B sont deux matrices permettant de définir le comportement thermique du bâtiment. Pour un bâtiment donné, cette équation peut être développée sous la forme suivante : U1(k) U2 (k) 5 [Ti(k + 1)1 [A1 A21 [ Ti(k) + [B1 B2 B3 B4 B5 B6 X U3 (k) [Tm (k + 1)] [A3 A4] [Tm (k)] [B7 B8 B9 B10 B11 B12J U4 (k) U5(k) U6(k)_ Les coefficients de la matrice A et de la matrice B, lorsqu'ils sont connus, permettent de calculer la réponse thermique d'un bâtiment. Ils correspondent donc aux paramètres physiques que l'on cherche à déterminer au moyen d'un procédé selon un premier aspect de l'invention. 10 Le procédé comprend en outre une étape 102 de création d'une pluralité de points d'explorations ci, chaque point d'exploration étant associé à un premier vecteur vi, dit vecteur vitesse, un deuxième vecteur xi, dit vecteur position et un scalaire pi, dit éloignement, chaque coordonnée du vecteur position xi étant associée à un paramètre physique à déterminer et l'éloignement étant fonction du vecteur 15 position, de la valeur de la première pluralité de grandeurs physiques à l'instant k, de la valeur de la deuxième pluralité de grandeurs physiques à l'instant k + 1 et de l'équation non linéaire. Plus particulièrement, l'éloignement est égal à pi = IIT(k + 1) - P"dei(k + 1)11 avec Tnwciei(k + 1) = iT(k) + 13' iU(k) où A' i et 13'i sont les matrices définissant 20 le comportement thermique dont les coefficients sont donnés par les coordonnées du vecteur position xi. Autrement dit, dans l'exemple de l'équation explicitée ci-dessus, le vecteur position xi du point d'exploration ci a pour coordonnées : xi = A'i2 A'i3 A'i4 13'i2 13'i3 13'i413'i513'i6 e 13'i8 1310 B11 B121. Le procédé comprend également une étape 103 de détermination d'une pluralité 3060814 18 de groupes de points d'exploration, chaque groupe comprenant un nombre de points d'exploration supérieur à 1, un point d'exploration pouvant appartenir à plusieurs groupes. Puis, pour chaque point d'exploration ci, le procédé comprend une étape 104 5 d'initialisation aléatoire des coordonnées des vecteurs vitesses et des vecteurs positions. Dans un mode de réalisation, le nombre de points d'exploration et le nombre de groupe sont fonction du nombre N de paramètres à déterminer. Dans l'exemple illustré plus haut, le nombre N de paramètres est égal à 16, ce qui correspond à 10 un nombre de points d'exploration I égal à 40 + -\/16 = 44 et un nombre de groupes M égal à4-. + -\/16 = 26. Lorsque cela est nécessaire , les valeurs obtenues sont arrondies. Le procédé comprend ensuite une pluralité d'itérations d'une étape 105 de détermination de nouvelles valeurs des vecteurs et de l'éloignement de chaque 15 point d'exploration ci de la pluralité de points d'exploration. Cette étape 105 de détermination comprend pour chaque itération t : une sous-étape 1051 de détermination du point d'exploration parmi la pluralité de points d'exploration dont l'éloignement à l'itération t est le moins élevé, dit point leader global ; 20 pour chaque groupe, une sous-étape 1052 de détermination du point d'exploration dont l'éloignement à l'itération t est le moins élevé au sein du groupe considéré, dit point d'exploration leader du groupe ; pour chaque point de la pluralité de point, une sous-étape 1053 de détermination des coordonnées des vecteurs et de l'éloignement dudit 25 point à l'itération t + 1, lesdites coordonnées et ledit éloignement étant déterminés en fonction des coordonnées des vecteurs dudit point, du ou des points leaders du ou des groupes auxquels appartient ledit point et du point leader global à l'itération t. De plus le procédé comprend, lorsque le critère d'arrêt est rempli (par exemple 3060814 19 lorsque l'éloignement pi du point leader global est inférieur à une valeur seuil) à l'issue de l'étape 105 de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration, une étape 106 de détermination des paramètres du système à partir des coordonnées du vecteur position dudit point leader global. 5 Autrement dit, si le point cs est le point leader lorsqu'il est mis fin au procédé, alors les coefficients des matrices A et B (c'est-à-dire les paramètres du système à déterminer) sont données par les coordonnées du vecteur position : xs = [AsA2 3 A's4 131 Blo BI BI] Dans un mode de réalisation, l'étape d'acquisition comprend l'acquisition de la 10 température T en différents points du bâtiment ainsi que les vecteurs énergies U qui affectent le bâtiment aux instants k à k + R et la température T en différents points du bâtiment aux instants k + 1 à k+(R + 1). De plus, l'éloignement pi du point ci à l'itération t est égal à la sommes des normes des différences entre la température T en différents points du bâtiment à un instant k + r + 1 et la valeur 15 de la température T en différents points du bâtiment calculée à partir de l'équation non linéaire et en fonction de la température T en différents points du bâtiment ainsi que les vecteurs énergies U qui affectent le bâtiment l'instant k + r pour tout r E [0, R], les paramètres de ladite équation étant égaux aux coordonnées de la position xi du point ci à l'itération t. Autrement dit, l'éloignement est défini par la 20 relation suivante : pi =11T(k + r + 1) - + r + 1)11 r=1 Un tel procédé peut donc être utilisé pour déterminer les paramètres physiques concernant le comportement thermique d'un bâtiment. . Comme précédemment, on notera que si l'instant k est choisi comme l'instant 0, la formule précédente 25 devient : R pi =11T(r + 1) - Tmodel + 1)11 r=1 3060814 20 Les paramètres ainsi déterminés pourront ensuite être utilisés dans un procédé de régulation du circuit de chauffage par exemple. Pour cela, un deuxième aspect de l'invention concerne un procédé de régulation thermique d'un bâtiment à partir d'un modèle dudit bâtiment, les paramètres dudit 5 modèle étant obtenus à l'aide d'un procédé selon un premier aspect de l'invention. Le modèle, une fois ces paramètres déterminés à l'aide d'un procédé de détermination selon un premier aspect de l'invention, pourra par exemple calculer précisément le moment qu'on doit allumer le chauffage et la quantité d'énergie qu'il faut mettre pour réduire la consommation du chauffage en tenant en compte 10 l'évolution de la météo. Afin de mettre en oeuvre un procédé selon un premier aspect ou un deuxième aspect de l'invention, un troisième aspect de l'invention concerne un dispositif de détermination des paramètres d'un système. Pour cela, le dispositif comporte des moyens d'acquisition d'une première pluralité de grandeurs physiques à un 15 premier instant k et d'une deuxième pluralité de grandeurs physique à un deuxième instant k + 1. Il comporte également des moyens de calculs reliés auxdits moyens d'acquisition. De manière générale, la nature des moyens d'acquisition et leur disposition au sein et/ou à proximité du système dépendent des paramètres que l'on cherche à 20 déterminer. Par exemple, dans le cas de la détermination des paramètres thermiques d'un bâtiment, ces moyens peuvent comprendre des capteurs de température disposés à l'intérieur et/ou à l'extérieur du bâtiment. La consommation de chauffage pourra être obtenue par l'intermédiaire d'un compteur électrique. Ledit dispositif pourra également comporter des moyens de connexion 25 à un réseau informatique, par exemple Internet, afin d'obtenir des informations météorologiques concernant la situation actuelle et futur de la zone dans laquelle se situe le bâtiment (par exemple des informations sur l'ensoleillement). Les moyens de calcul peuvent prendre la forme d'un processeur et d'une mémoire, d'un FPGA (pour Field-Programmable Gate Array en anglais, Circuit 30 Logique Programmable en français) ou d'une carte de type ASIC (pour 3060814 21 Application-Specific Integrated Circuit en anglais). Ces moyens de calcul sont configurés pour mettre en oeuvre un procédé selon un premier aspect de l'invention ou un second aspect de l'invention. Plus particulièrement, ces moyens de calcul sont configurés pour créer une 5 pluralité de points d'explorations ci, chaque point d'exploration étant associé à un premier vecteur vi, dit vecteur vitesse, un deuxième vecteur xi, dit vecteur position et à un scalaire pi, dit éloignement, chaque coordonnée du vecteur position xi étant associé à un paramètre physique à déterminer et l'éloignement étant fonction du vecteur position, de la valeur de la première pluralité de grandeurs 10 physiques à l'instant k, de la valeur de la deuxième pluralité de grandeurs physiques à l'instant k + 1 et de l'équation non linéaire. Les moyens de calcul sont également configurés pour créer une pluralité de groupes de points d'exploration, chaque groupe comprenant un nombre de points d'exploration supérieur à 1, un point d'exploration ci pouvant appartenir à plusieurs 15 groupes et, pour chaque point d'exploration ci, initialiser aléatoirement les coordonnées des vecteurs vitesse vi et position pi. Les moyens de calcul sont en outre configurés pour effectuer au moins une itération d'une étape 105 de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci de la pluralité de points 20 d'exploration comprenant, pour chaque itération t : une sous-étape 1051 de détermination du point d'exploration parmi la pluralité de point d'exploration dont l'éloignement à l'itération t est le moins élevé, dit point leader global ; pour chaque groupe, une sous-étape 1052 de détermination du point 25 d'exploration dont l'éloignement à l'itération t est le moins élevé, dit point d'exploration leader du groupe ; pour chaque point ci de la pluralité de point, une sous-étape 1053 de détermination des coordonnées des vecteurs dudit point à l'itération t + 1, lesdites coordonnées étant déterminées en fonction des 30 coordonnées des vecteurs dudit point ci, du ou des points leaders du ou 3060814 22 des groupes auxquels appartient ledit point ci et du point leader global à l'itération t. Enfin, les moyens de calculs sont configurés, pour lorsque le critère d'arrêt est rempli à l'issue de l'étape de détermination de nouvelles valeurs des vecteurs de 5 chaque point d'exploration, déterminer les paramètres du système à partir des coordonnées du vecteur position dudit point leader global, le critère d'arrêt étant vérifié, par le moyen de calcul, à chaque itération de l'étape 105 de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci. 10 REVENDICATIONS1. Procédé de détermination des paramètres physiques d'un système en fonction d'un critère d'arrêt caractérisé en ce qu'il comprend : une étape (101) d'acquisition d'une première pluralité de grandeurs physiques à un premier instant k et d'une deuxième pluralité de grandeurs physique à un deuxième instant k + 1, la valeur des grandeurs physiques de la deuxième pluralité à l'instant k + 1 pouvant être déduite de la valeur des grandeurs physiques de la première pluralité à l'instant k par l'intermédiaire d'une équation non linéaire dépendante desdits paramètres physiques ; une étape (102) de création d'une pluralité de points d'explorations ci, chaque point d'exploration étant associé à un premier vecteur vi, dit vecteur vitesse, un deuxième vecteur xi, dit vecteur position et à un scalaire pi, dit éloignement, chaque coordonnée du vecteur position xi étant associé à un paramètre physique à déterminer et l'éloignement étant fonction du vecteur position, de la valeur de la première pluralité de grandeurs physiques à l'instant k, de la valeur de la deuxième pluralité de grandeurs physiques à l'instant k + 1 et de l'équation non linéaire ; une étape (103) de création d'une pluralité de groupes de points d'exploration, chaque groupe comprenant un nombre de points d'exploration supérieur à 1, un point d'exploration ci pouvant appartenir à plusieurs groupes ; pour chaque point d'exploration ci, une étape (104) d'initialisation aléatoire des coordonnées des vecteurs vitesse vi et position pi ; au moins une itération d'une étape (105) de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci de la pluralité de points d'exploration comprenant, pour chaque itération t : ^ une sous-étape (1051) de détermination du point d'exploration 3060814 24 parmi la pluralité de point d'exploration dont l'éloignement à l'itération t est le moins élevé, dit point leader global ; ^ pour chaque groupe, une sous-étape (1052) de détermination du point d'exploration dont l'éloignement à l'itération t est le moins 5 élevé, dit point d'exploration leader du groupe ; ^ pour chaque point ci de la pluralité de point, une sous-étape (1053) de détermination des coordonnées des vecteurs dudit point à l'itération t + 1, lesdites coordonnées étant déterminées en fonction des coordonnées des vecteurs dudit 10 point ci, du ou des points leaders du ou des groupes auxquels appartient ledit point ci et du point leader global à l'itération t ; le procédé comprenant, lorsque le critère d'arrêt est rempli à l'issue de l'étape (105) de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration, une étape (106) de détermination des 15 paramètres du système à partir des coordonnées du vecteur position dudit point leader global, le critère d'arrêt étant vérifié à chaque itération de l'étape (105) de détermination de nouvelles valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci. 2. Procédé selon la revendication précédente caractérisé en ce que le critère 20 d'arrêt est rempli lorsqu'au moins une des conditions suivantes est vérifiée : l'éloignement pi du point leader global est inférieur à une première valeur seuil prédéterminée ; le nombre d'itérations de l'étape (105) de détermination de nouvelles 25 valeurs des vecteurs et de l'éloignement pi de chaque point d'exploration ci est supérieur à une deuxième valeur seuil prédéterminée ; la variation de l'éloignement pi du point leader entre J itérations consécutives de l'étape 105 de détermination de nouvelles valeurs des 30 vecteurs et de l'éloignement pi de chaque point d'exploration ci est inférieure à une troisième valeur seuil prédéterminée, J étant un entier 3060814 25 positif fixé en fonction de la précision souhaitée. 3. Procédé selon l'une des revendications précédentes caractérisé en ce que l'éloignement pi du point ci à l'itération t est égal à la norme de la différence entre la valeur de la deuxième pluralité de grandeurs physiques mesurée à 5 l'instant k + 1 et la valeur de la deuxième pluralité de grandeurs physiques calculée à partir de l'équation non linéaire et en fonction de la première pluralité de grandeur physique mesurée à l'instant k, les paramètres de ladite équation étant égaux aux coordonnées du vecteur position xi du point ci à l'itération t. 10 4. Procédé selon l'une des revendications précédentes caractérisé en ce que la détermination des coordonnées des vecteurs d'un point à l'itération t + 1 comprend : pour chaque groupe auquel appartient ledit point, une étape de réévaluation de la vitesse à l'itération t, cette dernière étant donnée par : 15 vi(t) = wxvi(t) + ai(t)xR(c1)x (49(0 - xi(t)) + R(c2)x(xi(t) - xi(t)) où w est un paramètre compris entre 0 et 1, R(c) est un paramètre aléatoire à partir d'une distribution uniforme en fonction de c et c = cl = c2 = 0.25, xf9(t) est la position du point leader du groupe relatif au point d'exploration considéré à l'itération t, xi(t) est la position du point leader 20 global à l'itération t, et bxi(t) est donnée par bxi(t) = xi(e) avec Mt') = min(pi(t"))Vt" E [0, t], la réévaluation de la vitesse dans un groupe donné tenant compte des réévaluations faites auparavant dans les autres groupes auxquels appartient ledit point ; une étape de détermination de la vitesse à l'itération t+ 1, cette dernière 25 étant donnée par vi(t + 1) = vi(t) + R(c3)x(bxi(t) - xi(t)) avec c3 = 0.5 ; une étape de détermination de la position à l'itération t + 1, cette dernière étant donnée par xi(t + 1) = xi(t) + vi(t + 1) ; une étape de calcul de l'éloignement à partir de la position nouvellement 30 déterminée. 3060814 26 5. Procédé selon l'une des revendications précédentes caractérisé en ce que le nombre de points d'exploration I est égal à 40 +/171, N étant le nombre de paramètres du système à déterminer. 6. Procédé selon l'une des revendications précédentes caractérisé en ce que 5 le nombre de groupes M est égal à -2i - V, I étant le nombre de points d'exploration et N le nombre de paramètres du système à déterminer. 7. Procédé selon l'une des revendications précédentes caractérisé en ce que l'étape (101) d'acquisition comprend l'acquisition d'une première pluralité de grandeurs physiques aux instants k à k + R et d'une deuxième pluralité 10 de grandeurs physiques aux instants k + 1 à k+(R + 1), la valeur des grandeurs physiques de la deuxième pluralité à l'instant k + (r + 1) pouvant être déduite de la valeur des grandeurs physiques de la première pluralité à l'instant k + r avec r E [0, R] par l'intermédiaire d'une équation non linéaire dépendante desdits paramètres physiques, et en ce que l'éloignement pi 15 du point ci à l'itération t est égal à la somme des normes des différences entre la valeur de la deuxième pluralité de grandeurs physiques mesurée à l'instant k + (r + 1) et la valeur de la deuxième pluralité de grandeurs physiques calculée à partir de l'équation non linéaire et en fonction de la première pluralité de grandeur physique mesurée à l'instant k + r pour tout 20 r E [0, R], les paramètres de ladite équation étant égaux aux coordonnées de la position xi du point ci à l'itération t. 8. Procédé selon l'une des revendications précédentes caractérisé en ce qu'un nombre de points d'exploration minimum Imin et maximum Imax ainsi qu'un nombre d'itérations maximal Tmax sont définis, et en ce que, lorsque 25 le nombre d'itérations T de l'étape (105) de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration est égal au nombre itérations maximal Tmax, le procédé comprend une étape (107) de suppression et de réévaluation comportant : une sous-étape (1071) de suppression du point d'exploration ayant 30 l'éloignement le plus élevé ; 3060814 27 - une sous-étape (1072) de réorganisation aléatoire des groupes ; - au moins une nouvelle itération de l'étape (105) de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration ci de la pluralité de points d'exploration, une nouvelle étape (107) de 5 suppression et de réévaluation étant mise en oeuvre lorsque le nombre d'itérations T de l'étape (105) de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration est égale au nombre d'itérations maximum Tmax ; le procédé comprenant, lorsque le nombre de points d'exploration I est égal 10 au nombre de points d'exploration minimum Imm, une étape (108) d'ajout et de réévaluation comportant : une sous-étape (1081) d'insertion d'un point d'exploration dont les coordonnées position correspondent à une région non encore explorée ; une sous-étape (1082) de réorganisation aléatoire des groupes ; 15 au moins une nouvelle itération de l'étape (105) de détermination de nouvelles valeurs des vecteurs de chaque point d'exploration ci de la pluralité de points d'exploration, une nouvelle étape (108) d'ajout et de réévaluation étant mise en oeuvre lorsque le nombre d'itérations T de l'étape (105) de détermination de nouvelles valeurs des vecteurs de 20 chaque point d'exploration est égale au nombre d'itérations maximum Tmax ; le procédé comprenant, lorsque le nombre de points d'exploration I est égal au nombre de points d'exploration maximum Im', une nouvelle étape (107) de suppression et de réévaluation. 25 9. Procédé de régulation thermique d'un bâtiment à partir d'un modèle dudit bâtiment caractérisé en ce que les paramètres dudit modèle sont obtenus à l'aide d'un procédé selon l'une des revendications précédentes. 10. Dispositif comprenant des moyens de calculs et des moyens d'acquisition, caractérisé en ce que lesdits moyens sont configurés pour mettre en oeuvre 30 un procédé selon l'une des revendications précédentes.