13.7

Приближенное интегрирование дифференциальных уравнений

[29/62%]
Показать
LaTeX
Задача 5342

Найти разложение в степенной ряд общего интеграла уравнения yxy=1y' - xy = 1.

?
Задача 5343

Найти разложение в степенной ряд общего интеграла уравнения (1+x)yny=0(1+x)y' - ny = 0.

?
Задача 5344

Найти разложение в степенной ряд общего интеграла уравнения y+y=0y'' + y = 0.

?
Задача 5345

Найти разложение в степенной ряд общего интеграла уравнения yxy=0y'' - xy = 0.

?
Задача 5346

Доказать, что соответствующий частный интеграл дифференциального уравнения

x(x1)y+((α+β+1)xγ)y+αβy=0 x(x-1)y'' + \left((\alpha +\beta +1)x-\gamma \right)y' + \alpha \beta y = 0

может быть представлен указанным степенным рядом:

y=1+αβ1γx+α(α+1)12β(β+1)γ(γ+1)x2+α(α+1)(α+2)123β(β+1)(β+2)γ(γ+1)(γ+2)x3+ y = 1 + \frac{\alpha \cdot \beta }{1\cdot \gamma }x + \frac{\alpha (\alpha +1)}{1\cdot 2}\cdot \frac{\beta (\beta +1)}{\gamma (\gamma +1)}x^{2} + \frac{\alpha (\alpha +1)(\alpha +2)}{1\cdot 2\cdot 3}\cdot \frac{\beta (\beta +1)(\beta +2)}{\gamma (\gamma +1)(\gamma +2)}x^{3} + \ldots

(гипергеометрический ряд Гаусса).

?
Задача 5347

Доказать, что соответствующий частный интеграл дифференциального уравнения

x2y+xy+(x2n2)y=0 x^{2}y'' + xy' + (x^{2}-n^{2})y = 0

может быть представлен указанным степенным рядом:

y=ν=0(1)ν(x2)n+2νν!Γ(n+ν+1) y = \sum _{\nu =0}^{\infty } \frac{(-1)^{\nu }\left(\frac{x}{2}\right)^{n+2\nu }}{\nu !\, \Gamma (n+\nu +1)}

(функция Бесселя Jn(x)J_{n}(x)).

?
Задача 5348

Доказать, что соответствующий частный интеграл дифференциального уравнения

x(x41)u(3x41)u+x3u=0 x(x^{4}-1)u'' - (3x^{4}-1)u' + x^{3}u = 0

может быть представлен указанным степенным рядом:

u=1+(14)2x4+(1548)2x8+(1594812)2x12+;x<1. u = 1 + \left(\frac{1}{4}\right)^{2}x^{4} + \left(\frac{1\cdot 5}{4\cdot 8}\right)^{2}x^{8} + \left(\frac{1\cdot 5\cdot 9}{4\cdot 8\cdot 12}\right)^{2}x^{12} + \ldots ; \qquad \left|x\right| < 1.
?
Задача 5349

Доказать, что соответствующий частный интеграл дифференциального уравнения

