Опубликован: 03.12.2012 | Уровень: для всех | Доступ: платный
Лекция 9:

Уравнения. Системы уравнений

8.5. Нахождение решений рекуррентных уравнений

Рекуррентными называются функции, значение которых в одной точке зависит от значения этой же функции в другой точке. Отношения между значениями функции в разных точках задаются рекуррентными уравнениями. Решает рекуррентные уравнения функция RSolve: в качестве результата она выдаёт зависимость значения функции от номера итерации. Задаётся функция в виде RSolve[eq,f[n],n], т.е. имеет, по крайней мере, три аргумента: первый — уравнение (систему уравнений) eq с начальными условиями, второй — функция (или список функций) f[n], которая зависит от номера итерации, счётчик которых n указываемой в третьем аргументе. Решение рекуррентного уравнения приведено в примерах In[1] и In[2] на рис. 8.21. Как мы видим, оба расчёта имеют одинаковый результат (ср. Out[1] и Out[2]) при том, что номера итераций указаны различным образом.

Функцию RSolve можно рассматривать как дискретный аналог DSolve. Именно поэтому специальные функции, генерируемые DSolve при решении дифференциальных уравнений, встречаются и при решении рекуррентных уравнений. В примере In[3] на рис. 8.21 при решении рекуррентного уравнения мы получили зависимость, содержащую функцию Бесселя.

О решении рекуррентных уравнений см. книгу В. П. Дьяконова [2, с. 263–264].

Решение рекуррентных уравнений при помощи функции RSolve

Рис. 8.21. Решение рекуррентных уравнений при помощи функции RSolve

Ключевые термины

Рекуррентными называются функции, значение которых в одной точке зависит от значения этой же функции в другой точке.

Трансцендентное уравнение — уравнение, не являющееся алгебраическим.

Краткие итоги

В данной лекции мы познакомились с основными встроенными функциями Mathematica для нахождения символьных и численных решений алгебраических, трансцендентных, рекуррентных и дифференциальных уравнений, а также систем уравнений. Мы познакомились с принципами работы некоторых встроенных функций решения уравнений. Также мы научились сочетать возможности Mathematica символьных вычислений, численных расчётов и графических средств для нахождения решений уравнений и систем уравнений.

Вопросы

  1. Какая функция Mathematica используется для нахождения символьных решений алгебраических уравнений? В каком виде задаются её аргументы? В каком виде она возвращает результат вычислений?
  2. Каким образом Mathematica представляет корни полиномиального уравнения в случае, если она не может представить полученное при помощи функции Solve решение в аналитическом виде?
  3. Каким образом Mathematica может отреагировать на попытку решить при помощи функции Solve трансцендентное уравнение?
  4. Какое действие выполняет функция FindInstance, применённая к алгебраическому уравнению? В каком виде задаются её аргументы?
  5. Для каких целей служит функция Eliminate?
  6. В каких случаях при работе с алгебраическими уравнениями функции Solve следует предпочесть функцию Reduce? В каком виде возвращается результат выполнения функции Reduce?
  7. Для чего используется функция SolveAlways? Какие обязательные аргументы она содержат?
  8. В каком виде должны быть заданы аргументы функции решения систем линейных уравнений LinearSolve?
  9. Какая функция Mathematica используется для нахождения численных решений алгебраических уравнений? В каком виде она возвращает результат вычислений?
  10. В каких случаях при решении уравнений используется функция FindRoot? Каким образом она ищет корни уравнений? Для чего используется опция MaxIterations?
  11. Для решения каких типов ДУ и систем ДУ используется функция DSolve? Какие обязательные аргументы она содержат?
  12. В чём заключается разница в результатах вычислений при задании второго аргумента функции DSolve как f и f[x]?
  13. В чём состоят отличия при задании аргументов функций аналитического DSolve и численного NDSolve решения ДУ?
  14. В каком виде возвращает результат решения функция NDSolve?
  15. Для каких целей используется функция RSolve? В каком виде задаются её обязательные аргументы?

