: WiFi, WiMAX, LTE, UMTS, . . .; ⇒ H´ et´ erog´ en´ eit´ e de l’environnement radio. Tirer profit de cette h´ et´ erog´ en´ eit´ e afin de satisfaire le concept de Always Best Connected (ABC); Qu’est ce que le ABC ? Une personne n’est pas seulement toujours connect´ ee, mais aussi connect´ ee en tout temps ` a un r´ eseau disponible offrant la qualit´ e de service la plus adapt´ ee ` a ses besoins 2/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
terminaux multistandards capables de migrer intelligemment d’un r´ eseau ` a un autre ` a la recherche de la QoS la plus ad´ equate; Ce processus de migration est appel´ e “Vertical Handover”. 3/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
` a un Horizontal Handover mais o` u la communication est prise en charge par un autre syst` eme; ´ Etapes d’un Vertical Handover : 1. Identification des syst` emes actifs dans l’environnent radio; 2. Analyse de la QoS disponible sur chaque r´ eseau d´ etect´ e; 3. D´ ecision et choix du r´ eseau le plus ad´ equat; 4. Ex´ ecution du Vertical Handover. 4/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
Triviale dans un contexte d’allocation statique du spectre; Probl´ ematique dans un contexte d’allocation dynamique du spectre (Radio Opportuniste); Les standards n’´ emettent plus sur des bandes connues. ⇒ N´ ecessit´ e d’identifier les r´ eseaux secondaires actifs. 5/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
eseau identifi´ e Contrairement au HHO, le RSB n’est pas une m´ etrique suffisante; Des informations et param` etres additionnels (“m´ etriques de handover“) doivent ˆ etre estim´ es : d´ ebit effectif, ´ etat de la batterie, . . . VHO transparent : estimation passive des m´ etriques sans se connecter, ´ economie en temps et en ´ energie. 6/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
plus ad´ equat Une fois les m´ etriques estim´ ees, une d´ ecision est prise en fonction de la QoS requise par les couches hautes et les informations collect´ ees. 4. Ex´ ecution du Vertical Handover Une fois le r´ eseau cible choisi, la connexion du terminal est transf´ er´ ee vers celui-ci. 7/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
travaux ? Concevoir des algorithmes de radio cognitive d´ edi´ es ` a : L’identification de syst` emes OFDM; L’estimation passive de m´ etriques de handover. 8/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
travaux ? Concevoir des algorithmes de radio cognitive d´ edi´ es ` a : L’identification de syst` emes OFDM; L’estimation passive de m´ etriques de handover. Pourquoi OFDM ? Sch´ ema de modulation tr` es r´ epandu de nos jours : WiFi, WiMAX, LTE, . . . 8/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
travaux ? Concevoir des algorithmes de radio cognitive d´ edi´ es ` a : L’identification de syst` emes OFDM; L’estimation passive de m´ etriques de handover. Pourquoi OFDM ? Sch´ ema de modulation tr` es r´ epandu de nos jours : WiFi, WiMAX, LTE, . . . Quels protocoles de multi-acc` es ? OFDMA (LTE, WiMAX Mobile); CSMA/CA (WiFi) 8/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
◮ Quelle signature ? Introduction N´ ecessit´ e : La fr´ equence porteuse n’est plus discriminante; Il faut trouver une autre signature; Exemple : Espacement entre sous-porteuses, cyclo-stationnarit´ e, motif de pilotes, . . . Approche retenue : motif de pilotes. 11/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
◮ Quelle signature ? Introduction N´ ecessit´ e : La fr´ equence porteuse n’est plus discriminante; Il faut trouver une autre signature; Exemple : Espacement entre sous-porteuses, cyclo-stationnarit´ e, motif de pilotes, . . . Approche retenue : motif de pilotes. freq time LTE Data tones Pilot tones IEEE 802.22 WiMAX DL-PUSC 11/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
◮ Quelle signature ? Pilot Induced Cyclostationnarity (PIC) Les auteurs montrent que si une corr´ elation est introduite entre paire de pilotes; Alors le signal rec ¸u pr´ esente une cylostationnarit´ e. freq time Data tones Pilot tones F.-X. Socheleau et al, “Cognitive OFDM System Detection Using Pilot Tones Second and Third-Order Cyclostationarity,” Elsevier Signal Processing, 2011 12/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
◮ Mod` ele du signal Mod` ele du signal Supposant qu’un signal OFDM est observ´ e sur une bande donn´ ee, apr` es FFT et extraction du pr´ efixe cyclique : Yn,k ∆ = Hkcn,k +Wn,k, o` u cn,k = pn,k if (n,k) ∈ P, dn,k, if (n,k) ∈ D. ⇒ Construire un test d’hypoth` ese Oularbi et al ≪ OFDM system identification using pilot tone signature ≫, WOSSPA 2011. Oularbi et al, ≪ Identification de syst` emes OFDM ` a partir de la signature des pilotes ≫, Gretsi 2011. Oularbi et al, ≪ GLR Test for OFDM System Identification Usning Pilot Tones Pattern ≫, article de revue soumis. 13/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
◮ Approche GLRT Les hypoth` eses H0 : Le motif de pilote P est absent; H1 : Syst` eme utilisant le motif P est actif. H0 : H ′ 0 : Seul du bruit est observ´ e; H ′′ 0 : Le syst` eme actif n’utilise pas le motif P. ⇒ N´ ecessit´ e d’effectuer une pr´ e-d´ etection 14/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
◮ Structure de pr´ e-d´ etection Structure de pr´ e-d´ etection H ′ 0 : Seul du bruit est observ´ e; H ′′ 0 : Le syst` eme actif n’utilise pas le motif P. f(Y;H ′ 0 ,σ2) = 1 |πσ2|MsNsc exp − Y 2 σ2 , f(Y;H ′′ 0 ,σ2 k ) = ∏Nsc −1 k=0 1 (πσ2 k )Ms exp − 1 σ2 k Y(k) 2 . Les inconnues σ2 et σ2 k sont estim´ ees au sens du M.V σ2 = 1 Ms Nsc tr(YYH) et σ2 k = 1 Ms Y(k) 2 15/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
◮ Structure de pr´ e-d´ etection Structure de pr´ e-d´ etection H ′ 0 : Seul du bruit est observ´ e; H ′′ 0 : Le syst` eme actif n’utilise pas le motif P. f(Y;H ′ 0 ,σ2) = 1 |πσ2|MsNsc exp − Y 2 σ2 , f(Y;H ′′ 0 ,σ2 k ) = ∏Nsc −1 k=0 1 (πσ2 k )Ms exp − 1 σ2 k Y(k) 2 . Les inconnues σ2 et σ2 k sont estim´ ees au sens du M.V σ2 = 1 Ms Nsc tr(YYH) et σ2 k = 1 Ms Y(k) 2 L(Y) = Ms Nsc log Y 2 −Ms ∑ k log ||Y(k)||2 H ′′ 0 ≷ H ′ 0 ηpre 15/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
◮ Structure d’identification Calcul du LLR de pr´ e-detection L(Y) L(Y) ? > ηpre Pfa,pre Observer la prochaine sous-bande Calcul du test statistique d’identification T ′(Y,P) P T ′(Y,P) ? > η Pfa Observer Y Un syst` eme utilisant le motif de pilotes P est pr´ esent oui non non oui FIGURE : Diagramme de l’identificateur propos´ e. 19/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
◮ R´ esultats de simulations Comparaison avec la technique PIC −13 −12 −11 −10 −9 −8 −7 −6 10−3 10−2 10−1 100 RSB (dB) probabilit´ e de non d´ etection GLRT RSB connu D´ etecteur PIC GLRT RSB estim´ e FIGURE : Probabilit´ e de non-d´ etection en fonction du RSB, Signaux WiMAX Mobile DL-PUSC, Nsc = 512, Ms = 24, Pfa = 0.02 21/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
des r´ eseaux OFDMA 1 Identification de syst` emes OFDM 2 M´ etriques de VHO pour des r´ eseaux OFDMA 3 M´ etriques de VHO pour des r´ eseaux CSMA/CA 4 Conclusion et perspectives 22/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
etriques de VHO pour des r´ eseaux OFDMA ◮ Introduction Introduction Variante multi-utilisateurs du sch´ ema OFDM; WiMAX Mobile, LTE (DownLink), IEEE 802.22; On s’int´ eresse ` a deux m´ etriques : • Nombre d’antennes de la station de base; • Taux d’allocation des slots temps-fr´ equence; Estimation passive effectu´ ee ` a partir d’observations sur la couche physique. 23/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux OFDMA ◮ Estimation du nombre d’antennes d’une B.S 1. Estimation du nombre d’antennes d’une B.S But : Mieux interpr´ eter le RSB; Utilisation de l’orthogonalit´ e entre motifs de pilotes : k n Antenne 1 Antenne 1 Antenne 1 Antenne 1 Antenne 1 Antenne 1 Antenne 1 Antenne 1 Antenne 1 Antenne 1 Antenne 1 Antenne 1 Antenne 2 Antenne 2 Antenne 2 Antenne 2 Antenne 2 Antenne 2 Antenne 2 Antenne 2 Antenne 2 Antenne 2 Antenne 2 Antenne 2 Tons pilotes Tons nuls Tons de donn´ ees Utiliser le GLRT ou le PIC pour v´ erifier la pr´ esence de chaque motif et d´ eduire le nombre d’antennes; Avantage : Une seule antenne est suffisante. 24/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux OFDMA ◮ Estimation du nombre d’antennes d’une B.S R´ esum´ e de la technique propos´ ee : Initialiser le Nom- bre d’antennes i = 1 P(i) Pr´ esent ? i = i + 1 i−1 antennes pr´ esentes oui non FIGURE : Diagramme de l’algorithme propos´ e. 25/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux OFDMA ◮ Estimation du nombre d’antennes d’une B.S Extension GLRT −10 −8 −6 −4 −2 0 2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 RSB (dB) Probabilit´ e de d´ etection 1 antenne 2 antennes 4 antennes FIGURE : Probabilit´ e de d´ etection du nombre d’antennes en fonction du RSB, en utilisant l’extension du GLRT, signaux LTE, Ms = 48 symboles observ´ es, Pfa = 0.02. Oularbi et al, ≪ GLR Test for BS number of antennas estimation ≫, IEEE com letters, accepted 2012.. 26/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux OFDMA ◮ Estimation du nombre d’antennes d’une B.S Extension PIC −10 −8 −6 −4 −2 0 2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 RSB (dB) Probabilit´ e de d´ etection 1 antenne 2 antennes 4 antennes FIGURE : Probabilit´ e de d´ etection du nombre d’antennes en fonction du RSB, en utilisant l’extension du PIC, signaux LTE, Ms = 48 symboles observ´ es, Pfa = 0.02. Oularbi et al, ≪ PIC detector for BS number of antennas estimation ≫, PIMRC 2012, Sydney. 27/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux OFDMA ◮ Estimation du taux d’allocation temps-fr´ equence 2. Estimation du taux d’allocation temps-fr´ equence Plus le nombre de sous porteuses actives est faible, plus le nombre d’utilisateurs et le trafic sur cette station de base le sont aussi. sous−porteuses temps Slots inactifs εk,n = 0 Slots actifs εk,n = 1 Il est pr´ ef´ erable pour un terminal multi-standards de se connecter ` a une BS o` u le taux d’activit´ e temps-fr´ equence est faible. ⇒ N´ ecessit´ e d’estimer ce taux. 28/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux OFDMA ◮ Estimation du taux d’allocation temps-fr´ equence 2. Estimation du taux d’allocation temps-fr´ equence Ce taux est diffus´ e sur certains syst` emes : N´ ecessit´ e d’intercepter les pr´ eambules des trames; Autre approche : estimation passive (Nombre de slots actifs divis´ e par le nombre total des slots observ´ es); Revient ` a un probl` eme de d´ etection; Oularbi et al.,≪ Blind estimation of the time-frequency activity rate of OFDMA signals ≫, ICUMT 2010. Oularbi et al.,≪ Estimation du taux d’activit´ e temps-fr´ equence d’un signal OFDMA pour des applications en radio- cognitive ≫, Gretsi 2011. 29/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux OFDMA ◮ Estimation du taux d’allocation temps-fr´ equence Mod` ele du signal Yn,k ∆ = εn,kHn,kcn,k +Wn,k. o` u εn,k sont des v.a i.i.d ∈ {0,1}. Si σ2 est connue alors le taux est donn´ e par : ˆ α = ∑n,k I(|Yn,k | > s(σ2)) Ms Nsc . Approche d´ ependante d’un seuil arbitraire; On propose une technique bas´ ee sur les moments d’ordre un et deux du signal observ´ e. 30/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux OFDMA ◮ Estimation du taux d’allocation temps-fr´ equence Algorithme propos´ e Le moment d’ordre deux s’´ ecrit : E[|Yn,k |2] = σ2 +αEs, ∀εn,k. Le moment d’ordre un s’´ ecrit : E[|Yn,k |] = αϕ(Es )+(1−α) √ π 2 σ. o` u ϕ(Es ) est une fonction qui d´ epend de la puissance du signal et de la distribution des symboles ´ emis. 31/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux OFDMA ◮ Estimation du taux d’allocation temps-fr´ equence Algorithme propos´ e L’estim´ ee du taux d’activit´ e temps-fr´ equence est la racine de l’´ equation suivante : α ϕ( µ2 −σ2 ˆ α )+(1− ˆ α) √ π 2 σ− ˆ µ1 = 0. µ1 et µ2 sont les moments d’ordre un et deux empiriques; R´ esolution par dichotomie. 33/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
des r´ eseaux CSMA/CA 1 Identification de syst` emes OFDM 2 M´ etriques de VHO pour des r´ eseaux OFDMA 3 M´ etriques de VHO pour des r´ eseaux CSMA/CA 4 Conclusion et perspectives 36/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
etriques de VHO pour des r´ eseaux CSMA/CA ◮ Introduction Introduction Le protocole CSMA/CA (Carrier Sense Multiple Access with Collision Avoidance) est d´ edi´ e aux r´ eseaux sans-fil; M´ ecanisme d’´ evitement de collisions, car liaison half duplex; R´ eseau le plus connu : WiFi; Nous nous int´ eressons ` a deux m´ etriques de handover : • Taux d’occupation du canal; • Taux de collisions. Oularbi et al, ≪ Physical layer IEEE 802.11 channel occupancy rate estimation ≫, ISIVC 2010. Oularbi et al, ≪ Two vertical handcover metrics toward an IEEE 802.11N network ≫, COGIS 2009. Oularbi et al, ≪ Vertical Handover Metrics Estimation for OFDM Based Systems ≫, Eurasip JWCN. 37/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux CSMA/CA ◮ Estimation du taux d’occupation du canal Qu’est ce que le taux de charge du canal ? Rapport entre le temps durant lequel le canal est d´ eclar´ e occup´ e et la taille de la fenˆ etre d’observation. 38/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux CSMA/CA ◮ Estimation du taux d’occupation du canal Qu’est ce que le taux de charge du canal ? Rapport entre le temps durant lequel le canal est d´ eclar´ e occup´ e et la taille de la fenˆ etre d’observation. Pourquoi cette m´ etrique? R´ ef : Q. Zhang et al, Efficient Mobility Management for vertical Handoff between WWAN and WLAN, IEEE communications magazine. 38/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux CSMA/CA ◮ Estimation du taux d’occupation du canal Particularit´ e d’une communication CSMA/CA Espaces inter-trame qui r´ egularisent les priorit´ es des transmissions; Couche Physique Couche MAC Trame 1 Trame 2 Signal + Bruit Signal + Bruit Bruit seulement Libre ... Bruit ... Espace Inter−trame ⇒ Taux de charge = Nbre d’´ ech signal+bruit Nbre total d’´ ech . 39/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux CSMA/CA ◮ Estimation du taux d’occupation du canal ◮ Technique propos´ ee Technique propos´ ee R´ ecepteur cognitif est dot´ e de N antennes; La variance du bruit est suppos´ ee connue ou estim´ ee; D´ efinition d’un crit` ere de vraisemblance J(u) J(u) = −(Ns −u)log(πσ2)− 1 Nσ2 N ∑ i=1 Ns ∑ m=u |yi (m)|2 . 40/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux CSMA/CA ◮ Estimation du taux d’occupation du canal ◮ R´ esultats exp´ erimentaux R´ esultats exp´ erimentaux FIGURE : Configuration : Un point d’acc` es WiFi, Nombre d’utilisateurs variable, d´ ebit variable. Couche PHY observ´ ee grˆ ace ` a l’USRP (Universal Software Radio Peripheral). ∗ Exp´ erimentations r´ ealis´ ees en collaboration avec S. Hadin (ing´ enieur de recherche, Dpt S&C). 45/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux CSMA/CA ◮ Estimation du taux d’occupation du canal ◮ R´ esultats exp´ erimentaux Influence du nombre d’utilisateurs 0 1 2 3 4 5 6 7 8 9 10 0 5 10 15 20 25 30 Longeur de la fenˆ etre d’observation (ms) Cor en % et variances 1 Client/Serveur 3 Clients/Serveurs FIGURE : Taux de charge mesur´ e en fonction de la taille de la fenˆ etre d’observation, d´ ebit fix´ e ` a 1Mbps. 46/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
pour des r´ eseaux CSMA/CA ◮ Estimation du taux d’occupation du canal ◮ R´ esultats exp´ erimentaux Influence du d´ ebit 0 1 2 3 4 5 6 7 8 9 10 −10 0 10 20 30 40 50 60 70 80 90 Longueur de la fenˆ etre d’observation Cor en % et variances (ms) 100kbps 1Mbps 10Mbps FIGURE : Taux de charge mesur´ e en fonction de la taille de la fenˆ etre d’observation, Nombre de couples C/S fix´ e ` a 3. 47/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
des r´ eseaux CSMA/CA ◮ Estimation du taux de collisions Pourquoi cette m´ etrique ? FIGURE : Probabilt´ e de collision en fonction du d´ ebit offert et du nombre de nœuds. R´ ef : K. Park, Perdu University, “CS 536 : Data Communication and Computer Networks”. 49/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
des r´ eseaux CSMA/CA ◮ Estimation du taux de collisions Comment faire ? Les trames issues d’une collision sont le m´ elange de deux sources ou plus; Estimer le nombre de sources constituant toute trame d´ etect´ ee; Taux de col = Nbre de collisions Total des trames Approche intuitive : Crit` eres d’information th´ eorique (AIC “Akaike information criterion” et MDL “Minimum Description Length”); 50/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
des r´ eseaux CSMA/CA ◮ Estimation du taux de collisions ◮ Mod` ele du signal Solution propos´ ee On consid` ere un r´ ecepteur dot´ e de N antennes; En consid´ erant une fenˆ etre d’observation de taille d : yd (m) = yT(m), ..., yT(m−d +1) T , Le rang r de la matrice d’autocorrelation de yd (m) vaut : r = Md +L On d´ emontre que si d = D (PC), alors M est le premier entier sup´ erieur ` a r 2d . 51/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
de syst` emes OFDM 2 M´ etriques de VHO pour des r´ eseaux OFDMA 3 M´ etriques de VHO pour des r´ eseaux CSMA/CA 4 Conclusion et perspectives 54/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes
propos´ e une nouvelle technique d’identification de syst` emes OFDM bas´ ee sur un GLRT : • Motifs de pilotes, estimation de certains param` etres (canal, variance du bruit, ...) Des techniques d’estimation de m´ etriques de vertical handover pour des syst` emes OFDMA : • Nombre d’antennes de la station de base; • Taux d’allocation temps fr´ equence. Des techniques d’estimation de m´ etriques de vertical handover pour des syst` emes CSMA/CA : • Taux d’occupation du canal; • Taux de collisions. 55/58 Oularbi Rabie (Dpt S&C) S´ eminaire SCEE —20/12/2012— Supelec Rennes