我有一个传递函数系数如下：
[b,a] = ellip(4,.2,40,[.41 .47]);
我想将此系统纳入其 PFE，使用[r,p,k] = residue(b,a). 然后，我想实现一个过滤器，它使用二阶子部分的并行组合。为此，我需要将每个复极点与其复共轭组合，以便整个二阶子部分具有实值系数。因此; 我选择了四对一阶项，并使用“残差”函数将它们重新组合成二阶小节。如何将生成的二阶子部分的系数存储在矩阵 c 和 d 中，以便每一行对应一个二阶系统？
我试图这样做，但我不知道它是否正确？
[c1 d1]=residue([r(1) r(2)],[p(1) p(2)],0);
[c2 d2]=residue([r(3) r(4)],[p(3) p(4)],0);
[c3 d3]=residue([r(5) r(6)],[p(5) p(6)],0);
[c4 d4]=residue([r(7) r(8)],[p(7) p(8)],k);
[c d]=[c1 c2 c3 c4 d1 d2 d3 d4]
有人可以给我一个想法吗？任何帮助，将不胜感激