3.3. Решение слау методом Крамера

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

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

Процедура поиска корней позволяет найти решение, если определитель основной матрицы коэффициентов A=(a_ij) не равен нулю. Для нахождения i-го корня ищут определитель, получаемый заменой i-того столбца матрицы А=(a_ij) на столбец значений правых частей уравнений:

для всех i = 1..n.

Тогда i-тый корень уравнения находят как

и

где ‑ алгебраическое дополнение a_ij.

Например, для матрицы размером 33 дополнение элемента a₁₂:

,

а определитель матрицы:

.

3.4. Итерационный метод решения системы линейных

алгебраических уравнений (метод Гаусса-Зейделя)

Приведенные выше методы решения СЛАУ носят конечный характер, что является их преимуществом. Однако во многих практических случаях системы уравнений предпочтительно решать итерационными методами, задав необходимую точность решения. Эти методы более удобны в их реализации с помощью вычислительных машин.

Рассмотрим систему из n уравнений c n неизвестными (3.1). По-прежнему полагают, что диагональные коэффициенты не равны нулю.

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

.

Итерационный процесс продолжается до тех пор, пока все не станут достаточно близки к. Критерием близости можно, например, задать условиеили, где определяют максимальное значение разности для всехi.

3.5. Решение системы нелинейных уравнений

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

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

Предполагается, что для системы уравнений вид итерационной формулы сохраняется:

, где ‑ вектор переменных,

‑ вектор значений функций,

‑ матрица Якоби,

‑обращенная матрица Якоби, которую можно получить описанным ранее методом Гаусса-Жордана. В свою очередь, произведение находят как произведение двух матрицС = AB у которых число строк первой равно числу столбцов второй: (n ‑ число строк матрицы A и число столбцов матрицы B).

ЗАДАНИЯ

1. Даны нормированные масс-спектры (фиктивные) четырех углеводородов.

	12	13	14	15	16	24	25	26	27	28	29	30
CH4	1	3	12	100	25
C2H6			2	80		3	4	13	37	24	100	70
C2H4	4	7	48	13		2	7	21	18	100
C2H2	5	23	57	2		40	45	100

Даны интенсивности четырех различных смесей этих четырех углеводородов:

Высота пика	Смесь 1	Смесь 2	Смесь 3	Смесь 4
Масса 14	119	169	326	60
Масса 15	195	284	501	113
Масса 26	134	168	433	21
Масса 28	124	148	124	100

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

2. Определить концентрации четырех компонентов смеси помощью закона Бера. Толщину поглощающего слоя принять равной 1,0 см.

Длина волны	Молярный коэффициент поглощения, л/(мольсм)
	n-ксилол	м-ксилол	о-ксилол	этил-бензол	общее поглощение
12,5	1,5020	0,0514	0	0,0408	0,10130
13,0	0,0261	1,1516	0	0,0820	0,09943
13,4	0,0342	0,0355	2,5320	0,2933	0,21940
14,3	0,0340	0,0684	0	0,3470	0,03396

3. Для определения концентрации компонентов газовой смеси можно использовать масс-спектроскопию, если наблюдаемая высота пика является линейной функцией парциального давления:

, где H_i – наблюдаемая высота пика, n – число компонентов, S_ij – чувствительность прибора к j – тому газу при определенном значении m/e, измеренная для чистого газа, P_j – парциальное давление j – того газа. Рассчитать мольную долю каждого газа в следующих смесях.

а) Смесь углеводородов:

m/e	Чувствительность				Общая высота пика
	этилцикло- пентан	цикло- гексан	цикло- гептан	метил- цикло- гексан
69	121,00	22,40	27,10	23,00	87,60
83	9,35	4,61	20,70	100,00	58,80
84	1,38	74,90	1,30	6,57	47,20
98	20,20	0,00	32,80	43,80	100,00

б) Смесь водорода и углеводородов:

m/e	Чувствительность							Общая высота пика
	H2	CH4	C2H4	C2H6	C3H6	C3H8	n-C5H12
2	16,87	0,165	0,202	0,317	0,234	0,182	0,110	17,1
16	0	27,70	0,962	0,062	0,730	0,131	0,120	65,1
26	0	0	22,35	13,05	4,420	6,001	3,043	186,0
30	0	0	0	11,28	0	1,110	0,371	82,7
40	0	0	0	0	9,850	1,684	2,108	84,2
44	0	0	0	0	0,299	15,98	2,107	63,7
72	0	0	0	0	0	0	4,670	119,7

Полное давление в системе равно 3,9910^‑2 мм.рт.ст.

4. Даны четыре емкости с растворами кислоты различной концентрации. Если смешать растворы в определенном соотношении, то получится кислота заданной процентной концентрации:

Концентрация конечного раствора, %	Доля кислоты в общей смеси
	1, %	2, %	3, %	4, %
25	1	1	1	1
20	4	3	2	1
25	4	1	1	4
22	4	1	4	1

Найдите концентрации кислот в каждом сосуде: 1, 2, 3 и 4.

5. Пусть в некоторой системе одновременно протекают две химические реакции. Они характеризуются выходами U₁ и U₂. Зависимость выхода, выраженного в %, в достаточно узкой области описывается двумя линейными относительно четырех параметров уравнениями:

U₁ = 10 + 1∙C₁ ‑ 2∙C₂ + 2.0∙pH + 0.5∙(T -300);

U₂ = 15 ‑ 1∙C₁ ‑ 1∙C₂ + 1.2∙pH + 0.9∙(T -300);

T – температура в градусах Кельвина; C₁ – концентрация вещества 1, моль/л; C₂ ‑ концентрация вещества 2, моль/л.

Найти pH 4 и температуру в градусах Кельвина при C₁ = C₂ = 2 моль/л, если выходы в этих условиях составляют U₁ = 36% и U₂ = 51,8%.

6. Найти мольные доли компонентов раствора (H₂SO₄, HSO₄^-, SO₄^2-) для заданного pH и исходной концентрации кислоты, решив систему уравнений. Например, для С₀=0,1 и H=2.76.

Диссоциация кислоты по первой ступени H₂SO₄  H⁺ + HSO₄^-

описывается константой диссоциации: .

Диссоциация кислоты по второй ступени HSO₄^-  H⁺ + SO₄^2-

описывается константой диссоциации: .

Уравнение материального баланса: .

7. Зависимость температуры плавления бинарной металлической системы от ее состава описывается соотношениями.

Найти состав эвтектической смеси и температуру плавления эвтектического сплава, если температуры плавления чистых компонентов T₁ = 316.15 К, Т₂ = 235,45 К, а энтальпии их плавления H₁ = 17200 Дж , H₂ = 18800 Дж.

8. Дана следующая схема химических реакций: .

Скорость расходования и накопления вещества описывается формулой:

.

Скорость реакции измерена при четырех различных концентрациях и температурах (см. табл.). Рассчитать четыре параметра уравнения Аррениуса k₀₁, E₁, k₀₂, E₂.

[A]	[B]	T	w
1	1	300	1.96610‑3
1	2	310	-0.0975
2	1	320	0.2435
1	1	330	-0.06219