meine früherenFrage (Integral der Intensitätsfunktion in Python)

Sie können das Beugungsmodell im Bild unten sehen:

Ich möchte das Integral der Intensität in jedem Pixel (Quadrat) berechnen, daher kann ich R und Theta nicht als Variable verwenden. Wie kann ich das in X-Y-Koordinaten machen?

Unsere Funktion:

anstelle von sin (theta) können wir verwenden:

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

Andere Konstanten:

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

Wenn Sie eine Funktion zeichnen, sieht das Ergebnis ungefähr so aus: (Das Bild oben ist eine Ansicht von oben)

Die Methode im vorherigen Thema: Berechnen Sie die Integration der Funktion in (0.0, dist) und danach * (2 * np.pix) welches x = ka * np.sin (theta), aber jetzt möchte ich in jedes pixel integrieren. was die vorherige Methode nicht funktioniert, da dies eine X-Y-Koordinate ist, nicht polar.