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} $










share|cite|improve this question






















  • Which are your known variables and unknown variables?
    – rzch
    Dec 5 at 12:23















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} $










share|cite|improve this question






















  • Which are your known variables and unknown variables?
    – rzch
    Dec 5 at 12:23













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} $










share|cite|improve this question













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






share|cite|improve this question













share|cite|improve this question











share|cite|improve this question




share|cite|improve this question










asked Dec 5 at 1:33









hipHopMetropolisHastings

1718




1718












  • 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




Which are your known variables and unknown variables?
– rzch
Dec 5 at 12:23















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
});


}
});














draft saved

draft discarded


















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
















draft saved

draft discarded




















































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.




draft saved


draft discarded














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





















































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







Popular posts from this blog

Bressuire

Cabo Verde

Gyllenstierna