x2(x41)v+3x(3x41)v+(19x41)v+8x3v=0 x^{2}(x^{4}-1)v''' + 3x(3x^{4}-1)v'' + (19x^{4}-1)v' + 8x^{3}v = 0

может быть представлен указанным степенным рядом:

v=1+(12)3x4+(1324)3x8+(135246)3x12+;x<1. v = 1 + \left(\frac{1}{2}\right)^{3}x^{4} + \left(\frac{1\cdot 3}{2\cdot 4}\right)^{3}x^{8} + \left(\frac{1\cdot 3\cdot 5}{2\cdot 4\cdot 6}\right)^{3}x^{12} + \ldots ; \qquad \left|x\right| < 1.
?
Задача 5350

Доказать равенство Гаусса:

(1+(14)2x4+(1548)2x8+(1594812)2x12+)2=1+(12)3x4+(1324)3x8+(135246)3x12+;x<1. \left(1 + \left(\frac{1}{4}\right)^{2}x^{4} + \left(\frac{1\cdot 5}{4\cdot 8}\right)^{2}x^{8} + \left(\frac{1\cdot 5\cdot 9}{4\cdot 8\cdot 12}\right)^{2}x^{12} + \ldots \right)^{2} = 1 + \left(\frac{1}{2}\right)^{3}x^{4} + \left(\frac{1\cdot 3}{2\cdot 4}\right)^{3}x^{8} + \left(\frac{1\cdot 3\cdot 5}{2\cdot 4\cdot 6}\right)^{3}x^{12} + \ldots ; \qquad \left|x\right|<1.
?
Задача 5351

Доказать, что, применяя способ Пикара к уравнению y=yy' = y, при a=0a=0, b=1b=1, получим:

y=1+x+x22!++xn1(n1)!+xnn!+ y = 1 + x + \frac{x^{2}}{2!} + \ldots + \frac{x^{n-1}}{(n-1)!} + \frac{x^{n}}{n!} + \ldots
?
Примечание.
?

Способ Пикара нахождения интеграла уравнения y=φ(x,y)y' = \varphi (x,y), который при x=ax=a обращается в bb, состоит в том, что составляется последовательно ряд функций y0,y1,y2,y_{0}, y_{1}, y_{2}, \ldots по формуле yn=b+axφ(x,yn1)dxy_{n} = b + \displaystyle \int_{a}^{x} \varphi (x,y_{n-1})\, dx; y0=by_{0}=b. Если ynyy_{n} \to y при nn \to \infty, то yy и есть искомый интеграл.

Задача 5352

Доказать, что в предыдущем примере (задача 5351) ynyy_{n} \to y, если nn \to \infty при произвольном выборе функции y0y_{0}.

?
Задача 5353

Доказать, что, заменив уравнение

y+p(x)y+q(x)y=0 y'' + p(x)y' + q(x)y = 0

системой уравнений

y=u;u=p(x)uq(x)y y' = u; \qquad u' = -p(x)u - q(x)y

и применяя способ Пикара к решению этих уравнений, получим сходящийся процесс, если радиус-вектор из начала в точку xx не проходит на комплексной плоскости переменной xx ни через одну особую точку функций p(x)p(x) и q(x)q(x).

?
Задача 5354

Решить уравнение yy=axyy' - y = axy по методу последовательных приближений, полагая y0=0y_{0}=0 и вычисляя величины yny_{n} по формуле ynyn=axyn1y_{n}' - y_{n} = axy_{n-1} (при x=0x=0 должно быть y=Cy=C).

?
Задача 5355

Таким же путем (см. задачу 5354) решить уравнение y+a2y=σxyy'' + a^{2}y = \sigma xy, где σ\sigma — малый параметр. При этом y=0y=0, y=ay'=a при x=0x=0. Полезно применить теорему: если y+a2y=f(x)y''+a^{2}y=f(x), то

y=C1sinax+C2cosax+1a0xf(t)sina(xt)dt. y = C_{1}\sin ax + C_{2}\cos ax + \frac{1}{a}\int _{0}^{x} f(t)\sin a(x-t)\, dt.
?
Задача 5356

Пользуясь подстановкой x=n+1nx = n + \dfrac {1}{\sqrt{n}}, показать, что интеграл уравнения y+xy=0y'' + xy = 0 в интервале (n,π+n)(n, \pi +n), где nn велико, имеет приблизительно n\sqrt{n} корней.

?
Задача 5357

Доказать, что интеграл уравнения Бесселя

x2y+xy+(x2n2)y=0 x^{2}y'' + xy' + (x^{2}-n^{2})y = 0

в весьма удаленном от начала интервале (a,a+b)(a, a+b) имеет приблизительно bπ\dfrac {b}{\pi } корней.

?
Задача 5358

Интеграл уравнения y+y=σy2y'' + y = \sigma y'^{2} имеет начальные условия: x=0x=0, y=1y=1, y=0y'=0. Найти три первых члена его разложения по степеням малого параметра σ\sigma.

?
Задача 5359

Многочлены Чебышева--Эрмита: Hn(x)=ex22n!dnex22dxnH_{n}(x) = \dfrac {e^{\frac{x^{2}}{2}}}{n!}\dfrac {d^{n}e^{-\frac{x^{2}}{2}}}{dx^{n}} удовлетворяют уравнению

Hn(x)xHn(x)+nHn(x)=0. H_{n}''(x) - xH_{n}'(x) + nH_{n}(x) = 0.

Подстановкой Hn(x)=f(x)yH_{n}(x) = f(x)\cdot y преобразовать это уравнение к виду

y+(x24+n+12)y=0. y'' + \left(-\frac{x^{2}}{4} + n + \frac{1}{2}\right)y = 0.
?
Задача 5360

Найти асимптотическое выражение при nn \to \infty для многочленов Чебышева--Эрмита Hn(x)=ex22n!dnex22dxnH_{n}(x) = \dfrac {e^{\frac{x^{2}}{2}}}{n!}\dfrac {d^{n}e^{-\frac{x^{2}}{2}}}{dx^{n}} (см. задачу 5359).

?
Задача 5361

Найти при 0<x<π0<x<\pi интеграл уравнения

y+4y=πxx2, y'' + 4y = \pi x - x^{2},

разложив его в ряд вида ansinnx\sum a_{n}\sin nx. Краевые условия: y=0y=0 при x=0x=0 и x=πx=\pi.

?
Задача 5362

Найти интеграл уравнения yIV+y=πxy^{\mathrm{IV}} + y = \pi - \left|x\right| при π<x<π-\pi <x<\pi и краевых условиях y=y=0y' = y''' = 0 для x=±πx = \pm \pi.

