数字信号处理上机实验答案(全)1

第十章上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。

上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。

本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。

实验一 系统响应及系统稳定性。

实验二 时域采样与频域采样。

实验三 用FFT 对信号作频谱分析。

实验四 IIR 数字滤波器设计及软件实现。

实验五 FIR 数字滤波器设计与软件实现实验六 应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。

建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR 数字滤波器设计及软件实现在。

学习完第六章进行；实验五在学习完第七章后进行。

实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。

10.1 实验一: 系统响应及系统稳定性1.实验目的（1）掌握 求系统响应的方法。

（2）掌握时域离散系统的时域特性。

（3）分析、观察及检验系统的稳定性。

2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。

已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。

在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。

也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。

系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。

重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。

系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。

或者系统的单位脉冲响应满足绝对可和的条件。

数字信号处理实验答案第十章上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。

上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。

本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。

实验一系统响应及系统稳定性。

实验二时域采样与频域采样。

实验三用FFT对信号作频谱分析。

实验四IIR数字滤波器设计及软件实现。

实验五FIR数字滤波器设计与软件实现实验六应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。

建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。

学习完第六章进行；实验五在学习完第七章后进行。

实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。

10.1 实验一: 系统响应及系统稳定性1.实验目的（1）掌握求系统响应的方法。

（2）掌握时域离散系统的时域特性。

（3）分析、观察及检验系统的稳定性。

2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。

已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。

在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB语言的工具箱函数filter函数。

也可以用MA TLAB语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。

系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。

重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。

系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。

或者系统的单位脉冲响应满足绝对可和的条件。

系统的稳定性由其差分方程的系数决定。

数字信号处理实验答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】实验一熟悉Matlab环境一、实验目的1.熟悉MATLAB的主要操作命令。

2.学会简单的矩阵输入和数据读写。

3.掌握简单的绘图命令。

4.用MATLAB编程并学会创建函数。

5.观察离散系统的频率响应。

二、实验内容认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。

在熟悉了MATLAB基本命令的基础上，完成以下实验。

上机实验内容：（1）数组的加、减、乘、除和乘方运算。

输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B并用stem语句画出A、B、C、D、E、F、G。

clear all;a=[1 2 3 4];b=[3 4 5 6];c=a+b;d=a-b;e=a.*b;f=a./b;g=a.^b;n=1:4;subplot(4,2,1);stem(n,a);xlabel('n');xlim([0 5]);ylabel('A');subplot(4,2,2);stem(n,b);xlabel('n');xlim([0 5]);ylabel('B');subplot(4,2,3);stem(n,c);xlabel('n');xlim([0 5]);ylabel('C');subplot(4,2,4);stem(n,d);xlabel('n');xlim([0 5]);ylabel('D');subplot(4,2,5);stem(n,e);xlabel('n');xlim([0 5]);ylabel('E');subplot(4,2,6);stem(n,f);xlabel('n');xlim([0 5]);ylabel('F');subplot(4,2,7);stem(n,g);xlabel('n');xlim([0 5]);ylabel('G');（2）用MATLAB实现下列序列：a) x(n)= 0≤n≤15b) x(n)=e+3j)n 0≤n≤15c) x(n)=3cosπn+π)+2sinπn+π) 0≤n≤15(n)=x(n+16),绘出四个周期。

第十章 上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。

上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。

本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。

实验一 系统响应及系统稳定性。

实验二 时域采样与频域采样。

实验三 用FFT 对信号作频谱分析。

实验四 IIR 数字滤波器设计及软件实现。

实验五 FIR 数字滤波器设计与软件实现实验六 应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。

建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR 数字滤波器设计及软件实现在。

学习完第六章进行；实验五在学习完第七章后进行。

实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。

10.1 实验一: 系统响应及系统稳定性1.实验目的（1）掌握 求系统响应的方法。

（2）掌握时域离散系统的时域特性。

（3）分析、观察及检验系统的稳定性。

2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。

已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。

在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。

也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。

系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。

重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。

系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。

或者系统的单位脉冲响应满足绝对可和的条件。

第十章 上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。

上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。

本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。

实验一 系统响应及系统稳定性。

实验二 时域采样与频域采样。

实验三 用FFT 对信号作频谱分析。

实验四 IIR 数字滤波器设计及软件实现。

实验五 FIR 数字滤波器设计与软件实现实验六 应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。

建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR 数字滤波器设计及软件实现在。

学习完第六章进行；实验五在学习完第七章后进行。

实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。

