Лаба 7 / realsignal

clc
clear all
close all

EEG=load('R7_04.txt');
N=length(EEG(:,1));
M=length(EEG(1,:));
Fs=185;
T=1/Fs;
tmax=N*T;
t=0:T:tmax-T;

figure
subplot(2,1,1)
dy=100;
for i=1:M
plot(t,EEG(:,i)-i*dy)
hold on
end
title('Графики 16-канальной ЭЭГ')

subplot(2,1,2)
R=corrcoef(EEG);
Ra=zeros(M+1,M+1);
Ra(1:M,1:M)=R;
pcolor(Ra)
colormap(gray)
axis ij
axis square
title('КВК ЭЭГ')

figure
subplot(1,2,1)
Limit=0.5;
for i=1:M
for j=1:M
if R(i,j)>Limit
R1(i,j)=1;
else
R1(i,j)=0;
end
end
end
Ra1=zeros(M+1,M+1);
Ra1(1:M,1:M)=R1;
pcolor(Ra1)
colormap(gray)
axis ij
axis square
title('КВК ЭЭГ, порог = 0.5')

subplot(1,2,2)
Limit=0.85;
for i=1:M
for j=1:M
if R(i,j)>Limit
R2(i,j)=1;
else
R2(i,j)=0;
end
end
end
Ra2=zeros(M+1,M+1);
Ra2(1:M,1:M)=R2;
pcolor(Ra2)
colormap(gray)
axis ij
axis square
title('КВК ЭЭГ, порог = 0.85')