目录
- 1.摘要
- 2.麻雀搜索算法SSA原理
- 3.改进策略
- 4.结果展示
- 5.参考文献
- 6.代码获取
 

1.摘要
为了应对暖通空调(HVAC)系统由于不当负荷分配导致的高能源消耗问题,本文提出了一种改进麻雀搜索算法(ISSA)。ISSA算法旨在在满足负载需求的同时,最小化系统的能源消耗,选取每个冷水机组的部分负载率作为优化变量。ISSA算法引入Circle混沌映射以初始化位置,有效提高了初始解的质量和多样性。同时,ISSA算法融合了灰狼优化算法的信息交换加强机制更新生产者位置,并通过混沌正弦余弦更新策略,从而增强了优化精度与算法的收敛速度。

PS:主要复现改进麻雀搜索算法~
2.麻雀搜索算法SSA原理
【智能算法】麻雀搜索算法(SSA)原理及实现
3.改进策略
Circle混沌映射
麻雀搜索算法因其初始种群随机生成,难以保证个体在搜索空间中的均匀分布,进而影响算法的搜索速度和优化性能。混沌映射利用混沌变量的遍历性和随机性,将待优化的PLR变量按照混沌映射规则映射至混沌变量的值区间,并通过线性变换将混沌序列转化为目标函数的搜索空间:
  
      
       
        
         
         
           X 
          
          
          
            i 
           
          
            , 
           
          
            j 
           
          
          
          
            t 
           
          
            + 
           
          
            1 
           
          
         
        
          = 
         
         
         
           m 
          
         
           o 
          
         
           d 
          
         
         
         
           ( 
          
          
          
            X 
           
           
           
             i 
            
           
             , 
            
           
             j 
            
           
          
            t 
           
          
         
           + 
          
         
           0.2 
          
         
           − 
          
          
          
            ( 
           
           
           
             0.5 
            
            
            
              2 
             
            
              π 
             
            
           
          
            ) 
           
          
          
          
            s 
           
          
            i 
           
          
            n 
           
          
          
          
            ( 
           
          
            2 
           
          
            π 
           
           
           
             X 
            
            
            
              i 
             
            
              , 
             
            
              j 
             
            
           
             t 
            
           
          
            ) 
           
          
         
           , 
          
         
           1 
          
         
           ) 
          
         
        
       
         X_{i,j}^{t+1}=\mathrm{mod}\left(X_{i,j}^t+0.2-\left(\frac{0.5}{2\pi}\right)\mathrm{sin}\left(2\pi X_{i,j}^t\right),1\right) 
        
       
     Xi,jt+1=mod(Xi,jt+0.2−(2π0.5)sin(2πXi,jt),1)
