scale.m

在matlab环境下

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;