Hi Stefani,
I understand that you looking for a sample code to resolve a system of your linear equations where all the variables are known except 
. The final equation you are trying to solve is 
.
. The final equation you are trying to solve is .For every j instead of calculating 
in 
for each f, it is possible to calculate 
(H x 1) at once using matrix operations in MATLAB. This in turn will reduce calculation at h and f level as it will be taken care internally using matrix operations. Simplify the given equation such that all the known variables are on right hand side of the equation as shown below.
in for each f, it is possible to calculate (H x 1) at once using matrix operations in MATLAB. This in turn will reduce calculation at h and f level as it will be taken care internally using matrix operations. Simplify the given equation such that all the known variables are on right hand side of the equation as shown below.
where 
and 
are (H x H) matrices and 
is a (H x 1) vector. The output 
will be a (H x 1) vector.
and are (H x H) matrices and is a (H x 1) vector. The output will be a (H x 1) vector.Refer to the pseudo code in the code snippet below to solve the equation for 
% Assuming that B_j, a_j, rho_j, ZX_j_rhoj variable are know and available,
% Xhat_j for every j calculated using the following code.
Xhat_j_rhoj_Plus1 = inv(a_j) * inv(B_j) * ((rho_j + 1) * (rho_j + 2)) * ZX_j_rhoj;
Xhat_j = Xhat_j_rhoj_Plus1.^(1/(rho_j +1)); % element wise power
Please refer to the following documentation link to learn more about Array and Matrix operations in MATLAB.
Hope this helps.
