WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     | 1 |   ...   | 3 | 4 || 6 | 7 |   ...   | 8 |

Таблица 6.Варианты заданий для самостоятельной работы № вар. Уравнение № вар. Уравнение 1 x3 - 2 x2 = x - 56 9 37 x - x3 = 2 2 x2 + x = 69 10 x3 + 73 x4 = x - 3 2 x4 - x3 = 0 11 x3 + 2 x = x4 5 x2 + 26 x = 6 12 27 x2 + 132 x = - 5 x3 - 2 x = 4 13 7 x2 - x = 6 2 x5 + x3 = 9 x2 14 x2 + 2 x = -7 24 x2 + x = 6,5 15 2 x4 - 37 x3 = 8 x2 + 26 x = 138 16 x2 - 26 x = Требования к отчёту Отчет о лабораторной работе должен включать цель работы, реферат первого раздела и протокол действий, самостоятельно выполненных студентом на компьютере. Рабочий документ выполнения лабораторной работы должен быть сохранён на ПЭВМ в личной папке студента.

Исходные данные:

W q 60000. T 1300.K 0.p max mW W.

3.5.10 35.

gkm.K4.K m2 m4..

Функция: f T T T T T q g gkm g p g p Построение графика для приближённого определения корня уравнения f Tg T 0.K, 100.K.. 2000.K g 5 2.5 f Tg 2.5 5 0 500 1000 1500 T g Принимаем в нулевом приближении T 1490.K gВычисление производной от f Tg 4..

T T T T q gkm g p g p.T 4.

gkm g Программный модуль для расчёта T методом Ньютона:

g T pog g T T g0 gwhile pog max 4..

T T T T q gkm g0 p g0 p T T gk g.T 4.

gkm gT T gk gpog 100.

T gk T T g0 gk T gk Результат: T = 1.45 103 K g Рис. 6.2. Решение задачи методом Ньютона При сдаче лабораторной работы студент должен ответить на следующие контрольные вопросы:

1. На чём основывается метод Ньютона 2. Напишите итерационную формулу метода Ньютона.

3. Каков геометрический смысл метода Ньютона 4. Какое условие сходимости ньютоновских итераций 5. Каковы преимущества метода Ньютона по сравнению с другими численными методами решения нелинейных уравнений 6. Зачем при нахождении решения нелинейного уравнения методом Ньютона строят график этого уравнения 7. Расскажите последовательность действий при решении нелинейного уравнения методом Ньютона.

8. Как найти производную функции в MathCAD 9. Как построить декартов график в MathCAD 10. Как приближённо найти корень функции на графике с использованием окна Trace 11. Объясните назначение кнопок ‘‘Add Line’’, ‘‘’’ и ‘‘while’’ панели программирования.

12. Что является условием окончания расчётов при использовании метода Ньютона 13. Какой тип индексов используется в примере (раздел 2): обычный или декоративный Чем они отличаются Лабораторная работа № Решение систем линейных алгебраических уравнений методом Гаусса Цель работы 1. Освоить метод Гаусса решения систем линейных алгебраических уравнений.

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

1. Общие сведения Система линейных уравнений записывается следующим образом:

a1,1 x1 + a1,2 x2 +... + a1,n xn = ba2,1 x1 + a2,2 x2 +... + a2,n xn = b, (1)........................................................

an,1 x1 + an,2 x2 +... + an,n xn = bn или в матричной форме:

A x = B, (2) где A = (ai, j ) i, j =1, 2,..., n – матрица коэффициентов при неизвестных, B = (b1, b2,..., bn ) – вектор свободных членов, x = (x1, x2,..., xn ) – вектор неизвестных.

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

Рассмотрим систему уравнений (1). Алгоритм исключения состоит из последовательных шагов.

Первый шаг. Исключение неизвестной х1 из всех уравнений системы (1), кроме первого. Пусть а1,1 0. Тогда из первого уравнения выражаем х x1 = (- a1,2 x2 -... - a1,n xn + b1) (3) a1,через остальные неизвестные и подставляем (3) во 2-е, 3-е, …, n-е уравнение системы (1). После подстановки получаем преобразованную исходную систему вида a1,1 x1 + a1,2 x2 + a1,3 x3 +... + a1,n xn = b[ [ [ [ a21] x2 + a21] x3 +... + a21] xn = b21],2,3,n, (4)....................................................

