# MIT 1805:  Studio 7 Script runsensorsim.r

# Setup the Signal/Sensor Simulation and 
# Run a nice sequence of calls to the function sensorsim


# 0. Simulation Setup ----
# 0.1 Define S to be 11 x 11 grid of sensors in row by column grid ----

source('sensorsim.r')

n = 11;
S = matrix(0,nrow=n*n,ncol=2);
for (j in c(1:n))
{
  for (k in c(1:n))
  {
    #row = j*(n-1) + k;
    row=(j-1)*n + k;
    S[row,1]=j;
    S[row,2]=k;
  }
} 

# 0.2 Define signal location ----
t0=c(3.5,3.5);

# 1. Uniform prior case ----
# 1.1. Plot Sensor Grid, Signal Location, and Prior density ----
msg = 'uniform prior';
sensorsim(S, t0,1,0,0,msg)

# 1.2. Plot Sensor Grid, Signal Location, and Posterior densities ----
#   updating  one-signal per step for 5 steps 
msg = 'posteriors:  1 signal per step (10 steps) ';
sensorsim(S, t0,1,1,5,msg)

# 1.3  5-signals per step for 10 steps
msg = 'posteriors:  5 signals per step (10 steps) ';
sensorsim(S, t0,1,5,10, msg)

# 2. Normal/Gaussian Prior ----

msg = 'Normal/Gaussian prior';
sensorsim(S, t0,0,0,0, msg)
  
msg = 'posteriors:  1 signal per step, 10 steps ';
sensorsim(S, t0,0,1,10,msg)
    
msg = 'posteriors:  5 signals per step 10 steps ';
sensorsim(S, t0,0,5,10,msg)
  
msg = 'posteriors:  10 signals per step, 5 steps ';
sensorsim(S, t0,0,10,5,msg)
  
t0 = c(5,15)
msg = 't0 outside sensors: posteriors:  1 signal per step, 10 steps ';
sensorsim(S, t0,0,1,10,msg)

t0 = c(5,20)
msg = 't0 outside sensors: posteriors:  1 signal per step, 10 steps ';
sensorsim(S, t0,0,1,10,msg)

t0 = c(12,20)
msg = 't0 outside sensors: posteriors:  1 signal per step, 10 steps ';
sensorsim(S, t0,0,1,10,msg)

# Other Cases to try:
#  t0 = c(2.1,8.1)