离散时间系统的时域特性分析(数字信号处理实验1答案)
数字信号处理(第三版)第1章习题答案
第 1 章 时域离散信号和时域离散系统
1.1.1
（1） 信号： 模拟信号、 时域离散信号、 数字信号三 者之间的区别； 常用的时域离散信号； 如何判断信号是周期 性的， 其周期如何计算等。
（2） 系统： 什么是系统的线性、 时不变性以及因果 性、 稳定性； 线性、 时不变系统输入和输出之 间的关系； 求解线性卷积的图解法（列表法）、 解析法， 以及用MATLAB工具箱函数求解； 线性常系数差分方程的递
x(n－n0)=x(n)*δ(n－n0)
（3）
Xˆ n ( j )
Байду номын сангаас
1 T
X a ( j
k
jks )
这是关于采样定理的重要公式， 根据该公式要求对
信号的采样频率要大于等于该信号的最高频率的两倍以上，
才能得到不失真的采样信号。
xa
(t
)
n
xa
(nt
)
sin[π(t nT ) / T π(t nT ) / T
第 1 章 时域离散信号和时域离散系统
第1章 时域离散信号和时域离散系统
1.1 学习要点与重要公式 1.2 解线性卷积的方法 1.3 例题 1.4 习题与上机题解答
第 1 章 时域离散信号和时域离散系统
1.1 学习要点与重要公式
本章内容是全书的基础。 学生从学习模拟信号分析与处 理到学习数字信号处理， 要建立许多新的概念。 数字信号 和数字系统与原来的模拟信号和模拟系统不同， 尤其是处理 方法上有本质的区别。 模拟系统用许多模拟器件实现， 数 字系统则通过运算方法实现。 如果读者对本章关于时域离散 信号与系统的若干基本概念不清楚， 则学到数字滤波器时， 会感到“数字信号处理”这门课不好掌握， 总觉得学习的不 踏实。 因此学好本章是极其重要的。
数字信号处理 实验作业：离散LSI系统的时域分析
实验2 离散LSI 系统的时域分析一、.实验目的：1、加深对离散系统的差分方程、单位脉冲响应、单位阶跃响应和卷积分析方法的理解。
2、初步了解用MA TLAB 语言进行离散时间系统时域分析的基本方法。
3、掌握求解离散时间系统的单位脉冲响应、单位阶跃响应、线性卷积以及差分方程的程序的编写方法，了解常用子函数的调用格式。
二、实验原理：1、离散LSI 系统的响应与激励由离散时间系统的时域分析方法可知，一个离散LSI 系统的响应与激励可以用如下框图表示：其输入、输出关系可用以下差分方程描述：[][]NMkk k k ay n k b x n m ==-=-∑∑2、用函数impz 和dstep 求解离散系统的单位脉冲响应和单位阶跃响应。
例2-1 已知描述某因果系统的差分方程为6y(n)+2y(n-2)=x(n)+3x(n-1)+3x(n-2)+x(n-3) 满足初始条件y(-1)=0，x(-1)=0，求系统的单位脉冲响应和单位阶跃响应。
解： 将y(n)项的系数a 0进行归一化，得到y(n)+1/3y(n-2)=1/6x(n)+1/2x(n-1)+1/2x(n-2)+1/6x(n-3)分析上式可知，这是一个3阶系统，列出其b k 和a k 系数： a 0=1， a ,1=0， a ,2=1/3， a ,3=0 b 0=1/6，b ,1=1/2， b ,2=1/2， b ,3=1/6程序清单如下： a=[1,0,1/3,0]; b=[1/6,1/2,1/2,1/6]; N=32; n=0:N-1; hn=impz(b,a,n); gn=dstep(b,a,n);subplot(1,2,1);stem(n,hn,'k');课程名称 数字信号处理 实验成绩 指导教师 ***实 验 报 告院系 班级学号 姓名 日期title('系统的单位序列响应'); ylabel('h(n)');xlabel('n');axis([0,N,1.1*min(hn),1.1*max(hn)]); subplot(1,2,2);stem(n,gn,'k'); title('系统的单位阶跃响应'); ylabel('g(n)');xlabel('n');axis([0,N,1.1*min(gn),1.1*max(gn)]); 程序运行结果如图2-1所示：102030系统的单位序列响应h (n )n1020300.20.30.40.50.60.70.80.911.11.2系统的单位阶跃响应g (n )n图2-13、用函数filtic 和filter 求解离散系统的单位序列响应和单位阶跃响应。
“数字信号处理”实验报告二
实验报告课程名称：数字信号处理实验任课教师：杨鉴实验名称：离散时间系统的时域分析年级、专业：2015级通信工程学号：**********姓名：***日期：2017 年10 月9 日云南大学信息学院一、实验1.通过MATLAB仿真一些简单的离散时间系统，并研究他们的时域特性。
2.掌握卷积在MATLAB的算法并理解滤波的概念。
二、实验内容1. 假定另一个系统为y[n]=x[n]x[n-1],修改程序P2.3，计算这个系统的输出序列y1[n],y2[n]和y[n]。
比较y[n]和yt[n]。
这两个序列是否相等？该系统是线性系统吗？2. 考虑另一个系统：y[n]=nx[n]+x[n-1],修改程序P2.4，以仿真上面的系统并确定该系统是否为时不变系统。
3.修改程序P2.7，计算长度为15的序列h[n]和长度为10的序列x[n]的卷积，重做问题Q2.28。
h[n]和x[n]的样本值你自己给定。
4.修改程序P2.9，将输入序列改变成扫频正弦序列（长度为301、最低频率为0、最高频率为0.5）。
那个滤波器能更好的抑制输入信号x[n]的高频分量？三、主要算法与程序Q2.11：clf;n = 0:40;a = 2;b = -3;f1=0.1;f2=0.4;x11=[0 cos(2*pi*f1*n) 0];x12=[0 0 cos(2*pi*f1*n)];x21=[0 cos(2*pi*f2*n) 0];x22=[0 0 cos(2*pi*f2*n)];x = a*x11 + b*x21;y1 = x11.*x12;y2 = x21.*x22;xd = a*x12+b*x22;y = x.*xd;yt = a*y1 + b*y2;d = y - yt; % Compute the difference output d[n]% Plot the outputs and the difference signalsubplot(3,1,1)stem([0 n 0],y);ylabel('Amplitude');title('Output Due to Weighted Input: a \cdot x_{1}[n] + b \cdot x_{2}[n]');subplot(3,1,2)stem([0 n 0],yt);ylabel('Amplitude');title('Weighted Output: a \cdot y_{1}[n] + b \cdot y_{2}[n]'); subplot(3,1,3)stem([0 n 0],d);xlabel('Time index n');ylabel('Amplitude');title('Difference Signal');Q2.17：clf;n = 0:40; D = 10;a = 3.0;b = -2;x = a*cos(2*pi*0.1*n) + b*cos(2*pi*0.4*n);xd = [zeros(1,D) x];nd=0:length(xd)-1;y=(n.*x)+[0 x(1:40)];yd=(nd.*xd)+[0 xd(1:length(xd)-1)];d = y - yd(1+D:41+D);subplot(3,1,1)stem(n,y);ylabel('振幅');title('输出 y[n]'); grid;subplot(3,1,2)stem(n,yd(1:41));ylabel('振幅');title('由于延时输入 x[n-10]的输出'); grid;subplot(3,1,3)stem(n,d);xlabel('时间序号 n'); ylabel('振幅');title('差值信号');grid;Q2.29：clf;h = [3 2 1 -2 1 0 -4 0 3 1 5 4 0 3 5]; % impulse responsex = [1 -2 3 -4 3 2 1 5 6 1]; % input sequencey = conv(h,x);n = 0:23;subplot(2,1,1);stem(n,y);xlabel('时间序号n');ylabel('振幅');title('用卷积得到的输出'); grid;x1 = [x zeros(1,14)];y1 = filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel('时间序号 n'); ylabel('振幅');title('由滤波生成的输出'); grid;Q2.35：f=w/2pi=(2a*n+b)/2pi=[0,0.5],所以b=0，a*n<=0.5*pi,当n=300时，a取pi/600。
数字信号处理课后习题答案(全)1-7章
x(n)=－δ(n+2)+δ(n－1)+2δ(n－3)
h(n)=2δ(n)+δ(n－1)+ δ(n－2)
由于
x(n)*δ(n)=x(n)
1
x(n)*Aδ(n－k)=Ax(n－k)
2
故
第 1 章 时域离散信号和时域离散系统
y(n)=x(n)*h(n)
=x(n)*［2δ(n)+δ(n－1)+ δ(n－2) 1 2
（5） 系统是因果系统， 因为系统的输出不取决于x(n)的未来值。 如果
|x(n)|≤M, 则|y(n)|=|ex(n)|≤e|x(n)|≤eM，
7． 设线性时不变系统的单位脉冲响应h(n)和输入序列x(n)如题7图所示，
要求画出y(n)输出的波形。
解： 解法（一）采用列表法。
y(n)=x(n)*h(n)=
0≤m≤3
－4≤m≤n
非零区间如下:
第 1 章 时域离散信号和时域离散系统
根据非零区间， 将n分成四种情况求解： ① n<0时， y(n)=0
② 0≤n≤3时， y(n)= ③ 4≤n≤7时， y(n)= ④ n>7时， y(n)=0
1=n+1
n
1=8－m n0
3
mn4
第 1 章 时域离散信号和时域离散系统
第 1 章 时域离散信号和时域离散系统
(3) 这是一个延时器， 延时器是线性非时变系统， 下面证明。 令输入为
输出为
x(n－n1)
y′(n)=x(n－n1－n0) y(n－n1)=x(n－n1－n0)=y′(n) 故延时器是非时变系统。 由于
T［ax1(n)+bx2(n)］=ax1(n－n0)+bx2(n－n0) =aT［x1(n)］+bT［x2(n)］
数字信号处理实验报告一二
数字信号处理课程实验报告实验一 离散时间信号和系统响应一. 实验目的1. 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解2. 掌握时域离散系统的时域特性3. 利用卷积方法观察分析系统的时域特性4. 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号及系统响应进行频域分析二、实验原理1. 采样是连续信号数字化处理的第一个关键环节。
对采样过程的研究不仅可以了解采样前后信号时域和频域特性的变化以及信号信息不丢失的条件，而且可以加深对离散傅里叶变换、Z 变换和序列傅里叶变换之间关系式的理解。
对连续信号()a x t 以T 为采样间隔进行时域等间隔理想采样，形成采样信号： 式中()p t 为周期冲激脉冲，()a x t 为()a x t 的理想采样。
()a x t 的傅里叶变换为()a X j Ω：上式表明将连续信号()a x t 采样后其频谱将变为周期的，周期为Ωs=2π/T 。
也即采样信号的频谱()a X j Ω是原连续信号xa(t)的频谱Xa(jΩ)在频率轴上以Ωs 为周期，周期延拓而成的。
因此，若对连续信号()a x t 进行采样，要保证采样频率fs ≥2fm ，fm 为信号的最高频率，才可能由采样信号无失真地恢复出原模拟信号ˆ()()()a a xt x t p t =1()()*()21()n a a a s X j X j P j X j jn T π∞=-∞Ω=ΩΩ=Ω-Ω∑()()n P t t nT δ∞=-∞=-∑计算机实现时，利用计算机计算上式并不方便，因此我们利用采样序列的傅里叶变换来实现，即而()()j j n n X e x n e ωω∞-=-∞=∑为采样序列的傅里叶变换2. 时域中，描述系统特性的方法是差分方程和单位脉冲响应，频域中可用系统函数描述系统特性。
已知输入信号，可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应。
数字信号处理实验报告
一、实验目的1. 理解数字信号处理的基本概念和原理。
2. 掌握离散时间信号的基本运算和变换方法。
3. 熟悉数字滤波器的设计和实现。
4. 培养实验操作能力和数据分析能力。
二、实验原理数字信号处理（Digital Signal Processing，DSP）是利用计算机对信号进行采样、量化、处理和分析的一种技术。
本实验主要涉及以下内容：1. 离散时间信号：离散时间信号是指时间上离散的信号，通常用序列表示。
2. 离散时间系统的时域分析：分析离散时间系统的时域特性，如稳定性、因果性、线性等。
3. 离散时间信号的变换：包括离散时间傅里叶变换（DTFT）、离散傅里叶变换（DFT）和快速傅里叶变换（FFT）等。
4. 数字滤波器：设计、实现和分析数字滤波器，如低通、高通、带通、带阻滤波器等。
三、实验内容1. 离散时间信号的时域运算（1）实验目的：掌握离散时间信号的时域运算方法。
（2）实验步骤：a. 使用MATLAB生成两个离散时间信号；b. 进行时域运算，如加、减、乘、除等；c. 绘制运算结果的时域波形图。
2. 离散时间信号的变换（1）实验目的：掌握离散时间信号的变换方法。
（2）实验步骤：a. 使用MATLAB生成一个离散时间信号；b. 进行DTFT、DFT和FFT变换；c. 绘制变换结果的频域波形图。
3. 数字滤波器的设计和实现（1）实验目的：掌握数字滤波器的设计和实现方法。
（2）实验步骤：a. 设计一个低通滤波器，如巴特沃斯滤波器、切比雪夫滤波器等；b. 使用MATLAB实现滤波器；c. 使用MATLAB对滤波器进行时域和频域分析。
4. 数字滤波器的应用（1）实验目的：掌握数字滤波器的应用。
（2）实验步骤：a. 采集一段语音信号；b. 使用数字滤波器对语音信号进行降噪处理；c. 比较降噪前后的语音信号，分析滤波器的效果。
四、实验结果与分析1. 离散时间信号的时域运算实验结果显示，通过MATLAB可以方便地进行离散时间信号的时域运算，并绘制出运算结果的时域波形图。
数字信号处理实验离散时间 LTI 系统的时域分析与 Z 域分析
实验一离散时间LTI系统的时域分析与Z域分析一、实验目的1、掌握用MATLAB求解离散时间系统的零状态响应、单位脉冲响应和单位阶跃响应；2、掌握离散时间系统系统函数零极点的计算方法和零极点图的绘制方法，并能根据零极点图分析系统的稳定性。
二、实验原理1、离散时间系统的时域分析（1）离散时间系统的零状态响应离散时间LTI系统可用线性常系数差分方程来描述，即MATLAB中函数filter可对式（1-1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。
函数filter的语句格式为：y=filter(b,a,x)其中，x为输入的离散序列；y为输出的离散序列；y的长度与x的长度一样；b与a分别为差分方程右端与左端的系数向量。
（2）离散时间系统的单位脉冲响应系统的单位脉冲响应定义为系统在 (n)激励下系统的零状态响应，用h(n)表示。
MATLAB求解单位脉冲响有两种方法：一种是利用函数filter；另一种是利用函数impz。
impz函数的常用语句格式为impz(b,a,n)，其中b和a的定义见filter，n表示脉冲响应输出的序列个数。
（3）离散时间系统的单位阶跃响应系统的单位阶跃响应定义为系统在ε(n)激励下系统的零状态响应。
MATLAB求解单位脉冲响应有两种方法：一种是利用函数filter，另一种是利用函数stepz。
stepz函数的常用语句格式为stepz(b,a,N)其中，b和a的定义见filter，N表示脉冲响应输出的序列个数。
2、离散时间系统的Z域分析(1)系统函数的零极点分析离散时间系统的系统函数定义为系统零状态响应的z变换与激励的z变换之比，即如果系统函数H(z)的有理函数表示式为那么，在MATLAB中系统函数的零极点就可通过函数roots得到，也可借助函数tf2zp得到。
roots的语法格式为：Z=roots(b)%计算零点b=[b1b2…bmbm+1]P=roots(a)%计算极点a=[a1a2…anan+1]tf2zp的语句格式为[Z,P,K]=tf2zp(b,a)其中，b与a分别表示H(z)的分子与分母多项式的系数向量。
数字信号处理课后习题答案
数字信号处理课后习题答案数字信号处理课后习题答案数字信号处理是一门重要的学科，它研究如何对数字信号进行处理和分析。
在学习过程中，我们经常会遇到一些习题，通过解答这些习题可以帮助我们更好地理解和掌握数字信号处理的知识。
本文将为大家提供一些数字信号处理课后习题的答案，希望对大家的学习有所帮助。
一、离散时间信号和系统1. 什么是离散时间信号？答：离散时间信号是在离散时间点上取值的信号，它可以用数学上的序列表示。
2. 什么是离散时间系统？答：离散时间系统是对离散时间信号进行处理的系统，它可以用差分方程或差分方程组来描述。
3. 离散时间信号和连续时间信号有何区别？答：离散时间信号是在离散时间点上取值的信号，而连续时间信号是在连续时间上取值的信号。
二、离散时间信号的表示和运算1. 如何表示离散时间信号？答：离散时间信号可以用数学上的序列表示，例如x(n)表示离散时间信号x在时间点n上的取值。
2. 离散时间信号的运算有哪些？答：离散时间信号的运算包括加法、减法、乘法和卷积等。
3. 什么是离散时间信号的卷积？答：离散时间信号的卷积是指两个离散时间信号之间的一种数学运算，它可以表示两个信号之间的线性叠加关系。
三、离散时间系统的性质和稳定性1. 离散时间系统有哪些常见的性质？答：离散时间系统常见的性质包括线性性、时不变性、因果性和稳定性等。
2. 什么是离散时间系统的稳定性？答：离散时间系统的稳定性是指当输入信号有界时，输出信号也有界。
3. 如何判断离散时间系统的稳定性？答：可以通过判断系统的冲激响应的绝对可和性来判断离散时间系统的稳定性。
四、离散傅里叶变换1. 什么是离散傅里叶变换（DFT）？答：离散傅里叶变换是将离散时间信号转换为离散频率信号的一种数学变换。
2. 离散傅里叶变换有何作用？答：离散傅里叶变换可以将时域的信号转换为频域的信号，从而方便对信号的频谱进行分析。
3. 如何计算离散傅里叶变换？答：可以通过对离散时间信号进行离散傅里叶变换公式的计算来得到离散傅里叶变换的结果。
数字信号处理实验1--5含代码
数字信号处理实验1--5含代码实验一离散时间信号的时域分析 1. 在MATLAB中利用逻辑关系式n,,0来实现序列，显示范围。
(产生如下，，,n,nn,n,n012图所示的单位脉冲信号的函数为impseq(n0,n1,n2)，程序如示例所示),3,n,10并利用impseq函数实现序列:; ，，，，，，yn,2,n,3，,n,6，，xn1nnnn120源代码:impseq.mfunction y=impseq(n0,n1,n2)n=[n1:n2]y=[(n-n0)==0]exp01-1.mfunction impseq(n0,n1,n2)n=-3:1:10y=2*impseq(3,-3,10)+impseq(6,-3,10);stem(n,y)n,,0，，2. 在MATLAB中利用逻辑关系式来实现序列，显示范围。
(自己编写un,nn,n,n012产生单位阶跃信号的函数，函数命名为stepseq(n0,n1,n2)) 并利用编写的stepseq函数实现序列: ，，，，，，yn,un，2，un,2,5,n,10源代码:stepseq.mfunction y=stepseq(n0,n1,n2)n=n1:1:n2y=[(n-n0)>=0]exp01-2.mfunction stepseq(n0,n1,n2)n=-5:1:20y=stepseq(-2,-5,20)+stepseq(2,-5,20)stem(n,y)3. 在MATLAB中利用数组运算符“.^”来实现一个实指数序列。
如: n ，，，，xn,0.30,n,15源代码:n=0:1:15;x=0.3.^nstem(n,x)4. 在MATLAB中调用函数sin或cos产生正余弦序列，如:π,, ，，，，xn,3sin0.4πn，，5cos0.3πn0,n,20,,5,,源代码:n=0:1:20x=11*sin(0.3*pi*n+pi/5)+5*cos(0.3*pi*n)stem(n,x)思考题:1.在MATLAB环境下产生单位脉冲序列和单位阶跃序列各有几种方法,如何使用,2.在MATLAB环境下进行序列的相乘运算时应注意什么问题,实验二离散时间系统的时域分析1. 在MATLAB中利用内部函数conv来计算两个有限长序列的卷积。
离散时间系统时域分析及稳定性实验
电子系电子信息工程实验报告 课程名称：《基于MA TLAB 的信号与系统及数字信号处理仿真实验》 实验项目名称：离散时间系统时域分析及稳定性实验实验时间：2012-6-13实验地点：信息学院四层机房班级：电子信息工程姓名：陆阿楠学号：2010117119一、实验目的1、掌握求系统响应的方法2、掌握时域离散系统的时域特性3、分析、观察及检验系统的稳定性二、实验原理在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。
已经输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应本实验仅在时域求解。
在计算机上适合用递推方求解差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。
也可以用MA TLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。
系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。
重点分析实验系统的稳定性，包括观察系统的暂态响应和稳态响应。
系统的稳定性是指对任意有节的输入信号，系统都能得到有界的系统响应。
或者系统的单位脉冲响应满足绝对可和的条件。
系统的稳定性由其差分方程的系数决定。
实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界的输出，或者检查系统的单位脉冲响应是否满足绝对可和的条件。
可行的方法是在系统的输入端加入单位阶跃序列，如果系统输出趋近一个常数（包括零），就可以断定系统是稳定的。
系统的稳态输出是指当n 趋向于无穷大时，系统的输出。
如果系统稳定信号加入系统后，系统输出的开始一段称为暂态响应，随n 的增大，幅度趋于稳定，达到稳态输出。
三、实验内容（1）编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter 函数或conv 函数求解系统输出响应的主程序。
程序中要有绘制信号波形的功能。
程序：A=input('输入：');B=input('单位脉冲响应序列');C=conv(A,B);stem(C);xlabel('n');ylabel('系统响应');输入：[1,3,5,7,9]单位脉冲响应序列：[2,4,6,8,10]图形如右：（2）给定一个低通滤波器的差分方程为)1(9.0)1(05.0)(05.0)(-+-+=n y n x n x n y 输入信号)()(81n R n x =)()(2n u n x =a) 分别求出系统对)()(81n R n x =和)()(2n u n x =的响应序列，并画出其波形。
数字信号处理实验MATLAB上机DOC
班级： 学号： 姓名： 日期： 实验一：离散时间信号的分析一、实验目的利用DFT 卷积实现系统的时域分析二、实验原理在离散时间、连续频率的傅里叶变换中，由于卷积性质知道，对系统输出的计算可以通过求x[n]和h[n]的DTFT ，将得到的X(e jw )和H(e jw )相乘就可以得到Y(e jw )，进而再通过反变换得到y[n]。
这就避免了在时域进行繁琐的卷积求解。
三、实验步骤（包括代码和波形）1-2（2）x[k]=g[k]=k+1,0<=k<=3;x[k]=g[k]=0,其他 编码如下：ak=1:4 gk=1:4Z=conv(ak,gk) stem(Z)波形如下：12345675101520251-3（1）已知序列x[k]={1,2,3,4;k=0,1,2,3},y[k]={-1,1,2,3;k=0,1,2,3},试计算x[k]的自相关函数以及序列x[k]与y[k]的互相关函数。
编码如下：x=[1,2,3,4];kx=0:3; y=[-1,1,-2,3];ky=0:3; xf=fliplr(x); s1=conv(x,xf); s2=conv(xf,y); yf=fliplr(y); s3=conv(yf,x);k1=kx(1)+ky(1):kx(end)+ky(end); kxf=-fliplr(kx);k2=kxf(1)+ky(1):kxf(end)+ky(end); kyf=-fliplr(ky);k3=kyf(1)+kx(1):kyf(end)+kx(end); subplot(2,2,1); stem(k1,s1);xlabel('k1');ylabel('s1'); subplot(2,2,2); stem(k2,s2);xlabel('k2');ylabel('s2'); subplot(2,2,3) stem(k3,s3);xlabel('k3');ylabel('s3');波形如下：0246102030k1s 1-4-2024-10-50510k2s 2-4-2024-10-50510k3s 3M-1已知g1[t]=cos(6*pi*t),g2=cos(14*pi*t),g3=cos(26*pi*t),以抽样频率f(max)=10HZ对上述三个信号进行抽样。
实验一 时域离散信号与系统分析(实验报告)-2015
《数字信号处理》 实验报告学院 专业 电子信息工程 班级 姓名 学号 时间实验一 时域离散信号与系统分析一、实验目的1、熟悉连续信号经理想采样后的频谱变化关系，加深对时域采样定理的理解。
2、熟悉时域离散系统的时域特性，利用卷积方法观察分析系统的时域特性。
3、学会离散信号及系统响应的频域分析。
4、学会时域离散信号的MATLAB 编程和绘图。
5、学会利用MATLAB 进行时域离散系统的频率特性分析。
二、实验内容1、序列的产生（用Matlab 编程实现下列序列（数组），并用stem 语句绘出杆图。
（要求标注横轴、纵轴和标题）（1）. 单位脉冲序列x(n)=δ（n ） （2）. 矩形序列x(n)=R N (n) ,N=10nδ(n )nR N (n )图1.1 单位脉冲序列 图1.2 矩形序列(3) . x(n)=e （0.8＋3j ）n ； n 取0－15。
4n|x (n )|201321111053 陈闽焜n<x (n )/R a d图1.3 复指数序列的 模 图1.4 复指数序列的 相角（4）. x(n)=3cos （0. 25πn ＋0.3π）＋2sin （0.125πn ＋0.2π） n 取0－15。
ny (n )图1.4 复合正弦实数序列（5）. 把第（3）小题的复指数x(n)周期化，周期20点，延拓3个周期。
4m|y (m )|201321111053 陈闽焜图1.5 第（3）的20点周期延拓杆图（6）. 假设x(n)= [1，－3，2，3，-2 ]， 编程产生以下序列并绘出杆图：y(n) y(n)= x(n)－2x(n+1)+x(n-1)+x(n-3)；5201321111053 陈闽焜图1.6 y(n)序列杆图（7）、编一个用户自定义matlab 函数，名为stepshf （n0，n1，n2）实现单位阶跃序列u[n －n1]。
其中位移点数n1在起点n0和终点n2之间任意可选。
自选3个入口参数产生杆图。
1、离散时间信号的时域分析
实验一离散时间信号的时域分析一、实验目的：1、学会用MATLAB在时域中产生一些基本的离散时间信号，并对这些信号进行一些基本的运算。
2、学会使用基本的MATLAB命令，并将它们应用到简单的数字信号处理问题中。
二、实验内容：1、学习并调试本章所给的例子。
2、回答书后给出的问题。
3、实验报告仅回答奇数信号的例子。
三、实验过程：（一）单位样本和单位阶跃序列：Q1.1运行程序P1.1以产生单位样本序列u[n]并显示它：clf;n = -10:20;u = [zeros(1,10) 1 zeros(1,20)];stem(n,u);xlabel('时间 n');ylabel('振幅');title('单位抽样序列');axis([-10 20 0 1.2]);Q1.3 修改程序P1.1以产生带有延时11个样本的延迟单位样本序列ud[n]，运行修改的程序并且显示产生的序列。
clf;n = -10:20;ud = [zeros(1,21) 1 zeros(1,9)];stem(n,ud);xlabel('时间 n');ylabel('振幅');title('单位抽样序列');axis([0 30 0 1.2]);Q1.5修改程序P1.1以产生带有超前7个样本的延时单位阶跃序列sd[n]，运行修改后的程序并且显示产生的序列：clf;n = -10:20;u = [zeros(1,3) 1 ones(1,27)];stem(n,u);xlabel('时间n');ylabel('振幅');title('单位抽样序列');axis([-10 20 0 1.2]);（二）指数信号：Q1.7 那个参数控制该序列的增长或衰减率？那个参数控制该序列的振幅？答：参数c控制该序列的增长或衰减率；参数k控制该序列的振幅。
数字信号处理实验指导
实验一 离散时间信号与系统的时域分析（基础验证型）1.实验目的（1）熟悉离散时间信号的产生与基本运算。
（2）熟悉离散时间系统的时域特性。
（3）利用卷积方法观察分析系统的时域特性。
2.实验原理（1）典型离散时间信号单位样本序列（通常称为离散时间冲激或单位冲激）用[]n δ表示，其定义为1,0[]0,0n n n δ=⎧=⎨≠⎩（1.1） 单位阶跃序列用[]n μ表示，其定义为1,0[]0,0n n n μ≥⎧=⎨<⎩ （1.2） 指数序列由 []n x n A α= （1.3）给定。
其中A 和α可以是任意实数或任意复数，表示为00(),j j e A A e σωφα+==式（1.3）可改写为 0000()00[]cos()sin()n j n n n x n A e A e n j A e n σωφσσωφωφ++==+++ （1.4） 带有常数振幅的实正弦序列形如0[]cos()x n A n ωφ=+ （1.5）其中A ，0ω和φ是实数。
在式（1.4）和（1.5）中，参数A ，0ω和φ分别称为正弦序列[]x n 的振幅、角频率和初始相位。
002f ωπ=称为频率。
（2）序列的基本运算长度N 的两个序列[]x n 和[]h n 的乘积，产生长度也为N 的序列[]y n[][][]y n x n h n =⋅ （1.6）长度为N 的两个序列[]x n 和[]h n 相加，产生长度也为N 的序列[]y n[][][]y n x n h n =+ （1.7）用标量A 与长度为N 的序列[]x n 相乘，得到长度为N 的序列[]y n[][]y n A x n =⋅ （1.8）无限长序列[]x n 通过时间反转，可得到无限长序列[]y n[][]y n x n =- （1.9）无限长序列[]x n 通过M 延时，可得到无限长序列[]y n[][]y n x n M =- （1.10）若M 是一个负数，式（1.10）运算得到序列[]x n 的超前。
数字信号处理实验指导书
《数字信号处理》实验指导书编写：刘梦亭审核：司玉娟阎维和适用专业：电子信息工程电子信息科学与技术通信工程等电子信息与工程系2009年9月目录实验一：离散时间信号分析 (1)实验二：离散时间系统分析 (3)实验三：离散系统的Z域分析 (6)实验四：FFT频谱分析及应用 (9)实验五：IIR数字滤波器的设计 (12)实验六：FIR数字滤波器的设计 (16)附录： MATLAB基本操作及常用命令 (20)实验一：离散时间信号分析实验学时：2学时 实验类型：验证 实验要求：必修 一、实验目的1） 掌握离散卷积计算方法； 2） 学会差分方程的迭代解法；3） 了解全响应、零输入响应、零状态响应和初始状态的物理意义和具体求法； 二、实验内容 1、信号的加数学描述 )()()(21n x n x n x += MATLAB 实现 21X X X +=设[ x10=[1 0.7 0.4 0.1 0]; x20=[0.1 0.3 0.5 0.7 0.9 1];]2、信号的乘数学描述 )()()(21n x n x n x *= MATLAB 实现 2.1X X X *=设[ x10=[1 0.7 0.4 0.1 0]; x20=[0.1 0.3 0.5 0.7 0.9 1];]3、计算卷积用MATLAB 计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。
首先用手工计算，然后用MATLAB 编程验证。
三、实验组织运行要求1、学生在进行实验前必须进行充分的预习，熟悉实验内容；2、学生根据实验要求，读懂并理解相应的程序；3、学生严格遵守实验室的各项规章制度，注意人身和设备安全，配合和服从实验室人员管理；4、教师在学生实验过程中予以必要的辅导，独立完成实验；5、采用集中授课形式。
四、实验条件1、具有WINDOWS 98/2000/NT/XP 操作系统的计算机一台； 2.、MATLAB 编程软件。
数字信号处理实验报告
实验报告课程名称：数字信号处理授课班级：学号：姓名：指导老师：实验一离散时间信号及系统的时域分析实验类别：基础性实验1实验目的：(1）了解MA TLAB 程序设计语言的基本特点，熟悉MA TLAB软件运行环境。
(2）熟悉MA TLAB中产生信号和绘制信号的基本命令，学会用MA TLAB在时域中产生一些基本的离散时间信号，并对这些信号进行一些基本的运算。
(3）通过MA TLAB仿真一些简单的离散时间系统，并研究它们的时域特性。
(4）通过MA TLAB进行卷积运算，利用卷积方法观察分析系统的时域特性。
2. 实验报告要求●简述实验原理及目的。
●结合实验中所得给定典型序列幅频特性曲线，与理论结果比较，并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法。
●记录调试运行情况及所遇问题的解决方法。
3.实验内容：思考题：9.2.1 运行程序P9.2.1，哪个参数控制该序列的增长或衰减：哪个参数控制该序列的振幅？若需产生实指数序列，应对程序作何修改？9.2.2运行程序P9.2.1，该序列的频率是多少？怎样改变它？哪个参数控制该序列的相位？哪个参数可以控制该序列的振幅？该序列的周期是多少？9.2.3 运行程序P9.2.3，对加权输入得到的y(n)与在相同权系数下输出y1(n)和y2(n)相加得到的yt(n)进行比较，这两个序列是否相等？该系统是线性系统吗？9.2.4 假定另一个系统为y(n)=x(n)x(n-1)修改程序，计算这个系统的输出序列y1(n),y2(n)和y(n)。
比较有y(n)和yt(n)。
这两个序列是否相等？该系统是线性系统吗？(提高部分)9.2.5运行程序P9.2.4，并比较输出序列y(n)和yd(n-10)。
这两个序列之间有什么关系？该系统是时不变系统吗？9.2.6 考虑另一个系统：修改程序,以仿真上面的系统并确定该系统是否为时不变系统。
（选做）n = 0:40; D = 10;a = 3.0;b = -2;x = a*cos(2*pi*0.1*n) + b*cos(2*pi*0.4*n);[x1,n1]=sigmult(n,n,x,n)[x2,n2]=sigshift(x,n,1)[y,ny1]= sigadd(x1,n1,x2,n2)[y1,ny11]= sigshift(y,ny1,D)[sx,sn]= sigshift(x,n,D)[sx1,sn1]=sigmult(n,n, sx,sn)[sx2,sn2]=sigshift(sx,sn,1)[y2,ny2]= sigadd(sx1,sn1,sx2,sn2)D= sigadd(y1,ny11,y2,ny22)六、实验心得体会：实验时间批阅老师实验成绩实验二 FFT 实现数字滤波实验类别：提高性实验 1.实验目的(1) 通过这一实验，加深理解FFT 在实现数字滤波（或快速卷积）中的重要作用，更好的利用FFT 进行数字信号处理。
实验一 离散时间信号的时域分析
实验一 离散时间信号的时域分析实验1 序列的产生1． 目的：熟悉C 语言产生和绘制，熟悉MATLAB 中产生信号和绘制信号的基本命令。
2． 具体实验：2.1 单位样本和单位阶跃序列。
Q1.1 运行程序P1.1 ，以产生单位样本序列u[n]并显示它。
答：如图1-1所示。
Q1.2 命令clf , axis , title , xlabel 和ylabel 的作用是什么？ 答：clf ：擦除当前图形窗口中的图形。
Axis ：调整坐标轴X 轴Y 轴的范围。
Title:给绘制的图形加上标题。
Xlabel:给X 轴加上标注。
Ylabel: 给Y 轴加上标注。
Q1.3 修改程序P1.1 以产生带有延时11个样本的延迟单位样本序列ud[n]。
运行修改的程序并显示产生的序列。
答：如图1-2所示。
Q1.4修改程序P1.1 以产生单位步长序列s[n]。
运行修改的程序并显示产生的序列。
答：如图1-3所示。
Q1.5修改程序P1.1 以产生带有超前7个样本的延迟单位样本序列sd[n]。
运行修改的程序并显示产生的序列。
答：如图1-4所示。
Figure 1-2 The unit sample sequence ud[n]Figure 1-1 The shifted unit sample sequence u[n]单位样本序列公式如下所示：Time index nA m p l i t u d eUnit Sample Sequence u[n]Time index nA m p l i t u d eShifted Unit Sample Sequence ud[n]1 , n=0 1 , n=k δ[n]= δ[n-k]=0 , 0≠0 0 , 0≠kFigure 1-3 The unit step sequence s[n] Figure 1-4 The shifted unit step sequence sd[n]单位阶跃序列公式如下所示：1 , n ≥0 1 , n ≥k μ[n]= μ[n-k]=0 , n ＜0 0 , n ＜k2.2 指数信号Q1.6 运行程序P1.2 ，以产生复数值的指数序列。
数字信号处理相关MATLAB实验内容--第1章
实验1 离散时间信号的时域分析一、实验目的（1）了解MATLAB 语言的主要特点及作用；（2）熟悉MATLAB 主界面，初步掌握MATLAB 命令窗和编辑窗的操作方法；（3）学习简单的数组赋值、数组运算、绘图的程序编写；（4）了解常用时域离散信号及其特点；（5）掌握MATLAB 产生常用时域离散信号的方法。
二、知识点提示本章节的主要知识点是利用MATLAB 产生数字信号处理的几种常用典型序列、数字序列的基本运算；重点是单位脉冲、单位阶跃、正（余）弦信号的产生；难点是MATLAB 关系运算符“==、>=”的使用。
三、实验内容1. 在MATLAB 中利用逻辑关系式0==n 来实现()0n n -δ序列，显示范围21n n n ≤≤。
（函数命名为impseq(n0,n1,n2)）并利用该函数实现序列：()()()632-+-=n n n y δδ；103≤≤-nn 0212. 在MATLAB 中利用逻辑关系式0>=n 来实现()0n n u -序列，显示范围21n n n ≤≤。
（函数命名为stepseq(n0,n1,n2)）并利用该函数实现序列：()()()20522≤≤--++=n n u n u n y3. 在MATLAB 中利用数组运算符“.^”来实现一个实指数序列。
如： ()()5003.0≤≤=n n x n4. 在MATLAB 中用函数sin 或cos 产生正余弦序列，如：()()2003.0cos 553.0sin 11≤≤+⎪⎭⎫ ⎝⎛+=n n n n x πππ5. 已知()n n x 102cos 3π=，试显示()()()3,3,+-n x n x n x 在200≤≤n 区间的波形。
6. 参加运算的两个序列维数不同，已知()()6421≤≤-+=n n u n x ，()()8542≤≤--=n n u n x ，求()()()n x n x n x 21+=。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
离散时间系统的时域特性分析(数字信号处理实验1答案)
二程序框图clf;
n=0:299;
x=cos(20*pi*n/256)+cos(200*pi*n/256);
num=[0.5 0.27 0.77];
den=[1 0 0];
y=filter(num,den,x);
subplot(2,1,1);
plot(n,x);
xlabel('时间信号n');
ylabel('信号幅度');
title('输入信号');
subplot(2,1,2);
plot(y);
xlabel('时间信号n');
ylabel('信号幅度');
title('输出信号');
clf;
n=0:299;
x=cos(20*pi*n/256)+cos(200*pi*n/256); num=[0.45 0.5 0.45];
den=[1 -0.53 0.46];
y=filter(num,den,x);
subplot(2,1,1);
plot(n,x);
xlabel('时间信号n'); ylabel('信号幅度'); title('输入信号'); subplot(2,1,2);
plot(y);
xlabel('时间信号n'); ylabel('信号幅度'); title('输出信号');
n=40;
num=[0.5 0.27 0.77]; den=[1 0 0];
y=impz(num,den,n); stem(y);
xlabel('时间信号n'); ylabel('信号幅度'); title('冲激响应'); grid;
n=40;
num=[0.45 0.5 0.45]; den=[1 -0.53 0.46]; y=impz(num,den,n); stem(y);
xlabel('时间信号n'); ylabel('信号幅度'); title('冲激响应'); grid;
n=0:60;
a=2;
b=3;
x1=cos(20*pi*n/256); x2=cos(200*pi*n/256); x=a*x1+b*x2;
num=[0.45 0.5 0.45]; den=[1 -0.53 0.46];
y1=filter(num,den,x1); y2=filter(num,den,x2); y=filter(num,den,x); yt=a*y1+b*y2; subplot(2,1,1);
stem(n,y);
ylabel('信号幅度'); subplot(2,1,2);
stem(n,yt);
ylabel('信号幅度');
n=0:60;
D=10;
a=1;
b=1;
x=a*cos(20*pi*n/256)+b*cos(200*pi*n/256); xd=[zeros(1,D) x];
num=[0.45 0.5 0.45];
den=[1 -0.53 0.46];
ic=[0 0];
y=filter(num,den,x,ic);
yd=filter(num,den,xd,ic);
N=length(y);
d=y-yd(1+D:N+D);
subplot(3,1,1);
stem(n,y);
ylabel('信号幅度');
title('输出y[n]');
grid;
subplot(3,1,2);
stem(n,yd(1:length(y))); ylabel('信号幅度'); subplot(3,1,3);
stem(n,d);
xlable('时间序号n'); ylable('信号幅度');
title('差值信号'); grid;。