头脑风暴算法MATLAB版本

头脑风暴算法的MATLAB实现，完美运行。运行成功后，显示曲线图

function  [dU] = analysis_calculation（L1dertHC）

%电磁校核计算程序
%变量公用区
global EFF Bt Bj J Sf;
      
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
                  %注释：双层、星接、径向、分数槽绕组、平行齿梨型槽%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%初始数据赋值%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

PN = 3;                                                                    %电机额定功率（kW）
UN = 231;                                                                  %电机额定相电压（V）
m = 3;                                                                     %电机相数
f = 50;                                                                    %额定频率（Hz）
NN = 1500;                                                                 %额定转速（rpm）
p = 2;                                                                     %电机极对数
dUN = 10;                                                                  %固有电压调整率（%100）
q = 3;                                                                     %每极每相槽数
                                                                           %额定功率因数为0.90
%Br20 = 1.24;                                                              %20°时永磁体剩余磁通密度（T）
%Hc20 = 900;                                                               %20°时永磁体的矫顽力（kA/m）
%HC = ;
%工作温度时的计算矫顽力     优化变量3
aBr = -0.12;                                                               %剩余磁密温度系数（%K-1）
tem = 75;                                                                  %预计工作时温度
mur = 1.05;                                                                %相对回复磁导率（H/m）
mu0 = 4*pi*10^（-7）;                                                        %真空磁导率（H/m）
ap = 0.785;                                                                %极弧系数
ai = 0.785;                                                                %计算极弧系数（表面凸出式结构）
KFe = 0.95;                                                                %定子铁心叠压系数
Q = 36;                                                                    %定子槽数
Ksk  =0.995;                                                               %斜槽因数
D1  =23.0;                                                                 %定子外径（cm）
Di1 = 15.6;                                                                %定子内径（cm）
%L1 = 10.0;
%%定子铁心长度（cm）  优化变量1
hs0 = 0.08;                                                                %槽口高0（cm）
hs1 = 0.12;                                                                %槽口高1（cm）
hs2 = 0.72;                                                                %槽高（cm）          
bs0 = 0.20;                                                                %槽口宽（cm）       
bs1 = 0.58;                                                                %槽宽1（cm）         
bs2 = 0.70;                                                                %槽宽2（cm）
Ns = 56;                                                                   %每槽导体数   

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      24322  2018-11-06 17:02  程序\analysis_calculation.m

     文件       3977  2018-11-06 17:03  程序\BSO_fun.m

     文件       1197  2018-11-06 17:02  程序\constraint.m

     文件      24708  2018-11-06 17:02  程序\fun.m

     文件      24760  2018-11-06 17:02  程序\fun1.m

     文件        423  2018-11-06 17:02  程序\getconstraints.m

     文件         92  2018-11-06 17:02  程序\getH.m

     目录          0  2019-03-30 09:42  程序

----------- ---------  ---------- -----  ----

                79479                    8