WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     | 1 |   ...   | 9 | 10 || 12 | 13 |

(Rm()) - m2Rm() = 0, (5.40) Rm(a) = am, Rm(0) <.

Условие при = 0 следует из естественного требования ограниченности функции u(, ). Уравнение в (5.40) является уравнением Эйлера, его частные решения можно искать в виде Rm r. Общие решения имеют вид R0() = A0 + B0 ln, Rm() = Amr|m| + Bmr-|m|, m = 0. (5.41) Из граничных условий определим константы:

am Rm(0) < Bm = 0, m; Rm(a) = am Am =. (5.42) a|m| Решение, записанное в виде ряда Фурье, имеет вид + |m| eim u(, ) = am. (5.43) a m=Ряд (5.43) можно просуммировать. Используя выражения (5.39) для коэффициентов Фурье функции f(), получим + + |m| eim 2 e-im 1 |m| f() d = f() eim(-) d.

a 2 a 2 0 m=- m=Преобразуем ряд в правой части равенства:

+ + + |m| m m eim(-) = eim(-) + e-im(-) = a a a m=- m=0 m=+ m = 2 Re eim(-) - 1.

a m= Введя обозначение q = ei(-), просуммируем ряд и выделим реальa ную часть:

+ 2 1+q 1 - |q|Re 2 qm - 1 = Re - 1 = Re =.

1-q 1-q 1-(q +q)+|q|m=Возвращаясь к исходным обозначениям, получим решение задачи (5.36) в виде интеграла Пуассона:

1 a2 - u(, ) = f() d. (5.44) 2 a2 - 2a cos( - ) + Отметим, что при решении конкретных задач бывает проще использовать формулу (5.43), так как вычислить коэффициенты Фурье am функции f() по формуле (5.39) зачастую проще, чем вычислить интеграл Пуассона.

5.13. Найти распределение потенциала электростатического поля u(x, y) внутри прямоугольника G = {0 < x < a, 0 < y < b}, если потенциал стороны этого прямоугольника, лежащей вдоль оси 0y, равен v0, а три другие стороны прямоугольника заземлены. Предполагается, что внутри G нет электрических зарядов.

5.14. Найти стационарное распределение температуры u(x, y) в прямоугольной однородной пластине G = {0 < x < a, 0 < y < b}, если ее стороны x = a и y = b покрыты тепловой изоляцией, две другие стороны (x = 0 и y = 0) поддерживаются при нулевой температуре, а в пластинке выделяется тепло с постоянной плотностью q.

5.15. Найти стационарное распределение температуры u(r, ) внутри бесконечного цилиндра радиуса r0, если на одной половине поверхности цилиндра (0 < < ) поддерживается температура -T0, а на другой половине ( < < 2) температура T0.

5.16. Найти решение краевой задачи для уравнения Пуассона в круге G = x2 + y2 < r0 :

u = -Axy, u|x +y2=r0 = 0.

5.17. Найти решение краевой задачи для уравнения Лапласа в кольце G = 1 < x2 + y2 < 4 (u1 и u2 константы):

u = 0, u|x +y2=1 = u1, u|x +y2=4 = u2.

5.18. Найти стационарное распределение температуры u(r) внутри цилиндра с радиусом основания r0 и высотой h, если к нижнему основанию подводится постоянный тепловой поток q, а боковая поверхность и верхнее основание поддерживаются при нулевой температуре.

5.19. Найти потенциал электростатического поля внутри цилиндра G = {x2 + y2 < a2, 0 < z < h}, если на его боковой поверхности поддерживается потенциал U0 z, U0 константа, а на торцах задано нулевое электрическое поле.

5.20. Найти потенциал электростатического поля внутри сферы радиуса a, если потенциал сферы имеет вид U0, 0 < < /2, f() = 0, /2 < <.

Ответы и указания к главе 5.1.

8 1 1) u(t, x) = e-(2k+1) t sin(2k + 1)x.

(2k + 1)k=2) u(t, x) = e4t - 1 + t cos 2x.

4 1 3) u(t, x) = e-t e-(2k+1) t - 1 cos(2k + 1)x.

(2k + 1)k=4) u(t, x) = et + sin t - cos t + e-3t cos 2x.

2 4 1 - e-(4m -1) t 5) u(t, x) = et - 1 + cos x - cos 2mx.

(4m2 - 1)m=5.2.

