Questions tagged [optics]

Questions on the use of Mathematica to model light propagation and interaction with matter. This tag should be used for questions specifically related to the physics of light--problems dealing with ray-tracing, radiative transfer, lens equations, and electromagnetic field equations all belong here.

Filter by
Sorted by
Tagged with
0 votes
0 answers
71 views

How to preserve my constants

\[Sigma] = Constant a = Constant b = Constant \[Epsilon] = Constant T[x_] := (\[Sigma]/(2*\[Epsilon]))[a/(a^2 + x^2)^(1/2) - b/(b^2 + x^2)^(1/2)] In mathematica ...
Mohd. Farhan Hassan's user avatar
0 votes
0 answers
125 views

Nonlinear Differential Equation (Nonlinear Helmholtz) - How to solve?

I am relatively new to Mathematica, but I know that it is a powerful tool. I am looking to solve a nonlinear equation that could be seen in nonlinear optics. I have tried in MATLAB, but was struggling ...
jgk5141's user avatar
4 votes
2 answers
441 views

Finding shadow area of a 3D object

Given a 3D object, how could one find the area of the shadow? In effect, this is a measure of optical cross section. The ultimate goal would be a table showing latitude, longitude, and, optical cross ...
dantopa's user avatar
  • 1,050
1 vote
0 answers
50 views

How to solve an Integral analytically using a predefined definition for Besselfunctions (phi part of angular spectrum representation)

I'd like to use mathematica to calculate an Integral that is dependent on phi and theta (to obtain the intensity distribution of a tightly focused TEM20 mode using the angular spectrum representation)....
Franz's user avatar
  • 21
1 vote
1 answer
139 views

Parabolic Mirror - proving focal length

I would like to prove that there is a focal length in the parabolic mirror. I wrote script, but it's not working, can you guys see what I did wrong? ...
XYZCODE123's user avatar
2 votes
2 answers
2k views

Planck black body radiation curve

I have a data set of intensities and corresponding wavelengths. I want to fit it to the Planck black body radiation curve so that I can estimate the temperature for the given data . I have used the ...
rabink's user avatar
  • 111
8 votes
4 answers
4k views

Circular Aperture Diffraction Pattern

I am unsure if this is even possible in Mathematica, but my optics professor assigned a project to me on circular apertures and Fourier transforms. I have found plenty on the Airy function, but I am ...
Madeline Statkewicz's user avatar
6 votes
1 answer
412 views

single slit diffraction numerical simulation

I'm trying to simulate the Fraunhoffer diffraction at slits(single,double,triple) with Mathematica. In the picture, the red one is analytical result and the green one is numerical result. The ...
Taeshin Kim's user avatar
3 votes
0 answers
667 views

Computation of a Fresnel Diffraction pattern with Discrete Hankel Transform

In the next link: Computation of Hankel Transform using FFT (Fourier) Rainer implemented a great solution given in the next reference: Manuel Guizar-Sicairos and Julio C. Gutiérrez-Vega, "Computation ...
Eduardo's user avatar
  • 31
2 votes
0 answers
149 views

How to do the Fourier transform of a picture? [duplicate]

I want to generate Fourier transforms of a picture. the transform can represent the diffraction pattern of the picture. Basically, if I have a picture. Then, the result should look like this: For a ...
BNHSX's user avatar
  • 929
1 vote
1 answer
768 views

Optics Diffraction Grating Plot [closed]

I was unsure whether to post this in physics stackexchange or mathematica stackexchange, so I posted it in both. I'm trying make an intensity plot for a diffraction grating that contains 100 lines/mm....
Rosie's user avatar
  • 83
4 votes
1 answer
159 views

Adding extinction factor to Monte Carlo analysis

I'm looking to add an absorption component to this Monte Carlo analysis, essentially I want there to be a 1% chance that a photon gets absorbed at each event as it travels through the medium, once the ...
AW1991's user avatar
  • 153
1 vote
0 answers
168 views

Change in direction when a line hits an ellipse? [duplicate]

I'm trying to simulate reflecting light off an elliptical mirror. I currently store the ellipse in the form {center,{radii},rotationAngle} and the ray of light is ...
Steve D's user avatar
  • 2,189
1 vote
0 answers
454 views

Numerically solving a 2D oscillating integral [closed]

I'm having trouble solving this integral numerically: ...
mike's user avatar
  • 11
6 votes
2 answers
653 views

How to improve visualisation of trajectories in 3D space

I have a list of trajectories from a bunch of photons that travel inside a scattering medium. Each trajectory is a list of $(x,y,z)$ pairs. I'd like to visualise them in 3D space. I came up with the ...
stathisk's user avatar
  • 3,044
67 votes
1 answer
2k views

Is this the most effective and efficient programming style for prototyping (a ray tracer) in Mathematica?

I am trying out Mathematica as a prototyping tool. As an initial exercise I have put together a brute force ray tracer, with a view to using Mathematica's built in probability distribution functions ...
jrp's user avatar
  • 771
1 vote
1 answer
315 views

Export Graphics using Move in Optica

I am modeling an eye and am using TurboPlot and DrawSystem. I would like to export the graph, but using the ...
user12162's user avatar
5 votes
1 answer
513 views

Volumetric Rendering of cylinder symmetric radial image

I would like to render a cylinder symmetric radial cross-section with Image3D. As an example one could use the propagation of a Bessel Beam through a Lens: Here the cylinder axis is the bottom of the ...
Rainer's user avatar
  • 2,861
16 votes
1 answer
7k views

Computation of Hankel Transform using FFT (Fourier)

To address circularly symmetric cases of 2-D Fourier Transformations, the so-called Hankel Transform can be applied (for a detailed derivation of the relation between the 2-D Fourier transform and the ...
Rainer's user avatar
  • 2,861
0 votes
1 answer
547 views

Finding resolution of an optical system using Rayleigh criterion

I want to use the Rayleigh criterion to find the resolution of an optical system. I have two functions f1 = Sinc[x]^2 and ...
Darwin's user avatar
  • 5
27 votes
2 answers
3k views

How can I create a glass distortion effect in an image?

I'd like to overlay a glass jar onto an image with realistic light bending. Can anyone think of a way to automate this effect? Perhaps this can be done with the raytracing package Rayica / Optica?
M.R.'s user avatar
  • 31.2k
4 votes
3 answers
1k views

Drawing graphics part by part

What I want to do is to draw graphics part by part, because I made a function that generates coordinates according to which I will draw graphics. The problem is that MMA does not operate as I expected....
balboa's user avatar
  • 683