Integral duplo na coordenada cartesiana em vez de (R, Theta)

meu anteriorPergunta, questão (Função Integral de Intensidade em python)

você pode ver o modelo de difração na imagem abaixo:

Eu quero calcular a integral de intensidade em cada pixel (quadrado), então não posso usar R e Theta como variável. Como posso fazer isso na coordenada X-Y.

Nossa função:

em vez de pecado (teta), podemos usar:

sintheta= (np.sqrt((x)**2 + (y)**2)/(np.sqrt((x)**2 + (y)**2 + d**2)))

Outras constantes:

lamb=550*10**(-9)
k=2.0*np.pi/lamb
a=5.5*2.54*10**(-2)
d=2.8

quando você plota a função, o resultado é mais ou menos assim: (A imagem acima é uma visão de cima)

o método no tópico anterior: calcule a integração da função em (0,0, dist) e depois * (2 * np.pix) qual x = ka * np.sin (theta), mas agora quero integrar em cada pixel. cujo método anterior não funciona, porque essa é a coordenada X-Y, não polar.

questionAnswers(1)

yourAnswerToTheQuestion