(Oszczędne próbkowanie)
Compressive sensing is a mathematical tool that creates hi-res data sets from lo-res samples.
It can be used to:
- resurrect old musical recordings,
- find enemy radio signals,
- and generate MRIs much more quickly.
Undersample
Fill in the dots
(l1)
Add shapes
(sparsity)
Add smaller shapes
Achieve clarity
What is the key?
Candès and Tao have shown mathematically that
the l1 minimization is all we need
library(R1magic)
N <- 100
# Sparse components
K <- 4
# Up to Measurements > K LOG (N/K)
M <- 40
# Measurement Matrix (Random Sampling Sampling)
phi <- GaussianMatrix(N,M)
# R1magic generate random signal
xorg <- sparseSignal(N, K, nlev=1e-3)
y <- phi %*% xorg ;# generate measurement
T <- diag(N) ;# Do identity transform
p <- matrix(0, N, 1) ;# initial guess
# R1magic Convex Minimization
ll <- solveL1(phi, y, T, p)
x1 <- ll$estimate
plot( 1:100, seq(0.011,1.1,0.011), type = "n",xlab="",ylab="")
title(main="Random Sparse Signal Recovery",
xlab="Signal Component",ylab="Spike Value")
lines(1:100, xorg , col = "red")
lines(1:100, x1, col = "blue", cex = 1.5)