help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: invfreqs -- specifying numerator without independant term


From: Doug Stewart
Subject: Re: invfreqs -- specifying numerator without independant term
Date: Wed, 22 Sep 2010 10:57:50 -0400

On Wed, Sep 22, 2010 at 10:17 AM, CdeMills <address@hidden> wrote:
>
> Hello,
>
> I'm trying to identify a frequency response of the form
>
> H(s) = (b2 s^2 + b1 s) / (1+a1 s)
>
> The specific problem is thus that the numerator has only two degrees of
> freedom. When I generate the frequency response using H(s) = freqs(B, A, w)
> then try to get back coefficients values through invfreqs, it fails as the
> equation system is singular. Is there some way to reformulate the problem to
> get around this missing term ?
>
> Regards
>
> Pascal
> --
> View this message in context: 
> http://octave.1599824.n4.nabble.com/invfreqs-specifying-numerator-without-independant-term-tp2550382p2550382.html
> Sent from the Octave - General mailing list archive at Nabble.com.
> _______________________________________________
> Help-octave mailing list
> address@hidden
> https://www-old.cae.wisc.edu/mailman/listinfo/help-octave
>



One idea that I use in control theory is to add a pole at a much
higher frequency. This makes the math work and doesn't change the
response at the low frequencies

try adding a pole say:

1000
_________
s + 1000

 I hope this helps.

Doug Stewart


reply via email to

[Prev in Thread] Current Thread [Next in Thread]