от
Я решил аналитическое уравнение Лапласа для двумерного квадратного квадрата и хотел бы построить векторы электрического поля как функцию от у и х. https://www.youtube.com/watch?v=ws3LIfWxfcM Это хорошая ссылка, чтобы показать конечное уравнение, которое я использовал. Я дифференцировал решение для потенциала, чтобы найти электрическое поле. Проблема, которую я получаю, состоит в том, что мой цикл суммирования, который я настроил, не сходится при увеличении n, следовательно, для n выше 50 я получаю ответы для электрического поля порядка 10 ^ 100. Поэтому моя программа не отображается, так как значения очень велики. Мой цикл суммирования до сих пор работал хорошо. Например, я построил график f (x, y) = sin (x) i ^ cos (y) j ^ и сравнил их с разложением по Тейлору sin (x) и cos (x). Я считаю, что проблема заключается в делении на sinh (npi). Когда я добавляю этот термин, он быстро становится бесконечно большим числом и отказывается от сюжета. импортировать matplotlib.pyplot как plt импортировать NumPy как NP импорт математики x = np.arange (-1,1,0.1) y = np.arange (-1,1,0.1) X, Y = np.meshgrid (x, y) a = np.zeros (len (X)) для n в диапазоне (50):     i = (np.cos (nnp.pi) -1) (- np.cos (nnp.piX) np.sinh (nnp.pi * Y)) * 1 / (np.sinh (np.pi)) #x Компонент электрического поля
a = np.add((i), a)
b = np.zeros (len (Y)) для n в диапазоне (50):     j = (np.cos (nnp.pi) -1) (np.sin (nnp.piX) np.cosh (nnp.pi * Y)) * 1 / (np.sinh (np.pi)) #y компонент электрического поля
b = np.add((i), a)
fig, ax = plt.subplots () когда я печатаю a или b, возвращается q = ax.quiver (X, Y, a, b) plt.show () Я ожидаю, что вектор электрического поля будет очень большим у «стен».              

Ваш ответ

Отображаемое имя (по желанию):
Конфиденциальность: Ваш электронный адрес будет использоваться только для отправки уведомлений.
Анти-спам проверка:
Чтобы избежать проверки в будущем, пожалуйста войдите или зарегистрируйтесь.
...