注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

FIR带阻滤波器的matlab实现  

2013-01-07 20:58:45|  分类: 编程 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
clf;N=300;Fs=1000;         %数据总数和采样频率
Order=100;    %滤波器阶数
n=[0:N-1];t=n/Fs;     %时间序列
wn=[100 200]/(Fs/2); %边界频率转换为归一化频率
b=fir1(Order,wn,'stop');          % 设计 100 阶的带阻滤波器
figure(1)
[H,f]=freqz(b,1,512,Fs);
subplot(2,1,1),plot(f,20*log10(abs(H)))
xlabel('频率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('频率/Hz');ylabel('相位/^o');grid on;
f1=50;f2=150;  %输入信号频率
x=sin(2*pi*f1*t)+sin(2*pi*f2*t);       %输入信号
y=fftfilt(b,x,256);% 对数据进行滤波
figure(2)
subplot(2,1,1);plot(t,x); grid on;                %绘制输入信号
title('输入信号');
subplot(2,1,2);plot(t,y);                   %绘制输出信号
hold on;plot(Order/2/Fs*ones(1,2),ylim, 'r')    %绘制延迟到的时刻
title('输出信号');
grid on;  xlabel('时间/s')
FIR带阻滤波器的matlab实现 - 阿英 - Mr.Right
 
  评论这张
 
阅读(645)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016