本科二年级电子信息工程专业《数字信号处理》课程设计教案

本科二年级电子信息工程专业《数字信号处理》课程设计教案

一、课程定位与目标

（一）课程定位

《数字信号处理》是电子信息工程专业本科二年级的核心专业必修课，是衔接《信号与系统》与后续《通信原理》、《数字图像处理》、《雷达信号处理》等方向的关键枢纽课程。本课程设计打破传统“理论讲授+验证实验”的割裂模式，以“复杂工程问题解决能力”为锚点，采用项目式学习范式，将离散傅里叶变换、快速算法、数字滤波器设计等核心知识模块统合于真实系统开发情境，旨在培养学生运用数学工具抽象物理世界、以算法实现系统功能的完整工程思维。

（二）课程目标

1.知识目标

【非常重要】【高频考点】系统阐释离散时间信号与系统的时域、频域及变换域分析框架，精准复述Z变换、DFT/FFT、IIR/FIR滤波器的数学原理与物理内涵；【重要】归纳DFT参数（采样率、点数、窗函数）对谱分析性能的制约规律；【一般】列举多速率信号处理的基本概念及工程应用场景。

2.能力目标

【非常重要】【热点】独立运用MATLAB/Python完成信号采集、算法仿真、性能评估全流程，具备针对具体指标（通带截止频率、阻带衰减、过渡带宽度）设计并实现数字滤波器的工程能力；【重要】能够对FFT运算量、滤波器阶数等资源约束进行定量估算与优化；【一般】通过小组协作完成一套实时频谱分析仪原型系统的开发，并以规范学术格式撰写技术报告。

3.素养目标

【重要】养成量化思维与误差意识，对仿真与理论结果的偏差能开展归因分析；【一般】理解信号处理技术在国防、通信、医疗等国家战略领域的支撑作用，增强科技报国的使命感；【一般】在项目迭代中锤炼抗挫折能力与团队沟通素养。

二、学情分析

（一）知识基础

学生已完成《高等数学》（级数、积分变换）、《复变函数》（留数定理）、《信号与系统》（连续/离散时间卷积、傅里叶变换）的学习，具备复频域分析的初步概念。但【难点】对离散域特有的“周期延拓”“栅栏效应”“有限字长效应”等非理想效应缺乏直观认知；【重要】混淆DTFT与DFT的适用场景，误将DFT视为连续傅里叶变换的简单数字化。

（二）能力水平

通过《高级语言程序设计》课程，学生具备C/Python基础语法能力，但【非常重要】将算法描述转化为可执行代码的工程实现经验严重不足，尤其对复数运算、矩阵索引、内存管理存在畏难情绪；【重要】MATLAB工具停留在“函数调用”层面，极少自主编写底层算法。

（三）学习风格

【一般】习惯于跟随教师板书进行纸笔推导，对开放性、非良构问题（如“设计一个滤波器，使…更好”）存在认知焦虑；【重要】对可视化、交互式仿真有较高兴趣，短视频平台的科普类信号处理实验（如“用声音烧毁蜡烛”）能迅速激发其探究欲望。

三、教材与参考资源

（一）主教材

【非常重要】奥本海姆（A.V.Oppenheim）《离散时间信号处理》（第3版），电子工业出版社。该书理论体系严谨，数学推导严密，为国内外顶尖高校通用经典教材，重点研读第3章（Z变换）、第8章（DFT）、第11章（FIR滤波）。

（二）实验指导书

【重要】自编《数字信号处理MATLAB/Python工程实践》，融合近五年全国大学生电子设计竞赛信号处理类真题及本院优秀课程设计案例，提供基础框架代码与调试模板。

（三）拓展资源

【一般】MITOpenCourseWare6.003课程视频（Oppenheim亲授）；IEEEXplore数据库“数字滤波器设计”综述性文献；TI公司《TMS320C6748DSP技术参考手册》；国家高等教育智慧教育平台《工程信号处理》虚拟仿真实验。

四、教学内容体系

（一）模块一：离散时间信号与系统基础（第1-2周）

1.离散时间信号

（1）序列表示与分类【重要】【基础】

（2）典型序列【一般】

[1]单位脉冲序列——时域表征基石，系统描述语言【非常重要】【高频考点】