[ [ [ an 1] x2 + an 1] x2 +... + an[1] xn = bn1],2,3,n [ где элементы матрицы ai, 1], 2 i, j n, и вектора bi[ 1], 2 i n, после первоj го шага получены по формулам 1 [ [ [ [ [ [ ai,k +1] = ai,k ] - ai,k ] akk ], bi[k +1] = bi[k ] - ai,k ] bkk ]. (5) j j [ [ akk ] k, j akk ] k,k,k где [k + 1] – номер шага исключения переменной, k = 0 (n – 1). Номер шага k = 0 отвечает записи системы (1).

Продолжая этот процесс исключения при условии, что [ [ a21] 0,..., ann-2] 0,,2 -1,n-после (n – 1)-го шага исключения получим преобразованную исходную систему в виде a1,1 x1 + a1,2 x2 + a1,3 x3 +... + a1,n xn = b[ [ [ [ a21] x2 + a21] x3 +... + a21] xn = b21],2,3,n................................................ (6) [ [ [ ann-2] xn-1 + ann-2] xn = bnn-2] -1,n-1 -1,n -[ [ ann-1] xn = bnn-1],n Преобразование исходной системы линейных уравнений к системе (6) – прямой ход исключения. На этом этапе мы ещё не вычислили ни одной компоненты вектора решения х, но эквивалентными преобразованиями привели систему к такой форме, для которой легко вычислить все компоненты решения х.

[ Пусть ann-1] 0. Тогда осуществляем обратный ход: вычисляем компо,n ненты вектора решения в обратном порядке.

Из (6) находим [ bnn -1] xn =, [ ann -1],n [ [ xn -1 = (bnn-2] - ann-2] xn ), [ (7) ann-2 ] -1 -1,n -1,n-.....................................................

x1 = (b1 - a1,n xn -... - a1,2 x2 ) a1,Число арифметических действий, выполняемых при решении системы из ‘‘n’’ линейных алгебраических уравнений методом Гаусса, в общем случае, равно n3.

2. Пример применения метода Гаусса Решить задачу:

Дано: куб (n = 6 см. рис. 7.1), внутренние поверхности которого имеют следующие температуры: Т1 = 800 К, Т2 = 900 К, Т3 = 1000 К, Т4 = 1100 К, Т5 = 1200 К, Т6 = 1300 К. Степень черноты всех поверхностей = 0.8;

коэффициент излучения абсолютно чёрного тела 0 = 5.668710 –8 Вт/(м2К4); угловой коэффициент излучения между всеми поверхностями равен = 0.2. Угловой коэффициент ‘‘сам на себя’’ равен нулю: ii = 0.

Найти: плотность эффективного излучения каждой поверхности.

Рис. 7.1. Схема рассРасчётные выражения:

читываемой системы Плотность потока эффективного излучения для i – й поверхности определяется по формуле:

эф qiэф = (1 - ) + Тi4 (8) q j ji j=Всего в кубе 6 поверхностей. Значит, имеем 6 уравнений. Их можно преобразовать к следующему виду (индекс ‘‘эф’’ для упрощения записи опущен):

1/[(1 - )] q1 - q2 - q3 - q4 - q5 - q6 = T14 /[(1 - )] - q1 + 1/[(1 - )] q2 - q3 - q4 - q5 - q6 = T24 /[(1 - )] - q1 - q2 + 1/[(1 - )] q3 - q4 - q5 - q6 = T34 /[(1 - )] -q1 -q2 - q3 +1/[(1-)]q4 -q5 - q6 = 0 T44 /[(1-)] - q1 - q2 - q3 - q4 + 1/[(1 - )] q5 - q6 = T54 /[(1 - )] - q1 - q2 - q3 - q4 - q5 + 1/[(1 - )] q6 = T64 /[(1 - )] Фрагмент рабочего документа с решением этой системы относительно q1, q2, q3, q4, q5 и q6 методом Гаусса представлен на рис. 7.2.

Указания к решению задачи:

1) Системная переменная ORIGIN определяет номер первого элемента массива (по умолчанию ORIGIN = 0). В нашем случае удобнее, чтобы все массивы начинались с элемента № 1;

2) Для облегчения решения задачи размерности лучше опустить;

