Материал: Matlab. Практический подход. Самоучитель

Глава 9. Символьные вычисления

Рис. 9.3. График функции, заданной на основе символьного выражения

Рис. 9.4. Замены в символьном выражении с несколькими переменными

351

Самоучитель Matlab

С вычислением значения такой своеобразной функции дела обстоят несколько хуже. Чтобы вычислить значение функции в точке, используем функцию subs(). Первым аргументом функции указывается переменная, определенная на основе символьного выражения. Второй аргумент – значение для символьной переменной. Это значение, кстати, тоже может быть символьным. В частности, командой subs(f,pi/2) вычисляется значение для числового аргумента, а командой subs(f,y) вычисляется выражение с заменой символьной переменной x на символьную переменную y.

Если в выражение входит несколько символьных переменных, то функцию subs() используют в немного ином формате. Первым аргументом, как и ранее, указывается переменная символьного выражения, а затем, через запятую, заменяемая в выражении символьная переменная и ее значение. На рис. 9.4 представлено продолжение предыдущего документа.

Команды здесь такие (жирным шрифтом выделен ввод пользователя):

>>z=cos(x*y)+x/y

z =

cos(x*y) + x/y

>>subs(z,x,1) ans =

cos(y) + 1/y

>>subs(z,y,x) ans =

cos(x^2) + 1

>>subs(subs(z,x,pi),y,0.5) ans =

6.2832

Командой z=cos(x*y)+x/y задается символьное выражение, в которое входят две символьные переменные x и y.

На заметку

Переменные x и y в качестве символьных объявлены ранее, в предыдущей части документа – см. рис. 9.3.

Результатом выполнения команды subs(z,x,1) является выражение, получающееся заменой в выражении z переменной x на 1. При этом выражение z не меняется. В результате выполнения команды subs(z,y,x) получаем выражение, которое получается из z заменой переменной y на x. Причем результат автоматически упрощается (насколько это возможно). Наконец, чтобы выполнить замену сразу нескольких символьных переменных, используем конструкцию из вложенных функций subs(), как в команде subs(subs(z,x,pi),y,0.5).

352

Глава 9. Символьные вычисления

Удобно то, что непосредственно в командном окне можно выполнять символьные операции и с числами. Для этого полезной станет функция sym(). Ее аргументом указывается (заключенное в одинарные кавычки) выражение, которое следует интерпретировать и обрабатывать как символьное. Таким символьным выражением может быть и число. Обратимся к примеру в документе, показанном на рис. 9.5.

Рис. 9.5. Символьные вычисления

Результатом команды sym('3/7')+sym('6/11') будет сумма дробей, причем вычисленная "точно", без округления. Результат представляется в виде правильной дроби. Далее, командой a=1+sqrt(sym('2')) опреде-

ляется переменная a с точным (символьным) значением 1 + 2 . Операции с этой переменной выполняются так, что их результат также вычисляется в символьном виде. Например, результат выражения a^2 символьный. Правда, он не очень компактный, но вызов команды simplify(a^2) упрощает соответствующее выражение.

На заметку

Для преобразования символьного числового выражения в формат числа с плавающей точкой используется функция vpa(). Ее аргументом указывают вычисляемое (преобразуемое) выражение. Второй необязательный аргумент определяет количество значащих цифр после десятичной запятой.

353

Самоучитель Matlab

Основные символьные операции

- Интересное решение. - Вы думаете?

- А что? Демократично, удобно.

- Честно говоря, мне виделось нечто более изящное.

К/ф "Чародеи"

В этом разделе кратко остановимся на некоторых задачах, которые часто приходится (или желательно) решать в символьном виде. К ним, безусловно, следует отнести дифференцирование, интегрирование, разложение функции в ряд Тейлора, а также решение дифференциальных и алгебраических уравнений. Для выполнения всех этих операций (а также многих других) в Matlab есть специальные встроенные функции, с которыми нам, собственно, и предстоит познакомиться. Начнем с дифференцирования.

Для вычисления производных (в том числе и частных) от символьных выражений используется функция diff(). У нее есть несколько форматов вызова. Самый простой случай – если необходимо вычислить производную от символьного выражения, в котором всего одна символьная переменная. В этом случае выражение указывается аргументом функции diff(). Этот и другие примеры вызова функции diff() представлены в документе на рис. 9.6.

Рис. 9.6. Вычисление производных

354

Глава 9. Символьные вычисления

Командой syms x y объявляются две символьные переменные x и y. Затем на основе этих переменных объявляем символьное выражение f=sin(x)/x. Для вычисления производной от этого выражения используем команду diff(f), в которой единственным аргументом функции diff() передано дифференцируемое выражение. Производная автоматически вычисляется по символьной переменной. Если переменных в выражении несколько, производная вычисляется по той, которая в алфавите "ближе" к букве x. Например, определим выражение F=sin(x)*y^2 и продифференцируем его с помощью команды diff(F). В данном случае производная вычисляется по переменной x. Чтобы вычислить производную по переменной y, используем команду diff(F,y). Здесь вторым аргументом явно указана переменная, по которой вычисляется производная.

Если вторым аргументом функции diff() указать целое число, то оно определит порядок производной от выражения, переданного первым аргументом функции. Обратимся к документу на рис. 9.7.

Рис. 9.7. Производная второго порядка

Здесь командой f=sym('(x^3-5*x^2+3*x)*sin(x)') определяется символьное выражение, в которое входит всего лишь одна символьная переменная x.

На заметку

Обратите внимание, как определяется символьное выражение: в одинарных кавычках оно указывается аргументом функции sym().

Вторую производную от созданного символьного выражения вычисляем с помощью команды diff(f,2). Можно также вызвать функцию diff() со следующим способом передачи аргументов: дифференцируемое выражение, переменная дифференцирования и порядок дифференцирования. Например, команда diff(f,x,2) означает вычисление второй производной от выражения f по переменной x.

355