[2]单位阶跃序列与矩形序列【一般】

[3]实指数序列收敛条件【重要】【易错】

[4]正弦序列周期性判定：2π/ω0为有理数【难点】【高频考点】

2.离散时间系统

（1）线性、时不变性、因果性、稳定性的判定方法【非常重要】【高频考点】

（2）差分方程迭代解法与系统框图（直接Ⅰ型、直接Ⅱ型）【重要】【一般】

（3）单位脉冲响应h(n)与卷积和：有限长/无限长序列卷积【非常重要】【高频考点】

3.Z变换

（1）Z变换定义：双边/单边，收敛域ROC与极点分布【非常重要】【高频考点】

（2）Z反变换：留数辅助定理、部分分式展开法、长除法适用场景【重要】【难点】

（3）Z变换性质：线性、时移、尺度变换、时域卷积、初值终值定理【非常重要】【高频考点】

（4）系统函数H(z)与频率响应H(e^{jω})：零极点图与幅频特性【非常重要】【高频考点】

（二）模块二：离散傅里叶变换与快速算法（第3-5周）

1.傅里叶变换的四种形式【重要】【基础】

2.离散傅里叶级数DFS与周期序列【一般】

3.离散傅里叶变换DFT

（1）DFT定义式：X(k)=Σx(n)W_N^{kn}，物理意义——频域采样【非常重要】【核心】

（2）DFT与DTFT、Z变换的关系：频域采样定理，内插恢复【非常重要】【高频考点】

（3）DFT性质：循环移位、循环卷积、共轭对称性、帕塞瓦尔定理【非常重要】【高频考点】

（4）频域采样造成时域周期延拓——栅栏效应、混叠与泄漏初探【重要】【难点】

4.快速傅里叶变换FFT

（1）基2时间抽取DIT-FFT算法：蝶形运算、旋转因子分布、原位计算【非常重要】【热点】

（2）基2频率抽取DIF-FFT算法：与DIT的转置关系【重要】【一般】

（3）FFT运算量分析：复数乘法次数N/2·log₂N，与直接DFT对比【重要】【工程应用】

（4）IFFT算法：通过FFT实现，缩放因子1/N【一般】

5.DFT/FFT参数选择

（1）采样率Fs与奈奎斯特频率【非常重要】【基础】

（2）频率分辨率Δf=Fs/N【非常重要】【高频考点】

（3）补零——提高视在分辨率，不改变物理分辨率【重要】【难点澄清】

（4）窗函数谱分析：主瓣宽度、旁瓣衰减、幅值恢复系数【重要】【热点】

（三）模块三：数字滤波器设计（第6-9周）

1.滤波器基本概念

（1）理想滤波器与非因果不可实现性【一般】

（2）技术指标：通带截止频率ωp、阻带截止频率ωs、通带波纹αp、阻带衰减αs【非常重要】【高频考点】

2.无限长单位脉冲响应IIR滤波器设计

（1）模拟滤波器原型：巴特沃斯（最大平坦）、切比雪夫Ⅰ/Ⅱ（等波纹）、椭圆（陡峭过渡带）【重要】【一般】

（2）脉冲响应不变法：时域采样，频域混叠【重要】【难点】

（3）双线性变换法：S平面到Z平面的非线性映射，预畸变校正【非常重要】【高频考点】【工程首选】

（4）IIR滤波器阶数估算公式及MATLAB（butter,cheby1,ellip）实现【重要】【热点】

3.有限长单位脉冲响应FIR滤波器设计

（1）线性相位条件：h(n)偶/奇对称，四种组合【非常重要】【高频考点】

（2）窗函数设计法：矩形、汉宁、汉明、布莱克曼、凯泽窗【非常重要】【工程常用】

（3）频率采样设计法：过渡带样本优化【重要】【一般】

（4）等波纹最佳逼近设计法（Parks-McClellan算法）【一般】【拓展】

（5）IIR与FIR比较：相位线性性、稳定性、阶数、实时性【重要】【高频考点】

4.滤波器实现结构

（1）IIR直接Ⅰ型、直接Ⅱ型（典范型）、级联型、并联型【重要】【一般】

（2）FIR直接型、级联型、频率采样型、快速卷积型【重要】【一般】

