针对大家评论区给出的很多问题，作者一直都有关注，因此在这里又写了一篇文章，而且思路与这篇文章有不同之处，至于具体的不同之处放在下一篇文章了，大家感兴趣的可以移步观看，下一篇文章可以说是作者的呕心力作。
(4条消息) 白鲸优化算法优化VMD参数，并提取特征向量，以西储大学数据为例，附MATLAB代码_今天吃饺子的博客-CSDN博客/woaipythonmeme/article/details/131226447?spm=1001..3001.5501好了，废话到此为止！接下来讲正文！
同样以西储大学数据集为例，选用105.mat中的X105_BA_time.mat数据。
首先进行VMD分解，采用麻雀优化算法（SSA）对VMD的两个关键参数（惩罚因子α和模态分解数K）进行优化，以最小包络熵为适应度值。其他智能优化算法同样适用，关键要学会最小包络熵代码的编写，网上的五花八门，代码中会对最小包络熵部分进行详细注释。
先上结果图：
实验过程中，会实时显示每次寻优后的最小包络熵值和VMD对应的两个最佳参数。本次寻优共100次（自己可以随意更改寻优次数）。
可以看到寻优100次后，最小包络熵为7.4036，对应两个vmd的最佳参数为122,8，其中惩罚因子为122，模态分解数为8。
收敛曲线如下所示：
代码包含一个脚本，三个函数！
主函数脚本：其中惩罚因子α的范围是[100-2500]，模态分解数K的范围是[3-10]。这里大家可以自行更改，但是一般建议是这样
clear all clcaddpath(genpath(pwd))CostFunction=@(x) Cost(x); % 适应度函数的调用，包络熵值，详情请看Cost%设置SSA算法的参数Params.nVar=2; % 优化变量数目Params.VarSize=[1 Params.nVar];% Size of Decision Variables MatrixParams.VarMin=[100 3];% 下限值，分别是a,kParams.VarMax=[2500 10]; % 上限值Params.MaxIter=30; % 最大迭代数目Params.nPop=30; % 种群规模[particle3, GlobalBest3,SD,GlobalWorst3,Predator,Joiner] = Initialization(Params,CostFunction,'SSA'); %初始化SSA参数disp(['***采用SSA算法开始寻优***'])[GlobalBest,SSA_curve] = SSA(particle3,GlobalBest3,GlobalWorst3,SD,Predator,Joiner,Params,CostFunction); %采用SSA参数优化VMD的两个参数fMin = GlobalBest.Cost;bestX = GlobalBest.Position;%画适应度函数图figureplot(1:Params.MaxIter,SSA_curve,'Color','r')title('Objective space')xlabel('Iteration');set(gca,'xtick',0:10:Params.MaxIter);ylabel('Best score obtained so far');legend('SSA优化VMD')display(['The best solution obtained by SSA is : ', num2str(round(bestX))]); %输出最佳位置display(['The best optimal value of the objective funciton found by SSA is : ', num2str(fMin)]); %输出最佳适应度值
关于VMD函数的详解和相关图像的绘制，可以参考这个文章：(5条消息) VMD分解，matlab代码，包络线，包络谱，中心频率，峭度值，能量熵，近似熵，包络熵，希尔伯特变换，包含所有程序MATLAB代码，-西储大学数据集为例_今天吃饺子的博客-CSDN博客
完整代码获取：下方卡片回复关键词：SSAVMD
觉着不错的给博主留个小赞吧！您的一个小赞就是博主更新的动力！谢谢！
