Comparer deux capteurs : correlation des séries réduites

Partie 1 · ★★★★

Deux capteurs mesurent le même phénomène à des échelles différentes. Écris comparer(A, B) qui centre-réduit les deux séries (écart-type population) puis renvoie le tuple (r, [c1, c2, ...]) où r est la moyenne des produits $A_i^*B_i^*$ et c_i = (A*_i + B*_i)/racine(2) est la projection sur $u=(1/\sqrt2,1/\sqrt2)$ . Sur A=[1,2,3], B=[10,20,30] : $r=1$ .

def comparer(A, B):
    def cr(xs):
        n = len(xs)
        m = sum(xs) / n
        s = (sum((x - m) ** 2 for x in xs) / n) ** 0.5
        return [(x - m) / s for x in xs]
    a, b = cr(A), cr(B)
    r = sum(x * y for x, y in zip(a, b)) / len(a)
    c = [(x + y) / ___ for x, y in zip(a, b)]
    return (r, c)

Ctrl/⌘+Entrée