numerical-methods

5熱度

2回答

我希望這個問題（以及它的可能答案）在範圍內足夠一般，這樣對其他人也是有用的。我試圖解決的數值問題，涉及到的功能double的產物，其形式爲其中是預定義的函數，以及i傳遞到積分。的複雜之處在於我的函數不是靜態的;整合是反覆做，並在每個整合■找不同的形式，例如，在第一循環也可能是第二變得和如此類推，其中表示函數的向量的成員。在積分的每次迭代中，我不知道的形式是什麼，只是它將是具有常係數的線性

0熱度

2回答

如何使高斯包在C中的方形屏障的數值模擬移動

我想用高斯包通過Trotter-Suzuki公式和快速傅里葉變換（FFT）研究傳輸概率，當遇到方形屏障時，只需如完成in this Quantum Python article。但我需要用C實現它。原則上，波函數在與方形障礙碰撞之前將保持其形狀。但是我發現波函數在與方形障礙碰撞之前會隨着時間顯着變平。有人在以下代碼中發現問題？這裏，兩個文件 - 結果和psi.txt - 創建存儲初始和演變波函數

0熱度

1回答

fsolve/fzero：沒有找到解決辦法，定期出現

我試圖做使用fsolve或fzero以下算法： K5=8.37e-2 P=1 Choose an A S2=(4*K5/A)^(2/3) S6=3*S2 S8=4*S2 SO2 = (5*P)/149 - (101*S2)/149 - (293*S6)/149 - (389*S8)/149 H2O = (40*P)/149 + (556*S

0熱度

1回答

半隱式Euler數值積分方法的時間可逆性

很明顯的是，semi-implicit Euler integration法辛，但我無法找到任何信息關於它的time-reversibility。所以這個問題：它是時間可逆的嗎？

0熱度

1回答

比較計算浮點數的準確性與x * 0.1和x/10.0

哪個表達式更準確，爲什麼？我認爲如果它不是10.0並且對於某個整數k是2^k，那麼我們只需要在指數中減去k即可。但有10個，我想不出一種比較它們的方法......

2熱度

2回答

Haskell中的數值積分：錯誤的符號

我已經寫了一些代碼，它是用梯形法則數值化的函數。它有效，但它產生的答案是錯誤的。爲什麼會這樣？的代碼是： integration :: (Double -> Double) -> Double -> Double -> Double integration f a b = h * (f a + f b + partial_sum) where h = (b - a)/10

5熱度

1回答

numpy.sum是否以避免數值錯誤的方式實現？

衆所周知，合計數字會導致數字錯誤（例如，如果第一個數字非常大，而其他許多小數字）。這可以通過以非直接方式累加數字來解決。請參閱：https://en.wikipedia.org/wiki/Kahan_summation_algorithm numpy.sum是否以避免數值錯誤的方式實現？

0熱度

1回答

如何繪製時間序列圖而不記住整個系列？

時間序列長無限期，但我希望能夠看看從開始時間t0到當前的時間階段t的情節。這意味着，該地塊的時間尺度是不斷萎縮。當然，我可以把所有的時間序列數據，並在每個時間步我只是重新繪製新規模的情節。但是如果內存也有限呢？假設編程被允許僅保留N個的存儲器陣列（他們可以存儲他們想要的陣列中的任何數目）。儘管使用有限的內存，程序員如何繪製一個保留真實圖形近似形狀的圖形？

-5熱度

1回答

程序評估多項式ax 3 + bx 2 + cx + d對於a，b，c和d的給定值的最小操作數

請建議子程序來評估多項式ax 3 + bx 2 + cx + d對於a，b，c和d的給定值，具有最少的操作次數。如果使用平分方法有什麼方法可以動態地猜測極限值。

8熱度

3回答

argsort爲多維的ndarray

我試圖讓索引按最後一個軸排序多維數組，例如， >>> a = np.array([[3,1,2],[8,9,2]]) 而且我想指數i使得 >>> a[i] array([[1, 2, 3], [2, 8, 9]]) 基礎的numpy.argsort我認爲應該做到這一點的文件上，但我得到的錯誤： >>> a[np.argsort(a)] IndexError: index 2