Диаграмма Боде корреляции двух сигналов

%%  Data upload
Data = xlsread('Data.xls');
time = Data(:,1);
A = Data(:,2);
B = Data(:,3);
 
%% Cross-correlation
e = fft(A);
f = fft(B);
 
g = e.*conj(f);
 
h = ifft(g); % cross correlation coefficiens
 
%% Plots
figure('Name','Corr Coef')
plot(time, h)
grid on
xlabel('Time')
ylabel('Amp')
 
% Amp = 10*log10(h); % Перевод в Db
 
F = (1:length(e))*(1/(time(2)-time(1)))/length(e); %Частотная сетка
 
figure('Name','Amp')
plot(F,abs(g))
% semilogx(F,abs(h))
% semilogx(F,Amp)
%xlim([0 10])      
grid on
xlabel('Frequency [Hertz]')
ylabel('Amp [???]')
 
phase = angle(g)*180/pi;
 
figure('Name','Phase')
plot(F,phase)
% xlim([0 10])
grid on
xlabel('Frequency [Hertz]')
ylabel('Phase [Degrees]')

h = xcorr(B,A,'coeff')

Data = xlsread('Data.xls');
time = Data(:,1);
A = Data(:,2);
B = Data(:,3);
tr = time(2)-time(1);
 
opts = bodeoptions('cstprefs');
opts.FreqScale = 'linear';
opts.MagScale = 'linear';
opts.PhaseUnits = 'deg';
opts.MagUnits = 'abs';
opts.Grid = 'on';
opts.Xlim = [0 10];
opts.Title.FontSize = 14;
opts.XLabel.FontSize = 14;
opts.YLabel.FontSize = 14;
 
DAT = iddata(A,B,tr);   %time domain
sys = tfest(DAT, 2);
figure('Name','Bode1 (tfest)')
bodeplot(sys, opts);
 
FDAT = fft(DAT);    %frequency domain
sys1 = etfe(FDAT, 128);
figure('Name','Bode1 (etfe)')
bodeplot(sys1, opts);