solver

1熱度

1回答

我想向矩陣中添加更多的數據來分析和解決問題，但是因爲它目前正在執行粗暴操作，所以如果我在分析中添加另一列，則會超出python的極限。有沒有解決方法可以找到類似的結果，而不是不得不通過組合來暴力行爲？ sample.csv也在下面列出。感謝您的任何建議。 import csv import itertools as it import numpy as np C = 2618.08 B

0熱度

1回答

VBA解算器每次粘貼結果

我已經搜索和搜索，我無法弄清楚這一點。很感謝任何形式的幫助。我正在運行一個基本的單純形求解器。根據預測（H列）和薪金（G列），我希望得出最高數額$ A $ 4。我的約束是穩定的，但我無法循環，然後在某處粘貼結果。二進制解決方案在列A中，總共有6個結果爲200.我在「優化程序1」中擁有所有這些數據，並且希望每次都將結果粘貼到「列表」中，從單元格A1-A6開始，然後第二個循環將是A7-A13等，我想

0熱度

1回答

Sympy在另一個求解器中重用求解的結果

下午好，我遇到了我的代碼問題。我將要展示的代碼是實際代碼的簡化版本，但想法是相同的。 Var1 = Symbol('Var1') Var2 = Symbol('Var2') A = 20 B = 30 Var1 = solve(12+A*B+Var1, Var1) Var2 = solve(Var1+A+B+Var2, Var2) print(Var1,Var2) 什麼問題是，例

1熱度

1回答

SymPy：如何使`nsolve`返回一個字典，找到的解決方案

我想用nsolve作爲後備，以solve，並喜歡用dict = True使solve回報與發現的解決方案和相應的變量的字典。但是，nsolve似乎沒有此選項。這是我使用的是什麼樣的解決方法： from sympy import * def nsolve(equations, variables, guesses, **flags): from sympy import nsolve

3熱度

1回答

SymPy：爲什麼如果公式列表中的「solve」有時會失敗？

按照SymPy documentation on solve，第一參數可以是：單個Expr的或聚必須是零，一個平等的關係表達式或一個布爾迭代或更多以上但是，有時它確實有所作爲，如果它是solve(equation, x)或。 [equation]應該是equation的迭代，對吧？我真的很困惑，爲什麼它只是一個問題。你可以看到我在下面MCVE問題： from sympy impor

0熱度

1回答

Excel VBA求解器循環花費時間過長

我創建了一個代碼，它通過求解器通過4個循環（每次15'迭代）完成60次優化。該代碼工作得很好，但它需要永久運行（超過一個小時）。每個優化都是一個簡單的線性模型（找到全局解決方案），只是改變我正在查看的月份。問題設置：我試圖儘量減少工人/測試平臺的數量，以滿足攝入量和容量不足的限制。我不知道我怎麼能讓它走得更快，但我不能把這個發送給其他人，並期望他們使用它。有人有任何建議嗎？下面是我的代碼：

0熱度

1回答

用數組中的符號代替數字值

我面臨的問題如下：我已經給出了一個包含2個NaN的6x1數組，代表2個未知值。我想要做的是用符號值替換這些值（例如，x_1和x_2）。我不知道該怎麼做，如果我給了一個12x1陣列，例如10 NaN，我必須改爲X_1，X_2 ... X_10，以便之後使用解算器。 array = [0; 0; NaN; 0; 0; 0; 0; NaN; 0]; 更改爲： array = [0; 0; x_1;

1熱度

1回答

添加函數時，Mathematica在步長相當大時忽略一個

首先，這個問題在視覺上更容易理解，我會張貼圖像，但我是新來的。我使用NDSolve求解兩個耦合微分方程（用於溫度和密度）。溫度方程以高斯形式添加了加熱功能，我希望減少高斯（降低方差），直到它接近delta函數，但是當我將方差降低到某個點時，NDSolve開始忽略加熱功能，可能與步長過大有關？下面是我使用的代碼：開始與一些行話： a = 1.99*10^-9; b = 0.24*10^-3;

1熱度

1回答

改善與涉及集成

我嘗試使用下面的代碼（不相關的部分去掉）來解決積分方程方程scipy.optimize.fsolve精度： def _pdf(self, a, b, c, t): pdf = some_pdf(a,b,c,t) return pdf def _result(self, a, b, c, flag): return fsolve(lambda t: flag - 1

0熱度

1回答

Python的解決給定的條件

大家是否有可能解決如下所示： x = np.matrix([[8.3,-20.6],[-20.6,65.8]]) y是P的函數： y = lambda P: np.matrix([[0.02P,-0.02P], [-0.02P,0.04P]]) 我想找到給出P值條件： P>0, det(x-y)==0; 有沒有方便的方法呢？非常感謝！ Shawn