Разрезать ультразвуковой сигнал между конкретными значениями с помощью октавы - PullRequest
Новая
       22

Разрезать ультразвуковой сигнал между конкретными значениями с помощью октавы

1 голос
/ 04 октября 2019

У меня ультразвуковая волна (ось графика: Вольт против микросекунды), и мне нужно разрезать сигнал / волну между двумя конкретными значениями для дальнейшего анализа этого отсечения. Моя идея - обрезать сигнал между 0,2 В (ось Y). Волна синусоидальной формы , как показано на рисунке, с желаемыми точками среза красным

В моем текущем коде я обрезаю сигнал в интервале от 1900 до 4000 мс (ось x). (Aa = A(1900:4000);), а затем я хочу сделать вышеупомянутое отсечение и продолжить работу с кодом.

Кто-нибудь знает, как я мог бы сделать это отсечение по оси Y?

Спасибо !! :) 

clear
clf
pkg load signal
for k=1:2
  w=1
  filename=strcat("PCB 2.1  (",sprintf("%01d",k),").mat") 
  load(filename)
  Lthisrun=length(A);
  Pico(k,1:Lthisrun)=A;


  Aa = A(1900:4000);



  Ah= abs(hilbert(Aa));

  step=100;
  hold on
i=1;
Ac=0;
      for index=1:step:3601
          Ac(i+1)=Ac(i)+Ah(i);  
          i=i+1
          r(k)=trapz(Ac)
      end
end

1 Ответ

0 голосов
/ 07 октября 2019

хорошо, вы хотите просто посмотреть на значения «выше уровня шума» в ваших данных. Или, в этом случае, «обрезать» все, что ниже 0,2 В. Самый простой способ сделать это с помощью логической индексации. Вы можете взять массив и создать подмассив, исключив все, что не соответствует определенному логическому условию. Смотрите этот пример: 

f = @(x) sin(x)./x;
x = [-100:.1:100];
y = f(x);
plot(x,y);

plot of sin(x)/x

figure;
x_trim = x(y>0.2);
y_trim = y(y>0.2);
plot(x_trim, y_trim);

image 0.2">

От вашего вопроса этоПохоже, вы хотите сделать отсечение после применения горизонтального окна от 1900-4000. (вы говорите, что это в миллисекундах, но ваше изображение показывает, что импульс гораздо раньше, чем 1900 мс). В любом случае, что-то вроде 

Ab = Aa(Aa > 0.2);

создаст другой массив Ab, который будет содержать только части Aa со значениями выше 0,2. Возможно, вам придется сделать нечто подобное (см. Пример) для горизонтальной оси, если ваши x-данные не просто индекс элемента.

...