From 19a477eddbfee6a1b0d0d6374278d82c669a0907 Mon Sep 17 00:00:00 2001
From: Diogo Cordeiro <diogo@fc.up.pt>
Date: Thu, 26 Nov 2020 15:14:39 +0000
Subject: [PATCH] First step towards TP3_2 and 3

---
 TP3a/DiogoEliseu_TP3_2.m | 54 ++++++++++++++++++++++++++++++++++++++++
 TP3a/DiogoEliseu_TP3_3.m | 25 +++++++++++++++++++
 2 files changed, 79 insertions(+)
 create mode 100644 TP3a/DiogoEliseu_TP3_3.m

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