Relatório Sudeste FINAL (+ calculo V)

367 544 658 534 688 321 301 199 142 
118 182 248 726 420 294 399 484 251 180 144 123 88 51 204 533]; 
TT=1:length(vaz); 
an=10; 
svaz0=0; 
u=0; 
for p=1:12; 
n=0; 
vaz0=0; 
u=0; 
for k=1:an; 
j=n*12; 
vaz0=vaz(p+j); 
svaz0(p)=vaz0+u; 
u=svaz0(p); 
n=n+1; 
end; 
end; 
n=0; 
for i=1:12 
svaz0m(i)= (svaz0(i))/an; 
end 
T=1:12; 
grau=3; 
z=polyfit(T,svaz0m,grau); 
86 
 
vcalc=polyval(z,T); 
figure(1); 
plot(TT,vaz); 
xlabel('Meses de 2004 a 2013') 
ylabel('Vazões [m^3/s]') 
figure(2); 
plot(T,svaz0m); 87 
xlabel('Meses do ano') 
ylabel('Vazões [m^3/s]') 
figure(3); 
plot(T,vcalc); 
xlabel('Meses do ano') 
ylabel('Vazões [m^3/s]') 
x = sym('x'); 
q=z(1)*x.^3+z(2)*x.^2+z(3)*x+z(4); 
a0 = 1/11*int(q,x,1,12); 
w0 = 2*pi/11; 
nn = 4; 
i=1; 
for n=1:nn; 
an(i)=(2/11).*int((q).*cos(n*w0.*x),1,12); 
bn(i)=(2/11).*int((q).*sin(n*w0.*x),1,12); 
cn(i)=sqrt((an(i)^2)+(bn(i)^2)); 
qn(i)=-atan(bn(i)/an(i)); 
i=i+1; 
end; 
figure(4); 
cnv=eval(cn) 
bar(cnv); %Plotagem do gráfico de barras das amplitudes de vazões em cada 
periodo estudado (anual, semestral, quadrimestral e trimestral) 
xlabel('n') 
ylabel('Variação da vazão [m^3/s]') 
FO=[]; 
for n=1:nn 
87 
 
aa=1; 
for x = 0:0.5:24; 
s(aa)=aa; 
fo(aa)=cn(n)*cos(n*w0*x+qn(n)); 
vfo(n,aa)=eval(fo(aa)); 
aa=aa+1; 
end; 
FO=[FO;fo]; %Armazenamento do valor encontrado para cada n 
end; 
figure (5) 
subplot (4,1,1); plot (s,FO(1,:),'r-'); 
ylabel('[m^3/s]') 
subplot (4,1,2); plot (s,FO(2,:),'r-'); 
ylabel('[m^3/s]') 
subplot (4,1,3); plot (s,FO(3,:),'r-'); 
ylabel('[m^3/s]') 
subplot (4,1,4); plot (s,FO(4,:),'r-'); 
xlabel('Periodo de dois anos') 
ylabel('[m^3/s]') 
a0=eval(a0) 
cn=eval(cn) 
qn=eval(qn)