【优化求解】基于School Based Optimization (SBO)求解单目标

1 简介

2 部分代码

登录后复制

          
          

           
           %% This function implements the basic School Based Optimization (SBO) algorithm for 10-bar truss optimization
          
          

• 1.

登录后复制

          
          

           
           %%
          
          

• 1.

登录后复制

          
          

           
           clc
          
          

• 1.

登录后复制

          
          

           
           clear all
          
          

• 1.

登录后复制

          
          

           
           close all
          
          

• 1.

登录后复制

          
          

           
           global D
          
          

• 1.

登录后复制

          
          

           
           % Specity SBO parameters
          
          

• 1.

登录后复制

          
          

           
           Itmax=300;                                                                 % Maximum number of iterations
          
          

• 1.

登录后复制

          
          

           
           NClass=5;                                                                  % Number of classes in the school
          
          

• 1.

登录后复制

          
          

           
           PopSize=15;                                                                % Population size of each class
          
          

• 1.

登录后复制

          
          

           
           % Optimization problem parameters
          
          

• 1.

登录后复制

          
          

           
           D=Data10;                                                                  % For truss function evaluate the functio to get the initial parameters
          
          

• 1.

登录后复制

          
          

           
           LB=D.LB;                                                                   % Lowerbound
          
          

• 1.

登录后复制

          
          

           
           UB=D.UB;                                                                   % Upperbound
          
          

• 1.

登录后复制

          
          

           
           FN='ST10';                                                                 % Name of analyzer function
          
          

• 1.

登录后复制

          
          

           
           %% Randomely generate initial designs between LB and UB
          
          

• 1.

登录后复制

          
          

           
           Cycle=1;
          
          

• 1.

登录后复制

          
          

           
           for I=1:PopSize
          
          

• 1.

登录后复制

          
          

           
           for NC=1:NClass
          
          

• 1.

登录后复制

          
          

           
           Designs{NC}(I,:)=LB+rand(1,size(LB,2)).*(UB-LB);                   % Row vector
          
          

• 1.

登录后复制

          
          

           
           end
          
          

• 1.

登录后复制

          
          

           
           end
          
          

• 1.

登录后复制

          
          

           
           % Analysis the designs
          
          

• 1.

登录后复制

          
          

           
           for NC=1:NClass
          
          

• 1.

登录后复制

          
          

           
           [PObj{NC},Obj{NC}]=Analyser(Designs{NC},FN);
          
          

• 1.

登录后复制

          
          

           
           Best{NC}=[];
          
          

• 1.

登录后复制

          
          

           
           end
          
          

• 1.

登录后复制

          
          

           
           %% SBO loop
          
          

• 1.

登录后复制

          
          

           
           for Cycle=2:Itmax
          
          

• 1.

登录后复制

          
          

           
           for NC=1:NClass
          
          

• 1.

登录后复制

          
          

           
           % Identify best designs and keep them
          
          

• 1.

登录后复制

          
          

           
           [Best{NC},Designs{NC},PObj{NC},Obj{NC},WMeanPos{NC}]=Specifier(PObj{NC},Obj{NC},Designs{NC},Best{NC});
          
          

• 1.

登录后复制

          
          

           
           TeachersPObj(NC,1)=Best{NC}.GBest.PObj;
          
          

• 1.

登录后复制

          
          

           
           TeachersDes(NC,:)=Best{NC}.GBest.Design;
          
          

• 1.

登录后复制

          
          

           
           end
          
          

• 1.

登录后复制

          
          

           
           for NC=1:NClass
          
          

• 1.

登录后复制

          
          

           
           % Select a teacher
          
          

• 1.

登录后复制

          
          

           
           SelectedTeacher=TeacherSelector(Best,NC,TeachersPObj);
          
          

• 1.

登录后复制

          
          

           
           % Apply Teaching
          
          

• 1.

登录后复制

          
          

           
           [Designs{NC},PObj{NC},Obj{NC}]=Teaching(LB,UB,Designs{NC},PObj{NC},Obj{NC},TeachersDes(SelectedTeacher,:),WMeanPos{NC},FN);
          
          

• 1.

登录后复制

          
          

           
           [Best{NC},Designs{NC},PObj{NC},Obj{NC},WMeanPos{NC}]=Specifier(PObj{NC},Obj{NC},Designs{NC},Best{NC});
          
          

• 1.

登录后复制

          
          

           
           % Apply Learning
          
          

• 1.

登录后复制

          
          

           
           [Designs{NC},PObj{NC},Obj{NC}]=Learning(LB,UB,Designs{NC},Obj{NC},PObj{NC},FN);
          
          

• 1.

登录后复制

          
          

           
           [Best{NC},Designs{NC},PObj{NC},Obj{NC},WMeanPos{NC}]=Specifier(PObj{NC},Obj{NC},Designs{NC},Best{NC});
          
          

• 1.

登录后复制

          
          

           
           end
          
          

• 1.

登录后复制

          
          

           
           % Find best so far solution and Mean
          
          

• 1.

登录后复制

          
          

           
           CumPObj=[];
          
          

• 1.

登录后复制

          
          

           
           for NC=1:NClass
          
          

• 1.

登录后复制

          
          

           
           ClassBestPObj(NC,1)=Best{NC}.GBest.PObj;
          
          

• 1.

登录后复制

          
          

           
           ClassMean(NC,1)=mean(PObj{NC});
          
          

• 1.

登录后复制

          
          

           
           CumPObj=[CumPObj;PObj{NC}];
          
          

• 1.

登录后复制

          
          

           
           end
          
          

• 1.

登录后复制

          
          

           
           [~,b]=min(ClassBestPObj);
          
          

• 1.

登录后复制

          
          

           
           OveralBestPObj=Best{b}.GBest.PObj;
          
          

• 1.

登录后复制

          
          

           
           OveralBestObj=Best{b}.GBest.Obj;
          
          

• 1.

登录后复制

          
          

           
           OveralBestDes=Best{b}.GBest.Design;
          
          

• 1.

登录后复制

          
          

           
           % Plot time history of the best solution vs. iteration and print the
          
          

• 1.

登录后复制

          
          

           
           % results
          
          

• 1.

登录后复制

          
          

           
           hold on;plot(Cycle,Best{b}.GBest.PObj,'b*');xlabel('Iteration');ylabel('Best solution value');pause(0.0001)
          
          

• 1.

登录后复制

          
          

           
           fprintf('Cycle: %6d, Best (Penalized): %6.4f, Objective: %6.4f\n',Cycle,OveralBestPObj,OveralBestObj);
          
          

• 1.

登录后复制

          
          

           
           end
          
          

• 1.

登录后复制

          
          

           
           Solution.PObj=OveralBestPObj;% Objective value for best non-penalized solution
          
          

• 1.

登录后复制

          
          

           
           Solution.Design=OveralBestDes;% Design for best non-penalized solution
          
          

• 1.

登录后复制

          
          

           
           img =gcf;  %获取当前画图的句柄
          
          

• 1.

登录后复制

          
          

           
           print(img, '-dpng', '-r600', './img.png')         %即可得到对应格式和期望dpi的图像
          
          

• 1.

登录后复制

          
          

           
           %% Save the results
          
          

• 1.

登录后复制

          
          

           
           save('SBO_Results.mat','Solution')
          
          

• 1.

3 仿真结果

4 参考文献

[1] Farshchin, M. , et al. "School based optimization algorithm for design of steel frames." Engineering Structures 171(2018):326-335.