158 CHAPTER 4. DEMONSTRATION EXAMPLES
Wdist,Wnoise,Wperf,Wactp,Wactv,piezo)
size(P)?
ans (a row vector) = 7 6 8
Now select a frequency grid for calculating the frequency responses. Some additional
points are included near the oscillatory modes.
omega = logspace(1,1000,15)’
omega = sort([omega; [5.275:0.05:5.625]’; [4.5:0.25:6.5]’])’
Examine the frequency response of the open-loop system.
vcmodelg = freq(vcmodel,omega)
gph1 = ctrlplot(vcmodelg,bode);
gph1 = plot(gph1,{title="voice coil model",!grid})?