求各位帮忙看看这个滤波程序哪有问题啊
是一个巴特沃斯滤波程序,下面是部分数据0 0 0 0 0 255 103 238 239 13
0 0 0 0 0 255 255 255 255 0
0 0 0 0 0 252 49 17 9 253
0 0 0 0 0 255 0 0 0 255
0 0 0 0 0 5 98 7 8 240
0 0 0 0 0 0 0 0 0 255
0 0 0 0 255 31 237 17 245 7
0 0 0 0 255 0 255 0 255 0
0 0 0 0 253 37 206 242 253 7
0 0 0 0 255 0 255 255 255 0
横着数从1到100
我编的滤波程序
clear all
fid=fopen('D:\学习\声学基础\衰减 成像实验数据\neck\sample_FB_00000.dat','r+');
shepindata=fread(fid,inf);
data=shepindata(1:100);
datah=reshape(data,10,10);
wn=;
=butter(3,wn,'s');
data_f=filter(B,A,data);
datah_f=reshape(data_f,10,10)
结果
datah_f =
1.0e+305 *
Columns 1 through 8
0 0 0 0 0 -0.0000 0.0000 0.0000
0 0 0 0 0 -0.0000 -0.0000 -0.0000
0 0 0 0 0 0.0000 0.0000 -0.0000
0 0 0 0 0 0.0000 0.0000 0.0000
0 0 0 0 0 -0.0000 -0.0000 0.0000
0 0 0 0 0 0.0000 -0.0000 -0.0000
0 0 0 0 0 0.0000 0.0000 -0.0000
0 0 0 0 0 -0.0000 0.0000 0.0000
0 0 0 0 0 -0.0000 -0.0000 0.0000
0 0 0 0 0.0000 0.0000 -0.0000 -0.0000
Columns 9 through 10
0.0000 NaN
0.0000 NaN
-0.0000 NaN
-0.0000 NaN
0.0000 NaN
1.8489 NaN
-Inf NaN
NaN NaN
NaN NaN
NaN NaN
我们要求中心频率7.5MHZ,带宽70%
为什么后面会出现那么多无穷大呢 我遍的哪有问题啊
[ 本帖最后由 eight 于 2007-7-23 11:44 编辑 ] 结果已经没有信号可言，考虑一下2种原因：
（1）巴特沃斯滤波器应该是针对频域信号，似乎不应该直接对原始数据直接操作
（2）滤波器是否有问题？
[ 本帖最后由 花如月 于 2007-7-21 14:55 编辑 ]
页:
[1]