%
% Written by:
% -- 
% John L. Weatherwax                2005-08-14
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

clear; close all; clc; 
clear functions; 

Htilde = [ 0, 1; 3, 4 ];

Z1 = [ 0, 1, 7, 8, 5, 7, 9, 10, 6, 4 ]; Z1 = Z1(:); 
Z2 = [ 10, 7, 4, 5, 5, 3, 0, 2, 2, 4 ]; Z2 = Z2(:); 
Z = [ Z1, Z2 ]; 
% All two by one measurement vectors stacked on top of each other: 
T = Z'; Zstacked = T(:);

nMeasurements = length(Z1);

% Use the matlab minimization routine to find the optimal xhat: 
%
my_fOptions = foptions;
my_fOptions(1) = 0; 
xhat = fmins( 'sect_1_6_min_fn', [1;1], my_fOptions, [], Zstacked )

% Given this value of xhat what is the disribution of the noise:
%