- Create a shape and save it by running defineShape.m. You can change the shape primarily by changing values of pvec. But keep in mind that first 3 values have to be same as last 3 values for curve to be closed. And length of the vector is also important. Rules: length(tvec) = length(pvec)+p+1 and length(pvec) = length(wvec)
- Call lightCurve_visual-function with these parameters:
- accuracy: accuracy of the measurement (How many light-rays)
- angle: angle between viewer's direction and light source
- folder: folder, where animation frames are saved
Example usage:
- Open defineShape.m
- Run it
- Create folder 'my_frames'
- Call a function lightCurve_visual(30, 3*pi/4, 'my_frames')
- Simulate MCMC shape reconstruction by running function inverse_problem_mcmc(N, folder, acc, frame_freq)
- Parameters:
- N: number of MCMC iterations
- folder: folder, where animation frames are saved
- acc: in light curve simulation, the number of light rays
- frame_freq: how frequently you want to record a frame (every frame_freq:th iteration will be recorded. In addition frame will be recorded whenever new maximum value likelihood occurs)
Example:
- Create folder for animation frames ie. 'mcmc_frames'
- Run inverse_problem_mcmc(1000, 'mcmc_frames', 10, 1)