这个滤波器怎么求带宽呢？

HolySaint

1-4/z^16+6/z^32-4/z^48+1/z^64
-------------------------------------------------------
         1-4/z+6/z^2-4/z^3+1/z^4

可以从图中看出来吗？

songzy41

一般对低通滤波器，带宽都是指比最大值下降3分贝处的频率。所以对楼主的滤波器可以先用最大值（0Hz处）归一，找出-3dB的频率，即是该滤波器的带宽。

HolySaint

这个是梳状滤波器的图
是按照多个带通的算还是怎么算？
只知道那个传递函数就可以算出来吧？

按楼上的算法是算到最后下降得还是第一次下降的呢？

yangzj

如果0频是x dB的话,找x-3 dB对应的值就是了.
-3 dB的点就只一个,后面的峰远不只降3 dB了

HolySaint

原帖由 yangzj 于 2008-5-15 12:08 发表
如果0频是x dB的话,找x-3 dB对应的值就是了.
-3 dB的点就只一个,后面的峰远不只降3 dB了

还有就是那个横轴的单位怎么读呢？？
x*2*pi  Hz？？

HolySaint

补充一下

这个是一个四阶的积分器，进行16倍下采样，在经过四阶差分器，最后输出

这个就是整个过程的传递函数

songzy41

原帖由 HolySaint 于 2008-5-15 12:42 发表

还有就是那个横轴的单位怎么读呢？？
x*2*pi  Hz？？

横轴是圆频率，单位是弧度/秒。楼主标出的rad/sample，好象不对，应为rad/second。

HolySaint

原帖由 songzy41 于 2008-5-15 16:37 发表

横轴是圆频率，单位是弧度/秒。楼主标出的rad/sample，好象不对，应为rad/second。

我的输入信号的采样是1024个点

1. t=0:1023;%点数
   
2. 
   
3. %%--------取样测试--------%%
   
4. 
   
5. %d0=fix((sin(2*pi*t/512)+abs(wgn(1,length(t),1))+1)*10000);
   
6. 
   
7. %d0=fix((sin(2*pi*t/512)+1)*10000);
   
8. 
   
9. d0=fix((sin(2*pi*t/512)+square(t*2*pi/12,0.5)+2)*10000);
   
10. 
    
11. %%--------画图比较--------%%
    
12. 
    
13. plot(d0);
    
14. 
    
15. title('周期取点数256');
    
16. 
    
17. %%--------转换数制--------%%
    
18. 
    
19. for i=1:1024
    
20. 
    
21.     b=dec2bin(d0(i),16);
    
22. 
    
23.     bout(i,:)=b;
    
24. 
    
25. end;
    
26. 
    
27. %%--------写入文本--------%%
    
28. 
    
29. fid = fopen('input.txt','wt');
    
30. 
    
31. for j=1:1024
    
32. 
    
33.     fprintf(fid,'%s\n',bout(j,:));
    
34. 
    
35. end;
    
36. 
    
37. fclose(fid);

复制代码

主要是我想证实一下这个滤波器的带宽是否大于24kHz？？

需要把输入改成24k的信号，，看结果吗？？

yangzj

图中是归一化的频率,1对应采样频率的一半.

HolySaint

我这个滤波器使用verilog写的
他的传函如一楼
输入时matlab那个程序生成的
我想验证带宽大于24k

我的输入要怎么给？
下图可以说明吗？？

附下图生成信号代码

1. t=0:47999;%点数
   
2. 
   
3. t1=0:.00001:.47999;
   
4. 
   
5. %%--------取样测试--------%%
   
6. 
   
7. d0=fix((sin(2*pi*t/24000)+1)*10000);
   
8. 
   
9. %%--------画图比较--------%%
   
10. 
    
11. plot(d0);
    
12. 
    
13. title('周期取点数256');
    
14. 
    
15. %%--------转换数制--------%%
    
16. 
    
17. for i=1:48000
    
18. 
    
19.     b=dec2bin(d0(i),16);
    
20. 
    
21.     bout(i,:)=b;
    
22. 
    
23. end;
    
24. 
    
25. %%--------写入文本--------%%
    
26. 
    
27. fid = fopen('input.txt','wt');
    
28. 
    
29. for j=1:48000
    
30. 
    
31.     fprintf(fid,'%s\n',bout(j,:));
    
32. 
    
33. end;
    
34. 
    
35. fclose(fid);

复制代码

[ 本帖最后由 HolySaint 于 2008-5-15 17:14 编辑 ]

HolySaint

滤波前是48000个点
滤波后是3000个点～

[ 本帖最后由 eight 于 2008-5-15 21:12 编辑 ]