Filter
Top Products
fitsys 241
real imaginary frequency damping
(rad/sec) ratio
-1.0000e+00 0.0000e+00 1.0000e+00 1.0000
5.0000e-01 3.1225e+00 3.1623e+00 -0.1581
5.0000e-01 -3.1225e+00 3.1623e+00 -0.1581
omega = logspace(0.001,100,200)
plantg = freq(plant,omega)
# Use complex cepstrum to fit minimum phase equivalent
# to the magnitude of the data. One of the principle
# uses of the fitsys function is fitting approximations
# to noisy data. To illustrate the concepts, no noise
# is added here.
cdata = mkphase(abs(plantg))
gph1 = ctrlplot([plantg,cdata],{phaseplot});
gph1 = plot(gph1,{title="Data and minimum phase fit",...
legend=["original data";"minimum phase fit"]})?