MATLAB高斯窗函数的源程序？急用

需要实现FIR数字低通滤波器技术指标：Wp= 0．2л， Ws=0．4л， Ap=0．25 dB ，As=50 dB 。对Gauss窗函数进行改进，引入函数w(n)=[1-(nT)*(nT)⼀10]*e-(nT)*(nT)⼀a cos(0.5nT)比较利用Gauss窗和改进窗设计的滤波器的性能。程序自己写的如下： ws=0.4*pi; %通带频率 wp=0.2*pi; wc=0.5*(wp+ws); %截至频率 N=32; %窗长度 n=-N⼀2:N⼀2; alpha=(N-1)⼀2; m=n-alpha+eps; hd=sin(wc*m).⼀(pi*m); %得到理想低通滤波器 T=6⼀32;a=2; %对T，a进行定义 wn=(1-(n.*T.*n.*T)⼀10).*exp((-n.*T.*n.*T).⼀a).*cos(0.5.*n.*T);%改进窗函数 h=hd.*wn; %得到FIR数字滤波器 [H,w1]=freqz(h,1); %求其频率响应 db=20*log10(abs(H)+eps); %得到幅值 plot(w1⼀pi,db);%绘图 axis([0,1,-100,0]); xlabel(✀归一化频率⼀pi✀); ylabel(✀db✀); title(✀ 高斯窗响应✀); 可是运行的图不对~能给出点思路看哪里错了吗