对于此特定表达式（但不仅限于此表达式）
>> n(1)
ans = (sym) R₁⋅(I₁⋅(C₂⋅R₁⋅R₂⋅s + R₁ + R₂) + V₂)
我想要八度符号包的partfrac做什么
>> partfrac(n(1),I1)
ans = (sym)
⎛ 2 2 ⎞
I₁⋅⎝C₂⋅R₁ ⋅R₂⋅s + R₁ + R₁⋅R₂⎠ + R₁⋅V₂
>>
，但输出的矢量分为两部分，第一个索引取决于 I1 ，第二个索引取决于 V2 （独立于 I1 >）。
对于“驱动”变量 [I1 V2] 等的向量，在最佳情况下，分别获得特定表达式的向量和零。
最佳输出可能如下所示
>> xyz = somefunc(n(1),[I1 V2])
xyz = (sym 2×2 matrix)
⎡ I₁ V₂⎤
⎢ ⎥
⎢ 2 2 ⎥
⎣C₂⋅R₁ ⋅R₂⋅s + R₁ + R₁⋅R₂ R₁⎦
>>
基于嵌入sympy代码的解决方案也将受到赞赏。
非常感谢