电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

MATLAB巴特沃兹滤波器VIP免费

MATLAB巴特沃兹滤波器_第1页
1/8
MATLAB巴特沃兹滤波器_第2页
2/8
MATLAB巴特沃兹滤波器_第3页
3/8
相 位 是 信 号 的 一 个 重 要 参 数 , 所 以 一 般 在 滤 波 中 希 望 实 现 零 相 位 偏 移 或 具 有 线性 相 位 偏 移 。 现 在 大 致 有 两 类 方 法 实 现 零 相 位 偏 移 : 1、将滤 波 器与一 个 全通滤 波 器级联, 用后者的 相 位 响应矫正前者。 需要 附 加全通滤 波 器的 设计运算量以 及额外的 矫正算法 运算量; 2、使用两 个 完全相 同的 滤 波 器, 将通过第一 个 滤 波 器后的 信 号 反转再通过 第二个 滤 波 器, 再将第二个 滤 波 器的 输出序列反转作为整个 系统的 输出。 综上我们采用第二种方 法 , 用 MATLAB 实 现 : clear clc close all %生成一个信号 Fs=1000; t=(0:1/Fs:1); y=sin(2*pi*5*t+pi/4)+sin(2*pi*50*t+pi/5)+sin(2*pi*60*t+pi/12); yy=sin(2*pi*5*t+pi/4); figure(1) subplot(211) plot(t,y),grid on title('输入信号波形') xlabel('t') ylabel('magnitude') subplot(212) plot(t,angle(y)) title('输入信号相位') xlabel('t') ylabel('angle') %设计一个低通滤波器 fp=10/500;%通带截止频率 fs=30/500;%阻带截止频率 wp=3;%通带最大衰减 ws=30;%阻带最大衰减 [n,fc]=buttord(fp,fs,wp,ws); [b,a]=butter(n,fc); figure(4) freqz(b,a); %将生成的信号通过低通滤波器 z=filter(b,a,y); figure(2) subplot(221) plot(t,z) title('通过第一个滤波器后的波形') xlabel('t') ylabel('output') axis([0 1 -1 1]) subplot(222) plot(t,angle(z)) title('通过第一个滤波器后的相位') xlabel('t') ylabel('angle') subplot(223) plot(t,yy) title('理论上过滤后的波形') xlabel('t') ylabel('output') axis([0 1 -1 1]) subplot(224) plot(t,angle(yy)); title('原相位') xlabel('t'); ylabel('angle') zf=fliplr(z); %对经过第一个滤波器后的信号序列反转 zfz=filter(b,a,zf); %对经过第一个滤波器后的信号序列反转再通过滤波器 zz=fliplr(zfz); %对经过两次滤波器后的输出反转作为最后输出 figure(3) subplot(221) plot(t,zz) title('通过两个滤波器后的波形') xlabel('t') ylabel('output') axis([0 1 -1 1]) subplot(222) plot(t,angle(zz)) title('通过两个滤波器后的相位') xlabel('t') ylabel('angle')...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

MATLAB巴特沃兹滤波器

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部