回 帖 发 新 帖 刷新版面

主题:求大家帮忙看看这个程序吧 谢谢

是一个巴特沃斯滤波程序,下面是部分数据
     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=[2*pi*4875000,2*pi*10125000];
[B,A]=butter(3,wn,'s');
data_f=filter(B,A,data);
datah_f=reshape(data_f,10,10)[/b]

结果 
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% 
为什么后面会出现那么多无穷大呢 我遍的哪有问题啊

回复列表 (共1个回复)

沙发

datah=reshape(data,10,10);
data_f=filter(B,A,data);
那这个datah有什么用??

我来回复

您尚未登录,请登录后再回复。点此登录或注册