Xlabel('Время (с)');
ylabel('X1(t)')
n = 10;%Постоянная времени фильтра ФНЧ1(Апертура)=10дискрет
Wn = 0.00001;
b = fir1(n,Wn);%Фильтр 10порядка
y1=filter(b,1,x1);
%subplot(2,1,2)
figure (2);
plot(t,y1),grid,set(gca,'FontName','Arial Cyr','FontSize',16)
title('Выходной процесс на выходе ФНЧ1');
Xlabel('Время (с)');
ylabel('Y1(t)')
r = zeros(1,100);
for k = 1:100 %Задержка сигнала Nz=0:99
y = [y1(k:end) zeros(1,k-1)];
y = y.*y1;
r(k) = sum(y)/length(y);
end
k =30;%КРАСНАЯ
y = [y1(k:end) zeros(1,k-1)];
y = y.*y1;
n = 100;%Постоянная времени фильтра ФНЧкорр(Апертура)=100дискрет
Wn = 0.00001;%частота среза
b = fir1(n,Wn);%Делаем КИХ фильтр(методом взвешивания)-Фильтр 100порядка
z=filter(b,1,r);
%subplot(3,1,3)
figure (3);
stem(r),grid,set(gca,'FontName','Arial Cyr','FontSize',16)
title('Автокорреляционая функция');
hold on
Xlabel('Время задержки tau');
ylabel('K(tau)')
xlim([0 100])
%ylim([-5 5])
Графики
Рис.1.
Рис.2.
Рис.3.
Рис.4.
Рис.5.
Вывод коэффициентов корреляции
r =
Columns 1 through 4
0.2509 0.2375 0.2017 0.1530
Columns 5 through 8
0.1031 0.0611 0.0314 0.0136
Columns 9 through 12
0.0046 0.0009 -0.0002 -0.0002
Columns 13 through 16
0.0006 0.0017 0.0029 0.0041
Columns 17 through 20
0.0050 0.0054 0.0053 0.0048
Columns 21 through 24
0.0039 0.0027 0.0016 0.0006
Columns 25 through 28
0.0002 0.0006 0.0018 0.0038
Columns 29 through 32
0.0063 0.0088 0.0109 0.0119
Columns 33 through 36
0.0118 0.0106 0.0087 0.0068
Columns 37 through 40
0.0052 0.0042 0.0036 0.0033
Columns 41 through 44
0.0028 0.0018 0.0003 -0.0018
Columns 45 through 48
-0.0044 -0.0071 -0.0096 -0.0112
Columns 49 through 52
-0.0117 -0.0108 0 0
Columns 53 through 56
0 0 0 0
Columns 57 through 60
0 0 0 0
Columns 61 through 64
0 0 0 0
Columns 65 through 68
0 0 0 0
Columns 69 through 72
0 0 0 0
Columns 73 through 76
0 0 0 0
Columns 77 through 80
0 0 0 0
Columns 81 through 84
0 0 0 0
Columns 85 through 88
0 0 0 0
Columns 89 through 92
0 0 0 0
Columns 93 through 96
0 0 0 0
Columns 97 through 100
0 0 0 0