Работы / ПЗ №8 Оифс
.docxФедеральное государственное бюджетное образовательное учреждение
высшего образования
«Петербургский государственный университет путей сообщения
Императора Александра I»
Кафедра «Методы и приборы неразрушающего контроля»
Отчет по практическому заданию № 8:
«Корреляционный анализ. Определение задержки между двумя сигналами»
Выполнил:
студент группы ПБ-713
Сабитов И.А.
Проверила:
Чурова В.В.
Санкт-Петербург
2020 г.
Цель работы: изучение метода вычисления относительной разности времени распространения ультразвуковых волн.
Исходные данные.
-
№ варианта
Значение file
8
file8.mat
Порядок выполнения работы.
В данной работе необходимо определить время задержки между двумя эхосигналами, который выполняется на основе взаимокорреляционной функции (ВКФ). Требуется ее вычислить.
Текст программы в MatLab представлен ниже:
load('file8.mat'); % загрузка исходных данных
c = xcorr( fx, fy ); % расчет ВКФ
[x, index] = max( abs(c) ); % поиск максимума ВКФ
offset = (index - (length(c)+1) / 2) * dt * 1e9; % расчет задержки
h = hilbert(c); % огибающая ВКФ
Полученные результаты максимума для модуля ВКФ и времени задержки:
С помощью преобразования Гильберта получили огибающую ВКФ и построили совместные графики двух сигналов и ВКФ с огибающей, используя код программы, представленный ниже:
clc;
clear;
load('file8'); % загрузка исходных данных
c = xcorr( fx, fy ); % расчет ВКФ
[x, index] = max( abs(c) ); % поиск максимума ВКФ
offset = (index - (length(c)+1) / 2) * dt * 1e9; % расчет задержки
h = hilbert(c); % огибающая ВКФ
subplot(1, 2, 1)
plot(c)
grid
hold on
plot(abs(h))
hold off
subplot(1, 2, 2)
plot(fx)
grid
hold on
plot(fy)
Результатом работы программы будет график (рисунок 1).
Рисунок 1