%% INFORME DEL PROCESADO FRECUENCIAL
%
%% Lectura de parámetros
%
load frecuencia

dens=Pyy;
densnor=PyyNorm;
denspor=PyyPorc;
n_ptos=num_puntos;
tm=tmuestreo;
n_canales=canales;
V=[];
cadena='Canal';

%% Representación de la Densidad Espectral de Potencia de los Canales
%
f=(0:n_ptos/2-1)/(tm*n_ptos);

for i=1:n_canales
    V=[V;strcat(cadena,num2str(i))];
end
plot(f,dens(1:n_ptos/2,3:2+n_canales));
title('Representación de la Densidad Espectral de Potencia')
xlabel('f (Hz)')
ylabel('Dens. Espectral de Pot.')
legend(V)

%% Representación de la Densidad Espectral de Potencia de la Temperatura
%
f=(0:n_ptos/2-1)/(tm*n_ptos);

plot(f,dens(1:n_ptos/2,1));
title('Representación de la Densidad Espectral de Potencia')
xlabel('f (Hz)')
ylabel('Dens. Espectral de Pot.')
legend('Temperatura')

%% Representación de la Dens. Espectral de Potencia de la Humedad
%
f=(0:n_ptos/2-1)/(tm*n_ptos);

plot(f,dens(1:n_ptos/2,2));
title('Representación de la Densidad Espectral de Potencia')
xlabel('f (Hz)')
ylabel('Dens. Espectral de Pot.')
legend('Humedad')


%% Representación de la Dens. Espectral de Pot. normalizada de los Canales
%
figure;
plot(f,densnor(1:n_ptos/2,3:2+n_canales));
title('Representación de la Densidad Espectral de Potencia Normalizada')
xlabel('f (Hz)')
ylabel('Dens. Espectral de Pot. Normalizada')
legend(V)

%% Representación de la Dens. Espectral de Pot. normalizada de Temperatura
%
figure;
plot(f,densnor(1:n_ptos/2,1));
title('Representación de la Densidad Espectral de Potencia Normalizada')
xlabel('f (Hz)')
ylabel('Dens. Espectral de Pot. Normalizada')
legend('Temperatura')


%% Representación de la Dens. Espectral de Pot. normalizada de la Humedad
%
figure;
plot(f,densnor(1:n_ptos/2,2));
title('Representación de la Densidad Espectral de Potencia Normalizada')
xlabel('f (Hz)')
ylabel('Dens. Espectral de Pot. Normalizada')
legend('Humedad')


%% Dens. Espectral de Pot. respecto al total de Potencia de los Canales
%
figure;
plot(f,denspor(1:n_ptos/2,3:2+n_canales));
title('Representación de la Densidad Espectral de Potencia respecto al total de Potencia')
xlabel('f (Hz)')
ylabel('Dens. Espectral de Pot (%)')
legend(V)


%% Dens. Espectral de Pot. respecto al total de Potencia de la Temperatura
%
figure;
plot(f,denspor(1:n_ptos/2,1));
title('Representación de la Densidad Espectral de Potencia respecto al total de Potencia')
xlabel('f (Hz)')
ylabel('Dens. Espectral de Pot (%)')
legend('Temperatura')

%% Dens. Espectral de Pot. respecto al total de Potencia de la Humedad
%
figure;
plot(f,denspor(1:n_ptos/2,2));
title('Representación de la Densidad Espectral de Potencia respecto al total de Potencia')
xlabel('f (Hz)')
ylabel('Dens. Espectral de Pot (%)')
legend('Humedad')


%% Cálculo de las principales estádisticas
%%
% 
% * *Potencia*
%%
%   Código que permite obtener la potencia para cada canal y sus
%   correspondientes valores
for i=1:n_canales+2
    if i==1
        fprintf('Potencia de la Temperatura : %7.3f (W)\n',Potencia(i))  
    elseif i==2
        fprintf('Potencia de la Humedad : %7.3f (W)\n',Potencia(i))
    else
        fprintf('Potencia del canal %d : %7.3f (W)\n',i-2,Potencia(i))
    end
end

%%
% 
% * *Valor máximo de la Densidad espectral de Potencia*
%%
%   Código que permite obtener el valor máximo para cada canal y sus
%   respectivos valores

for i=1:n_canales+2
    if i==1  
        fprintf('Valor máximo de la Dens. Espectral de Pot para la Temperatura: %7.3f\n',maximoPotencia(i))
    elseif i==2
        fprintf('Valor máximo de la Dens. Espectral de Pot para la Humedad: %7.3f\n',maximoPotencia(i))
    else
        fprintf('Valor máximo de la Dens. Espectral de Pot para el canal %d: %7.3f\n',i-2,maximoPotencia(i))
    end
end

%%
% 
% * *Valor mínimo de la Densidad espectral de Potencia*
%%
%   Código que permite obtener el valor mínimo para cada canal y sus
%   respectivos valores
for i=1:n_canales+2
    if i==1
      fprintf('Valor mínimo de la Dens. Espectral de Pot para la Temperatura: %7.3f\n',minimoPotencia(i))
    elseif i==2
        fprintf('Valor mínimo de la Dens. Espectral de Pot para la Humedad: %7.3f\n',minimoPotencia(i))
    else
        fprintf('Valor mínimo de la Dens. Espectral de Pot para el canal%d: %7.3f\n',i-2,minimoPotencia(i))
    end
end


%%
% 
% * *Máximos locales de la Densidad Espectral de Potencia*
%%
%   Código que permite obtener el valor de los máximos locales de la
%   densidad espectral de potencia
x1=[];
valor=[];

for i=1:n_canales+2
    x1=find(localmaxPot(:,i));
    valor=localmaxPot(:,i);
    if i==1
        fprintf('Valores máximos de la Dens. Esp. de Potencia para la Temperatura\n')
    elseif i==2
        fprintf('Valores máximos de la Dens. Esp. de Potencia para la Humedad\n')
    else
        fprintf('Valores máximos de la Dens. Esp. de Potencia para el canal %d\n',i-2)
    end
    for j=1:length(x1)
        fprintf('Valor: %8.4f\n',valor(x1(j)))
    end
    fprintf('\n')
end


%%
% 
% * *Mínimos locales de la Densidad Espectral de Potencia*
%%
%   Código que permite obtener el valor de los mínimos locales de la
%   densidad espectral de Potencia
x1=[];
valor=[];

for i=1:n_canales+2
    x1=find(localminPot(:,i));
    valor=localminPot(:,i);
    if i==1
        fprintf('Valores mínimos de la Dens. Esp. de Potencia para la Temperatura\n')
    elseif i==2
        fprintf('Valores mínimos de la Dens. Esp. de Potencia para la Humedad\n')
    else
        fprintf('Valores mínimos de la Dens. Esp. de Potencia para el canal %d\n',i-2)
    end
    for j=1:length(x1)
        fprintf('Valor: %8.6f\n',valor(x1(j)))
    end
    fprintf('\n')
end




       


 














