Лаба 7 / sinsoid

clc
close all
clear all

Fs=350;
tmax=0.8;
F=42;
T=1/Fs;

t=0:T:tmax-T;
M=5;
A(1:M)=[1 1 1 1 1];
C(1:M)=[0 0 0 0 0];
for i=1:M
Signal(:,i)=A(i)*sin(2*pi*F*t(:)-(i-1)*pi/4)+C(i);
end

figure
subplot(2,1,1)
dy=3;
for i=1:M
plot(t,Signal(:,i)-i*dy)
hold on
end
title('Графики тестовых сигналов')

subplot(2,1,2)
R=corrcoef(Signal);
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')