В статье рассматривается алгоритмизация процессов, наблюдаемых при описании вольт-амперных характеристик нелинейных элементов электротехнических устройств, при этом подбирается вид аналитической зависимости, позволяющий прогнозировать критический предел функционирования устройств.
УДК 519.654:519.657 ВАК 05.27.00 DOI: 10.22184/1992-4178.2017.169.8.136.140
В статье рассматривается алгоритмизация процессов, наблюдаемых при описании вольт-амперных характеристик нелинейных элементов электротехнических устройств, при этом подбирается вид аналитической зависимости, позволяющий прогнозировать критический предел функционирования устройств.
УДК 519.654:519.657 ВАК 05.27.00 DOI: 10.22184/1992-4178.2017.169.8.136.140
Существующие программные комплексы трудоемки в эксплуатации и предназначены для использования только высоко квалифицированными специалистами – программистами, электротехниками, электромеханиками. Разработка математической модели нелинейных процессов разного характера на основе единого подхода, достаточно компактной в реализации и обеспечивающей заданную количественную точность, является очень актуальной. Все нелинейные элементы (НЭ) подразделяют на электрические и магнитные. В настоящей работе рассматриваются электрические НЭ, причем для которых как в неуправляемых, так и в управляемых элементах наблюдаются процессы с насыщением, имеющие ВАХ определенного вида [2]. В качестве примера неуправляемых НЭ можно привести бареттер, вакуумный диод, а управляемых – биполярный и полевой транзисторы. Оптимальным для практики проектирования и эксплуатации электротехнических устройств является метод наименьших квадратов (МНК). К его достоинствам относятся: обеспечение требуемой точности для моделирования очень широкого класса функций, простота реализации, относительная независимость от количества точек исходных функций (чем больше точек, тем предпочтительнее), возможность уточнения модели объекта в процессе эксплуатации с учетом новых данных и, таким образом, уточнения прогноза на последующую эксплуатацию [3].
Анализ особенностей нелинейных характеристик наиболее распространенных элементов, ВАХ которых соответствует процессам с насыщением, позволяет выделить два наиболее характерных типа уравнений, которые используются для описания этих характеристик: дробно-рациональные, имеющие вид: или , а также степенные функции: f2 (x) = ахb. При степенные функции подходят для моделирования процессов с насыщением – ВАХ нелинейных резистивных элементов (в частности, бареттера). Для исследования качественного характера поведения ВАХ таких нелинейных элементов приведена таблица значений, отражающая закономерность процессов с насыщением (табл.1). Числовые значения в таблице приведены в качестве примера. Формализованная постановка задачи сводится к подбору вида аналитической зависимости у = f (х), приближающей график (рис.1). Для моделирования таких процессов ограничимся двумя параметрами в уравнениях. Определим параметры а и b методом наименьших квадратов (МНК), согласно этому методу наилучшими параметрами a и b считают те, для которых сумма квадратов отклонений минимальна: . Здесь n – число точек табличной функции (см. табл.1). Необходимо найти частные производные от этой суммы по параметрам a и b и решить уравнения:
I. Рассмотрим дробно-линейную аппроксимирующую функцию . Для нее система имеет вид:
Далее система уравнений преобразуется к виду:
Если ввести следующие обозначения: ; ; ; ; , то систему можно записать в нормальном виде по МНК:
Для вычисления коэффициентов нормальной системы нужно составить таблицу (см. табл.2), причем последние три столбца таблицы – расчетные значения , отклонения δ1i и квадраты отклонений δ1i2 – заполняются только после решения системы. С найденными числовыми коэффициентами система уравнений имеет вид:
Решение этой системы позволяет определить значения параметров: а = 0,0435; b = 0,1126. Тогда искомая функция f1 (x) будет представлена в виде: . Для анализа полученной функции f1 (x) при установленных для нее значениях аргумента вычисляем отклонения и дисперсию . Тогда среднеквадратическое отклонение (СКО): . II. Выполним анализ для степенной аппроксимирующей функции f2(x) = axb. Для нее система уравнений имеет вид:
Система нелинейна относительно неизвестных а и b. Ее решение в полученной форме невозможно. Здесь целесообразно логарифмирование функции f2 (x): y = a xb; ln y = ln a + b ln x. Обозначив z = ln y; A = ln a; p = ln x, можно записать: z = A + b p. Частные производные по А и b равны ; . Теперь систему уравнений можно записать в виде:
Систему уравнений несложно преобразовать к виду:
Если ввести следующие обозначения: то система уравнений будет представлена в виде:
Для вычисления коэффициентов системы нужно составить таблицу (табл.3), причем последние три столбца таблицы – расчетные значения f2 (xi), отклонения δ2i и квадраты отклонений δ22i – заполняются после решения системы. С числовыми коэффициентами систему уравнений запишем в виде (n = 4):
Решение системы позволяет определить значения параметров А = 1,9796; b = 0,4124, откуда а = eA = 7,2399. Тогда искомую функцию f2 (x) можно записать в виде: f2 (x) = 7,2399 x0,4124. Для анализа полученной функции при установленных для нее значениях аргумента вычислим отклонения и дисперсию . Тогда среднеквадратическое отклонение (СКО): . Очевидно, в рассмотренных в настоящей статье задачах аппроксимации по методу наименьших квадратов, погрешности, обеспечиваемые дробно-рациональной функцией f1 (x), и степенной функцией f2 (x) сравнимы по точности. Однако в каждом конкретном случае решение о применении той или иной зависимости следует принимать индивидуально на основе сравнительного анализа результатов аппроксимации. Следует также заметить, что для моделирования процессов с насыщением с горизонтальными асимптотами подходят также показательные (здесь – экспоненциальные) функции вида:
В ряде случаев, как следует из опыта, такие функции имеют преимущество в моделировании ВАХ, показывая бульшую точность приближения при сокращении алгоритмов и вычислительных операций по сравнению с функциями f1(x) и f2 (x). Таким образом, разработан обобщенный алгоритм решения задач аппроксимации процессов с насыщением различными аналитическими функциями. Для широкой номенклатуры нелинейных элементов электротехнических устройств их экспресс-диагностикой может служить анализ кривой ВАХ. Получив экспериментальные данные кривой, в каждом плановом контроле для нее можно установить исследуемую закономерность, анализируя значение СКО. Это позволит на основе достаточной статистики эксплуатации рассматриваемых устройств судить об их техническом состоянии и о продолжительности их последующей эксплуатации, а также прогнозировать критический предел функционирования и остаточное время жизни устройства [4]. ЛИТЕРАТУРА 1. Гольдберг О.Д. Испытания электрических машин. Учебное пособие. – М.: Высшая школа, 2000. 255 с.: ил. 2. Бессонов Л.А. Теоретические основы электротехники. Электрические цепи. Учебник. 10-е изд. – М.: Гардарики, 2002. 638 с. 3. Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа. Приближение функций, дифференциальные и интегральные уравнения. – СПб: Лань, 2010. 400 с. 4. Чернова Т.А., Лисов А.А., Кубрин П.В. Получение консолидированной информации о состоянии контролируемых объектов по совокупности их характеристических параметров // Наукоемкие технологии. 2015. № 6. Т. 16. 16–19 с.