diff --git a/TP3a/DiogoEliseu_TP3_2.m b/TP3a/DiogoEliseu_TP3_2.m index e69de29..60a9c68 100644 --- a/TP3a/DiogoEliseu_TP3_2.m +++ b/TP3a/DiogoEliseu_TP3_2.m @@ -0,0 +1,54 @@ +%% Inicialização do ambiente +clear ; close all; clc + +%% Exercício 2 +sintetico_data = fopen('sintetico.csv'); +sintetico = textscan(sintetico_data, '%f', 'Delimiter', ','); +sintetico = sintetico{1,1}; +fclose(sintetico_data); +cello = importdata('cello.mat'); +EEG = importdata('EEG.mat'); +EEG = EEG'; + +N = 256; + +figure(1) + +subplot(2,1,1) +plot(sintetico); +title('sinal sintetico'); + +subplot(2,1,2) +DFT_Xa = fft(sintetico, N+1); +DFT_Xa = DFT_Xa(1 : (N+1)/2); +plot(X((N+1)/2, 0, 1),abs(DFT_Xa)) +xlabel('\omega/\pi'), ylabel('|X(\omega)|'), title('DFT sintetico em [0, \pi]') + +figure(2) + +subplot(2,1,1) +plot(cello.x); +title('sinal cello'); + +subplot(2,1,2) +DFT_cello = fft(cello.x, N+1); +DFT_cello = DFT_cello(1 : (N+1)/2); +plot(X((N+1)/2, 0, 1),abs(DFT_cello)) +xlabel('\omega/\pi'), ylabel('|X(\omega)|'), title('DFT cello em [0, \pi]') + +figure(3) + +subplot(2,1,1) +plot(EEG); +title('sinal EEG'); + +subplot(2,1,2) +DFT_EEG = fft(EEG, N+1); +DFT_EEG = DFT_EEG(1 : (N+1)/2); +plot(X((N+1)/2, 0, 1),abs(DFT_EEG)) +xlabel('\omega/\pi'), ylabel('|X(\omega)|'), title('DFT EEG em [0, \pi]') + + +function f = X(N, inicio, fim) +f = (0:(1/N):(1 - 1/N))*(fim-inicio)+inicio; +end \ No newline at end of file diff --git a/TP3a/DiogoEliseu_TP3_3.m b/TP3a/DiogoEliseu_TP3_3.m new file mode 100644 index 0000000..619038b --- /dev/null +++ b/TP3a/DiogoEliseu_TP3_3.m @@ -0,0 +1,25 @@ +%% Inicialização do ambiente +clear ; close all; clc + +%% Exercício 3 +sintetico_data = fopen('sintetico.csv'); +sintetico = textscan(sintetico_data, '%f', 'Delimiter', ','); +sintetico = sintetico{1,1}; +fclose(sintetico_data); + +figure(1) + +subplot(2,1,1) +plot(sintetico); +title('sinal sintetico'); + +subplot(2,1,2) +DFT_Xa = fft(sintetico, N+1); +DFT_Xa = DFT_Xa(1 : (N+1)/2); +plot(X((N+1)/2, 0, 1),abs(DFT_Xa)) +xlabel('\omega/\pi'), ylabel('|X(\omega)|'), title('DFT sintetico em [0, \pi]') + + +function f = X(N, inicio, fim) +f = (0:(1/N):(1 - 1/N))*(fim-inicio)+inicio; +end \ No newline at end of file