normal rank calculation with tzero

22 次查看(过去 30 天)
Consider the following transfer matrix:
s = tf('s')
s = s Continuous-time transfer function.
G = [1/(s+2) 0; 0 1/((s+1)*(s+3))]
G = From input 1 to output... 1 1: ----- s + 2 2: 0 From input 2 to output... 1: 0 1 2: ------------- s^2 + 4 s + 3 Continuous-time transfer function.
the calcultation of the normal rank of G is done by:
[z, nrank] = tzero(G)
z = 0x1 empty double column vector
nrank = 1
and the result is
nrank = 1.
I was wondering, because in my understanding G will never lose rank for any and the normal rank should be equal to 2. So far, my observation is, that if the number of poles in the part transfer functions is equal,
i.e.
G = [1/(s+2) 0; 0 1/(s+1)]
G = From input 1 to output... 1 1: ----- s + 2 2: 0 From input 2 to output... 1: 0 1 2: ----- s + 1 Continuous-time transfer function.
[z, nrank] = tzero(G)
z = 0x1 empty double column vector
nrank = 2
the result meets with my expectations (nrank = 2).
So my question is, am I missing something when using tzero or is there a general problem of understanding regarding the normal rank of a transfer matrix?
Many thanks in advance.
  2 个评论
Paul
Paul 2025-3-27
Hi Andreas,
I took a look and I agree that the first case seems odd. I tried lots of variation in the tol input to tzero, but that had no effect. Unfortunately, the heavy lifting for tzero, at least for the first case, is in a .mex file, so I couldn't dig into the code.
If you open a case with Tech Support, would you mind posting back here with a summary of their response?
Andreas
Andreas 2025-3-28
编辑:Andreas 2025-4-2
Hi Paul,
thank you very much for your answer. sure, I'll do that.

请先登录,再进行评论。

采纳的回答

Christian
Christian 2025-4-7
We have identified that this is indeed a bug in "tzero".
We are actively working on a fix, and aim to provide it as soon as possible.

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Stability Analysis 的更多信息

产品


版本

R2024b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by