基于采样定理与混叠效应的动态演示系统

资 源 简 介

本项目旨在通过MATLAB编程实现采样定理的直观演示，涵盖时域连续信号采样以及频域频谱周期延拓的完整过程。系统能够生成特定频率的余弦信号或组合正弦信号作为原始信号，并根据用户设定的采样频率执行采样操作。程序将详细对比在满足奈奎斯特采样准则（即采样频率大于两倍信号最高频率）与不满足该准则（采样频率小于两倍最高频率）两种情况下的信号重建效果。在频域分析中，项目将展示信号频谱在采样后的周期性延拓现象，明确反映当采样频率过低时，相邻周期的频谱如何发生重叠并导致混叠效应。用户可以通过改变采样频率实时观察时域波形的变

详 情 说 明

基于MATLAB的采样定理与混叠效应动态演示系统

---

项目介绍

本项目是一个用于演示数学采样理论与数字信号处理核心概念的教学仿真系统。系统通过构造包含多个频率成分的复合余弦信号，完整模拟了信号从连续到离散的采样过程，并展示了信号如何通过理想插值算法进行无损重建。该系统重点展示了当采样频率高于和低于奈奎斯特频率（Nyquist Rate）时的不同表现，直观呈现了混叠效应在频域中的产生机理。

---

功能特性

1. 复合信号仿真：系统生成由10Hz和25Hz频率组合而成的原始信号，模拟真实世界中的多频率叠加场景。
2. 双模式对比实验：

* 满足采样准则模式：设置60Hz采样率（大于信号最高频率2倍的情况），演示信号的完美还原。 * 欠采样混叠模式：设置35Hz采样率（低于信号最高频率2倍的情况），演示信号失真与频率重叠现象。

1. 高精度频域分析：利用FFT技术展示信号频谱，并特别设计了周期延拓的可视化，展现主谱与旁谱之间的关系。
2. 理想重建算法：基于Sinc插值公式实现理想低通滤波提取，对比原始信号与重建信号的差异。
3. 量化评估指标：计算并输出均方误差（MSE）与最大偏差，为分析信号失真程度提供客观数据。

---

逻辑结构与实现方案

本演示系统主要由以下五个逻辑模块组成，各模块间严格遵循数字信号处理的时序逻辑：

1. 参数初始化与基准生成 逻辑中预设了10Hz和25Hz的测试频率，取其中最大值25Hz作为基准判定标准。程序首先以2000Hz的超高采样率渲染原始信号。这一高密度的时域轴被视为“连续模拟信号”，作为后续采样与重建对比的真值参考。

2. 离散采样模拟 系统针对两种不同的采样频率（60Hz和35Hz）遍历执行。根据采样频率计算采样间隔，从原始连续信号中提取对应的离散样本点。这一步骤完整模拟了模数转换（ADC）中的时间离散化过程。

3. 频域周期延拓可视化 这是本系统的核心演示逻辑。不仅计算当前信号的FFT，程序还手动模拟了频谱在频域上的周期性搬移。通过绘制主频谱及其左右相邻周期的副本，清晰展示了当采样频率降低时，相邻周期的频谱如何相互靠近并最终在重叠区产生混叠。

4. 基于Sinc函数的时域重建 系统并未直接使用内置滤波函数，而是采用了计算精度更高的Whittaker–Shannon插值公式。通过每个采样点与Sinc函数的乘积累加，模拟理想物理系统的无失真恢复逻辑。该过程能直观显示出当采样率不足时，重建信号如何偏离原始波形。

5. 误差量化分析与多维展示 系统实时计算重建信号与基准信号的偏差。在图形界面中，通过2x3的矩阵布局同时对比显示时域采样图、频域周期延拓图、以及时域重建对比图。当检测到欠采样状态时，界面会自动触发红色的“发生混叠”预警标志。

---

关键函数与算法说明

• Sinc插值算法：系统实现的核心算法为重建公式，它代表了时域中的卷积过程，等效于频域中的理想矩形低通滤波。这是证明采样定理能从离散点恢复连续包络的关键。
• 频谱搬移技术：通过fftshift对信号进行频谱中心化处理，并利用手动坐标偏移及填充绘图，展示了抽样后频谱以采样频率为周期进行循环延拓的数学事实。
• 计算指标：均方误差（MSE）的引入，使得用户不仅能从视觉上观察波形畸变，还能从数值上得到信号受混叠影响的损坏程度。

---

使用方法

1. 启动MATLAB开发环境。
2. 将系统关联的脚本文件保存至当前工作目录。
3. 在命令行窗口直接运行该程序主体。
4. 系统将自动弹出交互式图窗，左右对比展示不同采样率下的处理结果。
5. 观察控制台（Command Window）输出的性能指标，分析不同采样频率对信号还原精度的影响。

---

系统要求

• 运行环境：MATLAB R2016b 或更高版本。
• 所需工具箱：基础MATLAB组件（包含信号处理基本函数库）。
• 硬件建议：标准图形显示设备，建议分辨率不低于1280x800以获得最佳的多子图观察效果。