周圈圈 2019-05-05 17:19
浏览 543

matlab我用filter过滤谐波中的低次波,但出图就是有问题?有没有高手帮一下?

clear;
close all
wp=2*pi*2000;ws=2*pi*500;Rp=1;As=20;

[N,wc]=buttord(wp,ws,Rp,As,'s')%计算率波器的阶数和3dB截止频率

[B,A]=butter(N,wc,'high','s');%计算滤波器系统函数分子分母多项式

fk=0:800/512:8000;wk=2*pi*fk;

Hk=freqs(B,A,wk);

figure(1)
subplot(1,3,1);
plot(fk/1000,20*log10(abs(Hk)));
grid on;xlabel('频率(kHz)');ylabel('幅度(dB)');
title('巴特沃斯模拟高通滤波器');
axis([0,4,-35,5])

x=(0:0.000025:0.2);
Signal_Original=sin(2*pi*50*x)+0.1*sin(6000*pi*x)+0.05*sin(8000*pi*x);
subplot(1,3,2);plot(x,Signal_Original);xlabel('时间(s)');ylabel('安倍(A)');

Signal_Filter=filter(B,A,Signal_Original);
subplot(1,3,3);plot(x,Signal_Filter);

图片说明

就是不明白为什么第三个图
我用filter哪里有问题

  • 写回答

0条回答

    报告相同问题?

    悬赏问题

    • ¥65 永磁型步进电机PID算法
    • ¥15 sqlite 附加(attach database)加密数据库时,返回26是什么原因呢?
    • ¥88 找成都本地经验丰富懂小程序开发的技术大咖
    • ¥15 如何处理复杂数据表格的除法运算
    • ¥15 如何用stc8h1k08的片子做485数据透传的功能?(关键词-串口)
    • ¥15 有兄弟姐妹会用word插图功能制作类似citespace的图片吗?
    • ¥200 uniapp长期运行卡死问题解决
    • ¥15 latex怎么处理论文引理引用参考文献
    • ¥15 请教:如何用postman调用本地虚拟机区块链接上的合约?
    • ¥15 为什么使用javacv转封装rtsp为rtmp时出现如下问题:[h264 @ 000000004faf7500]no frame?