【求助】程序求助啊！！滤波器系数问题

[guanzhengting]

下面是我设计的滤波器系数程序，可在系数文件中没有系数，哪需要改下？？

fp=1500;fst=2250;fs=16000; %输入设计指标
wp=2*fp/fs; %求归一化数字通带截至频率
ws=2*fst/fs; %求归一化数字通带起始频率
deltaw=ws-wp; %求过渡带宽
NO=ceil(6.6/deltaw); %求窗口长度
N=NO+mod(NO+1,2); %确保窗口长度N为奇数
n=N-1; %求出滤波器的介数n
wn=(ws+wp)/2; %求滤波器的截止频率
b=fir1(n,wn); %利用fir1函数求出滤波器的系数
b1=round(b*10000); %将系数用蒸熟表示
fn=fopen('COEF_FIR.inc','w'); %将系数存放在文件COEF_FIR.inc中
fclose(fn)

[foxstar]

看看下面的行不?
fp=1500;fs=2250;Fs=16000; %输入设计指标
wp=2*fp/Fs; %求归一化数字通带截至频率
ws=2*fs/Fs; %求归一化数字通带起始频率
deltaw=ws-wp; %求过渡带宽
NO=ceil(6.6/deltaw); %求窗口长度
N=NO+mod(NO+1,2); %确保窗口长度N为奇数
n=N-1; %求出滤波器的介数n
wn=(ws+wp)/2; %求滤波器的截止频率
b=fir1(n,wn); %利用fir1函数求出滤波器的系数
b1=round(b*10000); %将系数用蒸熟表示
fn=fopen('COEF_FIR.inc','w'); %将系数存放在文件COEF_FIR.inc中
fprintf(fn,'the 70th coeff is: \n');

fprintf(fn,'%5d ',b1);

fclose(fn);