3) Температуры поверхностей заданы в виде матрицы – строки. Это удобно сделать для придания записи компактного вида. Шаблон матрицы можно вызвать нажатием кнопки с изображением квадратной матрицы на панели векторов и матриц 4) Надо иметь в виду, в расчётах MathCAD оперирует с одномерными массивами – столбцами. Поэтому выполняем транспонирование матрицы. Знак транспонирования матрицы Т можно вызвать нажатием комбинации клавиш [Ctrl + 1] или нажав кнопку ‘‘МТ’’ на панели векторов и матриц;

5) Проверку правильности решения можно произвести при помощи одной из встроенных функций решения систем линейных уравнений, например, lsolve, либо при помощи обратной матрицы. Стандартная функция lsolve имеет вид: lsolve(А, В), где А – матрица коэффициентов при неизвестных; В – матрица–столбец свободных членов.

Номера первых элементов массивов :

ORIGIN.10 Исходные данные: 0.8 5.6687 0.2 n Транспонирование матрицы-строки в матрицу-столбец для удобства ввода массива Ti :

T ( 800 900 1000 1100 1200 1300 )T Формирование матрицы коэффициентов при неизвестных и вектора свободных членов:

i 1.. n j 1.. n A i, j.

Ti 4 0.

A Bi i, i..

( 1 ) ( 1 ) 4.644 25 1 1 1 1 7.438 1 25 1 1 1 1 1 25 1 1 1 1.134 A = B = 1 1 1 25 1 1.66 1 1 1 1 25 2.351 1 1 1 1 1 3.238 Программный модуль для решения системы уравнений методом Гаусса:

q a A b B for k 1.. n for i k 1.. n for j k 1.. n ai, j ai, j.ai, k.ak, j ak, k bi bi.ai, k.bk ak, k bn xn an, n for k n 1.. s for i n.. k s s ak, i.xi xk. bk s ak, k x Проверка 1 (с использованием Проверка 2 (с использованием Результат:

стандартной функции lsolve): обратной матрицы):

3.63 104 3.63 104 3.63 4.705 104 4.705 104 4.705 6.205 104 6.205 104 6.205.B q = lsolve ( A, B) = A = 8.229 104 8.229 104 8.229 1.089 105 1.089 105 1.089 1.43 105 1.43 105 1.43 Рис. 7.2 Решение системы линейных уравнений методом Гаусса Задание для самостоятельной работы Решить систему линейных уравнений вида A x = B методом Гаусса.

Вариант задания взять из таблицы 7.1 согласно порядковому номеру студента по журналу.

Таблица 7.Варианты заданий для самостоятельной работы № вар. А В № вар. А В - 14 13 - - 58 - 0,5 200 9,5 326 0,5 560 3,5 1 9 - 9 9 300 4335 - 6 9 532 - 13 12 - - 23 - 1 400 9 946 3 526 6 2 10 - 8 8 600 7920 -1 8 400 -12 11 - -19 1, - 1,5 600 8,5 5456 1,3 126 3,5 3 11 - 7 7 900 5681 - 2 7 960 -11 10 - -1,1 - 2 800 7 5050 2,8 125 4.1 4 12 - 9 9 354 9520 - 23 1,7 - 22 11 - - 2,2 1, - 4,5 670 7.3 9126 25 45 7.3 5 13 - 5 9,7 300 1923 - 42 4,7 -15 65 - - 0,3 - 8 865 1 3541 8,3 821 1,2 6 14 - 5 21 35,4 9842 - 3 81 15,4 -18 53 - -1,8 - 6 630 9 3281 6,6 621 1 7 15 - 5 62 454 1562 - 3,6 92 - 2 38 - -18 - 3.5 370 4.3 2846 3 340 1 8 16 - 9 9,7 321 1973 - 5 5 154 Требования к отчёту Отчет о лабораторной работе должен включать цель работы, реферат первого раздела и протокол действий, самостоятельно выполненных студентом на компьютере. Рабочий документ выполнения лабораторной работы должен быть сохранён на ПЭВМ в личной папке студента.

При сдаче лабораторной работы студент должен ответить на следующие контрольные вопросы:

