try_quantize.m

%-----------------------------------------------------------------------------
%
% Copyright (C) 2016
% ASTRON (Netherlands Institute for Radio Astronomy) <http://www.astron.nl/>
% P.O.Box 2, 7990 AA Dwingeloo, The Netherlands
%
% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program.  If not, see <http://www.gnu.org/licenses/>.
%
%-----------------------------------------------------------------------------
% Author: E. Kooistra, 2016
%
% Purpose : Try quantize.m
% Description :

clear all;
close all;
fig=0;

% User settings
nof_bits = 3;

x_period = 2^nof_bits;
x_max = 2^(nof_bits-1)-1;
x_step = 0.5;
%x_step = 1;
x = -2*x_period:x_step:2*x_period-1;

y_half_away_wrap = quantize(x, x_max, nof_bits, 'half_away', 'wrap');
y_half_away_clip = quantize(x, x_max, nof_bits, 'half_away', 'clip');
y_half_up_wrap = quantize(x, x_max, nof_bits, 'half_up', 'wrap');
y_half_up_clip = quantize(x, x_max, nof_bits, 'half_up', 'clip');

% Plot
fig=fig+1;
figure('position', [300+fig*20 200-fig*20 1000 800]);
figure(fig);
plot(x, y_half_away_wrap, 'kx', x, y_half_away_clip, 'ro');
ylim([-x_period/2-0.5 x_period/2+0.5])
title(['Quantisation with half away and ', num2str(nof_bits), ' bits'] );
grid on;

fig=fig+1;
figure('position', [300+fig*20 200-fig*20 1000 800]);
figure(fig);
plot(x, y_half_up_wrap, 'kx', x, y_half_up_clip, 'ro');
ylim([-x_period/2-0.5 x_period/2+0.5])
title(['Quantisation with half up and ', num2str(nof_bits), ' bits'] );
grid on;