Курсовой ЭТФ

Вариант 9

В результате исследований, как должна меняться нагрузка при изменении напряжения, чтобы температура узла АД не превышала допустимого значения.

Таблица. Влияние тока холостого хода на величину нагрузки при постоянной температуре и изменении напряжения

Создайте текстовый файл ad.txt с исходными данными. Пользуясь его содержимым, определите предельную нагрузку для АД с I*0 0.5 при напряжении 0,77U1н, и для АД с I0* 0.6 при напряжении 0,9U1н и 0,85U1н. Функциональные зависимости β(U1) и β(I0) построить с помощью сплайн-интерполяции.

Впрограмме исходные данные рекомендуется считывать из файла ad.txt или

склавиатуры. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц проверить результаты работы программы, построить график экспериментальной зависимости и подобранной функции.

Вариант 10

В результате исследований, как должна меняться нагрузка при изменении напряжения, чтобы температура узла АД не превышала допустимого значения.

Таблица. Влияние тока холостого хода на величину нагрузки при постоянной температуре и изменении напряжения

Создайте текстовый файл ad.txt с исходными данными. Пользуясь его содержимым, определите предельную нагрузку для АД с I*0 0.5 при напряжении 0,77U1н, и для АД с I0* 0.6 при напряжении 0,9U1н и 0,85U1н. Функциональные зависимости β(U1) и β(I0) построить с помощью линейной интерполяции.

Впрограмме исходные данные рекомендуется считывать из файла ad.txt или

склавиатуры. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц проверить результаты работы программы, построить график экспериментальной зависимости и подобранной функции.

41

Вариант 11

В результате исследований, как должна меняться нагрузка при изменении напряжения, чтобы температура узла АД не превышала допустимого значения.

Таблица. Влияние тока холостого хода на величину нагрузки при постоянной температуре и изменении напряжения

Создайте текстовый файл ad.txt с исходными данными. Пользуясь его содержимым, определите предельную нагрузку для АД с I*0 0.5 при напряжении 0,77U1н, и для АД с I0* 0.6 при напряжении 0,9U1н и 0,85U1н. Функциональные зависимости β(U1) и β(I0) построить с помощью сплайн-интерполяции.

Впрограмме исходные данные рекомендуется считывать из файла ad.txt или

склавиатуры. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц проверить результаты работы программы, построить график экспериментальной зависимости и подобранной функции.

Вариант 12

В результате исследований, как должна меняться нагрузка при изменении напряжения, чтобы температура узла АД не превышала допустимого значения.

Таблица. Влияние тока холостого хода на величину нагрузки при постоянной

температуре и изменении напряжения

Создайте текстовый файл ad.txt с исходными данными. Пользуясь его содержимым, определите предельную нагрузку для АД с I*0 0.7 при напряжении

0,77U1н, и для АД с I0* 0.4 при напряжении 0,9U1н и 0,85U1н. Функциональные зависимости β(U1) и β(I0) построить с помощью линейной интерполяции.

Впрограмме исходные данные рекомендуется считывать из файла ad.txt или

склавиатуры. Результаты вычислений сохранить в текстовом файле rezult.txt. Средствами математического пакета или электронных таблиц проверить результаты работы программы, построить график экспериментальной зависимости и подобранной функции.

42

Вариант 13

Известная табличная зависимость пускового тока статора I1п(А) от входного напряжения U1 (В) для общепромышленного асинхронного двигателя АД АИУМ225М4, которая хранится в текстовом файле aium225m4.txt.

Построить интерполяционный полином Лагранжа и линейный сплайн зависимости I1(β). Найти значение пусковых токов при напряжениях 0.88U1н, 0.93U1н, 1.25U1н. Результаты вычислений сохранить в текстовом файле rezult.txt.

В программе исходные данные можно считывать с клавиатуры или из файла aium225m4.txt. Средствами математического пакета или электронных таблиц проверить результаты работы программы, построить график экспериментальной зависимости и подобранной функции.

43

Первый этап решения системы уравнений методом Гаусса, называемый прямым ходом метода Гаусса, заключается в приведении расширенной матрицы (3.2) к треугольному виду. Это означает, что все элементы матрицы (3.2) ниже главной диагонали должны быть равны нулю

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

2-я строка = 2-я строка – М 1-я строка

3-я строка = 3-я строка – М 1-я строка

………………………..

i-я строка = i-я строка – М 1-я строка

………………………….

n-я строка = n-я строка – М 1-я строка

Преобразование элементов второй строки будет происходить по формулам:

a21 a21 Ma11

a22 a22 Ma12

…

a2i a2i Ma1i …

a2n a2n Ma1n

b2 b2 Mb1 .

Так как целью данных преобразований является обнуление первого элемента строки, то М выбирается из условия a21 Ma11 0 .

a

Следовательно, M 21 .

a11

Элементы третьей строки и коэффициент М можно рассчитать аналогично:

a31 a31 Ma11

45

a32 a32 Ma12

…