?
Задача 5363

Составить таблицу значений yy через 0,050{,}05, если y+2xy=0y'+2xy=0 и y=1y=1 при x=0x=0. При этом значения y2,y1,y0,y1,y2y_{-2}, y_{-1}, y_{0}, y_{1}, y_{2} вычислить с помощью разложения в ряд, в формуле Адамса удерживать члены до третьих разностей, а вычисления вести до шестого знака после запятой.

?
Примечание.
?

Численное интегрирование дифференциального уравнения y=f(x,y)y'=f(x,y) состоит в том, что по данному начальному условию y=y0y=y_{0} при x=x0x=x_{0} составляют таблицу значений yy для значений xx, равных числам x0,x0+h,x0+2h,x_{0}, x_{0}+h, x_{0}+2h, \ldots, где hh — выбранное число, обычно довольно малое.

Формула Адамса--Штермера: Δyn=un+12Δun1+512Δ2un2+38Δ3un3+251720Δ4un4+13Δ5un5+\Delta y_{n} = u_{n} + \dfrac {1}{2}\Delta u_{n-1} + \dfrac {5}{12}\Delta^{2}u_{n-2} + \dfrac {3}{8}\Delta^{3}u_{n-3} + \dfrac {251}{720}\Delta^{4}u_{n-4} + \dfrac {1}{3}\Delta^{5}u_{n-5} + \ldots, где u=hf(x,y)u = hf(x,y); эта формула дает возможность вычислить yn+1y_{n+1}.

Задача 5364

Найти по методу Адамса--Штермера при x=2x=2 интеграл уравнения y=x+yy' = x+y, обращающийся в единицу при x=1x=1; при этом пользоваться формулой для Δy\Delta y, взятой до четвертых разностей включительно, и вести вычисления до пятого знака.

?
Задача 5365

Найти таблицу значений интеграла уравнения

y=x+yдляx=0; 0,1; 0,2;  1, y' = \sqrt{x} + \sqrt{y} \quad \text{для} \quad x = 0;\ 0{,}1;\ 0{,}2;\ \ldots \ 1,

если при x=0x=0, y=0y=0.

?
Примечание.
?

Разложение yy в начале имеет вид y=23x32+4723x74+17x2+14923x94+y = \dfrac {2}{3}x^{\frac32} + \dfrac {4}{7}\sqrt{\dfrac {2}{3}}x^{\frac74} + \dfrac {1}{7}x^{2} + \dfrac {1}{49}\sqrt{\dfrac {2}{3}}x^{\frac94} + \ldots Выгодно ввести новую переменную, положив x=ξ4x = \xi^{4}.

