模电习题7章 信号的运算和处理题解
第七章信号的运算和处理
自测题
一、判断下列说法是否正确，用“√”或“×”表示判断结果。
（1）运算电路中一般均引入负反馈。（）
（2）在运算电路中，集成运放的反相输入端均为虚地。（）
（3）凡是运算电路都可利用“虚短”和“虚断”的概念求解运算关系。
（）（4）各种滤波电路的通带放大倍数的数值均大于1。（）
解：（1）√（2）×（3）√（4）×
二、现有电路：
A. 反相比例运算电路
B. 同相比例运算电路
C. 积分运算电路
D. 微分运算电路
E. 加法运算电路
F. 乘方运算电路
选择一个合适的答案填入空内。
（1）欲将正弦波电压移相＋90O，应选用。
（2）欲将正弦波电压转换成二倍频电压，应选用。
（3）欲将正弦波电压叠加上一个直流量，应选用。
（4）欲实现A u＝－100的放大电路，应选用。
（5）欲将方波电压转换成三角波电压，应选用。
（6）欲将方波电压转换成尖顶波波电压，应选用。
解：（1）C （2）F （3）E （4）A （5）C （6）D
三、填空：
（1）为了避免50Hz电网电压的干扰进入放大器，应选用滤波电路。
（2）已知输入信号的频率为10kHz～12kHz，为了防止干扰信号的混入，应选用滤波电路。
（3）为了获得输入电压中的低频信号，应选用滤波电路。
（4）为了使滤波电路的输出电阻足够小，保证负载电阻变化时滤波特性不变，应选用滤波电路。
解：（1）带阻（2）带通（3）低通（4）有源
四、已知图T7.4所示各电路中的集成运放均为理想运放，模拟乘法器的乘积系数k 大于零。试分别求解各电路的运算关系。
图T7.4
解：图（a ）所示电路为求和运算电路，图（b ）所示电路为开方运算电
路。它们的运算表达式分别为
I
3
14
2O 2
O
4
3'O 43I 12O2
O1O I34
34
21f 2I21I1f O1 )b (d 1
)1()( )a (u R kR R R u ku R R u R R u R R u t u RC
u u R R R R R R R u R u R u ?=
?-=-=-=-
=?+?+++-=?∥
习题
本章习题中的集成运放均为理想运放。
7.1 分别选择“反相”或“同相”填入下列各空内。
（1）比例运算电路中集成运放反相输入端为虚地，而比例运算电路中集成运放两个输入端的电位等于输入电压。
（2）比例运算电路的输入电阻大，而比例运算电路的输入电阻小。
（3）比例运算电路的输入电流等于零，而比例运算电路的输入电流等于流过反馈电阻中的电流。
（4）比例运算电路的比例系数大于1，而比例运算电路的比例系数小于零。
解：（1）反相，同相（2）同相，反相（3）同相，反相
（4）同相，反相
7.2填空：
（1）运算电路可实现A u＞1的放大器。
（2）运算电路可实现A u＜0的放大器。
（3）运算电路可将三角波电压转换成方波电压。
（4）运算电路可实现函数Y＝aX1＋bX2＋cX3，a、b和c均大于零。
（5）运算电路可实现函数Y＝aX1＋bX2＋cX3，a、b和c均小于零。
（6）运算电路可实现函数Y＝aX2。
解：（1）同相比例（2）反相比例（3）微分（4）同相求和（5）反相求和（6）乘方
7.3 电路如图P7.3所示，集成运放输出电压的最大幅值为±14V ，填表。
图P 7.3
u I /V
0.1 0.5 1.0 1.5 u O 1/V u O 2/V
解：u O 1＝(－R f /R ) u I ＝－10 u I ，u O 2＝(1+R f /R ) u I ＝11 u I 。当集成运放工作到非线性区时，输出电压不是＋14V ，就是－14V 。
u I /V 0.1 0.5 1.0 1.5 u O 1/V －1 －5 －10 －14 u O 2/V
1.1
5.5
11
14
7.4 设计一个比例运算电路， 要求输入电阻R i ＝20k Ω， 比例系数为－100。
解：可采用反相比例运算电路，电路形式如图P7.3(a)所示。R ＝20k Ω，R f ＝2M Ω。
7.5 电路如图P7.5所示，试求： （1）输入电阻； （2）比例系数。
解：由图可知R i ＝50k Ω，u M ＝－2u I 。 342R R R i i i += 即 3
O
M 4M 2M R u u R u R u -+=-
输出电压 I M O 10452u u u -== 图P 7.5
7.6 电路如图P7.5所示，集成运放输出电压的最大幅值为±14V ，u I 为2V 的直流信号。分别求出下列各种情况下的输出电压。 （1）R 2短路；（2）R 3短路；（3）R 4短路；（4）R 4断路。
解：（1）V 4 2I 1
3
O -=-=-=u R R u （2）V 4 2I 1
2
O -=-=-
=u R R u （3）电路无反馈，u O ＝－14V （4）V 8 4I 1
3
2O -=-=+-=u R R R u
7.7 电路如图P7.7所示，T 1、T 2和T 3的特性完全相同，填空： （1）I 1≈ mA ，I 2≈ mA ；
（2）若I 3≈0.2mA ，则R 3≈ k Ω。
图P 7.7
解：（1）1，0.4；（2）10。
7.8 试求图P7.8所示各电路输出电压与输入电压的运算关系式。
图P 7.8
解：在图示各电路中，集成运放的同相输入端和反相输入端所接总电阻均相等。各电路的运算关系式分析如下：
（a ）13I2I1I33f I22f I11f O 522u u u u R R
u R R u R R u +--=?+?-?-= （b ）13I2I1I33
f I22f I11f O 1010u u u u R R
u R R u R R u ++-=?+?+?-= （c ）)( 8)(I1I2I1I21
f
O u u u u R R u -=-=
（d ）I44
f I33f I22f I11f O u R R
u R R u R R u R R u ?+?+?-?-
= 1413I2I1402020u u u u ++--=
7.9 在图P7.8所示各电路中，是否对集成运放的共模抑制比要求较高，为什么？
解：因为均有共模输入信号，所以均要求用具有高共模抑制比的集成运放。
7.10 在图P7.8所示各电路中，集成运放的共模信号分别为多少？要求写出表达式。
解：因为集成运放同相输入端和反相输入端之间净输入电压为零，所以它们的电位就是集成运放的共模输入电压。图示各电路中集成运放的共模信号分别为
（a ）I3IC u u = （b ）I3I2I3322I2323IC 111
1110u u u R R R u R R R u +=?++?+=
（c ）I2I2f 1f IC 9
8
u u R R R u =?+=
（d ）I4I3I4433I3434IC 41
1
4140u u u R R R u R R R u +=?++?+=
7.11 图P7.11所示为恒流源电路，已知稳压管工作在稳压状态，试求负载电阻中的电流。
图P 7.11
解：6.02
Z 2P L ===R U R u I mA
7.12 电路如图P7.12所示。
（1）写出u O 与u I 1、u I 2的运算关系式；
（2）当R W 的滑动端在最上端时，若u I 1＝10mV ，u I 2＝20mV ，则u O ＝？ （3）若u O 的最大幅值为±14V ，输入电压最大值 u I 1ma x ＝10mV ，u I 2ma x ＝20mV ，最小值均为0V ，则为了保证集成运放工作在线性区，R 2的最大值为多少？
图P 7.12
解：（1）A 2同相输入端电位 )( 10)(I1I2I1I2f
N2P2u u u u R
R u u -=-== 输出电压 ))(1(10)1(I1I212P212O u u R R
u R R u -+=?+= 或 )(10I1I21
W
O u u R R u -??
= （2）将u I 1＝10mV ，u I 2＝20mV 代入上式，得u O ＝100mV
（3）根据题目所给参数，)(I1I2u u -的最大值为20mV 。若R 1为最小值，则为保证集成运放工作在线性区， )(I1I2u u -＝20mV 时集成运放的输出电压应为＋14V ，写成表达式为
142010
10)(10min
1I1I2min 1W O =??=-??
=R u u R R u 故 R 1m i n ≈143Ω
R 2ma x ＝R W －R 1m i n ≈（10－0.143）k Ω≈9.86 k Ω
7.13 分别求解图P7.13所示各电路的运算关系。
图P 7.13
解：图（a ）所示为反相求和运算电路；图（b ）所示的A 1组成同相比例运算电路，A 2组成加减运算电路；图（c ）所示的A 1、A 2、A 3均组成为电压跟随器电路，A 4组成反相求和运算电路。
（a ）设R 3、R 4、R 5的节点为M ，则
))(( )(
2
I2
1I15434344M O 5
M
2I21I15342
I2
1I13M R u R u R R R R R R i u u R u R u R u i i i R u R u R u R R R R ++
+-=-=-+=
-=+-=
（b ）先求解u O 1，再求解u O 。
))(1(
)1()1( )1()1(I1I24
5
I245I11345I24
5O145O I11
3
O1u u R R u R R
u R R R R u R R
u R R u u R R u -+
=+++-=++-
=+=
（c ）A 1、A 2、A 3的输出电压分别为u I 1、u I 2、u I 3。由于在A 4组成的反相求和运算电路中反相输入端和同相输入端外接电阻阻值相等，所以 )( 10)(I3I2I1I3I2I11
4
O u u u u u u R R u ++=++=
7.14 在图P7.14（a ）所示电路中，已知输入电压u I 的波形如图（b ）所示，当t ＝0时u O ＝0。试画出输出电压u O 的波形。
图P 7.14
解：输出电压的表达式为 )(d 11O I O 2
1
t u t u RC u t t +-=?
当u I 为常量时
)
()(100 )()(10101
)()(1
1O 12I 1O 12I 7
5
112I O t u t t u t u t t u t u t t u RC
u O +-=+-?-=+--
=-－ 若t ＝0时u O ＝0，则t ＝5ms 时 u O ＝－100×5×5×10
－
3V ＝－2.5V 。
当t＝15mS时
u O＝[－100×(－5)×10×10－3＋（－2.5）]V＝2.5V。
因此输出波形如解图P7.14所示。
解图P7.14
7.15已知图P7.15所示电路输入电压u I的波形如图P7.4（b）所示，且当t＝0时u O＝0。试画出输出电压u O的波形。
图P7.15
解图P7.15
解：输出电压与输入电压的运算关系为u O＝100u I（t2－t1）＋u I－u C（t1），波形如解图P7.15所示。
7.16 试分别求解图P7.16所示各电路的运算关系。
图P 7.16
解：利用节点电流法，可解出各电路的运算关系分别为： （a ） t u u t u C R u R R u d 100d 1I I I 1I 12O ??--=--
= （b ） I I 3I 21I 1O 2d d 10d d u t
u
u C C t u RC u --=--=- （c ） t u t u RC
u d 10d 1I 3I O ??==
（d ） t u u t R u R u C u d )5.0(100d )(1I2I12
I21I1O +-=+-=??
7.17 在图P7.17所示电路中，已知R 1＝R ＝R '＝100k Ω，R 2＝R f ＝100k
Ω，C ＝1μF 。
图P 7.17
（1）试求出u O 与 u I 的运算关系。
（2）设t ＝0时u O ＝0，且u I 由零跃变为－1V ，试求输出电压由零上升到＋6V 所需要的时间。
解：（1）因为A 1的同相输入端和反相输入端所接电阻相等，电容上的电压u C ＝u O ，所以其输出电压
I O O 2
f I 1f O1u u u R R
u R R u -=?+?-= 电容的电流
R
u
R u u i I O O1C -=-= 因此，输出电压 t u t u RC t i C u d 10d 1d 1I I C O ??
?-=-==
（2）u O ＝－10u I t 1＝[－10×（－1）×t 1]V ＝6V ，故t 1＝0.6S 。即经0.6秒输出电压达到6V 。
7.18 试求出图P7.18所示电路的运算关系。
图P 7.18
解：设A 2的输出为u O 2。因为R 1的电流等于C 的电流，又因为A 2组成以u O 为输入的同相比例运算电路，所以
???
-=
=+
=-=-
=t
u u u u R R u t u t u C
R u d 2)1( d 2d 1
I O O O 3
2
O2I I 1O2
7.19 在图P7.19所示电路中，已知u I 1＝4V ，u I 2＝1V 。回答下列问题：
图P 7.19
（1）当开关S 闭合时，分别求解A 、B 、C 、D 和u O 的电位； （2）设t ＝0时S 打开，问经过多长时间u O ＝0？
解：（1）U A＝7V，U B＝4V，U C＝1V，U D＝－2V，u O＝2 U D＝－4V。
（2）因为u O＝2u D－u O3，2 u D＝－4V，所以u O3＝－4V时，u O才为零。即
mS
6.
28
4
7
10
10
50
1
1
6
3
A
1
O3
≈
-
=
?
?
?
?
-
=
?
?
-
=
-
t
t
t
u
C
R
u
7.20画出利用对数运算电路、指数运算电路和加减运算电路实现除法运算的原理框图。
解：答案如解图7.20所示。
解图7.20
7.21 为了使图P7.21所示电路实现除法运算，
（1）标出集成运放的同相输入端和反相输入端；
（2）求出u O和u I1、u I2的运算关系式。
图P7.21
解：（1）为了保证电路引入负反馈，A的上端为“－”，下端为“＋”。
（2）根据模拟乘法器输出电压和输入电压的关系和节点电流关系，可得
)1.0( I2O f
'
O
f
I1I2
O 'O u u R R R
u R R R u
u
ku u -?+=
+=
=
所以
I2
1
I f O )(10u u R R R u ?+-=
7.22 求出图P7.22所示各电路的运算关系。
图P 7.22
解：电路（a ）实现求和、除法运算，电路（b ）实现一元三次方程。它们的运算关系式分别为
)
( )(
)a (2
I21I1I33O I3O 2
I2
1I13'
O R u
R u ku R u u ku R u R u R u +-==+-=
I 1
4
3I 2342I 24O )b (u R R u k R R ku R R u ---
=
7.23 利用图7.3.9所示方框图的思路，分别设计5次方运算电路和5次幂运算电路。
解：方框图如图7.3.9所示，N ＝5时为5次方电路；N ＝0.2时为5次幂电路。
7.24 在下列各种情况下，应分别采用哪种类型（低通、高通、带通、带阻）的滤波电路。
（1）抑制50Hz 交流电源的干扰；
（2）处理具有1Hz 固定频率的有用信号； （3）从输入信号中取出低于2kHz 的信号； （4）抑制频率为100kHz 以上的高频干扰。
解：（1）带阻滤波器 （2）带通滤波器 （3）低通滤波器 （4）低通滤波器
7.25试说明图P7.25所示各电路属于哪种类型的滤波电路，是几阶滤波电路。
图P7.25
解：图（a）所示电路为一阶高通滤波器。
图（b）所示电路二阶高通滤波器。
图（c）所示电路二阶带通滤波器。
图（d）所示电路二阶带阻滤波器。
7.26 设一阶LPF 和二阶HPF 的通带放大倍数均为2，通带截止频率分别为2kHz 和100Hz 。试用它们构成一个带通滤波电路，并画出幅频特性。
解：低通滤波器的通带截止频率为2kHz ，高通滤波器的通带截止频率为100Hz 。将两个滤波器串联，就构成一个带通滤波电路。其通带放大倍数为
4p
=u A 通带增益为
12lg 20p
≈u A 幅频特性如解图P7.26所示。
解图P 7.26
7.27 分别推导出图P7.27所示各电路的传递函数，并说明它们属于哪种类型的滤波电路。
图P 7.27
解：利用节点电流法可求出它们的传递函数。 在图（a ）所示电路中 11)
( 1212C
sR C
sR sC
R R s A u +-
=+
-
=
故其为高通滤波器。 在图（b ）所示电路中
C
sR R R R sC R sC R s A u 21212211)
1
(1)( +?-=+?
-=
故其为低通滤波器。
7.28 试分析图P7.28所示电路的输出u O 1、u O 2和u O 3分别具有哪种滤波特性（LPF 、HPF 、BPF 、BEF ）？
图P 7.28
解：以u O 1为输出是高通滤波器，以u O 2为输出是带通滤波器，以u O 3为输出是低通滤波器。
数字信号处理实验(吴镇扬)答案-2
（1） 观察高斯序列的时域和幅频特性，固定信号)(n x a 中参数p=8,改变q 的 值，使q 分别等于2、4、8，观察他们的时域和幅频特性，了解当q 取不同值时，对信号序列的时域和幅频特性的影响；固定q=8，改变p,使p 分别等于8、13、14，观察参数p 变化对信号序列的时域和幅频特性的影响，注意p 等于多少时会发生明显的泄漏现象，混叠是否也随之出现？记录实验中观察到的现象，绘出相应的时域序列和幅频特性曲线。 ()() ?????≤≤=-其他0150,2n e n x q p n a 解：程序见附录程序一： P=8,q 变化时： t/T x a (n ) k X a (k ) t/T x a (n ) p=8 q=4 k X a (k ) p=8 q=4 t/T x a (n ) p=8 q=8 k X a (k ) p=8 q=8 幅频特性 时域特性
t/T x a (n ) p=8 q=8 k X a (k ) p=8 q=8 t/T x a (n ) 5 10 15 k X a (k ) p=13 q=8 t/T x a (n ) p=14 q=8 5 10 15 k X a (k ) p=14 q=8 时域特性幅频特性 分析： 由高斯序列表达式知n=p 为期对称轴； 当p 取固定值时，时域图都关于n=8对称截取长度为周期的整数倍，没有发生明显的泄漏现象；但存在混叠，当q 由2增加至8过程中，时域图形变化越来越平缓，中间包络越来越大，可能函数周期开始增加，频率降低，渐渐小于fs/2,混叠减弱； 当q 值固定不变，p 变化时，时域对称中轴右移，截取的时域长度渐渐地不再是周期的整数倍，开始无法代表一个周期，泄漏现象也来越明显，因而图形越来越偏离真实值， p=14时的泄漏现象最为明显，混叠可能也随之出现；
数字信号处理实验一
实验一 离散时间信号分析 班级 信息131班 学号 201312030103 姓名 陈娇 日期 一、实验目的 掌握两个序列的相加、相乘、移位、反褶、卷积等基本运算。 二、实验原理 1．序列的基本概念 离散时间信号在数学上可用时间序列)}({n x 来表示，其中)(n x 代表序列的第n 个数字，n 代表时间的序列，n 的取值范围为+∞<<∞-n 的整数，n 取其它值)(n x 没有意义。离散时间信号可以是由模拟信号通过采样得到，例如对模拟信号)(t x a 进行等间隔采样，采样间隔为T ，得到)}({nT x a 一个有序的数字序列就是离散时间信号，简称序列。 2．常用序列 常用序列有：单位脉冲序列（单位抽样）） （n δ、单位阶跃序列)(n u 、矩形序列)(n R N 、实指数序列、复指数序列、正弦型序列等。 3．序列的基本运算 序列的运算包括移位、反褶、和、积、标乘、累加、差分运算等。 4．序列的卷积运算 ∑∞ -∞==-= m n h n x m n h m x n y )(*)()()()( 上式的运算关系称为卷积运算，式中代表两个序列卷积运算。两个序列的卷积是一个序列与另一个序列反褶后逐次移位乘积之和，故称为离散卷积，也称两序列的线性卷积。其计算的过程包括以下4个步骤。 （1）反褶：先将)(n x 和)(n h 的变量n 换成m ，变成)(m x 和)(m h ，再将)(m h 以纵轴为对称轴反褶成)(m h -。
（2）移位：将)(m h -移位n ，得)(m n h -。当n 为正数时，右移n 位；当n 为负数时，左移n 位。 （3）相乘：将)(m n h -和)(m x 的对应点值相乘。 （4）求和：将以上所有对应点的乘积累加起来，即得)(n y 。 三、主要实验仪器及材料 微型计算机、Matlab6.5 教学版、TC 编程环境。 四、实验内容 （1）用Matlab 或C 语言编制两个序列的相加、相乘、移位、反褶、卷积等的程序； （2）画出两个序列运算以后的图形； （3）对结果进行分析； （4）完成实验报告。 五、实验结果 六、实验总结
7章 信号的运算和处理题解
（4）各种滤波电路的通带放大倍数的数值均大于1。（×） 二、现有电路： A. 反相比例运算电路 B. 同相比例运算电路 C. 积分运算电路 D. 微分运算电路 E. 加法运算电路 F. 乘方运算电路 选择一个合适的答案填入空内。 （1）欲将正弦波电压移相＋90O，应选用 C 。 （2）欲将正弦波电压转换成二倍频电压，应选用 F 。 （3）欲将正弦波电压叠加上一个直流量，应选用 E 。 （4）欲实现A u＝－100的放大电路，应选用A 。 （5）欲将方波电压转换成三角波电压，应选用 C 。 （6）欲将方波电压转换成尖顶波波电压，应选用 D 。 （1）为了避免50Hz电网电压的干扰进入放大器，应选用带阻滤波电路。 （2）已知输入信号的频率为10kHz～12kHz，为了防止干扰信号的混入，应选用带通滤波电路。 （3）为了获得输入电压中的低频信号，应选用低通滤波电路。 （4）为了使滤波电路的输出电阻足够小，保证负载电阻变化时滤波特性不变，应选用有源滤波电路。 四、已知图T7.4所示各电路中的集成运放均为理想运放，模拟乘法器的乘积系数k大于零。试分别求解各电路的运算关系。
图T7.4 解：图（a ）所示电路为求和运算电路，图（b ）所示电路为开方运算电路。它们的运算表达式分别为 I 3142O 2O 4 3'O 43I 12O2O1O I34 3421f 2I21I1f O1 )b (d 1 )1()( )a (u R kR R R u ku R R u R R u R R u t u RC u u R R R R R R R u R u R u ?= ?-=-=-=- =?+?+++-=?∥ 本章习题中的集成运放均为理想运放。 7.1 分别选择“反相”或“同相”填入下列各空内。 （1）反相 比例运算电路中集成运放反相输入端为虚地，而同相 比例运算电路中集成运放两个输入端的电位等于输入电压。 （2）同相比例运算电路的输入电阻大，而反相比例运算电路的输入电阻小。 （3）同相 比例运算电路的输入电流等于零，而 反相 比例运算电路的输入电流等于流过反馈电阻中的电流。 （4）同相 比例运算电路的比例系数大于1，而反相 比例运算电路的比例系数小于零。 （6） 乘方 运算电路可实现函数Y ＝aX 2。
数字信号处理实验作业
实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。
数字信号处理实验报告(实验1_4)
实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值； （2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算 （1）矩阵的乘法 已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B
（2）矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；
（5）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法 （1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；
（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。
5、MATLAB基本绘图命令 （1）绘制余弦曲线 y=cos(t)，t∈[0，2π]
（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求： （a）线形为点划线、颜色为红色、数据点标记为加号； （b）坐标轴控制：显示围、刻度线、比例、网络线 （c）标注控制：坐标轴名称、标题、相应文本； >> clear;
数字信号处理实验
实验一 离散傅里叶变换（DFT ）对确定信号进行谱分析 一．实验目的 1．加深对DFT 算法原理和基本性质的理解。 2．熟悉DFT 算法和原理的编程方法。 3．学习用DFT 对信号进行谱分析的方法，了解可能出现的误差及其原因，以便在实际中正确利用。 二．实验原理 一个连续信号)(t x a 的频谱可以用其傅里叶变换表示，即 dt e t x j X t j a a Ω-∞ ∞ -? = Ω)()( 若对)(t x a 进行理想采样可得采样序列 )(|)()(nT x t x n x a nT t a === 对)(n x 进行DTFT ，可得其频谱为： ∑∞ -∞ =-= n n j j e n x e X ωω )()( 其中数字频率ω与模拟频率Ω的关系为： s f T Ω = Ω=ω )(n x 的DFT 为∑∞ -∞ =-= n nk N j e n x k X π 2)()( 若)(t x a 是限带信号，且在满足采样定理的条件下，)(ω j e X 是)(Ωj X a 的周期延拓， )(k X 是)(ωj e X 在单位圆上的等间隔采样值，即k N j e X k X πωω2| )()(= =。 为在计算机上分析计算方便，常用)(k X 来近似)(ω j e X ，这样对于长度为N 的有限 长序列（无限长序列也可用有限长序列来逼近），便可通过DFT 求其离散频谱。 三．实验内容 1．用DFT 对下列序列进行谱分析。 （1）)()04.0sin(3)(100n R n n x π=
1 （2）]0,0,0,0,0,0,0,0,1,1,1,1[)(=n x 2．为了说明高密度频谱和高分辨率频谱之间的区别，考察序列 )52.0cos()48.0cos()(n n n x ππ+= （1）当0≤n ≤10时，确定并画出x(n)的离散傅里叶变换。 （2）当0≤n ≤100时，确定并画出x(n)的离散傅里叶变换。 四．实验结果 1. （1） （2）
数字信号处理实验作业
实验5 抽样定理 一、实验目的： 1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和插公式的编程方法。 二、实验原理： 1、时域抽样与信号的重建 （1）对连续信号进行采样 例5-1 已知一个连续时间信号sin sin(),1Hz 3 ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。 程序清单如下： %分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2; f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3; fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end 程序运行结果如图5-1所示：
原连续信号和抽样信号 图5-1 （2）连续信号和抽样信号的频谱 由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。 例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。 程序清单如下： dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm; t=-2:dt:2;N=length(t); f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 程序运行结果如图5-2所示。 由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的围，频谱出现了镜像对称的部分。
数字信号处理实验答案完整版
数字信号处理实验答案 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≤15 b) x(n)=e+3j)n 0≤n≤15 c) x(n)=3cosπn+π)+2sinπn+π) 0≤n≤15 d) 将c)中的x(n)扩展为以16为周期的函数x(n)=x(n+16),绘出四个周期。
数字信号处理实验1认识实验
实验1认识实验-MATLAB语言上机操作实践 一、实验目的 ㈠了解MATLAB语言的主要特点、作用。 ㈡学会MATLAB主界面简单的操作使用方法。 ㈢学习简单的数组赋值、运算、绘图、流程控制编程。 二、实验原理 ㈠简单的数组赋值方法 MATLAB中的变量和常量都可以是数组（或矩阵）,且每个元素都可以是复数。 在MATLAB指令窗口输入数组A=[1 2 3；4 5 6；7 8 9]，观察输出结果。然后，键入：A（4，2）= 11 键入：A (5，：) = [-13 -14 -15] 键入：A（4，3）= abs (A(5，1)) 键入：A ([2，5]，：) = [ ] 键入：A/2 键入：A (4，：) = [sqrt(3) (4+5)/6*2 –7] 观察以上各输出结果。将A式中分号改为空格或逗号，情况又如何？请在每式的后面标注其含义。 2．在MATLAB指令窗口输入B=[1+2i，3+4i；5+6i ，7+8i], 观察输出结果。 键入：C=[1，3；5，7]+[2，4；6，8]*i，观察输出结果。 如果C式中i前的*号省略，结果如何？ 键入：D = sqrt (2+3i) 键入：D*D 键入：E = C’, F = conj(C), G = conj(C)’ 观察以上各输出结果, 请在每式的后面标注其含义。 3．在MATLAB指令窗口输入H1=ones(3,2)，H2=zeros(2,3)，H3=eye(4)，观察输出结果。 ㈡、数组的基本运算 1．输入A=[1 3 5]，B= [2 4 6]，求C=A+B，D=A-2，E=B-A 2．求F1=A*3，F2=A.*B，F3=A./B，F4=A.\B, F5=B.\A, F6=B.^A, F7=2./B, F8=B.\2 *3．求B'，Z1=A*B’，Z2=B’*A 观察以上各输出结果,比较各种运算的区别，理解其含义。 ㈢、常用函数及相应的信号波形显示 例1：显示曲线f(t)=2sin(2πt)，（t>0） ⅰ点击空白文档图标（New M-file），打开文本编辑器。 ⅱ键入：t=0:0.01:3; (1) f=2*sin(2*pi*t); (2) plot(t,f); title(‘f（t）-t曲线’)； xlabel（‘t’），ylabel（‘f（t）’）；
数字信号处理实验(吴镇扬)答案-4
实验四 有限长单位脉冲响应滤波器设计 朱方方 0806020433 通信四班 (1) 设计一个线性相位FIR 高通滤波器，通带边界频率为0.6π，阻带边界频率为0.4π，阻 带衰减不小于40dB 。要求给出h(n)的解析式，并用MATLAB 绘出时域波形和幅频特性。 解： (1) 求数字边界频率: 0.6 , .c r ωπωπ== (2) 求理想滤波器的边界频率: 0.5n ωπ= (3) 求理想单位脉冲响应: []d s i n ()s i n [()] () ()1n n n n n n h n n παωαα παωα π?-- -≠??-=? ? -=?? (4) 选择窗函数。阻带最小衰减为-40dB ，因此选择海明窗(其阻带最小衰减为-44dB)；滤 波器的过渡带宽为0.6π-0.4π=0.2π，因此 6.21 0.231 , 152 N N N ππα-=?=== (5) 求FIR 滤波器的单位脉冲响应h(n): []31d sin (15)sin[0.5(15)] 1cos ()15()()()15(15)1 15 n n n R n n h n w n h n n n ππππ?---????-? ?≠? ???==-???? ? ?=? 程序： clear; N=31; n=0:N-1; hd=(sin(pi*(n-15))-sin(0.5*pi*(n-15)))./(pi *(n-15)); hd(16)=0.5; win=hanning(N); h=win'.*hd; figure; stem(n,h); xlabel('n'); ylabel('h(n)'); grid; title('FIR 高通滤波单位脉冲响应h(n)'); [H,w]=freqz(h,1); H=20*log10(abs(H)); figure;3 plot(w/pi,H); axis([0 1 -100 10]); xlabel('\omega/\pi'); ylabel('幅度/dB'); grid; title('FIR 高通滤波器，hanning 窗，N=31');
数字信号处理实验三
实验三：离散LSI 系统的频域分析 一、实验内容 2、求以下各序列的z 变换： 12030() ()sin() ()sin()n an x n na x n n x n e n ωω-=== 程序清单如下： syms w0 n z a; x1=n*a^n;X1=ztrans(x1) x2=sin(w0*n);X2=ztrans(x2) x3= exp(-a*n)*sin(w0*n);X3=ztrans(x3) 程序运行结果如下： X1 =z/(a*(z/a - 1)^2) X2 =(z*sin(w0))/(z^2 - 2*cos(w0)*z + 1) X3 =(z*exp(a)*sin(w0))/(exp(2*a)*z^2 - 2*exp(a)*cos(w0)*z + 1) 3、求下列函数的逆z 变换 0 312342 1 1() () () ()() 1j z z z z X z X z X z X z z a z a z e z ω---= = = = ---- 程序清单如下： syms w0 n z a; X1=z/(z-a);x1=iztrans(X1) X2= z/(a-z)^2;x2=iztrans(X2) X3=z/ z-exp(j*w0);x3=iztrans(X3) X4=(1-z^-3)/(1-z^-1);x4=iztrans(X4) 程序运行结果如下： x1 =a^n x2 =n*a^n/a 课程名称 数字信号 实验成绩 指导教师 实 验 报 告 院系 信息工程学院 班级 学号 姓名 日期
x3 =charfcn[0](n)-iztrans(exp(i*w0),w0,n) x4 =charfcn[2](n)+charfcn[1](n)+charfcn[0](n) 4、求一下系统函数所描述的离散系统的零极点分布图，并判断系统的稳定性 （1） (0.3)()(1)(1) z z H z z j z j -= +-++ z1=[0,0.3]';p1=[-1+j,-1-j]';k=1; [b1,a1]=zp2tf(z1,p1,k); subplot(1,2,1);zplane(z1,p1); title('极点在单位圆外); subplot(1,2,2);impz(b1,a1,20); 由图可见:当极点位于单位圆内，系统的单位序列响应随着频率的增大而收敛；当极点位于单位圆上，系统的单位序列响应为等幅振荡；当极点位于单位圆外，系统的单位序列响应随着频率的增大而发散。由此可知系统为不稳定系统。 -1 -0.5 00.51 -2 -1.5-1-0.500.511.5 2Real Part I m a g i n a r y P a r t 极点在单位圆外 n (samples) A m p l i t u d e Impulse Response
信号的运算和处理
第七章 信号的运算和处理 【本章主要内容】本章主要讲述基本运算电路和有源滤波电路。 【本章学时分配】本章分为2讲，每讲2学时。 第二十讲 运算电路概述和基本运算电路 一、主要内容 1、比例运算电路 分析方法，利用虚短、虚断的概念和基尔霍夫电流定理列出放大倍数表达式。 1) 反相比例运算电路 （1）电路的组成如图7.2.1所示。 （2）电路的放大倍数及特点 由分析得电路的放大倍数为 1 u R R A f -= 特点 ①输入信号接入反相输入端，u N 点虚地，其输出信号与输入信号反相。 ②电路不存在共模信号。 ③放大倍数可以大于1，可以小于1，也可以等于0。 ④因为电路引入电压并联负反馈，故电路的输入阻抗较低，即R i =R 1。 2) 同相比例运算电路 （1）电路的组成如图7.2.2所示。 （2）电路的放大倍数及特点 由分析得电路的放大倍数为 1 u R R 1A f += 特点 ①输入信号接入同相输入端，故其输出信号与输入信号同相。 ②电路存在共模信号，故应选用共模抑制比高的集成运放。 ③放大倍数只能大于或等于1。 ④因为电路引入电压串联负反馈，故其输入阻抗很高。 2、加减运算电路 分析方法，利用虚短、虚断的概念、结电电压法或叠加定理列出输出方程。 1) 反相求和运算电路 （1）电路的组成如下图所示
R u 1 u 2u o （2）电路的分析及特点 电路的输出表达式为 ??? ? ??+-=22 11 o u R R u R R u f f 电路的特点与反相比例运算电路的特点类似。 2) 同相求和运算电路 （1）电路的组成如下图所示 R 3u o （2）电路的分析及特点 电路的输出表达式为 ???? ??+++???? ? ? + =23 22132 31o u R R R u R R R R R 1u f 电路的特点与同相比例运算电路的特点类似。 3) 加减运算电路 （1）电路的组成如下图所示 R u 1u 2 u o （2）电路的分析及特点 电路的输出表达式为
数字信处理上机实验答案全
数字信处理上机实验答 案全 Document number【SA80SAB-SAA9SYT-SAATC-SA6UT-SA18】
第十章 上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一系统响应及系统稳定性。 实验二时域采样与频域采样。 实验三用FFT对信号作频谱分析。 实验四 IIR数字滤波器设计及软件实现。 实验五 FIR数字滤波器设计与软件实现 实验六应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱函数filter函数。也可以用MATLAB语言的工具箱函数conv函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。 系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。 系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。 实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的[19]。系统的稳态输出是指当∞ n时，系统的输出。如果系统稳定，信号加入 → 系统后，系统输出的开始一段称为暂态效应，随n的加大，幅度趋于稳定，达到稳态输出。 注意在以下实验中均假设系统的初始状态为零。 3．实验内容及步骤
数字信号处理实验4
数字信号处理实验四 第一题结果： （1）没有增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 %H(3,13) = 0.75;H(5,11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线
（2）增加过渡点 源码如下： N = 15; H = [1 1 1 0.5 zeros(1,7) 0.5 1 1 1]; %确定抽样点的幅度大小 H(3) = 0.75;H(13) = 0.75;H(5) = 0.25;H(11) = 0.25; %设置过渡点 k = 0:N-1; A = exp(-j*pi*k*(N-1)/N); %抽样点相位大小 HK = H.*A; %求抽样点的H(k) hn = ifft(HK,N); %求出FIR的单位冲激响应h(n) freqz(hn,1,256); %画出幅频相频曲线figure(2); stem(real(hn),'.'); %绘制单位冲激响应的实部 line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); 单位脉冲响应曲线 幅频和相频特性曲线 第二题结果：
章 信号的运算和处理题解 第四版模电答案
第七章信号的运算和处理 自测题 一、现有电路： A. 反相比例运算电路 B. 同相比例运算电路 C. 积分运算电路 D. 微分运算电路 E. 加法运算电路 F. 乘方运算电路 选择一个合适的答案填入空内。 （1）欲将正弦波电压移相＋90O，应选用。 （2）欲将正弦波电压转换成二倍频电压，应选用。 （3）欲将正弦波电压叠加上一个直流量，应选用。 （4）欲实现A u＝－100的放大电路，应选用。 （5）欲将方波电压转换成三角波电压，应选用。 （6）欲将方波电压转换成尖顶波波电压，应选用。 解：（1）C （2）F （3）E （4）A （5）C （6）D 二、填空： （1）为了避免50Hz电网电压的干扰进入放大器，应选用滤波电路。 （2）已知输入信号的频率为10kHz～12kHz，为了防止干扰信号的混入，应选用滤波电路。 （3）为了获得输入电压中的低频信号，应选用滤波电路。 （4）为了使滤波电路的输出电阻足够小，保证负载电阻变化时滤波特性不变，应选用滤波电路。 解：（1）带阻（2）带通（3）低通（4）有源 三、已知图T7.3所示各电路中的集成运放均为理想运放，模拟乘法器的乘积系数k大于零。试分别求解各电路的运算关系。 图T7.3 解：图（a）所示电路为求和运算电路，图（b）所示电路为开方运算电路。它
们的运算表达式分别为 习题 本章习题中的集成运放均为理想运放。 7.1填空： （1）运算电路可实现A u＞1的放大器。 （2）运算电路可实现A u＜0的放大器。 （3）运算电路可将三角波电压转换成方波电压。 （4）运算电路可实现函数Y＝aX1＋bX2＋cX3，a、b和c均大于零。 （5）运算电路可实现函数Y＝aX1＋bX2＋cX3，a、b和c均小于零。 （6）运算电路可实现函数Y＝aX2。 解：（1）同相比例（2）反相比例（3）微分（4）同相求和 （5）反相求和（6）乘方 7.2 电路如图P7.2所示，集成运放输出电压的最大幅值为±14V，填表。 图P7.2 u I/V 0.1 0.5 1.0 1.5 u O1/V u O2/V 解：u O1＝(－R f /R) u I＝－10 u I，u O2＝(1+R f /R ) u I＝11 u I。当集成运放工作到非线性区时，输出电压不是＋14V，就是－14V。 u I/V 0.1 0.5 1.0 1.5 u O1/V －1 －5 －10 －14 u O2/V 1.1 5.5 11 14 7.3设计一个比例运算电路，要求输入电阻R i＝20kΩ，比例系数为－100。 解：可采用反相比例运算电路，电路形式如图P7.2(a)所示。R＝20kΩ，R f＝2M Ω。 7.4电路如图P7.4所示，试求： （1）输入电阻； （2）比例系数。 解：由图可知R i＝50kΩ，u M＝－2u I。
数字信号处理第二章上机题作业
数字信号处理作业实验题报告 第一章16.(1) 实验目的： 求解差分方程所描述的系统的单位脉冲响应和单位阶跃响应。 实验要求： 运用matlab求出y(n)=0.6y(n-1)-0.08y(n-2)+x(n)的单位脉冲响应和单位阶跃响应的示意图。 源程序： B1=1;A1=[1, -0.6, 0.08]; ys=2; %设差分方程 xn=[1, zeros(1, 20)]; %xn=单位脉冲序列，长度N=31 xi=filtic(B1, A1, ys); hn1=filter(B1, A1, xn, xi); %求系统输出信号hn1 n=0:length(hn1)-1; subplot(2, 1, 1);stem(n, hn1, '.') title('单位脉冲响应'); xlabel('n');ylabel('h(n)') xn=ones(1, 20); sn1=filter(B1, A1, xn, xi); %求系统输出信号sn1 n=0:length(sn1)-1; Subplot(2, 1, 2); stem(n, sn1, '.') title('单位阶跃响应'); xlabel('n'); ylabel('s(n)')
运行结果： 实验分析： 单位脉冲响应逐渐趋于0，阶跃响应保持不变，由此可见，是个稳定系统。
第二章31题 实验目的： 用matlab判断系统是否稳定。 实验要求： 用matlab画出系统的极，零点分布图，输入单位阶跃序列u(n)检查系统是否稳定。 源程序： A=[2, -2.98, 0.17, 2.3418, -1.5147]; B=[0, 0, 1, 5, -50]; subplot(2,1,1); zplane(B,A); %求H(z)的极点 p=roots(A); %求H(z)的模 pm=abs(p); if max(pm)<1 disp('系统因果稳定'), else,disp('系统因果不稳定'),end un=ones(1,800); sn=filter(B, A, un); n=0:length(sn)-1; subplot(2, 1, 2);plot(n, sn) xlabel('n');ylabel('s(n)')
习题集-02 数字信号处理习题答案
§ Z 变换 ? Z 变换的定义及收敛域 【习题】 1. 假如)(n x 的z 变换代数表示式是下式，问)(z X 可能有多少不同的收敛域。 )83451)(411(411)(2122----+++- =z z z z z X 【分析】 ）要单独讨论，（环状、圆外、圆内：有三种收敛域：双边序列的收敛域为：特殊情况有：左边序列的收敛域为：因果序列的收敛域为：右边序列的收敛域为：特殊情况有：有限长序列的收敛域为 0 0 , , 0 0 , , 0 , 0 0 , 0 , 0 22 11 212 1∞==<<≤≤<≤<<≥≥∞≤<≥∞<<≤∞<≤≥∞≤<≤≤∞<<+ -++--z z R z R n n R z n n R z n n z R n n z R n z n z n n n z x x x x x x
解：对X (Z )的分子和分母进行因式分解得 )43 1 )(21 1)(211(2111111----+-+- =Z jZ jZ Z X (Z )的零点为：1/2，极点为：j/2，-j/2，-3/4 ∴ X (Z )的收敛域为： (1) 1/2 < | Z | < 3/4，为双边序列，见图一 (2) | Z | < 1/2，为左边序列，见图二 (3) | Z | > 3/4，为右边序列，见图三 图一 图二 图三 )431)(211)(411()211)(211()(11211-----++++- =Z Z Z Z Z Z X
? Z 反变换 【习题】 2. 有一右边序列 )(n x ，其 z 变换为)1)(211(1 )(11----=z z z X (a) 将上式作部分分式展开(用 1-z 表示)，由展开式求 )(n x 。 (b) 将上式表示成 z 的多项式之比，再作部分分式展开，由展开式求 )(n x ,并说明所得到的序列 与(a)所得的是一样的。 【注意】不管哪种表示法最后求出 x (n ) 应该是相同的。 解：(a) 因为11122 111)(---+--=z z z X 且x(n)是右边序列 所以 )()212()(n u n x n ?? ? ??-= (b) 122 1211 )1)(2 1(21231 )1)(2 1()(2 -+--+=---+=--=z z z z z z z z z X )()212( )1(2)1(21)()( n u n u n u n n x n n ??? ??-=-+-?? ? ??-=δ则
《数字信号处理》课程实验题目
计电学院《数字信号处理》课程实验 适用专业：电子通信工程专业；实验学时：9 学时 一、实验的性质、任务和基本要求 (一)本实验课的性质、任务 数字信号处理课程实验是数字信号处理课程的有效的补充部分，通过实验，使学生巩固和加深数字信号处理的理论知识的理解和掌握，在实验过程中了解简单但是完整的数字信号处理的工程实现方法和流程。通过实践进一步加强学生独立分析问题和解决问题的能力、实际动手能力、综合设计及创新能力的培养。 （二）基本要求 掌握数字信号处理基本理论知识和滤波器设计及应用。 （三）实验选项
二、实验教学内容 实验一 1、实验目的和要求 1）加深理解时域采样定理、体会使用MATLAB的离散FT函数fft( )来解决涉及模拟信号的问题； 2）加深理解对带通信号的采样特性，学会采用MATLAB解决该问题； 3）加深理解在频率采样法中，过渡点对所设计滤波器特性的影响。 2、实验要求 1）提供MATLAB程序，画出每个步骤的曲线图； 2）写实验报告，包含有对所得结果进行分析和说明。 第一组：张毅雷凌峰白法聪覃昱滔刘强何新文 第二组：邓志强林盛勇李日胜黎少锋梁聪杨晨 实验二 1、实验目的和要求 （1）加深理解采用数字信号处理方法对模拟信号处理的过程、掌握使用MATLAB处理的方法；对一段音乐信号进行处理和输出；要求画出滤波前后语音信号时域波形、信号和滤波器的幅度频率特性曲线、相位频率特性曲线； （2）加深对截断效应的理解； （3）掌握使用MATLAB设计滤波器，并对语音信号处理的方法。对一段音乐信号进行处理和输出；要求画出滤波前后语音信号时域波形、信号和滤波器的幅度频率特性曲线、相位频率特性曲线。 2、实验要求 1）提供MATLAB程序，画出每个步骤的曲线图； 2）写实验报告，包含有对所得结果进行分析和说明。 第九组:汪涛张汉毅巫金敏张经中柳泽举 第六组：罗涛梁乐杰黄乃生 实验三 1、实验目的和要求 掌握采用MATLAB数字滤波器设计软件编制方法。软件要求在界面内有不同类型（高通低通带通带阻）滤波器的选择、或者只对低通滤波器采用不同方法设
数字信号处理上机实验答案(全)1
第十章 上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一 系统响应及系统稳定性。 实验二 时域采样与频域采样。 实验三 用FFT 对信号作频谱分析。 实验四 IIR 数字滤波器设计及软件实现。 实验五 FIR 数字滤波器设计与软件实现 实验六 应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR 数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 10.1 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握 求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。 系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。 系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。 实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的[19]。系统的稳态输出是指当∞→n 时，系统的输出。如果系统稳定，信号加入系统后，系统输出的开始一段称为暂态效应，随n 的加大，幅度趋于稳定，达到稳态输出。 注意在以下实验中均假设系统的初始状态为零。 3．实验内容及步骤 （1）编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter 函数或conv 函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。 （2）给定一个低通滤波器的差分方程为