ВУЗ:
Составители:
Встроенная функция polyder предназначена для вычисления производных от полиномов. Вызов
polyder с аргументом – вектором, соответствующим полиному, приводит к вычислению вектора коэф-
фициентов производной полинома:
>> p = [1 0 1 0 0 1];
>> p1 = polyder(p)
p1 =
5 0 3 0 0
6.4 Интегрирование функций
Для вычисления интеграла используется quad, задавая первым аргументом имя файл-функции, от
которой вычисляется интеграл, а вторым и третьим – нижний и верхний предел интегрирования. Данная
функция вычисляет приближенное значение интеграла с точностью 10е-3.
>> I = quad('myfunction', a, b)
П р и м е р: Вычислить летучесть аммиака f как функцию давления при заданной температуре.
Наиболее точные значения коэффициента летучести получают из экспериментальных значений коэф-
фициента сжимаемости, определяемого из P-V-T измерений. Зависимость сжимаемости аммиака от дав-
ления определена экспериментально (табл. 9).
Для вычисления летучести как функции давления при заданной температуре используется зависи-
мость:
∫
−
=
P
dP
P
Z
P
f
0
1
ln
,
где Z – функция сжимаемости газа; P – давление; f – летучесть аммиака.
Вычислить летучесть и объем газа при заданных условиях.
9 Зависимость сжимаемости аммиака от давления
№
опыта
Давле-
ние
Функция
сжимае-
мости газа
Z
№ опы-
та
Давле-
ние
Функция
сжимае-
мости газа
Z
1 1 0,998 9 100 0,801
2 10 0,981 10 200 0,551
3 20 0,961 11 300 0,462
4 30 0,942 12 400 0,495
5 40 0,922 13 500 0,557
6 50 0,902 14 600 0,621
7 60 0,882 15 800 0,755
8 80 0,841
В функции ammiak1 записывается подынтегральная функция
function f = ammiak1(p, z)
% расчет зависимости летучести аммиака
% z – независимый параметр (летучесть)
f = (z-1)./p;
Процедура ammiak рассчитывает летучесть и объем газа, результаты расчета выводятся на графики
в одно окно.
Страницы
- « первая
- ‹ предыдущая
- …
- 45
- 46
- 47
- 48
- 49
- …
- следующая ›
- последняя »
