RU2727295C1 - Waveform analysis method - Google Patents

Waveform analysis method Download PDF

Info

Publication number
RU2727295C1
RU2727295C1 RU2019109299A RU2019109299A RU2727295C1 RU 2727295 C1 RU2727295 C1 RU 2727295C1 RU 2019109299 A RU2019109299 A RU 2019109299A RU 2019109299 A RU2019109299 A RU 2019109299A RU 2727295 C1 RU2727295 C1 RU 2727295C1
Authority
RU
Russia
Prior art keywords
array
elements
sample
max
oscillogram
Prior art date
Application number
RU2019109299A
Other languages
Russian (ru)
Inventor
Дмитрий Александрович Недорезов
Original Assignee
Акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва" filed Critical Акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва"
Priority to RU2019109299A priority Critical patent/RU2727295C1/en
Application granted granted Critical
Publication of RU2727295C1 publication Critical patent/RU2727295C1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00Methods or arrangements for processing data by operating upon the order or content of the data handled
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06KGRAPHICAL DATA READING; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K11/00Methods or arrangements for graph-reading or for converting the pattern of mechanical parameters, e.g. force or presence, into electrical signal

Abstract

FIELD: computing.SUBSTANCE: invention relates to calculation methods used for training of computer systems, and can be used for analysis of any waveforms, any types of signals, which have periodically repeating part, as well as one, and only one, upward transition from conventionally designated area of low amplitudes to conventionally designated region of high amplitudes within each individual period method allows detecting not only limited range of signal anomalies, such as wounds, glitches, pulses of certain duration, defined intervals between pulses, fronts, decays and the like, but any periods different from the automatically calculated sample, since the method is based on point-by-point comparison of a sample waveform with an oscillogram of the analyzed signal.EFFECT: technical result is simplification and acceleration of intelligent computer system training process.1 cl, 3 dwg

Description

Изобретение относится к способам вычислений, используемых для обучения компьютерных систем.The invention relates to computational methods used to train computer systems.

Известен способ поиска аномалий в осциллограмме WaveScan от компании Teledyne Lecroy (WaveScan [Электронный ресурс] // Сайт компании Teledyne Lecroy. 2018. URL: http://teledynelecroy.com/features/featureoverview.aspx?modelid=2107&capid=102&mid=556). Способ состоит в том, что воспроизводят осциллограмму целиком, либо частично на дисплее компьютерной системы в виде статического изображения. Затем проводят обучение компьютерной системы путем указания на часть осциллограммы, которую необходимо принять за графический образец для анализа на предмет наличия в осциллограмме таких же частей как образец. Либо проводят обучение компьютерной системы в ручном режиме, путем задания множества числовых параметров, описывающих искомые аномалии сигнала. При обнаружении данных аномалий компьютерная система сообщает об их наличии и сохраняет их изображения.There is a known method of searching for anomalies in the WaveScan waveform from Teledyne Lecroy (WaveScan [Electronic resource] // Teledyne Lecroy website. 2018. URL: http://teledynelecroy.com/features/featureoverview.aspx?modelid=2107&capid=102&mid=556) ... The method consists in the fact that the oscillogram is reproduced in whole or in part on the display of a computer system in the form of a static image. Then, the computer system is trained by pointing to a part of the waveform that must be taken as a graphic sample for analysis for the presence of the same parts as the sample in the waveform. Or, a computer system is trained in manual mode by setting a set of numerical parameters describing the desired signal anomalies. When these anomalies are detected, the computer system reports their presence and stores their images.

Недостатками данного способа являются: невозможность обнаружения аномалий сигнала не подпадающих под номенклатуру аномалий, которые способен обнаруживать данный способ, т.е анализ осциллограммы проводится только на предмет наличия определенных фронтов, нарушений монотонности, рантов, заданных измерений; необходимость ввода входных параметров поиска аномалий электрического сигнала, что замедляет процесс обучения компьютерной системы.The disadvantages of this method are: the impossibility of detecting signal anomalies that do not fall under the nomenclature of anomalies that this method is capable of detecting, that is, the analysis of the oscillogram is carried out only for the presence of certain fronts, violations of monotony, welts, specified measurements; the need to enter the input parameters for searching for electrical signal anomalies, which slows down the learning process of the computer system.

В рамках данной заявки, «Аномалия сигнала» - это период электрического сигнала или его часть, отличающиеся от периодов сигнала, принятых в качестве нормальных.Within the framework of this application, "Signal anomaly" is a period of an electrical signal, or a part of it, which differs from the periods of the signal taken as normal.

Наиболее близким (прототипом) является «Способ интеллектуальной синхронизации TriggerScan» от компании Teledyne Lecroy (TriggerScan [Электронный ресурс] // Сайт компании Teledyne Lecroy. 2018. URL: http://teledynelecroy.com/features/featureoverview.aspx?modelid=2108&capid=102&mid=556). Способ состоит в том, что воспроизводят осциллограмму целиком, либо частично на дисплее компьютерной системы в виде статического изображения. Затем проводят обучение компьютерной системы путем указания на часть осциллограммы, которую необходимо принять за графический образец. Далее компьютерная система на основе заданного графического образца автоматически рассчитывает десятки числовых параметров описывающих сигнал, принятый за образец и, на основе данных параметров, может обнаруживать некоторые части осциллограммы, отличные от заданного графического образца, а также сообщать о них специалисту, проводящему анализ, и сохранять их изображения.The closest (prototype) is the "TriggerScan intelligent synchronization method" from Teledyne Lecroy (TriggerScan [Electronic resource] // Teledyne Lecroy website. 2018. URL: http://teledynelecroy.com/features/featureoverview.aspx?modelid=2108&capid = 102 * mid = 556). The method consists in the fact that the oscillogram is reproduced in whole or in part on the display of a computer system in the form of a static image. Then, the computer system is trained by pointing to the part of the oscillogram that must be taken as a graphic sample. Further, the computer system, based on a given graphic sample, automatically calculates dozens of numerical parameters describing the signal taken as a sample and, based on these parameters, can detect some parts of the oscillogram that are different from the given graphic sample, and also inform the specialist conducting the analysis about them and save their images.

Недостатком прототипа является необходимость ввода входных параметров поиска аномалий электрического сигнала, что усложняет и замедляет процесс обучения компьютерной системы.The disadvantage of the prototype is the need to enter the input parameters for searching for electrical signal anomalies, which complicates and slows down the learning process of the computer system.

Для заявленного способа выявлены общие с прототипом существенные признаки: проводят обучение компьютерной системы, в процессе которого рассчитывают десятки числовых параметров описывающих образец сигнала и на основе данных параметров обнаруживают периоды осциллограммы отличные от рассчитанного образца, сообщают о них и сохраняют их изображения.For the claimed method, essential features common to the prototype have been revealed: the computer system is trained, during which tens of numerical parameters describing the signal sample are calculated and, based on these parameters, the oscillogram periods differ from the calculated sample are detected, reported and their images are saved.

Технической проблемой заявленного изобретения является упрощение и ускорение процесса обучения компьютерной системы.The technical problem of the claimed invention is to simplify and accelerate the learning process of a computer system.

Техническая проблема изобретения решается путем применения алгоритма расчета входных параметров поиска аномалий электрического сигнала, который состоит в следующем: подсчитывают общее количество точек S в осциллограмме, вычисляют коэффициент десятой части осциллограммы Sdc путем умножения S на 0.1 Sdc=S×0.1, округляют Sdc до целого значения, представляют начальный участок осциллограммы в виде массива числовых значений точек осциллограммы от первого ее элемента до элемента с индексом Sdc, присваивают данному массиву обозначение Ms, делят получившийся массив на 10 равных или приближенно равных частей, если в массиве имеется 10 или более элементов, если элементов в массиве менее 10, то делят Ms на количество частей равное количеству элементов Ms. Вычисляют минимальную и максимальную величину элемента по оси ординат в каждой части, затем вычисляют среднее арифметическое Ymin для минимальных значений всех частей массива Ms и среднее арифметическое Ymax для максимальных значений. Затем определяют максимальную амплитуду сигнала: Amax=Ymax-Ymin. Затем задают коэффициент средних игнорируемых амплитуд K, который может принимать значения от 0 до 0,5, либо оставляют значение по умолчанию 0,4, это делается для исключения из дальнейших расчетов, выполняемых для обнаружения переходов из условно обозначенной области (УОО) низких амплитуд в УОО высоких амплитуд, часто встречающихся паразитных колебаний в районе нулевой амплитуды при анализе переменного тока или средней амплитуды при анализе постоянного тока, величина K зависит от амплитуды помех. Далее определяют верхнюю границу диапазона УОО низких амплитуд: Nmax=Ymin+Amax×K и определяют нижнюю границу диапазона УОО низких амплитуд: Nmin=-∞. Определяют нижнюю границу диапазона УОО высоких амплитуд: Vmin=Ymax-Amax×K и определяют верхнюю границу диапазона УОО высоких амплитуд: Vmax=∞ (фиг. 1). Записывают все значения точек осциллограммы по оси ординат в виде чисел в массив MS0, состоящий из S элементов, причем сохраняют порядок следования элементов в MS0 таким же, как в осциллограмме-образце. Далее, для каждого элемента MS0 последовательно рассчитывают среднее значение из самого элемента и следующего сразу за ним, получившиеся средние значения записывают в массив MS1 в том же порядке, что в MS0, получают результат в виде массива такой же размерности - это делается для игнорирования в расчетах краткосрочных низкоамплитудных помех, которые могут усложнить обнаружение переходов из УОО низких амплитуд в УОО высоких. Для каждого элемента MS1 последовательно вычисляют к какому из диапазонов УОО он относится: низких амплитуд, высоких амплитуд или не относится ни к тому, ни к другому. По результатам вычислений составляют массив MS2, в котором заменяют каждый числовой элемент элементом логического типа, где все элементы из диапазона низких амплитуд заменяют константой «Ложь», все элементы из диапазона высоких амплитуд заменяют константой «Истина» значения, не вошедшие ни в тот, ни в другой диапазон, игнорируют и не переносят в массив MS2 - это делается для исключения из расчетов, выполняемых для обнаружения переходов из УОО низких амплитуд в УОО высоких амплитуд, часто встречающихся паразитных колебаний в районе нулевой амплитуды при анализе переменного тока или средней амплитуды при анализе постоянного тока. Вычисляют порядковые номера элементов массива MS2, которые имеют значение «Истина» и при этом предыдущий перед ними элемент имеет значение «Ложь», то есть обнаруживают места перехода от УОО низких к УОО высоких амплитуд. Выделяют область числового массива MS1 между двумя первыми переходами, получают точный образец осциллограммы одного периода, который считают окончательно рассчитанным образцом для последующих вычислений, вычисляют размерность этого массива. Таким же образом, как был вычислен период сигнала в образце, вычисляют все периоды во всей анализируемой осциллограмме, тем самым, получают числовой массив периода образца и множество числовых массивов всех периодов осциллограммы одинаковой размерности. Поочередно сравнивают все массивы периодов с образцом, причем первый элемент каждого числового массива сравнивают с первым элементом образца, второй со вторым, третий с третьим и так далее, поочередно, до последних элементов массива - данный шаг позволяет провести абсолютно полное сравнение периода образца с анализируемыми периодами, каждый минимально выделяемый элемент сравнивается с каждым («Поточечное сравнение»), что позволяет обнаружить абсолютно любое минимальное отличие от образца (аномалию). Затем вычисляют значения разностей элементов («Поточечного сравнения»), то есть отклонения от образца по оси ординат. Далее подсчитывают количество отклонений для каждого анализируемого периода осциллограммы по пяти критериям, значение разности по оси ординат входит в диапазон: от 5 до 10, от 10 до 20, от 20 до 30 и от 60 до ∞ процентов от Amax. Специалист, проводящий анализ, вводит критерии количества выходов элементов массивов периодов осциллограммы за пределы каждого из пяти диапазонов по оси ординат, в соответствие с которыми компьютерная система должна принять решение сообщать и формировать изображения анализируемого периода осциллограммы или не делать этого, то есть проигнорировать период, посчитав, что он недостаточно отличается от рассчитанного образца, причем компьютерная система сообщает о периоде при выходе, хотя бы из одного диапазона и обозначает какие именно диапазоны были превышены. Также имеется возможность оставить критерии по умолчанию и тем самым исключить работу специалиста, возложив ее, например, на ЭВМ.The technical problem of the invention is solved by applying the algorithm for calculating the input parameters for searching for electrical signal anomalies, which consists in the following: the total number of points S in the oscillogram is calculated, the coefficient of the tenth part of the oscillogram S dc is calculated by multiplying S by 0.1 S dc = S × 0.1, S dc is rounded to an integer value, represent the initial portion of the waveform as an array of numerical values of the waveform points from its first element to the element with index S dc , assign the designation M s to this array, divide the resulting array into 10 equal or approximately equal parts, if the array contains 10 or more elements, if there are less than 10 elements in the array, then divide M s by the number of parts equal to the number of elements M s . Calculate the minimum and maximum value of the element along the ordinate in each part, then calculate the arithmetic mean Y min for the minimum values of all parts of the array M s and the arithmetic mean Y max for the maximum values. Then, the maximum signal amplitude is determined: A max = Y max -Y min . Then the coefficient of the average ignored amplitudes K is set, which can take values from 0 to 0.5, or leave the default value of 0.4, this is done to exclude from further calculations performed to detect transitions from the conventionally designated area (ROO) of low amplitudes into RVO of high amplitudes, frequent parasitic oscillations in the region of zero amplitude when analyzing alternating current or average amplitude when analyzing direct current, the value of K depends on the amplitude of the noise. Next, determine the upper limit of the range of low amplitudes VAR: N max = Y min + A max × K and determine the lower limit of the range of low amplitudes VAR: N min = -∞. Determine the lower limit of the range of high amplitudes VAR: V min = Y max -A max × K and determine the upper limit of the range of high amplitudes of the VAR: V max = ∞ (Fig. 1). All the values of the points of the oscillogram along the ordinate axis are written in the form of numbers into the MS 0 array, consisting of S elements, and the order of the elements in MS 0 is kept the same as in the sample oscillogram. Further, for each element MS 0 , the average value is sequentially calculated from the element itself and the next one immediately after it, the resulting average values are written into the MS 1 array in the same order as in MS 0 , the result is obtained as an array of the same dimension - this is done for ignoring short-term low-amplitude noise in the calculations, which can complicate the detection of transitions from the low-amplitude ERR to the high-ERO. For each element of MS 1, it is sequentially calculated to which of the RVO ranges it belongs to: low amplitudes, high amplitudes, or does not apply to either. Based on the results of the calculations, an array MS 2 is compiled, in which each numeric element is replaced with an element of a logical type, where all elements from the range of low amplitudes are replaced by the constant "False", all elements from the range of high amplitudes are replaced by the constant "True" values that are not included in either to no other range, they are ignored and not transferred to the MS 2 array - this is done to exclude from the calculations performed to detect transitions from the low-amplitude CCD to the high-amplitude CCD, often occurring parasitic oscillations in the region of zero amplitude when analyzing the alternating current or the average amplitude at DC analysis. The ordinal numbers of the elements of the MS 2 array are calculated, which have the value "True" and the preceding element before them has the value "False", that is, the places of transition from the low to high amplitudes of the low to high amplitudes are detected. Allocate the area of the numerical array MS 1 between the first two transitions, obtain an exact sample of the oscillogram of one period, which is considered the final calculated sample for subsequent calculations, calculate the dimension of this array. In the same way as the period of the signal in the sample was calculated, all periods in the entire analyzed oscillogram are calculated, thereby obtaining a numeric array of the sample period and a set of numeric arrays of all periods of the oscillogram of the same dimension. All arrays of periods are compared in turn with the sample, and the first element of each numeric array is compared with the first element of the sample, the second with the second, the third with the third, and so on, in turn, until the last elements of the array - this step allows you to carry out an absolutely complete comparison of the period of the sample with the analyzed periods , each minimum selected element is compared with each ("Pointwise comparison"), which allows you to detect absolutely any minimum difference from the sample (anomaly). Then the values of the differences of the elements ("Pointwise comparison"), that is, the deviation from the sample along the ordinate axis, are calculated. Next, the number of deviations for each analyzed period of the oscillogram is calculated according to five criteria, the value of the difference along the ordinate is in the range: from 5 to 10, from 10 to 20, from 20 to 30 and from 60 to ∞ percent of A max . The analyst introduces the criteria for the number of outputs of the elements of the waveform period arrays beyond each of the five ranges along the ordinate axis, in accordance with which the computer system must decide to report and form images of the analyzed waveform period or not, that is, ignore the period, counting that it does not differ enough from the calculated sample, and the computer system informs about the period upon exiting, at least from one range, and indicates which ranges have been exceeded. It is also possible to leave the default criteria and thereby exclude the work of a specialist, assigning it, for example, to a computer.

Применение заявленного способа анализа осциллограмм позволяет обойтись без операции ввода входных параметров поиска аномалий электрического сигнала, что упрощает и ускоряет процесс.The application of the claimed method for analyzing oscillograms makes it possible to dispense with the operation of entering the input parameters for searching for electrical signal anomalies, which simplifies and speeds up the process.

В рамках данной заявки, «Период осциллограммы» - графическое отображение периода сигнала на осциллограмме.Within the framework of this application, "Waveform Period" is a graphical display of the signal period on the waveform.

Заявка поясняется изображениями:The application is illustrated by images:

Фиг. 1 - изображение осциллограммы синусоидальной формы, используемой для осуществления предлагаемого способа с обозначениями: УОО низких, средних и высоких амплитуд, единственного восходящего перехода от УОО низких к УОО высоких амплитуд и периода сигнала;FIG. 1 is an image of a sinusoidal oscillogram used for the implementation of the proposed method with the designations: low, medium and high amplitudes, the only upward transition from low to high low to high amplitudes and period of the signal;

Фиг. 2 - изображение осциллограммы треугольной формы, используемой для осуществления предлагаемого способа с обозначениями: УОО низких, средних и высоких амплитуд, единственного восходящего перехода от УОО низких к УОО высоких амплитуд и периода сигнала;FIG. 2 is an image of a triangular oscillogram used to implement the proposed method with the designations: low, medium and high amplitudes, the only upward transition from low to high low to high amplitudes and signal period;

Фиг. 3 - изображение осциллограммы прямоугольной формы, используемой для осуществления предлагаемого способа с обозначениями: УОО низких, средних и высоких амплитуд, единственного восходящего перехода от УОО низких к УОО высоких амплитуд и периода сигнала.FIG. 3 is an image of a rectangular oscillogram used for the implementation of the proposed method with the designations: low, medium and high amplitudes RVO, the only upward transition from the low RVR to the high RVR high amplitudes and signal period.

Способ применим для анализа любых типов сигналов, которые имеют в своем составе периодически повторяющуюся часть, а также один, и только один, восходящий переход от условно обозначенной области низких амплитуд к условно обозначенной области высоких амплитуд в пределах каждого отдельного периода. К такому типу относится большинство сигналов встречающихся в природе, в том числе сигналы синусоидальной (фиг. 1), треугольной (фиг. 2) и прямоугольной формы (фиг. 3).The method is applicable for analyzing any types of signals that include a periodically repeating part, as well as one, and only one, ascending transition from a conventionally designated area of low amplitudes to a conventionally designated area of high amplitudes within each separate period. This type includes most of the signals found in nature, including sinusoidal (Fig. 1), triangular (Fig. 2) and rectangular (Fig. 3) signals.

Способ осуществляют следующим образом.The method is carried out as follows.

Разрабатывают программу для ЭВМ в которой предусматривают функцию вывода данных из файлов, содержащих осциллограмму электрического сигнала. Подсчитывают общее количество точек в осциллограмме - S, вычисляют коэффициент десятой части осциллограммы Sdc путем умножения S на 0.1 Sdc=S×0.1, округляют Sdc до целого значения, представляют начальный участок осциллограммы в виде массива числовых значений точек осциллограммы от первого ее элемента до элемента с индексом Sdc, присваивают данному массиву обозначение - Ms, делят получившийся массив на 10 равных или приближенно равных частей, если в массиве имеется 10 или более элементов, если элементов в массиве менее 10, то делят Ms на количество частей равное количеству элементов Ms. Например, в массиве Ms 133333 элемента, тогда его возможно разделить на 10 частей Ms0…9 состоящих из 13333 элементов каждая, но 3 элемента останутся в остатке, их следует добавить в любые 3 части из полученных, т.к. больше одного остаточного элемента в одну часть добавлять нельзя. Например, их можно добавить в части Ms0, Ms1, Ms2 или в Ms3, Ms8, Ms9 и так далее. Или в массиве Ms 7 элементов, тогда его возможно разделить на 7 равных частей по 1 элементу в каждой.A computer program is being developed in which the function of outputting data from files containing an oscillogram of an electrical signal is provided. Count the total number of points in the oscillogram - S, calculate the coefficient of the tenth part of the oscillogram S dc by multiplying S by 0.1 S dc = S × 0.1, round S dc to an integer value, represent the initial section of the oscillogram as an array of numerical values of the oscillogram points from its first element to the element with index S dc , assign the designation to this array - M s , divide the resulting array into 10 equal or approximately equal parts, if the array has 10 or more elements, if the elements in the array are less than 10, then divide M s by the number of parts equal to the number of elements M s . For example, in the array M s 133333 elements, then it is possible to divide it into 10 parts M s0 ... 9 consisting of 13333 elements each, but 3 elements will remain in the remainder, they should be added to any 3 parts from the received ones, since more than one residual element cannot be added to one part. For example, they can be added to parts M s0 , M s1 , M s2 or to M s3 , M s8 , M s9, and so on. Or M s in the array elements 7, then it may be divided into 7 equal portions on each element 1.

Вычисляют минимальную и максимальную величину элемента по оси ординат в каждой части, затем вычисляют среднее арифметическое Ymin для минимальных значений всех частей массива Ms и среднее арифметическое Ymax для максимальных значений. Затем определяют максимальную амплитуду сигнала: Amax=Ymax-Ymin. Затем, специалист проводящий анализ, задает коэффициент средних игнорируемых амплитуд K, который может принимать значения от 0 до 0,5, либо оставляет значение по умолчанию 0,4, это делается для исключения из дальнейших расчетов, выполняемых для обнаружения переходов из условно обозначенной области (УОО) низких амплитуд в УОО высоких амплитуд, часто встречающихся паразитных колебаний в районе нулевой амплитуды при анализе переменного тока или средней амплитуды при анализе постоянного тока, величина K зависит от амплитуды помех. Далее определяют верхнюю границу диапазона УОО низких амплитуд: Nmax=Ymin+Amax×K и определяют нижнюю границу диапазона УОО низких амплитуд: Nmin=-∞. Определяют нижнюю границу диапазона УОО высоких амплитуд: Vmin=Ymax-Amax×K и определяют верхнюю границу диапазона УОО высоких амплитуд: Vmax=∞ (фиг. 1). Записывают все значения точек осциллограммы по оси ординат в виде чисел в массив MS0, состоящий из S элементов, причем сохраняют порядок следования элементов в MS0 таким же, как в осциллограмме-образце. Далее, для каждого элемента MS0 последовательно рассчитывают среднее значение из самого элемента и следующего сразу за ним, получившиеся средние значения записывают в массив MS1 в том же порядке, что в MS0, получают результат в виде массива такой же размерности - это делается для игнорирования в расчетах краткосрочных низкоамплитудных помех, которые могут усложнить обнаружение переходов из УОО низких амплитуд в УОО высоких. Для каждого элемента MS1 последовательно вычисляют к какому из диапазонов УОО он относится: низких амплитуд, высоких амплитуд или не относится ни к тому, ни к другому. По результатам вычислений составляют массив MS2, в котором заменяют каждый числовой элемент элементом логического типа, где все элементы из диапазона низких амплитуд заменяют константой «Ложь», все элементы из диапазона высоких амплитуд заменяют константой «Истина» значения, не вошедшие ни в тот, ни в другой диапазон, игнорируют и не переносят в массив MS2 - это делается для исключения из расчетов, выполняемых для обнаружения переходов из УОО низких амплитуд в УОО высоких амплитуд, часто встречающихся паразитных колебаний в районе нулевой амплитуды при анализе переменного тока или средней амплитуды при анализе постоянного тока. Вычисляют порядковые номера элементов массива MS2, которые имеют значение «Истина» и при этом предыдущий перед ними элемент имеет значение «Ложь», то есть обнаруживают места перехода от УОО низких к УОО высоких амплитуд. Выделяют область числового массива MS1 между двумя первыми переходами, получают точный образец осциллограммы одного периода, который считают окончательно рассчитанным образцом для последующих вычислений, вычисляют размерность этого массива. Таким же образом, как был вычислен период сигнала в образце, вычисляют все периоды во всей анализируемой осциллограмме, тем самым, получают числовой массив периода образца и множество числовых массивов всех периодов осциллограммы одинаковой размерности. Поочередно сравнивают все массивы периодов с образцом, причем первый элемент каждого числового массива сравнивают с первым элементом образца, второй со вторым, третий с третьим и так далее, поочередно, до последних элементов массива - данный шаг позволяет провести абсолютно полное сравнение периода образца с анализируемыми периодами, каждый минимально выделяемый элемент сравнивается с каждым («Поточечное сравнение»), что позволяет обнаружить абсолютно любое минимальное отличие от образца (аномалию). Затем вычисляют значения разностей элементов («Поточечного сравнения»), то есть отклонения от образца по оси ординат. Далее подсчитывают количество отклонений для каждого анализируемого периода осциллограммы по пяти критериям, значение разности по оси ординат входит в диапазон: от 5 до 10, от 10 до 20, от 20 до 30 и от 60 до ∞ процентов от Amax. Специалист, проводящий анализ, вводит критерии количества выходов элементов массивов периодов осциллограммы за пределы каждого из пяти диапазонов по оси ординат, в соответствие с которыми компьютерная система должна принять решение сообщать и формировать изображения анализируемого периода осциллограммы или не делать этого, то есть проигнорировать период, посчитав, что он недостаточно отличается от рассчитанного образца, причем компьютерная система сообщает о периоде при выходе, хотя бы из одного диапазона и обозначает какие именно диапазоны были превышены. Также имеется возможность оставить критерии по умолчанию и тем самым исключить работу специалиста, возложив ее, например, на ЭВМ.Calculate the minimum and maximum value of the element along the ordinate in each part, then calculate the arithmetic mean Y min for the minimum values of all parts of the array M s and the arithmetic mean Y max for the maximum values. Then, the maximum signal amplitude is determined: A max = Y max -Y min . Then, the specialist conducting the analysis sets the coefficient of the average ignored amplitudes K, which can take values from 0 to 0.5, or leaves the default value of 0.4, this is done to exclude from further calculations performed to detect transitions from the conventionally designated area ( UOO) low amplitudes in the UOO of high amplitudes, frequent parasitic oscillations in the region of zero amplitude when analyzing alternating current or average amplitude when analyzing direct current, the value of K depends on the amplitude of the noise. Next, determine the upper limit of the range of low amplitudes VAR: N max = Y min + A max × K and determine the lower limit of the range of low amplitudes VAR: N min = -∞. Determine the lower limit of the range of high amplitudes VAR: V min = Y max -A max × K and determine the upper limit of the range of high amplitudes of the VAR: V max = ∞ (Fig. 1). All the values of the points of the oscillogram along the ordinate axis are written in the form of numbers into the MS 0 array, consisting of S elements, and the order of the elements in MS 0 is kept the same as in the sample oscillogram. Further, for each element MS 0 , the average value is sequentially calculated from the element itself and the next one immediately after it, the resulting average values are written into the MS 1 array in the same order as in MS 0 , the result is obtained as an array of the same dimension - this is done for ignoring short-term low-amplitude noise in the calculations, which can complicate the detection of transitions from the low-amplitude ERR to the high-ERO. For each element of MS 1, it is sequentially calculated to which of the RVO ranges it belongs to: low amplitudes, high amplitudes, or does not apply to either. Based on the results of the calculations, an array MS 2 is compiled, in which each numeric element is replaced with an element of a logical type, where all elements from the range of low amplitudes are replaced by the constant "False", all elements from the range of high amplitudes are replaced by the constant "True" values that are not included in either to no other range, they are ignored and not transferred to the MS 2 array - this is done to exclude from the calculations performed to detect transitions from the low-amplitude CCD to the high-amplitude CCD, often occurring parasitic oscillations in the region of zero amplitude when analyzing the alternating current or the average amplitude at DC analysis. The ordinal numbers of the elements of the MS 2 array are calculated, which have the value "True" and the preceding element before them has the value "False", that is, the places of transition from the low to high amplitudes of the low to high amplitudes are detected. Allocate the area of the numerical array MS 1 between the first two transitions, obtain an exact sample of the oscillogram of one period, which is considered the final calculated sample for subsequent calculations, calculate the dimension of this array. In the same way as the period of the signal in the sample was calculated, all periods in the entire analyzed oscillogram are calculated, thereby obtaining a numeric array of the sample period and a set of numeric arrays of all periods of the oscillogram of the same dimension. All arrays of periods are compared in turn with the sample, and the first element of each numeric array is compared with the first element of the sample, the second with the second, the third with the third, and so on, in turn, until the last elements of the array - this step allows you to carry out an absolutely complete comparison of the period of the sample with the analyzed periods , each minimum selected element is compared with each ("Pointwise comparison"), which allows you to detect absolutely any minimum difference from the sample (anomaly). Then the values of the differences of the elements ("Pointwise comparison"), that is, the deviation from the sample along the ordinate axis, are calculated. Next, the number of deviations for each analyzed period of the oscillogram is calculated according to five criteria, the value of the difference along the ordinate is in the range: from 5 to 10, from 10 to 20, from 20 to 30 and from 60 to ∞ percent of A max . The analyst introduces the criteria for the number of outputs of the elements of the waveform period arrays beyond each of the five ranges along the ordinate axis, in accordance with which the computer system must decide to report and form images of the analyzed waveform period or not, that is, ignore the period, counting that it does not differ enough from the calculated sample, and the computer system informs about the period upon exiting, at least from one range, and indicates which ranges have been exceeded. It is also possible to leave the default criteria and thereby exclude the work of a specialist, assigning it, for example, to a computer.

Claims (1)

Способ анализа осциллограмм, при котором вычисляют минимальное значение амплитуды образца сигнала Ymin и максимальное значение амплитуды образца сигнала Ymax; определяют максимальную амплитуду сигнала: Amax=Ymax - Ymin; подсчитывают общее количество точек S в осциллограмме; затем задают коэффициент средних игнорируемых амплитуд K, который может принимать значения от 0 до 0,5, либо оставляют значение по умолчанию 0,4; далее определяют верхнюю границу диапазона условно обозначенной области (УОО) низких амплитуд: Nmax = Ymin + Amax × K и определяют нижнюю границу диапазона УОО низких амплитуд: Nmin = -∞; определяют нижнюю границу диапазона УОО высоких амплитуд: Vmin = Ymax - Amax × K и определяют верхнюю границу диапазона УОО высоких амплитуд: Vmax = ∞; записывают все значения точек осциллограммы по оси ординат в виде чисел в массив MS0, состоящий из S элементов, причем сохраняют порядок следования элементов в MS0 таким же, как в осциллограмме-образце; для каждого элемента MS0 последовательно рассчитывают среднее значение из самого элемента и следующего сразу за ним, получившиеся средние значения записывают в массив MS1 в том же порядке, что в MS0, получают результат в виде массива такой же размерности; для каждого элемента MS1 последовательно вычисляют к какому из диапазонов УОО он относится: низких амплитуд, высоких амплитуд или не относится ни к тому, ни к другому; по результатам вычислений составляют массив MS2, в котором заменяют каждый числовой элемент элементом логического типа, где все элементы из диапазона низких амплитуд заменяют константой «Ложь», все элементы из диапазона высоких амплитуд заменяют константой «Истина», значения, не вошедшие ни в тот, ни в другой диапазон, игнорируют и не переносят в массив MS2; вычисляют порядковые номера элементов массива MS2, которые имеют значение «Истина» и при этом предыдущий перед ними элемент имеет значение «Ложь», то есть обнаруживают места перехода от УОО низких к УОО высоких амплитуд; выделяют область числового массива MS1 между двумя первыми переходами, получают точный образец осциллограммы одного периода, который считают окончательно рассчитанным образцом для последующих вычислений, вычисляют размерность этого массива; таким же образом, как был вычислен период сигнала в образце, вычисляют все периоды во всей анализируемой осциллограмме, тем самым, получают числовой массив периода образца и множество числовых массивов всех периодов осциллограммы одинаковой размерности; поочередно сравнивают все массивы периодов с образцом, причем первый элемент каждого числового массива сравнивают с первым элементом образца, второй со вторым, третий с третьим и так далее, поочередно, до последних элементов массива; затем вычисляют значения разностей элементов, то есть отклонения от образца по оси ординат; далее подсчитывают количество отклонений для каждого анализируемого периода осциллограммы по пяти критериям, значение разности по оси ординат входит в диапазон: от 5 до 10, от 10 до 20, от 20 до 30 и от 60 до ∞ процентов от Amax. вводят критерии количества выходов элементов массивов периодов осциллограммы за пределы каждого из пяти диапазонов по оси ординат, в соответствии с которыми компьютерная система должна принять решение сообщать и формировать изображения анализируемого периода осциллограммы или не делать этого, то есть проигнорировать период, посчитав, что он недостаточно отличается от заданного образца, причем компьютерная система сообщает о периоде при выходе хотя бы из одного диапазона и обозначает, какие именно диапазоны были превышены; или оставляют критерии по умолчанию, отличающийся тем, что входные параметры Ymax и Ymin рассчитывают следующим образом: вычисляют коэффициент десятой части осциллограммы Sdc путем умножения S на 0.1 Sdc = S × 0.1, округляют Sdc до целого значения, представляют начальный участок осциллограммы в виде массива числовых значений точек осциллограммы от первого ее элемента до элемента с индексом Sdc, присваивают данному массиву обозначение Ms, делят получившийся массив на 10 равных частей, если в массиве имеется 10 или более элементов, если элементов в массиве менее 10, то делят Ms на количество частей равное количеству элементов Ms, если деление на 10 равных частей невозможно из-за неделимости числа элементов Ms без остатка, то остаточные элементы массива распределяют между произвольными частями, но не более одного элемента на часть; вычисляют минимальную и максимальную величину элемента по оси ординат в каждой части, затем вычисляют среднее арифметическое Ymin для минимальных значений всех частей массива Ms и вычисляют среднее арифметическое Ymax для максимальных значений всех частей массива Ms.A method for analyzing oscillograms, which calculates the minimum value of the amplitude of the signal sample Y min and the maximum value of the amplitude of the signal sample Y max ; determine the maximum signal amplitude: A max = Y max - Y min ; count the total number of points S in the oscillogram; then the coefficient of the average ignored amplitudes K is set, which can take values from 0 to 0.5, or leave the default value of 0.4; then determine the upper limit of the range of the conventionally designated area (ROO) of low amplitudes: N max = Y min + A max × K and determine the lower limit of the range of the RR of low amplitudes: N min = -∞; determine the lower limit of the range of high amplitudes VOO: V min = Y max - A max × K and determine the upper limit of the range of high amplitudes VOO: V max = ∞; write all the values of the points of the oscillogram along the ordinate axis in the form of numbers in the array MS 0 , consisting of S elements, and keep the order of the elements in MS 0 the same as in the sample oscillogram; for each element MS 0 , the average value is sequentially calculated from the element itself and the one immediately following it, the resulting average values are written into the MS 1 array in the same order as in MS 0 , the result is obtained in the form of an array of the same dimension; for each element of MS 1, it is sequentially calculated to which of the ROV ranges it belongs to: low amplitudes, high amplitudes, or does not belong to either one or the other; according to the results of calculations, an array MS 2 is compiled, in which each numeric element is replaced with an element of a logical type, where all elements from the range of low amplitudes are replaced by the constant "False", all elements from the range of high amplitudes are replaced by the constant "True", values that are not included in that , nor in any other range, are ignored and do not transfer to the MS array 2 ; calculate the ordinal numbers of the elements of the array MS 2 , which have the value "True" and the preceding element before them has the value "False", that is, they find the places of transition from the low to high amplitudes of the low to high amplitudes; allocate the area of the numerical array MS 1 between the first two transitions, obtain an accurate sample of the oscillogram of one period, which is considered the finally calculated sample for subsequent calculations, calculate the dimension of this array; in the same way as the period of the signal in the sample was calculated, all periods in the entire analyzed oscillogram are calculated, thereby obtaining a numerical array of the sample period and a plurality of numerical arrays of all periods of the oscillogram of the same dimension; alternately compare all arrays of periods with a sample, and the first element of each numeric array is compared with the first element of the sample, the second with the second, the third with the third, and so on, in turn, until the last elements of the array; then calculate the values of the differences of elements, that is, the deviation from the sample along the ordinate; then the number of deviations for each analyzed period of the oscillogram is calculated according to five criteria, the value of the difference along the ordinate is in the range: from 5 to 10, from 10 to 20, from 20 to 30 and from 60 to ∞ percent of A max . Criteria are introduced for the number of exits of the elements of the waveform period arrays beyond each of the five ranges along the ordinate, according to which the computer system must decide to report and form images of the analyzed waveform period or not, that is, ignore the period, considering that it does not differ enough from a given sample, and the computer system reports on the period when at least one range is exited and indicates which ranges have been exceeded; or leave the default criteria, characterized in that the input parameters Y max and Y min are calculated as follows: calculate the coefficient of the tenth part of the waveform S dc by multiplying S by 0.1 S dc = S × 0.1, round S dc to an integer value, represent the initial section oscillograms in the form of an array of numerical values of the oscillogram points from its first element to an element with index S dc , assign the designation M s to this array, divide the resulting array into 10 equal parts, if the array has 10 or more elements, if the elements in the array are less than 10, then divide M s into the number of parts equal to the number of elements M s , if division into 10 equal parts is impossible due to the indivisibility of the number of elements M s without a remainder, then the residual elements of the array are distributed between arbitrary parts, but no more than one element per part; calculate the minimum and maximum value of the element along the ordinate in each part, then calculate the arithmetic mean Y min for the minimum values of all parts of the array M s and calculate the arithmetic mean Y max for the maximum values of all parts of the array M s .
RU2019109299A 2019-03-29 2019-03-29 Waveform analysis method RU2727295C1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2019109299A RU2727295C1 (en) 2019-03-29 2019-03-29 Waveform analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2019109299A RU2727295C1 (en) 2019-03-29 2019-03-29 Waveform analysis method

Publications (1)

Publication Number Publication Date
RU2727295C1 true RU2727295C1 (en) 2020-07-21

Family

ID=71741185

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2019109299A RU2727295C1 (en) 2019-03-29 2019-03-29 Waveform analysis method

Country Status (1)

Country Link
RU (1) RU2727295C1 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1408443A1 (en) * 1986-12-08 1988-07-07 Всесоюзный Электротехнический Институт Им.В.И.Ленина Device for simulating pulsed interference
US20050049844A1 (en) * 2003-08-25 2005-03-03 Legend Design Technology, Inc. Glitch and metastability checks using signal characteristics
US20100052653A1 (en) * 2008-08-26 2010-03-04 Spx Corporation Digital Oscilloscope Module with Glitch Detection
RU2658101C1 (en) * 2017-03-29 2018-06-19 Николай Евгеньевич Моргаевский Oscillographic measuring device of amplitude characteristics of electric signals

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1408443A1 (en) * 1986-12-08 1988-07-07 Всесоюзный Электротехнический Институт Им.В.И.Ленина Device for simulating pulsed interference
US20050049844A1 (en) * 2003-08-25 2005-03-03 Legend Design Technology, Inc. Glitch and metastability checks using signal characteristics
US20100052653A1 (en) * 2008-08-26 2010-03-04 Spx Corporation Digital Oscilloscope Module with Glitch Detection
RU2658101C1 (en) * 2017-03-29 2018-06-19 Николай Евгеньевич Моргаевский Oscillographic measuring device of amplitude characteristics of electric signals

Similar Documents

Publication Publication Date Title
CN1230683C (en) Simultaneous display of primary measurement values and derived parameters
US20160034103A1 (en) Capacitance type touch panel and method of detecting input operation position therein
US10852323B2 (en) Measurement apparatus and method for analyzing a waveform of a signal
US6989833B2 (en) Waveform information display apparatus
Sewell et al. Attentional control in visual signal detection: Effects of abrupt-onset and no-onset stimuli.
CN105510664B (en) A kind of automatic setting method of digital oscilloscope
CN109558952A (en) Data processing method, system, equipment and storage medium
RU2727295C1 (en) Waveform analysis method
Mandal et al. Using eye movement data visualization to enhance training of air traffic controllers: A dynamic network approach
CN102227624B (en) System and method for displaying three-dimensional object scattergrams
Murray et al. Posterior probability matching and human perceptual decision making
CN100399033C (en) Surface mapping and 3-D parametric analysis
KR101908786B1 (en) Data similarity evaluation system
Schwartz et al. Probability maps for evaluating automated visual fields
RU2684203C1 (en) Method for intellectual analysis of oscillogram
CN104052541A (en) Display method and system for optical fiber detection and optical time domain reflectometer
RU2786156C1 (en) Method for detection of anomalies in shape of electrical signal
CN109145728B (en) Afterglow blanking method and device and computer readable storage medium
JP7458641B2 (en) Method for estimating electrode placement in biological tissue
CN108254604B (en) Waveform display method and device
RU2741762C1 (en) Method of training system for recognizing signals non-monotonicity
DE17170141T1 (en) Detection and visualization of correlations between measured correlation values and correlation of reference values of a path
RU2204167C2 (en) Method for on-line analysis of multiparametric object condition
JP2021071332A (en) Method and system for determining arrival of earthquake motion
RU2150742C1 (en) Method for color representation and analysis of dynamic state of object or process with multiple parameters