Align to sensors based on the quaternion

Hi all

I am collecting quaternion data from two IMU sensors. I fixed the sensors together aiming to capture both sensors orientation with quaternion presentation. Apparently, the initial quaternion of the sensors may not be the same. So I need to find the rotation matrix between to sensors then continuously apply it to one the sensor. It was frequently mentioned that rotation correcting quaternion can be acquired using the following code.

Rotation = quatmultiply(quatinv(...
    Sensor1_Quaternion(:,InitialPosition)),...
    Sensor2_Quaternion(:,InitialPosition));
Sensor1_Quaternion_Rotatated = quatmultiply(...
    Sensor1_Quaternion,...
    Rotation);

but it did not work as I expected. Would you mind please help me to find out how I can fix the problem?

Best Regards, Reza