ВУЗ:
Составители:
REZ2 =
1.0e+004 *
5.0000
3.0000
8.0000
6.3 Вычисление корней полинома и производных от полинома
Полином задается вектором его коэффициентов, число элементов вектора, т.е. число коэффициен-
тов полинома всегда на единицу больше его степени, нулевые коэффициенты должны содержаться в
векторе.
Вычисление значения полинома от некоторого аргумента:
>> p = [1 0 3.2 -5.2 0 0.5 1 -3];
>> polyval(p, 3)
ans =
2.5479e+003
Нахождение сразу всех корней полинома осуществляется при помощи функции roots, которая воз-
вращает вектор корней полинома, в том числе и комплексные.
>> r = roots(p)
r =
-0.5668 + 2.0698i
-0.5668 – 2.0698i
-0.6305 + 0.5534i
-0.6305 – 0.5534i
1.2149
0.5898 + 0.6435i
0.5898 – 0.6435i
Встроенная функция 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).
Страницы
- « первая
- ‹ предыдущая
- …
- 32
- 33
- 34
- 35
- 36
- …
- следующая ›
- последняя »