（3）量化效应：系数量化、运算舍入、极限环振荡【难点】【拓展】

（四）模块四：多速率信号处理与工程应用专题（第10周）

1.抽取与内插【重要】【热点】

2.多相分解与多相滤波器【一般】

3.采样率转换滤波器设计：半带滤波器、CIC滤波器【重要】【工程】

4.应用案例：软件无线电数字变频、音频重采样、雷达脉冲压缩【一般】【拓展】

五、教学实施过程

本课程总计48学时，理论讲授30学时，项目实践18学时。此处完整展开“基于FFT的实时频谱分析仪设计”项目（贯穿第3-5教学周），以此作为将数学变换转化为工程直觉的关键载体，彻底呈现教学实施全流程细节。

（一）阶段一：认知冲突与问题具象化（第3周第1次课，2学时）

1.课堂导入（0.25学时）

【教师行为】利用便携式声学测量套件现场演示：播放一段录制于嘈杂地铁环境的语音，时域波形完全淹没在噪声中；随后开启一台基于STM32的简易频谱分析仪，实时显示频域谱线，并用手持滤波器快速滤除50Hz工频及其谐波，语音信噪比显著提升。

【教师发问】“时域上看起来一团乱麻，为何在频域‘轻轻一刀’就干净了？那个将混乱变成有序的数学变换究竟做了什么？”【非常重要】

【学生反应】全体视觉聚焦于频谱跳动，产生强烈认知好奇。

【思政融入】介绍我国高速铁路轮轨噪声监测系统，正是通过车载频谱分析实时识别轴承故障，保障春运亿万旅客平安——科技的温度在于默默守护。

2.从连续到离散：傅里叶之路的四次跨越（0.75学时）

【教师行为】动画演示：傅里叶（热传导）→连续傅里叶变换→离散时间傅里叶变换（DTFT，频域连续）→计算机困境（无法存储无限长、连续频点）→DFT的诞生：频域离散化。

【板演核心】在黑板左侧写出DTFT公式X(e^{jω})=Σx(n)e^{-jωn}，右侧写出DFT公式X(k)=Σx(n)e^{-j2πkn/N}，邀请一名学生上台用彩色粉笔圈出两者的唯一区别：ω→2πk/N。

【教师强调】DFT是对DTFT频域的等间隔采样，采样点数N即为DFT点数；这一采样行为在时域引入了隐含的周期延拓。【非常重要】【高频考点】

【小组讨论】发放三组卡片，每组包含一个有限长序列及其周期延拓图，要求匹配对应的DFT幅谱形状。学生在移动白板上拼接并解释理由。

3.隐含的周期魔咒——从DFT公式到物理直觉（0.5学时）

【教师行为】运行MATLAB实时脚本：x=[1,1,1,1];N=4;X=fft(x);x_recon=ifft(X);重建序列与原序列完全一致。修改N=6，对x补零后ifft重建，得到[1,1,1,1,0,0]。

【教师追问】“重建结果证明了什么？”学生迟疑后答：“DFT对，IDFT也对。”教师升华：“这说明DFT自身是完备的数学对，它并不认为时域只有4个点，而是认为时域是这4个点的无限重复——所以重建时末尾补零，也是重复的一部分。”

【难点突破】将周期延拓类比为滚筒印花机：有限长图案印在无限布匹上，DFT系数就是滚筒的雕刻参数。【重要】

【课堂瞬间】学生出现豁然开朗的表情。

4.课后锚定任务（0.5学时）

【必做】阅读教材8.1-8.3节，完成习题8.2（DFT计算）、8.5（DFT与DTFT关系）。

【选做】查阅资料，简述科赫（Koch）噪声整形的原理，为后续FFT做心理铺垫。【一般】

（二）阶段二：DFT性质——从数学表达到工程权衡（第3周第2次课，2学时）

1.翻转课堂：性质的可视化演绎（0.5学时）

【学生展示】随机抽取三组，每组4分钟PPT+2分钟质疑。

第一组：演示循环移位对8点DFT幅度谱的影响，发现当移位长度等于周期整数倍时谱不变，非整数倍时相位谱旋转但幅度谱不变。

第二组：对比循环卷积与线性卷积，用卷积定理证明两者关系，并用条形图展示N≥L1+L2-1时两者相等。

