% setting up filter parameters and the svd of the input-to-state response
thetamid=1.325;
[A,B]=cjordan(5,0.88*exp(thetamid*1i));
% obtaining state statistics
R=dlsim_complex(A,B,y');
sv=Rsigma(A,B,th);
plot(th(:),sv(:));