How to decompose a MIMO system into a multiplication minimum and nonminimum phase part
$begingroup$
suppose the transfer function of m-input-m-output system is
$$Gleft( s right) = left[ {begin{array}{*{20}{c}}
{{g_{11}}left( s right)}& cdots &{{g_{1m}}left( s right)}\
vdots & ddots & vdots \
{{g_{m1}}left( s right)}& cdots &{{g_{mm}}left( s right)}
end{array}} right]$$
Each ${g_{ij}}(s)$ might contains nonminimum phase part.Then, how could the system be decomposed into the multiplication of minimum ${G_{min }}left( s right)$ and nonminimum phase part ${G_{nonmin }}left( s right)$
$$Gleft( s right) = {G_{min }}left( s right) * {G_{nonmin }}left( s right)$$
Is there any command in MATLAB to realize it?
matlab control-theory
$endgroup$
add a comment |
$begingroup$
suppose the transfer function of m-input-m-output system is
$$Gleft( s right) = left[ {begin{array}{*{20}{c}}
{{g_{11}}left( s right)}& cdots &{{g_{1m}}left( s right)}\
vdots & ddots & vdots \
{{g_{m1}}left( s right)}& cdots &{{g_{mm}}left( s right)}
end{array}} right]$$
Each ${g_{ij}}(s)$ might contains nonminimum phase part.Then, how could the system be decomposed into the multiplication of minimum ${G_{min }}left( s right)$ and nonminimum phase part ${G_{nonmin }}left( s right)$
$$Gleft( s right) = {G_{min }}left( s right) * {G_{nonmin }}left( s right)$$
Is there any command in MATLAB to realize it?
matlab control-theory
$endgroup$
$begingroup$
There is no unique way of doing this, it will only tell you how to split the zeros, but not the poles.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:09
$begingroup$
If $m=1$, then it is easy to decompose by simply finding the unstable zeros. However, for MIMO system, even if I can find unstable zeros for each $g_{ij}(s)$, it's still hard to obtain the final multiplication form.
$endgroup$
– user4756766
Jan 18 at 12:14
$begingroup$
So you do not care if one of the final transfer functions is not causal, otherwise you still also have to choose how to split up the poles, even for $m=1$.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:16
$begingroup$
What if I do not care about whether it's causal. Is there any simple method to realize the decomposition? Could you give me some hint about that @KwinvanderVeen
$endgroup$
– user4756766
Jan 18 at 12:22
add a comment |
$begingroup$
suppose the transfer function of m-input-m-output system is
$$Gleft( s right) = left[ {begin{array}{*{20}{c}}
{{g_{11}}left( s right)}& cdots &{{g_{1m}}left( s right)}\
vdots & ddots & vdots \
{{g_{m1}}left( s right)}& cdots &{{g_{mm}}left( s right)}
end{array}} right]$$
Each ${g_{ij}}(s)$ might contains nonminimum phase part.Then, how could the system be decomposed into the multiplication of minimum ${G_{min }}left( s right)$ and nonminimum phase part ${G_{nonmin }}left( s right)$
$$Gleft( s right) = {G_{min }}left( s right) * {G_{nonmin }}left( s right)$$
Is there any command in MATLAB to realize it?
matlab control-theory
$endgroup$
suppose the transfer function of m-input-m-output system is
$$Gleft( s right) = left[ {begin{array}{*{20}{c}}
{{g_{11}}left( s right)}& cdots &{{g_{1m}}left( s right)}\
vdots & ddots & vdots \
{{g_{m1}}left( s right)}& cdots &{{g_{mm}}left( s right)}
end{array}} right]$$
Each ${g_{ij}}(s)$ might contains nonminimum phase part.Then, how could the system be decomposed into the multiplication of minimum ${G_{min }}left( s right)$ and nonminimum phase part ${G_{nonmin }}left( s right)$
$$Gleft( s right) = {G_{min }}left( s right) * {G_{nonmin }}left( s right)$$
Is there any command in MATLAB to realize it?
matlab control-theory
matlab control-theory
asked Jan 18 at 12:00
user4756766user4756766
357
357
$begingroup$
There is no unique way of doing this, it will only tell you how to split the zeros, but not the poles.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:09
$begingroup$
If $m=1$, then it is easy to decompose by simply finding the unstable zeros. However, for MIMO system, even if I can find unstable zeros for each $g_{ij}(s)$, it's still hard to obtain the final multiplication form.
$endgroup$
– user4756766
Jan 18 at 12:14
$begingroup$
So you do not care if one of the final transfer functions is not causal, otherwise you still also have to choose how to split up the poles, even for $m=1$.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:16
$begingroup$
What if I do not care about whether it's causal. Is there any simple method to realize the decomposition? Could you give me some hint about that @KwinvanderVeen
$endgroup$
– user4756766
Jan 18 at 12:22
add a comment |
$begingroup$
There is no unique way of doing this, it will only tell you how to split the zeros, but not the poles.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:09
$begingroup$
If $m=1$, then it is easy to decompose by simply finding the unstable zeros. However, for MIMO system, even if I can find unstable zeros for each $g_{ij}(s)$, it's still hard to obtain the final multiplication form.
$endgroup$
– user4756766
Jan 18 at 12:14
$begingroup$
So you do not care if one of the final transfer functions is not causal, otherwise you still also have to choose how to split up the poles, even for $m=1$.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:16
$begingroup$
What if I do not care about whether it's causal. Is there any simple method to realize the decomposition? Could you give me some hint about that @KwinvanderVeen
$endgroup$
– user4756766
Jan 18 at 12:22
$begingroup$
There is no unique way of doing this, it will only tell you how to split the zeros, but not the poles.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:09
$begingroup$
There is no unique way of doing this, it will only tell you how to split the zeros, but not the poles.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:09
$begingroup$
If $m=1$, then it is easy to decompose by simply finding the unstable zeros. However, for MIMO system, even if I can find unstable zeros for each $g_{ij}(s)$, it's still hard to obtain the final multiplication form.
$endgroup$
– user4756766
Jan 18 at 12:14
$begingroup$
If $m=1$, then it is easy to decompose by simply finding the unstable zeros. However, for MIMO system, even if I can find unstable zeros for each $g_{ij}(s)$, it's still hard to obtain the final multiplication form.
$endgroup$
– user4756766
Jan 18 at 12:14
$begingroup$
So you do not care if one of the final transfer functions is not causal, otherwise you still also have to choose how to split up the poles, even for $m=1$.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:16
$begingroup$
So you do not care if one of the final transfer functions is not causal, otherwise you still also have to choose how to split up the poles, even for $m=1$.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:16
$begingroup$
What if I do not care about whether it's causal. Is there any simple method to realize the decomposition? Could you give me some hint about that @KwinvanderVeen
$endgroup$
– user4756766
Jan 18 at 12:22
$begingroup$
What if I do not care about whether it's causal. Is there any simple method to realize the decomposition? Could you give me some hint about that @KwinvanderVeen
$endgroup$
– user4756766
Jan 18 at 12:22
add a comment |
1 Answer
1
active
oldest
votes
$begingroup$
If you define the zeros of a transfer function as the values for $s$ at which it loses rank (also known as transmission zeros) then you could make use of the Smith McMillan form. Namely by multiplying any MIMO transfer function matrix by its pole polynomial one can obtain a polynomial matrix $P(s)$. Any polynomial can be written as
$$
P(s) = U_1(s),D(s),U_2(s)
$$
where $U_1(s)$ and $U_2(s)$ are unimodular (does not lose rank for any $s$) and $D(s)$ is diagonal (does loses rank at $s$ equal to a zero). This diagonal matrix can be split up into two diagonal matrices, one with only right half plane zeros and the other with only left half plane zeros.
For example
$$
G(s) =
begin{bmatrix}
frac{s^4+s+1}{(s^2+1)(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+3)} \
frac{s-1}{(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+1)(s+3)}
end{bmatrix}
$$
can also be written as
begin{align}
G(s) &= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix} \
&= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & 1
end{bmatrix}
begin{bmatrix}
1 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix}
end{align}
Multiplying the left unimodular matrix with the left diagonal matrix and the right diagonal matrix with right unimodular matrix and distribute the terms of the pole polynomial gives
$$
G(s) =
begin{bmatrix}
frac{s+2}{s+1} & 1 \
0 & frac{1}{s+1}
end{bmatrix}
begin{bmatrix}
frac{1}{(s^2+1)(s+3)} & 0 \
frac{s-1}{s+3} & frac{1-s}{(s^2+1)(s+3)}
end{bmatrix}
$$
Of course you could use a different distribution of the pole polynomial and you could throw any two unimodular matrices, which are each others inverse, in-between the two two diagonal matrices. So splitting a transfer function matrix into a minimum phase and non-minimum phase part will in general not be unique.
$endgroup$
add a comment |
Your Answer
StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3078152%2fhow-to-decompose-a-mimo-system-into-a-multiplication-minimum-and-nonminimum-phas%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
$begingroup$
If you define the zeros of a transfer function as the values for $s$ at which it loses rank (also known as transmission zeros) then you could make use of the Smith McMillan form. Namely by multiplying any MIMO transfer function matrix by its pole polynomial one can obtain a polynomial matrix $P(s)$. Any polynomial can be written as
$$
P(s) = U_1(s),D(s),U_2(s)
$$
where $U_1(s)$ and $U_2(s)$ are unimodular (does not lose rank for any $s$) and $D(s)$ is diagonal (does loses rank at $s$ equal to a zero). This diagonal matrix can be split up into two diagonal matrices, one with only right half plane zeros and the other with only left half plane zeros.
For example
$$
G(s) =
begin{bmatrix}
frac{s^4+s+1}{(s^2+1)(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+3)} \
frac{s-1}{(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+1)(s+3)}
end{bmatrix}
$$
can also be written as
begin{align}
G(s) &= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix} \
&= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & 1
end{bmatrix}
begin{bmatrix}
1 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix}
end{align}
Multiplying the left unimodular matrix with the left diagonal matrix and the right diagonal matrix with right unimodular matrix and distribute the terms of the pole polynomial gives
$$
G(s) =
begin{bmatrix}
frac{s+2}{s+1} & 1 \
0 & frac{1}{s+1}
end{bmatrix}
begin{bmatrix}
frac{1}{(s^2+1)(s+3)} & 0 \
frac{s-1}{s+3} & frac{1-s}{(s^2+1)(s+3)}
end{bmatrix}
$$
Of course you could use a different distribution of the pole polynomial and you could throw any two unimodular matrices, which are each others inverse, in-between the two two diagonal matrices. So splitting a transfer function matrix into a minimum phase and non-minimum phase part will in general not be unique.
$endgroup$
add a comment |
$begingroup$
If you define the zeros of a transfer function as the values for $s$ at which it loses rank (also known as transmission zeros) then you could make use of the Smith McMillan form. Namely by multiplying any MIMO transfer function matrix by its pole polynomial one can obtain a polynomial matrix $P(s)$. Any polynomial can be written as
$$
P(s) = U_1(s),D(s),U_2(s)
$$
where $U_1(s)$ and $U_2(s)$ are unimodular (does not lose rank for any $s$) and $D(s)$ is diagonal (does loses rank at $s$ equal to a zero). This diagonal matrix can be split up into two diagonal matrices, one with only right half plane zeros and the other with only left half plane zeros.
For example
$$
G(s) =
begin{bmatrix}
frac{s^4+s+1}{(s^2+1)(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+3)} \
frac{s-1}{(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+1)(s+3)}
end{bmatrix}
$$
can also be written as
begin{align}
G(s) &= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix} \
&= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & 1
end{bmatrix}
begin{bmatrix}
1 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix}
end{align}
Multiplying the left unimodular matrix with the left diagonal matrix and the right diagonal matrix with right unimodular matrix and distribute the terms of the pole polynomial gives
$$
G(s) =
begin{bmatrix}
frac{s+2}{s+1} & 1 \
0 & frac{1}{s+1}
end{bmatrix}
begin{bmatrix}
frac{1}{(s^2+1)(s+3)} & 0 \
frac{s-1}{s+3} & frac{1-s}{(s^2+1)(s+3)}
end{bmatrix}
$$
Of course you could use a different distribution of the pole polynomial and you could throw any two unimodular matrices, which are each others inverse, in-between the two two diagonal matrices. So splitting a transfer function matrix into a minimum phase and non-minimum phase part will in general not be unique.
$endgroup$
add a comment |
$begingroup$
If you define the zeros of a transfer function as the values for $s$ at which it loses rank (also known as transmission zeros) then you could make use of the Smith McMillan form. Namely by multiplying any MIMO transfer function matrix by its pole polynomial one can obtain a polynomial matrix $P(s)$. Any polynomial can be written as
$$
P(s) = U_1(s),D(s),U_2(s)
$$
where $U_1(s)$ and $U_2(s)$ are unimodular (does not lose rank for any $s$) and $D(s)$ is diagonal (does loses rank at $s$ equal to a zero). This diagonal matrix can be split up into two diagonal matrices, one with only right half plane zeros and the other with only left half plane zeros.
For example
$$
G(s) =
begin{bmatrix}
frac{s^4+s+1}{(s^2+1)(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+3)} \
frac{s-1}{(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+1)(s+3)}
end{bmatrix}
$$
can also be written as
begin{align}
G(s) &= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix} \
&= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & 1
end{bmatrix}
begin{bmatrix}
1 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix}
end{align}
Multiplying the left unimodular matrix with the left diagonal matrix and the right diagonal matrix with right unimodular matrix and distribute the terms of the pole polynomial gives
$$
G(s) =
begin{bmatrix}
frac{s+2}{s+1} & 1 \
0 & frac{1}{s+1}
end{bmatrix}
begin{bmatrix}
frac{1}{(s^2+1)(s+3)} & 0 \
frac{s-1}{s+3} & frac{1-s}{(s^2+1)(s+3)}
end{bmatrix}
$$
Of course you could use a different distribution of the pole polynomial and you could throw any two unimodular matrices, which are each others inverse, in-between the two two diagonal matrices. So splitting a transfer function matrix into a minimum phase and non-minimum phase part will in general not be unique.
$endgroup$
If you define the zeros of a transfer function as the values for $s$ at which it loses rank (also known as transmission zeros) then you could make use of the Smith McMillan form. Namely by multiplying any MIMO transfer function matrix by its pole polynomial one can obtain a polynomial matrix $P(s)$. Any polynomial can be written as
$$
P(s) = U_1(s),D(s),U_2(s)
$$
where $U_1(s)$ and $U_2(s)$ are unimodular (does not lose rank for any $s$) and $D(s)$ is diagonal (does loses rank at $s$ equal to a zero). This diagonal matrix can be split up into two diagonal matrices, one with only right half plane zeros and the other with only left half plane zeros.
For example
$$
G(s) =
begin{bmatrix}
frac{s^4+s+1}{(s^2+1)(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+3)} \
frac{s-1}{(s+1)(s+3)} & frac{1-s}{(s^2+1)(s+1)(s+3)}
end{bmatrix}
$$
can also be written as
begin{align}
G(s) &= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix} \
&= frac{1}{(s^2+1)(s+1)(s+3)}
begin{bmatrix}
1 & s+1 \
0 & 1
end{bmatrix}
begin{bmatrix}
s+2 & 0 \
0 & 1
end{bmatrix}
begin{bmatrix}
1 & 0 \
0 & s-1
end{bmatrix}
begin{bmatrix}
1 & 0 \
s^2+1 & -1
end{bmatrix}
end{align}
Multiplying the left unimodular matrix with the left diagonal matrix and the right diagonal matrix with right unimodular matrix and distribute the terms of the pole polynomial gives
$$
G(s) =
begin{bmatrix}
frac{s+2}{s+1} & 1 \
0 & frac{1}{s+1}
end{bmatrix}
begin{bmatrix}
frac{1}{(s^2+1)(s+3)} & 0 \
frac{s-1}{s+3} & frac{1-s}{(s^2+1)(s+3)}
end{bmatrix}
$$
Of course you could use a different distribution of the pole polynomial and you could throw any two unimodular matrices, which are each others inverse, in-between the two two diagonal matrices. So splitting a transfer function matrix into a minimum phase and non-minimum phase part will in general not be unique.
answered Jan 18 at 16:30
Kwin van der VeenKwin van der Veen
5,4952828
5,4952828
add a comment |
add a comment |
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3078152%2fhow-to-decompose-a-mimo-system-into-a-multiplication-minimum-and-nonminimum-phas%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
$begingroup$
There is no unique way of doing this, it will only tell you how to split the zeros, but not the poles.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:09
$begingroup$
If $m=1$, then it is easy to decompose by simply finding the unstable zeros. However, for MIMO system, even if I can find unstable zeros for each $g_{ij}(s)$, it's still hard to obtain the final multiplication form.
$endgroup$
– user4756766
Jan 18 at 12:14
$begingroup$
So you do not care if one of the final transfer functions is not causal, otherwise you still also have to choose how to split up the poles, even for $m=1$.
$endgroup$
– Kwin van der Veen
Jan 18 at 12:16
$begingroup$
What if I do not care about whether it's causal. Is there any simple method to realize the decomposition? Could you give me some hint about that @KwinvanderVeen
$endgroup$
– user4756766
Jan 18 at 12:22