8 1) u(t, x) = - cos [(2k + 1)t] sin [(2k + 1)x].

3 (2k + 1)k=2) u(t, x) = sh2 t - t2 cos 2x.

x 2 1 1 3) u(t, x) = cos + cos k+ t - 1 cos k+ x.

2 2 (k + )k= 1 4) u(t, x) = - t cos 2t sin x + sin[2 3t] sin 3x.

2 5.3.

1) u(t, x) = x + t sin x + 1 - e-8t sin 3x.

2) u(t, x) = x2 + 2e9t + (2t - sin 2t) cos 3x.

1 3) u(t, x) = x + t2 + e5t - 1 cos x + 1 - e-3t cos 3x.

5 2 1 1 e[ 6-(2k+1)2] t-4) u(t, x) = tx2+x- - cos(2k+1)x.

2k-1 2k+3 6 - (2k + 1)k=5.4.

1) u(t, x) = 2tx + 2et - e-t - 3te-t cos x.

1 t sin µkt 2) u(t, x) = t (1-x) + e- 2 cos µkt + -2 sin kx, (k)3 µk k=где µk = (k)2 -.

3) u(t, x) = 3 + t + t2 x + 8 + 4t - 8et + 5tet sin x.

1 1 2 (-1)k 4) u(t, x) = tx + t3 sin x + sin µkt - t sin kx, 3 kµ2 µk k k= где µk = k2 - 1.

5.5. Указание. Для нахождения температуры u(t, r) необходимо решить смешанную задачу ut = a2u, u|G = 0, u|t=0 = xy(1 - x)(1 - y).

Ответ.

2 e-a 2[(2m+1)2+(2n+1)2]t u(t, r) = sin[(2m+1) x] sin[(2n+1) y].

(2m+1)3 (2n+1)m,n=5.6. Указание. Направим ось z по оси цилиндра. Условия задачи предполагают, что зависимости температуры от координаты z нет и граничные условия на торцах цилиндра, удаленных на бесконечность, не оказывают влияния на решение. Переходя к полярным координатам (r, ) : x = r cos, y = r sin, получаем смешанную задачу для температуры u = u(t, r, ):

ut = a2u, ur|r=r = q, u|t=0 = u0.

Ответ.

r 2 J0 k r0 ak t 1 r2 u(t, r) = u0 + 2qr0 a + - - e-( r0 ) t, 2 r0 4 r0 2 k J0(k) k=где k k-й корень уравнения J0() = 0.

5.7. Указание. Условие конвективного теплообмена на поверхности S имеет вид u -k = h (u - )|S, n S где k коэффициент теплопроводности; h коэффициент теплообмена; температура окружающей среды. С учетом этого смешанная задача для определения температуры u = u(t, r) в цилиндре имеет вид ut = a2u, ur|r=1 = 0, uz|z=0 = 0, (uz + u)|z=l = 0, = h/k, u|t=0 = u0 1 - r2.

Ответ.

J0(kr) 2 u(t, r, z) = u0 1 - 8 e-a kt kJ0(k) k= an l n + 2l2 nz e- l t cos, n (n + l + 2l2) l n=где k k-й корень уравнения J0() = 0; n n-й корень уравнения l tg =.

5.8. Указание. Смешанная задача для определения температуры u(t, r) имеет вид ut = a2u, u|r=r = 0, u|t=0 = 0.

r j1 xk r0 axk r Ответ. u(t, r, ) = A cos + 2 e-( r0 ) t, r0 k=1 xk j1(xk) где xk k-й корень уравнения j1(x) = 0.

5.9. u(t, x, y) = 3 cos( 5at) sin x sin 2y + cos(5at) sin 3x sin 4y.

5.10. Указание. Смешанная задача имеет вид p utt = a2u + sin t, u|r=r = 0, u|t=0 = 0, ut|t=0 = 0, где поверхностная плотность мембраны.

axk Ответ. Нерезонансный случай: пусть = k =, где xk k-й rкорень уравнения J0(x) = 0, тогда r p0 (k sin t - sin kt) J0 xk ru(t, r) = 2.

k (k - 2) xkJ1(xk) k=Резонансный случай: при = n получим r p0 (t n cos nt - sin nt) J0 xnru(t, r) = - + n xn J1(xn) (5.45) r p0 (k sin nt - n sin kt) J0 xk r+ 2.

k (k - n) xkJ1(xk) k=k =n В этом случае амплитуда колебаний возрастает линейно по t.

