假定m是一个数，你可以把-PBR^(-1)B'P+mPDD'P合并成-P(BR^(-1)B'-mDD')P，也就是把B看成I，R看成(BR^(-1)B'-mDD')^{-1}，这样解虽然不见得稳定但至少可以利用一下现成的care
