Лаба 3 / eeg

%eeg
clc;
clear all;
close all;

X0=50;
Y0=80;
W0=1200;
H0=600;

x1=50;
y1=430;
w1=200;
h1=150;
dy=30;
dx=20;
w2=40;
h2=10;

figure('Color','[1 1 1]','Position',[X0,Y0,W0,H0]);
hAxes11=axes('Units','pixels','Position',[x1,y1,w1,h1]);
hAxes21=axes('Units','pixels','Position',[x1,y1-h1-dy,w1,h1]);
hAxes31=axes('Units','pixels','Position',[x1,y1-2*h1-2*dy,w1,h1]);
hAxes12=axes('Units','pixels','Position',[x1+w1+dx,y1,2*w1,h1]);
hAxes22=axes('Units','pixels','Position',[x1+w1+dx,y1-h1-dy,2*w1,h1]);
hAxes32=axes('Units','pixels','Position',[x1+w1+dx,y1-2*h1-2*dy,2*w1,h1]);
hAxes13=axes('Units','pixels','Position',[x1+3*w1+2*dx,y1,2*w1,h1]);
hAxes23=axes('Units','pixels','Position',[x1+3*w1+2*dx,y1-h1-dy,2*w1,h1]);
hAxes33=axes('Units','pixels','Position',[x1+3*w1+2*dx,y1-2*h1-2*dy,2*w1,h1]);

%параметры сигнала
Fs=250;
S=load('R3_07.txt');
N=length(S);
S1=S(1:N,1);
S2=S(1:N,2);
S3=S(1:N,3);

T=1/Fs;
tmax=N*T;

t=0:T:tmax-T;

%пределы
Amax(1)=max(S1);
Amax(2)=max(S2);
Amax(3)=max(S3);

Amin(1)=min(S1);
Amin(2)=min(S2);
Amin(3)=min(S3);

%графики
axes(hAxes11);
plot(t,S1);
grid minor;
ylim([min(Amin) max(Amax)]);
hold on;

axes(hAxes21);
plot(t,S2);
grid minor;
ylim([min(Amin) max(Amax)]);
hold on;

axes(hAxes31);
plot(t,S3);
grid minor;
ylim([min(Amin) max(Amax)]);
hold on;

%АКФ
tcf=-tmax+T:T:tmax-T;

acf1=xcorr(S1,'coeff');
axes(hAxes12);
hold on;
plot(tcf,acf1);
grid minor;
%ylim([min(Amin) max(Amax)]);
hold on;

acf2=xcorr(S2,'coeff');
axes(hAxes22);
hold on;
plot(tcf,acf2);
grid minor;
%ylim([min(Amin) max(Amax)]);
hold on;

acf3=xcorr(S3,'coeff');
axes(hAxes32);
hold on;
plot(tcf,acf3);
grid minor;
%ylim([min(Amin) max(Amax)]);
hold on;

%ВКФ
ccf13=xcorr(S1,S2,'coeff');
axes(hAxes13);
hold on;
plot(tcf,ccf13);
grid minor;
hold on;

ccf23=xcorr(S2,S3,'coeff');
axes(hAxes23);
hold on;
plot(tcf,ccf23);
grid minor;
hold on;

ccf31=xcorr(S3,S1,'coeff');
axes(hAxes33);
hold on;
plot(tcf,ccf31);
grid minor;
hold on;