Двойной интеграл в декартовой координате вместо (R, Theta)
мой предыдущийвопрос (Функция интеграла интенсивности в питоне)
Вы можете увидеть дифракционную модель на изображении ниже:
Я хочу вычислить интеграл интенсивности в каждом пикселе (квадрате), поэтому я не могу использовать R и Theta в качестве переменной. Как я могу сделать это в координате X-Y.
Наша функция:
вместо греха (тэта) мы можем использовать:
sintheta= (np.sqrt((x)**2 + (y)**2)/(np.sqrt((x)**2 + (y)**2 + d**2)))
Другие константы:
lamb=550*10**(-9)
k=2.0*np.pi/lamb
a=5.5*2.54*10**(-2)
d=2.8
когда вы строите график, результат выглядит примерно так: (изображение выше - вид сверху)
метод в предыдущем разделе: вычислить интегрирование функции в (0.0, dist) и после этого * (2 * np.pi).х) который х = кa * np.sin (theta), но теперь я хочу интегрировать в каждый пиксель. который предыдущий метод не работает, потому что это координата X-Y, а не полярная.