Главная > Л.И.МАНДЕЛЬШТАМ. ТОМ IV
<< Предыдущий параграф Следующий параграф >>
Пред.
След.
Макеты страниц

Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше

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

Если задана внешняя сила
\[
g(x, t)=g(x) \cos \sqrt{\lambda} t
\]
( $\lambda$ здесь дано) и смещение стержня
\[
y(x, t)=\varphi(x) \cos \sqrt{\lambda} t,
\]

то, как мы видели, для $\varphi(x)$ получается интегральное уравнение
\[
\varphi(x)=\lambda \int_{0}^{l} V(x, \xi) q(\xi) \varphi(\xi) d \xi+f(x) \text {, }
\]

где
\[
f(x)=\int_{0}^{l} V(x, \xi) g(\xi) d \xi
\]
-заданная функџия $x$. Если нет внешней силы, интегральное уравнение принимает вид
\[
\varphi(x)=\lambda \int_{0}^{l} V(x, \xi) q(\xi) \varphi(\xi) d \xi .
\]

Здесь $\lambda$ должно быть определено из самого интегрального уравнения.
$V(x, \xi)$ – функџия Грина нашей задачи. Функџия $V(x, \xi) q(\xi)$ называется ядром интегрального уравнения. На основании физического смысла функџии Грина ее часто называют функџией влияния.

Мы нашли $V(x, \xi)$ для двух частных случаев: для стержня с закрепленными кондами и стержня с одним закрепленным концом. Если оба конда свободны, то функџию Грина построить нельзя. Это очень типичное и интересное обстоятельство. Постараемся его связать с уже известными нам свойствами.

Если стержень свободен, производная $\frac{\partial y}{\partial x}$ равна нулю на обоих конџах. Кроме того, при равновесии $\frac{\partial^{2} y}{\partial x^{2}}$ равно нулю всюду, кроме точки приложения силы. Таким образом, $y$ должно быть постоянно на протяжении всего стержня, и нельзя удовлетворить условию скачка производной в точке приложения силы. Подчеркнем, что несуществование функџии Грина связано не с дифферендиальным уравнением, а с граничными условиями.

В случае свободных конџов одно из собственных значений есть нуль. Уравнение
\[
\varphi^{\prime \prime}+\lambda q \varphi=0
\]

превращается при $\lambda=0$ в уравнение
\[
\varphi^{\prime \prime}=0
\]

с непрерывным решением
\[
\varphi=\text { const. }
\]

Если существует непрерывная функџия, удовлетворяющая уравнению и граничным условиям, то не существует удовлетворяющей ему разрывной фунхции. Наоборот, если не существует непрерывной функџии, удовлетворяющей дифференџиальному уравнению и граничным условиям, то существует разрывная.

Существование или несуществование функџии Грина имеет ясный физический смысл. Если стержень свободен на обоих конџах, то при действии на него силы не существует состояния равновесия. Физически этот случай является исключительным: стержень может не только колебаться, но и двигаться равномерно.

Однако для случая свободных конџов можно построить функџию Грина в обобщенном смысле. Для этого, кроме единичной сосредоточенной силы, нужно приложить к стержню еше компенсирующую ее распределенную силу.

Можно рассматривать стержень как совокупность масс, соединенных „пружинами“ (довольно трудно ответить на вопрос: что лучше отражает – действительный стержень, сплошная или дискретная модель). Мы уже рассматривали предельный переход от дискретной модели к сплошной и показали, что решение, получаемое для дискретной системы, в пределе переходит в решение некоторого дифферендиального уравнения в частных производных ${ }^{1}$. Интересно провести рассуждение так, чтобы от уравнений дискретной системы перейти к интегральному уравнению.

Возьмем „стержень“, закрепленный на одном конџе и свободный на другом (рис. 180), и поставим сначала статическую задачу (здесь также можно перейти от статической задачи к динамической с помощью принџипа Даламбера). Пусть на каждую точку действует некоторая сила. Обозначим через $\bar{g}_{k}$ силу, действующую на $k$-ую точку. Нужно определить смещения $y_{1}, y_{2}, \ldots, y_{n}$ наших $n$ частиџ.

