[МУЗЫКА] Вот,
а теперь, наконец, перейдем к задаче об определении
сопротивления между двумя разными узлами на сетке.
Сетка сопротивлений.
[БЕЗ_ЗВУКА] Значит,
еще раз, имеется сетка сопротивлений,
сопротивление между двумя узлами
это R, сопротивление вот такого звена это R.
И мы хотим найти, вообще говоря,
сопротивление между двумя произвольно заданными узлами сетки.
Главное — выбрать правильные переменные, и в этой задаче я предлагаю в качестве
переменных использовать электростатический потенциал в точках x и y.
Будем считать, что для простоты размер сетки равен единице,
чтобы с лишними буквами жизнь себе не усложнять, тогда x и y,
они принимают целочисленные значения.
И вот на этот потенциал нужно написать уравнение и его решить.
Ну, поехали.
Какое уравнение?
Естественно, это уравнения Кирхгофа, которые говорят,
что ток у нас в каждом узле сохраняется.
Давайте мы теперь напишем, что мы ожидаем для тока,
который в каком-то узле собирается.
Значит, это будет уравнение в левой части и в правой части.
Уравнение будет длинное.
Давайте начнем с левой части.
Вот я нахожусь здесь, в точке x, y.
И я смотрю, какой ток у меня протекает по 4-м ребрам и я его суммирую.
Значит, нужно, скажем, вот это ребро, да?
Здесь x меняется с x до x + 1, поэтому,
чтобы посчитать разность потенциалов вот здесь,
я пишу x + 1 y − вот здесь и так по всем четырем ребрам.
Поэтому в результате получается такое выражение: φ от x − 1, y.
Это я работал этот узел и этот узел,
потом я добавляю узлы по вертикали.
Я оставляю x, меняю y вначале на + 1,
потом на −
1 и вычитаю
4 значения потенциала в точке x, y.
И это у меня тот ток,
который перетекает, втекает в этот узел по ребрам.
Теперь как ставится задача?
Задача ставится так, что в точку, есть две выделенные точки.
Есть точка, ну, скажем, пусть эта точка будет 0, 0,
куда приделана
приделан один контактометр И имеется точка,
скажем, X, Y, куда приделан второй контактометр.
Соответственно, вот этот вот ток, подходящий к вершине,
он равен 0 по всем узлам, кроме вот этих двух, значит,
я должен дописать сюда, в правую часть выражения,
вклад от этих двух точек, где вот этот вот ток не равен 0.
и я напишу так.
Это есть I * R, слева стоит потенциал, поделенный на R, дает мне ток,
мне удобней написать сюда, и есть два вклада, первый говорит,
что вот есть ток отсюда, и это есть дельта-функция, которая говорит,
что она отлична от нуля, точнее, символ Кронекера, x = X,
а y = Y.
И минус вот отсюда,
символ Кронекера говорит, что x = 0 и y = 0.
Вот так выглядит исходное уравнение, которое требуется решить.
В принципе, если я его смогу решить,
то я найду распределение потенциала вообще по всей сетке.
Это такое дискретное уравнение.
И будем его решать.
Ну, понятно, что у нас тут все узлы зацеплены, и просто так не решишь,
но можно воспользоваться аналогией с предыдущей задачей,
если сюда посмотреть, то вы с легкостью обнаружите,
что здесь написано такое дискретное дискретный оператор Лапласа,
который обобщает вторую производную, на случай функций, заданных на решетке.
Поэтому теперь я говорю, что раз так, то я покажу,
перехожу в представление Фурье, и смотрите,
что у меня происходит с левой частью.
Значит, я перешел представление Фурье,
ввел букву Φ, то есть я теперь написал,
что у меня давайте я его еще раз повторю: я написал, что Φ в точке x и y
есть интеграл, по зона Бриллюэна, давайте посмотрим, как устроена зона Бриллюэна.
Давайте посмотрим.
Значит, теперь у нас двумерная задача, поэтому у меня имеются qx и qy,
это qx, это qy, и зона Бриллюэна,
размер ячейки 1, зона Бриллюэна это такой квадрат,
вот это вот π это π, это − π и это − π.
И я по такому квадрату должен проинтегрировать.
Так устроена зона Бриллюэна.
То есть я должен проинтегрировать по зоне Бриллюэна.
[БЕЗ_ЗВУКА] Здесь я напишу Φ от qx,
qy, сюда я должен написать e в степени i
qxX + i qyY
и здесь стоит d2q / 2π в
квадрате потому что двумерная задача.
Я так сделал, теперь я подставляю эту формулу сюда и смотрите, что получается.
Как и в предыдущей задачей с цепочкой спиннов, здесь x сдвигается на 1,
здесь на − 1, и это генерит нам следующие члены.
Значит, справа мы имеем e в степени из первого i qx вот отсюда.
Здесь сдвигается в другую сторону, получается e в степени − iqx.
Здесь то же самое, только в направлении y,
поэтому я должен написать e в степени i qy e в степени
− iqy, а отсюда − 4.
И все это на Φ от q.
Это левая часть.
Теперь что стало с правой частью.
С правой частью, смотрите, у нас здесь стоит символ Кронекера,
поэтому, когда я буду считать преобразование Фурье,
получится следующее: здесь будет стоять
I R на e
в степени − iqx
/ X, − iqyY,
отсюда будет просто − 1.
Так переписывается исходное уравнение
Кирхгофа в терминах Фурье-образа.
Ну тут надо немножко поработать и записать ответ
уже в более читаемом виде, он устроен так: Ф
(q) есть и это никуда не девается I R.
Здесь дополнительный минус есть, и здесь минус этот я сокращу,
1 − e в степени − iqx
* X − iqy * Y,
а в знаменателе вот отсюда
стоит 2 1 − cos qx, лучше 2,
наверное, написать, − cos qy.
Это то, что пришло из левой части.
Ну вот смотрите, после этой процедуры мы видим, что, в принципе,
мы знаем Фурье-образ потенциала.
То есть задача почти решена.
Имеется явное выражение, которое, в зависимости от того,
чему равен вот этот x и y, где мы ток снимаем, показывает нам Ф от q.
Ну, хорошо, это замечательно, но это еще не совсем то, что нам нужно,
а хочется нам посчитать Ф,
вернуться в реальное пространство и нам Ф везде не нужен,
нам нужно посчитать разность потенциалов между вот этой точкой и вот этой точкой.
И если это мы найдем, тогда мы найдем эффективность сопротивления сетки между
этими двумя точками.
Ну хорошо, поехали.
Считаем эффективное напряжение, как Ф 0,
0 − Ф от X,Y.
Здесь минус здесь.
Хорошо.
Что для этого надо сделать?
Вот написано выражение.
Нужно написать интеграл по зоне Бриллюэна с вот такими экспонентами и с Ф,
которое вот мы здесь уже нашли.
Теперь смотрите, когда мы пишем Ф 0, 0 в начале координат,
вот здесь тогда экспонента вообще не работает, потому что x y = 0.
Теперь, когда мы считаем в этом месте, у нас получается, где X и Y, они большие.
Значит, что у нас в результате получается?
Ну, давайте мы напишем предварительное выражение.
Имеется интеграл по зоне Бриллюэна,
Ф (q)
1 −
eiqx *
X + iqy *
Y d2q
/ 2π2 [МУЗЫКА]
[МУЗЫКА]