02356数字信号处理
数字信号处理(吴镇扬)课后习题答案(比较详细的解答过程)chap6
x ( m) x1 (m) = 0
或
m = 0,± M ,±2M ,⋯
其它
∞
(6.2a)
(6.2b .2b) x1 (m) = x(m) p(m) = x(m) ∑ δ (m − Mi) (6.2b)
i =−∞
是一脉冲串序列， 式中 p(m) 是一脉冲串序列， 它在 M 的整数倍处的值 其余皆为零。 表示将采样率减少 为 1，其余皆为零。令 ↓M 表示将采样率减少 M 倍 的抽取， 6.1.1） 6.1.2 式的含意如图 6.1.1 （6.1.1 和 .2） 6.1. 的抽取， 6.1.1） （6.1.2） （ 所示， M=3。 所示，图中 M=3。
1 p( n ) = M 数展开。 数展开。
M −1 k =0
e j 2πnk / M 为周期序列 p(n) 的付里叶级 p(n)的付里叶级 ∑
所以
1 M −1 j (ω − 2πk ) / M ′(e ) = X ) （6.4） .4） ∑ X (e M k =0
jω
′(e jω ) ， X (e jω ) 分 别 是 x ′(n) 和 x (n) 的 式中 X DTFT。这样， DTFT。这样， X ′(e jω ) 是原信号频谱 X (e jω ) 先作 的移位叠加 位叠加， M 倍的扩展再在 ω 轴上每隔 2π / M 的移位叠加，
而 X 1 (e ) =
jω n = −∞
∞
∑ x ( n ) p ( n)e
− jωn
1 M −1 j 2πnk / M − jωn = ∑ [ x ( n) ]e ∑e n = −∞ M k =0 1 M −1 = X (e j (ω − 2πk / M ) ) (6.3b (6.3b) ∑ M k =0
数字信号处理(姚天任江太辉第三版)课后习题答案
第二章2.1 判断下列序列是否是周期序列。
若是，请确定它的最小周期。
（1）x(n)=Acos(685ππ+n ) （2）x(n)=)8(π-ne j（3）x(n)=Asin(343ππ+n ) 解 (1)对照正弦型序列的一般公式x(n)=Acos(ϕω+n )，得出=ω85π。
因此5162=ωπ是有理数，所以是周期序列。
最小周期等于N=)5(16516取k k =。
（2）对照复指数序列的一般公式x(n)=exp[ωσj +]n,得出81=ω。
因此πωπ162=是无理数，所以不是周期序列。
（3）对照正弦型序列的一般公式x(n)=Acos(ϕω+n )，又x(n)=Asin(343ππ+n )＝Acos(-2π343ππ-n )＝Acos(6143-n π)，得出=ω43π。
因此382=ωπ是有理数，所以是周期序列。
最小周期等于N=)3(838取k k =2.2在图2.2中，x(n)和h(n)分别是线性非移变系统的输入和单位取样响应。
计算并列的x(n)和h(n)的线性卷积以得到系统的输出y(n)，并画出y(n)的图形。
(a)1111(b)(c)111110 0-1-1-1-1-1-1-1-1222222 33333444………nnn nnnx(n)x(n)x(n)h(n)h(n)h(n)21u(n)u(n)u(n)a n ===22解 利用线性卷积公式y(n)=∑∞-∞=-k k n h k x )()(按照折叠、移位、相乘、相加、的作图方法，计算y(n)的每一个取样值。
(a) y(0)=x(O)h(0)=1y(l)=x(O)h(1)+x(1)h(O)=3y(n)=x(O)h(n)+x(1)h(n-1)+x(2)h(n-2)=4,n ≥2 (b) x(n)=2δ(n)-δ(n-1)h(n)=-δ(n)+2δ(n-1)+ δ(n-2)y(n)=-2δ(n)+5δ(n-1)= δ(n-3) (c) y(n)=∑∞-∞=--k kn k n u k u a)()(=∑∞-∞=-k kn a=aa n --+111u(n)2.3 计算线性线性卷积 (1) y(n)=u(n)*u(n) (2) y(n)=λnu(n)*u(n)解：(1) y(n)=∑∞-∞=-k k n u k u )()(=∑∞=-0)()(k k n u k u =(n+1),n ≥0即y(n)=(n+1)u(n) (2) y(n)=∑∞-∞=-k k k n u k u )()(λ=∑∞=-0)()(k kk n u k u λ=λλ--+111n ,n ≥0 即y(n)=λλ--+111n u(n)2.4 图P2.4所示的是单位取样响应分别为h 1(n)和h 2(n)的两个线性非移变系统的级联，已知x(n)=u(n), h 1(n)=δ(n)-δ(n-4), h 2(n)=a nu(n),|a|<1,求系统的输出y(n).解 ω(n)=x(n)*h 1(n) =∑∞-∞=k k u )([δ(n-k)-δ(n-k-4)]=u(n)-u(n-4)y(n)=ω(n)*h 2(n) =∑∞-∞=k kk u a )([u(n-k)-u(n-k-4)]=∑∞-=3n k ka,n ≥32.5 已知一个线性非移变系统的单位取样响应为h(n)=an-u(-n),0<a<1 用直接计算线性卷积的方法，求系统的单位阶跃响应。
自学考试_浙江省2014年10月高等教育自学考试数字信号处理试题(02356)
A.{3,3,3,3}
B.{1,3,3,3}
C.{0,3,3,3}
D.{1,2,3,4}
8.4 点 序 列 x(n)和 5 点 序 列 y(n)的 线 性 卷 积 长 度 为
A.10
B.8
C.5
D.4
9.
fs N
N DFT
X0
A.0
B.fs/N
C.2fs/N
D.(N-1)fs/N
10.滤 波 器 的 一 个 极 点 为z= -0.95,没 有 非 原 点 处 的 零 点 ,则 该 滤 波 器 的 类 型 是
02356# 数字信号处理试题 第1页 (共3页 )
6.已知线性移不变因 果 系 统 的 系 统 函 数 为 H (z)= (1-1.3z-1)2(1+0.7z-1),则 该 系 统 的 系 统函数收敛域为
A.|z|>1.3
B.|z|>0.7
C.0.7<|z|<1.3 D.|z|<0.7
7.已 知 4 点 序 列 x(n)= {1,2,3,4},其 圆 周 共 轭 对 称 分 量 为
{x(n),0≤n≤N-1
y(n)=
的 DFT。
0,N≤n≤2N-1
22.已知某离散线性移不变系统的单位冲激响应h(n)={3,2,1,1;n=0,1,2,3},若 系 统 输 入
序 列 为 x(n)= {1,2,1,1;n=0,1,2,3},
(1)试 计 算 系 统 的 输 出 序 列 y(n);
非选择题部分
注意事项: 用 黑 色 字 迹 的 签 字 笔 或 钢 笔 将 答 案 写 在 答 题 纸 上 ,不 能 答 在 试 题 ,每 空 2 分 ,共 20 分 )
《数字信号处理》研究生复试大纲
《数字信号处理》研究生复试大纲南京理工大学研究生复试大纲学院（系）：电光学院课程名称：数字信号处理执笔人：宋耀良修（制）订日期： 2003年3月一、课程的考试目的与基本要求本课程的考试目的主要在于考查学生掌握离散时间信号与系统的分析方法、离散傅里叶变换的概念以及离散傅里叶变换的快速算法、掌握数字滤波器的各种设计方法，以及了解有限字长效应对系统性能的影响等方面内容的情况；通过考查同时检验学生分析问题，解决实际问题的能力。
要求学生较全面深刻地掌握数字信号变换的基本理论及其在数字系统分析的应用，全面掌握数字系统的基本结构和系统设计和实现的基本方法。
二、考查内容以教学大纲为依据，具体内容包括:1 绪论数字信号和数字信号处理的基本概念2离散时间信号与系统2.1 离散时间信号与系统的频域表示2.1.1 系统的频率响应；2.1.2 系统频率响应的两个性质；2.1.3 系统频率响应与单位取样响应的关系；2.1.4 序列的频域表示法；2.1.5 输出序列与输入序列的傅氏变换间的关系。
2.2 傅里叶变换的对称性质2.2.1 序列的共轭对称和共轭反对称2.2.2 偶序列与奇序列2.2.3 傅里叶变换的共轭对称和共轭反对称2.2.4 傅里叶变换的对称性质3连续时间信号的采样3.1 周期采样3.2 采样的频域表示3.2.1 奈奎斯特采样定理(△)3.3 由样本重构带限信号（★）3.4 连续时间信号的离散时间处理3.4.1 线性时不变离散时间系统3.4.2 冲激响应不变3.5 离散时间信号的连续时间处理3.5.1 非整数延时3.5.2 滑动平均3.6 利用离散时间处理改变采样率（★）3.6.1 采样率按整数因子减小3.6.2 采样率按整数因子增加3.6.3 采样率按非整数因子变化3.7 实际应用考虑的问题3.7.1 消除混叠的预滤波3.7.2 模拟到数字(A/D)转换3.7.3 量化误差分析3.7.4 D/A转换3.7.5 A/D和D/A转换中抽取和内插的应用4 z 变换4.1 z变换的定义和性质4.2 z 变换的存在条件及系统的收敛域4.3 逆z变换5线性时不变 ( LTI ) 系统的变换分析5.1 LTI系统的频率响应5.1.1 理想频率选择性滤波器5.1.2 相位失真和延迟5.2 用线性常系数差分方程所表征的系统函数5.2.1 稳定性和因果性；5.2.2 逆系统（★）5.2.3 有理系统函数的冲激响应5.3 有理系统函数的频率响应5.3.1 单个零点或极点的频率响应5.3.2 多重零极点的例子5.4 幅度和相位特性之间的关系5.5 全通系统5.6 最小相位系统5.6.1 频率响应的补偿5.6.2 最小相位系统的性质5.7 广义线性相位的线性系统（△）5.7.1 线性相位系统5.7.2 广义线性相位5.7.3 因果广义线性相位系统5.7.4线性相位系统与最小相位系统的关系6 离散时间系统结构6.1 线性常系数差分方程的方框图表示6.2 线性常系数差分方程的信号流图表示6.3 IIR系统的基本结构6.3.1 直接型(△)6.3.2 级联型(△)6.3.3 并联型(△)6.3.4 IIR系统中的反馈6.4 转置形式6.5 FIR系统的基本结构6.5.1 直接型(△)6.5.2 级联型(△)6.5.3 线性相位FIR系统的结构(△)6.6 有限精度数值效应概述6.6.1 数的表示法6.6.2 在系统中实现的量化6.8 系数量化的影响6.7.1 IIR系统中系数量化的影响(△)6.7.2 FIR系统中系数量化的影响6.7.3 其它灵敏度分析方法6.8 数字滤波器中舍入噪声的影响6.8.1 直接型IIR结构分析6.8.2 IIR系统定点实现中的幅度加权6.8.3 级联和并联IIR结构分析6.8.4 直接型FIR结构分析6.8.5 离散时间系统的浮点实现6.9 IIR数字滤波器定点实现中的零输入极限环6.9.1 由于舍入和截尾引起的极限环6.9.2 由于溢出而出现的极限环7滤波器设计方法7.1 由连续时间波器设计离散时间IIR滤波器7.1.1 滤波器设计的冲激响应不变法（△）7.1.2 双线性变换法（△）7.1.3 双线性变换法设计举例7.2 低通IIR滤波器的频率变换法7.3 用窗函数法设计FIR滤波器7.3.1 常用窗函数的性质（△）7.3.2 广义线性相位的合并7.3.3 凯泽窗滤波器设计法（△）7.3.4 凯泽窗与其它窗之间的关系7.4 凯泽窗法设计FIR滤波器举例7.4.1 高通滤波器举例7.4.2 离散时间微分器7.6 FIR滤波器的最佳逼近（★）7.6.1 最佳I型低通滤波器7.6.2 最佳II型低通滤波器7.6.3 帕克斯-麦克莱伦算法7.6.4 最佳FIR滤波器的特性7.7 FIR等波纹逼近举例7.7.1 低通滤波器7.7.2 零阶保持器的补偿7.7.3 带通滤波器8离散傅里叶变换8.1 周期序列的表示：离散傅里叶级数8.2 离散傅里叶级数的性质（△）8.2.1 线性；8.2.2 序列的移位8.2.3 对偶性8.2.4 对称性8.2.5 周期卷积8.3 周期序列DFS表示的性质汇总8.4 周期信号的傅里叶变换8.5 对傅里叶变换采样8.5 有限长序列的傅里叶表示：离散傅里叶变换（△）8.7 离散傅里叶变换的性质（△）8.7.1 线性8.7.2 序列的循环移位（★）8.7.3 对偶性8.7.4对称性8.7.5 循环卷积（★）8.8 离散傅里叶变换的性质汇总8.9 用离散傅里叶变换实现线性卷积（△）8.9.1 两个有限长序列的线性卷积8.9.2 循环卷积作为带有混叠的线性卷积8.9.3 用DFT实现线性时不变系统9 离散傅里叶变换的计算9.1 离散傅里叶变换的高效计算9.2 戈泽尔算法9.3 按时间抽取的FFT算法9.3.1 算法原理(△)9.3.2 算法流图(△)9.3.3 算法特点(△)9.3.4 变形的算法流图(★)9.4 按频率抽取的FFT算法9.4.1 算法原理(△)9.4.2 算法流图(△)9.4.3 算法特点(△)9.4.4 变形的算法流图(★)9.5 FFT算法的实现9.5.1 标号排列9.5.2 系数9.5.3 实现FFT的算法流程图9.5.4 实序列FFT算法9.5 复合数N的FFT算法（★）9.5.1 库利-图基算法9.5.2 素因子算法9.6 用卷积实现DFT线性调频变换算法（△）9.7 离散傅里叶变换计算中有限寄存器长度的影响9.7.1 在FFT的直接法计算中量化的分析9.7.2 定点FFT算法中量化效应的分析9.7.3 浮点FFT算法中量化效应小结9.7.4 在FFT中系数量化的影响Alan V. Oppenheim & Ronald W. Schafer,三. 主要参考教材Discrete-Time Signal Processing,(2and ) Prince Hall. （教材）1997.注:课程复试大纲填写说明:重点△，难点★。
山东自学考试法律专业科目表
山东省高等教育自学考试2012年10月份考试专业科目一览表
类 别
时
10月20日
03709马克思主义基本原理概论
本 00199中外民俗 03709马克思主义基本原理概论 00798商务交流
本 03709马克思主义基本原理概论
00844日语阅读（二） 07128日语写作基础 00654新闻采访写作 00657新闻心理学
00068外国财政
00076国际金融
00101外经贸经营与管理
概论
主考院校 山东财经大学
烟台大学 山东省广播电视大学 中国海洋大学 山东财经大学
山东大学
山东大学
青岛大学 山东轻工业学院 山东财经大学
030301 行政管理
040101 学前教育 040103 小学教育
00147人力资源管理（一） 专 00163管理心理学
00386幼儿文学 专 00389学前教育科学研究
理论和“三个代表”重要思
想概论
04729大学语文
03707毛泽东思想、邓小平
理论和“三个代表”重要思想 05372国际物流导论
概论
05362物流英语
00055企业会计学
07522运营管理
03707毛泽东思想、邓小平
理论和“三个代表”重要思
想概论
04729大学语文
03707毛泽东思想、邓小平
理论和“三个代表”重要思想
9：00-11：30
00075证券投资与管理
2013年四月份自考课程考试安排(定稿)
环境微生物学 环境信息系统 环境管理学 建筑结构 02655 建筑施工（二） 工程招标与合同管 06393 土木工程概论 00139 西方经济学 理 06397 工程造价管理与案 02391 工程力学（二） 例分析 06217 人机工程学 02051 物理化学（二） 02539 化学基础 06706 细胞工程技术
03285 00053 00144 00183 00144 00158 00207
990091 市场营销
990092 企业财务管理
990097 建筑工程
02391 工程力学（二） 02439 结构力学（二） 04270 西方现代心理学派 05615 心理健康教育概论 05626 变态心理学（一） 别
附件二：
湖北省2013年4月份高等教育自学考试全日制助学班及本科衔接教育考试科目及时间安排表
本
4月19日（星期五） 课程 上午(9：00~11： 课程 下午(2：30~5： 代号 30) 代号 00) 课程 代号 03707
科
4月20日（星期六） 4月21日（星期日） 课程 课程 课程 上午(9：00~11：30) 下午(2：30~5：00) 上午(9：00~11：30) 代号 代号 代号 毛泽东思想、邓小平 03706 思想道德修养与法律基础 02197 概率论与数理统计 00012 （二） 03709 马克思主义基本原理概论 理论和“三个代表” 00015 04183 概率论与数理统计（经 00020 00009 政治经济学（财经类） 重要思想概论 管类） 00420 物理（工） 中国近现代史纲要 00022 管理系统中计算机应 02198 线性代数 00023 用 04729 04184 线性代数（经管类） 健康教育学 护理学基础 内科护理学（二） 美学 现代汉语 新闻学概论 传播学概论 市场信息学 互联网数据库 电子商务法概论 市政学 公共政策 宪法学 03000 营养学 03179 生物化学（三） 03200 预防医学（二） 00530 中国现代文学作品选 00661 中外新闻作品研究 00892 商务交流（二） 00898 互联网软件应用与开发 00908 网络营销与策划 00054 管理学原理 00315 当代中国政治制度 02864 微生物学与免疫学基础 02899 03201 护理学导论 03007 00040 00651 00889 00902 00913 00226 00923
数字信号处理第三版西科大课后答案第3和4章
采用数字信号处理技术对生物医学信号进行分析与处理，如心电图、 脑电图等信号的处理与识别。
04
重点难点总结与复习指导
第三章重点难点总结
离散时间信号与系统的时域分析
掌握离散时间信号的定义、性质及分类，理解离散时间系统的描述方式，掌握卷积和的计 算方法。
离散时间信号的频域分析
理解周期信号的傅里叶级数展开，掌握离散时间信号的傅里叶变换及其性质，了解频域采 样理论。
内部奇点的留数和。这种方法适用于X(z)在复平面上有奇点的情况。
系统函数H(z)求解方法
直接法
根据系统差分方程，直接写出系统函 数H(z)的表达式。这种方法简单直接， 但需要注意差分方程的初始条件和边 界条件。
间接法
先求出系统的单位冲激响应h(n)，然 后根据h(n)求出H(z)。这种方法需要 先确定系统的单位冲激响应，计算量 相对较大。
课后习题解答与技巧
熟练掌握z变换的定义和性质，能够灵活运用这些 性质进行信号处理和系统分析。
理解系统函数H(z)的物理意义，掌握其求解方法 ，并能够根据H(z)分析系统的稳定性和频率响应 特性。
掌握z反变换的计算方法，能够根据具体情况选择 合适的方法进行求解。
在解答课后习题时，注意审题和理解题意，明确 题目要求和已知条件，选择合适的公式和方法进 行求解。同时，注意计算过程和结果的准确性， 避免出现计算错误或遗漏重要步骤的情况。
时不变性质
系统对输入信号的响应不随时间推移而改变，即 输入信号延迟或提前一定时间后，输出信号也相 应延迟或提前相同的时间。
稳定性判定
系统对任意有界输入信号的响应也是有界的，即 输出信号的幅度不会无限制地增长。
课后习题解答与技巧
自学考试_浙江省2016年4月高等教育自学考试数字信号处理试题(02356)
-57
6
加窗后滤波器性能指标
过渡带宽度 阻带最小衰减
/×2π/N
/dB
0.9
-21
2.1
-25
3.1
-44
3.3
-53
5.5
-74
02356# 数字信号处理试题 第4页 (共4页 )
个蝶形。
A.112
B.224
C.448
D.896
02356# 数字信号处理试题 第1页 (共4页 )
6.下 列 关 于 双 线 性 变 换 法 设 计 数 字 滤 波 器 描 述 错 误 的 是
A.设 计 数 字 滤 波 器 时 存 在 相 位 失 真
B.从S 平面映射到Z 平面时,引入了中间平面S1
。
题 16 图
17.线 性 移 不 变 系 统 h(n)=δ(n)
(是 / 不 是 )因 果 系 统 ,
(是 / 不 是 )稳 定 系 统 。
18.已知线性移不变系统的冲激响应为h(n)=1,0≤n≤ N-1,则|H(ejω)|=
,Φ(ω)=
,群 时 延 为
。
19.利用 DFT 分 析 时 域 连 续 信 号 频 谱 时,可 能 出 现
在 每 小 题 列 出 的 四 个 备 选 项 中 只 有 一 个 是 符 合 题 目 要 求 的 ,请 将 其 选 出 并 将 “答 题 纸 ”的 相
应 代 码 涂 黑 。 错 涂 、多 涂 或 未 涂 均 无 分 。
1.离散时间序列x(n) =ej511π(n+1) 的周期是
A.11
B.22/5
11.y(n)=x(n)+u(n+1)是 移 不 变 系 统 。
12.X(z)和 X* (z)的收敛域相同。
数字信号处理实验五 谱分析 哈工程
实验五谱分析一，试验目的：1 研究不同类型的窗函数，研究一些不同的方法来测试窗的性能；2 专注于有关窄带信号的几个不同的情形。
二，实验原理信号是无限长的，而在进行信号处理时只能采用有限长信号，所以需要将窄带信号“截断”。
在信号处理中，“截断”被看成是用一个有限长的“窗口”看无限长的信倍号，或者从分析的角度是无限长的信号x(t)乘以有限长的窗函数w(t),由傅立叶变换性质可知：x(t)w(t)==1/2πX(jw)W(jw)如果x(t)是频带有限信号，而w(t)是频带无限函数，截断后的信号也必是频带无限信号，从而产生所谓的频谱泄露。
频谱泄露是不可避免的，但是尽量减小，因此设计了不同的窗函数满足不同的要求。
从能量的角度，频谱泄露也是能量泄露，因为加窗后，是原来的信号集中在宅频带内的能量分散到无限的频带范围。
三，实验内容1．用MATLAB编程绘制各种窗函数的形状（1）矩形窗程序N=20;w=boxcar(N);nn=0:N-1;plot(nn,w)（2）汉宁窗程序N=20;w=hanning(N);nn=0:N-1;plot(nn,w)（3）汉明窗程序N=20;w=hamming(N);nn=0:N-1;plot(nn,w)（4）巴特利特窗程序N=20;w=bartlett(N);nn=0:N-1;plot(nn,w)（5）布莱克曼窗程序N=20;w=blackman(N);nn=0:N-1;plot(nn,w)（6）Triang窗程序N=20;w=triang(N); nn=0:N-1;plot(nn,w)（7）Kaiser窗程序N=20;w=kaiser(N,10); nn=0:N-1;plot(nn,w)（8）切比雪夫窗程序N=20;w=chebwin(N,30);nn=0:N-1;plot(nn,w)2，用MATLAB编程绘制各种窗函数的形状及其幅度响应。
(1)矩形窗N=20;w=boxcar(N);[H,W]=dtft(w,512);subplot(111),plot(W/2/pi,20*log(abs(H)));grid,title('MAGNITUDE RESPONSE')xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|/dB')（2）汉宁窗N=20;w=hanning(N);[H,W]=dtft(w,512);subplot(111),plot(W/2/pi,20*log(abs(H)));grid,title('MAGNITUDE RESPONSE')xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|/dB')（3）汉明窗N=20;w=hamming(N);[H,W]=dtft(w,512);subplot(111),plot(W/2/pi,20*log(abs(H)));grid,title('MAGNITUDE RESPONSE')xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|/dB')（4）巴特利特窗N=20;w=bartlett(N);[H,W]=dtft(w,512);subplot(111),plot(W/2/pi,20*log(abs(H)));grid,title('MAGNITUDE RESPONSE')xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|/dB')（5）布莱克曼窗N=20;w=blackman(N);[H,W]=dtft(w,512);subplot(111),plot(W/2/pi,20*log(abs(H)));grid,title('MAGNITUDE RESPONSE')xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|/dB')（6）Triang窗N=20;w=triang(N);[H,W]=dtft(w,512);subplot(111),plot(W/2/pi,20*log(abs(H)));grid,title('MAGNITUDE RESPONSE')xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|/dB')（7）Kaiser窗N=20;w=kaiser(N,10);[H,W]=dtft(w,512);subplot(111),plot(W/2/pi,20*log(abs(H)));grid,title('MAGNITUDE RESPONSE')xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|/dB')（8）切比雪夫窗N=20;w=chebwin(N,30);[H,W]=dtft(w,512);subplot(111),plot(W/2/pi,20*log(abs(H)));grid,title('MAGNITUDE RESPONSE')xlabel('NORMALIZED FREQUENCY'),ylabel('|H(w)|/dB')3，绘制矩形窗的幅频响应，窗长度分别为：N=10，N=20，N=50，N=100。
《数字信号处理》第三版课后实验答案_西安电子科技大学出版社
程序清单及波形显示： clc;close all;clear all;%======内容1：调用filter 解差分方程，由系统对u(n)的响应判断稳定性======A=[1,-0.9];B=[0.05,0.05]; %系统差分方程系数向量B 和A x1n=[1 1 1 1 1 1 1 1 zeros(1,50)]; %产生信号x1(n)=R8(n) x2n=ones(1,128); %产生信号x2(n)=u(n)y1n=filter(B,A,x1n); %求系统对x1(n)的响应y1(n) n=0:length(y1n)-1;subplot(2,2,1);stem(n,y1n,'.'); title('(a) 系统对R8(n)的响应y1(n)');xlabel('n');ylabel('y1(n)');y2n=filter(B,A,x2n); %求系统对x2(n)的响应y2(n) n=0:length(y2n)-1;subplot(2,2,2);stem(n,y2n,'.'); title('(b) 系统对u(n)的响应y2(n)');xlabel('n');ylabel('y2(n)');hn=impz(B,A,58); %求系统单位脉冲响应h(n) n=0:length(hn)-1;subplot(2,2,3);y=hn;stem(n,hn,'.'); title('(c) 系统单位脉冲响应h(n)');xlabel('n');ylabel('h(n)');(a) 系统对R8(n)的响应y1(n)ny 1(n )(b) 系统对u(n)的响应y2(n)ny 2(n )(c) 系统单位脉冲响应h(n)nh (n )%===内容2：调用conv 函数计算卷积============================x1n=[1 1 1 1 1 1 1 1 ]; %产生信号x1(n)=R8(n) h1n=ones(1,10); h2n=[1 2.5 2.5 1 ];y21n=conv(h1n,x1n); y22n=conv(h2n,x1n); figure(2)n=0:length(h1n)-1;subplot(2,2,1);stem(n,h1n); title('(d) 系统单位脉冲响应h1n');xlabel('n');ylabel('h1(n)'); n=0:length(y21n)-1;subplot(2,2,2);stem(n,y21n); title('(e) h1(n)与R8(n)的卷积y21n');xlabel('n');ylabel('y21(n)'); n=0:length(h2n)-1;subplot(2,2,3);stem(n,h2n); title('(f) 系统单位脉冲响应h2n');xlabel('n');ylabel('h2(n)'); n=0:length(y22n)-1;subplot(2,2,4);stem(n,y22n); title('(g) h2(n)与R8(n)的卷积y22n');xlabel('n');ylabel('y22(n)');nh 1(n )ny 21(n )(f) 系统单位脉冲响应h2nnh 2(n)(g) h2(n)与R8(n)的卷积y22nny 22(n )%=========内容3：谐振器分析======================== un=ones(1,256); %产生信号u(n) n=0:255;xsin=sin(0.014*n)+sin(0.4*n); %产生正弦信号A=[1,-1.8237,0.9801];B=[1/100.49,0,-1/100.49]; %系统差分方程系数向量B 和Ay31n=filter(B,A,un); %谐振器对u(n)的响应y31(n) y32n=filter(B,A,xsin); %谐振器对u(n)的响应y31(n) figure(3)n=0:length(y31n)-1;subplot(2,1,1);stem(n,y31n,'.'); title('(h) 谐振器对u(n)的响应y31n');xlabel('n');ylabel('y31(n)'); n=0:length(y32n)-1;subplot(2,1,2);stem(n,y32n,'.'); title('(i) 谐振器对正弦信号的响应y32n');xlabel('n');ylabel('y32(n)');050100150200250300(h) 谐振器对u(n)的响应y31nny 31(n )(i) 谐振器对正弦信号的响应y32nny 32(n )程序清单及波形显示：% DTMF 双频拨号信号产生6位电话号码 %clear all;clc;tm=[1,2,3,65;4,5,6,66;7,8,9,67;42,0,35,68]; % DTMF 信号代表的16个数N=205;K=[18,20,22,24,31,34,38,42];f1=[697,770,852,941]; % 行频率向量 f2=[1209,1336,1477,1633]; % 列频率向量 TN=input('键入6位电话号码= '); % 输入6位数字TNr=0; %接收端电话号码初值为零for l=1:6;d=fix(TN/10^(6-l)) TN=TN-d*10^(6-l); for p=1:4; for q=1:4;if tm(p,q)==abs(d); break,end % 检测码相符的列号q endif tm(p,q)==abs(d); break,end % 检测码相符的行号p endn=0:1023; % 为了发声，加长序列x = sin(2*pi*n*f1(p)/8000) + sin(2*pi*n*f2(q)/8000);% 构成双频信号sound(x,8000); % 发出声音 pause(0.1)% 接收检测端的程序X=goertzel(x(1:205),K+1); % 用Goertzel 算法计算八点DFT样本val = abs(X); % 列出八点DFT向量subplot(3,2,l);stem(K,val,'.');grid;xlabel('k');ylabel('|X(k)|') % 画出DFT(k)幅度axis([10 50 0 120])limit = 80; %for s=5:8;if val(s) > limit, break, end % 查找列号endfor r=1:4;if val(r) > limit, break, end % 查找行号endTNr=TNr+tm(r,s-4)*10^(6-l);enddisp('接收端检测到的号码为：') % 显示接收到的字符disp(TNr)显示结果：键入6位电话号码= 123456d = 1d = 2d = 3d = 4d = 5d = 6接收端检测到的号码为：123456050100|X (k )||X (k )|050100k|X (k )|050100k|X (k )|050100k|X (k )|050100k|X (k )|% DTMF 双频拨号信号产生8位电话号码 %clear all;clc;tm=[1,2,3,65;4,5,6,66;7,8,9,67;42,0,35,68]; % DTMF 信号代表的16个数 N=205;K=[18,20,22,24,31,34,38,42];f1=[697,770,852,941]; % 行频率向量 f2=[1209,1336,1477,1633]; % 列频率向量 TN=input('键入8位电话号码= '); % 输入8位数字TNr=0; %接收端电话号码初值为零for l=1:8;d=fix(TN/10^(8-l)) TN=TN-d*10^(8-l); for p=1:4; for q=1:4;if tm(p,q)==abs(d); break,end % 检测码相符的列号q endif tm(p,q)==abs(d); break,end % 检测码相符的行号p endn=0:1023; % 为了发声，加长序列 x = sin(2*pi*n*f1(p)/8000) + sin(2*pi*n*f2(q)/8000);% 构成双频信号 sound(x,8000); % 发出声音 pause(0.1)% 接收检测端的程序X=goertzel(x(1:205),K+1); % 用Goertzel 算法计算八点DFT 样本val = abs(X); % 列出八点DFT 向量 subplot(4,2,l);stem(K,val,'.');grid;xlabel('k');ylabel('|X(k)|') % 画出DFT(k)幅度axis([10 50 0 120])limit = 80; % for s=5:8;if val(s) > limit, break, end % 查找列号 endfor r=1:4;if val(r) > limit, break, end % 查找行号 endTNr=TNr+tm(r,s-4)*10^(8-l); enddisp('接收端检测到的号码为：') % 显示接收到的字符 disp(TNr) 显示结果：键入8位电话号码= 12345678 d = 1 d = 2 d = 3 d = 4 d = 5 d = 6 d = 7 d = 8接收端检测到的号码为：12345678|X (k )||X (k )||X (k )||X (k )||X (k )||X (k )|k|X (k )|k|X (k )|程序清单及波形显示： % 时域采样理论验证程序Tp=64/1000; %观察时间Tp=64微秒Fs=1000;T=1/Fs; M=Tp*Fs;n=0:M-1;t=n*T;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xat=A*exp(-alph*t).*sin(omega*t);Xk=T*fft(xat,M); %M点FFT[xat)]subplot(3,2,1); stem(n,xat,'.'); xlabel('n');ylabel('x1(n)'); title('(a) Fs=1000Hz');k=0:M-1;fk=k/Tp;subplot(3,2,2);plot(fk,abs(Xk));title('(a) T*FT[xa(nT)],Fs=1000Hz');xlabel('\omega/hz');ylabel('(H1(ejw))');axis([0,Fs,0,1.2*max(abs(Xk))]);Fs=300;T=1/Fs; M=Tp*Fs;n=0:M-1;t=n*T;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xat=A*exp(-alph*t).*sin(omega*t);Xk=T*fft(xat,M); %M点FFT[xat)]subplot(3,2,3); stem(n,xat,'.'); xlabel('n');ylabel('x2(n)'); title('(b)Fs=300Hz');k=0:M-1;fk=k/Tp;subplot(3,2,4);plot(fk,abs(Xk));title('(a) T*FT[xa(nT)],Fs=300Hz');xlabel('\omega/hz');ylabel('(H2(ejw))');axis([0,Fs,0,1.2*max(abs(Xk))]);Fs=200;T=1/Fs; M=Tp*Fs;n=0:M-1;t=n*T;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xat=A*exp(-alph*t).*sin(omega*t);Xk=T*fft(xat,M); %M点FFT[xat)]subplot(3,2,5); stem(n,xat,'.'); xlabel('n');ylabel('x3(n)'); title('(c) Fs=200Hz');k=0:M-1;fk=k/Tp;subplot(3,2,6);plot(fk,abs(Xk));title('(a) T*FT[xa(nT)],Fs=200Hz');xlabel('\omega/hz');ylabel('(H3(ejw))');axis([0,Fs,0,1.2*max(abs(Xk))])nx 1(n )(a) Fs=1000Hz(a) T*FT[xa(nT)],Fs=1000Hzω/hz(H 1(e j w ))nx 2(n )(b) Fs=300Hz(a) T*FT[xa(nT)],Fs=300Hzω/hz(H 2(e j w ))nx 3(n )(c) Fs=200Hz(a) T*FT[xa(nT)],Fs=200Hzω/hz(H 3(e j w ))%频域采样理论验证程序 clc;clear;close all; M=27;N=32;n=0:M; xa=0:(M/2); xb= ceil(M/2)-1:-1:0; xn=[xa,xb]; %产生M 长三角波序列x(n) Xk=fft(xn,1024); %1024点FFT[x(n)], 用于近似序列x(n)的TF X32k=fft(xn,32) ;%32点FFT[x(n)] x32n=ifft(X32k); %32点IFFT[X32(k)]得到x32(n) X16k=X32k(1:2:N); %隔点抽取X32k 得到X16(K) x16n=ifft(X16k,N/2); %16点IFFT[X16(k)]得到x16(n) subplot(3,2,2);stem(n,xn,'.'); title('(b) 三角波序列x(n)');xlabel('n');ylabel('x(n)');axis([0,32,0,20]) k=0:1023;wk=2*k/1024; subplot(3,2,1);plot(wk,abs(Xk));title('(a)FT[x(n)]'); xlabel('\omega/\pi');ylabel('|X(e^j^\omega)|');axis([0,1,0,200]) k=0:N/2-1; subplot(3,2,3);stem(k,abs(X16k),'.'); title('(c) 16点频域采样');xlabel('k');ylabel('|X_1_6(k)|');axis([0,8,0,200]) n1=0:N/2-1; subplot(3,2,4);stem(n1,x16n,'.'); title('(d) 16点IDFT[X_1_6(k)]');xlabel('n');ylabel('x_1_6(n)');axis([0,32,0,20]); k=0:N-1; subplot(3,2,5);stem(k,abs(X32k),'.'); title('(e) 32点频域采样');xlabel('k');ylabel('|X_3_2(k)|');axis([0,16,0,200]);n1=0:N-1; subplot(3,2,6);stem(n1,x32n,'.');box on title('(f) 32点IDFT[X_3_2(k)]');xlabel('n');ylabel('x_3_2(n)');axis([0,32,0,20])(b) 三角波序列x(n)nx (n )0100200(a)FT[x(n)]ω/π|X (e j ω)|(c) 16点频域采样k|X 16(k)|102030(d) 16点IDFT[X 16(k)]nx 16(n)(e) 32点频域采样k|X 32(k )|(f) 32点IDFT[X 32(k)]nx 32(n )程序清单及波形显示：% 用FFT 对信号作频谱分析 clear all;close all %实验内容(1)=================================================== x1n=[ones(1,4)]; %产生序列向量x1(n)=R4(n) M=8;xa=1:(M/2); xb=(M/2):-1:1; x2n=[xa,xb]; %产生长度为8的三角波序列x2(n) x3n=[xb,xa]; X1k8=fft(x1n,8); %计算x1n 的8点DFT X1k16=fft(x1n,16); %计算x1n 的16点DFT X2k8=fft(x2n,8); %计算x1n 的8点DFT X2k16=fft(x2n,16); %计算x1n 的16点DFT X3k8=fft(x3n,8); %计算x1n 的8点DFT X3k16=fft(x3n,16); %计算x1n 的16点DFT %以下绘制幅频特性曲线 subplot(1,2,1);stem(X1k8,'.'); %绘制8点DFT 的幅频特性图 title('(1a) 8点DFT[x_1(n)]');xlabel('ω/π');ylabel('幅度'); subplot(1,2,2);stem(X1k16,'.'); %绘制16点DFT 的幅频特性图 title('(1b)16点DFT[x_1(n)]');xlabel('ω/π');ylabel('幅度');figure(2)subplot(2,2,1);stem(X2k8,'.'); %绘制8点DFT 的幅频特性图 title('(2a) 8点DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度');subplot(2,2,2);stem(X2k16,'.'); %绘制16点DFT 的幅频特性图 title('(2b)16点DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度'); subplot(2,2,3);stem(X3k8,'.'); %绘制8点DFT 的幅频特性图 title('(3a) 8点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');subplot(2,2,4);stem(X3k16,'.'); %绘制16点DFT 的幅频特性图 title('(3b)16点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');-2.5-2-1.5-1-0.500.511.522.5(1a) 8点DFT[x 1(n)]ω/π幅度-2.5-2-1.5-1-0.500.511.522.5(1b)16点DFT[x 1(n)]ω/π幅度2468-4-2024(2a) 8点DFT[x 2(n)]ω/π幅度5101520-20-1001020(2b)16点DFT[x 2(n)]ω/π幅度2468-4-2024(3a) 8点DFT[x 3(n)]ω/π幅度-10-50510(3b)16点DFT[x 3(n)]ω/π幅度%实验内容(2) 周期序列谱分析==================================N=8;n=0:N-1; %FFT 的变换区间N=8 x4n=cos(pi*n/4);x5n=cos(pi*n/4)+cos(pi*n/8);X4k8=fft(x4n); %计算x4n 的8点DFT X5k8=fft(x5n); %计算x5n 的8点DFT N=16;n=0:N-1; %FFT 的变换区间N=16 x4n=cos(pi*n/4);x5n=cos(pi*n/4)+cos(pi*n/8);X4k16=fft(x4n); %计算x4n 的16点DFT X5k16=fft(x5n); %计算x5n 的16点DFT figure(3)subplot(2,2,1);stem(X4k8,'.'); %绘制8点DFT 的幅频特性图 title('(4a) 8点DFT[x_4(n)]');xlabel('ω/π');ylabel('幅度');subplot(2,2,3);stem(X4k16,'.'); %绘制16点DFT 的幅频特性图 title('(4b)16点DFT[x_4(n)]');xlabel('ω/π');ylabel('幅度');subplot(2,2,2);stem(X5k8,'.'); %绘制8点DFT 的幅频特性图 title('(5a) 8点DFT[x_5(n)]');xlabel('ω/π');ylabel('幅度');subplot(2,2,4);stem(X5k16,'.'); %绘制16点DFT 的幅频特性图title('(5b)16点DFT[x_5(n)]');xlabel('ω/π');ylabel('幅度');2468-1012x 10-15(4a) 8点DFT[x 4(n)]ω/π幅度5101520-4-2024x 10-15(4b)16点DFT[x 4(n)]ω/π幅度2468-4-2024(5a) 8点DFT[x 5(n)]ω/π幅度5101520-4-2024x 10-15(5b)16点DFT[x 5(n)]ω/π幅度%实验内容(3) 模拟周期信号谱分析=============================== figure(4) Fs=64;T=1/Fs; N=16;n=0:N-1; %FFT 的变换区间N=16 x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); %对x6(t)16点采样 X6k16=fft(x6nT); %计算x6nT 的16点DFT X6k16=fftshift(X6k16); %将零频率移到频谱中心 Tp=N*T;F=1/Tp; %频率分辨率F k=-N/2:N/2-1;fk=k*F; %产生16点DFT 对应的采样点频率（以零频率为中心） subplot(3,1,1);stem(fk,abs(X6k16),'.');box on %绘制8点DFT 的幅频特性图 title('(6a) 16点|DFT[x_6(nT)]|');xlabel('f(Hz)');ylabel('幅度'); axis([-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k16))])N=32;n=0:N-1; %FFT的变换区间N=16x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); %对x6(t)32点采样X6k32=fft(x6nT); %计算x6nT的32点DFTX6k32=fftshift(X6k32); %将零频率移到频谱中心Tp=N*T;F=1/Tp; %频率分辨率Fk=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心）subplot(3,1,2);stem(fk,abs(X6k32),'.');box on %绘制8点DFT的幅频特性图title('(6b) 32点|DFT[x_6(nT)]|');xlabel('f(Hz)');ylabel('幅度');axis([-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k32))])N=64;n=0:N-1; %FFT的变换区间N=16x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); %对x6(t)64点采样X6k64=fft(x6nT); %计算x6nT的64点DFTX6k64=fftshift(X6k64); %将零频率移到频谱中心Tp=N*T;F=1/Tp; %频率分辨率Fk=-N/2:N/2-1;fk=k*F; %产生16点DFT对应的采样点频率（以零频率为中心）subplot(3,1,3);stem(fk,abs(X6k64),'.'); box on%绘制8点DFT的幅频特性图title('(6a) 64点|DFT[x_6(nT)]|');xlabel('f(Hz)');ylabel('幅度');axis([-N*F/2-1,N*F/2-1,0,1.2*max(abs(X6k64))])-30-20-1001020300510(6a) 16点|DFT[x 6(nT)]|f(Hz)幅度010(6b) 32点|DFT[x 6(nT)]|f(Hz)幅度020(6a) 64点|DFT[x 6(nT)]|f(Hz)幅度程序清单及波形显示： clc;clear all;close all;fc1=250; fm1=15; fc2=500; fm2=50; fc3=1000; fm3=100;N=800; Fs=10000;Ts=1/Fs; n=[0:N-1];t=n*Ts;x11=cos(2*pi*fc1*t); x12=cos(2*pi*fm1*t); x1=x11.*x12; subplot(3,1,1);plot(t,x11,'g');plot(t,x12,'r');plot(t,x1,'b'); legend('载波','调制波 ','已调 ');xlabel('t/s');ylabel('波形')x=cos(2*pi*fc1*t).*cos(2*pi*fm1*t)+cos(2*pi*fc2*t).*cos(2*pi*fm2*t)+cos(2*pi*fc3*t).*cos(2*pi*fm3*t); subplot(3,1,2);plot(t,x);X=fft(x)subplot(3,1,3)k=[0:(N-1)/2]stem(k*2/N,abs(X(k+1))/max(abs(X(k+1))),'.');axis([0,0.3,0,1]);xlabel('\omeg a/\pi');ylabel('幅度');wp=[0.04,0.06];ws=[0.03,0.07];rp=0.1;rs=60;[N1,wp]=ellipord(wp,ws,rp,rs);[B,A]=ellip(N1,rp,rs,wp);y1=filter(B,A,x);figuresubplot(3,1,1);plot(t,x);xlabel('t/s');title('3路混合信号波形')subplot(3,1,2);[H1,w1]=freqz(B,A,N);plot(w1/pi,20*log10(abs(H1)));axis([0,0.5,-80,1]);xlabel('\omega/\pi');ylabel('|H(e^j\omega)|');title('中心频率为250Hz的频率响应');subplot(3,1,3);plot(t,y1);xlabel('t/s');ylabel('y1(t)');title('中心频率为250H的滤波信号')wp=[0.08,0.12];ws=[0.07,0.13];rp=0.1;rs=60;[N1,wp]=ellipord(wp,ws,rp,rs);[B,A]=ellip(N1,rp,rs,wp);y1=filter(B,A,x);figuresubplot(3,1,1);plot(t,x);xlabel('t/s');title('3路混合信号波形');subplot(3,1,2)[H1,w1]=freqz(B,A,N);plot(w1/pi,20*log10(abs(H1)));axis([0,0.5,-90,2]);xlabel('\omega/\pi');ylabel('|H(e^j\omega)|');title('中心频率为500Hz的频率响应')subplot(3,1,3);plot(t,y1);xlabel('t/s');ylabel('y1(t)');title('中心频率为500H的滤波信号')wp=[0.17,0.23];ws=[0.16,0.24];rp=0.1;rs=60;[N1,wp]=ellipord(wp,ws,rp,rs);[B,A]=ellip(N1,rp,rs,wp);y1=filter(B,A,x);figuresubplot(3,1,1);plot(t,x);xlabel('t/s');title('3路混合信号波形');subplot(3,1,2);[H1,w1]=freqz(B,A,N);plot(w1/pi,20*log10(abs(H1)));axis([0,0.5,-100,10]);xlabel('\omega/\pi');ylabel('|H(e^j\omega)|');title('中心频率为1000Hz的频率响应')subplot(3,1,3);plot(t,y1);xlabel('t/s');ylabel('y1(t)');title('中心频率为1000H的滤波信号')00.010.020.030.040.050.060.070.08-101t/s波形00.010.020.030.040.050.060.070.08-505ω/π幅度-505t/s3路混合信号波形-80-60-40-200ω/π|H (e j ω)|中心频率为250Hz 的频率响应-202t/sy 1(t )中心频率为250H 的滤波信号0.010.020.030.040.050.060.070.08-505t/s3路混合信号波形0.050.10.150.20.250.30.350.40.450.5-80-60-40-200ω/π|H (e j ω)|中心频率为500Hz 的频率响应0.010.020.030.040.050.060.070.08-202t/sy 1(t )中心频率为500H 的滤波信号0.010.020.030.040.050.060.070.08-505t/s3路混合信号波形0.050.10.150.20.250.30.350.40.450.5-100-500ω/π|H (e j ω)|中心频率为1000Hz 的频率响应0.010.020.030.040.050.060.070.08-202t/sy 1(t )中心频率为1000H 的滤波信号程序清单及波形显示： clc;clear;clear allN=1000; Fs=1000;T=1/Fs;Tp=N*T; t=0:T:(N-1)*T;fc=Fs/10;f0=fc/10; %载波频率fc=Fs/10，单频调制信号频率为f0=Fc/10; mt=cos(2*pi*f0*t); %产生单频正弦波调制信号mt ，频率为f0 ct=cos(2*pi*fc*t); %产生载波正弦波信号ct ，频率为fc xt=mt.*ct; %相乘产生单频调制信号xt nt=2*rand(1,N)-1; %产生随机噪声ntfp=150; fs=200;Rp=0.1;As=70; % 滤波器指标fb=[fp,fs];m=[0,1]; % 计算remezord 函数所需参数f,m,devdev=[10^(-As/20),(10^(Rp/20)-1)/(10^(Rp/20)+1)];[n,fo,mo,W]=remezord(fb,m,dev,Fs); % 确定remez 函数所需参数hn=remez(n,fo,mo,W); % 调用remez 函数进行设计,用于滤除噪声nt 中的低频成分yt=filter(hn,1,10*nt); %滤除随机噪声中低频成分，生成高通噪声yt xt=xt+yt; %噪声加信号 fst=fft(xt,N);k=0:N-1;f=k/Tp;subplot(2,1,1);plot(t,xt);grid;xlabel('t/s');ylabel('x(t)');axis([0,Tp/5,min(xt),max(xt)]);title('(a) 信号加噪声波形')subplot(2,1,2);plot(f,abs(fst)/max(abs(fst)));grid;title('(b) 信号加噪声的频谱')axis([0,Fs/2,0,1.2]);xlabel('f/Hz');ylabel('幅度')-10-50510t/sx (t )(a) 信号加噪声波形0501001502002503003504004505000.51(b) 信号加噪声的频谱f/Hz幅度%==调用xtg 产生信号xt, xt 长度N=1000,并显示xt 及其频谱,========= N=1000;xt=xtg(N);fp=120; fs=150;Rp=0.2;As=60;Fs=1000; % 输入给定指标% (1) 用窗函数法设计滤波器wc=(fp+fs)/Fs; %理想低通滤波器截止频率(关于pi归一化）B=2*pi*(fs-fp)/Fs; %过渡带宽度指标Nb=ceil(11*pi/B); %blackman窗的长度Nhn=fir1(Nb-1,wc,blackman(Nb));Hw=abs(fft(hn,1024)); % 求设计的滤波器频率特性ywt=fftfilt(hn,xt,N); %调用函数fftfilt对xt滤波subplot(4,1,1);plot(f,Hw);xlabel('f/Hz');ylabel('幅度'); title('(a)低通滤波器幅频特性');subplot(4,1,2);plot(t,ywt); title('(b)滤除噪声后的信号波形');xlabel('t/s');ylabel('ywt');% (2) 用等波纹最佳逼近法设计滤波器fb=[fp,fs];m=[1,0]; % 确定remezord函数所需参数f,m,dev dev=[(10^(Rp/20)-1)/(10^(Rp/20)+1),10^(-As/20)];[Ne,fo,mo,W]=remezord(fb,m,dev,Fs); % 确定remez函数所需参数hn=remez(Ne,fo,mo,W); % 调用remez函数进行设计Hw=abs(fft(hn,1024)); % 求设计的滤波器频率特性yet=fftfilt(hn,xt,N); % 调用函数fftfilt对xt滤波subplot(3,1,1);plot(f,Hw);xlabel('f/Hz');ylabel('幅度'); title('(c)低通滤波器幅频特性');subplot(3,1,2);plot(t,yet);title('(d)滤除噪声后的信号波形');xlabel('t/s');ylabel('yet');。
数字信号处理西安电子高西全课后答案
因果系统
因果系统是指系统的输出仅与输入的时间点有关，与输入的时间点无关。
信号与系统的关系
01
系统对信号的作用
系统对信号的作用可以改变信号 的幅度、频率和相位等基本属性 。
02
信号在系统中的传 播
信号在系统中传播时，会受到系 统的特性影响，从而改变信号的 基本属性。
03
系统对信号的响应
系统对信号的响应可以反映系统 的特性，从而可以用来分析和设 计系统。
02 离散傅里叶变换的定义
离散傅里叶变换是针对离散时间信号和系统的傅 里叶变换，它将离散时间信号分解成不同频率的 正弦波的叠加。
03 离散傅里叶变换的性质
离散傅里叶变换具有周期性、对称性和Parseval 等重要性质。
快速傅里叶变换算法
1 2 3
快速傅里叶变换算法的定义
快速傅里叶变换是一种高效计算离散傅里叶变换 的算法，它利用了循环卷积和分治的思想来降低 计算的复杂度。
03
数字信号处理技术能够提高通信系统的抗干扰性能、
传输效率和可靠性。
数字信号处理在通信中的应用
调制解调技术
调制是将低频信号转换为适 合传输的高频信号，解调是 将高频信号还原为原始的低
频信号。
通过调制解调技术，可以实 现信号的多路复用和高效传 输。
数字信号处理在通信中的应用
01
信道编码技术
02
信道编码是在发送端对信号进行编码，以增加信号的冗余 度，提高信号的抗干扰能力。
FIR数字滤波器的优 点
FIR数字滤波器具有稳定性好、易 于实现、没有递归运算等优点， 因此在一些需要稳定的系统中得 到广泛应用。
08
信号处理的应用
数字信号处理在通信中的应用
2024江苏省自学考试1月份考试科目资料
建筑工程
2080806
02446 建筑设备
28887 土木工程概论
06006 地基处理技术
30456 建筑工程事故分析
环境工程
2081102
28447 大气污染限制工程
28490 固体废弃物处理处置工程
28525 环境分析与监测
29804 工程制图（二）
28529 环境化学
27092 财务管理学
国际经济与贸易
2024243
271871 国际运输与保险
27332 当代中国经济运行
27186 制单结汇与报关实务
27183 国际经济学
会计
2024204
11002 公司法与企业法
27332 当代中国经济运行
27350 企业会计准则与制度
市场营销
2024208
11002 公司法与企业法
2024282
03615 选购绩效管理
27332 当代中国经济运行
03614 选购法务与合同管理
03618 选购项目管理
社会工作与管理
2030203
00281 社区社会工作
00279 团体社会工作
00280 西方社会学理论
00278 社会统计学
27051 社会保险学
27052 社会福利思想
学前教化
2040102
02189 机械制造基础
房屋建筑工程
1080801
00170 建筑工程定额与预算
02394 房屋建筑学
化工工艺
108120102481 物理化学（三）来自02175 分析化学（一）
27056 化工仪表及自动化
03146 化工原理（二）
数字信号处理第三版西科大课后答案第5章
第 4 章 时域离散系统的网络结构及数字信号处理的实现
［例4.4.2］ 假设系统函数如下式， 画出它的并联型 结构。
H (z) (2 0.379z 1)(4 1.24z 1 5.264z 2 ) (1 0.5z 1)(1 z 1 0.5z 2 )
解： 上式的分子分母是因式分解形式， 再写成下式:
［例4.4.3］ 为了保证滤波器的因果稳定性， 其系统 函数的极点必须保证全部集中在单位圆内。 如果有极点在 单位圆上， 则可以形成一个正弦波发生器。 利用这一原理 试设计正弦波发生器。
解： 假设有两个系统函数
H1(z)
Y1 ( z) X (z)
1
(sin0 )z 1 2(cos0 )z 1
z 2
4.4 例 题
［例4.4.1］ 设FIR滤波器的系统函数为 H (z) 1 (1 0.9z 1 2.1z 2 0.9z 3 z 4 ) 10
求出其单位脉冲响应， 判断是否具有线性相位， 画出直 接型结构和线性相位结构（如果存在）。
第 4 章 时域离散系统的网络结构及数字信号处理的实现 图4.4.1
N
H(z) C
Ak
k1 1 pk z 1
式中, pk是极点l， C是常整数， Ak是展开式中的系数。 一 般pk、 Ak都是复数。 为了用实数乘法， 将共轭成对的极点 放在一起， 形成一个二阶网络， 公式为
Hk
(z)
1
bk0 bk1z 1 ak1z 1 ak 2
z 2
第 4 章 时域离散系统的网络结构及数字信号处理的实现
第 4 章 时域离散系统的网络结构及数字信号处理的实现
该公式是频率采样结构的基本公式， 但它是一个不考虑 稳定性， 又可以应用复数乘法器的公式。为了稳定， 且使 用实数乘法器，
综改试点专业考试安排
2014年10月自考综改课程安排表
1、黑色为与自考共卷课程
2、红色为江西科技学院单独申报的安排课程
3、紫色为南昌理工学院单独申报的安排课程
4、绿色为江西渝州职业学院单独申报的安排课程
5、蓝色为江西科技学院、南昌理工学院申报相同的安排课程
6、紫红色为江西科技学院、江西渝州职业学院申报相同的安排课程
7.橙色为江西科技学院和江西服装学院申报相同的安排课程
8、褐色为江西科技学院、江西服装学院、江西渝州职业学院申报相同的安排课程
9.天蓝色为江西科技学院、南昌理工学院、江西服装学院申报相同的安排课程
二0一四年十月民办综改统考课程表。
山东省高等教育自学考试2011年10月份考试专业科目一览表
00051管理系统中计算机应 用
山东经济学院
青岛大学
03709马克思主义基本原理 概论
00102世界市场行情
00098国际市场营销学
04183概率论与数理统计 （经管类）
04184线性代数（经管类）
本
03709马克思主义基本原 理概论
04127文化市场与营销
00853广告学（二）
04129人类文明史
02797家畜饲养管理学
青岛农业大学
00182公共关系学 03005护理教育导论 00263外国法制史 05680婚姻家庭法
03628社区精神卫生护理 山东大学
00169房地产法
山东大学 青岛大学
00228环境与资源保护法学
00229证据法学
00228环境与资源保护法学 山东师范大学 山东财政学院
00924婚姻家庭法原理与实 00925公证与基层法律服务
00281社区社会工作 00285中国福利思想 00182公共关系学 00316西方政治制度 00321中国文化概论 00399学前游戏论
00449教育管理原理
06057生理心理学
山东财政学院
00228环境与资源保护法学 青岛大学 山东经济学院
00278社会统计学
00287发展社会学
00322中国行政史 00323西方行政学说史
05565物业管理法规
本
03709马克思主义基本原理 概论
04183概率论与数理统计
（经管类）
05374物流企业财务管理
00147人力资源管理（一） 03325劳动关系学
本 03709马克思主义基本原理
概论
03323劳动经济学
本
03709马克思主义基本原理 概论
《数字信号处理》第三版高西全版课后习题答案详解
数字信号处理课后答案 高西全、丁美玉版1.2 教材第一章习题解答1. 用单位脉冲序列()n δ及其加权和表示题1图所示的序列。
解：()(4)2(2)(1)2()(1)2(2)4(3) 0.5(4)2(6)x n n n n n n n n n n δδδδδδδδδ=+++-+++-+-+-+-+-2. 给定信号：25,41()6,040,n n x n n +-≤≤-⎧⎪=≤≤⎨⎪⎩其它（1）画出()x n 序列的波形，标上各序列的值；（2）试用延迟单位脉冲序列及其加权和表示()x n 序列； （3）令1()2(2)x n x n =-，试画出1()x n 波形； （4）令2()2(2)x n x n =+，试画出2()x n 波形； （5）令3()2(2)x n x n =-，试画出3()x n 波形。
解：（1）x(n)的波形如题2解图（一）所示。
（2）()3(4)(3)(2)3(1)6() 6(1)6(2)6(3)6(4)x n n n n n n n n n n δδδδδδδδδ=-+-+++++++-+-+-+-（3）1()x n 的波形是x(n)的波形右移2位，在乘以2，画出图形如题2解图（二）所示。
（4）2()x n 的波形是x(n)的波形左移2位，在乘以2，画出图形如题2解图（三）所示。
（5）画3()x n 时，先画x(-n)的波形，然后再右移2位，3()x n 波形如题2解图（四）所示。
3. 判断下面的序列是否是周期的，若是周期的，确定其周期。
（1）3()cos()78x n A n ππ=-，A 是常数；（2）1()8()j n x n e π-=。
解：（1）3214,73w w ππ==，这是有理数，因此是周期序列，周期是T=14； （2）12,168w wππ==，这是无理数，因此是非周期序列。
5. 设系统分别用下面的差分方程描述，()x n 与()y n 分别表示系统输入和输出，判断系统是否是线性非时变的。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
高纲1197江苏省高等教育自学考试大纲02356数字信号处理南京理工大学编江苏省高等教育自学考试委员会办公室I、课程性质与设置目的要求数字信号处理课程是江苏省高等教育自学考试电子工程专业的必修课，是为培养和检验自学应考者具有数字信号处理系统的基本理论和实际应用的一门专业基础课。
数字信号处理是随着计算机技术的发展而迅速发展起来的一门新兴而古老的学科，它在新的领域如生物医学工程、声学、雷达、地震不、语音通信、数据通信、核科学等学科发挥着重要的作用，而它所采用的各种方法及众多应用已有悠久的历史；同时也是一门具有很强的理论性与实践性，且理论和技术发展都十分迅速的前沿性学科。
随着数字化时代的来临，科学技术的进步而生产发展需求的与日俱增，促进了数字信号处理学科的发展，产生了各种巧妙的信号处理算法；特别是计算机技术的飞速发展，为数字信号处理增添了巨大的生命力。
数字信号处理主要是研究用数字或符号序列表示和处理信号。
处理的目的可以是削弱信号中的多余内容，滤除混杂的噪声和干扰，或者是将信号变换为容易分析和识别的形式，便于估计和选择它的特征参数。
例如通过分析和运算，可以估计脑电图或心电图中的某种特征参数，帮助医生查找病因和分析病情，确定合理的治疗方案；又如，信号在传输时，要受到各种干扰，包括失真、衰落和混入的背景噪声，信号处理要排除这些干扰。
设置本课程具体的目的要求是：使自学应考者理解数字信号处理系统的基本组成和基本原理；掌握线性时不变离散时间系统的三种分析方法——卷积和方法、频域分析方法和Z变换分析法；离散时间信号的频谱——傅里叶变换；离散傅里叶变换（DFT）及其快速算法（FFT）；数字滤波器的基本结构；IIR和FIR数字滤波器设计的常用方法，并对新的信号处理方法——小波变换及常用的数字信号处理仿真软件Matlab有所了解。
通过系统而专门的学习，培养电子工程专业学生应用数字信号处理基本理论的基本方法处理实际信号的能力。
Ⅱ、考试目标（考核知识点、考核要求）绪论一、考试知识点（一）信号、系统和信号处理的基本概念（二）数字信号处理的基本组成及实现方法（三）数字信号处理系统的优点及其应用二、考核要求（一）信号、系统和信号处理的基本概念1、理解:(1)信号的基本概念及四类信号的划分；（2）系统的基本概念和四类系统的划分；（3）信号处理的本质和目的。
（二）数字信号处理的基本组成及实现方法1、理解：（1）数字信号处理系统的基本组成及各部分作用；（2）数字信号处理系统的实现方法。
（三）数字信号处理系统的优点及其应用1、理解：（1）数字信号处理系统的优点；（2）数字信号处理的应用领域。
第1章离散时间信号与系统一、考核知识点(一)离散时间信号——序列（二）连续时间信号的采样（三）离散时间系统时域分析（四）离散时间系统频域及复频域分析二、考核要求(一)离散时间信号——序列1、理解：（1）离散时间序列的图形表示；（2）序列能量的含义。
2、掌握：（1）常用的离散时间序列；（2）周期序列的概念。
3、应用：（1）序列的运算（移位，翻褶，加法，乘法，标乘）；（2）应用单位采样序列表示任意序列。
（二）连续时间信号的采样1、理解：（1）理想采样的概念；（2）带限信号重构的概念。
2、掌握：（1）理想采样信号的频谱特点；（2）采样信号恢复为连续时间信号的原理与方法。
3、应用：（1）奈奎斯特采样定理。
（2）离散时间系统处理连续时间信号。
（三）离散时间系统时域分析1、理解：（1）线性系统的概念；（2）时不变系统的概念；（3）常系数线性差分方程的概念。
2、掌握：（1）系统级联和并联的概念；（2）因果系统的概念；（3）稳定系统的概念。
3、应用：（1）单位脉冲响应，包括据此判断系统的因果性和稳定性；（2）线性卷积。
（四）离散时间系统频域及复频域分析1、理解：（1）z变换及其逆变换的概念；（2）序列的傅里叶变换的概念；（3）极点和零点；（4）IIR系统和FIR系统。
2、掌握：（1）z变换及其逆变换的求法；（2）z变换收敛域的确定；（3）z变换的基本性质；（4）拉氏变换和z变换的映射关系；（5）系统函数；（6）系统的频率响应。
3、应用：（1）应用部分分式展开法求解逆z变换；（2）z变换性质；（3）系统函数，包括：由收敛域判断系统因果性和稳定性，由零、极点分布画系统的幅频响应曲线，与单位脉冲响应的关系，与差分方程的关系等。
第2章离散傅里变换一、考核知识点（一）离散傅里叶级数（DFS）（二）离散傅里叶变换（DFT）（三）频域采样理论二、考核要求（一）离散傅里叶级数1、理解：（1）周期离散序列的离散傅里叶级数的意义；（2）主值序列的含义。
2、掌握：（1）四种形式信号的傅里叶变换形式；（2）DFS与IDFS；（3）DFS的性质。
（二）离散傅里叶变换1、理解：（1）取模值运算（又称取余运算）和周期延拓的概念。
2、掌握：（1）有限长序列的DFT与IDFT；（2）有限长序列的DFT与周期序列DFS的关系；（3）离散傅里叶变换与离散序列的傅里叶变换和Z变换的关系；（4）DFT的性质。
3、应用：（1）应用DFT的性质计算有限长序列的DFT；（2）应用基本的数学关系及DFS 的性质证明DFT的性质。
（三）频域采样理论1、理解：（1）频域采样含义；（2）频域采样的条件。
第3章快速傅里叶变换（FFT）一、考核知识点（一）FFT算法（二）线性调频Z变换算法（三）FFT应用二、考核要求（一）FFT算法1、理解：（1）FFT的意义；（2）基-2的含义；（3）蝶形运算的含义；（4）信号流图的概念；（5）同址运算的概念；（6）倒位序的含义。
2、掌握：（1）按时间抽取的基-2 FFT运算原理和运算特点；（2）按频率抽取的基-2FFT 运算原理和运算特点；（3）按时间抽取的基-2 FFT蝶形运算流图；（4）按频率抽取的基-2FFT 运算流图。
3、应用：（1）DFT及基-2FFT的复数运算量的计算；（2）应用FFT算法和运算流图确定各个参量，如倒位序，蝶形运算两节点的距离，权系数Wr N，以及程序框图等。
（二）线性调频Z变换算法1、理解：（1）线性调频Z变换的概念；（2）算法原理。
2、掌握：（1）线性调频Z变换的实现步骤。
（三）FFT应用1、理解：（1）连续时间信号离散频谱分析的基本概念（框图结构）及其所造成的误差；（2）频率分辨率的物理意义；（3）重叠相加法及重叠保留法的概念；（4）信号消噪的概念。
2、掌握：（1）连续时间信号离散频谱分析所需满足的条件；（2）快速卷积法及其实现框图。
3、应用：（1）应用FFT分析实际连续信号的频谱；（2）应用FFT分析实信呈的DFT；（3）应用FFT计算线性时不变离散时间系统的输出。
第4章数字滤波器的基本结构一、考核知识点（一）数字滤波器的基本概念（二）IIR数字滤波器结构（三）FIR数字滤波器结构二、考核要求（一）数字滤波器的基本概念1、理解：（1）数字滤波器的基本概念；（2）数字滤波器的基本运算单元；（3）数字滤波器的运算结构（又可用信号流图表示）。
（二）IIR数字滤波器结构1、理解：（1）IIR数字滤波器的含义；（2）递归结构、反馈支路、节点的概念。
2、掌握：（1）直接型（Ⅱ型）结构；（2）级联型结构；（3）并联型结构。
3、应用：（1）根据描述系统的差分方程、系统函数或单位脉冲冲应画出不同形式的IIR系统信号流图；（2）根据IIR系统信号流图写出描述系统的差分方程、系统函数或单位脉冲响应。
（三）FIR数字滤波器结构1、理解：（1）FIR数字滤波器的含义；（2）非递归结构的概念；（3）横截型结构的意义。
2、掌握：（1）直接型结构；（2）级联型结构：（3）频率采样型结构；（4）快速卷积型结构。
3、应用：（1）根据描述FIR系统的差分方程、系统函数或单位脉冲应画出不同形式的系统信号流图；（2）根据系统信号流图写出描述FIR系统的差分方程、系统函数或单位脉冲响应。
第5章无限长单位脉冲响应（IIR）数字滤波器的设计方法一、考核知识点（一）数字滤波器设计的基本概念（二）常用模拟低通滤波器的设计方法（三）用脉冲响应不变变换法设计IIR数字滤波器（四）用双线性变换法设计IIR数字滤波器（五）设计IIR滤波器的频率变换法（六）Z平面变换法二、考核要求（一）数字滤波器设计的基本概念1、理解：（1）选频滤波器的概念；（2）滤波器的设计指标；（3）阻带衰减和过渡带的概念；（4）IIR和FIR滤波器设计特点。
2、掌握：（1）高通滤波器；（2）低通滤波器；（3）带通滤波器；（4）带阻滤波器；（5）数字滤波器的设计步骤；（6）IIR和FIR滤波器。
3、应用：（1）根据系统函数零极点分布或幅度频率响应曲线判断系统的滤波器类型。
（二）常用模拟低通滤波器的设计方法1、理解：（1）模拟滤波器的基本概念；（2）巴特活沃思滤波器、切比雪夫滤波器的幅度频率响应特点；（3）3dB截止频率和3dB带宽；（4）模拟低通、高通、带通和带阻滤波器的概念。
2、掌握：（1）由幅度平方函数确定系统的原理和方法；（2）巴特沃思模拟低通滤波器的设计方法；（3）切比雪夫横拟通滤波器的设计方法。
3、应用：（1）应用幅度平方函数确定巴特沃思模拟低通滤波器。
（三）脉冲响应不变变换法设计IIR数字滤波器1、理解：（1）脉冲响应不变变换法的概念；（2）混叠失真的概念。
2、掌握：（1）脉冲响应不变变换法设计原理；（2）模拟波波器的数字化方法；（3）脉冲响应不变变换法设计IIR数字滤波器的设计步骤；（4）适用范围3、应用：（1）用脉冲响应不变变换法设计IIR数字滤波器。
（四）双线性变换法设计IIR数字滤波器1、理解：（1）双线性变换法的概念；（2）频率畸变的概念；（3）非线性失真的概念。
2掌握：（1）双线性变换法变换原理；（2）预畸变的方法；（3）模拟滤波器的数字化方法；（4）双线性变换法设计IIR数字滤波器的设计步骤；（5）适用范围。
3、应用：用双线性变换法设计IIR数字滤波器（五）设计IIR滤波器的频率变换法1、理解：（1）模拟原型的概念：（2）频率变换的方法及意义。
（六）Z平面变换法1、理解：（1）Z平面变换的概念2、掌握：（1）Z平面变换满足的基本原则；（2）数字低通-数字低通变换的基本特性。
第6章有限长单位脉冲响应（FIR）数字滤波器的设计方法一、考核知识点（一）线性相位FIR滤波器（二）窗函数设计FIR滤波器（三）频率采样法设计FIR滤波器（四）等波纹线性相位滤波器（五）数字滤波器的应用二、考核要求（一）线性相位FILR滤波器1、理解：（1）线性相位数字滤波器的概念；（2）群延时的概念。
2、掌握：（1）线性相位FIR数字滤波器的条件；（2）四种形式的线性相位FIR数字滤波器的幅小辈特性和相频特性特点；（3）线性相位FIR数字滤波器的零点、极点分布特点。