Changes

Page history
notebooks authored by Lee McCuller's avatar Lee McCuller
Hide whitespace changes
Inline Side-by-side
notebooks/cplx/digest.md
View page @ d8d102ff
# fit_sequence_v1
# fit_sequence version 1
> `fit_sequence_v1`
> `v1.fit_sequence`
Version 1 smart fitter in IIRrational library. Uses SVD method with high order over fitting, then switches to nonlinear fits with heuristics to remove poles and zeros down to a reasonable system order.
......@@ -34,18 +34,18 @@ Performs the SVD for a rough initial guess
![Fitttt](plot-1-3.png)
### `1.4` choose poles
### `1.4` choose direct
> `if`
Chose the Poles SVD fitter as it had the smaller residual of 1.84e+02 vs. 2.06e+04 for the zeros
Chose the direct (non-SVD) fitter as it had the smallest residual of 4.22e+02
### `1.5` seq_iter_1
> `IIRrational.fit_poles, IIRrational.fit_zeros`
First iterations, enforcing mindelay and stabilized poles residual of 8.41e+03
First iterations, enforcing mindelay and stabilized poles residual of 4.05e+02
![Fitttt](plot-1-5.png)
......@@ -54,7 +54,7 @@ First iterations, enforcing mindelay and stabilized poles residual of 8.41e+03
> `IIRrational.fit_poles, IIRrational.fit_zeros`
First iterations, enforcing stabilized poles residual of 9.17e+02
First iterations, enforcing stabilized poles residual of 3.93e+02
![Fitttt](plot-1-6.png)
......@@ -63,7 +63,7 @@ First iterations, enforcing stabilized poles residual of 9.17e+02
> `IIRrational.fit_poles, IIRrational.fit_zeros`
First iterations, enforcing stabilized poles residual of 2.46e+02
First iterations, enforcing stabilized poles residual of 3.87e+02
![Fitttt](plot-1-7.png)
......@@ -72,7 +72,7 @@ First iterations, enforcing stabilized poles residual of 2.46e+02
> `IIRrational.fit_poles, IIRrational.fit_zeros`
First iterations, enforcing stabilized poles residual of 1.68e+02
First iterations, enforcing stabilized poles residual of 3.82e+02
![Fitttt](plot-1-8.png)
......@@ -112,16 +112,16 @@ ID Pairs for order reduction
![Fitttt](plot-3-1.png)
### `3.2` seq_order_reduce_opt_1
### `3.2` duals_ID_2
> `MultiReprFilter.optimize`
Optimize again after order reduction
ID Pairs for order reduction
![Fitttt](plot-3-2.png)
### `3.3` duals_ID_2
### `3.3` duals_ID_3
> `MultiReprFilter.optimize`
......@@ -130,6 +130,24 @@ ID Pairs for order reduction
![Fitttt](plot-3-3.png)
### `3.4` duals_ID_4
> `MultiReprFilter.optimize`
ID Pairs for order reduction
![Fitttt](plot-3-4.png)
### `3.5` duals_ID_5
> `MultiReprFilter.optimize`
ID Pairs for order reduction
![Fitttt](plot-3-5.png)
## `4` flip mindelay
## `5` order_reduce 2
......@@ -143,3 +161,21 @@ ID Pairs for order reduction
![Fitttt](plot-5-1.png)
### `5.2` duals_ID_2
> `MultiReprFilter.optimize`
ID Pairs for order reduction
![Fitttt](plot-5-2.png)
### `5.3` duals_ID_3
> `MultiReprFilter.optimize`
ID Pairs for order reduction
![Fitttt](plot-5-3.png)