Рассмотрим сначала частный случай, когда внешняя сила, равная единиџе, действует только на $k$-ую частиду.
${ }^{1}$ [См. 32-ю лекуию части I и 1-ю лекџию части II.]

ЛЕКЦИИ ПО КОЛЕБАНИЯМ. ЧАСТЬ ВТОРАЯ
Пусть $\alpha$-коэффиџиент упругости каждой пружины. Тогда результирующая упругая сила, действующая на $i$-ую частиџу $(i
eq k)$, есть
\[
\alpha\left(y_{i}-y_{i-1}\right)-\alpha\left(y_{i+1}-y_{i}\right)
\]
(на все частиџы, кроме $k$-ой, другие силы не действуют). При равновесии эта результирующая сила равна нулю, т. е.
\[
-y_{i-1}+2 y_{i}-y_{i+1}=0
\]

Напишем уравнение для первой частиџы:
\[
y_{1}-y_{2}=0
\]

Рис. 180.
и для последней:
\[
y_{n}-y_{n-1}=0 .
\]

Нужно еше написать уравнение для $k$-ой частијы. Здес результирующая упругая сила уравновешивается внешней силой, равной единице, т. е.
\[
-y_{k-1}+2 y_{k}-y_{k+1}=\frac{1}{\alpha} .
\]

Для $i<k$ смешения $y_{i}$ образуют арифметическую прогрессию:
\[
y_{i}=i r \quad(i<k) .
\]

Все частиды с $i \geqslant k$ имеют одинаковое смещение. Так как џепочка не имеет разрыва, то
\[
y_{i}=k r \quad(i \geqslant k) .
\]