1. Как записывается система линейных уравнений в матричной форме 2. В чём суть метода Гаусса 3. Какой вычислительный процесс называется прямым ходом исключения в методе Гаусса 4. Какой вычислительный процесс называется обратным ходом в методе Гаусса 5. Прокомментируйте формулу преобразования матрицы коэффициентов при неизвестных во время прямого хода исключения. Что в неё входит и как получена 6. Прокомментируйте формулу преобразования вектора свободных членов во время прямого хода исключения. Что в неё входит и как получена 7. Прокомментируйте формулу нахождения неизвестных во время обратного хода. Что в неё входит и как получена 8. Какую информацию содержит системная переменная ORIGIN 9. Что происходит при транспонировании матрицы Как осуществить транспонирование матрицы в MathCAD 10. Как можно проверить правильность решения системы линейных уравнений методом Гаусса 11. Объясните назначение кнопок ‘‘Add Line’’, ‘‘’’ и ‘‘while’’ панели программирования.

12. Как транспонировать матрицу–столбец в матрицу–строку Лабораторная работа № Применение метода наименьших квадратов для аппроксимации табличных данных Цель работы 1. Освоить метод наименьших квадратов.

2. Получить практические навыки аппроксимации табличных данных методом наименьших квадратов.

1. Общие сведения Метод наименьших квадратов – основной метод статистической обработки результатов исследования, позволяющий решить, какое из произвольных уравнений даёт наилучшее приближение к фактической зависимости.

Пусть для описания исследуемой зависимости выбрано уравнение y = f(x). Например, y = a0 + a1 x + a2 x2 + a3 x3 + K. Согласно выбранному уравнению y = f(x), значению аргумента xi должно соответствовать значение функции yi. Как правило, получается отклонение фактических значений от расчётных. Сущность метода наименьших квадратов заключается в том, что наилучшее приближение к истинной зависимости даёт такое уравнение, для которого сумма квадратов отклонений экспериментальных и расчётных данных имеет минимальное значение, т.е.

n S = - f (xi )] min. (1) [yi Для сравнения обычно выбирают уравнения определённого типа, но с неизвестными коэффициентами. Величину S можно рассматривать как функцию коэффициентов а1, а2, …, аn искомого уравнения. Задача состоит в том, чтобы найти значения коэффициентов уравнения, соответствующие минимуму S.

Известно, что необходимым условием минимума дифференцируемой функции является равенство нулю первых производных. В данном случае S S S = 0, = 0,..., = 0. (2) a1 a2 an Эти равенства можно рассматривать как систему нормальных уравнений относительно коэффициентов а1, а2, …, аn, которая имеет единственное решение, минимизирующее величину S.

Рассмотрим применение метода наименьших квадратов на примере линейной зависимости типа y = a0+a1x. Отклонение линейного уравнения от искомой прямой в отдельных точках может быть записано в виде системы уравнений:

y1 = y - y1 = a0 + a1 x1 - y1;

y2 = y - y2 = a0 + a1 x2 - y2 ;

(3)...............................................

yn = y - yn = a0 + a1 xn - yn.

Если в каждом частном уравнении системы (3) левую и правую части возвести в квадрат и сложить все уравнения, получим сумму квадратов отклонений n n 2 + a1 xi - yi ). (4) (y - yi ) = (ai=1 i=Приравняем к нулю частные производные:

n yin i= = 2 + a1 xi - yi )= 0;

(a a i= (5) n yi n i= = 2 + a1 xi - yi ) xi ]= 0.

[(a ai= После ряда простейших преобразований получим систему линейных уравнений:

n n yi = n a0 + a1 xi ;

i=1 i= (6) n n n xi yi = a0 xi + a1 xi2.

i=1 i=1 i=Решая (6) любым подходящим методом (методом Гаусса, методом простых итераций и т.д.), находим значения коэффициентов а0 и а1.

2. Пример применения метода наименьших квадратов Решить задачу:

Дано: экспериментально полученная зависимость содержания СО2 в продуктах полного горения топлива от содержания О2 в дутье (дутьё – окислитель топлива на основе воздуха).

Содержание 21 24 27 36 30 42 70 48 О2 в дутье, % СО2, % 11,8 13,6 15,6 22 17,7 26,6 54 31,6 Найти: линейное уравнение вида y = a0 + a1x, описывающее эту зависимость.

Фрагмент рабочего документа с решением данной задачи представлен на рис. 8.1.

Номера первых элементов массивов :

ORIGIN Исходные данные:

Pages:     | 1 |   ...   | 3 | 4 || 6 | 7 |   ...   | 8 |



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.