function mainfun()
Wp=0.2*pi;
Ws=0.4*pi;
tr_width=Ws-Wp; %过渡带宽度
N=ceil(6.6*pi/tr_width)+1 %滤波器长度
n=0:1:N-1;
Wc=(Ws+Wp)/2; %理想低通滤波器的截至频率
hd=ideal_lp(Wc,N)'; %理想低通滤波器的单位冲激响应
w_ham=(hamming(N)) %海明窗
h=hd.*w_ham; %截取得到实际的单位脉冲响应
[db,mag,pha,w]=freqz_m4(h,[1]); %计算实际滤波器的幅度响应
delta_w=2*pi/1000;
Ap=-(min(db(1:1:Wp/delta_w+1))) %实际通带纹波
As=-round(max(db(Ws/delta_w+1:1:501))) %实际阻带纹波
subplot(221)
stem(n,hd) %火柴棍图
title('理想单位脉冲响应hd(n)')
subplot(222)
stem(n,w_ham)
title('海明窗w(n)')
subplot(223)
stem(n,h)
title('实际单位脉冲响应h(n)')
subplot(224)
plot(w/pi,db)
title('幅度响应(dB)')
axis([0,1,-100,10])
msgbox('MATLAB编程答疑,请加QQ: 993878382','MATLAB答疑','help')
web http://993878382.qzone.qq.com -browser
end
%本程序用到的自定义函数
function hd=ideal_lp(Wc,N)
alpha= (N-1)/2;
n=0:1:N-1;
m=n-alpha+eps;
hd=sin (Wc*m)./(pi*m);
end
function[db,mag,pha,w]=freqz_m4(b,a)
[H,w]=freqz(b,a,1000,'whole');
H=(H(1:1:501));
w=(w(1:1:501));
mag=abs(H);
db=20*log10((mag+eps)/max(mag));
pha=angle(H);
end
没有合适的资源?快使用搜索试试~ 我知道了~
matlab源码之设计的简单滤波器程序源代码.zip
共1个文件
m:1个
需积分: 0 0 下载量 14 浏览量
2023-07-14
12:31:43
上传
评论
收藏 871B ZIP 举报
温馨提示
matlab源码之设计的简单滤波器程序源代码 matlab源码之设计的简单滤波器程序源代码 matlab源码之设计的简单滤波器程序源代码 matlab源码之设计的简单滤波器程序源代码 matlab源码之设计的简单滤波器程序源代码
资源推荐
资源详情
资源评论
收起资源包目录
matlab源码之设计的简单滤波器程序源代码.zip (1个子文件)
mainfun.m 1KB
共 1 条
- 1
资源评论
大学生资源网
- 粉丝: 137
- 资源: 1334
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 【Unity尾巴动画插件】Tail Animator 轻松创建和控制角色的尾巴动画
- RabbitMQ 延时执行的功能插件
- Java数组反转技巧:保持元素原始类型与代码实现
- 【Unity纹理生成和材料编辑工具】Surforge
- 【日常办公必须工具】文件管理+批量移动文件+实用工具+软件开发+windows必备
- 基于stm32的六轴机械臂控制+openmv颜色识别-识别不同的物块分放(源码+文档说明)
- 深入浅出Pandas:利用Python进行数据处理与分析 (李庆辉)
- 【Unity烟雾特效插件】VFX Graph - Stylized Smoke - Vol. 1 高质量的烟雾特效
- 2019-灵活就业数据集.dta
- 2017-灵活就业数据集.dta
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功