Implement function for efficient calculation and visualization of right-most spectrum envelope.
Envelope is defined via
$$
\begin{aligned}
s \geq c_d \\
s \in \mathcal{Z} \left(h\left(s\right)\right)
\end{aligned}
\implies |s| \leq \rho(s)\sum\limits_{k=0}^{m_A} \lVert A_k\rVert_2 e^{-\Re(s)\tau_k}
$$
where
$$
\rho(s) = \max\limits_{\vec{\theta} \in [0,2\pi]^n}\left\lVert \sum\limits_{i=1}^{m_H} H_i~e^{-\Re(s)\tau_i}~e^{j\theta_i} \right\lVert_2^{-1}
$$
for neutral system, for retarded $c_d = -\infty$ and $\rho(s)=1$
Re-use gridding functions for theta.
Implement function for efficient calculation and visualization of right-most spectrum envelope.
Envelope is defined via
where
for neutral system, for retarded$c_d = -\infty$ and $\rho(s)=1$
Re-use gridding functions for theta.