guanzhengting
2007-06-16, 16:10
我设计了数字滤波器的系数求法，但是系数文件自己建了，可就是没系数在里面，哪要修改啊？？？
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)
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)