Brevet FR3000246A1 — publié le 20140627

Générateur de séquences chaotiques La présente invention concerne un générateur de séquences chaotiques de valeurs entières et un procédé de génération de séquences chaotiques associé. Plus particulièrement, l'invention se rapporte au domaine de la sécurité des données partagées, transmises et stockées dans les réseaux de transmission d'informations. Le transfert des données confidentielles (documents d'entreprise, informations médicales, résultats de recherche, informations personnelles de types photos et vidéos, etc.) dans un environnement ouvert utilisant les canaux usuels de communication (câbles, Internet, Radio-mobiles, satellites,...) doit être réalisé avec une sécurité maximale et un débit suffisant. A cet effet, les crypto-systèmes basés sur les signaux chaotiques sont adéquats pour atteindre les objectifs cités. Un élément déterminant dans tout cryptosystème basé chaos est le générateur des séquences chaotiques, qui sert à la génération des clés secrètes et au processus de chiffrement/déchiffrement des données dans les opérations de substitution et de permutation. La confidentialité des données dépendra entre autres du degré du chaos (c'est-à-dire de l'aléatoire) des séquences produites par le générateur des séquences chaotiques utilisé. Les séquences chaotiques utilisables dans les crypto-systèmes sont caractérisées par des propriétés de longueur de période et de diffusion. La demande de brevet W02011/121218 de Safwan El Assad et Hassan Noura, présente un générateur de séquences chaotiques qui permet de générer des séquences chaotiques de longueur très grande, donc de périodicité négligeable, utilisables dans des applications cryptographiques concrètes avec un grand niveau de sécurité. La présente invention propose un générateur qui améliore la propriété de « confusion-diffusion » des générateurs de séquences chaotiques, tout en ayant une longueur de période suffisamment grande pour une application cryptographique avec une bonne sécurité. Selon un premier aspect, l'invention concerne un générateur de séquences chaotiques de nombres de valeurs entières représentées sur un nombre de bits prédéterminé, lesdites séquences étant destinées à être utilisées notamment dans des systèmes de cryptage d'informations à base de clé, ledit générateur comportant un nombre m de moyens de mise en oeuvre de fonctions non linéaires. Le générateur comporte, en sortie desdits moyens de mise en oeuvre de fonctions non linéaires, des moyens de combinaison aptes à combiner les sorties desdits moyens de mise en oeuvre de fonctions non linéaires par application d'une matrice de diffusion binaire donnée, permettant d'obtenir un nombre m de sorties, chaque valeur de sortie étant calculée par une combinaison binaire de sorties desdits moyens de mise en oeuvre de fonctions non linéaires. Avantageusement, le générateur proposé améliore la caractéristique de « confusion-diffusion » des séquences chaotiques générées, tout en préservant une bonne propriété de périodicité. De plus, un générateur selon l'invention est aisé à mettre en oeuvre et permet la génération de séquences chaotiques en un temps de calcul assez rapide comparé aux générateurs standard proposés par le NIST (« National lnstitute of Standards and Technology »). Le générateur de séquences chaotiques selon l'invention peut présenter une ou plusieurs des caractéristiques ci-dessous, prises indépendamment ou en combinaison : - il comporte en outre des moyens de perturbation raccordés en sortie des moyens de combinaison ; - il comporte m moyens de perturbation, un dit moyen de perturbation étant raccordé à chaque sortie des moyens de combinaison ; - un dit moyen de perturbation comprend un registre à décalage à réaction ; - lesdits moyens de mise en oeuvre de fonctions non linéaires comportent des moyens de mise en oeuvre de cartes chaotiques d'un premier type et des moyens de mise en oeuvre de cartes chaotiques d'un deuxième type ; - les moyens de mise en oeuvre de cartes chaotiques de premier type et les moyens de mise en oeuvre de cartes chaotiques de deuxième type sont alternés ; - les cartes chaotiques de premier type sont des cartes chaotiques de type « skew tent » et les cartes chaotiques de deuxième type sont des cartes chaotiques cartes chaotiques linéaires par morceaux PWLCM ; - les nombre d'entrées et de sorties des moyens de combinaison est égal à l'un des nombres suivants : 4, 8, 32, 64. - les moyens de combinaison comportent, pour l'obtention d'une sortie, un nombre m d'interrupteurs branchés en sortie des moyens de mise en oeuvre de fonctions non linéaires, un interrupteur fermé correspondant à la présence d'un élément égal à un dans la matrice de diffusion binaire et un interrupteur ouvert correspondant à la présence d'un élément égal à zéro dans la matrice de diffusion binaire, les sorties desdits interrupteurs étant raccordées à une porte « ou exclusif ». Selon un deuxième aspect, l'invention propose un procédé de génération de séquences chaotiques de nombres de valeurs entières représentées sur un nombre de bits prédéterminé, lesdites séquences étant destinées à être utilisées notamment dans des systèmes de cryptage d'informations à base de clé, le procédé étant caractérisé en ce qu'il comporte les étapes de : - application, de m fonctions non linéaires sur m valeurs initiales données, permettant d'obtenir m valeurs chaotiques, - combinaison des m valeurs chaotiques par application d'une matrice de diffusion binaire, permettant d'obtenir m valeurs chaotiques combinées et -application d'une perturbation sur au moins une sous-partie des m valeurs chaotiques combinées. D'autres caractéristiques et avantages de l'invention ressortiront de la description qui en est donnée ci-dessous, à titre indicatif et nullement limitatif, en référence aux figures annexées, parmi lesquelles : -la figure 1 est un synoptique d'un générateur de séquences chaotiques selon l'invention ; - la figure 2 est un synoptique illustrant une mise en oeuvre des moyens de combinaison pour la génération d'une sortie combinée, et - la figure 3 est un organigramme illustrant les principales étapes d'un mode de réalisation d'un procédé de génération de séquences chaotiques selon l'invention. La figure 1 illustre un générateur de séquences chaotiques 2 selon l'invention, permettant de générer m séquences de valeurs entières Xt(n), 1 i in, où n est un entier et chaque valeur X (n) est représentée sur un nombre N de bits prédéterminé, par exemple N=32, ce qui permet d'obtenir une précision finie de calcul, quelque soit la plateforme matérielle sur laquelle le générateur est implémenté. De préférence, N est une puissance de deux. Les m séquences générées sont combinables, si nécessaire, en une seule séquence de nombre chaotiques S {X1(1),..., Xm (1), Xi(2),..., Les séquences générées sont utilisables dans divers crypto-systèmes utilisant des nombres pseudo-aléatoires, notamment pour des processus de chiffrement/déchiffrement appliquées pour sécuriser des données partagées, transmises et stockées, mais également pour des applications de stéganographie, de tatouage numérique et de génération de clés secrètes. Le générateur 2 de nombres chaotiques prend en entrée initialement un ensemble de m conditions initiales Xi (0), X2 Xm (0), définissant les valeurs initiales ou graines pour n=0. Le générateur 2 comporte m moyens de mise en oeuvre d'une fonction non- linéaire notés 4, 6, 8, 10. Les moyens de mise en oeuvre d'une fonction non-linéaire comprennent des circuits mettant en oeuvre la fonction xlog(x) , ou xexp[cos(x)], ou une carte de Chebyshev, ou une carte « skew tent » ou encore une carte chaotique linéaire par morceaux ou PWLCM (« piecewise linear chaotic map »). Les fonctions non-linéaires () associées à chaque entrée X (n) sont détaillées ci-dessous dans un mode de réalisation préféré. De préférence, les moyens de mise en oeuvre de fonctions non linéaires de différents types sont utilisés, afin de maximiser la propriété de diffusion de la séquence chaotique générée. En particulier, selon un mode de réalisation préféré, les moyens de mise en oeuvre de fonction non-linéaire sont à base de cartes chaotiques de deux types alternées d'une entrée à l'entrée suivante. Il a été constaté qu'une alternance entre des cartes chaotiques « skew tent » de fonction non linéaire F10 et des cartes chaotiques de type PWLCM de fonction non linéaire F2 () donne de très bons résultats. Sans perte de généralité, la fonction F1 () est appliquée pour les X (n -1) où i est impair, et la fonction F2 0 est appliquée pour les X (n -1) où i est pair. La fonction non linéaire F1 () basée sur une carte chaotique de type « skew tent », appelée carte chaotique de premier type, est la suivante : Fi (X (n -1)) = 2N x X(n -1) si() X(n-1" P N 2N - X (n -1) 2N - P +1 siP En variante, une autre alternance des cartes chaotiques de type « skew tent » et de type PWLCM est utilisée. Ainsi, chaque fonction non linéaire utilisée est définie, dans ce mode de réalisation, par son type et par un paramètre P associé. Selon encore une autre variante, plus de deux types de cartes chaotiques sont utilisés, par exemple trois ou quatre cartes chaotiques alternées successivement sur les entrées du générateur 2 de séquences chaotiques. Les sorties X' , (n) = F, (X t(n -1)) après mise en oeuvre des fonctions non- linéaires choisies sont fournies en entrée des moyens de combinaison 20, aptes à appliquer une matrice de diffusion binaire de dimension mxm. Chaque sortie Y(n)desmoyens de combinaison 20 est calculée à partir des X (n), j = 1,...,m et de la matrice de diffusion D.1 comme suit : 17,(n) = D ,j (n) , j = 1, . , m , où le symbole CD désigne l'opération de combinaison décrite plus en détail ci-après. Les sorties des moyens de combinaison 20, en même nombre que les entrées (m) sont perturbées par des moyens de génération de séquences de perturbation notés 22, 24, 26 et 28, qui sont des registres à décalage à réaction (« linear feedback shift reg ister »). Dans ce mode de réalisation, chaque sortie i des moyens de combinaison 20 est perturbée par addition, au moyen de fonctions « ou exclusif » 30, 32, 34, 36, d'une séquence binaire de perturbation fournie par chaque registre à décalage 22, 24, 26, 28, permettant d'obtenir une valeur chaotique X , (n) . Les registres à décalage 22, 24, 26, 28 sont aptes à mettre en oeuvre m polynômes primitifs de degrés respectifs k,, i=1,...,M. Chaque registre à décalage génère une séquence de perturbation Q(n) représentée sur k, bits. Chaque registre à décalage est caractérisé par une bonne fonction d'auto- corrélation, une distribution presque uniforme, un cycle de longueur maximale égale à 24 -1 et une implémentation logicielle ou matérielle facile. En variante, une perturbation est appliquée seulement à un sous-ensemble des sorties Y(n) des moyens de combinaison 20. Les valeurs de sortie X , (n) sont fournies à nouveau sur les entrées des moyens de mise en oeuvre de fonctions non-linéaires, afin de générer une nouvelle série de m valeurs chaotiques. La figure 2 illustre une mise en oeuvre des moyens de combinaison permettant de générer une valeur de sortie Y, (n) . Comme illustré à la figure 2, la génération d'une valeur de sortie Y, (n) est obtenue potentiellement à partir de l'ensemble des valeurs d'entrée Les références des éléments communs avec ceux de la figure 1 sont reprises. Ainsi, le générateur 2 applique d'abord les moyens de mise en oeuvre de fonctions non-linéaires à chacune des valeurs d'entrée X 1(n -1),..., X ,,,(n -1) . Les sorties X' (n) sont fournies à des interrupteurs 40, 42, 44, 46, dont la position ouverte ou fermée est fonction de la valeur de la matrice de diffusion binaire mise en oeuvre par les moyens de combinaison 20. Par exemple, dans le cas où m=4, la matrice de diffusion appliquée est par exemple la matrice suivante : 1 1 1 0 1 1 0 1 MD4 1 0 1 1 0 1 1 1 Pour le calcul appliqué à X'3(n), les interrupteurs 40, 44 et 46 sont fermés, et l'interrupteur 42 est ouvert. De manière plus générale, la matrice de diffusion MDm de taille mxm étant composée d'éléments un interrupteur appliqué à la sortie des moyens de mise en oeuvre de la fonction non-linéaire pour l'entrée X j (n -1) est ouvert si D,,J=0 et fermé si D,,J=1. Une sous-partie déterminée par les valeurs Di,j de la matrice de diffusion est combinée par un « ou exclusif » 48 appliqué bit à bit, permettant d'obtenir une valeur combinée Y1(n) représentée sur N bits. Ainsi, il apparaît clairement que les moyens de combinaison 20 sont faciles à implémenter par des interrupteurs et des portes « ou exclusif », ce qui permet de générer des calculs rapides. Ainsi, en sortie, les moyens de combinaison 20 génèrent des valeurs d'échantillons combinées : Y1(n) = D.1 0 X' j (n), j =1,..., m Où D, est un élément de la matrice de diffusion binaire comme expliqué ci-dessus, X:(n) est la sortie des moyens de mise en oeuvre de fonction non linéaire Ft() appliqués, et le symbole 0 représente l'opération de combinaison. De préférence, comme illustré dans l'exemple de matrice de diffusion donné ci- dessus pour m=4, la matrice de diffusion MDm est une matrice telle que MDm-1 = MDm , ce qui facilite l'implémentation de l'opération de déchiffrement dans une application d'un générateur de séquences chaotiques dans un crypto-système à base de clé dynamique. Dans un mode de réalisation, m=8 et la matrice de diffusion MD8 suivante est de préférence utilisée : 1 0 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 1 1 0 1 1 0 1 0 1 1 1 1 1 1 0 MD8 = 1 1 0 0 0 1 1 1 0 1 1 0 1 0 1 1 0 0 1 1 1 1 0 1 1 0 0 1 1 1 1 0 Dans un autre mode de réalisation, m=16 et la matrice de diffusion MD16 suivante est de préférence utilisée : 0 0 0 1 1 0 1 0 1 1 0 0 0 1 1 0 0 0 1 0 0 1 0 1 1 1 0 0 1 0 0 1 0 1 0 0 1 0 1 0 0 0 1 1 1 0 0 1 1 0 0 0 0 1 0 1 0 0 1 1 0 1 1 0 1 0 1 0 0 1 0 0 1 0 0 1 0 0 1 1 0 1 0 1 1 0 0 0 0 1 1 0 0 0 1 1 1 0 1 0 0 0 0 1 0 1 1 0 1 1 0 0 0 1 0 1 0 0 1 0 1 0 0 1 1 1 0 0 MD16 1 1 0 0 1 0 0 1 0 0 1 0 0 1 0 1 1 1 0 0 0 1 1 0 0 0 0 1 1 0 1 0 0 0 1 1 0 1 1 0 1 0 0 0 0 1 0 1 0 0 1 1 1 0 0 1 0 1 0 0 1 0 1 0 0 1 1 0 0 0 1 1 0 1 0 1 1 0 0 0 1 0 0 1 0 0 1 1 1 0 1 0 0 1 0 0 1 0 0 1 1 1 0 0 0 1 0 1 0 0 1 0 0 1 1 0 1 1 0 0 1 0 1 0 0 0 0 1 Il est à noter que pour une application par calcul algorithmique de la matrice de diffusion binaire MD16 ci-dessus, il est possible de réaliser des calculs intermédiaires de termes communs à plusieurs sorties. Ce regroupement permet également de diminuer le nombre de portes « ou exclusif » à mettre en oeuvre dans une implémentation matérielle. Cependant, une implémentation matérielle simple, à base de 96 portes « ou exclusif » est suffisamment rapide pour la mise en oeuvre des moyens de combinaison 20. Dans un autre mode de réalisation, m=32 et la matrice de diffusion MD32 proposée est construite en appliquant la démarche décrite par Koo et al. dans l'article « On constructing of a 32x32 binary matrix as a diffusion layer for a 256-bit block cipher », publié dans Proceedings ICISC 2006, édition Springer Verlag LNCS 4296, pages 51-64. Dans cet article, il est proposé de générer une matrice de diffusion pour des applications cryptographiques à partir de l'équation suivante : MD32 = Mod(Li xMm x4,2) Où ri = Les matrices L et M sont définies comme suit : B1,1 B1,2 B1,8 B2,1 B2,2 B2,8 . . B8 B8,2 - - - B 4x4 si h1 = 1 4x4 si h1 0, hi,1 e H1 B.= où Avec /4x4 la matrice identité de taille 4x4, 04x4 la matrice zéro de taille 4x4, et : H1 La matrice M, entrant dans la construction de la matrice de diffusion binaire est : Mm= H1 °8x8 °8x8 °8x8 °8x8 H2 °8x8 °8x8 °8x8 °8x8 H3 °8x8 Hm4 _ _°8x8 °8x8 °8x8 Où 08,8 est la matrice zéro de taille 8x8 (matrice ne contenant que des 0) et les matrices Elm, à Hm4 sont les suivantes : 0 1 1 1 1 1 1 1 1 0 1 1 0 0 0 0 1 0 1 1 1 1 1 1 0 1 1 1 0 0 0 0 1 1 0 1 1 1 1 1 1 1 1 0 0 0 0 0 1 1 1 0 1 1 1 1 1 1 0 1 0 0 0 0 1 1 1 1 0 1 1 1 ; H',2 0 0 0 0 1 0 1 1 1 1 1 1 1 0 1 1 0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 1 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 1 1 0 1 1 1 0 1 0 0 0 0 1 1 1 1 1 1 1 0 1 1 1 0 0 0 0 0 1 1 1 1 1 1 0 1 0 1 1 1 0 0 0 0 1 1 1 1 1 0 1 1 1 0 1 1 0 0 0 0 1 1 1 1 0 1 1 1 H,,,3 0 0 0 0 1 1 0 1 ; H',4 1 1 1 0 1 1 1 1 0 0 0 0 1 1 1 0 1 1 0 1 1 1 1 1 0 0 0 0 0 1 1 1 1 0 1 1 1 1 1 1 0 0 0 0 1 0 1 1 0 1 1 1 1 1 1 1 Comme déjà expliqué ci-dessus, quelque soit le choix de m parmi les valeurs 4, 8, 16, 32, les m sorties des moyens de combinaison 20 sont perturbées via des séquences de perturbation générées par des registres à décalage. Un registre à décalage 26 est représenté sur la figure 2. La séquence de perturbation est appliquée à la sortie Y, (n) par la porte « ou exclusif » 34. Les sorties Y, (n) pour i=1,...,m sont perturbées à l'instant initial et tous les A, , où A, est la valeur de l'orbite chaotique sans perturbation. La mesure de l'orbite chaotique est effectuée de préférence selon la technique décrite dans la demande de brevet W02011/121218. Alternativement, d'autres techniques connues de mesure d'orbite chaotique peuvent être utilisées. La perturbation appliquée peut donc s'écrire comme suit, pour chaque i, de 1 à m : X i (n) Y (n) si mod(n, Ai ) #0 si mod(n, Ai ) =0 ( ( Y, (n) CD mod Q. -n+ 1 ,24 A ) ) Où le symbole 0 représente l'opération « ou exclusif ». L'équation ci-dessus montre que la perturbation est appliquée sur les k, derniers bits de poids faibles de Y(n) aux instants n = lxA, ,1 = 0,1,..... Pour la valeur m=32, les polynômes g suivants sont mis en oeuvre, un polynôme g étant mis en oeuvre par le registre à décalage mis en oeuvre pour perturber la séquence (n) : g (x) = X15 + X13 + Xi° + X1 + 1 g 2(x) =x17 +x3 +x2 + +1 g 3(x) = x19 + x5 + x2 + + 1 g 4 (x) =x23 +x12 +x5 +x4 +1 g 5 (x) =x21 + x 2 +1 g 6(x) = x17 +x7 +x4 +x3 +1 g 7 (x) = x15 + x9 + x4 + + 1 g 8 (x) = x19 + x9 + x8 + x7 + x6 + x3 + 1 g 9 (x) + x21 +x14 +x7 +x2 +1 g10 (x) = x15 + x7 + x4 + x1 + 1 g11 (x) = x17 16 3 1 g12 (x) = x21 + X + X + X +1 +x13+ X5 + X2 +1 g13 =X 15 + x14 + x12 + x2 +1 g14 = X15 + X13 + X111 + X6 + 1 + x23 + x5 + x4 + x1 + 1 g15 17 12 6 3 2 g16 =X +X +X +X +X +X+ 1 g17 = X15 + X13 + X6 + X6 + 1 g18 x 21 + x14 + x7 + x6 + x3 + x 2 + 1 g19 + x21 + x10 + x6 + x4 + x3 + x2 + 1 g20 = X23 + X12 + X" + X5 + g21 + x15 + x14 + x9 + x2 + 1 g22 + x17 + x9 + x8 + x6 + x4 + x1 + 1 g23 =x15 +x13 +x12 + x10 +1 g24 + x21 + x8 + x7 + x4 + x3 + x 2 + g25 x15 + x12 + x3 + x1 + 1 g26 = x17 + x' +x7 +x6 +x4 +x3 +1 g22 = x19 +X13 +x' +X5 +x4 +x3 +1 = x15 + x13 + x7 + x4 +1 g29 x21 +X14 +X12 +X7 +X6 +X4 +X3 +x2 +1 g30 x19 + x12 + x10 + x9 + x7 + x3 + 1 g31 x17 + x16 + x14 + x10 + x3 + x2 + 1 g32 = x15 +x1 +1 Dans les cas où m=4,8 ou 16, un nombre m de polynômes primitifs est sélectionné dans la liste ci-dessus, cette sélection pouvant être effectuée de façon arbitraire. Selon une variante non représentée, afin d'améliorer encore les performances du générateur chaotique selon l'invention, il est possible d'appliquer la technique dite de sous échantillonnage, qui consiste à itérer l'application des moyens de mise en oeuvre de fonction non linéaire s fois, avec s un entier positif supérieur à 3. Néanmoins, le temps calculatoire est augmenté si cette variante est appliquée. Le générateur de séquences chaotiques, décrit ci-dessus selon différents modes de réalisation, permet de générer des séquences chaotiques ayant de bonnes performances cryptographiques, notamment des probabilités de bits égaux à 0 et 1 très proches de 0,5 pour toutes les sorties X (n), i = 1,...,m. Comme expliqué ci-dessus, dans le mode de réalisation préféré, une alternance des types cartes chaotiques définissant les fonctions non linéaires mises en oeuvre par les moyens de mise en oeuvre de fonctions est appliquée. Ainsi, un générateur 2 de séquence chaotiques est entièrement défini par une clé secrète qui comporte : m conditions initiales sur N bits, m/2xN bits pour représenter les paramètres des cartes chaotiques de premier type, m/2x(N-1) bits pour représenter les paramètres des cartes chaotiques de deuxième type, m conditions initiales de registres à décalage ce qui correspond à minimum mx15bits. Ainsi, pour ce mode de réalisation, les tailles de clé secrète sont les suivantes : TK4=314 bits pour m=4, TK8=628 bits pour m=8, TK18=1256 bits pour m=16 et TK32=2512 bits pour m=32. Comme on peut le constater, les tailles de clés sont grandes, même pour un nombre d'entrées égal à 4. La figure 3 illustre les principales étapes d'un procédé de génération de séquences chaotiques de valeurs entières selon un mode de réalisation de l'invention. Un tel procédé peut être implémenté par un code logiciel, exécuté par une unité central d'un processeur, par exemple d'un ordinateur. En entrée 50 du procédé sont fournies les m conditions initiales, qui sont les valeurs initiales X1(0),. ..,X,(0), chaque valeur étant codée sur N bits. De plus, les paramètres permettant de définir les fonctions non linéaires et les perturbations à appliquer sont également fournis en entrée du procédé, ainsi que le nombre Ns d'échantillons à générer pour chacune des m sorties. Par exemple, une clé secrète de format déterminé contient toutes les conditions initiales et tous les paramètres dans un ordre prédéfini. Une première étape d'initialisation 52 consiste à initialiser une variable n entière à 1. Ensuite, une fonction non linéaire Fi() est appliquée à X,(n-1) à l'étape 54 de génération de nombres chaotiques, pour i=1,...,m. Comme expliqué ci-dessus, les fonctions non linéaires sont de différents types selon les indices i, par exemple une fonction de premier type pour les indices impairs et une fonction de deuxième type pour les indices pairs. De manière optionnelle, chaque étape de génération 54 est itérée un nombre s de fois à l'étape 56. Un ensemble de m valeurs chaotiques {X'1(n),..., X','(n)} est obtenu à l'issue de l'étape 56 lorsque toutes les valeurs X,(n-1) ont été traitées. Ensuite, une étape 58 d'application de matrice de diffusion binaire est mise en oeuvre. Une matrice de diffusion binaire de taille mxm est mémorisée, et l'application de la matrice consiste en une multiplication suivie de l'application d'une opération « ou exclusif » bit à bit, pour chaque i de i=1 à m, comme expliqué ci-dessus en référence à la figure 2. Un ensemble de valeurs chaotiques combinées {Yi(n),...,Y,(n)} est obtenu à l'issue de cette étape. Ensuite, une perturbation telle que décrite ci-dessus est appliquée à l'étape 60 sur chaque Y,(n) pour obtenir une valeur XJ(n), pour j=1,...,m, et les valeurs X(n) obtenues sont mémorisées. Ensuite il est vérifié à l'étape de test 62 si n est inférieur à Ns, qui est le nombre total d'échantillons à générer par séquence, à l'étape 62. Si n est inférieur à Ns, l'étape 62 est suivie de l'étape 64 d'incrémentation de n, suivie de l'étape 54 déjà décrite. Si n est supérieur à Ns, le procédé prend fin, les séquences chaotiques {X, (m)} pour i=1,..., m et n=1,...Ns étant fournies en sortie de ce procédé. Il est à noter qu'en variante, la technique de perturbation n'est appliquée qu'a un sous-ensemble des valeurs chaotiques combinées Y,(n). Les inventeurs ont mis en oeuvre divers tests définis par le NIST, dits tests de NIST, permettant de mesurer les propriétés de confusion-diffusion des séquences chaotiques générées, et ont pu montrer de très bonnes propriétés. En particulier, les probabilités des valeurs 0 et 1 mesurées sur 10000 échantillons générés sont très proches de 0,5, donc la distribution est quasi uniforme et l'inter-corrélation des séquences générées est très faible. Le tableau 1 ci-dessous illustre la corrélation p(X ' X j) calculée pour Ns=10000 échantillons, pour m=4 : Indices(i,j) 1 2 3 4 1 1 0,00271649 -0,0020199 -0,01118557 2 0,00271649 1 -0,01352033 -0,00860343 3 -0,0020199 _0,01352033 1 0,0083838 4 -0,01118557 _0,00860343 0,0083838 1 Tableau 1 Le tableau 2 illustre les probabilités d'avoir des bits égaux à « 0 » et à « 1 », notées respectivement Pr _0 et Pr 1, dans une séquence chaotique de Ns=10000 échantillons, pour m=4 : Pr _0 Pr _1 X1 0,50049688 0,49950313 X2 0,50073125 0,49926875 X3 0,49981875 0,50018125 X4 0,50096250 0,4990375 Tableau 2 Pour montrer la sensibilité à la clé secrète, le test suivant a été effectué. Pour chaque entrée d'indice i, la conditions initiale X(0) a été modifiée par changement du bit de poids faible : X:(0) = X(0) 01, et les sorties respectives X(n) ont été générées, pour tous les i allant de 1 à m. Ensuite, la corrélation p(X i, X: ) a été calculée, ainsi que la distance de Hamming PDH entre les séquences X i (n) et X: (n) : PDH (X t, Nsx2N X i (1) X(l) 1=1 X % Ns x 2N Où Ns est le nombre d'échantillons testés. Ns=10000 pour les résultats donnés dans le tableau 3 ci-dessous, et m=4. La distance de Hamming est exprimée en pourcentage. PDH P X1 50,0846875 0,00085691 X2 50,0190625 0,01004383 X3 49,971875 -0,00973299 X4 49,9796875 -0,0053712 Tableau 3 Comme on peut le constater à partir des tableaux 1 à 3, de très bons résultats sont obtenus avec m=4. Des résultats encore améliorés ont été obtenus pour les valeurs supérieures de m, à savoir m=8, 16 et 32. Ainsi, le générateur de séquences chaotiques selon l'invention permet d'obtenir des séquences chaotiques avec de très bonnes performances cryptographiques de confusion-diffusion, protégées par des clés secrètes de plus de 300 bits, qui trouvent des applications dans tous les domaines nécessitant des séquences de nombres pseudoaléatoires : chiffrement par bloc et par flux, stéganographie, tatouage numérique robuste et génération de clés secrètes. REVENDICATIONS1.- Générateur (2) de séquences chaotiques de nombres de valeurs entières représentées sur un nombre de bits prédéterminé, lesdites séquences étant destinées à être utilisées notamment dans des systèmes de cryptage d'informations à base de clé, ledit générateur comportant un nombre m de moyens (4, 6, 8, 10) de mise en oeuvre de fonctions non linéaires, caractérisé en ce qu'il comporte, en sortie desdits moyens (4, 6, 8, 10) de mise en oeuvre de fonctions non linéaires, des moyens de combinaison (20) aptes à combiner les sorties desdits moyens (4, 6, 8, 10) de mise en oeuvre de fonctions 10 non linéaires par application d'une matrice de diffusion binaire donnée, permettant d'obtenir un nombre m de sorties, chaque valeur de sortie étant calculée par une combinaison binaire de sorties desdits moyens (4, 6, 8, 10) de mise en oeuvre de fonctions non linéaires. 15 2.- Générateur de séquences chaotiques selon la revendication 1, caractérisé en ce qu'il comporte en outre des moyens (22, 24, 26, 28) de perturbation raccordés en sortie des moyens de combinaison (20). 3.- Générateur de séquences chaotiques selon la revendication 2, caractérisé en 20 ce qu'il comporte m moyens (22, 24, 26, 28) de perturbation, un dit moyen de perturbation étant raccordé à chaque sortie des moyens de combinaison (20). 4.- Générateur de séquences chaotiques selon l'une des revendications 2 ou 3, caractérisé en ce qu'un dit moyen (22, 24, 26, 28) de perturbation comprend un registre à 25 décalage à réaction. 5.- Générateur de séquences chaotiques selon l'une des revendications 1 à 4, caractérisé en ce que lesdits moyens (4, 6, 8, 10) de mise en oeuvre de fonctions non linéaires comportent des moyens de mise en oeuvre de cartes chaotiques d'un premier 30 type et des moyens de mise en oeuvre de cartes chaotiques d'un deuxième type. 6.- Générateur de séquences chaotiques selon la revendication 5, caractérisé en ce que les moyens de mise en oeuvre de cartes chaotiques de premier type et les moyens de mise en oeuvre de cartes chaotiques de deuxième type sont alternés. 35 7.- Générateur de séquences chaotiques selon l'une des revendications 5 ou 6, caractérisé en ce que les cartes chaotiques de premier type sort des cartes chaotiques de type « skew tent » et les cartes chaotiques de deuxième type sont des cartes chaotiques cartes chaotiques linéaires par morceaux PWLCM. 8.- Générateur de séquences chaotiques selon l'une des revendications précédentes, caractérisé en ce que les nombre d'entrées et de sorties des moyens de combinaison est égal à l'un des nombres suivants : 4, 8, 32, 64. 10 9.- Générateur de séquences chaotiques selon l'une des revendications précédentes, caractérisé en ce que les moyens de combinaison (20) comportent, pour l'obtention d'une sortie, un nombre m d'interrupteurs (40, 42. 44, 46) branchés en sortie des moyens (4, 6, 8, 10) de mise en oeuvre de fonctions non linéaires, un interrupteur fermé correspondant à la présence d'un élément égal à un dans la matrice de diffusion 15 binaire et un interrupteur ouvert correspondant à la présence d'un élément égal à zéro dans la matrice de diffusion binaire, les sorties desdits interrupteurs étant raccordées à une porte « ou exclusif » (48). 10.- Procédé de génération de séquences chaotiques de nombres de valeurs 20 entières représentées sur un nombre de bits prédéterminé, mis en oeuvre par un générateur de séquences chaotiques selon la revendication 1, lesdites séquences étant destinées à être utilisées notamment dans des systèmes de cryptage d'informations à base de clé, le procédé étant caractérisé en ce qu'il comporte les étapes de : - application (54), de m fonctions non linéaires sur m valeurs initiales données, 25 permettant d'obtenir m valeurs chaotiques, - combinaison (62) des m valeurs chaotiques par application d'une matrice de diffusion binaire, permettant d'obtenir m valeurs chaotiques combinées. 11. Procédé de génération de séquences chaotiques selon la revendication 10, 30 caractérisé en ce qu'il comporte en outre une étape d'application (66) d'une perturbation sur au moins une sous-partie des m valeurs chaotiques combinées.