- 1、本文档共51页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
关于几个常见数字滤波器的设计
实验一 信号、系统及系统响应
编写常用信号的程序:
(1):单位采样序列δ(n)。
程序如下:
function [x,n]=impseq(n0,ns,nf)
n0=0;
ns=-5;
nf=5;
n=[ns:nf];
x=[(n-n0)==0]
stem(n,x)
grid
title(单位采样序列)
xlabel(n)
ylabel(x(n))
运行结果如下:
:单位阶跃序列u(n)。
程序如下:
function [x,n]=stepseq(n0,ns,nf)
n0=0;
ns=-5;
nf=5;
n=[ns:nf];
x=[(n-n0)=0]
stem(n,x)
grid
title(单位阶跃序列)
xlabel(n)
ylabel(u(n))
运行结果如下:
(3):单位矩形序列。
程序如下:
function [x,n]=juxingseq(n0,ns,nf)
n0=4;
ns=0;
nf=9;
n=[ns:nf];
x=[((n-ns)=0)((n-nf)=0)]
stem(n,x)
grid
title(单位矩形序列)
xlabel(n)
ylabel(R(n))
end
运行结果如下:
正弦序列。
程序如下:
function [x,n]=sin(w,n)
w=0.2*pi;
n=(0:0.2:10);
x=sin(w*n);
stem(n,x)
title(正弦序列)
xlabel(n)
ylabel(x(n))
运行结果如下:
(5):。
程序如下:
function x=exp(n)
n=(0:0.0001:0.01);
x=exp(-1000*abs(n))
stem(n,x,.)
运行结果如下:
(6):。
程序如下:
function x=rs(n)
n=(0:0.1:10);
A=1;c=2;w=2*pi;
x=A*exp(-c*n).*sin(w*n)
stem(n,x,.)
运行结果如下:
2、认真复习采样理论、 离散信号与系统、 线性卷积、 序列的傅里叶变换及性质等有关内容, 阅读本实验原理与方法。
观察并分析采用不同频率时,对函数的频谱影响。
(a):以,对其进行采样得到。
(b):以,对其进行采样得到。
(c):对、采用理想内插函数重建原始信号。
程序如下:
function x=exp(n)
n1=(0:0.0002:0.01);
x=exp(-1000*abs(n1))
subplot(2,1,1)
stem(n1,x,.)
grid on
n2=(0:0.001:0.01);
x=exp(-1000*abs(n2))
subplot(2,1,2)
stem(n2,x,.)
grid on
运行结果如下:
实现两序列的卷积。利用y=conv(x,h).例如:求解
试求:.
程序如下:
x=[5,9,3,6,-8];
h=[18,7,5,20,11,14,9];
y=conv(x,h)
n=-4:6;
stem(n,y,b)
hold on
grid on
运行结果如下:
(3)系统单位脉冲响应序列产生子程序。 本实验要用到两种FIR系统。
a. ha(n)=R10(n);
b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3)
hb(n)信号的程序如下:
function [x,n]=impseq(n0,ns,nf)
n0=0;
ns=-5;
nf=5;
n=[ns:nf];
x0=[(n-n0)==0]
x1=2.5*[(n-1)==0]
x2=2.5*[(n-2)==0]
x3=[(n-3)==0]
x=x0+x1+x2+x3
stem(n,x)
grid
title(单位采样序列)
xlabel(n)
ylabel(x(n))
运行结果如下:
3.时域离散信号、 系统和系统响应分析。
观察信号xb(n)和系统hb(n)的时域和频域特性; 利用线性卷积求信号xb(n)通过系统hb(n)的响应y(n), 比较所求响应y(n)和hb(n)的时域及频域特性, 注意它们之间有无差别, 绘图说明, 并用所学理论解释所得结果。
b. 观察系统ha(n)对信号xc(n)的响应特性。
程序如下:
w=(-3*pi:0.001:3*pi)+eps;
Xa=1;
Xb=1+2.5*exp(-j*w)+2.5*exp(-2*j*w)+exp(-3*j*w);
subplot(2,1,1)
plot(w/pi,abs(Xa),b)
xlabel(w/pi)
ylabel(|Xa|)
title(Xa频域特性)
hold on
grid on
subpl
文档评论(0)