第三组：帕塞瓦尔定理验证——时域能量与频域能量完全一致，用音频dropout能量泄漏反向验证。

【教师精讲】针对第一组，补充“循环移位在时域对应频域乘旋转因子，幅度谱确实不变，但若移出主值区间发生序列截断，则另当别论。”该细节为历年高频考点盲区。【非常重要】

2.工程困境：无限长信号如何用DFT做实时卷积（0.8学时）

【情境引入】现场连接麦克风，学生说一句长句子“信号处理是美丽的学科”，计算机实时采集。教师提问：“计算机必须等你说完整句话才能开始滤波吗？那直播延时将无法忍受。”

【分组探究】全班分两大阵营，左区研究重叠保留法，右区研究重叠相加法。

教师下发半成品MATLAB代码，已搭建好主循环框架，核心分块、补零、拼接部分留空。

【学生活动】每组2-3人围坐，在代码空白处填写正确索引，运行并测试一段chirp信号，对比两种方法输出差异与运行时间。

【成果对比】邀请两组代表将算法流程图画于黑板两侧，全班发现重叠保留法需丢弃每块前若干点，而重叠相加法需交叠相加——两者本质均为将无限长切割为有限长，再用DFT加速循环卷积，最后拼回线性卷积。【非常重要】

3.学术伦理与工程精神（0.2学时）

【案例】讲述我国“魂芯二号A”DSP研发团队平均年龄32岁，在无任何源代码参考下，从FFT硬件架构开始重新设计，最终运算效率达到国际同类产品1.35倍。【一般】

4.微测验（0.5学时）

【测验1】已知x(n)是8点序列，其6点DFT与8点DFT有何区别？【高频考点】

【测验2】两个4点序列，线性卷积长度7，若用8点循环卷积替代，哪些点与线性卷积不同？【重要】

（三）阶段三：FFT——从数学构造到工程复杂度革命（第4周第1次课，2学时）

1.科学史叙事（0.2学时）

【讲述】1965年，IBM研究员Cooley与普林斯顿大学Tukey在《MathematicsofComputation》发表论文，将DFT计算复杂度从N²降至Nlog₂N。彼时正值美苏核试验监测关键期，快速谱分析使次声波传感器数据处理时间从周缩至时。

2.蝶形算法的发生学建构（1.0学时）

【教师行为】从N=2点DFT开始，手算F(0)=x(0)+x(1)，F(1)=x(0)-x(1)——这是最简单的蝶形。

【追问】N=4时，能否将4点序列按奇偶分成两个2点DFT？

【板演】全黑板推演N=4DIT-FFT流图：第一级旋转因子均为W_2^0=1；第二级出现W_4^0和W_4^2=-j。

【学生行为】每人发一张A4打印的N=8空白流图，用蓝笔标数据流向，红笔填旋转因子W_N^k，绿笔写输出位序。教师巡堂发现共性问题：W_8^2与W_4^1混淆。

【规律总结】旋转因子W_N^k=e^{-j2πk/N}，在各级之间满足递推关系；原位运算特性使得存储空间只需N个复数单元。【非常重要】【高频考点】

3.代码显形：20行Python还原FFT本质（0.6学时）

【教师行为】JupyterNotebook逐行编写：

deffft(x):

n=len(x)

ifn<=1:returnx

even=fft(x[0::2])

odd=fft(x[1::2])

