clear;

sin_params = [0.3 0 0.001; 0.5 0 0.002];
sig = genSumSin(sin_params(:,1),sin_params(:,2),sin_params(:,3),1000);

[sig_sd,e] = sigmadelta(sig);


figure;
hold on
    plot(sig)
    plot(sig_sd)
hold off

% figure
% plot(e)

% sig_rec - signal reconstruction
for i = 1:1:length(sig_sd)
    sig_rec(i) = sum(sig_sd(1:i));
end

% skalowanie do podobnej amplitudy
sig_rec = ( sig_rec - mean(sig_rec) )/length(sig_sd)*10;

figure;
hold on
    plot(sig_rec)
    plot(sig)
hold off