Identify this type of matrix system of equation: $M_i = C_i + sum_{k=1}^K A_i M_k A_i^{prime}$
up vote
0
down vote
favorite
I'm trying to find more information on how to solve this linear system of equations with a funky form.
$M_i = C_i + sum_{k=1}^K pi_k A_i M_k A_i^{prime}$,
for all $i in {1,...,K}$ s.t. $M_i$ is symmetric and PSD, and $sum_{k=1}^K pi_k =1$. $M_i$ is a covariance matrix. $M_i$, $A_i$, and $C_i$ are all $(n times n)$ matrices, and $pi_i$ is a scalar.
I've found on my own that this may be a special type of a DARE equation but the Wikipedia article was hard to understand.
I've found a way to write the equation in block-diagonal form, but I still can't solve it analytically, due to the off-diagonal elements of the first-term after the equality being a function of $M_1$ and $M_2$. Here it is for the case where $K=2$.
$begin{bmatrix} M_1 & 0 \ 0 & M_2 end{bmatrix} = begin{bmatrix} C_1 & -pi_1 A_1M_1A_2^{prime} - pi_2 A_1 M_2 A_2^{prime} \ -pi_1 A_2M_1A_1^{prime} - pi_2 A_2 M_2 A_1^{prime} & C_2 end{bmatrix} + begin{bmatrix} pi_1 A_1 & pi_2 A_1 \ pi_1 A_2 & pi_2 A_2 end{bmatrix} begin{bmatrix} M_1 & 0 \ 0 & M_2 end{bmatrix} begin{bmatrix} A_1^{prime} & A_2^{prime} \ A_1^{prime} & A_2^{prime} end{bmatrix} $
statistics systems-of-equations matrix-equations
add a comment |
up vote
0
down vote
favorite
I'm trying to find more information on how to solve this linear system of equations with a funky form.
$M_i = C_i + sum_{k=1}^K pi_k A_i M_k A_i^{prime}$,
for all $i in {1,...,K}$ s.t. $M_i$ is symmetric and PSD, and $sum_{k=1}^K pi_k =1$. $M_i$ is a covariance matrix. $M_i$, $A_i$, and $C_i$ are all $(n times n)$ matrices, and $pi_i$ is a scalar.
I've found on my own that this may be a special type of a DARE equation but the Wikipedia article was hard to understand.
I've found a way to write the equation in block-diagonal form, but I still can't solve it analytically, due to the off-diagonal elements of the first-term after the equality being a function of $M_1$ and $M_2$. Here it is for the case where $K=2$.
$begin{bmatrix} M_1 & 0 \ 0 & M_2 end{bmatrix} = begin{bmatrix} C_1 & -pi_1 A_1M_1A_2^{prime} - pi_2 A_1 M_2 A_2^{prime} \ -pi_1 A_2M_1A_1^{prime} - pi_2 A_2 M_2 A_1^{prime} & C_2 end{bmatrix} + begin{bmatrix} pi_1 A_1 & pi_2 A_1 \ pi_1 A_2 & pi_2 A_2 end{bmatrix} begin{bmatrix} M_1 & 0 \ 0 & M_2 end{bmatrix} begin{bmatrix} A_1^{prime} & A_2^{prime} \ A_1^{prime} & A_2^{prime} end{bmatrix} $
statistics systems-of-equations matrix-equations
Which are your known variables and unknown variables?
– rzch
Dec 5 at 12:23
add a comment |
up vote
0
down vote
favorite
up vote
0
down vote
favorite
I'm trying to find more information on how to solve this linear system of equations with a funky form.
$M_i = C_i + sum_{k=1}^K pi_k A_i M_k A_i^{prime}$,
for all $i in {1,...,K}$ s.t. $M_i$ is symmetric and PSD, and $sum_{k=1}^K pi_k =1$. $M_i$ is a covariance matrix. $M_i$, $A_i$, and $C_i$ are all $(n times n)$ matrices, and $pi_i$ is a scalar.
I've found on my own that this may be a special type of a DARE equation but the Wikipedia article was hard to understand.
I've found a way to write the equation in block-diagonal form, but I still can't solve it analytically, due to the off-diagonal elements of the first-term after the equality being a function of $M_1$ and $M_2$. Here it is for the case where $K=2$.
$begin{bmatrix} M_1 & 0 \ 0 & M_2 end{bmatrix} = begin{bmatrix} C_1 & -pi_1 A_1M_1A_2^{prime} - pi_2 A_1 M_2 A_2^{prime} \ -pi_1 A_2M_1A_1^{prime} - pi_2 A_2 M_2 A_1^{prime} & C_2 end{bmatrix} + begin{bmatrix} pi_1 A_1 & pi_2 A_1 \ pi_1 A_2 & pi_2 A_2 end{bmatrix} begin{bmatrix} M_1 & 0 \ 0 & M_2 end{bmatrix} begin{bmatrix} A_1^{prime} & A_2^{prime} \ A_1^{prime} & A_2^{prime} end{bmatrix} $
statistics systems-of-equations matrix-equations
I'm trying to find more information on how to solve this linear system of equations with a funky form.
$M_i = C_i + sum_{k=1}^K pi_k A_i M_k A_i^{prime}$,
for all $i in {1,...,K}$ s.t. $M_i$ is symmetric and PSD, and $sum_{k=1}^K pi_k =1$. $M_i$ is a covariance matrix. $M_i$, $A_i$, and $C_i$ are all $(n times n)$ matrices, and $pi_i$ is a scalar.
I've found on my own that this may be a special type of a DARE equation but the Wikipedia article was hard to understand.
I've found a way to write the equation in block-diagonal form, but I still can't solve it analytically, due to the off-diagonal elements of the first-term after the equality being a function of $M_1$ and $M_2$. Here it is for the case where $K=2$.
$begin{bmatrix} M_1 & 0 \ 0 & M_2 end{bmatrix} = begin{bmatrix} C_1 & -pi_1 A_1M_1A_2^{prime} - pi_2 A_1 M_2 A_2^{prime} \ -pi_1 A_2M_1A_1^{prime} - pi_2 A_2 M_2 A_1^{prime} & C_2 end{bmatrix} + begin{bmatrix} pi_1 A_1 & pi_2 A_1 \ pi_1 A_2 & pi_2 A_2 end{bmatrix} begin{bmatrix} M_1 & 0 \ 0 & M_2 end{bmatrix} begin{bmatrix} A_1^{prime} & A_2^{prime} \ A_1^{prime} & A_2^{prime} end{bmatrix} $
statistics systems-of-equations matrix-equations
statistics systems-of-equations matrix-equations
asked Dec 5 at 1:33
hipHopMetropolisHastings
1718
1718
Which are your known variables and unknown variables?
– rzch
Dec 5 at 12:23
add a comment |
Which are your known variables and unknown variables?
– rzch
Dec 5 at 12:23
Which are your known variables and unknown variables?
– rzch
Dec 5 at 12:23
Which are your known variables and unknown variables?
– rzch
Dec 5 at 12:23
add a comment |
active
oldest
votes
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',
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%2f3026473%2fidentify-this-type-of-matrix-system-of-equation-m-i-c-i-sum-k-1k-a-i-m%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
active
oldest
votes
active
oldest
votes
active
oldest
votes
active
oldest
votes
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.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- 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.
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%2f3026473%2fidentify-this-type-of-matrix-system-of-equation-m-i-c-i-sum-k-1k-a-i-m%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
Which are your known variables and unknown variables?
– rzch
Dec 5 at 12:23