Расчет и оптимизация процессов и аппаратов химических и пищевых производств в среде MatLab. Дворецкий Д.С - 43 стр.

UptoLike

го процесса
‘final’ Выводится информация
только о завершении вычис-
лительного процесса
‘Max-
FunEvals
Целое число Максимальное количество
вызовов исследуемой функ-
ции
Рис. 6.1 Решение системы дифференциальных уравнений
6.2 Решение уравнений
Корни произвольного уравнения по заданному приближению позволяет найти встроенная функ-
ция fzero, которая вызывается следующим образом:
>> x = fzero(‘myfunction’, x0)
где ‘myfunction’ имя файл-функции, вычисляющей левую часть уравнения; x0 начальное приближе-
ние к корню либо вектор, содержащий концы отрезка, на котором ищутся корни; x найденное при-
ближенное значение корня. Данная функция находит решение с точностью до шестнадцатого знака по-
сле запятой.
Функция, содержащая выражение для расчета pH, сохранена под именем concentr, а входным аргу-
ментом ее является значение концентрации ионов водорода.
function f = concentr(x)
% определение pH растворов слабых кислот
kw = 1.85e-5;
ka = 1.0e-14;
ha = 0.1;
f = (x^2-kw)/ka+x-(kw/x)-ha;
В программе aim производится построение графика функции concentr для определения начального
приближения, а затем решение нелинейного уравнения и определение pH.
% программа aim определения pH раствора
% построение графика функции для определения
% начального приближения