Jak uzyskać kształty trybu znormalizowanego masy?

Mam system MDOF o 4 stopniach swobody z wymuszonym wejściem na jednym z dof (4. dof) (jest to fala sinusoidalna o określonej częstotliwości, która ma być generowana za pomocą shakera); odpowiedź przesunięcia jest mierzona dla tego układu. Chciałbym eksperymentalnie wyodrębnić kształty trybów dla tego systemu, ale aby osiągnąć ten cel, podjąłem teoretyczny przypadek zrozumienia różnych komplikacji i problemów. Więc teraz mój system może być reprezentowany jako:

{ Y 1 Y 2 Y 3 Y 4 } = [ H 11 H 12 H 13 H 14 H 21 H 22 H 23 H 24 H 31 H 32 H 33 H 34 H 41 H 42 H 43 H

$$\begin{equation} G_{YY}=H G_{YX} \end{equation}$$ $$\begin{Bmatrix} Y_{1}\\ Y_{2}\\ Y_{3}\\ Y_{4} \end{Bmatrix}$$ { 0 0 0 X 4 $$\begin{bmatrix} H_{11} & H_{12} & H_{13} & H_{14}\\ H_{21} & H_{22} & H_{23} & H_{24}\\ H_{31} & H_{32} & H_{33} & H_{34}\\ H_{41} & H_{42} & H_{43} & H_{44}& \end{bmatrix}$$ $$\begin{Bmatrix} 0\\ 0\\ 0\\ X_{4} \end{Bmatrix}$$

Tutaj H jest funkcją odpowiedzi częstotliwościowej (FRF) i można zauważyć, że tylko jedna kolumna byłaby potrzebna do uzyskania kształtów modalnych systemu. Teraz, pisząc kod na Matlabie, daj mi:

> FRF(i,:)=GYY./GYX
> plot(fr,abs(real(GYY./GYX)))

$$\begin{equation} \phi^{T}\left[M_A\right]\phi=I \end{equation}$$

Oto, co powinienem uzyskać, co potwierdziła analiza własna macierzy masy M i sztywności K:

   -0.1442    0.3651    0.4152   -0.2710
   -0.2710    0.3651   -0.1442    0.4152
   -0.3651   -0.0000   -0.3651   -0.3651
   -0.4152   -0.3651    0.2710    0.1442

Ale z frf wyodrębniany jest kształt pierwszego trybu:

0.0098
0.0186
0.0254
0.0296

Wygląda porportalnie z rzeczywistą, ale nie jest skalowana do rzeczywistych wartości w celu spełnienia normalizacji masy. Gdzie się mylę?