Формирование кривой намагничивания

В качестве кривой намагничивания использована кривая по первым гармоникам, заданная в виде таблицы (табл. 5.1). Весьма удобно задать в программе кривую намагничивания в виде матрицы, полностью повторяющей табл. 5.1. Первый столбец — магнитная напряженность, второй — индукция. Назовем эту матрицу НВ.

Таблица 5.7.

Кривая намагничивания стали.

Окончание табл. 5.7.

Кривую намагничивания создаем с помощью функции интерполяции interp, которая в качестве одного из аргументов использует функцию pspline. Для обеих этих функций требуются табличные данные в виде векторов-столбцов. Эти столбцы надо вырезать из матрицы НВ:

Обратить внимание! В пакете Mathcad все элементы матриц нумеруются с нулевого индекса!

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

Поскольку при нулевом значении индукции магнитная напряженность также равна нулю, в начальной точке задана магнитная проницаемость, равная 50р_о.