RU2787309C1 - Способ идентификации мультисинусоидальных цифровых сигналов - Google Patents
Способ идентификации мультисинусоидальных цифровых сигналов Download PDFInfo
- Publication number
- RU2787309C1 RU2787309C1 RU2022112232A RU2022112232A RU2787309C1 RU 2787309 C1 RU2787309 C1 RU 2787309C1 RU 2022112232 A RU2022112232 A RU 2022112232A RU 2022112232 A RU2022112232 A RU 2022112232A RU 2787309 C1 RU2787309 C1 RU 2787309C1
- Authority
- RU
- Russia
- Prior art keywords
- signal
- model
- multisinusoidal
- block
- values
- Prior art date
Links
- 238000001914 filtration Methods 0.000 claims abstract description 10
- 238000001514 detection method Methods 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 15
- 239000000126 substance Substances 0.000 abstract 1
- 230000001131 transforming Effects 0.000 description 19
- 239000011159 matrix material Substances 0.000 description 16
- 238000005259 measurement Methods 0.000 description 9
- 238000000844 transformation Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 210000004544 DC2 Anatomy 0.000 description 3
- NLKNQRATVPKPDG-UHFFFAOYSA-M Potassium iodide Chemical compound [K+].[I-] NLKNQRATVPKPDG-UHFFFAOYSA-M 0.000 description 3
- 230000002530 ischemic preconditioning Effects 0.000 description 3
- 230000000737 periodic Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000007374 clinical diagnostic method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000000875 corresponding Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral Effects 0.000 description 1
Images
Abstract
Предлагаемым изобретением ставится задача идентификации мультисинусоидального цифрового сигнала, позволяющая осуществить фильтрационную обработку сигнала с помощью структурного преобразователя, в идентификаторе непрерывной цепной С-дроби вести подсчет числа гармонических составляющих сигнала, оценивать параметры сигнала (амплитуда, частота) поэтапно по мере убывания значений их амплитуд в восстановителе модели гармоники с последующим восстановлением остатка сигнала, определять наличие смещения и восстанавливать его значение, расширяя тем самым область применения метода идентификации для различных типов мультисинусоидальных сигналов, автоматически определяя структуру и параметры полезного сигнала. Техническим результатом при реализации заявленного решения является повышение точности, достоверности и быстродействия процесса идентификации, что позволяет использовать данный метод как способ текущей идентификации в реальном масштабе времени, определять изменение структуры модели и тем самым повышать точность, надежность и качество результатов. 8 ил.
Description
Изобретениие относится к области обработки информации и измерительной техники и может быть использовано для контроля работоспособности электротехнических и электромеханических устройств. Способ может быть применен для выделения полезного мультисинусоидального сигнала и для определения его математической модели на основе цифровых отсчетов, и реализован с использованием ЭВМ в автоматическом режиме, в реальном масштабе времени.
Известен СПОСОБ ВЫДЕЛЕНИЯ ПОЛЕЗНОГО СИГНАЛА РЕАЛИЗУЕМОГО ПРОЦЕССА (патент РФ №2394216, G01M 7/00, G01R 23/16, опубл. 10.07.2010 http://www.freepatent.ru/patents/2394216), основанный на узкополостной фильтрации сигнала. Узкополосную фильтрацию этого сигнала путем автоматической частотной настройки соответствующего фильтра проводят с последующим детектированием отфильтрованного сигнала и его сравнением с пороговой величиной, после чего формируют сигнал управления компьютерной обработкой процесса, а затем осуществляют широкополосную фильтрацию реализуемого процесса.
Данный способ имеет следующие недостатки:
- о спектральном составе сигнала можно судить только с помощью визуального анализа графиков, что не позволяет произвести автоматизацию расчетов и снижает быстродействие и точность метода;
- после фильтрации шумовой составляющей определяются только пороговые значения периодического сигнала, а число гармоник и их параметры (амплитуда, частота) остаются неизвестными.
Известен способ определение наличия гармонических составляющих и их частот в дискретных сигналах (патент РФ №2498324, МПК G01R 23/00, опубл. 10.11.2013 http://www.freepatent.ru/patents/2498324), использующий быстрое преобразование Фурье входных сигналов. Определяют комплексно-сопряженные значения прямого преобразования Фурье анализируемого дискретного сигнала, попарно умножают полученные комплексные сигналы прямого преобразования Фурье анализируемого дискретного сигнала, тем самым формируют m сигналов и подвергают их обратному преобразованию Фурье, определяют частотно-временную автокорреляционную функцию и строят ее график, по которому судят о наличии гармонических составляющих и их частотах. Данный способ авторы предлагают использовать для анализа сигналов различного происхождения в задачах неразрушающего контроля и диагностики оборудования, однако такие сигналы имеют в своем составе шумовые составляющие. Данный способ не предусматривает предварительную фильтрацию шумов, что может впоследствии привести к ошибочным оценкам структуры и параметров цифровых сигналов. Кроме того, данный способ имеет следующие недостатки, снижающие быстродействие и точность метода:
- большое количество операций для осуществления быстрого преобразования Фурье, формирование большого количества дополнительных сигналов (более 1000), вследствие чего наблюдается громоздкость расчетов и неточность вычислений;
- для реализации метода необходимо большое количество измерений сигнала (более 16000) и малый шаг дискретизации (менее 10-4);
- визуальный анализ графиков частотно-временной корреляционной функции для определения количества гармоник и частот сигнала, не позволяющий произвести автоматизацию расчетов;
- нет возможности оценить амплитуды гармонических составляющих.
Известен способ гармонического анализа периодического многочастотного сигнала (патент РФ №2435168, МПК G01R 23/16, опубл. 27.11.2011, http://www.freepatent.ru/patents/2435168), использующий итерационный процесс определения необходимой формы весовой функции. В результате многократных оценок частот периодического сигнала форма весовой функции подбирается такой, чтобы на частоте любого из анализируемых составляющих сигнала слагаемые спектра других составляющих сигнала по амплитуде были равны нулю, причем кратность нуля, определяемая порядком производных модуля спектра, может быть заданной степени. Амплитуду и фазу каждого составляющего сигнала оценивают после определения точных значений частот.
Данный способ имеет следующие недостатки, снижающие быстродействие и точность метода:
- итерационная процедура подбора весовой функции и проведение дополнительных циклов гармонического анализа требуют значительного числа операций, поэтому обладают недостаточной быстродействием и наличием погрешности расчетов;
- громоздкость и неточность нахождения Фурье-образа спектра;
- визуальный анализ лепестков спектра значительно затрудняет процедуру оценки частот сигнала и не позволяет автоматизировать расчеты;
- для реализации метода необходимо большое количество измерений сигнала (более 100 периодов колебаний сигнала).
Наиболее близким к предлагаемому способу является способ идентификации мультисинусоидальных цифровых сигналов (патент РФ 2703933, МПК G01R 23/16, опубл. 22.10.2019), включающий использование непрерывных цепных дробей путем измерения сигнала в равноотстоящие промежутки времени, подачи их на идентификатор непрерывной С-дроби с последующим восстановлением числа гармоник, круговых частот, амплитуд, прогнозирующей модели и модельных значений сигнала, отличающийся тем, что при обработке значений сигнала в идентификаторе непрерывной С-дроби на основе правила останова судят о наличии смещения сигнала, а при его выявлении выполняют последующее восстановление его значения и далее с учетом полученных результатов судят о начальных фазовых сдвигах в гармониках сигнала.
Этот способ может применяться только для не зашумленных мультисинусоидальных сигналов, что значительно снижает качество идентификации, поскольку в реальных системах контроля работоспособности электротехнических и электромеханических устройств измерительная информация обычно искажена помехами. Для мультисинусоидальных сигналов, содержащих шумовые составляющие, способ приводит к получению ошибочной модели полезного сигнала.
Предлагаемым изобретением ставится задача идентификации мультисинусоидального цифрового сигнала, позволяющая осуществить фильтрационную обработку сигнала с помощью структурного преобразователя, в идентификаторе непрерывной цепной С-дроби вести подсчет числа гармонических составляющих сигнала, оценивать параметры сигнала (амплитуда, частота) поэтапно по мере убывания значений их амплитуд в восстановителе модели гармоники с последующим восстановлением остатка сигнала, определять наличие смещения и восстанавливать его значение, расширяя тем самым область применения метода идентификации для различных типов мультисинусоидальных сигналов, автоматически определяя структуру и параметры полезного сигнала, значительным образом повышая степень автоматизации, точность, достоверность и быстродействие процесса идентификации, что позволяет использовать данный метод как способ текущей идентификации в реальном масштабе времени, определять изменение структуры модели и, тем самым, повышать точность, надежность и качество результатов моделирования.
Техническим результатом предлагаемого способа идентификации мультисинусоидальных цифровых сигналов является то, что обеспечивается быстродействие, универсальность реализации способа, простота и высокая точность вычислений в условиях наличия помех в сигналах.
Изобретение поясняется на фигурах 1-8.
Структурная схема системы, реализующая предлагаемый способ, изображена на фиг. 1. Она содержит датчик анализируемого сигнала (ДАС), с выхода которого сигнал x(kΔt) поступает на вход последовательно соединенного блока 1 - структурного преобразователя, в котором осуществляется предварительная фильтрационная обработка значений сигнала на основе трех последовательных структурных преобразований. Результаты третьего структурного преобразования сигнала последовательно передаются на вход блока 2 - идентификатора непрерывной цепной С-дроби, который включает в себя 2 подблока. В подблоке 2.1 счетчика гармоник рассчитывается матрица-идентификатор, в которой 4-ая строка будет являться нулевой, и далее запускается подсчет числа гармоник i=i+1 (первоначально i=0). В подблоке 2.2 аппроксиматора модели определяется модель гармоники в форме дискретной передаточной функции (ДПФ) объекта-идентификатора. Далее в блоке 3 - восстановителя модели гармоники определяют параметры гармоники - круговую частоту и амплитуду. Полученную модель гармоники подают на вход блока 4 - восстановителя сигнала остатка, в котором из исходного сигнала вычитают модельные значения сигнала выявленной гармоники, и далее результаты сигнала остатка подаются на вход блока 1. Поскольку число гармоник исходного сигнала неизвестно, все действия повторяют до выполнения правила останова в блоке 2 - появления 4-ой ненулевой строки в матрице-идентификаторе. Выявление гармоник в блоке 2 осуществляется последовательно по мере убывания значений их амплитуд. Результат выполнения правила останова поступает на вход блока 5 - анализатора смещения, в котором определяют наличие смещения и выполняют восстановление его значения. Далее в блоке 6 - восстановителя модели сигнала определяют прогнозирующую модель мультисинусоидального сигнала. Затем прогнозирующая модель поступает на вход блока 7 - восстановителя модельных значений сигнала, в котором определяется модельный сигнал xм(kΔt).
Предлагаемый способ осуществляется следующим образом: с ДАС результаты измерений сигнала x(kΔt) (k=0, 1, 2, …) в равноотстоящие промежутки времени с шагом дискретизации Δt поступают на вход блока 1, в котором производится расчет первого структурного преобразования
второго структурного преобразования
и третьего структурного преобразования
где N1 - число измерений сигнала x(kΔt); N2, N3 - число расчетных значений второго и третьего структурных преобразований соответственно.
Результаты расчетов значений третьего структурного преобразования (3) передают на вход блока 2, где заносят их в 0-ую строку матрицы-идентификатора:
Последующие строки матрицы-идентификатора определяются с помощью формулы
Если анализируемый сигнал является мультисинусоидальным и представим в виде модели
где μ - постоянное смещение, n - число гармоник, ωi - круговая частота i-ой гармоники, Ci - амплитуда i-ой гармоники, a(kΔt) - шумовая составляющая, k=0, 1, 2, …, то в матрице (4) 4-ая строка будет являться нулевой. Запускается подсчет числа гармоник i=i+1 (первоначально i=0).
Первый столбец матрицы-идентификатора позволяет построить непрерывную цепную С-дробь в виде:
сворачивая которую, определяют модель гармоники в форме ДПФ объекта-идентификатора. Полученную модель передают на вход блока 3, где приступают к нахождению круговой частоты ω1 данной гармоники, имеющей максимальную амплитуду:
где z1=u±iv - полюс ДПФ (7). Тогда максимальную амплитуду гармоники находят по формуле:
где k=1, …, N3. Выявленная гармоника с максимальной амплитудой будет иметь вид:
Поскольку число гармоник исходного сигнала неизвестно, необходимо провести следующую процедуру. Полученную модель гармоники (10) передают на вход блока 4, где значения (10) вычитаются из исходного сигнала, получая на выходе блока сигнал остатка x1(kΔt)
После этого с сигналом (11) следует повторить всю процедуру для выявления остальных гармоник в нем. Результаты сигнала остатка (11) поступают на вход блока 1. Следующей гармоникой будет являться та, у которой значение амплитуды является наибольшим из оставшихся. Процедура повторяется до тех пор, пока в матрице-идентификаторе (4) элементы 4-й строки будут обращаться в нуль. Таким образом, критерием останова процедуры выявления гармонических компонент в исходном сигнале служит появление ненулевых элементов в 4-й строке матрицы-идентификатора.
После выполнения правила останова модель сигнала будет иметь вид
Она подается на вход блока 5, где определяют значение смещения
Если μ≠0, то фиксируют наличие смещения и принимают его равным (13). Результаты поступают на вход блока 6, где окончательно формируют прогнозирующую модель мультисинусоидального сигнала в виде (14)
если μ≠0, или в виде (12), если μ=0.
Идентифицированная модель мультисинусоидального сигнала поступает на вход блока 7, который восстанавливает значения модельного сигнала
Таким образом, предлагаемый способ идентификации мультисинусоидальных цифровых сигналов, отличается от известного тем, что восстановление гармоник в идентификаторе непрерывной дроби осуществляется последовательно по мере убывания значений их амплитуд до выполнения правила останова с предварительной фильтрационной обработкой значений сигнала, а после выполнения правила останова судят о наличии смещения и выполняют восстановление его значения.
Способ позволяет структурно преобразовать сигнал с целью фильтрации помех, автоматически оценить неизвестную заранее структуру математической модели мультисинусоидального сигнала, а также все его параметры. Вследствие использования метода процесс идентификации сигнала упрощается, увеличивается быстродействие, повышается точность, качество и достоверность результатов исследований. Универсальность и работоспособность метода продемонстрирована на примере различного рода мультисинусоидальных сигналов. В конечном итоге, предлагаемый способ позволяет достоверно прогнозировать значения анализируемых физических процессов различных типов, принимать адекватные решения по их контролю, управлению и диагностике.
Пример 1
График выходного сигнала спирального дозатора при подаче на вход йодида калия приведен на фиг. 2. К закономерным колебаниям выходного сигнала дозатора могут добавляться особые условия (зависание материала в бункере, очаговые уплотнения материала, влажность и т.п.), приводящие к стохастическим изменениям выходного сигнала дозатора. Пусть модель выходного сигнала имеет вид
С выхода ДАС поступают измерения анализируемого сигнала (шаг дискретизации Δt=0.5, число измерений N1=240), в блоке 1 над сигналом берутся три последовательных структурных преобразования. График третьего структурного преобразования приведен на фиг. 3.
Результаты расчетов значений третьего структурного преобразования (3) поступают на вход блока 2 и заполняют 0-ую строку матрицы-идентификатора:
Ее 4-ая строка обращается в нулевую. В подблоке 2.2 находят ДПФ объекта-идентификатора
Полученную модель передают на вход блока 3, где приступают к нахождению круговой частоты ω1 гармоники с максимальной амплитудой по формуле (8) ω1=2.00330. Сама амплитуда согласно (9) получается равна C1=0.49973.
На вход блока 4 подают полученную модель гармоники s1(kΔt)=0.49973sin(2.00330⋅kΔt), получая на выходе блока 4 сигнал остатка x1(kΔt)=x(kΔt)-0.49973sin(2.00330⋅kΔt). Измерения сигнала остатка поступают на вход блока 1, где с сигналом осуществляются три последовательных структурных преобразования. Результаты расчетов значений третьего структурного преобразования (3) поступают на вход блока 2 в 0-ую строку матрицы-идентификатора:
В 4-й строке матрицы-идентификатора появляются ненулевые элементы, что служит критерием останова процедуры выявления гармонических компонент в сигнале. Тогда модель сигнала будет иметь вид
Она подается на вход блока 5, где фиксируют наличие смещения и определяют его значение μ=8.02655 согласно (13). Результаты поступают на вход блока 6, где окончательно формируют прогнозирующую модель сигнала в виде
Максимальная абсолютная погрешность оценки параметров имеет порядок 10-2, структура сигнала (15) восстановлена верно. Идентифицированная модель сигнала поступает на вход блока 7, который восстанавливает значения модельного сигнала
Результаты расчетов модельных значений сигнала в блоке 7 приведены в Таблице 1, где - модельные значения сигнала, рассчитанные с помощью предлагаемого способа, - погрешности модельных значений. Сравнивая экспериментальные и модельные значения сигнала на выходе спирального дозатора и принимая во внимание наличие стохастических помех и погрешностей, можно сделать заключение о высокой точности идентификации выходного сигнала (с точностью восстановления до вычислительных погрешностей).
Пример 2.
Выходной вибросигнал роторного оборудования горнодобывающего производства имеет вид (график сигнала приведен на фиг. 4):
С выхода ДАС поступают измерения анализируемого сигнала (шаг дискретизации Δt=0.5, число измерений N1=100), в блоке 1 над сигналом берутся три последовательных структурных преобразования. Далее рассчитывается матрица-идентификатор, в которой 0-ая строка заполняется значениями третьего структурного преобразования (график третьего структурного преобразования приведен на фиг. 5), поступающих с выхода блока 1:
Ее 4-ая строка содержит нулевые значения. В подблоке 2.2 находят ДПФ объекта-идентификатора
Далее модель ДПФ объекта-идентификатора передают на вход блока 3, где приступают к нахождению параметров гармоники с максимальной амплитудой: круговая частота ω1=4.002, амплитуда C1=7.778.
На вход блока 4 подают полученную модель гармоники получая на выходе блока 4 сигнал остатка Сигнала остатка поступает на вход блока 1, где с сигналом осуществляются три последовательных структурных преобразования. Результаты расчетов значений третьего структурного преобразования (фиг. 6) поступают на вход блока 2 и заполняется матрица-идентификатор:
В подблоке 2.2 определяют ДПФ объекта-идентификатора
Далее модель ДПФ объекта-идентификатора передают на вход блока 3, где находят параметры второй гармоники: круговая частота ω2=1.411, амплитуда С2=1.888.
На вход блока 4 передают модель второй гармоники s2(kΔt)=1.888sin(1.411⋅kΔt), получая на выходе сигнал остатка x2(kΔt)=x1(kΔt)-1.888sin(1.411⋅kΔt), который предается на вход блока 1. Результаты расчетов третьего структурного преобразования (фиг. 7) поступают на вход блока 2 и заполняется матрица-идентификатор:
В подблоке 2.1 счетчик гармоник фиксирует наличие третьей гармоники и в подблоке 2.2 восстанавливают ее модель в форме ДПФ объекта-идентификатора
ДПФ объекта-идентификатора передают на вход блока 3, где находят параметры третьей гармоники: круговая частота ω3=2.017, амплитуда С3=1.026.
В блок 4 поступает модель третьей гармоники получая на выходе сигнал остатка Сигнал остатка поступает на вход блока 1, где производится расчет третьего структурного преобразования (фиг. 8). В блоке 2 заполняется матрица-идентификатор:
В 4-й строке матрицы-идентификатора появляются ненулевые элементы, что служит критерием останова процедуры выявления гармонических компонент в сигнале. Тогда модель сигнала будет иметь три гармоники:
Она подается на вход блока 5, где фиксируют наличие смещения и определяют его значение μ=10.047 согласно (13). Результаты поступают на вход блока 6, где окончательно формируют прогнозирующую модель мультисинусоидального сигнала
Максимальная абсолютная погрешность оценки параметров имеет порядок 10-2, структура модели мультисинусоидального сигнала (16) восстановлена верно. Идентифицированная модель сигнала поступает на вход блока 7, который восстанавливает значения модельного сигнала Результаты расчетов модельных значений сигнала в блоке 7 приведены в Таблице 2, где - модельные значения сигнала, рассчитанные с помощью предлагаемого способа, - погрешности модельных значений. Сравнивая экспериментальные и модельные значения вибросигнала на выходе роторного оборудования горнодобывающего производства, и принимая во внимание наличие стохастических помех и погрешностей, можно сделать заключение о высокой точности идентификации мультисиносидального выходного сигнала (с точностью восстановления до вычислительных погрешностей).
Claims (1)
- Способ идентификации мультисинусоидальных цифровых сигналов, включающий использование непрерывных цепных дробей путем измерения сигнала в равноотстоящие промежутки времени, подачи их на идентификатор непрерывной С-дроби и с последующим восстановлением числа гармоник, круговых частот, амплитуд, прогнозирующей модели и модельных значений сигнала, отличающийся тем, что выявление гармоник в идентификаторе непрерывной дроби осуществляется последовательно по мере убывания значений их амплитуд до выполнения правила останова с предварительной фильтрационной обработкой значений сигнала, а после выполнения правила останова судят о наличии смещения и выполняют восстановление его значения.
Publications (1)
Publication Number | Publication Date |
---|---|
RU2787309C1 true RU2787309C1 (ru) | 2023-01-09 |
Family
ID=
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2229725C1 (ru) * | 2002-11-12 | 2004-05-27 | Томский политехнический университет | Способ спектрального анализа периодического многочастотного сигнала |
RU2435168C1 (ru) * | 2010-04-09 | 2011-11-27 | ООО предприятие "КОНТАКТ-1" | Способ гармонического анализа периодического многочастотного сигнала |
US10080026B2 (en) * | 2015-03-20 | 2018-09-18 | Dolby Laboratories Licensing Corporation | Signal reshaping approximation |
RU2703933C1 (ru) * | 2018-11-08 | 2019-10-22 | федеральное государственное бюджетное образовательное учреждение высшего образования "Кемеровский государственный университет" (КемГУ) | Способ идентификации мультисинусоидальных цифровых сигналов |
RU2728832C2 (ru) * | 2013-02-05 | 2020-07-31 | Телефонактиеболагет Л М Эрикссон (Пабл) | Способ и устройство для управления маскировкой потери аудиокадров |
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2229725C1 (ru) * | 2002-11-12 | 2004-05-27 | Томский политехнический университет | Способ спектрального анализа периодического многочастотного сигнала |
RU2435168C1 (ru) * | 2010-04-09 | 2011-11-27 | ООО предприятие "КОНТАКТ-1" | Способ гармонического анализа периодического многочастотного сигнала |
RU2728832C2 (ru) * | 2013-02-05 | 2020-07-31 | Телефонактиеболагет Л М Эрикссон (Пабл) | Способ и устройство для управления маскировкой потери аудиокадров |
US10080026B2 (en) * | 2015-03-20 | 2018-09-18 | Dolby Laboratories Licensing Corporation | Signal reshaping approximation |
RU2703933C1 (ru) * | 2018-11-08 | 2019-10-22 | федеральное государственное бюджетное образовательное учреждение высшего образования "Кемеровский государственный университет" (КемГУ) | Способ идентификации мультисинусоидальных цифровых сигналов |
Non-Patent Citations (1)
Title |
---|
Статья: "Гибридный алгоритм идентификации частот мультисинусоидального сигнала", Ж. Научно-технический вестник информационных технологий, механики и оптики, номер 1 (83) 2013. * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8285498B2 (en) | Method for detection and automatic identification of damage to rolling bearings | |
Sedlak et al. | New automatic localization technique of acoustic emission signals in thin metal plates | |
RU2514061C2 (ru) | Способ определения событий вибрации с резонансной частотой в узле вращающихся лопаток | |
US7031873B2 (en) | Virtual RPM sensor | |
CN112326245B (zh) | 一种基于变分希尔伯特黄变换的滚动轴承故障诊断方法 | |
Minda et al. | A review of interpolation methods used for frequency estimation | |
RU2787309C1 (ru) | Способ идентификации мультисинусоидальных цифровых сигналов | |
JP6605855B2 (ja) | 絶縁異常診断装置 | |
CN104156339B (zh) | 一种利用二次排列熵识别周期微弱脉冲信号的方法 | |
KR102061616B1 (ko) | 비선형 스펙트럼 상관관계를 이용한 구조물 진단 방법 및 구조물 진단 시스템 | |
CN116298765A (zh) | 一种高温高湿反偏试验*** | |
RU2256950C2 (ru) | Способ идентификации линеаризованного динамического объекта | |
RU2703933C1 (ru) | Способ идентификации мультисинусоидальных цифровых сигналов | |
RU124411U1 (ru) | Система обработки и анализа данных цифровых регистраторов для мониторинга переходных режимов в энергообъединении | |
Henry et al. | Sensor validation via ultrasonic signal processing analysis | |
Hejtmanek et al. | Comparison of two denoising techniques to improve UHF partial discharge localization | |
Hou et al. | Adaptive maximum cyclostationarity blind deconvolution and its application on incipient fault diagnosis of bearing | |
Lesieutre et al. | A system response persistence measure for use in ambient data monitoring | |
RU2306592C1 (ru) | Способ активной идентификации линейных объектов управления | |
Herasimov et al. | Method for Assessing Meter Error Characteristics of Random Signals | |
RU2399060C1 (ru) | Способ анализа многочастотных сигналов, содержащих скрытые периодичности | |
CN117076932B (zh) | 高灵敏度电容变化检测方法、***、电子设备和存储介质 | |
Gillich et al. | ABOUT THE ACCURACY OF ESTIMATED FREQUENCIES WITH THE PyFEST SOFTWARE | |
CN109711036B (zh) | 飞行控制***试验结果的评估方法 | |
Totis et al. | Wavelet-like analysis in the frequency-damping domain for modal parameters identification |