givens.m

This is mat lab code for adaptive lattice filters.

function [c,s] = givens(a,b) % returns cosine and sine parameters


if b == 0
   c=1;
   s=0;
elseif  a == 0
   c=0;
   s=1;
else
   rho=b/a; 
   c=1/sqrt(1+abs(rho)^2); 
   s=rho*c;  
end