MATLAB中绘制冲击响应谱

介绍了利用快速傅里叶变换实现有限冲击响应低通数字滤波器的方法，并使用MATLAB语言进行程序设计。

这个程序从指定的属性表中提取正常冲击的性质。

在Matlab中，频率响应可以通过离散傅里叶变换在单位圆上的Z变换来表达，这反映了内插函数的频率特性。

MATLAB中的DFT源代码演示-功率谱和功率谱密度（PSD）这是展示如何实时计算功率谱（PS）和功率谱密度（PSD）的演示。我们直接使用DFT（或FFT）计算功率谱。使用MATLAB函数周期图来计算功率谱和功率谱密度。观察到：（1）功率谱等于DFT绝对值的平方。 （2）功率谱中所有功率谱线的总和等于输入信号的功率。 （3）PSD的积分等于输入信号的功率。通过周期图获得的PSD是等效噪声功率谱密度（ENPSD）。可以看到ENPSD与PS的相关系数为1 / T，等于频率分辨率或频率间隔。应当注意，功率谱是离散序列或离散连续自变量函数，而ENPSD是非离散连续自变量函数。为了强调这一点，我将茎用于功率谱，将图用于ENPSD。在本演示中，我们从具有各种参数的正弦信号开始。然后，我们继续处理实际的音频信号。

冲击滤波器的核心思想是在局部区域应用膨胀或侵蚀过程，取决于像素是否属于最大值或最小值的影响区域。基于Kramer-Bruckner的拉普拉斯算子集合{-1, 0, +1}，使用符号函数来决定膨胀和侵蚀的优先顺序。这种方法产生的效果是对输入图像进行增强和锐化。参考文献包括Guichard和Morel在尺度空间和形态学中的研究，以及Aubert和Kornprobst在图像处理中的数学问题的探讨。

绘制多条曲线的方法是通过plot(x)函数实现，其中x的列数决定了曲线的数量。如果需要同时绘制多条曲线，可以使用plot(x, y)，确保x和y具有相同维度。例如，plot(x, sin(x),'r:', x, cos(x),'-bo')。

通过链效应，Lissajouse曲线在Matlab / Octave中生成。函数调用示例： LissajouseCurveTable() % 展示6x6的演示表格或者LissajouseCurveTable(10,1) % 展示10x10表格，并保存每一帧。生成的gif文件用于展示动态效果，为用户带来视觉乐趣。

本练习展示了倒谱提升器在语音信号频谱平滑中的应用。通过计算窗口加权的实际倒谱并保存对数幅度谱，随后使用低频率提升器来有效平滑频谱，截止频率从20Hz逐步增至100Hz，步长为20Hz。详细操作指南请参阅“3.14Cepstral Liftering.pdf”。

本研究介绍了一种方法，用于交换两幅图像的幅度谱和相位谱，并利用交换后的谱实现双谱重构。该方法包括幅度谱和相位谱的交换算法、双谱计算方法以及实验验证。实验结果表明，该方法能够有效地交换图像谱，并实现双谱重构，为图像处理和分析提供了新的可能性。