Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement