Корреляционная функция

[agr0sha]

Здравствуйте.
решил проанализировать корреляционную функцию записанных частот (wav файл прилагаю) РДНС-20 "Альфа", использовал MATLAB для наиболее быстрого способа получения результата, написал код и получил график (идет во вложении), подскажите почему сигнал на частоте 14881 ГГц имеет такой всплеск?

Вложения

lr2Y5gS6hHU.jpg

(72.01 KiB) Скачиваний: 317

[КРАМ]

Может речь идет об АВТОкорреляционной функции?
ЗЫ. ГГц полагаю опечаткой...?

[agr0sha]

Нет, функция корреляционная, но суть не в этом
С частотой и вправду недоразумение вышло, опечатка
Вопрос всё так же остается актуальным

[HochReiter]

Нет, функция корреляционная

Кого с кем?

[mickbell]

Вероятно, синей и зелёной кривулин.

[КРАМ]

agr0sha писал(а):

Нет, функция корреляционная, но суть не в этом

Тогда речь идет о том, что две первых частоты коррелированы с третьей, а третья сама с собой....

[agr0sha]

Вот программный код

clc
clear
close all
[y, Fs, nbits] = wavread('sample.wav');
D = y(:,1);
L = length(D);

f1 = 11905;
f2 = 12649;
f3 = 14881;

t(:,1)=1/Fs:1/Fs:((L)/Fs);
D1 = sin(2*pi*f1*t).*(y(:,1));
D2 = sin(2*pi*f2*t).*(y(:,1));
D3 = sin(2*pi*f3*t).*(y(:,1));

CORR = zeros(L,3);
z1 = 0;
z2 = 0;
z3 = 0;
for i=1:L
z1 = z1 + D1(i);
z2 = z2 + D2(i);
z3 = z3 + D3(i);

CORR(i,1) = abs(z1);
CORR(i,2) = abs(z2);
CORR(i,3) = abs(z3);
end
plot(t,CORR)
xlabel('t,c')
ylabel('Корреляционнонная функция')
legend('11905','12649','14881',0)

[КРАМ]

Ну и нахрена нам эта портянка?
Просто формулы написать лень? Или Вы плохо понимаете что считает этот код?
Тем более, что никто не знает что там в этом sample.wav