result = ZerosAndPoles.Analysis.evaluate(zp, p, den_min=0)
Function Analysis.evaluate evaluates the ZerosAndPoles transfer function at a given (complex) value of p and returns the value G(p)=N(p)/D(p). The optional argument den_min with default 0 is used to guard against a division by zero.
if |(D(p))| >= den_min then G(p) = N(p) / D(p); elseif D(p).re >= 0.0 then G(p) = N(p) / den_min else G(p) = -N(p) / den_min end if;
import Modelica.ComplexMath.j; ZerosAndPoles p = Modelica_LinearSystems2.ZerosAndPoles.p(); Modelica_LinearSystems2.ZerosAndPoles zp=(p+1)/(p^2+p+1); Complex result; algorithm result := Modelica_LinearSystems2.ZerosAndPoles.Analysis.evaluate(zp, j+1); // result = 0.538462 - 0.307692*j
Math.Polynomial.evaluateComplex
encapsulated function frequencyEvaluate extends Modelica.Icons.Function; import Modelica; import Modelica_LinearSystems2.Internal; input Real gain "Gain of transfer function"; input Real Zeros[:, 2] "Zeros as Real matrix (first column: real, second column imaginary values)"; input Real Poles[:, 2] "Poles as Real matrix (first column: real, second column imaginary values)"; input Real re "Real value of complex number"; input Real im "Imaginary value of complex number"; output Real A "Amplitude of complex value"; output Modelica.Units.SI.Angle phi "Angle of complex value"; output Integer info "= 0/1/2 success/infinity(A is a large value)/indefinite (A=0/0; A=0 returned)"; end frequencyEvaluate;
Date | Author | Comment |
---|---|---|
2010-05-31 | Marcus Baur, DLR-RM | Realization |