j1(xkr) 5.11. u(t, r, ) = 2A cos (cos xkt - 1), x3 j0(xk) k k=где xk k-й корень уравнения j1(x) = 0.

5.12. Указание. Возникающие в сосуде акустические колебания газа описываются потенциалом скоростей u(t, r): u(t, r) = v(t, r). В момент времени t = 0 все частицы имеют одинаковую скорость v0. Направляя ось z параллельно v0, получим u|t=0 = v0z. На границе сосуда радиальная компонента скорости равна нулю. Следовательно, смешанная задача для u(t, r) имеет вид utt = a2u, ur|r=1 = 0, u|t=0 = v0r cos, ut|t=0 = 0.

j1(xkr) Ответ. u(t, r, ) = 2v0 cos cos xkt, (x2 - 2) j1(xk) k k=где xk k-й корень уравнения j1(x) = 0.

2v0 sh k(a - x) sh ky (2k + 1) 5.13. u(x, y) =, где k =.

b k sh ka b k=5.14. Указание. Краевая задача для определения u(x, y) имеет вид u = -q/k, u|x=0 = 0, u|y=0 = 0, ux|x=a = 0, uy|y=b = 0.

Здесь k коэффициент теплопроводности внутри пластины. Для выбора базиса возможны три варианта: из задачи на собственные значения по переменной x, по переменной y или по обеим переменным.

Ответ.

2q ch m(y - b) sin mx u(x, y) = 1 - = ka ch mb m m= 2q ch n(x - a) sin ny = 1 - = kb ch na n n= 4q sin mx sin ny =.

2 kab mn (m + n) m,n=5.15. Указание. Функцию u(r, ) можно найти либо в виде ряда, используя формулы (5.39) и (5.43), либо по формуле Пуассона (5.44).

Ответ.

2m+4T0 r sin(2m + 1) u(r, ) = - = r0 2m + m=2T0 2rr= - arctg sin.

r0 - rAr2 5.16. u(r, ) = r0 - r2 sin 2.

ln r 5.17. u(r) = u1 + (u2 - u1).

ln 5.18. Указание. Функция u(r) является решением краевой задачи u = 0, -kuz|z=0 = q, u|r=r = 0, u|z=h = 0, где k коэффициент теплопроводности.

k(h - z) kr 2qr0 sh r0 J0 rОтвет. u(r, z) =, kh k k J1(k) ch k=rгде k k-й корень уравнения J0() = 0.

5.19. Указание. Потенциал u(r) является решением краевой задачи u = 0, u|r=a = U0 z, uz|z=0 = 0, uz|z=h = 0.

U0h 4U0 I0(kr) cos kz (2k + 1) Ответ. u(r, z) = -, где k =.

2 h I0(ka) k h k=5.20. Указание. Разлагая неизвестную функцию u = u(r,, ) в ряд Фурье по сферическим гармоникам и оставляя только ограниченные при r = 0 радиальные функции, получим l l r u(r,, ) = Alm Ylm(, ).

a l=0 m=-l Константы Alm необходимо определить, используя граничные условия.

В силу аксиальной симметрии задачи (граничные условия не зависят от угла ), отличны от нуля только Al0:

/ 2l + Al0 = f() Yl0(, ) sin dd = 2 U0 Pl(cos ) sin d.

0 0 Далее следует использовать рекуррентное соотношение для полиномов Лежандра:

Pl(x) = (Pl+1(x) - Pl-1(x)), P-1(x) = 0, 2l + и частные значения полиномов Лежандра:

(2l)! Pl(1) = 1, P2l+1(0) = 0, P2l(0) = (-1)l.

22l(l!)Ответ.

2l+U0 4l + 3 (2l)! r u(r, ) = + U0 (-1)l P2l+1(cos ).

2 l + 1 22l+2(l!)2 a l=Глава Обобщенные функции В математической физике обобщенные функции (ОФ) возникают совершенно естественным образом. Действительно, в основе аналитических методов решения задач математической физики лежит теорема о разложении в ряд Фурье, которая позволяет свести задачу о нахождении искомого решения к нахождению его коэффициентов Фурье. Фактически наряду с известными способами задания функции явным, неявным и параметрическим теорема Фурье дает еще один способ, а именно, посредством задания набора ее коэффициентов Фурье, которые очень просто выражаются через исходную функцию. Так, например, для функции f(x) L2[a, b] коэффициенты Фурье b ck = (f, k) = f(x) k(x) dx, (6.1) a где {k(x)} ортонормированный базис. При этом функции f(x) однозначно соответствует счетный набор чисел ck (обратное утверждение, как известно, неверно). Аналог этого в конечномерном случае задание вектора посредством указания его декартовых координат (проекций вектора на базисные орты). Элементарные алгебраические операции (сложение и умножение на числа) совершенно естественно переносятся с функций на коэффициенты Фурье.

