You use
for W0 = rand (m,r);
for is defined as proceeding along the columns of whatever is on the right hand side, so W0 will be set to columns of length m x 1, and that will happen r times.
However, it turns out that you never use W0 in your code, so the overall effect is as if you had just set up a loop to happen r times.
On the next line you have
H0 = h0.* (V*w0')./((w0'*w0)*h0 + (10^-9));
This makes use of a bunch of variables that we as readers do not know the definition of. Your variable w0 there is not the same as W0 because MATLAB is case sensitive.
We know that V is either 2 or 3 dimensional, and is an integer data type, as it was read in by imread(). You use it with a * operation with w0 . We can tell by context that w0 is intended to be a vector or a 2D array. If w0 is not a scalar, then you have a integer data class * a vector or array. But the * operation is not fully implemented for integer data classes: if you have an integer data class, the only permitted operations are:
- integer scalar * integer scalar
- integer scalar * integer array
- integer scalar * double scalar
- integer scalar * double array
- integer array * integer scalar
- integer array * double scalar
- double scalar * integer scalar
- double scalar * integer array
- double array * integer scalar
or to put it another way, when one side of the * is an integer array (non-scalar), then the other side must be a scalar integer or scalar double, and it is not permitted to have a non-scalar integer array on one side and a non-scalar array of any kind on the other side.
The easiest fix:
V = double( imread('circuit.tif') );
