I've recently played with the IMU filter in MATLAB. When using their examples, they always plot the rotations by stating something alike this:

plot(t,eulerd(orientation,'ZYX','frame'))

( Note: orientation variable would be the collection of quaternion results )

As a result, I did the same, as I wanted to observe my roll, pitch, and yaw values. However, I ended up using the solution from this post. ( Which gave me the same results )

With that in mind, when I do a stack query for roll, yaw, and pitch from quaternions, I am almost always taken to a tait bryan angle representation. Does this mean that capturing these values is not possible through quaternions? Or is more difficult to visualize? Obviously my quaternion knowledge is small, as most of it is significantly over my head, but shouldn't I be able to calculate rotations and show them without relying on euler rotations?