Question

matlab低通滤波程序,怎么用

看到很多类似的程序，比如下面这个
Fp=3500;Fs=4500;
Wp=2*pi*Fp;Ws=2*pi*Fs;
[N,Wn]=buttord(Wp,Ws,0.5,30,'s');
[b,a]=butter(N,Wn,'s');
wa=0:(3*Ws)/511:3*Ws;
h=freqs(b,a,wa);
plot(wa/(2*pi),20*log10(abs(h)));grid;
第一个问题是freqs是做什么用的？第二是这个程序的输入是什么（哪个量），滤波器不是要有输入吗？

Answer 1

h = freqs(b, a, w) 根据系数向量计算返回滤波器的复频域响应。

调用这个低通滤波器时，使用下面的函数

sf=filter(a,b,s); %s为需滤波的数据，sf经过你设计的低通滤波器以后的新数据

追问

T是我的需处理数据。我输入命令后提示错误
Fp=2100;Fs=2800;
Wp=2*pi*Fp;Ws=2*pi*Fs;
[N,Wn]=buttord(Wp,Ws,10,30,'s');
[b,a]=butter(N,Wn,'s');
sf=filter(a,b,T);
英文的意思我记得是滤波器的第一个系数需要是非零。怎么回事？
（补充一个情况：我把b和a调换了一下，倒是能运行了，但是处理结果是一堆NaN）

追答

你滤波器设计的有问题。
[N,Wn]=buttord(Wp,Ws,10,30,'s');
在matlab中输入 doc buttord
看下这个函数，wp和ws的范围是0到1之间，需要wp和ws要归一化的

Answer 2

freqs() 是把时域的资料转到频域
wa 就是你的输入资料
[b, a] 是低通滤波器参数
h 是输出
程式中以 20*log10的方式绘图

Answer 3

h = freqs(b, a, w) 根据系数向量计算返回滤波器的复频域响应。

调用这个低通滤波器时，使用下面的函数

sf=filter(a,b,s);  %s为需滤波的数据，sf经过你设计的低通滤波器以后的新数据

Fp=2100;Fs=2800;
Wp=2*pi*Fp;Ws=2*pi*Fs;
[N,Wn]=buttord(Wp,Ws,10,30,'s');
[b,a]=butter(N,Wn,'s');
sf=filter(a,b,T)