ousta_fod.m

来自「高等引用数学的matlab求解的配套程序」· M 代码 · 共 7 行

   function G=ousta_fod(r,N,w_L,w_H)
   mu=w_H/w_L; k=-N:N; 
   w_kp=(mu).^((k+N+0.5-0.5*r)/(2*N+1) )*w_L;
   w_k=(mu).^((k+N+0.5+0.5*r)/(2*N+1) )*w_L;
   K=(mu)^(-r/2)*prod(w_k./w_kp); 
   G=tf(zpk(-w_kp',-w_k',K));