Определеение максимума и интеревала между максимумами

@Sergei_89 · Регистрация: 18.03.2013

Студворк — интернет-сервис помощи студентам

Добрый день всем!

Решается следующая задача.
Имеются несколько результатов измерений переведенные из xml в txt. Оттуда вы извлекаем эту информацию путем import data... Далее присваиваем имена для данных из файла Probe1_1 соответственно probe1_1, а из файла Probe2_1 соответственно probe2_1.

Далее воспользовавшись след кодом:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
clc;
n=0; Selection1=0
for i=1:length(probe1_1)
    if probe1_1(i)>0
        n=n+1;
        Selection1(n)=probe1_1(i);
    else
        n=n+1;
        Selection1(n)=probe1_1(i)*(-1);
    end
end
[b,a]=butter(3,0.0002);
y1=filter(b,a, Selection1);
plot(y1);

У нас есть некий график... точнее три графика, но это не столь важно. Если взять данные из второго файла с именем Probe2_1 то примерно такие же графики получаем, только с двумя пиками.

Вопрос собственно заключается вот в чем:
Каким образом осуществить поиски пика??? Обычный перебор не особо помогает. Смысл в том, что данные у нас в форме <Nx1> представлены, где N - количество элементов, порядка 60000. К сожалению txt фаил не получается загрузить, ибо он слишком велик. Могу выслать на почту или еще каким-то способом.

Знаю, что есть еще вариант определения перегиба функции по второй производной, но тут мне кажется это не самый лучший вариант. Я нашел некую функцию которая позволяет находить вторую производную, но как-то она не внушает доверия, если честно.

Может быть я чего-то просто не замечаю, глаз, так сказать, замылился.

@vital792 · 18.03.2013, 14:24

Сообщение от Sergei_89

есть еще вариант определения перегиба функции по второй производной

не пойму, зачем вам перегиб по второй производной. Судя по описанию вам скорее надо найти максимум по первой производной, приравнивая ее значение нулю. Но так как мы все используем дискретные сигналы, понятие первой производной тут сильно упрощается, вам достаточно пройти в цикле по выборке проверить условия

Matlab M
1
y1(k)>y1(k-1) && y1(k)>y1(k+1)

и собрать в массив пары [значение положение] точек удовлетворяющих условию. Если при этом появятся лишние точки маленького уровня, можно срезать их задав пороговое значение, либо уменьшить полосу фильтра. Кстати не забываете, что когда вы фильтруете, имеется групповая(или фазовая, не помню как правильно) задержка. Так что полученные максимумы будут смещены на длину импульсной характеристики фильтра.

@Sergei_89 · 19.03.2013, 12:02 **[ТС]**

Спасибо за совет.

Но на самом деле все оказалось в принципе еще проще

Matlab M
1
2
3
4
5
6
7
8
maximum=0.0002; %задаем пороговое значение ниже которого результаты не анализруются
 for i=1:length(y1)
     if y1(i)>0.0002
         if y1(i)>maximum
            maximum=y1(i);
         end
     end
 end

Новые блоги и статьи Все статьи Все блоги /
Оттенки серого Argus19 18.03.2026 Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .	SDL3 для Desktop (MinGW): Рисуем цветные прямоугольники с помощью рисовальщика SDL3 на Си и C++ 8Observer8 17.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-rectangles-sdl3-c. zip finish-rectangles-sdl3-cpp. zip	Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая ссылка» (hard link),. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .