Способ анализа осциллограмм Российский патент 2020 года по МПК G06F17/00 

Описание патента на изобретение RU2727295C1

Изобретение относится к способам вычислений, используемых для обучения компьютерных систем.

Известен способ поиска аномалий в осциллограмме WaveScan от компании Teledyne Lecroy (WaveScan [Электронный ресурс] // Сайт компании Teledyne Lecroy. 2018. URL: http://teledynelecroy.com/features/featureoverview.aspx?modelid=2107&capid=102&mid=556). Способ состоит в том, что воспроизводят осциллограмму целиком, либо частично на дисплее компьютерной системы в виде статического изображения. Затем проводят обучение компьютерной системы путем указания на часть осциллограммы, которую необходимо принять за графический образец для анализа на предмет наличия в осциллограмме таких же частей как образец. Либо проводят обучение компьютерной системы в ручном режиме, путем задания множества числовых параметров, описывающих искомые аномалии сигнала. При обнаружении данных аномалий компьютерная система сообщает об их наличии и сохраняет их изображения.

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

В рамках данной заявки, «Аномалия сигнала» - это период электрического сигнала или его часть, отличающиеся от периодов сигнала, принятых в качестве нормальных.

Наиболее близким (прототипом) является «Способ интеллектуальной синхронизации TriggerScan» от компании Teledyne Lecroy (TriggerScan [Электронный ресурс] // Сайт компании Teledyne Lecroy. 2018. URL: http://teledynelecroy.com/features/featureoverview.aspx?modelid=2108&capid=102&mid=556). Способ состоит в том, что воспроизводят осциллограмму целиком, либо частично на дисплее компьютерной системы в виде статического изображения. Затем проводят обучение компьютерной системы путем указания на часть осциллограммы, которую необходимо принять за графический образец. Далее компьютерная система на основе заданного графического образца автоматически рассчитывает десятки числовых параметров описывающих сигнал, принятый за образец и, на основе данных параметров, может обнаруживать некоторые части осциллограммы, отличные от заданного графического образца, а также сообщать о них специалисту, проводящему анализ, и сохранять их изображения.

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

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

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

Техническая проблема изобретения решается путем применения алгоритма расчета входных параметров поиска аномалий электрического сигнала, который состоит в следующем: подсчитывают общее количество точек 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. Специалист, проводящий анализ, вводит критерии количества выходов элементов массивов периодов осциллограммы за пределы каждого из пяти диапазонов по оси ординат, в соответствие с которыми компьютерная система должна принять решение сообщать и формировать изображения анализируемого периода осциллограммы или не делать этого, то есть проигнорировать период, посчитав, что он недостаточно отличается от рассчитанного образца, причем компьютерная система сообщает о периоде при выходе, хотя бы из одного диапазона и обозначает какие именно диапазоны были превышены. Также имеется возможность оставить критерии по умолчанию и тем самым исключить работу специалиста, возложив ее, например, на ЭВМ.

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

В рамках данной заявки, «Период осциллограммы» - графическое отображение периода сигнала на осциллограмме.

Заявка поясняется изображениями:

Фиг. 1 - изображение осциллограммы синусоидальной формы, используемой для осуществления предлагаемого способа с обозначениями: УОО низких, средних и высоких амплитуд, единственного восходящего перехода от УОО низких к УОО высоких амплитуд и периода сигнала;

Фиг. 2 - изображение осциллограммы треугольной формы, используемой для осуществления предлагаемого способа с обозначениями: УОО низких, средних и высоких амплитуд, единственного восходящего перехода от УОО низких к УОО высоких амплитуд и периода сигнала;

Фиг. 3 - изображение осциллограммы прямоугольной формы, используемой для осуществления предлагаемого способа с обозначениями: УОО низких, средних и высоких амплитуд, единственного восходящего перехода от УОО низких к УОО высоких амплитуд и периода сигнала.

Способ применим для анализа любых типов сигналов, которые имеют в своем составе периодически повторяющуюся часть, а также один, и только один, восходящий переход от условно обозначенной области низких амплитуд к условно обозначенной области высоких амплитуд в пределах каждого отдельного периода. К такому типу относится большинство сигналов встречающихся в природе, в том числе сигналы синусоидальной (фиг. 1), треугольной (фиг. 2) и прямоугольной формы (фиг. 3).

Способ осуществляют следующим образом.

Разрабатывают программу для ЭВМ в которой предусматривают функцию вывода данных из файлов, содержащих осциллограмму электрического сигнала. Подсчитывают общее количество точек в осциллограмме - 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 элементу в каждой.

Вычисляют минимальную и максимальную величину элемента по оси ординат в каждой части, затем вычисляют среднее арифметическое 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. Специалист, проводящий анализ, вводит критерии количества выходов элементов массивов периодов осциллограммы за пределы каждого из пяти диапазонов по оси ординат, в соответствие с которыми компьютерная система должна принять решение сообщать и формировать изображения анализируемого периода осциллограммы или не делать этого, то есть проигнорировать период, посчитав, что он недостаточно отличается от рассчитанного образца, причем компьютерная система сообщает о периоде при выходе, хотя бы из одного диапазона и обозначает какие именно диапазоны были превышены. Также имеется возможность оставить критерии по умолчанию и тем самым исключить работу специалиста, возложив ее, например, на ЭВМ.

Похожие патенты RU2727295C1

название год авторы номер документа
СПОСОБ ИНТЕЛЛЕКТУАЛЬНОГО АНАЛИЗА ОСЦИЛЛОГРАММ 2017
  • Недорезов Дмитрий Александрович
RU2684203C1
СПОСОБ ОБУЧЕНИЯ СИСТЕМЫ РАСПОЗНАВАНИЯ НЕМОНОТОННОСТИ СИГНАЛОВ 2020
  • Недорезов Дмитрий Александрович
RU2741762C1
СПОСОБ ОБНАРУЖЕНИЯ АНОМАЛИЙ ФОРМЫ ЭЛЕКТРИЧЕСКОГО СИГНАЛА 2021
  • Недорезов Дмитрий Александрович
RU2786156C1
ОБЪЕДИНЕННОЕ СООБЩЕНИЕ HARQ 2008
  • Бориу Эдриан
  • Махешвари Шашикант
  • Саифуллах Юсуф
RU2456754C2
Способ обнаружения импульсных помех 2019
  • Недорезов Дмитрий Александрович
RU2731320C1
Способ клонального микроразмножения in vitro сортового хмеля 2021
  • Хлебова Любовь Петровна
  • Бровко Елена Сергеевна
  • Мироненко Ольга Николаевна
  • Бычкова Ольга Владимировна
  • Хлыновский Михаил Данилович
RU2777200C1
СПОСОБ ИССЛЕДОВАНИЯ ПРОЦЕССА РАЗРУШЕНИЯ ТВЕРДОГО ТЕЛА ПОД ДЕЙСТВИЕМ НАГРУЗКИ ПО ЕГО ЭЛЕКТРОМАГНИТНОМУ ИЗЛУЧЕНИЮ 2007
  • Опарин Виктор Николаевич
  • Кулаков Геннадий Иванович
RU2343282C1
ОПТИЧЕСКОЕ УСТРОЙСТВО, ОБЕСПЕЧИВАЮЩЕЕ ПО МЕНЬШЕЙ МЕРЕ ДВЕ ОПТИЧЕСКИЕ ФУНКЦИИ 2018
  • Борискин, Артём
  • Шрамкова, Оксана
RU2755257C1
СПОСОБ И СЕТЬ, ПОДХОДЯЩИЕ ДЛЯ УВЕЛИЧЕНИЯ SINR КАНАЛА ПЕРЕДАЧИ ДАННЫХ 2008
  • Ким Дзее Хиун
  • Шульц Эгон
  • Цирвас Вольфганг
RU2452093C2
СПОСОБ ОБНАРУЖЕНИЯ АНОМАЛЬНЫХ ИЗМЕРЕНИЙ БЕЗ ОЦЕНКИ ФУНКЦИИ ТРЕНДА И УСТРОЙСТВО, ЕГО РЕАЛИЗУЮЩЕЕ 2005
  • Марчук Владимир Иванович
  • Шерстобитов Александр Иванович
  • Воронин Вячеслав Владимирович
  • Токарева Светлана Викторовна
RU2302655C1

Иллюстрации к изобретению RU 2 727 295 C1

Реферат патента 2020 года Способ анализа осциллограмм

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

Формула изобретения RU 2 727 295 C1

Способ анализа осциллограмм, при котором вычисляют минимальное значение амплитуды образца сигнала 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.

Документы, цитированные в отчете о поиске Патент 2020 года RU2727295C1

Способ обработки целлюлозных материалов, с целью тонкого измельчения или переведения в коллоидальный раствор 1923
  • Петров Г.С.
SU2005A1
Приспособление для суммирования отрезков прямых линий 1923
  • Иванцов Г.П.
SU2010A1
ОСЦИЛЛОГРАФИЧЕСКИЙ ИЗМЕРИТЕЛЬ АМПЛИТУДНЫХ ХАРАКТЕРИСТИК ЭЛЕКТРИЧЕСКИХ СИГНАЛОВ 2017
  • Моргаевский Николай Евгеньевич
RU2658101C1
Устройство для моделирования импульсных помех 1986
  • Евликов Александр Александрович
  • Лазутин Геннадий Тимофеевич
  • Матвеев Владимир Васильевич
SU1408443A1

RU 2 727 295 C1

Авторы

Недорезов Дмитрий Александрович

Даты

2020-07-21Публикация

2019-03-29Подача