a3i a3i Ma1i …

a3n a3n Ma1n

b3 b3 Mb1,

a31 Ma11 0 M a31 .

a11

Таким образом, преобразование элементов i–й строки будет происходить следующим образом:

ai1 ai1 Ma11

ai2 ai2 Ma12

…

aii aii Ma1i

…

ain ain Ma1n

bi bi Mb1 .

Коэффициент М для i–й строки выбирается из условия ai1 Ma11 0, и

равен M aai1 .

11

После проведения подобных преобразований для всех строк матрица (3.2) примет вид:

Очевидно, что если повторить описанный выше алгоритм для следующих столбцов матрицы (3.2), причем преобразование второго столбца начнется с третьего элементы, преобразование третьего столбца – с четвертого элемента и т. д., то в результате будет получена матрица (3.3).

Если в матрице (3.2) на главной диагонали встретится элемент akk, равный

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

46

В результате выполнения прямого хода метода Гаусса матрица (3.2) преобразуется в матрицу (3.3), а система уравнений (3.1) будет иметь вид

...

Решение системы (3.4) называют обратным ходом метода Гаусса.

Последнее n–е уравнение системы (3.4) имеет вид:

(3.4), а следовательно, и система (3.1) имеют бесконечное множество решений. При ann 0 и bn 0 система (3.4), а, значит, и система (3.1) решения не

имеет. Предпоследнее (n–1)–е уравнение системы (3.4) имеет вид:

an 1n 1xn 1 an 1n xn bn 1.

Следовательно,

xn 1 bn 1 an 1n xn .

Следующее, (n–2)–е уравнение системы (3.4) будет иметь вид:

an 2n 2 xn 2 an 2n 1xn 1 an 2n xn bn 2 .

Отсюда имеем

Таким образом, формула для вычисления i–го значения x будет иметь вид:

n

bi aij xj

47

Метод LU–разложения

Метод LU-разложения является точным методом решения линейных систем. Решение систем линейных уравнений методом LU–разложения заключается в следующем. Пусть дана система

Ax=b. ( 3.5)

Представим матрицу A={aij} в виде произведения треугольных матриц

L={lij} и U={uij}:

где элементы lij и uij определяются по формулам:

Формулы (3.6) и (3.7) дают нам LU–разложение матрицы А. Итак, если А=L U, то можно записать уравнение эквивалентное (3.5):

L U x=b. ( 3.8)

Введем вектор вспомогательных переменных y=(y1, y2, …, yn)T и представим уравнение (3.8) в виде системы

L x b,

U x y

Таким образом, решение системы (3.1) с квадратной матрицей коэффициентов, сводится к решению двух систем с треугольными матрицами коэффициентов:

и

48

Из (3.9) и (3.10) понятно, что все yi могут быть вычислены по формуле

Итак, решить систему (3.1) методом LU–разложения можно, выполнив следующие действия:

Определить элементы матриц L и U по формулам (3.6) и (3.7).

Получить решение заданной системы по формулам (3.11) и (3.12).

Метод QR–разложения

Метод QR–разложения можно считать более сложной версией LU–разложения, так как в нем главную матрицу A линейной системы (3.5) так же представляют в виде произведения двух матриц – ортогональной матрицы Q и верхней треугольной матрицы R. Далее поступают так же, как и в случае LU– разложения, то есть вводят дополнительные переменные y и записывают систему уравнений

Q R x=b (Q y=b и R x=y).

Следовательно, решение системы уравнений (3.1) сводится к решению двух систем, матрица коэффициентов первой – ортогональная, второй – верхняя треугольная. А поскольку ортогональная матрица обладает свойством QТ=Q–1, то у можно вычислить по формуле у =QT b, что значительно проще чем решать уравнение у =Q–1 b. Вычислить x из второй системы так же не составит труда, это можно сделать, например, по формулам (3.12), заменив элементы матрицы U на R

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

Метод ортогонализации

Методом ортогонализации является точным методом решения системы линейных алгебраических уравнений. Заключается он в том, что систему (3.1) ортонормируют, то есть приводят ее к виду

где Q – ортогональная матрица. Для этого разделим первое уравнение системы

n

(3.1) на выражение a12k

k 1

49

Получим

Если система (3.1) совместна, то, воспользовавшись формулами (3.14) и (3.15) можно прийти к новой системе (3.13), где матрица Q будет ортогональной, а это значит, что она обладает свойством

QT Q=E,

где Е – единичная матрица. Значит, решение системы (3.13) можно записать в виде

QT Q x=QT d,

или

Метод Зейделя

Метод Зейделя – это приближенный метод решения систем уравнений,

представляющий собой модификацию метода простой итерации. Основная идея этого метода заключается в том, что при вычислении (к+1)–го приближения неизвестной xi учитываются уже вычисленные к этому времени (к+1)–е предшествующих неизвестных x1, x2, …, xi–1.

Пусть дана приведенная линейная система

Выберем в качестве начального приближения корней вектор

50