這篇文章給大家介紹Python中怎么實(shí)現(xiàn)線性規(guī)劃,內(nèi)容非常詳細(xì),感興趣的小伙伴們可以參考借鑒,希望對(duì)大家能有所幫助。
小店網(wǎng)站制作公司哪家好,找創(chuàng)新互聯(lián)!從網(wǎng)頁(yè)設(shè)計(jì)、網(wǎng)站建設(shè)、微信開(kāi)發(fā)、APP開(kāi)發(fā)、響應(yīng)式網(wǎng)站開(kāi)發(fā)等網(wǎng)站項(xiàng)目制作,到程序開(kāi)發(fā),運(yùn)營(yíng)維護(hù)。創(chuàng)新互聯(lián)成立于2013年到現(xiàn)在10年的時(shí)間,我們擁有了豐富的建站經(jīng)驗(yàn)和運(yùn)維經(jīng)驗(yàn),來(lái)保證我們的工作的順利進(jìn)行。專注于網(wǎng)站建設(shè)就選創(chuàng)新互聯(lián)。
運(yùn)籌學(xué)
運(yùn)籌學(xué)是一種科學(xué)的決策方法,它通常是在需要分配稀缺資源的條件下,尋求系統(tǒng)的優(yōu)秀設(shè)計(jì)??茖W(xué)的決策方法需要使用一個(gè)或多個(gè)數(shù)學(xué)模型(優(yōu)化模型)來(lái)做出最優(yōu)決策。
優(yōu)化模型試圖在滿足給定約束的決策變量的所有值的集合中,找到優(yōu)化(最大化或最小化)目標(biāo)函數(shù)的決策變量的值。 它的三個(gè)主要組成部分是:
鴻蒙官方戰(zhàn)略合作共建——HarmonyOS技術(shù)社區(qū)
目標(biāo)函數(shù):要優(yōu)化的函數(shù)(最大化或最小化)。
決策變量:影響系統(tǒng)性能的可控變量。
約束:決策變量的一組約束(即線性不等式或等式)。非負(fù)性約束限制了決策變量取正值。
優(yōu)化模型的解稱為最優(yōu)可行解。
建模步驟
對(duì)運(yùn)籌學(xué)問(wèn)題進(jìn)行準(zhǔn)確建模是很重要的任務(wù),也是很困難的任務(wù)。錯(cuò)誤的模型會(huì)導(dǎo)致錯(cuò)誤的解決方案,從而不能解決原來(lái)的問(wèn)題。團(tuán)隊(duì)成員應(yīng)按照以下步驟進(jìn)行建模:
鴻蒙官方戰(zhàn)略合作共建——HarmonyOS技術(shù)社區(qū)
問(wèn)題定義:定義項(xiàng)目的范圍,并確定三個(gè)要素:決策變量、目標(biāo)和限制(即約束)。
模型構(gòu)建:將問(wèn)題定義轉(zhuǎn)化為數(shù)學(xué)關(guān)系。
模型求解:使用標(biāo)準(zhǔn)優(yōu)化算法。在獲得解后,需要進(jìn)行靈敏度分析,以找出由于某些參數(shù)的變化而導(dǎo)致的解的行為。
模型有效性:檢查模型是否按預(yù)期工作。
實(shí)現(xiàn):將模型和結(jié)果轉(zhuǎn)換為解決方案。
線性規(guī)劃
線性規(guī)劃(Linear Programming,也稱為L(zhǎng)P)是一種運(yùn)籌學(xué)技術(shù),當(dāng)當(dāng)所有的目標(biāo)和約束都是線性的(在變量中)并且當(dāng)所有的決策變量都是連續(xù)的時(shí)使用。線性規(guī)劃是最簡(jiǎn)單的運(yùn)籌學(xué)方法。
Python的SciPy庫(kù)包含用于解決線性編程問(wèn)題的linprog函數(shù)。在使用linprog時(shí),編寫(xiě)代碼要考慮的兩個(gè)注意事項(xiàng):
鴻蒙官方戰(zhàn)略合作共建——HarmonyOS技術(shù)社區(qū)
這個(gè)問(wèn)題必須表述為一個(gè)最小化問(wèn)題。
不等式必須表示為≤。
最小化問(wèn)題
讓我們考慮以下要解決的最小化問(wèn)題:
讓我們看一下Python代碼:
# Import required libraries import numpy as np from scipy.optimize import linprog # Set the inequality constraints matrix # Note: the inequality constraints must be in the form of <= A = np.array([[-1, -1, -1], [-1, 2, 0], [0, 0, -1], [-1, 0, 0], [0, -1, 0], [0, 0, -1]]) # Set the inequality constraints vector b = np.array([-1000, 0, -340, 0, 0, 0]) # Set the coefficients of the linear objective function vector c = np.array([10, 15, 25]) # Solve linear programming problem res = linprog(c, A_ub=A, b_ub=b) # Print results print('Optimal value:', round(res.fun, ndigits=2), '\nx values:', res.x, '\nNumber of iterations performed:', res.nit, '\nStatus:', res.message)
輸出結(jié)果:
# Optimal value: 15100.0 # x values: [6.59999996e+02 1.00009440e-07 3.40000000e+02] # Number of iterations performed: 7 # Status: Optimization terminated successfully.
最大化問(wèn)題
由于Python的SciPy庫(kù)中的linprog函數(shù)是用來(lái)解決最小化問(wèn)題的,因此有必要對(duì)原始目標(biāo)函數(shù)進(jìn)行轉(zhuǎn)換。通過(guò)將目標(biāo)函數(shù)的系數(shù)乘以-1(即通過(guò)改變其符號(hào)),可以將最小化問(wèn)題轉(zhuǎn)化為一個(gè)最大化問(wèn)題。
讓我們考慮下面需要解決的最大化問(wèn)題:
讓我們看一下Python的實(shí)現(xiàn):
# Import required libraries import numpy as np from scipy.optimize import linprog # Set the inequality constraints matrix # Note: the inequality constraints must be in the form of <= A = np.array([[1, 0], [2, 3], [1, 1], [-1, 0], [0, -1]]) # Set the inequality constraints vector b = np.array([16, 19, 8, 0, 0]) # Set the coefficients of the linear objective function vector # Note: when maximizing, change the signs of the c vector coefficient c = np.array([-5, -7]) # Solve linear programming problem res = linprog(c, A_ub=A, b_ub=b) # Print results print('Optimal value:', round(res.fun*-1, ndigits=2), '\nx values:', res.x, '\nNumber of iterations performed:', res.nit, '\nStatus:', res.message)
上述代碼的輸出結(jié)果為:
# Optimal value: 46.0 # x values: [5. 3.] # Number of iterations performed: 5 # Status: Optimization terminated successfully.
關(guān)于Python中怎么實(shí)現(xiàn)線性規(guī)劃就分享到這里了,希望以上內(nèi)容可以對(duì)大家有一定的幫助,可以學(xué)到更多知識(shí)。如果覺(jué)得文章不錯(cuò),可以把它分享出去讓更多的人看到。