I have a matrix (phi), say of 1502 rows (harmonic) by 10 columns (observations), where the i,j th element represents the phase of the ith harmonic of the jth observation; say the name of the matrix is X;
I need to find the modeshape corresponding to each harmonic; say y = sin(x + phi(i,j)), where x is a position vector, say x = 0:1:11 (11 position elements);
But since each of these modeshapes has its own corresponding 11 'y' values, what I need is a 3D matrix with size 1502x10x11 with all the modeshape values. How can I use the 2D matrix phi to get the 3D matrix (Modeshapes) that I want?
I tried the following code, which gives a better indication of what I need but it didn't work.
modeshapes = zeros(1502,10,11)
x = 0:1:11
for i = 1:1502
for j = 1:10
modeshapes(i,j,:) = sin(phi(i,j)+x)