むだ時間モデルの近似値を求める
tfm = TransferFunctionModel[{{{-20.}}, 30. - 20./E^s + 0.1*s + s^2}, s, SamplingPeriod -> None];delaypoles = Flatten@TransferFunctionPoles[tfm, {{-10, 10}, {-50, 50}}];
ListPlot[Table[{Re[p], Im[p]}, {p, delaypoles}], PlotMarkers -> {Automatic, 8}]SystemsModelDelayApproximateを使って,むだ時間系の遅延のない近似値を求めることができる.次数が増えるにつれて,近似値はもとのモデルに収束する:
Manipulate[With[{approxtfm = SystemsModelDelayApproximate[tfm, order]}, Row[{BodePlot[{tfm, approxtfm}, Ticks -> None, Frame -> False, ImageSize -> 250][[1, 1, 1]], ListPlot[{Table[{Re[p], Im[p]}, {p, delaypoles}], Table[{Re[p], Im[p]}, {p, Flatten[TransferFunctionPoles[approxtfm]]}]}, PlotMarkers -> {Automatic, 8}, Ticks -> None, ImageSize -> 250, PlotRange -> {{-10, 10}, {-50, 50}}]}]], {order, 0, 15, 1}, SaveDefinitions -> True]