Задача 5366

Составить таблицу значений функций sinx\sin x и cosx\cos x через 0,10{,}1 от 00 до 11, пользуясь уравнениями: (sinx)=cosx(\sin x)' = \cos x; (cosx)=sinx(\cos x)' = -\sin x и применив метод Адамса--Штермера.

?
Задача 5367

Составить таблицу значений функций Бесселя J0(x)J_{0}(x) и J0(x)J_{0}'(x) для значений xx через одну десятую от нуля до единицы, при этом уравнение для J0(x)J_{0}(x), имеющее вид

xJ0(x)+J0(x)+xJ0(x)=0, xJ_{0}''(x) + J_{0}'(x) + xJ_{0}(x) = 0,

заменить системой: xv+v+xu=0xv' + v + xu = 0; u=vu' = v.

?
Примечание.
?

Для уравнений второго порядка вида y=f(x,y,y)y'' = f(x,y,y') можно пользоваться формулами:

Δ2yn1=h2fn+112h2Δ2fn2+112h2Δ3fn3+19240h2Δ4fn4+340h2Δ5fn5,Δyn=hyn+h22fn+h26Δfn1+h28Δ2fn2+19h2180Δ3fn3+3h232Δ4fn4,Δyn=hfn+h2Δfn1+512hΔ2fn2+38hΔ3fn3+251720hΔ4fn4. \begin{align} \Delta ^{2}y_{n-1} & = h^{2}f_{n} + \frac{1}{12}h^{2}\Delta ^{2}f_{n-2} + \frac{1}{12}h^{2}\Delta ^{3}f_{n-3} + \frac{19}{240}h^{2}\Delta ^{4}f_{n-4} + \frac{3}{40}h^{2}\Delta ^{5}f_{n-5}, \\ \Delta y_{n} & = hy_{n}' + \frac{h^{2}}{2}f_{n} + \frac{h^{2}}{6}\Delta f_{n-1} + \frac{h^{2}}{8}\Delta ^{2}f_{n-2} + \frac{19h^{2}}{180}\Delta ^{3}f_{n-3} + \frac{3h^{2}}{32}\Delta ^{4}f_{n-4}, \\ \Delta y_{n}' & = hf_{n} + \frac{h}{2}\Delta f_{n-1} + \frac{5}{12}h\Delta ^{2}f_{n-2} + \frac{3}{8}h\Delta ^{3}f_{n-3} + \frac{251}{720}h\Delta ^{4}f_{n-4}. \end{align}
Задача 5368

Составить таблицу значений через 0,10{,}1 в интервале (0,1)(0,1) функции y=chxy = \operatorname {ch}x, удовлетворяющей уравнению y=yy'' = y и начальным условиям: при x=1x=1 y=1y=1, y=0y'=0.

?
Задача 5369

Составить таблицу значений функции J0(x)J_{0}(x) для значений xx, идущих через 0,10{,}1, зная, что J0(x)J_{0}(x) удовлетворяет дифференциальному уравнению xJ0(x)+J0(x)+xJ0(x)=0xJ_{0}''(x) + J_{0}'(x) + xJ_{0}(x) = 0, и начальные ее значения при x=0x=0: J0(0)=1J_{0}(0)=1, J0(0)=0J_{0}'(0)=0.

?
Задача 5370

Зная величины: cos1=0,54030\cos 1 = 0{,}54030; cos1,1=0,45360\cos 1{,}1 = 0{,}45360; cos1,2=0,36236\cos 1{,}2 = 0{,}36236; cos1,3=0,26750\cos 1{,}3 = 0{,}26750; cos1,4=0,16997\cos 1{,}4 = 0{,}16997, с помощью уравнения y=yy'' = -y найти значения cosx\cos x при x=1,5x=1{,}5; 1,61{,}6; 1,71{,}7; 1,81{,}8; 1,91{,}9; 2,02{,}0.

?