Re: zp2sos error: cplxpair: could not pair all complex numbers

[Nicholas Jankowski]

On Thu, Jan 5, 2017 at 12:30 PM, Nicholas Jankowski <address@hidden> wrote:
>> If you can help identify what Matlab does in this situation, and what
>> should be fixed and where, then maybe something can be done about it.
>
> Running the same script in Matlab 2016b produces no errors and the
> following output.
> Note Octave output from cheby2 is a row vector, Matlab is a column. I
> don't know if this matters. Matlab output follows. Octave 4.2.0,
> signal pkg v.1.3.2 ouput shown below that.

Oops, I left off the important part of my scripts as well.

-----------------------
Matlab 2016b:
>>  [sos, g] = zp2sos (z, p, k)
sos =
  Columns 1 through 5
   1.000000000000000   1.999811353006249   0.999999999999999
1.000000000000000   1.899456851014267
   1.000000000000000   1.998592402557042   1.000000000000000
1.000000000000000   1.923500197249950
   1.000000000000000   1.997374287874268   0.999999999999998
1.000000000000000   1.967506243219755
   1.000000000000000  -1.999999259461443   0.999999999999999
1.000000000000000  -1.993570273311965
   1.000000000000000  -1.999994472762730   1.000000000000001
1.000000000000000  -1.995281755784758
   1.000000000000000  -1.999989686441154   1.000000000000001
1.000000000000000  -1.998255341271938
  Column 6
   0.902325817320760
   0.927580075649315
   0.972881387742377
   0.993582084086443
   0.995298350815657
   0.998276741178964
g =
   0.896584733013338

-----------
Octave 4.2.0:

>> [sos, g] = zp2sos (z, p, k)
error: cplxpair: could not pair all complex numbers
error: called from
    cplxpair at line 118 column 9
    cplxreal at line 52 column 9
    zp2sos at line 77 column 10

%% just to verify it's not as simple as the orientation of z and p:

>> [sos, g] = zp2sos (z', p', k)
error: cplxpair: could not pair all complex numbers
error: called from
    cplxpair at line 118 column 9
    cplxreal at line 52 column 9
    zp2sos at line 77 column 10