function eval_gainControl()

close all;

% 1
G_Target = [40*ones(1,200) 20*ones(1,200) 40*ones(1,200)];
G_PA = 43+5*sin(1/400*2*pi*(1:600));

gainControl(G_Target, G_PA, 3, 0.2);

% 2
figure
G_Target = [40*ones(1,600)];
G_PA = [43*ones(1,200) 46*ones(1,200) 23*ones(1,200)];

gainControl(G_Target, G_PA, 3, 0.2);

% 2
figure
G_Target = [40*ones(1,600)];
G_PA = [43*ones(1,200) 40*ones(1,200) 83*ones(1,200)];

gainControl(G_Target, G_PA, 3, 0.2);