Упражнения

  1. Решите средствами Mathematica следующие уравнения и системы уравнений. Для каждого уравнения или системы уравнений используйте по два наиболее подходящих способа решения, если это возможно. Сравните результаты решения каждым из способов друг с другом.

    x^2+3x-7=0

    (2x+5)^2=(3x-1)^2

    x^3-8x^2+4x+48=0

    x^5-2bx+4=0 (относительно x)

    x^2+y^2=3 (относительно x)

    6x-10/(x+1)=9

    ((x-5)/(x+2))^{1/2}+((x-4)/(x+3))^{1/2}=(7/(x+2))((x+2)/(x+3))^{1/2}

    (x+x^{1/2})^{1/2}-(x-x^{1/2})^{1/2}=1.5(x/(x+x^{1/2}))^{1/2}

    cos(2t)=4t

    sin(2x)=a (относительно x)

    14-3x=-7y, 7y^2-28=46x-7x^2

    x^{1/2}+y^{1/3}=4, (x+2y)^{1/2}=4

    ax+2y=0, 2x+(1-2a)y=1 (относительно x и y)

    x^2+2y^2=3, xy=4

    x^2+y^2+z^2=1, 2x+3y-4z=5,xyz=6

    (x^2+y^2)^{1/2}-(x^2-y^2)^{1/2}=y, x^4+y^4=144a^4 (относительно x и y)

    ax^2+by^2=c, x^2=2y

    exp(x-1)=y, y=x^{1/2}

  2. Из приведённых ниже систем уравнений исключите указанные переменные.

    x+2y+3z=4, 5x+6y+7z=8 (исключить переменную z)

    x+2y+3z=4, 5x+6y+7z=8 (переменную y)

    x^2+y^2+z^2=1, x+y+z=2, x-y^3=z (переменную z)

    x^2+y^2+z^2=1, x+y+z=2, x-y^3=z (переменные y и z)

    x^2+y^2+z^2=1, x+y+z=2, x-y^3=z (переменные x, y и z)

    b^3+(a^2-4)b^2+(2a-3)b+4, 0=(1-a^2)b^4+ab+7a^2 (переменную b)

    sin(x+y)=x, cos(x-y)=y (переменную y)

    z=x^5+y^5,a=x^2+y^2,b=x+y (переменные x и y)

    exp(x-1)=y, y=x^{1/2} (переменную y)

  3. Решите средствами Mathematica следующие ДУ и системы ДУ. Проиллюстрируйте полученные решения графически, если это возможно.

    y'(x)+y(x)=x cos(x)

    y'(x)+y(x)=x cos(x), y(0)=0

    y'(x)+x (y(x)-sin(x))=0, y(0)=x

    y''(x)-2y(x)-y'(x)/2=0

    y''(x)-2y(x)-y'(x)/2=0, y(0)=1, y(3)=-11

    y''(x)+2y'(x)+3y(x)=4, y(5)=6, y(7)=8

    3x^3 y'''(x)+10x^2 y''(x)+(x-x^2)y'(x)+(2-4x)y(x)=0

    y'(x)-(y(x)/x)^{1/3}=y(x)/x

    3df(u,v)/du+5df(u,v)/dv=u

    df(u,v)/du+4df(u,v)/dv=8, f(u,0)=u

    x''(t)=x(t)y(t), y'(t)=2-x(t), x(0)=x(2)=4, y(0)=1

    df(u,v)/du+d^3f(u,v)/dv^3+6f(u,v) df(u,v)/dv=0

    y'(t)=t^2y(t), x'(t)=5x(t)

    y'(t)=cos(\pi/4)y(t)-sin(\pi/4)x(t)+1/2, x'(t)=sin(\pi/4)y(t)+cos(\pi/4)x(t), y(0)=1/2, x(0)=1/2

    y'(t)-x(t)=3cos(t), y(t)-x(t)=3/2, y(\pi)=1/2

    y'(t)=x(t)^{-1/2}, x'(t)=y(t)

    x'(t)=-2(x(t)-y(t)), y'(t)=50x(t)-y(t)-x(t)z(t), z'(t)=x(t)y(t)-z(t), x(0)=z(0)=0, y(0)=1

    d^2f(u,v)/du^2=d^2f(u,v)/dv^2, f(0,v)=exp(-v^2), df(0,v)/du=0, f(u,-10)=f(u,10)

Наталья Егорова
Наталья Егорова
Россия, Новосибирск, НГПУ, 1994
Максим Радунцев
Максим Радунцев
Россия