欧美三区_成人在线免费观看视频_欧美极品少妇xxxxⅹ免费视频_a级毛片免费播放_鲁一鲁中文字幕久久_亚洲一级特黄

模擬退火算法Python實現

系統 1985 0

模擬退火算法Python實現

  • 瞎BB
  • 代碼
    • 導入庫以及參數設置
    • 目標函數
    • 主函數

瞎BB

代碼

導入庫以及參數設置

            
              
                import
              
               matplotlib
              
                .
              
              pyplot 
              
                as
              
               plt

              
                import
              
               math

              
                import
              
               random
T_init 
              
                =
              
              
                100
              
              
                # 初始最大溫度
              
              
alpha 
              
                =
              
              
                0.95
              
              
                # 降溫系數
              
              
T_min 
              
                =
              
              
                1e
              
              
                -
              
              
                3
              
              
                # 最小溫度,即退出循環條件
              
            
          

目標函數

            
              
                def
              
              
                obj
              
              
                (
              
              x
              
                )
              
              
                :
              
              
    y 
              
                =
              
              
                10
              
              
                *
              
               math
              
                .
              
              sin
              
                (
              
              
                5
              
              
                *
              
               x
              
                )
              
              
                +
              
              
                7
              
              
                *
              
               math
              
                .
              
              cos
              
                (
              
              
                4
              
              
                *
              
               x
              
                )
              
              
                return
              
              
                -
              
              y

            
          

主函數

            
              
                def
              
              
                SA
              
              
                (
              
              T_init
              
                ,
              
              alpha
              
                ,
              
              T_min
              
                )
              
              
                :
              
              
    T 
              
                =
              
               T_init
    x_new 
              
                =
              
               random
              
                .
              
              random
              
                (
              
              
                )
              
              
                *
              
              
                10
              
              
                #初解
              
              
    x_current 
              
                =
              
               x_new
    y_current 
              
                =
              
              
                float
              
              
                (
              
              
                'inf'
              
              
                )
              
              
    x_best 
              
                =
              
               x_new
    y_best 
              
                =
              
              
                float
              
              
                (
              
              
                'inf'
              
              
                )
              
              
                while
              
               T 
              
                >
              
               T_min
              
                :
              
              
                for
              
               i 
              
                in
              
              
                range
              
              
                (
              
              
                100
              
              
                )
              
              
                :
              
              
            delta_x 
              
                =
              
               random
              
                .
              
              random
              
                (
              
              
                )
              
              
                -
              
              
                0.5
              
              
                # 自變量變化后仍要求在[0,10]之間
              
              
                if
              
              
                0
              
              
                <
              
              
                (
              
              x_new 
              
                +
              
               delta_x
              
                )
              
              
                <
              
              
                10
              
              
                :
              
              
                x_new 
              
                =
              
               x_new 
              
                +
              
               delta_x
            
              
                else
              
              
                :
              
              
                x_new 
              
                =
              
               x_new 
              
                -
              
               delta_x
            y_new 
              
                =
              
               obj
              
                (
              
              x_new
              
                )
              
              
                if
              
              
                (
              
              y_new
              
                <
              
              y_current
              
                )
              
              
                :
              
              
                y_current
              
                =
              
              y_new
                x_current
              
                =
              
              x_new
                
              
                if
              
              
                (
              
              y_new
              
                <
              
              y_best
              
                )
              
              
                :
              
              
                    y_best
              
                =
              
              y_new
                    x_best
              
                =
              
              x_new
            
              
                else
              
              
                :
              
              
                if
              
               random
              
                .
              
              random
              
                (
              
              
                )
              
              
                <
              
               math
              
                .
              
              exp
              
                (
              
              
                -
              
              
                (
              
              y_new 
              
                -
              
               y_current
              
                )
              
              
                /
              
               T
              
                )
              
              
                :
              
              
                    y_current
              
                =
              
              y_new
                    x_current
              
                =
              
              x_new
                
              
                else
              
              
                :
              
              
                    x_new
              
                =
              
              x_current
        T 
              
                *=
              
               alpha
    
              
                print
              
              
                (
              
              
                '最優解'
              
              
                ,
              
              x_best
              
                ,
              
              obj
              
                (
              
              x_best
              
                )
              
              
                )
              
              
SA
              
                (
              
              T_init
              
                ,
              
              alpha
              
                ,
              
              T_min
              
                )
              
            
          

更多文章、技術交流、商務合作、聯系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦!!!

發表我的評論
最新評論 總共0條評論
主站蜘蛛池模板: 999久久久精品视频在线观看 | 成片在线看一区二区草莓 | 奇米激情 | 天天透天天干 | 国产品久久| 日本中文字幕在线播放 | 国产黄色2| 久久久久久久免费看 | h网站国产 | 欧美精品成人 | 一区二区在线看 | 欧美精品网站 | 精品国产三级 | 亚洲国产片高清在线观看 | 99久久久精品国产一区二区 | 中国黄色一级生活片 | 日韩深夜视频 | 亚洲精品午夜国产va久久成人 | 激情九月| 无码又黄又爽又舒服的A片 综合久久网 | 亚洲国产国产综合一区首页 | 神马九九 | 亚洲欧美视频一区 | 国产精品爱久久久久久久电影 | 丝袜 亚洲 另类 欧美 综合 | 台湾av在线 | 久久久久久久国产精品毛片 | 久草精品在线 | 国产精品视频网 | 欧美三区在线观看 | 亚洲电影一区二区 | 亚洲精品乱码久久久久久蜜桃91 | 美女露尿口视频 | 精品国产成人 | 色接久久| 欧美日韩国产在线人成dvd | 日韩不卡高清 | 日本视频网 | 日本精品久久久久护士 | www.youjizz.com亚洲 | 欧美精品一区二区三区在线 |