demo_biosensor
Demo of denoising using SASS
Biosensor data from Stephen Arnold
Ivan Selesnick December 2011. Revised August 2016
Contents
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)')