关于几个常见数字滤波器的设计.docVIP

实验一 信号、系统及系统响应 编写常用信号的程序： （１）：单位采样序列δ（ｎ）。 程序如下： 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)) 运行结果如下： ：单位阶跃序列ｕ（ｎ）。 程序如下： 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)) 运行结果如下： （３）：单位矩形序列。 程序如下： 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、认真复习采样理论、 离散信号与系统、 线性卷积、 序列的傅里叶变换及性质等有关内容， 阅读本实验原理与方法。 观察并分析采用不同频率时，对函数的频谱影响。 （ａ）：以，对其进行采样得到。 　　（ｂ）：以，对其进行采样得到。 （ｃ）：对、采用理想内插函数重建原始信号。 程序如下： 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 运行结果如下： 实现两序列的卷积。利用ｙ＝ｃｏｎｖ（ｘ，ｈ）．例如：求解 试求：． 程序如下： 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