This repository contains the supporting code for the article Interferometric Focal Planes (under review). The code implements Monte Carlo simulations to estimate parameters in several prototype objects, including:
- Two point sources
- A line source
- An edge source
- A right-angle source
These simulations explore fixed and adaptive measurement strategies for mutual intensity measurement and compare them with the estimation results from intensity measurements. Jacobian matrices of the forward models are also calculated to determine the degrees of freedom of the system.
SimulationCode/
βββ Corner/
β βββ CornerEstimation.m
β βββ JacobianCorner.m
β βββ MeasureCorner.m
β
βββ Edge/
β βββ EdgeEstimation.m
β βββ JacobianEdge.m
β βββ MeasureEdge.m
β
βββ Line/
β βββ Jacobian1Line.m
β βββ LineEstimation.m
β βββ MeasureLine.m
β
βββ twopoints/
β βββ Jacobian2P.m
β βββ Measurel2P2D.m
β βββ two_pts_estimation.m
β
βββ README.md
The code is organized into different folders, each corresponding to a prototype object.
Each folder contains three MATLAB scripts:
-
Main Program (
*_Estimation.m)- Runs Monte Carlo simulations for parameter estimation.
- Implements fixed and adaptive measurement strategies.
- Uses a forward model to simulate measurements.
- Includes visualization code to generate comparison plots.
-
Forward Model (
Measure*.m)- Defines the physical measurement process for the given object.
- Used in the Monte Carlo simulations to generate synthetic data.
-
Jacobian Calculation (
Jacobian*.m)-
Computes the Jacobian matrix of the forward model.
-
Not used in Monte Carlo simulations, but provides theoretical insights.
-
Used for visualizing degrees of freedom in parameter estimation.
-
-
Clone or download this repository.
-
Open MATLAB in the
SimulationCodedirectory. -
Choose a geometry and run its main script. For example, for the corner (right-angle) source:
run Corner/CornerEstimation.m
This code is written in MATLAB R2022b and requires the following toolboxes:
- Parallel Computing Toolbox (for
parfor, used in Monte Carlo simulations) - Optimization Toolbox (for
fminunc, used in Maximum Likelihood Estimation)
If running the code without a Parallel Computing Toolbox, replace parfor loops with for loops to ensure compatibility.
If you use this code, please cite our article Interferometric Focal Planes (once published).
License: This code is provided for research purposes only.