# Declip Saturated Signals Using Your Own Function

Sensors can return clipped readings if the data are larger than a given saturation point. To reconstruct the readings, you can fit a polynomial through the points adjacent to the saturated intervals. Write a function that performs the reconstruction and integrate it into Signal Analyzer.

Generate a two-channel signal sampled at 1 kHz for 14 seconds. The signal has several peaks of varying sizes and shapes. A sensor that reads the signal saturates at 0.1 V.

fs = 1000;
t = 0:1/fs:14-1/fs;

0.85*besselj(0,5*(sin(2*pi*(t+1.5).^2/20).^2)).*sin(2*pi*t/9)]';

sigsat = sig;
stv = 0.1;
sigsat(sigsat >= stv) = stv;

Open Signal Analyzer and drag the original signal and the saturated signal to the Signal table. Drag the first channel of each signal to the top display, and the second channel of each signal to the bottom display.

Write a function that uses a polynomial to reconstruct the signal peaks:

• The first input argument, x, is the input signal. This argument must be a vector and is treated as a single channel.

• The second input argument, tIn, is a vector of time values. The vector must have the same length as the signal. If the input signal has no time information, the function reads this argument as an empty array.

• Use varargin to specify additional input arguments. If you do not have additional input arguments, you can omit varargin. Enter the additional arguments as an ordered comma-separated list in the Function Parameters panel inside the preprocessing mode.

• The first output argument, y, is the preprocessed signal.

• The second output argument, tOut, is a vector of output time values. If the input signal has no time information, tOut is returned as an empty array.

• To implement your algorithm, you can use any MATLAB® or Signal Processing Toolbox™ function.

function [y,tOut] = declip(x,tIn,varargin)
% Declip saturated signal by fitting a polynomial

% Initialize the output signal

y = x;

% For signals with no time information, use sample numbers as abscissas

if isempty(tIn)
tOut = [];
t = (1:length(x))';
else
t = tIn;
tOut = t;
end

% Specify the degree of the polynomial as an optional input argument
% and provide a default value of 4

if nargin<3
ndx = 4;
else
ndx = varargin{1};
end

% To implement your algorithm, you can use any MATLAB or Signal
% Processing Toolbox function

% Find the intervals where the signal is saturated and generate an
% array containing the interval endpoints

idx = find(x==max(x));
fir = [true;diff(idx)~=1];
ide = [idx(fir) idx(fir([2:end 1]))];

% For each interval, fit a polynomial of degree ndx over the ndx+1 points
% before the interval and the ndx+1 points after the interval

for k = 1:size(ide,1)
bef = ide(k,1); aft = ide(k,2);
intv = [bef-1+(-ndx:0) aft+1+(0:ndx)];
[pp,~,mu] = polyfit(t(intv),x(intv),ndx);
y(bef:aft) = polyval(pp,t(bef:aft),[],mu);
end

end

Add the function to Signal Analyzer as a custom preprocessing function. Select sigsat in the Signal table and on the Analyzer tab, click Preprocess to enter the preprocessing mode. In the Functions gallery, select Add Custom Function. Input the function name and description. Paste the text of your function in the editor window that appears. Save the file. The function appears in the Custom Functions list in the Functions gallery.

Demonstrate that the function you created reconstructs the saturated regions.

1. Select the first channel of the saturated signal in the Signal table.

2. In the Functions gallery, select declip.

3. In the Function Parameters panel, click Apply.

4. Click Accept All.

1. Select sig and sigsat in the Signal table. Do not select individual channels.

2. On the Analyzer tab, click Time Values. Select Sample Rate and Start Time and specify fs as the sample rate.

Check that the function works when you specify optional inputs.

1. Select the second channel of the saturated signal in the Signal table.

2. Click Preprocess, and select declip from the Functions gallery. In the Function Parameters panel, enter 8 in the Arguments field and click Apply. The preprocessing function uses a polynomial of degree 8 to reconstruct the saturated regions.

3. Click Accept All.