I have two subsystems which the transfer functions you could see below as g(s) and h(s). According to some references show below, I found that the bode plot of a cascade system should be the sum of subsystems. Then, to comfirm that I implemented code as below in Matlab, but the result seems violate with the reference. So, where I went wrong?
Subsystem transfer functions:
Matlab code:
den = [0.3352 0 209.8621 0];
gs = tf(num, den);
num = [1 3];
den = [21 2 10 40];
hs = tf(num, den);
c = gs*hs;
figure
bode(hs,gs,c)
grid on
legend('subsystem1','subsystem2','cascade system')
Result:



