%% 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
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
