blob: fa5fa8af20ea8b45cef4c6aff3871116b2f93826 (
plain)
1
2
3
4
5
6
7
8
9
10
|
// function [phi,psi] = ch_pol(N,epsilon)
// Returns desired characteristic polynomial and numerator
// N = rise time in number of sample times
// epsilon = overshoot as a fraction of ss gain
function [phi,psi] = ch_pol(N,epsilon)
omega = %pi/2/N; r = epsilon^(omega/%pi);
phi = [1 -2*r*cos(omega) r^2];
psi = [1-r*cos(omega) (r^2-r*cos(omega))];
endfunction;
|