2.4 运行JSOP
在上一节中,我们将适应度函数设定为Rastrigin函数,现在我们就通过JSOP来寻找该函数的全局最小值(JSOP默认寻找函数的最小值而非最大值)。
打开usr.m,相关参数设置如下
%仿真设置 model_name='test.cst'; %CST模型路径 result_route='1D Results\S-Parameters'; %仿真结果在Navigation Tree中的路径 result_name='SZmax(1),Zmax(1)'; %仿真结果的名称(参见rtdata.txt) plot_mode='phase'; %仿真结果的模式,包括'magnitude','magnitudedb', 'phase', 'polar', 'real', 'imaginary', 'smith', 'smithy' solver_type='FDSolver'; %频域求解器:FDSolver 时域求解器:Solver start_cst=false; %是否调用CST %变量设置 var_name={'d','w'}; %变量名称 var_bound=[[-10,10];[-10,10]]; %变量变化范围 [ [变量1的下界,变量1的上界] ; [变量2的下界,变量2的上界] ; ...] inserts=[ ]; %插入指定个体 var_constrain=''; %变量的限制条件,若无限制请输入'' %种群设置 MAXGEN=100; %最大代数 NIND=40; %一代中个体的数量 PRECI=20; %单个变量的二进制位数 MTR=1; %变异率(默认值的倍数) %其它设置 plot_graph=2; %0:不绘制进化曲线 1:每进化一代就更新进化曲线 2:完成最后一代后再绘制进化曲线 plot_ave=false; %是否显示种群平均值 %---------------------------------------------------------------------------------------- main_function( model_name, result_route, result_name, plot_mode, solver_type, ... var_name, var_bound', inserts, var_constrain, ... MAXGEN, NIND, PRECI, MTR, ... plot_graph, plot_ave, start_cst); %JSOP1.2
由于不需要调用CST,因此“仿真设置”中除了start_cst,其它的变量可任意设置,只要没有语法错误即可。
运行usr.m,几秒钟后程序即运行完毕并显示出进化曲线图,如下图所示。
每一代的最优个体的适应度值输出在命令行窗口,可以看到,当种群进化到第40代,就基本找到了全局最小值0。
打开JSOP\interface\result.txt,可以看到每一代的每个个体的参数和适应度值,如下图所示,图中展示了某一代的第35-40号个体的信息。在每一代信息的最下方还有每一代的最优个体和平均适应度值,在图中,最优个体为17号个体,平均适应度值为8.0588
如果JSOP没有优化出期望的结果,可以重新运行JSOP进行优化,运行时JSOP会先将result.txt中的数据备份到同一文件夹下的result.old.txt中,再将result.txt清空。需要注意的是,result.old.txt只对上一次的运行结果进行备份,因此若要保留运行结果,建议对result.txt进行手动备份。