randsys 331
real imaginary frequency damping
(rad/sec) ratio
-1.0408e+00 0.0000e+00 1.0408e+00 1.0000
-1.1390e+00 0.0000e+00 1.1390e+00 1.0000
-2.8050e-01 -1.4858e+00 1.5120e+00 0.1855
-2.8050e-01 1.4858e+00 1.5120e+00 0.1855
-1.5337e+00 0.0000e+00 1.5337e+00 1.0000
-5.5566e-01 -1.6958e+00 1.7845e+00 0.3114
-5.5566e-01 1.6958e+00 1.7845e+00 0.3114
-3.1454e+00 0.0000e+00 3.1454e+00 1.0000
-5.6798e+00 0.0000e+00 5.6798e+00 1.0000
-6.1039e+00 0.0000e+00 6.1039e+00 1.0000
Zeros:
real imaginary frequency damping
(rad/sec) ratio
-8.6788e-02 -1.0161e+00 1.0198e+00 0.0851
-8.6788e-02 1.0161e+00 1.0198e+00 0.0851
-1.0395e+00 0.0000e+00 1.0395e+00 1.0000
-1.5152e+00 0.0000e+00 1.5152e+00 1.0000
-6.3796e-01 -1.9633e+00 2.0643e+00 0.3090
-6.3796e-01 1.9633e+00 2.0643e+00 0.3090
-3.1425e+00 0.0000e+00 3.1425e+00 1.0000
-5.8916e+00 1.5472e-01 5.8937e+00 0.9997
-5.8916e+00 -1.5472e-01 5.8937e+00 0.9997
fHz = logspace(0.01,100,100)
sys1g = freq(sys1,fHz)
sys2g = freq(sys2,fHz)
gph1 = ctrlplot([sys1g,sys2g],{bode});
gph1 = plot(gph1,{legend=["sys1 (input 1)";...
"sys1 (input 2)";"sys2"]})?