《DSP using MATLAB》示例Example7.11
2024-08-29 19:44:38
代码:
M = 45; As = 60; n = [0:1:M-1];
beta = 0.1102*(As - 8.7)
%beta = 0.1102*(As - 8.7) + 0.3 w_kai = (kaiser(M, beta))'; wc1 = pi/3; wc2 = 2*pi/3; hd = ideal_lp(wc1, M) + ideal_lp(pi, M) - ideal_lp(wc2, M); h = hd .* w_kai; [db, mag, pha, grd, w] = freqz_m(h, [1]); %Plot figure('NumberTitle', 'off', 'Name', 'Exameple 7.11')
set(gcf,'Color','white'); subplot(2,2,1); stem(n, hd); axis([0 M-1 -0.2 0.8]); grid on;
xlabel('n'); ylabel('hd(n)'); title('Ideal Impulse Response'); subplot(2,2,2); stem(n, w_kai); axis([0 M-1 0 1.1]); grid on;
xlabel('n'); ylabel('w(n)'); title('Kaiser Window'); subplot(2,2,3); stem(n, h); axis([0 M-1 -0.2 0.8]); grid on;
xlabel('n'); ylabel('h(n)'); title('Actual Impulse Response'); subplot(2,2,4); plot(w/pi, db); axis([0 1 -80 10]); grid on;
xlabel('frequency in \pi units'); ylabel('Decibels'); title('Magnitude Response in dB');
运行结果:
参数β =5.6533,最小的阻带衰减小于60dB。很明显,我们增加β 就能使得衰减达到60dB,此时β =5.9533。
将上面代码的第2个β 公式前面注释去掉,再次运行,结果如下:
此时满足设计要求。
最新文章
- html 超文本标记语言
- JS中变量名和函数名重名
- Swift3.0都有哪些变化
- 2016-11-10:win7下VMware虚拟机中CentOS6.5网络配置
- java中会存在内存泄漏吗,请简单描述。
- [codevs2070]爱情之路
- Eclipse开发过程中个VM Arguments的设置
- Javascript 基础知识学习--javascript中的参数传递都是按值传递的
- BZOJ3630 : [JLOI2014]镜面通道
- Mtk Android编译命令
- LeetCode49 Group Anagrams
- 高健壮性css---Float详细
- Errore HTTP 404.2 - Not Found"; IIS 7.5 请求的内容似乎是脚本,因而将无法由静态文件处理程序来处理
- 【STL】算法 — partial_sort
- (转载)ANDROID STRINGS.XML 中的特殊字符转义
- .NET Core微服务之开源项目CAP的初步使用
- android apk打包编译好的so
- The Little Prince-12/01
- Informatica学习:3、用户创建与权限管理
- innodb next-key lock引发的死锁