demo_biosensor

Demo of denoising using SASS

Biosensor data from Stephen Arnold

Ivan Selesnick
December 2011. Revised August 2016

Contents

Load data
Define low-pass filter
SASS - run algorithm
SASS - plot denoised signal
SASS - sparse signal component

Load data

clear

y = load('data/wgm_data.txt');           % load WGM sensor data

N = length(y);          % N : 1500
n = 1:N;

figure(1)
clf
plot(n, y)
title('Biosensor data');
print -dpdf figures/demo_biosensor_data

Define low-pass filter

d = 1;
fc = 0.005;

SASS - run algorithm

K = 1;
lam = 60.0;

[x, u, cost, v] = sass_L1(y, d, fc, K, lam);

txt = sprintf('SASS (K = %d, lam = %.2f, d = %d, fc = %.3f)', K, lam, d, fc);

0 incorrectly locked zeros detected.

SASS - plot denoised signal

figure(1)
clf
plot(n, x);
title(txt)
print -dpdf figures/demo_biosensor_SASS

SASS - sparse signal component

figure(1)
clf
subplot(2, 1, 1)
plot(u)
xlim([0 N])
title('Sparse signal (determined via SASS)')