Попытаемся посмотреть на коэффициенты Фурье с другой, более общей точки зрения. Фактически формула (6.1) для вычисления коэффициентов Фурье определяет некоторое правило, по которому каждой функции базиса однозначно соответствует определенное число:

b f k(x) - ck = (f, k) = f(x) k(x) dx. (6.2) a Такое правило в математике принято называть функционалом. Таким образом, функцию можно задать посредством задания некоторого функционала, в данном случае линейного, на функциях базиса. Это утверждение есть просто изменение терминологии, и само по себе не дает ничего нового.

Принципиально важен следующий шаг. На данном этапе областью определения функционала является множество функций, определенных на фиксированном интервале и образующих там ортонормированный базис. Оба ограничения весьма неудобны, так как в различных задачах приходится иметь дело с разными интервалами и разными базисами. Формально снять эти ограничения несложно, нужно просто расширить область определения функционала, т. е., например, определить его тем же выражением на всей вещественной оси и на достаточно широком множестве функций:

+ f (x) - (f, ) = f(x) (x) dx. (6.3) При этом, однако, возникает проблема существования интеграла от произведения функций (попросту говоря, чем хуже свойства функции f(x), тем лучше должны быть функций (x)). К счастью, компромисс возможен, и множество функций, образующих область определения функционала, называют пространством основных функций (точное определение будет дано далее).

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

f (x) - (f, ), (6.4) однако потребуем, чтобы, как и в случае интеграла, функционал был линейным и непрерывным (эти понятия будут строго определены далее). Теперь появляется возможность наравне с обычными функциями рассматривать и новые объекты. Обозначение для функционала (f, ), совпадающее с обозначением скалярного произведения, просто дань традиции, и по этой же причине множество всех линейных непрерывных функционалов называют пространством ОФ.

В итоге мы перешли от обычных функций к ОФ. Здесь уместна аналогия: переход от вещественных x R чисел к комплексным z C.

Каждому x естественно соответствует z, а именно то, для которого Re z = x и Im z = 0. Обратное неверно: числу z с Im z = 0 нет соот ветствия в R. Далее, для комплексных чисел можно определить целый ряд действий, привычных для вещественных чисел, причем так, что в частном случае вещественных чисел, т. е. если Im z = 0, старые правила остаются неизменными. При этом некоторые действия (например, сложение) выглядят вполне естественно, а другие (например, умножение и деление) по форме выглядят весьма непривычно. Более того, некоторые понятия (например, понятие неравенства чисел) вообще не имеют аналога для комплексных чисел.

То же происходит и при переходе от обычных функций к ОФ: некоторые действия (сложение и умножение на числа) переносятся вполне естественным образом, другие же (например, операция дифференцирования) вводится весьма необычным образом, а третьи (например, умножение), вообще говоря, не определены для ОФ. В итоге, используя эти правила действий, удается перейти от УМФ для обычных функций к уравнениям для ОФ и, самое главное, указать новые общие методы решения краевых задач математической физики. В настоящее время теория ОФ общепризнанно является основой современной математической физики.

Важно отметить, что зачатки понятия ОФ появились в самом начале создания математического аппарата физики. Действительно, некоторые активно используемые в физических теориях понятия, такие, как точечный заряд, материальная точка, мгновенный импульс и т. п., не могут быть описаны обычными функциями.

Рассмотрим, например, материальную точку с массой m = 1 в начале координат, r = 0, и поставим задачу: найти функцию (r), которая описывает это распределение массы. Для (r) естественно потребовать, чтобы 1, если 0 V, (r) dr = (6.5) 0, если 0 V.

/ V Попробуем, учитывая определение материальной точки, такой вариант:

C, если r = 0, (r) = (6.6) 0, если r = 0.

Но для такой функции (r) dr =0 при любом V и любом выборе C.

Pages:     | 1 |   ...   | 9 | 10 || 12 | 13 |



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.