T=[np.exp(-2jnp.pi

k/n)*odd[k]forkinrange(n//2)]

return[even[k]+T[k]forkinrange(n//2)]+[even[k]-T[k]forkinrange(n//2)]

运行后对比np.fft.fft，误差10^-15量级。

【学生惊叹】“原来课本上密密麻麻的流图，核心代码不过十行！”学生对递归思想与分治策略产生审美共鸣。【非常重要】

4.工程隐喻：FFT与数字示波器的帧率（0.2学时）

【讲解】某品牌高精度示波器FFT刷新率达50fps，秘诀之一是将点数固定为1024/2048等2的幂，使硬件流水线无分支预测惩罚。

（四）阶段四：项目式沉浸——实时频谱分析仪原型开发（第4周第2次课至第5周，6学时）

本阶段为课程设计能力转化核心引擎，采用“目标拆解—支架提供—协作建构—对抗验证”四阶模型。

1.任务发布与指标分级（0.5学时）

【项目任务】基于计算机声卡/手机麦克风，开发实时频谱分析仪。

【必选指标】（1）实时采集音频，显示频率-幅度谱线；（2）可调采样率8k/16k/44.1kHz；（3）可调FFT点数256/512/1024；（4）可加矩形、汉宁、汉明窗之一；（5）在谱峰处标注频率值。

【进阶指标】（1）瀑布图（三维频谱堆叠）；（2）峰值保持轨迹；（3）总谐波失真（THD）估算。

【技术约束】禁止直接调用matplotlib.mlab.specgram或scipy.signal.spectrogram；FFT核允许调用np.fft.fft，但加窗、分帧、频率映射必须手写。

2.知识脚手架搭建——参数选择的物理制约（1.0学时）

【教师微讲座】“频谱仪不是魔镜，它有分辨率极限。”

（1）频率分辨率Δf=Fs/N，要分辨10Hz与11Hz，需Δf≤1Hz。

（2）采样定理：Fs≥2f_max，否则混叠。

（3）窗函数：矩形窗主瓣最窄但旁瓣高（-13dB），汉宁窗主瓣略宽（-31dB），用于强窄带信号；平顶窗幅值精度高，用于校准。【非常重要】【高频考点】

【仿真对比】同一段含有50Hz、120Hz成分的信号，分别加矩形窗和汉宁窗，矩形窗泄漏谱线掩盖了120Hz附近的弱小信号，汉宁窗清晰分辨。

3.小组协作开发——分模块迭代（4.0学时）

【第1学时】数据采集与分帧

教师演示pyaudio或sounddevice回调函数写法，强调缓冲区大小与帧移的平衡（帧移=帧长/4为常用值）。学生分组调试，常见错误：数据类型为int16而非float32，导致FFT结果异常。

【第2学时】加窗与频域变换

学生完成hamming窗序列生成，逐帧加窗，调用np.fft.fft，取模，转换dB（20log10|X|），生成频率轴向量。巡堂发现共性问题：只对实部模值，未处理复数；dB公式漏加eps致log0。

【第3学时】GUI与动态刷新

利用matplotlib.animation或pyqtgraph，设置blit加速。教师展示对比：无blit时CPU占用80%，有blit时降至25%。

【第4学时】峰值检测与标注

利用scipy.signal.find_peaks或简单差分阈值法，标注前三大峰值频率。学生自主测试纯净正弦波，精度可达±0.1Hz。

4.中期对抗性评审（0.5学时）

【现场挑战】随机交换小组代码，运行并输入同一声源——手机播放A4音（440Hz）。若某组频谱主峰偏移超过2Hz，全体研讨归因。

【原因聚焦】某组采样率设置44.1k，点数1024，Δf≈43Hz，远大于1Hz，故无法精确测量440Hz。学生主动提议降低采样率至8k或增加点数至8192。教师引导：点数增加虽改善分辨率，但牺牲实时性；工程即权衡。

5.成果升华与跨媒介测试（第5周课后+课堂展示）

【测试1】用函数发生器播放扫频信号100-1000Hz，频谱仪实现“追频”动画。

【测试2】手机播放吉他调音App，所有小组频谱仪均在196Hz（G3）、246.9Hz（B3）、293.7Hz（D4）附近出现峰值。

【测试3】播放一段电子琴和弦，频谱仪清晰显示基频及三次、五次谐波。

【教师升华】“贝多芬失聪后如何作曲？他用想象‘看见’了频谱。今天，你们用代码让声音变得可见——傅里叶变换是人类的第七感。”【思政自然融入】

六、评价与考核设计

（一）形成性评价（40%）

1.随堂测验（10%）：3次闭卷，题型为DFT参数计算、滤波器指标换算、FFT流图填空。【非常重要】

2.实验报告（15%）：共4次，采用“数据-图表-代码-误因”四段式结构，凡结论与理论不符但能合理解释者不扣分，鼓励质疑精神。【重要】

3.课堂表现与互评（15%）：基于翻转课堂主讲质量、组内贡献度调查表、Git提交记录（如采用）。【一般】

（二）项目终结性评