怎么通过频率响应求FIR滤波器的参数?

snowlion1111

就是已知一个FIR滤波器的频率响应,是DB表示的,
怎么通过这些条件求FIR滤波器的传递函数.

liuzhanbin

第一种思路：通过对频率响应的分析，得出通带临界频率Wp，通带内最大衰减Rp，阻带临界频率Ws，阻带内最小衰减Rs，[n,Wn]=butterd(Wp,Ws,Rp,Rs,)，[b,a]=butter(n,Wn)
第二种是还原出输入信号
用下列matlab函数：
A,G,r,a]=autolpc(x,p)
x=x(:);
L=length(x);
r=zeros(p+1,1);
for i=0:p
    r(i+1)=x(1:L-i)'*x(1+i:L);
end
R=toeplitz(r(1:p));
a=R\r(2:p+1);
l=length(a);
i=[1,l];
i1=[2,l+1];
A(1)=1;
A(i1)=-a(i);
%A=[1,-a];
%G=sqrt(sum(A.*r));
G=sqrt(A*r);
%p为阶次，x为输入信号序列，a为预测滤波器系数矢量，A＝[1,-a],r为x的自相关矢量，G为增益常数。

不知道能否可行 一起探讨一下