# How to transform the diffraction image into the function of light intensity with the order of diffraction? [closed] In our pre-experiment, I got the image. It’s obvious that the order of diffraction from the top down is 3,2,1,0,-1,-2,-3, the background of optics is so easy.

Then I have to get the light intensity of the diffraction spot, denote as I(n),n in -3 to 3(in this experiment).

How to think about the problem and manipulate the processes with Mathematica? Thanks so much!

If it’s a little hard, in the future experiment, we will use CCD, with which we can get the direct image as follows: You can also just analyze this simple case.

(Forgive me for my poor English)

• The Mathematica StackExchange is not a free coding service. You need help from an expert consultant. Apr 25, 2021 at 7:50
• The first image is atrocious. It should have been taken flat on, from further away, and zoomed in because it has bad lens distortion / curvature, and will need to be inverse-perspective transformed too. Apr 25, 2021 at 13:37
• @Roman so sorry for that when I was new that day. Thank you for your advice. Sep 23, 2021 at 4:51

Consider a single slit illuminated by a parallel light beam so that at time t all the points in the slit have the same phase. The light from the slit that falls on a screen parallel to the slit is considered. The near filed is complicated, but the field far away, where the angle to some point on the screen is approx. the same for all points in the slit, is much easier. We may then also neglect the effect of changing intensity due to different way length. Further, at time t al the points in the slit have the same phase. But due to different way lengths, the phase of light coming from different slits position does have different phases. This leads to interference.

Here is an simple example for relative intensities. k denotes the wave length of the light, d the distance between slit and screen, x are coordinates of the slit, y of the screen:

k = 10 ;
dist = 1000;
intensity[y_] :=
Norm[NIntegrate[
Exp[2 Pi I Sqrt[dist^2 + (x - y)^2] k], {x, -.5, .5}]]^2
Plot[intensity[y], {y, 0, dist/2}] 