完善资料让更多小伙伴认识你，还能领取20积分哦， 立即完善>
3天内不再提示
16267
|
本帖最后由 hanl 于 2011-7-21 10:07 编辑 我用C语言实现一个butterworth低通滤波器，现在遇到问题，please help me！ 实现过程： 1、采样频率10KHz，通带边界频率2KHz，通带波纹小于1dB，阻带边界频率3KHz，阻带衰减大于40dB。 2、用matlab求出滤波器系数b，a
3、差分方程： y(0) = b(0)*x(0) + b(1)*x(1) +...+ b(9)*x(9) - {a(1)*y(1) + a(2)*y(2) +...+ a(9)*y(9)} 4、程序：
5、测试结果 a、无输入 
（单位：0.1ms） b、输入1KHz正弦波，峰峰值 4V
（单位：0.1ms）
到此，不知道哪里出了问题，请高人赐教 
0
|
|
相关推荐
6个回答
|
|
|
|
|
|
|
|
火奴小子 发表于 2013-5-17 14:41 我那个就是AY=BX的形式，系数没搞反，请看仔细 |
|
|
|
|
|
|
|
|
|
|
|
你好，问题找到吗？在学习数字滤波。请指教。
|
|
|
2016-7-28 14:43:01
评论
|
|
|
差分方程是Z域的，你要求的是时域的。你如果要用差分方程就要Z逆变换，要不就用卷积。
|
|
|
|
|
|
滤波函数 0.15
H(z) = ---------------------- 1 - 0.8*z^-1 1. b = 0.15 ; 2. a = [1 -0.8]; 3. 4. n = [0:100]; 5. x = 2*sin(0.05*pi*n) + 2*randn(1, 101); %滤波前序列 6. 7. imp = [1; zeros(100, 1)]; 8. h = filter(b, a, imp); % h为冲激函数的响应（时域） 9. w = conv(x, h); % 做卷积 10. yc = w(1:101); %滤波后序列 11. z= filter(b, a, x); % 进行滤波 12. xlabel('n'); 13. ylabel('x y yc'); 14. plot(n, x,'r',n,z ,'b', n,yc,'m'); 15. grid; 思想很简单，首先是用冲击函数 Imp = [1 ; zeros(100,1)] 取求出系统的冲激响应，然后再与输入进行卷积 运算结果： |
|
|
|
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问，请使用“评论”功能。
基于 DSP5509 进行数字图像处理中 Sobel 算子边缘检测的硬件连接电路图
1713 浏览 0 评论
631 浏览 0 评论
普中科技F28335开发板中，如何使用aic23播放由代码生成的正弦波
2132 浏览 0 评论
2931 浏览 1 评论
1107 浏览 1 评论
/7 
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 （电路图） 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证：合字B2-20210191
工商网监
湘ICP备2023018690号
淘帖
3055
