Материал: Sb98338

ЛАБОРАТОРНАЯ РАБОТА № 1

АППРОКСИМАЦИЯ ОБРАТНОЙ КРИВОЙ НАМАГНИЧИВАНИЯ ЭЛЕКТРИЧЕСКОЙ МАШИНЫ НА ОСНОВЕ МЕТОДА НАИМЕНЬШИХ КВАДРАТОВ

Цель работы: аппроксимировать нелинейную зависимость F(Ф) заданную таблично, в промежуточных точках; аппроксимирующую функцию найти в виде полинома заданной степени; оценить зависимость точности аппроксимации от степени полинома.

Постановка задачи

В последующих лабораторных работах понадобится использование функции Ф–1(F), т. е. функции F(Ф), заданной таблично для некоторого числа точек N. При отыскании статических режимов, при численном интегрировании СНЛДУ, а также при линеаризации необходимо иметь значения F при любом значении Ф из диапазона [0, Фmax], а не только в конечном числе

точек. Поэтому нужно найти некоторую функцию р(Ф), определенную на всем отрезке [0, Фmax], обладающую свойством p(Фi ) F(Фi ) в каждой i-й точке

таблицы и аппроксимирующую таблицу в промежуточных точках. В данной работе использована функция р(Ф) в виде полинома заданной степени n:

p(x) c0 c1x cn xn ,

где x – переменная.

Значение n выбирают из условия обеспечения требуемой точности. Чем больше n, тем выше точность аппроксимации. При заданном значении п полином р определяют его коэффициенты c0 , ,cn . Метод наименьших квадра-

тов (МНК) позволяет аналитически определить коэффициенты полинома, когда критерием точности является функционал вида

i 1

где N – количество точек в таблице. В выражении (1.1) разность

i Fi p Фi

называется невязкой (i -й точки), поэтому значение функционала I есть сумма квадратов значений невязок по всем N точкам. Метод наименьших квадратов дает формулу для определения коэффициентов полинома, при которых значение функционала I будет наименьшим.

6

Содержание работы

1. Рассмотреть таблицу F(Ф), которую следует переписать в виде F (Ф) , где Ф – нормированный магнитный поток Ф Ф/ Фн, F – нормированная МДС F F / Fн. Такое преобразование необходимо с целью обеспечения хорошей обусловленности матриц в формуле МНК. Аппроксимации подлежит таблица F (Ф) . Ее аппроксимирующий полином обозначить через p . После нахождения этого полинома коэффициенты р легко определить путем перехода от Ф к Ф и от F к F.

2. Написать и отладить программу расчета коэффициентов полинома p , т. е. с0, ,cn , затем рассчитать коэффициенты полинома р.

Предварительно нужно выбрать значение степени полинома п. Следует учитывать, что график F (Ф) симметричен относительно начала координат ( F (Ф) – нечетная функция), следовательно, полином p будет иметь нулевые коэффициенты при нечетных степенях переменной, а также с0 0 . Поэтому рассчитывают только коэффициенты с1,с3, ,сn, где n должно быть нечетным.

Алгоритм расчета коэффициентов полинома p :

1) построить матрицу G G1,G3, ,Gn , в которой каждый столбец Gi имеет вид

Ф1i Gi ,ФNi

где показатель степени i =1, 3, 5, ... , п ;

2) рассчитать вектор C GTG 1 GT F , где F – вектор-столбец из таблицы F (Ф) . Компоненты вектора C есть коэффициенты полинома p при

нечетных степенях, причем последний элемент вектора С есть коэффициент при старшей степени. В языке МAТLAВ полином представлен вектором, где первый элемент – это коэффициент при старшей степени. Поэтому при п = 3 полином p формируется в программе так:

pc 2 ,0,c 1 ,0 ;

3)сформировать полином р.

4)построить графики F(Ф) и р(Ф) и оценить качество аппроксимации.

7

Рекомендации по программированию

Понадобятся две функции языка MATLAB:

x inv y – обращает квадратную невырожденную матрицу у, результат помещает в х;

z polyval p, x – вычисляет значение полинома р (заданного как век-

тор коэффициентов) на каждом элементе вектора х, результат помещает в соответствующий элемент вектора у, т. е. реализует расчет значений функции у = р(х) для всех точек из х.

Содержание отчета

1.Формулировка цели работы.

2.Таблицы F(Ф)и р(Ф).

3.Алгоритм расчета и соответствующий ему текст программы.

4.Значения коэффициентов аппроксимирующих полиномов p и р для

n = 3 и n = 5.

5. Графики F(Ф) и р(Ф) и выводы о качестве аппроксимации.

8

ЛАБОРАТОРНАЯ РАБОТА № 2

ИССЛЕДОВАНИЕ ПЕРЕХОДНЫХ ПРОЦЕССОВ ДИНАМИЧЕСКОЙ СИСТЕМЫ

Цель работы: исследовать характер переходных процессов, используя численное интегрирование СНДУ объекта; построить модель динамической системы в среде SIMULINK.

Постановка задачи

Исходные уравнения (Вl) описывают полностью динамику объекта. Исследование переходных процессов подразумевает расчет изменения во времени всех переменных в уравнениях (Вl), что может быть осуществлено любым численным методом интегрирования. Однако для решения численным методом необходимо: во-первых, избавиться от алгебраических уравнений (исключив промежуточные переменные); во-вторых, записать дифференциальные уравнения в канонической форме Коши. Принято представлять уравнения динамической системы в виде векторных дифференциального уравнения состояния и алгебраического уравнения выхода:

где x Rn , u Rl ; функция g(x, и) действует из Rn l в Rm и является в общем случае нелинейной. Переменная y Rm есть вектор выходных коорди-

нат, т. е. переменных, значения которых доступны измерению. Непосредственно измерить все компоненты вектора х в большинстве случаев невозможно. Измерениям доступна, как правило, лишь некоторая функция переменных состояния и управления. Значение этой функции в любой момент времени t, т. е. y(t), есть вектор измеряемых переменных – информация, по которой можно судить о том, что происходит в системе.

Пример. Для приведения модели (Вl) к обобщенной форме необходимо

полином р(Ф), полученный в лабораторной работе № 1. Далее ток возбуждения iв выразим через uв и rв. Затем произведем нормирование переменных (в данном примере – относительно номинальных значений параметров):

9

Ф Ф/ Фн, iг iг / iгн, / н, uв uв / uвн, Mв Mв / Mвн.

В итоге преобразованная форма записи уравнений имеет вид:

можно перейти к описанию СНДУ в унифицированном виде (2.1). Стандартные методы и пакеты программ для исследования динамиче-

ских систем ориентированы на такое представление. При этом функция g(x,y) может быть сформирована как угодно.

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

ными во времени значениями x(t) x(0) и u(t) u(0) . Переходный процесс в системе возникает, когда вектор входов принимает новое постоянное значение u(1) u(0) . В этом случае происходит изменение вектора х от значения x(0) к значению x(1) . Переходный процесс в статический режим (u(1) , x(1) ) практически заканчивается за конечное время. Поэтому построение кривых переходного процесса заключается в следующем: для СНДУ (2.1) вектор x(0) принимают в качестве начальных условий x(t0 ) x(0) , далее численным ме-

10