改进生产者位置更新
在处理最优冷水机组加载(OCL)问题时,至关重要的一点是精确控制并联冷水机组的运作,以便在满足多样化负载需求的同时,最大限度地降低电力消耗。然而,在麻雀搜索算法(SSA)的初始搜索阶段,生产者倾向于快速聚焦于全局最优解,这种行为减少了种群的多样性,并可能导致算法过早地收敛到局部最优解。为解决这一问题,本文引入了灰狼优化算法中的信息交换加强机制,用于更新生产者的位置:
  
      
       
        
         
         
           X 
          
          
          
            i 
           
          
            , 
           
          
            j 
           
          
          
          
            t 
           
          
            + 
           
          
            1 
           
          
         
        
          = 
         
         
         
           X 
          
          
          
            i 
           
          
            , 
           
          
            j 
           
          
         
           t 
          
         
        
          + 
         
        
          A 
         
        
          ⋅ 
         
         
         
           D 
          
          
          
            g 
           
          
            b 
           
          
            e 
           
          
            s 
           
          
            t 
           
          
         
        
          + 
         
        
          B 
         
        
          ⋅ 
         
         
         
           D 
          
          
          
            s 
           
          
            u 
           
          
            b 
           
          
            b 
           
          
            e 
           
          
            s 
           
          
            t 
           
          
         
        
       
         X_{i,j}^{t+1}=X_{i,j}^t+A\cdot D_{gbest}+B\cdot D_{subbest} 
        
       
     Xi,jt+1=Xi,jt+A⋅Dgbest+B⋅Dsubbest
 其中,参数表述为:
  
      
       
        
         
          
           
            
           
          
          
           
            
             
            
              A 
             
            
              = 
             
            
              2 
             
            
              d 
             
            
              ξ 
             
            
              − 
             
            
              d 
             
            
           
          
         
         
          
           
            
           
          
          
           
            
             
            
              d 
             
            
              = 
             
             
             
               { 
              
              
               
                
                 
                  
                  
                    1 
                   
                  
                    + 
                   
                   
                    
                    
                      [ 
                     
                    
                      cos 
                     
                    
                       
                     
                     
                     
                       ( 
                      
                      
                      
                        t 
                       
                      
                        M 
                       
                      
                     
                       π 
                      
                     
                       ) 
                      
                     
                    
                      ] 
                     
                    
                   
                     ′ 
                    
                   
                  
                    , 
                   
                  
                    0 
                   
                  
                    < 
                   
                  
                    t 
                   
                  
                    ≤ 
                   
                   
                   
                     M 
                    
                   
                     2 
                    
                   
                  
                 
                
               
               
                
                 
                  
                 
                
               
               
                
                 
                  
                  
                    1 
                   
                  
                    − 
                   
                   
                    
                    
                      [ 
                     
                    
                      − 
                     
                    
                      cos 
                     
                    
                       
                     
                     
                     
                       ( 
                      
                      
                      
                        t 
                       
                      
                        M 
                       
                      
                     
                       π 
                      
                     
                       ) 
                      
                     
                    
                      ] 
                     
                    
                   
                     ′ 
                    
                   
                  
                    , 
                   
                   
                   
                     M 
                    
                   
                     2 
                    
                   
                  
                    < 
                   
                  
                    t 
                   
                  
                    ≤ 
                   
                  
                    M 
                   
                  
                 
                
               
              
             
            
           
          
         
        
       
         \left. \begin{aligned} & A=2d\xi-d \\ & d=\left\{ \begin{array} {ll}1+\left[\cos\left(\frac{t}{M}\pi\right)\right]^{\prime},0<t\leq\frac{M}{2} \\ \\ 1-\left[-\cos\left(\frac{t}{M}\pi\right)\right]^{\prime},\frac{M}{2}<t\leq M \end{array}\right. \end{aligned}\right. 
        
       
     A=2dξ−dd=⎩ 
                       ⎨ 
                       ⎧1+[cos(Mtπ)]′,0<t≤2M1−[−cos(Mtπ)]′,2M<t≤M
改进拾荒者位置更新
在处理高维优化问题时,麻雀搜索算法中的拾荒者位置更新过程高度依赖于先前个体。如果这些个体的位置处于局部最优,拾荒者便可能陷入这一局部最优解并导致优化过程停滞,这限制了算法的整体优化效率。为了克服这一问题,本文引入了正弦余弦机制和混沌算子到拾荒者位置更新:
  
      
       
        
         
          
           
            
             
             
               X 
              
              
              
                i 
               
              
                , 
               
              
                j 
               
              
              
              
                t 
               
              
                + 
               
              
                1 
               
              
             
            
              = 
             
             
             
               X 
              
              
              
                i 
               
              
                , 
               
              
                j 
               
              
             
               ′ 
              
             
            
              + 
             
             
             
               r 
              
             
               1 
              
             
            
              ⋅ 
             
            
              sin 
             
            
               
             
            
              ( 
             
             
             
               r 
              
             
               2 
              
             
            
              ) 
             
            
              ⋅ 
             
             
             
               ∣ 
              
              
              
                r 
               
              
                3 
               
              
             
               ⋅ 
              
              
              
                X 
               
               
               
                 r 
                
               
                 a 
                
               
                 n 
                
               
                 d 
                
               
               
               
                 t 
                
               
                 + 
                
               
                 1 
                
               
              
             
               − 
              
              
              
                X 
               
               
               
                 i 
                
               
                 , 
                
               
                 j 
                
               
              
                ′ 
               
              
             
               ∣ 
              
             
             
             
             
               r 
              
             
               4 
              
             
            
              < 
             
            
              0.5 
             
            
           
          
         
         
          
           
            
             
             
               X 
              
              
              
                i 
               
              
                , 
               
              
                j 
               
              
              
              
                t 
               
              
                + 
               
              
                1 
               
              
             
            
              = 
             
             
             
               X 
              
              
              
                i 
               
              
                , 
               
              
                j 
               
              
             
               ′ 
              
             
            
              + 
             
             
             
               r 
              
             
               1 
              
             
            
              ⋅ 
             
            
              cos 
             
            
               
             
            
              ( 
             
             
             
               r 
              
             
               2 
              
             
            
              ) 
             
            
              ⋅ 
             
             
             
               ∣ 
              
              
              
                r 
               
              
                3 
               
              
             
               ⋅ 
              
              
              
                X 
               
               
               
                 p 
                
               
                 b 
                
               
                 e 
                
               
                 s 
                
               
                 t 
                
               
               
               
                 t 
                
               
                 + 
                
               
                 1 
                
               
              
             
               − 
              
              
              
                X 
               
               
               
                 i 
                
               
                 , 
                
               
                 j 
                
               
              
                t 
               
              
             
               ∣ 
              
             
             
             
             
               r 
              
             
               4 
              
             
            
              ≥ 
             
            
              0.5 
             
            
           
          
         
        
       
         \begin{gathered} X_{i,j}^{t+1}=X_{i,j}^{\prime}+r_{1}\cdot\sin(r_{2})\cdot\left|r_{3}\cdot X_{rand}^{t+1}-X_{i,j}^{\prime}\right|\quad r_{4}<0.5 \\ X_{i,j}^{t+1}=X_{i,j}^{\prime}+r_{1}\cdot\cos(r_{2})\cdot\left|r_{3}\cdot X_{pbest}^{t+1}-X_{i,j}^{t}\right|\quad r_{4}\geq0.5 \end{gathered} 
        
       
     Xi,jt+1=Xi,j′+r1⋅sin(r2)⋅ 
                      r3⋅Xrandt+1−Xi,j′ 
                      r4<0.5Xi,jt+1=Xi,j′+r1⋅cos(r2)⋅ 
                      r3⋅Xpbestt+1−Xi,jt 
                      r4≥0.5
流程图

4.结果展示

 

 
 
5.参考文献
[1] Xue Z, Yu J, Zhao A, et al. Optimal chiller loading by improved sparrow search algorithm for saving energy consumption[J]. Journal of Building Engineering, 2023, 67: 105980.







![[服务器][教程]Ubuntu24.04 Server开机自动挂载硬盘教程](https://i-blog.csdnimg.cn/img_convert/b96a9ceaa03e976263bb4c54dda7d586.png)











