Skip to content

Instantly share code, notes, and snippets.

@larsoner
Created January 3, 2014 15:03
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save larsoner/8239281 to your computer and use it in GitHub Desktop.
Save larsoner/8239281 to your computer and use it in GitHub Desktop.
SOS tests for scipy
========================== BESSELF =========================
z =
Empty matrix: 0-by-1
p =
-3.000000000000000e-01
k =
3.000000000000000e-01
z =
Empty matrix: 0-by-1
p =
-3.893784325301967e-01 + 2.082923181459433e-02i
-3.893784325301967e-01 - 2.082923181459433e-02i
-3.862735176153423e-01 + 6.252441204484148e-02i
-3.862735176153423e-01 - 6.252441204484148e-02i
-3.800063969264777e-01 + 1.043324153265292e-01i
-3.800063969264777e-01 - 1.043324153265292e-01i
-3.704569665070659e-01 + 1.463376911059390e-01i
-3.704569665070659e-01 - 1.463376911059390e-01i
-3.574300381579267e-01 + 1.886403949711383e-01i
-3.574300381579267e-01 - 1.886403949711383e-01i
-3.406329048958464e-01 + 2.313670251522233e-01i
-3.406329048958464e-01 - 2.313670251522233e-01i
-3.196358682519445e-01 + 2.746876212958544e-01i
-3.196358682519445e-01 - 2.746876212958544e-01i
-2.938003295524644e-01 + 3.188464058882734e-01i
-2.938003295524644e-01 - 3.188464058882734e-01i
-2.621375273973739e-01 + 3.642225746223757e-01i
-2.621375273973739e-01 - 3.642225746223757e-01i
-2.229943267172552e-01 + 4.114690805661358e-01i
-2.229943267172552e-01 - 4.114690805661358e-01i
-1.731977157734973e-01 + 4.619039345030320e-01i
-1.731977157734973e-01 - 4.619039345030320e-01i
-1.046396975135688e-01 + 5.191384340047482e-01i
-1.046396975135688e-01 - 5.191384340047482e-01i
k =
1.596772093453531e-09
z =
Empty matrix: 0-by-1
p =
-1.398710433237015e-01 + 6.735753852075907e-01i
-1.398710433237015e-01 - 6.735753852075907e-01i
-2.311112505805679e-01 + 5.965841249062884e-01i
-2.311112505805679e-01 - 5.965841249062884e-01i
-2.970756579733747e-01 + 5.286265803850322e-01i
-2.970756579733747e-01 - 5.286265803850322e-01i
-5.055436388355750e-01 + 5.658992268447834e-02i
-5.055436388355750e-01 - 5.658992268447834e-02i
-4.989198616630246e-01 + 1.132893832531636e-01i
-4.989198616630246e-01 - 1.132893832531636e-01i
-4.877302860874728e-01 + 1.702165916211008e-01i
-4.877302860874728e-01 - 1.702165916211008e-01i
-4.717331331160929e-01 + 2.275088451191541e-01i
-4.717331331160929e-01 - 2.275088451191541e-01i
-4.505514519448975e-01 + 2.853371310741720e-01i
-4.505514519448975e-01 - 2.853371310741720e-01i
-4.236209682271996e-01 + 3.439293121342233e-01i
-4.236209682271996e-01 - 3.439293121342233e-01i
-3.900940978977325e-01 + 4.036111169859529e-01i
-3.900940978977325e-01 - 4.036111169859529e-01i
-3.486505808113257e-01 + 4.648872649574286e-01i
-3.486505808113257e-01 - 4.648872649574286e-01i
-5.077370186741994e-01
k =
1.615565688961578e-06
b =
Columns 1 through 3
1.000000000000000e+00 4.279007703722472e-13 7.139162763562802e-11
Columns 4 through 5
-9.521656951460538e-09 -3.326167160852700e-06
a =
Columns 1 through 3
1.000000000000000e+00 9.603257618831046e+02 3.952395295441562e+05
Columns 4 through 5
8.434637829684693e+07 8.099999999999998e+09
========================== BUTTER =========================
z =
-1
p =
3.249196962329063e-01
k =
3.375401518835469e-01
z =
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
p =
2.068257195514592e-01 + 9.238294351481734e-01i
2.068257195514592e-01 - 9.238294351481734e-01i
1.874933103892023e-01 + 8.269455076775277e-01i
1.874933103892023e-01 - 8.269455076775277e-01i
1.717435567330153e-01 + 7.383078571194629e-01i
1.717435567330153e-01 - 7.383078571194629e-01i
1.588266870755982e-01 + 6.564623730651094e-01i
1.588266870755982e-01 - 6.564623730651094e-01i
1.481881532502603e-01 + 5.802343458081779e-01i
1.481881532502603e-01 - 5.802343458081779e-01i
1.394122576319697e-01 + 5.086609000582009e-01i
1.394122576319697e-01 - 5.086609000582009e-01i
1.321840881809715e-01 + 4.409411734716436e-01i
1.321840881809715e-01 - 4.409411734716436e-01i
1.262633413354405e-01 + 3.763990035551881e-01i
1.262633413354405e-01 - 3.763990035551881e-01i
1.214660449478046e-01 + 3.144545234797277e-01i
1.214660449478046e-01 - 3.144545234797277e-01i
1.104868766650320e-01 + 2.771505404367791e-02i
1.104868766650320e-01 - 2.771505404367791e-02i
1.111768629525075e-01 + 8.331369153155753e-02i
1.111768629525075e-01 - 8.331369153155753e-02i
1.125740630842972e-01 + 1.394219509611784e-01i
1.125740630842972e-01 - 1.394219509611784e-01i
1.147138487992747e-01 + 1.963932363793666e-01i
1.147138487992747e-01 - 1.963932363793666e-01i
1.176516491045901e-01 + 2.546021573417188e-01i
1.176516491045901e-01 - 2.546021573417188e-01i
k =
1.446671081817286e-06
z =
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
p =
-1.772572785680147e-01 + 9.276431102995948e-01i
-1.772572785680147e-01 - 9.276431102995948e-01i
-1.600766565322114e-01 + 8.264026279893268e-01i
-1.600766565322114e-01 - 8.264026279893268e-01i
-1.461948419016121e-01 + 7.341841939120078e-01i
-1.461948419016121e-01 - 7.341841939120078e-01i
-1.348975284762046e-01 + 6.493235066053785e-01i
-1.348975284762046e-01 - 6.493235066053785e-01i
-1.256628210712206e-01 + 5.704921366889227e-01i
-1.256628210712206e-01 - 5.704921366889227e-01i
-1.181038235962314e-01 + 4.966120551231630e-01i
-1.181038235962314e-01 - 4.966120551231630e-01i
-1.119304913239356e-01 + 4.267938916403775e-01i
-1.119304913239356e-01 - 4.267938916403775e-01i
-1.069237739782691e-01 + 3.602914879527338e-01i
-1.069237739782691e-01 - 3.602914879527338e-01i
-1.029178030691416e-01 + 2.964677964142126e-01i
-1.029178030691416e-01 - 2.964677964142126e-01i
-9.978747500816100e-02 + 2.347687643085738e-01i
-9.978747500816100e-02 - 2.347687643085738e-01i
-9.743974496324025e-02 + 1.747028739092479e-01i
-9.743974496324025e-02 - 1.747028739092479e-01i
-9.580754551625957e-02 + 1.158246860771989e-01i
-9.580754551625957e-02 - 1.158246860771989e-01i
-9.484562207782568e-02 + 5.772118357151691e-02i
-9.484562207782568e-02 - 5.772118357151691e-02i
-9.452783117928215e-02
k =
9.585686688851069e-09
z =
1
1
1
1
1
1
1
1
-1
-1
-1
-1
-1
-1
-1
-1
p =
4.979909925436156e-01 + 8.367609424799387e-01i
4.979909925436156e-01 - 8.367609424799387e-01i
4.913338722555539e-01 + 7.866774509868817e-01i
4.913338722555539e-01 - 7.866774509868817e-01i
5.035229361778706e-01 + 7.401147376726750e-01i
5.035229361778706e-01 - 7.401147376726750e-01i
5.307617160406101e-01 + 7.029184459442954e-01i
5.307617160406101e-01 - 7.029184459442954e-01i
5.680556159453138e-01 + 6.788228792952775e-01i
5.680556159453138e-01 - 6.788228792952775e-01i
6.100962560818854e-01 + 6.693849403338664e-01i
6.100962560818854e-01 - 6.693849403338664e-01i
6.904694312740631e-01 + 6.930501690145245e-01i
6.904694312740631e-01 - 6.930501690145245e-01i
6.521767004237027e-01 + 6.744414640183752e-01i
6.521767004237027e-01 - 6.744414640183752e-01i
k =
3.398854055800844e-08
z =
-1.594474531383421e-02 + 9.998728744679880e-01i
-1.594474531383421e-02 - 9.998728744679880e-01i
-1.594474531383421e-02 + 9.998728744679880e-01i
-1.594474531383421e-02 - 9.998728744679880e-01i
-1.594474531383421e-02 + 9.998728744679880e-01i
-1.594474531383421e-02 - 9.998728744679880e-01i
-1.594474531383421e-02 + 9.998728744679880e-01i
-1.594474531383421e-02 - 9.998728744679880e-01i
-1.594474531383421e-02 + 9.998728744679880e-01i
-1.594474531383421e-02 - 9.998728744679880e-01i
-1.594474531383421e-02 + 9.998728744679880e-01i
-1.594474531383421e-02 - 9.998728744679880e-01i
-1.594474531383421e-02 + 9.998728744679880e-01i
-1.594474531383421e-02 - 9.998728744679880e-01i
p =
-1.766850742887729e-01 + 9.466951258673900e-01i
-1.766850742887729e-01 - 9.466951258673900e-01i
1.467897662432886e-01 + 9.515917126462422e-01i
1.467897662432886e-01 - 9.515917126462422e-01i
-1.370083529426906e-01 + 8.880376681273993e-01i
-1.370083529426906e-01 - 8.880376681273993e-01i
1.086774544701390e-01 + 8.915240810704319e-01i
1.086774544701390e-01 - 8.915240810704319e-01i
-7.982704457700891e-02 + 8.506056315273435e-01i
-7.982704457700891e-02 - 8.506056315273435e-01i
5.238812787110331e-02 + 8.524011102699969e-01i
5.238812787110331e-02 - 8.524011102699969e-01i
-1.357545000491310e-02 + 8.382287744986582e-01i
-1.357545000491310e-02 - 8.382287744986582e-01i
k =
4.577122512960063e-01
b =
Columns 1 through 3
0 0 0
Columns 4 through 6
0 1.599999999999999e+09 0
Columns 7 through 9
0 0 0
a =
Columns 1 through 3
1.000000000000000e+00 5.226251859505511e+02 2.565685424949238e+05
Columns 4 through 6
6.794127417357160e+07 1.519411254969542e+10 2.038238225207147e+12
Columns 7 through 9
2.309116882454312e+14 1.411088002066486e+16 8.099999999999991e+17
z =
0
0
0
0
p =
-4.179137760733086e+00 + 1.095935899082837e+02i
-4.179137760733086e+00 - 1.095935899082837e+02i
-9.593598668443835e+00 + 1.034745398029734e+02i
-9.593598668443835e+00 - 1.034745398029734e+02i
-8.883991981781929e+00 + 9.582087115567160e+01i
-8.883991981781929e+00 - 9.582087115567160e+01i
-3.474530886568715e+00 + 9.111599925805801e+01i
-3.474530886568715e+00 - 9.111599925805801e+01i
k =
1.600000000000001e+05
========================== CHEBY1 =========================
z =
-1
p =
-5.390126972799615e-01
k =
7.695063486399808e-01
z =
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
p =
-6.136558509657073e-01 + 2.700091504942893e-01i
-6.136558509657073e-01 - 2.700091504942893e-01i
-3.303348340927516e-01 + 6.659400861114254e-01i
-3.303348340927516e-01 - 6.659400861114254e-01i
8.779713780557169e-03 + 8.223108447483040e-01i
8.779713780557169e-03 - 8.223108447483040e-01i
2.742361123006911e-01 + 8.356666951611864e-01i
2.742361123006911e-01 - 8.356666951611864e-01i
4.562984557158206e-01 + 7.954276912303594e-01i
4.562984557158206e-01 - 7.954276912303594e-01i
5.777335494123628e-01 + 7.435821817961783e-01i
5.777335494123628e-01 - 7.435821817961783e-01i
6.593260977749194e-01 + 6.955390907990932e-01i
6.593260977749194e-01 - 6.955390907990932e-01i
7.149590948466562e-01 + 6.559437858502012e-01i
7.149590948466562e-01 - 6.559437858502012e-01i
7.532432388188739e-01 + 6.256158042292060e-01i
7.532432388188739e-01 - 6.256158042292060e-01i
7.794365244268271e-01 + 6.042099234813333e-01i
7.794365244268271e-01 - 6.042099234813333e-01i
7.967253874772997e-01 + 5.911966597313203e-01i
7.967253874772997e-01 - 5.911966597313203e-01i
8.069756417293870e-01 + 5.862214589217275e-01i
8.069756417293870e-01 - 5.862214589217275e-01i
k =
6.190427617192018e-04
z =
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
p =
-9.574019834435322e-01
-9.541085838336029e-01 + 8.482104798298945e-02i
-9.541085838336029e-01 - 8.482104798298945e-02i
-9.445488056314215e-01 + 1.673388481770735e-01i
-9.445488056314215e-01 - 1.673388481770735e-01i
-9.296386567980504e-01 + 2.454488552711194e-01i
-9.296386567980504e-01 - 2.454488552711194e-01i
-9.107657310168502e-01 + 3.174021886467137e-01i
-9.107657310168502e-01 - 3.174021886467137e-01i
-8.896168944779259e-01 + 3.818926727829247e-01i
-8.896168944779259e-01 - 3.818926727829247e-01i
-8.068787133559521e-01 + 5.861636377350449e-01i
-8.068787133559521e-01 - 5.861636377350449e-01i
-8.089131305095553e-01 + 5.742050583785998e-01i
-8.089131305095553e-01 - 5.742050583785998e-01i
-8.169256222996423e-01 + 5.534810333010162e-01i
-8.169256222996423e-01 - 5.534810333010162e-01i
-8.301876822600938e-01 + 5.239306947572915e-01i
-8.301876822600938e-01 - 5.239306947572915e-01i
-8.679949999243093e-01 + 4.380664730669522e-01i
-8.679949999243093e-01 - 4.380664730669522e-01i
-8.476563955221054e-01 + 4.854659489348472e-01i
-8.476563955221054e-01 - 4.854659489348472e-01i
k =
-5.544055466212737e-17
z =
1
1
1
1
1
1
1
1
-1
-1
-1
-1
-1
-1
-1
-1
p =
3.077784854851463e-01 + 9.453307017592942e-01i
3.077784854851463e-01 - 9.453307017592942e-01i
3.280567400654425e-01 + 9.272377218689016e-01i
3.280567400654425e-01 - 9.272377218689016e-01i
3.677912763284301e-01 + 9.038008865279966e-01i
3.677912763284301e-01 - 9.038008865279966e-01i
4.194425632520948e-01 + 8.769407159656157e-01i
4.194425632520948e-01 - 8.769407159656157e-01i
4.740921994669189e-01 + 8.496508528630974e-01i
4.740921994669189e-01 - 8.496508528630974e-01i
5.234866481897429e-01 + 8.259608422808477e-01i
5.234866481897429e-01 - 8.259608422808477e-01i
5.844717632289875e-01 + 8.052901363500210e-01i
5.844717632289875e-01 - 8.052901363500210e-01i
5.615189063336070e-01 + 8.100667803850766e-01i
5.615189063336070e-01 - 8.100667803850766e-01i
k =
5.007028718074307e-09
z =
-1.583844403245361e-01 + 9.873775210440450e-01i
-1.583844403245361e-01 - 9.873775210440450e-01i
-1.583844403245361e-01 + 9.873775210440450e-01i
-1.583844403245361e-01 - 9.873775210440450e-01i
-1.583844403245361e-01 + 9.873775210440450e-01i
-1.583844403245361e-01 - 9.873775210440450e-01i
-1.583844403245361e-01 + 9.873775210440450e-01i
-1.583844403245361e-01 - 9.873775210440450e-01i
-1.583844403245361e-01 + 9.873775210440450e-01i
-1.583844403245361e-01 - 9.873775210440450e-01i
-1.583844403245361e-01 + 9.873775210440450e-01i
-1.583844403245361e-01 - 9.873775210440450e-01i
-1.583844403245361e-01 + 9.873775210440450e-01i
-1.583844403245361e-01 - 9.873775210440450e-01i
p =
-8.942974551472813e-02 + 3.482480481185926e-01i
-8.942974551472813e-02 - 3.482480481185926e-01i
1.293775154041798e-01 + 8.753499858081858e-01i
1.293775154041798e-01 - 8.753499858081858e-01i
3.399741945062013e-02 + 9.690316022705607e-01i
3.399741945062013e-02 - 9.690316022705607e-01i
4.167225522796539e-04 + 9.927338161087488e-01i
4.167225522796539e-04 - 9.927338161087488e-01i
-3.912966549550960e-01 + 8.046122859255742e-01i
-3.912966549550960e-01 - 8.046122859255742e-01i
-3.307805547127368e-01 + 9.133455018206508e-01i
-3.307805547127368e-01 - 9.133455018206508e-01i
-3.072658345097743e-01 + 9.443589759799366e-01i
-3.072658345097743e-01 - 9.443589759799366e-01i
k =
3.619438310405028e-01
b =
Columns 1 through 3
1.000000000000006e+00 0 3.255000000000020e+05
Columns 4 through 6
0 4.238010000000026e+10 0
Columns 7 through 9
2.758944510000017e+15 0 8.980364380050052e+19
Columns 10 through 11
0 1.169243442282517e+24
a =
Columns 1 through 3
1.000000000000000e+00 4.630555945694342e+02 4.039266454794788e+05
Columns 4 through 6
1.338060988610237e+08 5.844333551294591e+10 1.357346371637638e+13
Columns 7 through 9
3.804661141892782e+15 5.670715850340080e+17 1.114411200988328e+20
Columns 10 through 11
8.316815934908471e+21 1.169243442282517e+24
z =
0
0
0
0
0
0
0
0
0
0
p =
-3.144743169501551e+03 + 3.511680029092744e+03i
-3.144743169501551e+03 - 3.511680029092744e+03i
-5.633065604514602e+02 + 2.023615191183945e+03i
-5.633065604514602e+02 - 2.023615191183945e+03i
-1.946412183352025e+02 + 1.372309454274755e+03i
-1.946412183352025e+02 - 1.372309454274755e+03i
-7.987162953085479e+01 + 1.105207708045358e+03i
-7.987162953085479e+01 - 1.105207708045358e+03i
-2.250315039031946e+01 + 1.001723931471477e+03i
-2.250315039031946e+01 - 1.001723931471477e+03i
k =
8.912509381337453e-01
========================== CHEBY2 =========================
z =
-1
p =
9.967826460175649e-01
k =
1.608676991217512e-03
z =
9.981088955489852e-01 + 6.147058341984388e-02i
9.981088955489852e-01 - 6.147058341984388e-02i
9.832702870387426e-01 + 1.821525257215483e-01i
9.832702870387426e-01 - 1.821525257215483e-01i
9.550760158089112e-01 + 2.963609353922882e-01i
9.550760158089112e-01 - 2.963609353922882e-01i
9.162054748821922e-01 + 4.007087817803773e-01i
9.162054748821922e-01 - 4.007087817803773e-01i
8.700619897368064e-01 + 4.929423232136168e-01i
8.700619897368064e-01 - 4.929423232136168e-01i
5.889791753434985e-01 + 8.081482110427953e-01i
5.889791753434985e-01 - 8.081482110427953e-01i
5.984900456570295e-01 + 8.011302423760501e-01i
5.984900456570295e-01 - 8.011302423760501e-01i
6.172880888914629e-01 + 7.867371958365343e-01i
6.172880888914629e-01 - 7.867371958365343e-01i
6.448899971038180e-01 + 7.642754030030161e-01i
6.448899971038180e-01 - 7.642754030030161e-01i
6.804845629637927e-01 + 7.327624168637228e-01i
6.804845629637927e-01 - 7.327624168637228e-01i
8.202619107108660e-01 + 5.719881098737678e-01i
8.202619107108660e-01 - 5.719881098737678e-01i
7.228410452536148e-01 + 6.910143437705678e-01i
7.228410452536148e-01 - 6.910143437705678e-01i
7.702121399578629e-01 + 6.377877856007792e-01i
7.702121399578629e-01 - 6.377877856007792e-01i
p =
7.365546198286450e-01 + 4.842085129329526e-02i
7.365546198286450e-01 - 4.842085129329526e-02i
7.292038510962885e-01 + 1.442201672097581e-01i
7.292038510962885e-01 - 1.442201672097581e-01i
7.151293788040354e-01 + 2.369925800458584e-01i
7.151293788040354e-01 - 2.369925800458584e-01i
6.955051820787286e-01 + 3.250341363856910e-01i
6.955051820787286e-01 - 3.250341363856910e-01i
6.719122956045220e-01 + 4.070475750638047e-01i
6.719122956045220e-01 - 4.070475750638047e-01i
6.461722130611300e-01 + 4.821965916689270e-01i
6.461722130611300e-01 - 4.821965916689270e-01i
5.528045062872224e-01 + 8.162920513838372e-01i
5.528045062872224e-01 - 8.162920513838372e-01i
5.464847782492791e-01 + 7.869899955967304e-01i
5.464847782492791e-01 - 7.869899955967304e-01i
5.488033111260949e-01 + 7.520442354055579e-01i
5.488033111260949e-01 - 7.520442354055579e-01i
6.201874719022955e-01 + 5.500894392527353e-01i
6.201874719022955e-01 - 5.500894392527353e-01i
5.586478152536709e-01 + 7.112676877332921e-01i
5.586478152536709e-01 - 7.112676877332921e-01i
5.958145844148228e-01 + 6.107074340842115e-01i
5.958145844148228e-01 - 6.107074340842115e-01i
5.747812938519067e-01 + 6.643001536914696e-01i
5.747812938519067e-01 - 6.643001536914696e-01i
k =
9.932997786497189e-02
z =
9.690690376586687e-01 + 2.467897896011971e-01i
9.690690376586687e-01 - 2.467897896011971e-01i
9.999999999999492e-01
8.835111277191199e-01 + 4.684101698261429e-01i
8.835111277191199e-01 - 4.684101698261429e-01i
7.613142857900539e-01 + 6.483830335935022e-01i
7.613142857900539e-01 - 6.483830335935022e-01i
6.232625173626231e-01 + 7.820126817709752e-01i
6.232625173626231e-01 - 7.820126817709752e-01i
4.864456563413621e-01 + 8.737108351316745e-01i
4.864456563413621e-01 - 8.737108351316745e-01i
3.618368136816749e-01 + 9.322414495530347e-01i
3.618368136816749e-01 - 9.322414495530347e-01i
2.549486883466794e-01 + 9.669545833752675e-01i
2.549486883466794e-01 - 9.669545833752675e-01i
1.676175432109457e-01 + 9.858520980390212e-01i
1.676175432109457e-01 - 9.858520980390212e-01i
1.975218468277521e-03 + 9.999980492540941e-01i
1.975218468277521e-03 - 9.999980492540941e-01i
1.786959496651858e-02 + 9.998403260399917e-01i
1.786959496651858e-02 - 9.998403260399917e-01i
9.967933660557139e-02 + 9.950196127985684e-01i
9.967933660557139e-02 - 9.950196127985684e-01i
5.013970951219547e-02 + 9.987422137518890e-01i
5.013970951219547e-02 - 9.987422137518890e-01i
p =
4.218866331906864e-01
4.120110200127552e-01 + 1.361290593621978e-01i
4.120110200127552e-01 - 1.361290593621978e-01i
3.835890113632530e-01 + 2.664910809911026e-01i
3.835890113632530e-01 - 2.664910809911026e-01i
3.399195570456499e-01 + 3.863983538639875e-01i
3.399195570456499e-01 - 3.863983538639875e-01i
2.855977834508353e-01 + 4.929444399540688e-01i
2.855977834508353e-01 - 4.929444399540688e-01i
2.255765441339322e-01 + 5.851631870205766e-01i
2.255765441339322e-01 - 5.851631870205766e-01i
1.644087535815792e-01 + 6.637356937277153e-01i
1.644087535815792e-01 - 6.637356937277153e-01i
-7.293633845273095e-02 + 9.739218252516307e-01i
-7.293633845273095e-02 - 9.739218252516307e-01i
1.058259206358626e-01 + 7.304739464862978e-01i
1.058259206358626e-01 - 7.304739464862978e-01i
-5.703971947785402e-02 + 9.291057542169088e-01i
-5.703971947785402e-02 - 9.291057542169088e-01i
5.263875132656864e-02 + 7.877974334424453e-01i
5.263875132656864e-02 - 7.877974334424453e-01i
-3.007943405982616e-02 + 8.846331716180016e-01i
-3.007943405982616e-02 - 8.846331716180016e-01i
6.857277464483946e-03 + 8.383275456264492e-01i
6.857277464483946e-03 - 8.383275456264492e-01i
k =
6.507068761705037e-03
z =
-9.999999999999999e-01
3.676588029658514e-01 + 9.299607543341383e-01i
3.676588029658514e-01 - 9.299607543341383e-01i
7.009689684982283e-01 + 7.131917730894889e-01i
7.009689684982283e-01 - 7.131917730894889e-01i
7.815697973765858e-01 + 6.238178033919218e-01i
7.815697973765858e-01 - 6.238178033919218e-01i
8.063793628819866e-01 + 5.913986160941200e-01i
8.063793628819866e-01 - 5.913986160941200e-01i
1.000000000000001e+00
9.944493019920448e-01 + 1.052168511576739e-01i
9.944493019920448e-01 - 1.052168511576739e-01i
9.854674703367308e-01 + 1.698642543566085e-01i
9.854674703367308e-01 - 1.698642543566085e-01i
9.762751735919308e-01 + 2.165335665157851e-01i
9.762751735919308e-01 - 2.165335665157851e-01i
9.792277171575134e-01 + 2.027636011479496e-01i
9.792277171575134e-01 - 2.027636011479496e-01i
p =
8.143803410489621e-01 + 5.411056063397541e-01i
8.143803410489621e-01 - 5.411056063397541e-01i
7.650769827887418e-01 + 5.195412242095543e-01i
7.650769827887418e-01 - 5.195412242095543e-01i
6.096241204063443e-01 + 3.568440484659796e-01i
6.096241204063443e-01 - 3.568440484659796e-01i
6.918192770246239e-01 + 4.770463577106911e-01i
6.918192770246239e-01 - 4.770463577106911e-01i
6.986241085779207e-01 + 1.146512226180060e-01i
6.986241085779207e-01 - 1.146512226180060e-01i
8.654645923909734e-01 + 1.604208797063147e-01i
8.654645923909734e-01 - 1.604208797063147e-01i
9.164831670444591e-01 + 1.969181049384918e-01i
9.164831670444591e-01 - 1.969181049384918e-01i
9.630425777594550e-01 + 2.317513360702271e-01i
9.630425777594550e-01 - 2.317513360702271e-01i
9.438104703725529e-01 + 2.193509900269860e-01i
9.438104703725529e-01 - 2.193509900269860e-01i
k =
9.345352824659604e-03
z =
6.230544895101009e-01 + 7.821784343111114e-01i
6.230544895101009e-01 - 7.821784343111114e-01i
9.086608545660115e-01 + 4.175349702471991e-01i
9.086608545660115e-01 - 4.175349702471991e-01i
9.478129721465802e-01 + 3.188268649763867e-01i
9.478129721465802e-01 - 3.188268649763867e-01i
-6.230544895100982e-01 + 7.821784343111109e-01i
-6.230544895100982e-01 - 7.821784343111109e-01i
-9.086608545660116e-01 + 4.175349702472088e-01i
-9.086608545660116e-01 - 4.175349702472088e-01i
-9.478129721465784e-01 + 3.188268649763897e-01i
-9.478129721465784e-01 - 3.188268649763897e-01i
p =
-9.464094036167638e-01 + 1.720048695084344e-01i
-9.464094036167638e-01 - 1.720048695084344e-01i
-8.715844103386737e-01 + 1.370665039509297e-01i
-8.715844103386737e-01 - 1.370665039509297e-01i
-8.078751204586425e-01 + 5.729329866682983e-02i
-8.078751204586425e-01 - 5.729329866682983e-02i
9.464094036167665e-01 + 1.720048695084332e-01i
9.464094036167665e-01 - 1.720048695084332e-01i
8.078751204586447e-01 + 5.729329866683007e-02i
8.078751204586447e-01 - 5.729329866683007e-02i
8.715844103386721e-01 + 1.370665039509331e-01i
8.715844103386721e-01 - 1.370665039509331e-01i
k =
2.917823332763358e-03
b =
Columns 1 through 3
1.000000000000000e+00 6.683253076978249e-12 2.111512500000000e+07
Columns 4 through 6
1.134325604589552e-04 1.782966433781250e+14 7.216787944356781e+02
Columns 7 through 9
7.525901316990656e+20 2.039829265789886e+09 1.587960565565748e+27
Columns 10 through 11
2.161236218626134e+15 1.339913493808585e+33
a =
Columns 1 through 3
1.000000000000000e+00 1.849550755473371e+02 2.113222918998538e+07
Columns 4 through 6
3.125114149732283e+09 1.785133457155609e+14 1.979158697776348e+16
Columns 7 through 9
7.535048322653831e+20 5.567966191263037e+22 1.589246884221346e+27
Columns 10 through 11
5.871210648525566e+28 1.339913493808590e+33
========================== ELLIP =========================
z =
-9.999999999999998e-01
p =
-6.660721153525525e-04
k =
5.003330360576763e-01
z =
9.761875332501075e-01 + 2.169283290099910e-01i
9.761875332501075e-01 - 2.169283290099910e-01i
8.413503353963494e-01 + 5.404901600661900e-01i
8.413503353963494e-01 - 5.404901600661900e-01i
7.160082576305009e-01 + 6.980918098681732e-01i
7.160082576305009e-01 - 6.980918098681732e-01i
6.456533638965329e-01 + 7.636306264739803e-01i
6.456533638965329e-01 - 7.636306264739803e-01i
6.127321820971366e-01 + 7.902906256703928e-01i
6.127321820971366e-01 - 7.902906256703928e-01i
5.983607817490196e-01 + 8.012267936512676e-01i
5.983607817490196e-01 - 8.012267936512676e-01i
5.922577552594799e-01 + 8.057485658286990e-01i
5.922577552594799e-01 - 8.057485658286990e-01i
5.896952092563588e-01 + 8.076258788449631e-01i
5.896952092563588e-01 - 8.076258788449631e-01i
5.886248765538837e-01 + 8.084063054565607e-01i
5.886248765538837e-01 - 8.084063054565607e-01i
5.881802711123132e-01 + 8.087298490066037e-01i
5.881802711123132e-01 - 8.087298490066037e-01i
5.879995719101164e-01 + 8.088612386766461e-01i
5.879995719101164e-01 - 8.088612386766461e-01i
5.879354086709576e-01 + 8.089078780868164e-01i
5.879354086709576e-01 - 8.089078780868164e-01i
p =
-3.184805259081650e-01 + 4.206951906775851e-01i
-3.184805259081650e-01 - 4.206951906775851e-01i
1.417279173459985e-01 + 7.903955262836452e-01i
1.417279173459985e-01 - 7.903955262836452e-01i
4.042881216964651e-01 + 8.309042239116594e-01i
4.042881216964651e-01 - 8.309042239116594e-01i
5.128964442789670e-01 + 8.229563236799665e-01i
5.128964442789670e-01 - 8.229563236799665e-01i
5.569614712822724e-01 + 8.155957702908510e-01i
5.569614712822724e-01 - 8.155957702908510e-01i
5.750478870161392e-01 + 8.118633973883931e-01i
5.750478870161392e-01 - 8.118633973883931e-01i
5.825314018170804e-01 + 8.101960910679270e-01i
5.825314018170804e-01 - 8.101960910679270e-01i
5.856397379751872e-01 + 8.094825218722543e-01i
5.856397379751872e-01 - 8.094825218722543e-01i
5.869326035251949e-01 + 8.091827531557583e-01i
5.869326035251949e-01 - 8.091827531557583e-01i
5.874697218855733e-01 + 8.090593298213502e-01i
5.874697218855733e-01 - 8.090593298213502e-01i
5.876904783532237e-01 + 8.090127161018823e-01i
5.876904783532237e-01 - 8.090127161018823e-01i
5.877753105317594e-01 + 8.090050577978136e-01i
5.877753105317594e-01 - 8.090050577978136e-01i
k =
4.918081266957108e-02
z =
9.999999999998661e-01
6.603717261750994e-01 + 7.509388678638675e-01i
6.603717261750994e-01 - 7.509388678638675e-01i
2.788635267510325e-01 + 9.603307416968041e-01i
2.788635267510325e-01 - 9.603307416968041e-01i
1.070215532544218e-01 + 9.942567008268131e-01i
1.070215532544218e-01 - 9.942567008268131e-01i
4.049427369978163e-02 + 9.991797705105507e-01i
4.049427369978163e-02 - 9.991797705105507e-01i
1.531059368627931e-02 + 9.998827859909265e-01i
1.531059368627931e-02 - 9.998827859909265e-01i
5.808061438534933e-03 + 9.999831330689181e-01i
5.808061438534933e-03 - 9.999831330689181e-01i
2.224277847754599e-03 + 9.999975262909676e-01i
2.224277847754599e-03 - 9.999975262909676e-01i
8.731857107534554e-04 + 9.999996187732845e-01i
8.731857107534554e-04 - 9.999996187732845e-01i
3.649057346914968e-04 + 9.999999334218996e-01i
3.649057346914968e-04 - 9.999999334218996e-01i
1.765538109802615e-04 + 9.999999844143768e-01i
1.765538109802615e-04 - 9.999999844143768e-01i
1.143655290967426e-04 + 9.999999934602630e-01i
1.143655290967426e-04 - 9.999999934602630e-01i
p =
-6.322017026545028e-01
-4.648423756662754e-01 + 5.852407464440732e-01i
-4.648423756662754e-01 - 5.852407464440732e-01i
-2.249233374627773e-01 + 8.577853017985717e-01i
-2.249233374627773e-01 - 8.577853017985717e-01i
-9.234137570557621e-02 + 9.506548198678851e-01i
-9.234137570557621e-02 - 9.506548198678851e-01i
-3.585663561241373e-02 + 9.821494736043981e-01i
-3.585663561241373e-02 - 9.821494736043981e-01i
-1.363917242312723e-02 + 9.933844128330656e-01i
-1.363917242312723e-02 - 9.933844128330656e-01i
-5.131505238923029e-03 + 9.975221173308673e-01i
-5.131505238923029e-03 - 9.975221173308673e-01i
-1.904937999259502e-03 + 9.990680819857982e-01i
-1.904937999259502e-03 - 9.990680819857982e-01i
-6.859439885466834e-04 + 9.996492201426826e-01i
-6.859439885466834e-04 - 9.996492201426826e-01i
-2.269936267937089e-04 + 9.998686250679161e-01i
-2.269936267937089e-04 - 9.998686250679161e-01i
-5.687071588789117e-05 + 9.999527573294513e-01i
-5.687071588789117e-05 - 9.999527573294513e-01i
-6.948417068525226e-07 + 9.999882737700173e-01i
-6.948417068525226e-07 - 9.999882737700173e-01i
k =
1.220910020289434e-02
z =
-9.999999999999991e-01
6.856610961780020e-01 + 7.279209168501619e-01i
6.856610961780020e-01 - 7.279209168501619e-01i
7.850346167691289e-01 + 6.194518952058737e-01i
7.850346167691289e-01 - 6.194518952058737e-01i
7.999038743173071e-01 + 6.001281461922627e-01i
7.999038743173071e-01 - 6.001281461922627e-01i
9.999999999999999e-01
9.862938983554124e-01 + 1.649980183725925e-01i
9.862938983554124e-01 - 1.649980183725925e-01i
9.788558330548762e-01 + 2.045513580850601e-01i
9.788558330548762e-01 - 2.045513580850601e-01i
9.771155231720003e-01 + 2.127093189691258e-01i
9.771155231720003e-01 - 2.127093189691258e-01i
p =
8.063992755498643e-01 + 5.858071374778874e-01i
8.063992755498643e-01 - 5.858071374778874e-01i
8.050395347071724e-01 + 5.639097428109795e-01i
8.050395347071724e-01 - 5.639097428109795e-01i
8.113124936559144e-01 + 4.855241143973142e-01i
8.113124936559144e-01 - 4.855241143973142e-01i
8.665595314082394e-01 + 3.334049560919331e-01i
8.665595314082394e-01 - 3.334049560919331e-01i
9.412369011968871e-01 + 2.457616651325908e-01i
9.412369011968871e-01 - 2.457616651325908e-01i
9.679465190411238e-01 + 2.228772501848216e-01i
9.679465190411238e-01 - 2.228772501848216e-01i
9.747235066273385e-01 + 2.178937926146544e-01i
9.747235066273385e-01 - 2.178937926146544e-01i
k =
8.354782670263239e-03
z =
3.528578094286510e-01 + 9.356769561794296e-01i
3.528578094286510e-01 - 9.356769561794296e-01i
3.769716042264783e-01 + 9.262248159096587e-01i
3.769716042264783e-01 - 9.262248159096587e-01i
4.406101783111199e-01 + 8.976985411420985e-01i
4.406101783111199e-01 - 8.976985411420985e-01i
5.539386470258847e-01 + 8.325574907062760e-01i
5.539386470258847e-01 - 8.325574907062760e-01i
6.748464963023645e-01 + 7.379581332490555e-01i
6.748464963023645e-01 - 7.379581332490555e-01i
7.489887970285254e-01 + 6.625826604475596e-01i
7.489887970285254e-01 - 6.625826604475596e-01i
7.913118471618432e-01 + 6.114127579150699e-01i
7.913118471618432e-01 - 6.114127579150699e-01i
7.806804740916381e-01 + 6.249303940216475e-01i
7.806804740916381e-01 - 6.249303940216475e-01i
p =
-1.025299146693730e-01 + 5.662682444754943e-01i
-1.025299146693730e-01 - 5.662682444754943e-01i
1.698463595163031e-01 + 8.926678667070186e-01i
1.698463595163031e-01 - 8.926678667070186e-01i
2.750532687820631e-01 + 9.351020170094005e-01i
2.750532687820631e-01 - 9.351020170094005e-01i
3.070095178909486e-01 + 9.457373499553291e-01i
3.070095178909486e-01 - 9.457373499553291e-01i
7.695332312152288e-01 + 2.792567212705257e-01i
7.695332312152288e-01 - 2.792567212705257e-01i
8.083818999225620e-01 + 4.990723496863960e-01i
8.083818999225620e-01 - 4.990723496863960e-01i
8.066158014414928e-01 + 5.649811440393374e-01i
8.066158014414928e-01 - 5.649811440393374e-01i
8.062787978834571e-01 + 5.855780880424964e-01i
8.062787978834571e-01 - 5.855780880424964e-01i
k =
2.068622545291259e-01
b =
Columns 1 through 3
1.000000000000000e+00 1.743506051190569e-13 2.426561778314366e+05
Columns 4 through 6
3.459426536825722e-08 2.348218683400168e+10 2.559179747299313e-03
Columns 7 through 9
1.132780692872241e+15 8.363229375535731e+01 2.724038554089566e+19
Columns 10 through 11
1.018700994113120e+06 2.612380874940186e+23
a =
Columns 1 through 3
1.000000000000000e+00 1.337266601804649e+02 2.486725353510667e+05
Columns 4 through 6
2.628059713728125e+07 2.436169536928770e+10 1.913554568577315e+12
Columns 7 through 9
1.175208184614438e+15 6.115751452473410e+16 2.791577695211466e+19
Columns 10 through 11
7.241811142725384e+20 2.612380874940182e+23
z =
-5.583607317695175e-14 + 1.433755965989225e+02i
-5.583607317695175e-14 - 1.433755965989225e+02i
5.740106416459296e-14 + 1.261678754570291e+02i
5.740106416459296e-14 - 1.261678754570291e+02i
-2.199676239638652e-14 + 6.974861996895196e+01i
-2.199676239638652e-14 - 6.974861996895196e+01i
-3.372595657044283e-14 + 7.926145989044531e+01i
-3.372595657044283e-14 - 7.926145989044531e+01i
0
p =
-8.814960004852743e-01 + 1.104124501436066e+02i
-8.814960004852743e-01 - 1.104124501436066e+02i
-2.477372459140184e+00 + 1.065638954516534e+02i
-2.477372459140184e+00 - 1.065638954516534e+02i
-3.072156842945799e+00 + 9.995404870405324e+01i
-3.072156842945799e+00 - 9.995404870405324e+01i
-2.180456023925693e+00 + 9.379206865455268e+01i
-2.180456023925693e+00 - 9.379206865455268e+01i
-7.230484977485752e-01 + 9.056598800801140e+01i
-7.230484977485752e-01 - 9.056598800801140e+01i
k =
3.774571622827070e-02
========================== BUTTORD =========================
n =
16
Wn =
2.034240543183123e-01
n =
18
Wn =
2.994397774822765e-01
n =
18
Wn =
1.988019498910527e-01 5.020435331649425e-01
n =
20
Wn =
1.481914223819211e-01 5.986892392286634e-01
========================== CHEB1ORD =========================
n =
8
Wn =
2.000000000000000e-01
n =
9
Wn =
3.000000000000000e-01
n =
9
Wn =
2.000000000000000e-01 5.000000000000000e-01
n =
10
Wn =
1.000000000000000e-01 6.000000000000000e-01
========================== CHEB2ORD =========================
n =
8
Wn =
3.000000000000000e-01
n =
9
Wn =
2.000000000000000e-01
n =
9
Wn =
1.000000000000000e-01 6.000000000000000e-01
n =
10
Wn =
2.000000000000000e-01 5.000000000000000e-01
========================== ELLIPORD =========================
n =
5
Wn =
2.000000000000000e-01
n =
6
Wn =
3.000000000000000e-01
n =
6
Wn =
2.000000000000000e-01 5.000000000000000e-01
n =
7
Wn =
1.000000000000000e-01 6.000000000000000e-01
========================== ANALOG =========================
n =
7
Wn =
2.236556391943067e+02
n =
4
Wn =
700
n =
11
Wn =
1.673740595370124e+01 5.974641487254268e+01
n =
9
Wn =
1000 6000
========================== CHEB2ORD =========================
n =
4
Ws =
3.000000000000000e-01
n =
7
Ws =
1.000000000000000e-01 5.000000000000000e-01
>>
clear all; close all; clc
format long e
format compact
disp('========================== BESSELF =========================')
%[z, p, k] = besself(0, 0.3) % degenerate (octave doesn't allow this)
[z, p, k] = besself(1, 0.3) % first order is same for all types
[z, p, k] = besself(24, 0.43) % high even order
[z, p, k] = besself(23, 0.56) % high odd order
% with transfer function conversion, without digital conversion
[b, a] = besself(4, 300, 's')
disp('========================== BUTTER =========================')
%[z, p, k] = butter(0, 0.3) % degenerate (octave doesn't allow this)
[z, p, k] = butter(1, 0.3) % first order is same for all types
[z, p, k] = butter(28, 0.43, 'high') % high even order
[z, p, k] = butter(27, 0.56, 'high') % high odd order
[z, p, k] = butter(8, [0.25 0.33]) % bandpass
[z, p, k] = butter(7, [0.45 0.56], 'stop') % bandstop
% with transfer function conversion, without digital conversion
[b, a] = butter(4, [100 300], 'bandpass', 's')
[z, p, k] = butter(4, [90.5 110.5], 's') % Paarmann Example 9.11
disp('========================== CHEBY1 =========================')
%[z, p, k] = cheby1(0, 1, 0.6) % degenerate (octave doesn't allow this)
[z, p, k] = cheby1(1, 0.1, 0.3) % first order is same for all types
[z, p, k] = cheby1(24, 0.7, 0.2, 'high') % high even order
[z, p, k] = cheby1(23, 0.8, 0.8, 'high') % high odd order
[z, p, k] = cheby1(8, 1, [0.3 0.4]) % bandpass
[z, p, k] = cheby1(7, 1, [0.5 0.6], 'stop') % bandstop
% with transfer function conversion, without digital conversion
[b, a] = cheby1(5, 0.9, [210 310], 'stop', 's')
[z, p, k] = cheby1(10, 1, 1000, 'high', 's') % Paarmann Example 9.6
disp('========================== CHEBY2 =========================')
%[z, p, k] = cheby2(0, 40, 0.3) % degenerate (octave doesn't allow this)
[z, p, k] = cheby2(1, 50, 0.3) % first order is same for all types
[z, p, k] = cheby2(26, 60, 0.3, 'high') % high even order
[z, p, k] = cheby2(25, 80, 0.5, 'high') % high odd order
[z, p, k] = cheby2(9, 40, [0.07 0.2]) % bandpass
[z, p, k] = cheby2(6, 55, [0.1 0.9], 'stop') % bandstop
% with transfer function conversion, without digital conversion
[b, a] = cheby2(5, 20, [2010 2100], 'stop', 's')
disp('========================== ELLIP =========================')
%[z, p, k] = ellip(0, 1, 45, 0.1) % degenerate (octave doesn't allow this)
[z, p, k] = ellip(1, 1, 55, 0.3) % first order is same for all types
[z, p, k] = ellip(24, 1, 80, 0.3, 'high') % high even order
[z, p, k] = ellip(23, 1, 70, 0.5, 'high') % high odd order
[z, p, k] = ellip(7, 1, 40, [0.07 0.2]) % bandpass
[z, p, k] = ellip(8, 1, 65, [0.2 0.4], 'stop') % bandstop
% with transfer function conversion, without digital conversion
[b, a] = ellip(5, 1, 40, [201 240], 'stop', 's')
[z, p, k] = ellip(5, 1, 75, [90.5 110.5], 's') % Paarmann Example 9.12
disp('========================== BUTTORD =========================')
% minimum order calculations (digital):
[n, Wn] = buttord(0.2, 0.3, 3, 60) % lowpass
[n, Wn] = buttord(0.3, 0.2, 3, 70) % highpass
[n, Wn] = buttord([0.2, 0.5], [0.1, 0.6], 3, 80) % bandpass
[n, Wn] = buttord([0.1, 0.6], [0.2, 0.5], 3, 90) % bandstop
disp('========================== CHEB1ORD =========================')
[n, Wn] = cheb1ord(0.2, 0.3, 3, 60) % lowpass
[n, Wn] = cheb1ord(0.3, 0.2, 3, 70) % highpass
[n, Wn] = cheb1ord([0.2, 0.5], [0.1, 0.6], 3, 80) % bandpass
[n, Wn] = cheb1ord([0.1, 0.6], [0.2, 0.5], 3, 90) % bandstop
disp('========================== CHEB2ORD =========================')
[n, Wn] = cheb2ord(0.2, 0.3, 3, 60) % lowpass
[n, Wn] = cheb2ord(0.3, 0.2, 3, 70) % highpass
[n, Wn] = cheb2ord([0.2, 0.5], [0.1, 0.6], 3, 80) % bandpass
[n, Wn] = cheb2ord([0.1, 0.6], [0.2, 0.5], 3, 90) % bandstop
disp('========================== ELLIPORD =========================')
[n, Wn] = ellipord(0.2, 0.3, 3, 60) % lowpass
[n, Wn] = ellipord(0.3, 0.2, 3, 70) % highpass
[n, Wn] = ellipord([0.2, 0.5], [0.1, 0.6], 3, 80) % bandpass
[n, Wn] = ellipord([0.1, 0.6], [0.2, 0.5], 3, 90) % bandstop
disp('========================== ANALOG =========================')
% analog
[n, Wn] = buttord(200, 600, 3, 60, 's') % lowpass
[n, Wn] = cheb1ord(700, 100, 3, 70, 's') % highpass
[n, Wn] = cheb2ord([20, 50], [10, 60], 3, 80, 's') % bandpass
[n, Wn] = ellipord([1000, 6000], [2000, 5000], 3, 90, 's') % bandstop
disp('========================== CHEB2ORD =========================')
% octave does not match matlab's documentation:
[n, Ws] = cheb2ord(40/500, 150/500, 3, 60)
[n, Ws] = cheb2ord([60/500, 200/500], [50/500, 250/500], 3, 40)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment