%%
clear all
[DN,DT]=xlsread( 'Nivel_Sistema_Cantareira.xls', '2003_2014' );
ordem = 15*(0:(size(DN,1)-1)).';
nivel = DN(:,5);
%%
figure
subplot( 3, 1, 1:2 )
plot( ordem, nivel, '*' )
set( gca, 'fontsize', 16 )
ylabel( 'Nível do sistema Cantareira (%)' )
% xlabel( 'Dia desde 1/1/2003 (dias)' )

Nm = 50;
AA = 10;
fA = -pi/2;
AS = 40;
fS = pi/2;
PAR = [ Nm AA fA AS fS ];
Nfit = funcao_nivel( PAR, ordem );
hold on
% plot( ordem, Nfit, '-r' )

PAR_aj = nlinfit( ordem, nivel, 'funcao_nivel', PAR )
N_aj = funcao_nivel( PAR_aj, ordem );
plot( ordem, N_aj, '-k' )
xlim( [0 max(ordem)] )

subplot( 3, 1, 3 )
plot( ordem, nivel - N_aj, '*r' )
xlim( [0 max(ordem)] )
set( gca, 'fontsize', 16 )
ylabel( 'Diferença (%)' )
xlabel( 'Dia desde 1/1/2003 (dias)' )