Не трудно определить, что $r=1 / \alpha$. Таким образом, значения $y_{i}$ $(i=1,2, \ldots)$ являются функџиями номера $i$ и номера $k$ той частиџы, к которой приложена сила, а именно:
\[
V_{i k}=\left\{\begin{array}{ll}
\frac{i}{\alpha} & (i<k), \\
\frac{k}{\alpha} & (i \geqslant k) .
\end{array}\right.
\]

Если сила равна не единиџе, а $\bar{g}$, то $y_{i}$ просто умножается на $\bar{g}$. Если теперь силы действуют на все точки и $\bar{g}_{k}$ – сила, приложенная к $k$-ой точке, то
\[
y_{i}=\sum_{k=1}^{n} V_{i k} \bar{g}_{k} .
\]

Если
\[
\begin{array}{l}
\bar{g}_{i}=g_{i} \cos \sqrt{\lambda} t, \\
y_{i}=\varphi_{i} \cos \sqrt{\lambda} t,
\end{array}
\]

то, применяя принџип Даламбера, получаем окончательно для динамической задачи:
\[
\varphi_{i}=\lambda \sum_{k} V_{i k} m_{k} \varphi_{k}+\sum_{k} V_{i k} g_{k} .
\]

Здесь уместно поставить вопрос о том, как поступить, если мы хотим вычислить период колебаний распределенной системы, апроксимируя ее посредством дискретной системы. Подход, основанный на дифферендиальных уравнениях, проще. Имея в виду только простоту, можно было бы не переходить к анализу посредством интегральных уравнений.
Для перехода от (5) к интегральному уравнению напишем
\[
\alpha=\frac{p}{a},
\]

где $p$-макроскопический модуль упругости, $a$ – длина отдельной „пружины“. Тогда
\[
V_{i k}=\left\{\begin{array}{ll}
\frac{a i}{\alpha a} & (i<k), \\
\frac{a k}{\alpha a} & (i \geqslant k) .
\end{array}\right.
\]

Пусть теперь $i a=x$, причем
\[
a \rightarrow 0, \quad i \rightarrow \infty,
\]

и, кроме того,
\[
\frac{m_{k}}{a}=q(\xi), \quad a=\Delta \xi .
\]

Мы получаем в пределе
\[
\varphi(x)=\lambda \int_{0}^{l} V(x, \xi) q(\xi) \varphi(\xi) d \xi+\int_{0}^{l} V(x, \bar{\xi}) g(\xi) d \xi .
\]

Итак, исходя из дискретной модели и перейди к пределу (мы не вдаемся в вопрос о законности перехода), мы снова получили интегральное уравнение (3). Интегральное уравнение устанавливается на основании физической картины явления так же хорошо, как и дифференџиальное.

Пусть дана краевая задача, которая формулируется на языке дифференџиальных уравнений следующим образом:
\[
\begin{array}{c}
\frac{d}{d x}\left(p \frac{d \varphi}{d x}\right)+\lambda q \varphi=-g(x), \\
\left(\alpha_{1} \varphi-\alpha_{2} \frac{d \varphi}{d x}\right)_{0}=0, \\
\left(\beta_{1} \varphi+\beta_{2} \frac{d \varphi}{d x}\right)_{l}=0
\end{array}
\]
(общие граничные условия Штурма-Лиувилля). Для нее существуют определенные собственные числа и собственные функџии. Можно доказать, что все собственные функции и собственные числа краевой задачи удовлетворяют интегральному уравнению для той же задачи, и наоборот. Совокупность собственных чисел и собственных функций в обеих задачах одинакова.

Разумеется, то, что мы из физики получили для обоих случаев обе картины, в значительной степени предрешает этот результат, но мы постараемся его теперь доказать. Мы докажем, что решения интегрального уравненяя (3) удовлетворяют дифференџиальному уравнению (6) при краевых условиях (7) и (8), и наоборот.
Функџия $V(x, \xi)$ непрерывна и удовлетворяет уравнению
\[
\frac{d}{d x}\left(p \frac{d V}{d x}\right)=0 .
\]

Ее производная непрерывна всюду, кроме точки $x=\xi$, где она претерпевает скачок
\[
V^{\prime}(x, \xi)_{x=\xi-0}-V^{\prime}(x, \xi)_{x=\xi+0}=\frac{1}{p} .
\]

Уравнение (9) можно проинтегрировать в квадратурах. Мы найдем одно решение для левой части системы $(x<\xi)$, другое- для правой $(x>\xi)$. Обозначим эти решения соответственно через ? и б. Имеем:
\[
\frac{d}{d x}\left(p \rho^{\prime}\right)=0, \quad \frac{d}{d x}\left(p \sigma^{\prime}\right)=0,
\]

откуда
\[
\frac{d}{d x}\left[p\left(\sigma \rho^{\prime}-\rho \sigma^{\prime}\right)\right]=0,
\]

или, интегрируя,
\[
p\left(\sigma \rho^{\prime}-\rho c^{\prime}\right)=\text { const. }
\]

Для случая, когда $\lambda$ не есть собственное значение, всегда можно построить функџию Грина из следуюших соображений: $\rho$ и $\sigma$ выбраны так, чтобы $\rho$ удовлетворяло краевому условию (7), а $\sigma$-краевому условию (8). При этом $\rho$ и $\sigma$ линейно независимы $\left(\sigma \rho^{\prime}-\rho \sigma^{\prime}
eq 0\right)$ и при подходящей нормировке
\[
p\left(\sigma \rho^{\prime}-\rho \sigma^{\prime}\right)=1 .
\]

Возьмем теперь
\[
\begin{array}{ll}
V(x, \xi)=p(x) \sigma(\xi) & (0<x<\xi), \\
V(x, \xi)=p(\xi) \sigma(x) & (\xi<x<l)
\end{array}
\]

и проверим, что построенная таким способом функџия $V(x, \xi)$ удовлетворяет всем поставленным условиям.
Имеем, дифференџируя (10) и (11) по $x$ :
\[
\begin{aligned}
V^{\prime}(x, \xi)_{x=\xi-0} & =\rho^{\prime} \sigma, \\
V^{\prime}(x, \xi)_{x=\xi+0} & =\rho \sigma^{\prime}
\end{aligned}
\]

и, следовательно,
\[
V^{\prime}(x, \xi)_{x=\xi-0}-V^{\prime}(x, \xi)_{x=\xi+0}=\frac{1}{p} .
\]

Я утверждаю, что составленное с таким $V(x, \xi)$ интегральное уравнение эквивалентно дифференџиальной схеме (6) – (8). Можно показать обычным способом, что функџия $V(x, \xi)$ – единственная функция, удовлетворяющая всем поставленным условиям.

Покажем, что функџия $\varphi$, удовлетворяющая дифференџиальному уравнению и краевым условиям, удовлетворяет и найденному интегральному уравнению.
Имеем:
\[
\begin{array}{c}
\frac{d}{d x}\left(p \varphi^{\prime}\right)=-\lambda q \varphi-g(x) ; \\
\frac{d}{d x}\left(p V^{\prime}\right)=0 .
\end{array}
\]

Умножив (13) на $V$, (14) – на ф и вычитая, мы получаем:
\[
\frac{d}{d x}\left[p\left(\varphi^{\prime} V-V^{\prime} \varphi\right)\right]=-\lambda V q \varphi-V g .
\]

Будем теперь интегрировать по $x$ в интервале $(0, l)$. Так как при: этом нужно учитывать разрывность $V^{\prime}$, мы разделим интервал $(0, l)$ на интервалы $(0, \xi)$ и $(\xi, l)$ и будем интегрировать отдельно по каждому интервалу:
\[
\begin{array}{l}
\left.p\left(\varphi^{\prime} V-V^{\prime} \varphi\right)\right|_{0} ^{\xi-0}=-\lambda \int_{0}^{-0} V q \varphi d x-\int_{0}^{\xi-0} V g d x ; \\
\left.p\left(\varphi^{\prime} V-V^{\prime} \varphi\right)\right|_{\xi+0} ^{l}=-\lambda \int_{\xi+0}^{l} V q \varphi d x-\int_{\xi+0}^{l} V g d x .
\end{array}
\]

Левые части обращаются в нуль соответственно при $x=0, x=l$, так что после сложения (15) и (16) остается:
\[
-p \varphi(\xi)\left[V_{\xi-0}^{\prime}-V_{\xi+0}^{\prime}\right]=-\lambda \int_{0}^{l} V(x, \xi), q(x) \varphi(x) d x-\int_{0}^{l} V(x, \xi) g(x) d x .
\]

Учитывая (12), получаем далее:
\[
\varphi(\xi)=\lambda \int_{0}^{l} V(x, \xi) q(x) \varphi(x) d x+\int_{0}^{l} V(x, \xi) g(x) d x,
\]

откуда, после замены $x$ на $\xi$ и $\xi$ на $x$, следует уравнение
\[
\varphi(x)=\lambda \int_{0}^{l} V(\xi, x) q(\xi) \varphi(\xi) d \xi+\int_{0}^{l} V(\xi, x) g(\xi) d \xi .
\]

Это уравнение отличается от (3) лишь тем, что вместо $V(x, \xi)$ стоит $V(\xi, x)$, но так как функция Грина симметрична:
\[
V(x, \xi)=V(\xi, x)
\]
(и это здесь существенно), уравнение (17) совпадает с (3).
Итак, мы убедились, что всякая функция, удовлетворяющая схеме (6) – (8), одновременно удовлетворяет и интегральному уравнению (3). Можно также путем простого дифференџирования показать и обратное.

Дифференџиальная схема нашей задачи эквивалентна интегральному уравнению с функџией Грина, определенной выше. В этой эквивалентности заключается, с одной стороны, интерес, который представляют интегральные уравнения, с другой стороны, 一их слабость, если бы дело этим ограничивалось.

Из того, что было сказано, следует, что интегральное уравнение (3) имеет бесконечное множество собственных значений и каждому собственному значению соответствует одна собственная функџия.

Мы вывели интегральное уравнение, соответствующее определенному узкому классу дифференџиальных уравнений. Вид ядра связан определенным образом с задачей Штурма-Лиувилля. Для спеџиального типа интегральных уравнений, ядра которых соответствуют дифреренциальной схеме (6)-(6), мы доказали некоторые свойства собственных чисел и собственных функций. Однако необходима более общая теория, пригодная при любом ядре. Такая теория существует. Теория интегральных уравнений гораздо шире, чем проблема Штурма-Лиувилля. Ряд вопросов (например, вопрос о разложимости функџий по собственным функџиям) легче, проще и изящнее решается с помощью интегральных уравнений. Даже ради одного этого стоит их изучить.

Мне хочется изложить задачу, где ссылка на дифференџильные уравнения не поможет. Мы получим в ней интегральное уравнение другого типа, чем (3). Эта задача относится к теории оптического изображения ${ }^{1}$.

Пусть имеется система линз, дающая изображение (рис. 181). Каждой точке объекта соответствует определенная точка изображения. Предположим, что система не дает увеличения (или уменьшения). Согласно геометрической оптике яркость в точке $M$ определяет освещенность в точке $M$. Но в действительности геометрическая оптика несправедлива. Вследствие диффракџии каждой точке предмета соответствует в плоскости изображения некоторое распределение интенсивности. Максимум находится в точке, где было бы изображение согласно геометрической оптике, но с удалением от этой точки интенсивность спадает постепенно.

Пусть одна точка предмета дает в плоскости изображения распределение амплитуды
\[
K(x, \xi)
\]
${ }^{1}$ [Cр. том I, стр. 229.]

(амплитуда есть функџия координаты точки объекта $\xi$ и координаты точки изображения $x$ ). Если диафрагма имеет вид щели, то
\[
K(x, \xi)=\frac{\sin c(x-\xi)}{x-\xi},
\]

где
\[
c=\frac{2 \pi d}{L f}
\]
( $d$ – половина ширины щели; $L-$ длина волны; $f$ – фокусное расстояние). Если точка объекта перемещается, то перемещается вся картина; поэтому $K(x, \xi)$ есть функџия от $x-\xi$.
Амплитуда в точке $x$ от элемента предмета длины $d \xi$ есть
\[
z(\xi) K(x, \xi) d \xi,
\]

где $z(\xi)$ характеризует распределение амплитуды по предмету. Картины от отдельных Рис. 181. элементов предмета накладываются, и результирующая амплитуда в плоскости наблюдения имеет распределение
\[
A(x)=\int_{a}^{b} K(x, \xi) z(\xi) d \xi
\]
(здесь предполагается, что предмет освещен так, что отдельные его элементы посылают когерентные колебания).

Конечно, распределение $A(x)$, вообще говоря, не подобно распределению $z(\xi)$. Они только связаны между собой соотношением (18). Вообще говоря, при изображении предмет искажается, распределение света в изображении иное, чем в предмете. Возникают вопросы: как нужно строить оптический прибор для того, чтобы искажение было возможно меньше? Существуют ли такие предметы, которые не искажаются заданной оптической системой? Если бы оказалось, что существуют такие функџии $z(x)$, для которых
\[
A(x)=\frac{1}{\lambda} z(x),
\]

где $\lambda$-постоянная, то это означало бы, что распределения вида $z(x)$ не искажаются.

Подставляя (19) в (18), получаем:
\[
z(x)=\lambda \int_{a}^{0} K(x, \xi) z(\xi) d \xi .
\]

Неискажаемые структуры должны удовлетворять интегральному уравнению (20), ядро которого задано свойствами оптической системы.

Ядро уравнения (20) симметрично, но оно совсем другого типа, чем в задаче Штурма-Лиувилля. Оно не имеет разрыва производной.
Существует ли структура, удовлетворяющая уравнению (20)?
К сожалению, ответ неутешителен. Интегральное уравнение (20) имеет решения, но эти решения – бесконечные синусоиды. Только такие распределения изображаются точно. Таким образом, решение существует, но физически оно мало интересно.

Но можно поставить вопрос иначе: насколько можно приблизиться к точному изображению? Оказывается, можно указать такие структуры, которые мало искажаются. Можно высказать на этот счет ряд общих теорем.

Таким образом, области дифференџиальных и интегральных уравнений далеко не перекрываются. Интегральные уравнения представляют интерес и независимо от дифференциальных уравнений.

Categories

1
Оглавление
email@scask.ru