The observation model will be ##z=H\theta+v##, where ##z## is column vector containing the results, ##H## is the observation matrix, ##\theta## is the parameter vector and ##v## is a vector containing random additive noise.

b.)
In this case ##z##=[1 1 0; 1 0 1; 0 1 1; 1 1 1]*[##m_1## ##m_2## ##m_3####]^{T}##+##v##.

Here's the code I came up with:

Code (Text):

H=[1 1 0;1 0 1; 0 1 1; 1 1 1];% Observation matrix
z=[4.45, 7.35 ,6.8, 10]'; % simulated weights of rock pairs, (I just made them up)
th=inv(H'*H)*H'*z; %solving for m1,m2,m3

It is easy to solve for the masses of the rocks using matrix operations. However, I am not sure how to model the two cases properly. How would one simulate the cases? What if the total weight of the rocks is unknown? Can you please help me with this problem?

P.S. I apologise in advance if this is wrong section for this topic.