说明：MMATLAB四阶龙格库塔法 求解微分方程数值解 源程序代码，可以直接运行，对于时间步长是稳定的
说明：%基于VQ的说话人识别系统，矢量量化起着双重作用。在训练阶段，把每一个说话者所提取的特征参数进行分类，产生不同码字所组成的码本。在识别(匹配)阶段，我们用VQ方法计算平均失真测度(本系统在计算距离d时，采用欧氏距离测度)，从而判断说话人是谁。本段代码是对音频进行训练，也就是提取特征参数
说明：数据模型归一化，模态振动，使用高阶累积量对MPSK信号进行调制识别，一种噪声辅助数据分析方法,加入重复控制，本程序的性能已经达到较高水平，应用小区域方差对比，程序简单。
说明：求频响函数的实频、虚频、幅频、相频及Nyqiust图！在读取实验所得振动数据后，通过程序提取模态特性，利用拟合圆法，求三阶固有频率、阻尼比。计算各阶主振型。