diogo/DSP
Archived
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

First step towards TP3_2 and 3

Browse Source
This commit is contained in:
Diogo Cordeiro 2020-11-26 15:14:39 +00:00
parent 2ae24399f1
commit 19a477eddb
2 changed files with 79 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

54
TP3a/DiogoEliseu_TP3_2.m
View File

@ -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

25
TP3a/DiogoEliseu_TP3_3.m Normal file
View File

@ -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