www.pudn.com > bfunction.rar > scale.m

function  f=scale(m,w) 
 
%a=1-exp(-j*w); 
%a=a./(-j*w); 
%a=sin(w/2)./(w/2); 
b=w./2; 
x=sin(b)./b; 
y=exp(-i*b); 
a=x.*y; 
 
f1=a.^m; 
%caculater the denominator 
N=B(2*m); 
 
d=0; 
 for k=1-m:-1 
     d=d+2*N(m+k+1)*cos(k*w); 
 end  
 d=sqrt(d+N(m+1)); 
 f=f1./d;