请帮忙看下这段程序 谢谢各位

ambelmm

clear all
fid=fopen('D:\学习\声学基础\衰减 成像实验数据\neck\sample_FB_00000.dat','r');
shepindata=fread(fid,inf);
data=shepindata(1:52800);
n=3;
B=zeros(1:(n+1));
A=zeros(1:(n+1));
Wn=[7800000,16200000];
[B,A]=BUTTER(n,Wn);
data_f=filter(B,A,data)

倒数第二行总有问题
这个滤波器是三阶的巴特沃斯带通滤波器 按理说系数应该是4维向量 不知道我前面的定义是不错了
麻烦各位帮忙看下

yangzj

看看butter 的帮助里对参数的说明

ambelmm

看过了 B A表示系数  应是n+1维的向量 是不我的wn和我的数据不符啊
不知道程序本身有什么问题

yangzj

The cutoff frequency Wn must be 0.0 < Wn < 1.0, with 1.0 corresponding to half the sample rate。