FR 2913403 A1 20080912 FR 0701580 A 20070305 COMMANDE D'ACTIONNEUR REDUISANT LE NIVEAU DE VIBRATIONS D'UNE STRUCTURE FLEXIBLE ASSOCIEE La présente invention concerne un procédé et dispositif de commande d'au moins un système mécanique présentant au moins un élément de structure flexible et au moins un actionneur ou groupe d'actionneurs, lesdits procédé et dispositif étant destinés à réduire le niveau de vibrations dudit élément de structure flexible tout en commandant ledit actionneur ou groupe d'actionneurs, de sorte à atteindre au moins un objectif assigné audit système mécanique, tel que réaliser un déplacement. Quelques exemples issus de l'état de la technique et des solutions connues qui visent à répondre à ce problème sont considérés ci-dessous. Une première illustration concerne le contrôle d'orbite d'un engin spatial, par exemple un satellite de télécommunication de quelques tonnes que l'on souhaite positionner précisément au voisinage d'un point donné de l'orbite géostationnaire, une fois à poste. Dans le cas d'un satellite à propulsion chimique, des corrections d'orbite doivent être effectuées à intervalles de temps réguliers (typiquement une fois tous les quinze jours), sous forme d'incréments de vitesse (typiquement de 1 à 2 m/s) devant être imprimés au satellite dans chacune des directions tangente à l'orbite et perpendiculaire au plan de l'orbite. Dans le cas considéré ici, ces incréments de vitesse sont réalisés physiquement grâce au système de propulsion monté à bord du satellite, par expulsion de gaz dans des tuyères, l'incrément de vitesse total devant être réalisé étant de 50 m/s par an. Le système de contrôle du satellite calcule, selon un algorithme pré-défini, les ordres d'ouverture à envoyer aux vannes d'éjection de gaz des différentes tuyères afin de réaliser l'incrément de vitesse désiré. Un incrément de vitesse selon une direction donnée peut être réalisé, suivant le principe d'action-réaction, en commandant l'ouverture de la ou des vannes d'une ou plusieurs tuyères éjectant du gaz dans la direction opposée à l'incrément de vitesse désiré. Dans le cas général couvrant la plupart des missions spatiales, les incréments de vitesse à imprimer aux satellites et sondes spatiales pour le contrôle d'orbite s'échelonnent typiquement de quelques cm/s à quelques centaines de m/s.. Dans la très grande majorité des applications spatiales, la commande d'ouverture des vannes est une commande en tout ou rien et non une commande proportionnelle. De ce fait, la poussée résultant de l'ouverture d'une vanne est approximativement constante (de l'ordre de 1 à 20 Newton dans le cas du contrôle d'orbite de satellites de télécommunication, de quelques dixièmes à quelques centaines de Newton dans le cas général). Le contrôle de l'incrément de vitesse réalisé s'effectue en boucle ouverte, à partir d'un modèle de poussée des tuyères, en contrôlant la durée totale d'ouverture des vannes de chaque tuyère active. Le plus simple consiste à commander l'ouverture des vannes de façon continue pendant le temps juste nécessaire pour réaliser la totalité de l'incrément de vitesse désiré, typiquement quelques centaines de secondes dans le cas considéré. Cependant, l'ouverture continue d'une ou plusieurs vannes pendant une telle durée n'est pas recommandée, car les couples perturbateurs résultant des mésalignements des jets des tuyères par rapport à leur orientation nominale pourraient dépasser la limite admissible et engendrer des dépointages transitoires du satellite et de ses charges utiles allant au-delà de la limite autorisée par la mission (typiquement 0.05 à 0.1 degré pour l'application considérée). Suivant l'état de l'art, afin de remédier à ce problème, l'incrément de vitesse est réalisé non pas d'un coup, en totalité, en ouvrant en continu les tuyères concernées pendant le laps de temps juste nécessaire, mais plutôt par une succession de petits incréments répartis sur un laps de temps plus important, typiquement deux à quatre fois la durée juste nécessaire dans le cas d'une poussée continue. L'effet produit est en moyenne équivalent à une diminution d'un même facteur de la poussée équivalente imprimée au satellite pendant la manoeuvre (réduisant d'autant l'effet transitoire néfaste des couples perturbateurs), ce qui ne peut être produit par une ouverture partielle des vannes du système de propulsion qui, rappelons-le, travaillent selon nos hypothèses en tout ou rien. De façon pratique, suivant l'état de l'art, les ordres de commande d'ouverture des tuyères sont échantillonnés à une période d'échantillonnage T. Dans chaque intervalle d'échantillonnage compris entre deux instants d'échantillonnage successifs, chaque vanne de tuyère active est commandée périodiquement à la période T en position ouverte dans un intervalle de temps de durée AT strictement inférieure à la période d'échantillonnage T, et typiquement égale à une fraction entière de T. Pour des raisons qui apparaîtront ci-après, l'ordre d'ouverture des vannes est de préférence centré au milieu de l'intervalle d'échantillonnage. Ainsi, si on considère conventionnellement que la commande envoyée au système d'ouverture de chaque vanne est égale à 1 pour commander une ouverture totale, et 0 pour commander une fermeture totale, le signal de commande envoyé à chaque vanne du système de propulsion a la forme d'une succession périodique, de période T, de petits incréments de valeur 1 et de durée AT centrés au milieu des intervalles d'échantillonnage, la commande étant nulle en dehors de ces durées AT (voir ci-après l'exemple de la figure 2a). Il est à noter que la période d'échantillonnage T doit être choisie la plus grande possible (typiquement quelques secondes) afin de minimiser, pendant la durée de vie du satellite qui est de l'ordre de quinze ans, le nombre total de transitoires ouverture/fermeture des vannes de tuyères qui sont sensibles à ce paramètre. Par ailleurs, cette période d'échantillonnage doit être suffisamment grande pour que les incréments à réaliser soient d'un ordre de grandeur plus grand que la durée minimum d'ouverture des tuyères (ce qu'on appelle le Minimum Impulse Bit), en dessous de laquelle il y a une perte notable d'efficacité et une surconsommation de carburant importante. L'avantage d'un tel procédé connu est donc d'étaler la réalisation de l'incrément de vitesse sur une durée plus grande d'un facteur T/AT que la durée juste nécessaire dans le cas d'une ouverture en continu des tuyères, ce qui a pour effet moyen de diminuer du même facteur la poussée équivalente imprimée au satellite, et donc de réduire d'autant les transitoires indésirables dus aux couples perturbateurs. Cependant, en contre partie de cet avantage décisif qu'il s'agit de conserver, ce procédé suivant l'état de l'art présente une limitation importante, qui peut affaiblir, voire annuler le bénéfice escompté. Cette limitation est due au caractère périodique de la commande d'ouverture des vannes ainsi créé par l'étalement de la manœuvre tel qu'expliqué précédemment. Par exemple, pour des valeurs T = 1 seconde et AT = 0.2 seconde, la figure 4a (décrite ci-après) présente le spectre de Fourier du profil de commande selon l'état de l'art présenté ci-dessus. Ce spectre montre que l'énergie de la commande est concentrée sur les multiples entiers k/T de la fréquence d'échantillonnage, entraînant un risque de résonance entre les incréments de poussée imprimés au satellite pendant la manoeuvre et les modes de vibration d'éléments flexibles, comme, par exemple, les grands générateurs solaires, d'autant que la fréquence d'échantillonnage 1/T est relativement basse, pour les raisons expliquées plus haut. Il faut insister sur le fait qu'une telle résonance entre une commande en boucle ouverte de réalisation d'incréments de vitesse et des modes flexibles de grands appendices comme des générateurs solaires peut être extrêmement néfaste, du fait que ces modes flexibles sont très faiblement amortis (facteur d'amortissement de l'ordre de 0,001). Suivant l'état de l'art, on cherche à compenser cet effet en sélectionnant urne fréquence d'échantillonnage 1/T telle que ses multiples entiers soient tous suffisamment éloignés des fréquences des modes flexibles des générateurs solaires afin d'éviter toute coïncidence ou rapprochement néfaste. Cependant, ces fréquences de modes flexibles sont connues avec une incertitude importante de l'ordre de 10 à 20 surtout pour les hautes fréquences, ce qui en pratique ne permet pas de réduire le risque de coïncidence. De plus, la densité des modes flexibles peut être telle qu'il soit impossible de trouver un bon emplacement entre ces modes. Enfin, il serait avantageux d'augmenter la fréquence d'échantillonnage afin de repousser le problème de coïncidence vers des fréquences plus hautes, là où l'énergie des modes flexibles est moins importante et leurs effets plus limités, mais cette augmentation est contrainte par les effets indiqués précédemment (nombre de cycles ouverture/fermeture des tuyères, perte d'efficacité induite par des incréments trop courts). Suivant l'état de l'art, on peut également chercher à augmenter la bande passante du contrôle d'attitude du satellite afin de contrôler au moins les premiers modes des appendices flexibles et ainsi réduire l'effet néfaste des éventuelles résonances. Cependant, cette solution est également limitée par plusieurs facteurs : niveau d'autorité de contrôle disponible (problème de saturation), risque de déstabilisation du contrôle d'attitude qui doit gérer un grand nombre de modes flexibles mal connus et très peu amortis, risque subsistant de résonance avec des modes flexibles situés en dehors de la bande passante du contrôle d'attitude. Un autre exemple d'approche conventionnelle concerne le pointage de grande charge utile présentant des modes flexibles d'assez basse fréquence. C'est le cas, par exemple, de la mise en rotation des grands générateurs solaires de satellites de télécommunication comme ceux précédemment évoqués, qui doivent exécuter un tour en 24 heures par rapport à la plate-forme du satellite afin de rester pointés vers le soleil tout au long de l'orbite, alors que le satellite reste pointé vers la Terre. Ces générateurs solaires sont mis en rotation par des moteurs pas à pas. La commande de ces moteurs est échantillonnée à une période T. Toutes les T secondes, une commande d'incrément angulaire est envoyée au moteur afin de réaliser la rotation souhaitée (comme décrit ci-après en référence à la figure 3a). Les incréments ont nominalement la même valeur. Là encore, du fait de la périodicité des incréments, il y a un risque éventuel de résonance entre la commande et les modes flexibles de l'organe mécanique commandé. En réponse aux limitations de l'état de l'art mentionnées ci-dessus, l'objet de la présente invention est de fournir un procédé de commande d'au moins un système mécanique comportant au moins un élément de structure flexible et au moins un actionneur ou groupe d'actionneurs commandé au moyen d'un signal de commande de façon à atteindre au moins un objectif assigné audit système mécanique, ledit procédé étant destiné à réduire le niveau de vibration dudit élément de structure flexible provoqué par l'activation dudit actionneur ou groupe d'actionneurs. A cet effet, l'invention propose un procédé de commande d'au moins un système mécanique présentant au moins un élément de structure flexible et au moins un actionneur ou groupe d'actionneurs, ledit procédé étant destiné à réduire le niveau de vibrations dudit élément de structure flexible tout en commandant ledit actionneur ou groupe d'actionneurs de sorte à atteindre au moins un objectif assigné audit système mécanique, tel que réaliser un déplacement, ladite commande étant constituée d'incréments devant être réalisés à une période T dite période d'échantillonnage, chacun desdits incréments étant constant en amplitude dans un intervalle de temps dit intervalle d'échantillonnage s'étendant entre deux instants d'échantillonnage successifs et chacun desdits incréments s'appliquant pendant un intervalle d'application d'incrément d'une durée inférieure ou égale à ladite période d'échantillonnage T, ledit procédé étant caractérisé en ce que, dans chaque intervalle d'échantillonnage, l'un au moins des instants initial et final d'application dudit incrément est modifié par l'ajout d'un écart temporel dT variable, aléatoire ou pseudoaléatoire, dit écart aléatoire. Ce qu'on entend précisément par écart aléatoire, en relation avec un intervalle d'échantillonnage respectif, est donc une valeur respective d'une suite discrète de valeurs aléatoires ayant avantageusement la même loi de probabilité, et étant avantageusement indépendantes statistiquement les unes des autres. Cependant, étant donnée une loi de probabilité définie mathématiquement, la réalisation par un programme informatique connu de valeurs aléatoires obéissant à cette loi de probabilité idéale conduit à une suite de nombres qui ont des propriétés statistiques se rapprochant de ladite loi de probabilité idéale, sans toutefois y obéir parfaitement. Une telle suite de nombres est dite pseudoaléatoire , et on attribue à ces nombres, par abus de langage, la loi de probabilité idéale dont ils sont censés être issus. L'intérêt du procédé selon l'invention est le suivant. La commande nominale (avant application du procédé de l'invention) est une répétition d'incréments (qui, dans le cas général, peuvent ne pas être tous parfaitement égaux) à une période T. Même si ce signal de commande n'est pas forcément parfaitement périodique (incréments non égaux, imperfections de réalisation, etc), il a néanmoins un fort contenu fréquentiel centré sur la fréquence d'échantillonnage 1/T et ses harmoniques, ce qui conduit aux risques mentionnés ci--dessus. L'ajout selon l'invention d'un écart aléatoire , selon la définition ci-dessus, à au moins un instant d'application des incréments de la commande permet de modifier le contenu fréquentiel de la commande en atténuant son caractère périodique. Suivant la loi de probabilité choisie, tout ou partie des pics d'énergie centrés sur les multiples entiers de la fréquence d'échantillonnage 1/T seront atténués, voire supprimés, comme cela est décrit ci-après pour des exemples particuliers. La loi de probabilité de l'écart aléatoire est choisie à bon escient, pour minimiser les pics d'énergie en priorité là où les résonances pourraient être les plus préjudiciables, c'est-à-dire proches des fréquences de vibration des éléments flexibles du système mécanique. L'homme de l'art peut aisément généraliser le procédé de l'invention lorsque, non pas la commande elle-même, mais une de ses dérivées temporelles possède le profil temporel objet du procédé (suite périodique d'incréments). Dans ce cas, le procédé peut être appliqué à ladite dérivée possédant ces caractéristiques, et la commande est calculée comme l'intégrale d'ordre correspondant de cette dérivée modifiée. Selon une réalisation particulière de l'invention, pour chaque intervalle d'échantillonnage, la durée d'application desdits incréments est égale à la période d'échantillonnage T, et un écart aléatoire selon l'invention est ajouté à l'instant initial d'application desdits incréments (l'instant final d'un incrément étant l'instant initial du suivant). Selon une autre réalisation particulière de l'invention, pour chaque intervalle d'échantillonnage, la durée d'application desdits incréments est strictement inférieure à la période d'échantillonnage T, et un écart aléatoire est ajouté soit à l'instant initial, soit à l'instant final, soit aux deux instants initial et final d'application desdits incréments, avec des valeurs pouvant être différentes dans ce dernier cas. En variantes de la réalisation précédente, dans chaque intervalle d'échantillonnage, le même écart aléatoire peut être ajouté à la fois à l'instant initial et à l'instant final d'application desdits incréments, ce qui a pour effet de translater dans le temps l'ensemble de l'incrément ou créneau de commande devant être réalisé, d'une valeur égale à l'écart aléatoire appliqué. Comme expliqué précédemment, il est avantageux de définir la loi de probabilité choisie pour déterminer ledit écart aléatoire de sorte à modifier au moins un des pics d'énergie de la commande, en diminuant, déplaçant ou supprimant au moins un des pics d'énergie indésirables, et en particulier ceux qui coïncident, ou sont proches des fréquences propres d'éléments flexibles du système mécanique considéré. Avantageusement, la valeur absolue maximum de l'écart aléatoire est inférieure à une valeur fixée, de préférence une fraction entière de la période d'échantillonnage, de sorte que l'incrément modifié par l'écart aléatoire puisse être réalisé entièrement à l'intérieur de son intervalle d'échantillonnage, sans déborder sur le suivant. Dans la plupart des cas, l'écart aléatoire est à moyenne nulle, de sorte que l'effet moyen de la commande modifiée selon le procédé de l'invention soit le même que celui de la commande initiale ou nominale. Pour des raisons pratiques d'implantation informatique, on peut sélectionner les valeurs de l'écart aléatoire parmi les multiples entiers d'une valeur pré- définie T', celle-ci pouvant être un diviseur entier de la période d'échantillonnage T. Dans ce cas, la fréquence 1/T' peut être choisie de sorte à limiter les éventuelles résonances avec les modes de vibrations dudit élément de structure flexible, et de préférence en ne coïncidant pas avec les fréquences propres desdits modes de vibrations. Afin de simplifier la réalisation pratique de l'invention, une loi de probabilité choisie pour déterminer ledit écart aléatoire est entièrement définie par un nombre fini de paramètres calculés de sorte à modifier au moins un des pics d'énergie de la commande. Dans une réalisation particulière de l'invention suivant ce principe général, la loi de probabilité déterminant ledit écart aléatoire est une loi de probabilité discrète à m valeurs, ces M valeurs et leurs probabilités étant calculées de sorte à modifier au moins un des pics d'énergie de la commande. Par exemple, la loi de probabilité déterminant ledit écart aléatoire est une loi binomiale, prenant des valeurs extrêmes d'écart aléatoire +dTmax avec une probabilité de 1/2, et -dTmax avec une probabilité de 1/2. Dans le cas d'une translation aléatoire des incréments, l'effet est de diminuer l'énergie de la commande aux fréquences k/T, multiples entières de la fréquence d'échantillonnage 1/T, dans un rapport cos(GJ.dTmax) où 0) = 2. 7t.k/T, et k est un nombre entier naturel. Dans une autre réalisation, la loi de probabilité déterminant l'écart aléatoire est une loi discrète à m valeurs équiprobables en progression arithmétique entre une valeur minimale dTmax et une valeur maximale +dTmax des écarts aléatoires, avec dTmax = T.(M-1)/2M. Dans le cas d'une translation aléatoire des incréments, l'effet est que le spectre de la commande ne peut pas comporter de raies aux fréquences 1/T, 2/T, ... (M-1)/T, ainsi que (Mp+l)/T, (Mp+2)/T, ... (Mp+M-1)/T où p est un entier naturel. On peut également choisir comme loi de probabilité pour déterminer l'écart aléatoire selon le procédé de l'invention, une loi uniforme entre des valeurs extrêmes d'écart aléatoires dTmin et dTmax, les valeurs dTmin et dTmax étant choisies pour modifier au moins un pic d'énergie de la commande. Plus particulièrement, on peut choisir, comme loi de probabilité déterminer l'écart aléatoire, une loi uniforme entre des valeurs maximales -dTmax et dTmax. Dans le cas d'une translation aléatoire des incréments, l'effet est de diminuer l'énergie de la commande aux fréquences k/T multiples entières de la fréquence d'échantillonnage 1/T, dans un rapport sin(().dTmax) / (G).dTmax) où O) = 2. 7L.k/T. Dans ce dernier cas, on peut choisir dTmax = T/2, de sorte que le spectre fréquentiel de la commande ne comporte pas de raie à la fréquence d'échantillonnage 1/T. Un autre choix peut être dTmax = T/4, de sorte que le spectre fréquentiel de la commande ne comporte pas de raie aux fréquences multiples impaires de la fréquence d'échantillonnage 1/T. A des fins d'implantation numérique du procédé, les écarts aléatoires selon l'invention peuvent prendre les valeurs de nombres successifs d'au moins une séquence aléatoire ou pseudo-aléatoire de nombres produite à chaque période d'échantillonnage par au moins un programme informatique générateur de nombres aléatoires ou pseudoaléatoires. La mise en oeuvre du procédé est encore plus simple si les écarts aléatoires selon le procédé prennent les valeurs de nombres successifs d'une séquence longue mais finie de nombres présentant des caractéristiques pseudoaléatoires ou aléatoires proches de caractéristiques désirées, ladite séquence étant pré-calculée et appliquée suivant le procédé, puis reproduite en boucle une fois qu'elle est achevée. L'invention a simultanément pour objet un produit programme d'ordinateur comprenant des instructions de code adaptées à la mise en oeuvre du procédé selon l'invention et tel que défini ci-dessus, lorsqu'il est chargé et exécuté sur des moyens informatiques. De plus, un dispositif de mise en oeuvre du procédé selon l'invention comprend des moyens de calcul et informatiques permettant l'exécution du produit programme selon l'invention et tel que défini ci-dessus, et des moyens électriques et/ou électroniques permettant de transformer des ordres de commande numériques élaborés par les moyens de calcul et informatiques en signaux électriques envoyés à au moins un actionneur, ledit dispositif étant tel que : - au moins un programme informatique dudit produit programme est mis en oeuvre par lesdits moyens de calcul et informatiques pour produire lesdits incréments de ladite commande d'au moins un desdits actionneur et groupe d'actionneurs dudit système mécanique de sorte à atteindre au moins un objectif assigné audit système mécanique, tel que réaliser un déplacement ; - ledit programme informatique est mis en oeuvre pour produire, à chaque période d'échantillonnage T, ledit écart aléatoire suivant une loi de probabilité se rapprochant au mieux d'une loi de probabilité désirée - ledit programme informatique est mis en oeuvre pour modifier, à chaque période d'échantillonnage, le ou les instants d'application dudit incrément en lui ajoutant ledit écart aléatoire selon le procédé de l'invention et produire ainsi la commande numérique à appliquer effectivement audit système mécanique, et - lesdits moyens électriques et/ou électroniques transforment ladite commande numérique produite par ledit programme informatique en signaux électriques permettant à au moins un desdits actionneur et groupe d'actionneurs de réaliser physiquement la commande assignée. Dans le cas où la commande selon le procédé de l'invention est modifiée en utilisant une séquence longue mais finie d'écarts aléatoires ayant les caractéristiques aléatoires désirées, ladite sequence est calculée préalablement, puis avantageusement stockée dans une mémoire desdits moyens de calcul et informatiques dudit dispositif, ces écarts aléatoires étant, successivement et suivant un ordre établi, ajoutés à au moins un instant d'application desdits incréments de ladite commande, ladite séquence étant reproduite en boucle une fois qu'elle est achevée. Il peut être avantageux qu'au moins une partie des paramètres intervenant dans le calcul des écarts aléatoires par ledit programme informatique dudit dispositif, tels que les paramètres d'une loi de probabilité, ou une séquence de valeurs aléatoires, peut être modifiée au cours du fonctionnement dudit système mécanique, afin par exemple d'optimiser l'efficacité du procédé (si on bénéficie d'une connaissance a posteriori des modes de vibration des éléments flexibles du système, et en particulier leur fréquence propre). Par exemple, le dispositif de mise en oeuvre du procédé selon l'invention peut être tel qu'il comprend au moins un capteur dédié de mesure des vibrations d'au moins un élément de structure flexible présent dans le système mécanique, des moyens de calcul du spectre fréquentiel de ces vibrations, et des moyens de modification de paramètres intervenant dans le calcul desdits écarts aléatoires afin de réduire le niveau desdites vibrations en modifiant au moins un des pics d'énergie du spectre fréquentiel de la commande. Comme déjà évoqué, le procédé de l'invention s'applique à la commande d'ouverture d'au moins une tuyère d'un système de contrôle d'orbite, de stabilisation et/ou d'orientation d'un engin se déplaçant dans l'espace, tel que lanceur, satellite artificiel, sonde interplanétaire, etc, ledit engin comportant au moins un élément de structure flexible, tel de grands appendices, antenne, panneau solaire et structure primaire, etc. Plus précisément, cette application du procédé de l'invention consiste à : - commander l'ouverture d'au moins une tuyère afin 30 d'imprimer audit engin un incrément de vitesse permettant de corriger la trajectoire dudit engin ; -réaliser nominalement cet incrément de vitesse en une suite de petits incréments effectués périodiquement à une période T ; - réaliser lesdits petits incréments par la commande d'ouverture d'au moins une tuyère, ladite commande étant conventionnellement égale à 0 pour un ordre de fermeture de tuyère, et à 1 pour un ordre d'ouverture de tuyère ; - avant la mise en oeuvre dudit procédé, calculer une commande dite nominale de façon à être périodique de période T dans l'intervalle de temps [Ti, Ti + N.T], où N est un nombre entier, de réalisation d'une impulsion totale, et, dans cet intervalle de temps, égale à un incrément de 1 sur tous les intervalles de temps [Ti + k.T + Tl, Ti + k.T + T2], où k est un nombre entier inférieur à N, Tl et T2 étant des constantes telles que 0 Généralement, cet écart temporel supplémentaire est petit (typiquement d'un facteur 10) par rapport à la durée T2-Tl d'application de l'incrément. Dans un autre type d'application du procédé selon l'invention, on commande la rotation d'au moins un moteur pas à pas pour le pointage d'au moins une charge utile d'un engin spatial, telle qu'au moins un panneau solaire, une antenne, un instrument, etc, ladite charge utile présentant au moins un élément de structure flexible. Dans ce cas, l'application du procédé de l'invention est telle que : - une commande a de l'angle entre le stator et le rotor dudit moteur est échantillonnée à une période d'échantillonnage T, et s'incrémente d'une certaine valeur à chaque période d'échantillonnage, - une commande nominale a se présente sous la forme : ak = k_1 + (d a) k où a k représente la valeur constante de la commande entre les instants k.T et (k+l).T, où k est un nombre entier, et (d a) k est la valeur de l'incrément à réaliser à la date k.T, les incréments 30 s'effectuant par convention aux instants d'échantillonnage, - ladite commande est modifiée selon le procédé de l'invention en ajoutant, aux instants k.T de réalisation des incréments, un écart temporel aléatoire ayant des caractéristiques désirées dans le but de minimiser des résonances avec des modes de vibration dudit élément de structure flexible. Un cas particulier de cette application est tel que les incréments (d OC) k sont tous égaux, comme dans le cas du maintien en rotation à vitesse constante d'un panneau solaire. De même que précédemment, on peut superposer un signal de commande en boucle fermée au signal de commande en boucle ouverte que l'on veut modifier selon l'invention. Dans ce cas, l'application du procédé de l'invention est telle que les incréments (d OC)k nominalement égaux sont légèrement modifiés, de manière connue, pour réaliser le contrôle d'au moins un phénomène physique résultant de la rotation du moteur, comme par exemple le pilotage d'un satellite à voile solaire dans le cas de la commande d'un panneau solaire muni d'au moins un volet. Dans le cas de la commande d'un moteur pas à pas, on peut appliquer le procédé de l'invention à d'autres profils de rotation plus complexes que de simples rotations à vitesse constante, comme par exemple des basculements alternatifs à basse fréquence. Dans ce cas, l'application du procédé selon l'invention est telle que le spectre fréquentiel de la commande nominale à appliquer au moteur se situe principalement à des fréquences très inférieures à la fréquence d'échantillonnage, en particulier pour réaliser un mouvement de balayage à faible vitesse. D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante de quelques modes de réalisation, donnés à titre d'exemples non limitatifs en regard des dessins joints : - les figures la et lb représentent respectivement en perspective et de manière schématique un satellite à panneaux solaires déployés tels qu'en orbite, et dans lequel est mis en oeuvre l'invention ; - les figures 2a, 2b et 2c sont, respectivement, un chronogramme représentant l'évolution temporelle d'un signal de commande conventionnel dans le cas d'un contrôle d'orbite du satellite des figures la et lb, et deux chronogrammes représentant l'évolution temporelle du même signal de commande modifié selon deux formes de réalisation de l'invention ; - les figures 3a et 3b sont respectivement un chronogramme représentant l'évolution temporelle d'un signal de commande conventionnel dans le cas du pointage d'un générateur solaire par un moteur pas à pas, et un chronogramme représentant l'évolution temporelle du signal de commande modifié selon une forme de réalisation de l'invention ; et - les figures 4a, 4b et 4c sont des représentations spectrales respectivement du signal de commande conventionnelle correspondant à celui représenté sur la figure 2a, et du même signal de commande modifié selon la forme de réalisation de l'invention correspondant à la figure 2c, avec un écart aléatoire suivant une loi de probabilité discrète à 3 valeurs équiprobables -dTmax, 0 et +dTmax, avec dTmax pour la figure 4c. Les figures la et lb représentent de manière schématique un satellite 1 qui comprend une pluralité de générateurs ou panneaux solaires 2, utiles pour l'alimentation électrique, et des charges utiles comme des antennes telles que 5. Ce satellite 1 peut contrôler son attitude et/ou son orbite au moyen de tuyères 3 (voir figures la) qui peuvent, dans leur globalité, exercer une pour la figure 4b, et dTmax T 4 T 3 poussée dans une ou plusieurs directions. Certains éléments du satellite 1 comme les générateurs solaires 2 ou des charges utiles, telles que les antennes 5, peuvent être de très grandes dimensions, avec des modes de vibration à très basse fréquence qui peuvent être excités par le contrôle d'attitude et d'orbite du satellite 1. Les tuyères 3 sont actionnées par des vannes d'alimentation 6 (voir figure lb) qui permettent l'ouverture des tuyères 3 afin d'exercer une poussée dans une direction prédéterminée. Cela permet d'ajuster l'orbite du satellite 1 afin de maintenir une altitude correcte, ou bien d'améliorer les conditions de pointage du satellite, par exemple dans le cadre d'une communication optique avec un autre satellite, ou d'une communication radio avec une station terrestre. La commande de chaque tuyère 3 est effectuée par ouvertures successives pendant une durée prédéterminée de la vanne d'alimentation 6 correspondante. Ceci a pour effet de générer une poussée constante pendant la durée d'ouverture de la vanne 6. On peut donc moduler l'intensité de la force appliquée par les tuyères 3 uniquement par la durée d'ouverture des vannes 6. L'orientation des panneaux solaires 2 est contrôlée par des moteurs pas à pas 4. Il s'agit de moteurs électriques dont le signal d'entrée de chacun n'est pas un simple courant ou tension analogique, mais une valeur représentative de l'angle du rotor du moteur 4 par rapport à un point de référence. Le pas minimal de rotation de ces moteurs 4 dépend du type de moteur considéré. Ces moteurs 4 permettent une régulation simple de l'orientation des panneaux solaires 2 en fonction de l'orientation du soleil vis-à-vis du satellite 1. Les moteurs pas à pas 4 et les vannes d'alimentation 6 des tuyères 3 sont des actionneurs commandés par un dispositif 7 d'application du signal de commande. Ce dispositif (DRV) 7 est contrôlé par une unité de traitement centrale (CPU) 8 électronique, qui calcule et génère le signal de commande qu'applique le dispositif d'application 7 aux différents actionneurs 4 et 6. L'unité centrale 8 peut, selon un mode de réalisation, générer une séquence d'écarts temporels aléatoires (dTi) 28 ou 35 (voir figures 2b, 2c et 3b) répondant par exemple à une loi de probabilité prédéterminée, afin de générer un signal de commande, tel que 27 sur les figures 2b et 2c, à appliquer par le dispositif 7 d'application de signal de commande. Ce signal de commande est formé par une succession d'incréments 23 ou 33 (figures 2b, 2c et 3b). Dans le cas de la commande de vannes 6 de tuyères 3, ces incréments 23 (voir figure 2b) ont la forme d'une impulsion d'amplitude et de durée sensiblement constantes. Dans le cas de commande des moteurs pas à pas 4 (figures 3a et 3b), chaque commande correspond à un échelon angulaire d'un incrément 33. Le signal de commande, dont le profil temporel de commande nominale (avant application du procédé de l'invention) est représenté en 21 ou 31 en traits pleins sur les figures 2a et 3a (qui représentent les commandes nominales correspondantes), et en pointillés sur les figures 2b, 2c et 3b (qui représentent les commandes modifiées par l'application du procédé de l'invention), comprend une succession d'intervalles d'échantillonnage 22 ou 32, répétés périodiquement, à la période T, et donc avec une fréquence f égale à 1/T, déterminée. A l'intérieur de chacun des intervalles d'échantillonnage 22 ou 32, un incrément 23 ou 33 est généré, et maintenu appliqué, soit pendant un intervalle 24 d'application des incréments 23 qui est de durée OT strictement inférieure à la durée T de l'intervalle d'échantillonnage 22 dans le cas de la commande des vannes 6 des tuyères 3, soit pendant l'intervalle d'échantillonnage 32 de durée T dans le cas de la commande des moteurs pas à pas 4. Cependant, selon un premier exemple de mise en oeuvre de l'invention, l'instant initial 25 ou 34 ou éventuellement final 26 d'application de l'incrément 23 ou 33 est modifié, comme en 29 sur les figures 2b et 2c, de manière aléatoire par l'unité centrale 8. A cet effet, l'unité centrale 8 permet de retarder ou d'avancer d'un écart temporel aléatoire (dTi) 28 ou 35 l'instant initial 25 ou 34 ou éventuellement final 26 d'application de chaque incrément 23 ou 33, par rapport à un instant de référence dans son intervalle d'échantillonnage 22 ou 32 respectif. En général, pour une commande nominale selon la figure 2a, on prend comme instant de référence, le milieu de l'intervalle d'échantillonnage 22. Selon un autre mode de réalisation, l'unité centrale 8 est connectée à une base de données (DB) 9 formée par des moyens de stockage de données. A l'intérieur de cette base de données 9, il existe des séquences d'écarts temporels dTi pseudo aléatoires. Ces séquences peuvent être considérées comme des suites de variables aléatoires, si le nombre d'éléments dans la séquence est suffisant, et que l'on répète cette séquence en boucle. De telles séquences ou suites de variables pseudo aléatoires peuvent être générées par des programmes informatiques connus. Par ailleurs, le satellite 1 peut de plus comprendre des capteurs de vibrations 10, disposés au niveau des éléments flexibles comme par exemple les générateurs solaires 2, et capables de mesurer avec précision leurs fréquences de vibrations. Les capteurs 10 sont reliés à l'unité centrale 8 afin de calculer le spectre fréquentiel desdites vibrations de telle sorte que l'unité centrale 8 génère les séquences d'écarts dTi aléatoires ou pseudo aléatoires en fonction du spectre de fréquence et du niveau de vibration détecté par les capteurs 10. De manière plus détaillée, sur la figure 2a est représenté le profil temporel 21 d'un signal de commande nominale (conventionnelle) utilisé, par exemple, pour le contrôle d'orbite d'un satellite 1. Dans ce cas, le signal 21 représente les ordres d'ouverture et de fermeture de vannes 6 de tuyères 3. Par convention, ce signal vaut 1 pour commander une ouverture, et 0 pour une fermeture. Ce signal comprend des incréments 23 (ordre d'ouverture) à réaliser périodiquement à une période T, dite période d'échantillonnage. Dans chaque intervalle d'échantillonnage 22, l'incrément 23 se réalise sur un intervalle 24 d'application d'incrément de durée AT qui, dans ce cas particulier, a une durée strictement inférieure à la période d'échantillonnage T. L'intervalle 24 d'application de l'incrément est, de préférence, centré sur l'intervalle d'échantillonnage 22, mais la figure 2a présente le cas général où l'instant initial 25 d'application de l'incrément 23 et l'instant final 26 d'application de l'incrément 23 ne sont pas symétriques par rapport au centre de l'intervalle d'échantillonnage 22. Sur la figure 2b est représenté le profil temporel 27 d'un signal de commande modifié à partir de celui de la figure 2a, en appliquant une première réalisation particulière du procédé selon l'invention, qui consiste à ajouter à l'instant initial 25 d'application de l'incrément 23 un écart temporel (dTi) 28 aléatoire, ou pseudo aléatoire, pouvant être positif ou négatif, et donc correspondre respectivement à un retard ou à une avance dans le temps, de sorte que le nouvel instant initial ou instant initial modifié 29 se produise aléatoirement autour de la valeur de l'ancien instant initial 25. Une réalisation alternative, mais préférée, du procédé est visualisée sur la figure 2c, où, dans ce cas, le même écart aléatoire (dTi) 28, positif ou négatif, est appliqué sur l'instant initial 25 et sur l'instant final 26 d'application de l'incrément 23, de sorte que le nouvel intervalle d'application de l'incrément s'étendant entre les nouveaux instants initial 29 et final 30 (modifiés) soit juste le translaté de l'ancien intervalle d'application 24 d'un écart temporel aléatoire (dTi) 28. Sur les figures 2b et 2c, dTi est positif lorsque dTi+l et dTi+2 sont négatifs. Dans une autre réalisation, applicable en particulier à la commande d'un moteur pas à pas 4 d'entraînement d'un générateur solaire 2, le profil temporel 31 du signal de commande nominale (avant application du procédé de l'invention) représenté sur la figure 3a est constitué d'incréments 33 qui sont successivement appliqués sur l'ensemble d'intervalles d'échantillonnage 32, successifs de durée T, les instants d'application de ces incréments 33 étant conventionnellement les instants d'échantillonnage 34. La figure 3b montre comment le signal de commande nominal de la figure 3a est modifié selon une réalisation du procédé de l'invention. Un écart temporel aléatoire (dTi) 35 (qui peut être positif ou négatif et correspondre respectivement à un retard ou a une avance dans le temps), et variable aléatoirement d'un intervalle d'échantillonnage 32 au suivant, est appliqué à chaque instant d'application 34 de l'incrément 33, sachant que, dans ce cas, l'instant final d'application d'un incrément correspond à l'instant initial 34 d'application de l'incrément suivant. Sur la figure 3b, dTi-1, dTi+l et dTi+2 sont positifs, et dTi est négatif, mais la somme des écarts aléatoires peut être nulle, c'est-à-dire que les écarts (dTi) sont à moyenne nulle. Afin de comprendre l'influence de l'écart aléatoire dTi sur le spectre du signal de commande, il convient de comparer le spectre du signal de commande conventionnel et celui du signal de commande modifié après application du procédé selon l'invention. En reprenant le cas du signal de commande conventionnel de la figure 2a, noté f(t) , l'énergie de ce signal périodique est concentrée aux fréquences multiples de la fréquence d'échantillonnage 1/T, comme le montre son spectre fréquentiel représenté sur la figure 4a. A la fréquence k/T, où k est un nombre entier naturel, le spectre du signal est précisément donné par le coefficient de Fourier Fk de rang k du signal où, par définition : Fk = ÏT+I)l f(t)e'~`dt avec w=2zk/T Dans le cas considéré : F _ATsin(wOT/2)e_'"''~ où w=21tk/T k T wAT / 2 et on retrouve les valeurs numériques de la figure 4a, où est représentée la norme des coefficients Fk dans le cas où T = 1 seconde et AT = 0.2 seconde (et dTi=0 puisqu'il s'agit du signal de commande nominale). Afin de généraliser au cas du signal modifié par le procédé suivant l'invention, on note que le coefficient de Fourier de rang k peut s'écrire également, du fait de la périodicité de la fonction f(t): F = 1 Io' f(t)e Fk (2N+1)T o' Un calcul montre qu'on peut passer à la limite: Fk = lim 1n+> ~f(t)e-'0'dt N-9-(2N+ 1)T fi! 25 Cette formule peut s'appliquer à toute réalisation particulière du signal modifié, selon le procédé, par l'ajout d'un écart aléatoire (nous noterons f(t) un tel signal modifié). On fait, tout d'abord l'hypothèse qu'à chaque intervalle d'échantillonnage 22, le même écart aléatoire dT(n) est appliqué à l'instant initial 25 et à l'instant final 26 d'application de l'incrément. Dans ce cas, qui est celui de la figure 2c, la transformée de Fourier du signal modifié à la fréquence k/T, noté Fk, est 10 donnée par la formule : Fk = hm 1 J-T n f(t) e-'~'dt N-* (2N + 1)T n = lim AT sin(COOT/2) ' e-rwTl2 m(nT+dT(n)) Nom (2N + 1)T T wLT /2 ä__N où w =27rk/T. La valeur moyenne de ce coefficient pour différentes réalisations de l'écart aléatoire dT(n) 28 est 15 donnée par : 1 A T sin(wAT /2) _1 mTl2 n=N -i cotir md'f (n) E[Fk 1= lim e ~ e Me- Nù (2N + 1)T T (OOT / 2 n=-N où on reconnaît le coefficient de Fourier du signal d'origine modifié par un terme venant de l'écart aléatoire introduit suivant l'invention : 20 E[Fk ]= Fk E[e-'"d''(n) Si, par exemple, on prend comme loi de probabilité l'écart aléatoire dT, une loi discrète à m valeurs dTm même probabilité, dans ce cas, l'harmonique de rang k du signal d'origine sera multiplié par la valeur : 1 M ( d' \k 25 M m=I i Sachant que de de M k (e 2ivrm/M) =0 pour k=1,...,M-1 m=1 on en déduit qu'une distribution équi-répartie des valeurs dTm entre deux valeurs -dTmax et dTmax, à savoir : dTm = ù dTnax + 2 dTmax pour m =1,..., M -1 i m-1 ~M-1, permet d'éliminer les M-1 premières harmoniques du signal d'origine. Par exemple, pour M=3, considérons un écart aléatoire prenant les valeurs -dTmax, 0, et dTmax avec chacune une probabilité de 1/3. Dans le cas où dTmax/T = 1/3, les deux premières harmoniques du signal d'origine, aux fréquences 1/T et 2/T, sont éliminées, comme le montre la figure 4c qui représente le spectre du signal dans ce cas. Il est à noter qu'un écart aléatoire moins important, correspondant à dTmax/T = 1/4 , sans éliminer strictement les premières harmoniques, les aurait tout de même atténuées fortement, comme le montre la figure 4b qui représente le spectre du signal dans ce cas. De façon générale, on calcule, selon le procédé une loi de probabilité afin de générer des écarts temporels (dTi) aléatoires vis-à-vis de l'instant de référence du signal d'application dans le but de modifier le spectre fréquentiel du signal de commande. On peut, par exemple, utiliser une loi de probabilité uniforme, c'est-à-dire que l'écart temporel est distribué sur un intervalle compris entre deux valeurs extrémales, par exemple [-dTmax, dTmax]. Chacune des valeurs de cet intervalle a une probabilité égale d'être appliquée. Dans un tel cas, on peut montrer que le coefficient d'atténuation du spectre du signal de commande initial, est donné aux fréquences multiples de la fréquence d'échantillonnage par l'expression : sin(wdTrax ) où w=2Jrk/T wdTäa,, En jouant sur la valeur de dTmax, on peut plus ou moins atténuer le spectre du signal de commande. Plus on augmente la valeur de dTmax, plus on distribue la probabilité de l'instant d'application de l'incrément, et plus l'atténuation des pics d'énergie du signal d'origine aux fréquences multiples de la fréquence d'échantillonnage augmente. Un autre exemple de loi de probabilité utilisable est l'exemple de la loi de probabilité binomiale dans laquelle la valeur -dTmax a une probabilité d'être appliquée égale à p, tandis que la valeur d'écart temporel +dTmax a une probabilité d'être appliquée égale à 1-P. Dans un tel cas, on utilise en général P=1/2, c'est-à-dire que le retard négatif -dTmax et le retard positif + dTmax ont une probabilité égale d'être appliqués. Dans ce cas, le coefficient d'atténuation est donné par : cos(w dT nax ) où w = 2n' k / T En particulier, dans le cas où dTmax vaut T/4, les 20 harmoniques impaires du signal d'origine disparaissent totalement. D'autres lois de probabilité peuvent être utilisées afin d'établir la commande. Par ailleurs, on rappelle qu'il est possible d'utiliser des séquences de nombres 25 pseudo aléatoires générées au préalable. Cette génération peut être effectuée avant le lancement d'un satellite 1 et inscrite dans la mémoire 9, ou en cours de fonctionnement grâce à un générateur de nombres aléatoires qui génère une séquence dont les propriétés spectrales permettent 30 d'atténuer des harmoniques proches d'une ou de plusieurs fréquences de résonance. Par ailleurs, il convient de noter que, dans le cas de la commande de tuyères 3, il est possible de modifier la durée des impulsions de telle sorte qu'on puisse assurer, en plus de la consigne de déplacement souhaitée, un contrôle de pointage ou de positionnement du satellite 1 simultanément. De même, dans le cas de l'orientation des panneaux solaires 2, il est possible de modifier l'incrément au cours du temps de telle sorte qu'on puisse, par exemple, piloter un satellite à voile solaire. REVENDICATIONS 1. Procédé de commande d'au moins un système mécanique (1) présentant au moins un élément de structure flexible (2, 5) et au moins un actionneur ou groupe d'actionneurs (3, 4), ledit procédé étant destiné à réduire le niveau de vibrations dudit élément de structure flexible (2, 5) tout en commandant ledit actionneur ou groupe d'actionneurs (3, 4) de sorte à atteindre au moins un objectif assigné audit système mécanique (1), tel que réaliser un déplacement, ladite commande étant constituée d'incréments (23, 33) devant être réalisés à une période T dite période d'échantillonnage, chacun desdits incréments (23, 33) étant constant en amplitude dans un intervalle de temps dit intervalle d'échantillonnage (22, 32) s'étendant entre deux instants d'échantillonnage successifs et chacun desdits incréments (23, 33) s'appliquant pendant un intervalle d'application d'incrément (24) d'une durée (AT) inférieure ou égale à ladite période d'échantillonnage T, ledit procédé étant caractérisé en ce que, dans chaque intervalle d'échantillonnage (22, 32), l'un au moins des instants initial (25, 34) et final (26) d'application dudit incrément (23, 33) est modifié par l'ajout d'un écart temporel (dT) variable, aléatoire ou pseudo-aléatoire, dit écart aléatoire (28, 35). 2. Procédé selon la revendication 1, caractérisé en ce que la durée d'application desdits incréments (33) est égale à la période d'échantillonnage T, et un écart aléatoire (35) est ajouté à l'instant initial (34) d'application desdits incréments (33). 3. Procédé selon la revendication 1, caractérisé en ce que la durée d'application (AT) desdits incréments (23) est strictement inférieure à la périoded'échantillonnage T, et un écart aléatoire (28) est ajouté soit à l'instant initial (25), soit à l'instant final (26), soit aux deux instants initial (25) et final (26) d'application desdits incréments (23), avec des valeurs pouvant être différentes dans ce dernier cas. 4. Procédé selon la revendication 3, caractérisé en ce que dans chaque intervalle d'échantillonnage (22), le même écart aléatoire (28) est ajouté à la fois à l'instant initial (25) et à l'instant final (26) d'application desdits incréments (23). 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'une loi de probabilité choisie pour déterminer ledit écart aléatoire (28, 35) est définie de sorte à modifier au moins un des pics d'énergie de la commande. 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que la valeur absolue maximum de l'écart aléatoire (28, 35) est inférieure à une valeur fixée, de préférence une fraction entière de la période d'échantillonnage T. 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que l'écart aléatoire (28, 35) est à moyenne nulle. 8. Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce que les valeurs de l'écart aléatoire (28, 35) sont sélectionnées parmi les multiples entiers d'une valeur pré-définie T' 9. Procédé selon la revendication 8, caractérisé en ce que ladite valeur pré-définie T' est un 30 diviseur entier de la période d'échantillonnage T 10. Procédé selon l'une quelconque des revendications 8 et 9, caractérisé en ce que la fréquence 1/T' est choisie de sorte à limiter les éventuellesrésonances avec les modes de vibrations dudit élément de structure flexible (25), et de préférence ne coïncide pas avec les fréquences propres desdits modes de vibrations. 11. Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce qu'une loi de probabilité choisie pour déterminer ledit écart aléatoire (28, 35) est entièrement définie par un nombre fini de paramètres calculés de sorte à modifier au moins un des pics d'énergie de la commande. 12. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce qu'une loi de probabilité choisie pour déterminer ledit écart aléatoire (28, 35) est une loi de probabilité discrète à m valeurs, lesdites M valeurs et leurs probabilités étant calculées de sorte à modifier au moins un des pics d'énergie de la commande. 13. Procédé selon la revendication 12, caractérisé en ce que ladite loi de probabilité est une loi binomiale prenant des valeurs extrêmes d'écart aléatoire (28, 35) +dTmax avec une probabilité de 1/2, et -dTmax avec une probabilité de 1/2. 14. Procédé selon la revendication 12, caractérisé en ce que ladite loi de probabilité est une loi discrète à m valeurs équiprobables en progression arithmétique entre une valeur minimale -dTmax et une valeur maximale +dTmax des écarts aléatoires (28, 35), avec dTmax = T.(M-1)/2M. 15. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce qu'une loi de probabilité choisie pour déterminer l'écart aléatoire (28, 35) est une loi uniforme entre des valeurs extrêmes d'écart aléatoire dTmin et dTmax, les valeurs dTmin et dTmax étant choisies pour modifier au moins un pic d'énergie de la commande. 16. Procédé selon la revendication 15, caractérisé en ce que ladite loi de probabilité est une loi uniforme entre des valeurs maximales -dTmax et dTmax. 17. Procédé selon la revendication 16 telle que rattachée à la revendication 4, caractérisé en ce que dTmax = T/2, de sorte que le spectre de la commande ne comporte pas de raie à la fréquence d'échantillonnage 1/T. 18. Procédé selon la revendication 16 telle que rattachée à la revendication 4, caractérisé en ce que dTmax = T/4, de sorte que le spectre de la commande ne comporte pas de raie aux fréquences multiples impaires de la fréquence d'échantillonnage 1/T. 19. Procédé selon l'une quelconque des revendications 1 à 18, caractérisé en ce que les écarts aléatoires (28, 35) prennent les valeurs de nombres successifs d'au moins une séquence aléatoire ou pseudoaléatoire de nombres produite à chaque période d'échantillonnage par au moins un programme informatique générateur de nombres aléatoires ou pseudo-aléatoires 20. Procédé selon l'une quelconque des revendications 1 à 18, caractérisé en ce que les écarts aléatoires (28, 35) prennent les valeurs de nombres successifs d'une séquence longue mais finie de nombres présentant des caractéristiques aléatoires proches de caractéristiques désirées, ladite séquence étant pré-calculée et appliquée suivant le procédé, puis reproduite en boucle une fois qu'elle est achevée. 21. Produit programme d'ordinateur comprenant des instructions de code adaptées à la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 20, lorsqu'il est chargé et exécuté sur des moyens informatiques. 22. Dispositif de mise en oeuvre du procédé selonl'une quelconque des revendications 1 à 20, caractérisé en ce qu'il comprend des moyens de calcul et informatiques permettant l'exécution du produit programme selon la revendication 21, et des moyens électriques et/ou électroniques permettant de transformer des ordres de commande numériques élaborés par les moyens de calcul et informatiques en signaux électriques envoyés à au moins un actionneur (3, 4), ledit dispositif étant tel que - au moins un programme informatique dudit produit programme est mis en oeuvre par lesdits moyens de calcul et informatiques pour produire lesdits incréments (23, 33) de ladite commande d'au moins un desdits actionneur et groupe d'actionneurs (3, 4) dudit système mécanique (1) de sorte à atteindre au moins un objectif assigné audit système mécanique (1), tel que réaliser un déplacement ; - ledit programme informatique est mis en oeuvre pour produire, à chaque période d'échantillonnage T, ledit écart aléatoire (28, 35) suivant une loi de probabilité se rapprochant au mieux d'une loi de probabilité désirée ; - ledit programme informatique est mis en oeuvre pour modifier le ou les instants (25, 26, 34) d'application dudit incrément (23, 33) en lui ajoutant ledit écart aléatoire (28, 35), selon le procédé de l'une quelconque des revendications 1 à 20, et produire ainsi la commande numérique à appliquer effectivement audit système mécanique (1) ; et - lesdits moyens électriques et/ou électroniques transforment ladite commande numérique produite par ledit programme informatique en signaux électriques permettant à au moins un desdits actionneur et groupe d'actionneurs (3, 4) de réaliser physiquement la commande assignée. 23. Dispositif selon la revendication 22 pour la mise en oeuvre du procédé selon la revendication 20, caractérisé en ce qu'une séquence longue mais finied'écarts aléatoires (28, 35) ayant les caractéristiques aléatoires désirées est calculée préalablement, puis stockée dans une mémoire (9) desdits moyens de calcul et informatiques (8), lesdits écarts aléatoires (28, 35) étant, successivement et suivant un ordre établi, ajoutés à au moins un instant d'application (25, 26, 34) desdits incréments (23, 33) de ladite commande, ladite séquence étant reproduite en boucle une fois qu'elle est achevée. 24. Dispositif selon la revendication 22, caractérisé en ce qu'au moins une partie des paramètres intervenant dans le calcul des écarts aléatoires (28, 35) par ledit programme informatique, tels que les paramètres d'une loi de probabilité, ou une séquence de valeurs aléatoires, peut être modifiée au cours du fonctionnement dudit système mécanique (1), afin notamment d'optimiser l'efficacité du procédé. 25. Dispositif selon la revendication 24, caractérisé en ce qu'il comprend au moins un capteur (10) dédié de mesure des vibrations d'au moins un élément de structure flexible (2, 5) présent dans le système mécanique (1), des moyens de calcul du spectre fréquentiel de ces vibrations, et des moyens de modification de paramètres intervenant dans le calcul desdits écarts aléatoires (28, 35) afin de réduire le niveau desdites vibrations en modifiant au moins un des pics d'énergie du spectre fréquentiel de la commande. 26. Application du procédé selon l'une quelconque des revendications 1 à 20, caractérisée en ce qu'elle consiste à commander l'ouverture d'au moins une tuyère (3) d'un système de contrôle d'orbite, de stabilisation et/ou d'orientation d'un engin (1) se déplaçant dans l'espace, tels que lanceur, satellite artificiel, sonde interplanétaire, ledit engin (1) comportant au moins un élément de structure flexible, telde grands appendices, antenne (5), panneau solaire (2), et structure primaire. 27. Application selon la revendication 26, caractérisée en ce qu'elle consiste à - commander l'ouverture d'au moins une tuyère (3) afin d'imprimer audit engin (1) un incrément de vitesse permettant de corriger la trajectoire dudit engin (1) ; - réaliser nominalement ledit incrément de vitesse en une suite de petits incréments effectués périodiquement 10 à une période T ; - réaliser lesdits petits incréments par la commande (7) d'ouverture d'au moins une tuyère (3), ladite commande étant conventionnellement égale à 0 pour un ordre de fermeture de tuyère (3), et à 1 pour un ordre 15 d'ouverture de tuyère (3) ; - avant la mise en oeuvre dudit procédé, calculer une commande nominale de façon à être périodique de période T dans l'intervalle de temps [Ti, Ti + N.T], où N est un nombre entier, de réalisation d'une impulsion totale, et, 20 dans ledit intervalle de temps, égale à un incrément de 1 sur tous les intervalles de temps [Ti + k.T + Ti, Ti + k.T + T2] où k est un nombre entier inférieur à N, Tl et T2 étant des constantes telles que 0 28. Application selon la revendication 27, caractérisée en ce qu'en plus desdits écarts aléatoires (28, 35), on ajoute à un seul des instants initial (25) et final (26) d'application de l'incrément un autre écart temporel, calculé de manière connue pour réaliser un contrôle d'attitude ou de pointage de l'engin spatial (1) en même temps que l'on réalise la commande d'incrément de vitesse pour le contrôle d'orbite. 29. Application du procédé selon l'une quelconque des revendications 1 à 20, caractérisée en ce qu'elle consiste à commander la rotation d'au moins un moteur pas à pas (4) pour le pointage d'au moins une charge utile d'un engin spatial (1), telle qu'au moins un panneau solaire (2), une antenne (5), un instrument, ladite charge utile présentant au moins un élément de structure flexible (2). 30. Application selon la revendication 29, caractérisée en ce que - une commande î de l'angle entre le stator et le rotor dudit moteur (4) est échantillonnée à une période d'échantillonnage T, et s'incrémente d'une certaine valeur à chaque période d'échantillonnage, - une commande nominale î se présente sous la forme : OC k = k_1 + (d U) k Oë a k représente la valeurconstante de la commande entre les instants k.T et (k+l).T, où k est un nombre entier, et (da) k est la valeur de l'incrément à réaliser à la date k.T, les incréments s'effectuant par convention aux instants d'échantillonnage (34), ladite commande est modifiée selon le procédé de l'une quelconque des revendications 1 à 20 en ajoutant, aux instants k.T de réalisation (34) des incréments, un écart temporel aléatoire (35) ayant des caractéristiques désirées dans le but de minimiser des résonances avec des modes de vibration dudit élément de structure flexible (2). 31. Application selon la revendication 30, caractérisée en ce que les incréments (dia) k sont tous 15 égaux, comme dans le cas du maintien en rotation à vitesse constante d'un panneau solaire (2). 32. Application selon la revendication 30, caractérisée en ce que les incréments (d OC) k nominalement égaux sont légèrement modifiés, de manière connue, pour 20 réaliser le contrôle d'au moins un phénomène physique résultant de la rotation du moteur (4), comme le pilotage d'un satellite (1) à voile solaire dans le cas de la commande d'un panneau solaire (2) muni d'au moins un volet. 33. Application selon l'une quelconque des 25 revendications 31 et 32, caractérisée en ce que le spectre fréquentiel de la commande nominale à appliquer au moteur (4) se situe principalement à des fréquences très inférieures à la fréquence d'échantillonnage, en particulier pour réaliser un mouvement de balayage à faible 30 vitesse.