Confusion au sujet du problème de détection compressé

13

J'ai lu quelques références dont celle-ci .

Je suis un peu confus quel problème d'optimisation compressé détecte les builds et essaie de résoudre. Est-ce

minimizex1subject toAx=b

ou et

minimizex0subject toAx=b

ou / et autre chose?

Tim
la source

Réponses:

18

Brian est sur place. Mais je pense qu'il est utile d'ajouter un contexte de détection compressé.

Tout d'abord, notez que la soi-disant norme 0 - la fonction de cardinalité, ou le nombre de valeurs non nulles dans x - n'est pas une norme . Il est probablement préférable de l'écrire comme quelque chose comme la carte ( x ) dans tout sauf les contextes les plus décontractés. Ne vous méprenez pas, vous êtes en bonne compagnie lorsque vous utilisez le x 0x0xcard(x)x0 raccourci , mais je pense que cela a tendance à semer la confusion.

Les gens savent depuis longtemps que minimiser la norme x 1 tend à produire des solutions éparses. Il existe des raisons théoriques à cela liées à la complémentarité linéaire. Mais ce qui était le plus intéressant, ce n'était pas que les solutions étaient rares, mais qu'elles étaient souvent les plus rares possibles . C'est-à-dire que minimiser x 1 vous donne vraiment la solution de cardinalité minimale dans certains cas utiles. (Comment ont-ils compris cela, alors que le problème de cardinalité minimale est NP-difficile? En construisant des problèmes artificiels avec des solutions rares connues.) Ce n'était pas quelque chose que la théorie de la complémentarité linéaire pouvait prédire.1x1x1

Le domaine de la détection compressée est né lorsque les chercheurs ont commencé à identifier des conditions sur la matrice qui leur permettraient de garantir à l'avance que la solution 1 était également la plus rare. Voir, par exemple, les premiers articles de Candés, Romberg et Tao , et d'autres discussions sur la propriété d'isométrie restreinte, ou RIP . Un autre site Web utile si vous voulez vraiment plonger dans une théorie est la page de détection compressée de Terence Tao .A1

Michael Grant
la source
12

Nous serions ravis de pouvoir résoudre

minx0

st

Ax=b

mais ce problème est un problème d'optimisation combinatoire NP-difficile qui n'est pas pratique à résoudre dans la pratique lorsque , x et b sont de tailles typiques dans la détection en compression. Il est possible de résoudre efficacementAxb

minx1

st

Ax=b

à la fois en théorie (cela peut être fait en temps polynomial) et en pratique informatique pour des problèmes même assez importants qui surviennent dans la détection compressive. Nous utilisons comme "substitut" pour le x 0 . Cela a une justification intuitive (la minimisation à une norme préfère les solutions avec moins d'entrées non nulles dans x ), ainsi que des justifications théoriques beaucoup plus sophistiquées (théorèmes de la forme "Si A x = b a une solution k-clairsemée, minimisant alors x 1x1x0xAx=bx1 trouvera cette solution avec une forte probabilité « .

Ax=bAxb2δ

minAxb22+λx1

Brian Borchers
la source
8

Je n'ai rien à ajouter aux explications de Brians et Michaels 1 contre. 0. Mais comme la question semble concerner la détection compressée, j'aimerais ajouter mon point de vue: la détection compressée ne consiste pas à résoudre

minX0stUNEX=b
ni sur
minX1stUNEX=b.
La détection comprimée est plus un paradigme , qui peut être énoncé très grossièrement comme

Il est possible d'identifier des signaux clairsemés à partir de quelques mesures.

La détection compressée consiste en réalité à prendre le moins de mesures possible pour identifier un signal dans une certaine classe de signaux.

Une phrase accrocheuse est:

Pourquoi votre appareil photo 5 mégapixels devrait-il vraiment mesurer 15 millions de valeurs (trois pour chaque pixel) qui vous coûtent 15 mégaoctets de données alors qu'il ne stocke qu'environ 2 mégaoctets (après compression)?
Serait-il possible de mesurer les 2 mégaoctets tout de suite?

Il existe des cadres assez différents:

  • mesures linéaires
  • non linéaires (par exemple "sans phase")
  • données vectorielles, données matricielles / tensorielles
  • la rareté comme juste le nombre de non-zéros
  • clarté comme «bas rang» ou même «faible complexité»).

Et il existe également plus de méthodes pour calculer des solutions clairsemées comme les poursuites d'appariement (des variantes comme la poursuite d'appariement orthogonal (OMP), la poursuite d'appariement orthogonal régularisé (ROMP), CoSaMP) ou les méthodes plus récentes basées sur des algorithmes de transmission de messages .

Si l'on identifie la détection comprimée avec de simples 1- ou 0-minimisation, on manque de beaucoup de flexibilité face à des problèmes pratiques d'acquisition de données.

Cependant, si l'on ne s'intéresse qu'à l'obtention de solutions éparses aux systèmes linéaires, on fait quelque chose que j'appellerais une reconstruction éparse .

Poignard
la source
Merci! Pouvez-vous reformuler ce qui suit dans la formulation mathématique: "Il est possible d'identifier des signaux clairsemés à partir de quelques mesures. La détection compressée consiste en réalité à prendre le moins de mesures possible pour identifier un signal dans une certaine classe de signaux."
Tim
1
Non, je ne peux pas, car la Compressed Sensing n'est pas une théorie mathématique mais plutôt un concept d'ingénierie.
Dirk
1
Je pense que cette réponse est une bonne contribution et je l'ai votée. Quant à la phrase accrocheuse, cependant, j'ai toujours eu un problème avec elle. Cela suggère que CS est si puissant que vous pouvez simplement jeter 13 millions de pixels et récupérer l'image de toute façon. Mais non, vous ne devez jamais jeter les données au hasard, même dans un système CS --- un bon algorithme de récupération peut toujours utiliser plus de données. La promesse de CS est le potentiel de développer des capteurs qui collectent moins de données en premier lieu en échange d'importantes économies pratiques: économies d'énergie, collecte plus rapide, etc.
Michael Grant
@MichaelGrant Je suis d'accord: ne jetez pas les données que vous avez déjà mesurées et utilisez la date que vous pouvez mesurer avec un minimum d'effort.
Dirk