资源描述：
《实验1离散系统时域分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验1离散系统时域分析一.离散卷积计算二.程序如下：n=0:19;x=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1];%x(n)的输入n=0:19;h=0.5*sin(0.5*n);%h(n)的输入[y]=conv(x,h);%公式法求卷积ny=length(y);n=0:ny-1;stem(n,y);xlabel('n');ylabel('y(n)');%图形显示二.理解卷积y(n)=x(n)*h(n)，并采用conv命令求解的差分方程y(n)=x(n)+6x(
2、n-1)-3x(n-2)+5x(n-3)-8x(n-4)-7x(n－5)其中，程序如下：1.单位取样序列Impseq：function[x,n]=impseq(n0,n1,n2)n=[n1:n2];x=[(n-n0)==0];2.差分方程求解：n=[0:10];x=impseq(0,0,10);%单位取样序列n=0:10;h=[1,6,-3,5,-8,-7];%差分方程系数的输入[y]=conv(x,h);%x与h做卷积运算求yny=length(y);n=0:ny-1;stem(n,y);xlab
3、el('n');ylabel('y(n)');%图像显示三．采用conv命令求解的差分方程y(n)=x(n)+6x(n-1)-3x(n-2)+5x(n-3)-8x(n-4)-7x(n－5)其中，x(n)=s(n-2)+2s(n-1)+s(n),y(n)=0当n<0程序如下：1.单位取样序列Impseq：同上2.差分方程求解：n=[0:10];x=impseq(2,0,10)+2*impseq(1,0,10)+impseq(0,0,10);%x(n)的输入n=[0:10];h=[1,6,-3,5,-8
4、,-7];%差分方程系数的输入[y]=conv(x,h);%x与h做卷积运算求yny=length(y);n=0:ny-1;stem(n,y);xlabel('n');ylabel('y(n)');%图像显示实验2利用FFT和IFFT计算线性卷积利用FFT实现以下两序列的线性卷积程序如下：n=0:19;x=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1];%x(n)的输入h=0.5*sin(0.5*n);%h(n)的输入;L=length(x)+length(h)-
5、1;X=fft(x,L);H=fft(h,L);y=ifft(X.*H);subplot(3,1,1);stem(n,x);ylabel('序列x(n)');xlabel('n');subplot(3,1,2);stem(n,h);ylabel('序列h(n)');xlabel('n');n=0:L-1;subplot(3,1,3);stem(n,real(y));ylabel('卷积结果y(n)');xlabel('n');实验3IIR滤波器设计实验内容(1)设计一个数字低通滤波器．设计指标如下：
6、通带内幅度衰减人大于1dB阻带内幅度衰减大于15dB(2)编制计算设计的数字滤波器幅度特性与相位特性程序，并进行实验验证；(3)编制实现该数字滤波器程序并且模拟实现数字滤波①分别让满足所设计滤波器的通带、过滤带、阻带频率特性的正弦波通过滤波器，验证滤波器性能。②改变正弦抽样时间，验证数字低通滤波器的模拟截止频率是抽样时间的函数。双线性变换设计数字巴特沃斯低通滤波器1.主程序如下：wp=0.2*pi;ws=0.35*pi;Rp=1;As=15;%数字滤波器指标T=1;Fs=1/T;%转换成模拟域指标o
7、megap=(2/T)*tan(wp/2);%预畸变截止频率omegapomegas=(2/T)*tan(ws/2);%预畸变截止频率omegas[cs,ds]=afd_buttap(omegap,omegas,Rp,As);%模拟巴特沃斯滤波器的计算%双线性变换[b,a]=bilinear(cs,ds,Fs);%(模拟到数字滤波变换的双线性变化方法)[C,B,A]=dir2cas(b,a);%直接型到级联型的转换函数[db,mag,pha,grd,w]=freqz_m(b,a);%计算离散系统频率
8、响应的函数%图形显示subplot(2,2,1);plot(w/pi,mag);ylabel('
9、幅度
10、');xlabel('以pi为单位的频率');title('幅度响应');grid;axis([0101.5]);subplot(2,2,3);plot(w/pi,db);title('幅度响应(dB)');grid;xlabel('以pi为单位的频率');ylabel('对数幅度/dB');axis([01-800]);subplot(2,2,2);plot(
