原创力文档- 1、本文档共20页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
word
word
PAGE / NUMPAGES
word
数字信号处理
班级：13通信1班
同组人：
骆燕 1311011320
X逸媛 1310200732
鲁鑫 1310200a33
李秀婷1310211221
何思雨 1310313113
数字信号处理课程设计
1、试设计一个模拟低通Butterworth滤波器
解：经计算可知低通Butterworth滤波器的阶数为3
用Matlab软件计算的程序如下：
OmegaP=0.2*pi;
OmegaS=0.3*pi;
Rp=7;
As=16;
N=ceil(log10((10.^(0.1*abs(Rp))-1)./(10.^(0.1*abs(As))-1))/(2*log10(OmegaP/OmegaS)))
运行程序后出来的结果如下图:
使用Matlab软件仿真的程序如下：
Wp=input(Normalized passband edge =);
Ws=input(Normalized stopband edge =);
Rp=input(Passband ripple in dB =);
Rs=input(Minimum stopband attenuation in dB =);
[N,Wn]=buttord (Wp,Ws,Rp,Rs);
[b,a]=butter(N,Wn);
[h,omega]=freqz(b,a,512);
plot(omega/pi,20*log10(abs(h)));
grid;
xlabel(\omega/\pi);
ylabel(Gain dB);
title(IIR Butterworth Lowpass Filter);
得出的仿真波形如下：
2、设Chebyshev模拟滤波器的技术指标?c=2??2kHz, ?s=2??4kHz,在?=?s处的衰减小于-15dB,通带波动参数，?2 =0.2，求其阶数，并绘出幅频响应分贝图。
解：经计算可知切比雪夫模拟滤波器的阶数为4
用matlab计算的程序如下：
clc;
clear all;
OmegaP=2*pi*2000;
OmegaS=2*pi*4000;
Rp=0.17;
As=15;
g=sqrt((10^(As/10)-1)/(10^(Rp/10)-1));
OmegaR=OmegaS/OmegaP;
N=ceil(log10(g+sqrt(g*g-1))/log10(OmegaR+sqrt(OmegaR*OmegaR-1)))
其计算结果显示如下：
使用matlab软件仿真程序如下：
clc;
clear all;
OmegaP=2*pi*2000;
OmegaS=2*pi*4000;
Rp=0.17;
As=15;
g=sqrt((10^(As/10)-1)/(10^(Rp/10)-1));
OmegaR=OmegaS/OmegaP;
N=ceil(log10(g+sqrt(g*g-1))/log10(OmegaR+sqrt(OmegaR*OmegaR-1)));
OmegaC=OmegaS;
[z0,p0,k0]=cheb2ap(N,As);
a0=real(poly(p0));
aNn=a0(N+1);p=p0*OmegaC;a=real(poly(p));
aNu=a(N+1);b0=real(poly(z0));M=length(b0);
bNn=b0(M);z=z0*OmegaC;b=real(poly(z));
bNu=b(M);k=k0*(aNu*bNn)/(aNn*bNu);
b=k*b
w0=[OmegaP,OmegaS];
[H,w]=freqs(b,a);
Hx=freqs(b,a,w0);
dbHx=-20*log10(abs(Hx)/max(abs(H)))
plot(w/(2*pi)/1000,20*log10(abs(H)));
xlabel(f(kHz));ylabel(dB);axis([-1,12,-55,1]);
set(gca,xtickmode,manual,xtick,[0,1,2,3,4,5,6,7,8,9]);
set(gca,ytickmode,manual,ytick,[-50,-40,-30,-20,-10,0]);
grid;
显示结果如下图：
3、用矩形窗、汉宁窗和布莱克曼窗设计FIR低通滤波器。
程序为：
clear all;
N=11;
n=1:N;
wdhn=hamming(N);
wdbl=blackman(N);w=boxcar(N);
plot(n,[wdhn,wdbl,w])
显示结果如下图：
4、用matlab设计一个LPF FIR wp=0.24π ws=0.28
文档评论（0）