clear all;clc;
[sinal,Fs] = audioread('sinal_fsk.wav'); %% as informações do sinal que esta contido no arquivo.wav
N = length(sinal); %lenght é o tamanho do sinal
ts = 1/Fs; %%tempo de amostragem

%%Calculando a FFT do sinal (função computacional mais eficiente para o
%%cálculo da DFT)

y = fft(sinal);
y = y/N;

%%Cálculo do eixo das frequências
m = 0:(N-1);
f = m*Fs/N;%%resolução
y = y(1:N/2); %pegando metade dos pontos
f = f(1:N/2);

%%Espectros em escala linear e log
magnitude = abs(y); %Magnitude em escala linear
fase = angle(y);
f_ref = max(magnitude); %escolhe o maior valor para ser a referencia para normalizacao
y_db = -20*log10(magnitude/f_ref); %Magnitude em db (- indicativo para inverter)
fase_db = -20*log10(fase/(2/pi));

%%Plotando
figure(1);
subplot(4,1,1); stem(f,magnitude); title('Magnitude do espectro escala linear'); xlabel('freq(Hz)');ylabel('Amplitude');
subplot (4,1,2); stem(f,fase,'r'); title('Fase');
subplot(4,1,3);stem(f,y_db); title('Potencia do espectro'); xlabel('m'); ylabel('dB');
subplot(4,1,4);stem(f,fase_db,'r'); title('Fase em dB');