10.1 实验一: 系统响应及系统稳定性1.实验目的（1）掌握 求系统响应的方法。

（2）掌握时域离散系统的时域特性。

（3）分析、观察及检验系统的稳定性。

2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。

已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。

在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。

也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。

系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。

重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。

系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。

或者系统的单位脉冲响应满足绝对可和的条件。

第十章上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验.上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力.本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验.实验一系统响应及系统稳定性。

实验二时域采样与频域采样。

实验三用FFT对信号作频谱分析。

实验四IIR数字滤波器设计及软件实现。

实验五FIR数字滤波器设计与软件实现实验六应用实验—-数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。

建议自学的读者在学习完第一章后作实验一;在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在.学习完第六章进行；实验五在学习完第七章后进行。

实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行.function tstem(xn，yn)%时域序列绘图函数％ xn：信号数据序列，yn：绘图信号的纵坐标名称（字符串)n=0:length（xn)—1；stem（n,xn,’.’);box onxlabel(’n'）；ylabel(yn)；axis（[0,n(end），min(xn）,1.2*max(xn）］）10.1 实验一: 系统响应及系统稳定性1。

实验目的（1）掌握求系统响应的方法。

(2）掌握时域离散系统的时域特性.(3）分析、观察及检验系统的稳定性.2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。

已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应,本实验仅在时域求解。

在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB语言的工具箱函数filter函数.也可以用MATLAB语言的工具箱函数conv函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。

数字信号处理实验答案第十章上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。

上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。

本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。

实验一系统响应及系统稳定性。

实验二时域采样与频域采样。

实验三用FFT对信号作频谱分析。

实验四IIR数字滤波器设计及软件实现。

实验五FIR数字滤波器设计与软件实现实验六应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。

建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。

学习完第六章进行；实验五在学习完第七章后进行。

实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。

实验一: 系统响应及系统稳定性1.实验目的（1）掌握求系统响应的方法。

（2）掌握时域离散系统的时域特性。

（3）分析、观察及检验系统的稳定性。

2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。

已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。

在计算机上适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱函数filter函数。

也可以用MATLAB语言的工具箱函数conv函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。

系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。

重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。

系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。

或者系统的单位脉冲响应满足绝对可和的条件。

系统的稳定性由其差分方程的系数决定。

数字信号处理上机答案（含程序及图片）第三版高西全丁玉美著数字信号处理实验一内容一a=0.8;ys=0;A=[1,-0.9];B=[0.05,0.05];xn=[1,zeros(1,50)];x1n=[1 1 1 1 1 1 1 1 zeros(1,50)];x2n=ones(1,128);xi=filtic(B,A,ys);hn=filter(B,A,xn,xi)n=0:length(hn)-1;subplot(2,2,1);stem(n,yn,'.')title('(a) 系统单位脉冲响应h(n)');xlabel('n');ylabel(hn);y1n=filter(B,A,x1n,xi);n=0:length(y1n)-1;subplot(2,2,2);y='y1(n)'; stem(n,y1n,'.')title('(b) 系统对R8(n)的响应y1(n)');xlabel('n');ylabel(yn);y2n=filter(B,A,x2n,xi);n=0:length(y2n)-1;subplot(2,2,4);y='y2(n)'; stem(n,y2n,'.')title('(c) 系统对u(n)的响应y2(n)');xlabel('n');ylabel(yn);20400.020.040.060.080.1nh (n )(a) 系统单位脉冲响应h(n)020400.20.40.6ny 1(n )(b) 系统对R8(n)的响应y1(n)501000.20.40.60.81ny 2(n )(c) 系统对u(n)的响应y2(n)内容二x1n=[1 1 1 1 1 1 1 1 ];h1n=[ones(1,10) zeros(1,10)]; h2n=[1 2.5 2.5 1 zeros(1,10)]; y21n=conv(h1n,x1n); y22n=conv(h2n,x1n); M1=length(y21n)-1; M2=length(y22n)-1; n1=0:1:M1; n2=0:1:M2;n11=0:length(h1n)-1; n22=0:length(h2n)-1;subplot(2,2,1); tstem(n11,h1n); title('(d) 系统单位脉冲响应h1(n)'); xlabel('n');ylabel(h1(n));subplot(2,2,2); stem(n1,y21n,'fill'); title('(e) h1(n)与R8(n)的卷积y21(n)'); xlabel('n');ylabel(y21(n));subplot(2,2,3); tstem(n22,h2n); title('(f) 系统单位脉冲响应h2(n)'); xlabel('n');ylabel(h2(n));subplot(2,2,4); stem(n1,y22n,'fill'); title('(g) h2(n)与R8(n)的卷积y22(n)'); xlabel('n');ylabel(y22(n));5101500.51nh 1(n )(d) 系统单位脉冲响应h1(n)010202468ny 21(n )(e) h1(n)与R8(n)的卷积y21(n)510123nh 2(n )(f) 系统单位脉冲响应h2(n)510152002468ny 22(n )(g) h2(n)与R8(n)的卷积y22(n)内容三谐振器对u(n)的响应a=0.8;ys=0;xn=[1,zeros(1,250)];B=[1/100.49,-1/100.49];A=[1,-1.8237,0.9801]; xi=filtic(B,A,ys); yn=filter(B,A,xn,xi) n=0:length(yn)-1;subplot(1,1,1);stem(n,yn,'.')谐振器对正弦信号的响应a=0.8;ys=0;xsin=sin(0.014*n)+sin(0.4*n);B=[1/100.49,-1/100.49];A=[1,-1.8237,0.9801]; xi=filtic(B,A,ys); yn=filter(B,A,xsin,xi) n=0:length(yn)-1;subplot(1,1,1);stem(n,yn,'.')50100150200250-0.01-0.008-0.006-0.004-0.00200.0020.0040.0060.0080.0150100150200250-0.5-0.4-0.3-0.2-0.100.10.20.30.40.5数字信号处理实验三实验（1）x1n=[ones(1,4)]; X1k8=fft(x1n,8); X1k16=fft(x1n,16); N=8;f=2/N*(0:N-1); figure(1);subplot(1,2,1);stem(f,abs(X1k8),'.'); title('(la) 8µãDFT[x_1(n)]');xlabel('\omega/\pi¡¯);ylabel(¡®|(e^j^\omega)|'); N=16;f=2/N*(0:N-1);subplot(1,2,2);stem(f,abs(X1k16),'.'); title('(la) 16µãDFT[x_1(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|');实验（1-2，1-3）M=8;xa=1:(M/2);xb=(M/2):-1:1;x2n=[xa,xb];x3n=[xb,xa];X2k8=fft(x2n,8);X2k16=fft(x2n,16);X3k8=fft(x3n,8);X3k16=fft(x3n,16);figure(2);N=8;f=2/N*(0:N-1);subplot(2,2,1);stem(f,abs(X2k8),'.');title('(2a) 8µãDFT[x_2(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|'); subplot(2,2,3);stem(f,abs(X3k8),'.');title('(3a) 8µãDFT[x_3(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|'); N=16;f=2/N*(0:N-1);subplot(2,2,2);stem(f,abs(X2k16),'.');title('(2a) 16µãDFT[x_2(n)]');xlabel('');ylabel('');subplot(2,2,4);stem(f,abs(X3k16),'.');title('(3a) 16µãDFT[x_3(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|');实验（2-1，2-2）N=8;n=0:N-1;x4n=cos(pi*n/4);x5n=cos(pi*n/4)+cos(pi*n/8);X4k8=fft(x4n,8);X4k16=fft(x4n,16);X5k8=fft(x5n,8);X5k16=fft(x5n,16);figure(3);N=8;f=2/N*(0:N-1);subplot(2,2,1);stem(f,abs(X4k8),'.');title('(4a) 8µãDFT[x_4(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|'); subplot(2,2,3);stem(f,abs(X5k8),'.');title('(5a) 8µãDFT[x_5(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|'); N=16;f=2/N*(0:N-1);subplot(2,2,2);stem(f,abs(X4k16),'.');title('(4a) 16µãDFT[x_4(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|'); subplot(2,2,4);stem(f,abs(X5k16),'.');title('(5a) 16µãDFT[x_5(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|');实验（3）Fs=64;T=1/Fs;N=16;n=0:N-1;nT=n*T;x8n=cos(8*pi*nT)+cos(16*pi*nT)+cos(20*pi*nT); X8k16=fft(x8n,16);N=16;f=2/N*(0:N-1);figure(4);subplot(2,2,1);stem(f,abs(X8k16),'.');title('(8a) 16µãDFT[x_8(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|'); N=32;n=0:N-1;nT=n*T;x8n=cos(8*pi*nT)+cos(16*pi*nT)+cos(20*pi*nT); X8k32=fft(x8n,32);N=32;f=2/N*(0:N-1);figure(4);subplot(2,2,2);stem(f,abs(X8k32),'.');title('(8a) 32µãDFT[x_8(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|'); N=64;n=0:N-1;nT=n*T;x8n=cos(8*pi*nT)+cos(16*pi*nT)+cos(20*pi*nT); X8k64=fft(x8n,64);N=64;f=2/N*(0:N-1);figure(4);subplot(2,2,3);stem(f,abs(X8k64),'.');title('(8a) 64µãDFT[x_8(n)]');xlabel('\omega/\pi');ylabel('|(e^j^\omega)|');数字信号处理实验四内容一function st=mstgN=800Fs=10000;T=1/Fs;Tp=N*T; t=0:T:(N-1)*T;k=0:N-1;f=k/Tp; fc1=Fs/10; fm1=fc1/10; fc2=Fs/20; fm2=fc2/10; fc3=Fs/40; fm3=fc3/10;xt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); st=xt1+xt2+xt3; fxt=fft(st,N); subplot(3,1,1)plot(t,st);grid;xlabel('t/s');ylabel('s(t)');axis([0,Tp/8,min(st),max(st)]);title('(a) s(t)的波形') subplot(3,1,2)stem(f,abs(fxt)/max(abs(fxt)),'.');grid;title('(b) s(t)的频谱') axis([0,Fs/5,0,1.2]);xlabel('f/Hz');ylabel('幅度')0.0010.0020.0030.0040.0050.0060.0070.0080.0090.01-10123t/ss (t )(a) s(t)的波形20040060080010001200140016001800200000.51(b) s(t)的频谱f/Hz幅度内容二Fs=10000;T=1/Fs;st=mstg;%低通滤波器设计与实现fp=280;fs=450;wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60;[N,wp]=ellipord(wp,ws,rp,rs);[B,A]=ellip(N,rp,rs,wp);y1t=filter(B,A,st);figure(2);subplot(3,1,1);[H,w]=freqz(B,A,1000);m=abs(H);plot(w/pi,20*log(m/max(m)));grid on;title('低通滤波损耗函数曲线'); xlabel('w/pi ');ylabel('幅度'); axis([0,1,0,1.2*max(H)])yt='y1(t)'; subplot(3,1,2); plot(t,y1t);title('低通滤波后的波形');xlabel('t/s');ylabel(y1(t));%带通滤波器设计与实现fpl=440;fpu=560;fsl=275;fsu=900;wp=[2*fpl/Fs,2*fpu/Fs];ws=[2*fsl/Fs,2*fsu/Fs];rp=0.1;rs=60;[N,wp]=ellipord(wp,ws,rp,rs);[B,A]=ellip(N,rp,rs,wp);y2t=filter(B,A,st);figure(3);subplot(3,1,1);[H,w]=freqz(B,A,1000);m=abs(H);plot(w/pi,20*log(m/max(m)));grid on;title('带通滤波损耗函数曲线'); xlabel('w/pi ');ylabel('幅度'); axis([0,1,0,1.2*max(H)])yt='y2(t)'; subplot(3,1,2); plot(t,y2t);title('带通滤波后的波形');xlabel('t/s');ylabel(y2(t));%高通滤波器设计与实现fp=890;fs=600;wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60;[N,wp]=ellipord(wp,ws,rp,rs);[B,A]=ellip(N,rp,rs,wp,'high');y3t=filter(B,A,st);figure(4);subplot(3,1,1);[H,w]=freqz(B,A,1000);m=abs(H);plot(w/pi,20*log(m/max(m)));grid on;title('高通滤波损耗函数曲线'); xlabel('w/pi ');ylabel('幅度'); axis([0,1,0,1.2*max(H)])yt='y3(t)'; subplot(3,1,2); plot(t,y3t);title('高通滤波后的波形');xlabel('t/s');ylabel(y3(t));低通滤波器损耗函数及其分离出的调幅信号y1(t)带通滤波器损耗函数及其分离出的调幅信号y2(t)高通滤波器损耗函数及其分离出的调幅信号y3(t)数字信号处理实验五1、function xt=xtg(N)Fs=1000;T=1/Fs;Tp=N*T;t=0:T:(N-1)*T;fc=Fs/10;f0=fc/10;mt=cos(2*pi*f0*t);ct=cos(2*pi*fc*t);xt=mt.*ct;nt=2*rand(1,N)-1;fp=150; fs=200;Rp=0.1;As=60;fb=[fp,fs];m=[0,1];dev=[10^(-As/20),(10^(Rp/20)-1)/(10^(Rp/20)+1)];[n,fo,mo,W]=remezord(fb,m,dev,Fs);hn=remez(n,fo,mo,W);yt=filter(hn,1,10*nt);xt=xt+yt;fst=fft(xt,N);k=0:N-1;f=k/Tp;subplot(3,1,1);plot(t,xt);grid;xlabel('t/s');ylabel('x(t)');axis([0,Tp/5,min(xt),max(xt)]);title('(a) 信号加噪声波形')subplot(3,1,2);plot(f,abs(fst)/max(abs(fst)));grid;title('(b) 信号加噪声的频谱')axis([0,Fs/2,0,1.2]);xlabel('f/Hz');ylabel('幅度')2、xt=xtg;N=1000;Fs=1000;T=1/Fs;Tp=N*T;k=0:N-1;f=k/Tp;t=0:T:(N-1)*T;fp=120;fs=150;Rp=0.1;As=60;Fs=1000;wc=(fp+fs)/Fs;B=2*pi*(fs-fp)/Fs;M=ceil(11*pi/B);hn=fir1(M-1,wc,blackman(M));Hw=abs(fft(hn,N));ywt=fftfilt(hn,xt,N);figure;subplot(2,1,1);plot(f,20*log10(Hw)/max(Hw));grid onxlabel('f/Hz');ylabel('幅度（dB ）');title('(a)低通滤波器的幅频特性')axis([0,500,-160,5]);subplot(2,1,2);plot(t,ywt);grid onxlabel('t/s');ylabel('y_1(t)');title('(b)滤除噪声后的信号波形')050100150200250300350400450500-150-100-500f/Hz幅度（d B ）(a)低通滤波器的幅频特性00.10.20.30.40.50.60.70.80.91-1-0.50.51t/s y 1(t )(b)滤除噪声后的信号波形。

数字信处理上机实验答案GE GROUP system office room 【GEIHUA16H-GEIHUA GEIHUA8Q8-第十章上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。

上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。

本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。

实验一系统响应及系统稳定性。

实验二时域采样与频域采样。

实验三用FFT对信号作频谱分析。

实验四 IIR数字滤波器设计及软件实现。

实验五 FIR数字滤波器设计与软件实现实验六应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。

建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。

学习完第六章进行；实验五在学习完第七章后进行。

实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。

functiontstem(xn,yn)%时域序列绘图函数%xn:信号数据序列，yn:绘图信号的纵坐标名称（字符串）n=0:length(xn)-1;stem(n,xn,'.');boxonxlabel('n');ylabel(yn);axis([0,n(end),min(xn),1.2* max(xn)])10.1 实验一: 系统响应及系统稳定性1.实验目的（1）掌握求系统响应的方法。

（2）掌握时域离散系统的时域特性。

（3）分析、观察及检验系统的稳定性。

2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。

已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。

第十章上机实验数字信号处理是一门理论和实际密切结合的课程.为深入掌握课程容.最好在学习理论的同时.做习题和上机实验。

上机实验不仅可以帮助读者深入的理解和消化基本理论.而且能锻炼初学者的独立解决问题的能力。

本章在第二版的基础上编写了六个实验.前五个实验属基础理论实验.第六个属应用综合实验。

实验一系统响应及系统稳定性。

实验二时域采样与频域采样。

实验三用FFT对信号作频谱分析。

实验四 IIR数字滤波器设计及软件实现。

实验五 FIR数字滤波器设计与软件实现实验六应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度.安排学生上机进行实验。

建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。

学习完第六章进行；实验五在学习完第七章后进行。

实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验.在学习完本课程后再进行。

function tstem(xn,yn)%时域序列绘图函数% xn:信号数据序列.yn:绘图信号的纵坐标名称（字符串）n=0:length(xn)-1;stem(n,xn,'.');box onxlabel('n');ylabel(yn);axis([0,n(end),min(xn),1.2*max(xn)])10.1 实验一: 系统响应及系统稳定性1.实验目的（1）掌握求系统响应的方法。

（2）掌握时域离散系统的时域特性。

（3）分析、观察及检验系统的稳定性。

2.实验原理与方法在时域中.描写系统特性的方法是差分方程和单位脉冲响应.在频域可以用系统函数描述系统特性。

已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应.本实验仅在时域求解。

在计算机上适合用递推法求差分方程的解.最简单的方法是采用MATLAB语言的工具箱函数filter函数。

也可以用MATLAB语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积.求出系统的响应。