numerical-methods

0熱度

1回答

我試圖用Matlab編寫歐拉方法的'啞'版本，但我總是沒有想到什麼。我的代碼是垃圾:-( 這種方法請參見本僞代碼： ‘set integration range xi = 0 xf = 0 ‘initialize variables x = xi y = 1 ‘set step size and determine ‘number of calculation steps dx =

7熱度

4回答

JavaScript數字圖書館 - 向量，矩陣，確定，反演

純JavaScript中是否有合理的數值庫？我想要基於陣列的向量，矩陣，確定和矩陣反演。這隻適用於小問題，不適用於大問題。我知道在Python，Java，C++中有選擇，或者哦，也許是FORTRAN，更不用說R.嗯...... FORTRAN ...... 我會很高興得到基本矢量和矩陣運算，而無需在FORTRAN中重建輪或重新編碼一兩天。

0熱度

3回答

增量縮放功能

如果您知道'input1'嚴格介於0和1之間，或者通常'min'和'max'（其中min和max已知在介於0之間但不是嚴格介於0和1之間）你會得到'input1'來通過'input2'給出的數值跳躍來增加或減少，保證新值嚴格在最小值和最大值之間，並且永遠不會達到最小值或最大值？

0熱度

2回答

修改rk4方法的C實現

坦率地說，我的問題是，我不確定這是如何工作的。我需要修改雙f（）函數來解決任意微分方程d2θ/ dt2 =-ω2sinθ，但因爲它是我不確定如何繼續。 rk4函數runge4（）本身我明白;我不明白的是f（）函數如何爲諧振子返回正確的值。有人請至少解釋f（）函數背後的邏輯嗎？原始代碼如下。 /* **********************************************

1熱度

3回答

如何在c語言中將正在程序中計算的變量作爲指針（使用例如梯形法則）進行數值積分

我有一個代碼，並非由我自己編寫。在這個複雜的代碼中，許多規則被應用於計算數量d（x）。在代碼中正在使用一個指針來計算它。我想計算一個積分，例如： W = Int_0^L d（x）dx？我這樣做： #define DX 0.003 void WORK(double *d, double *W) { double INTE5=0.0; int N_X_POINTS=333

3熱度

1回答

最喜歡的/最好的A.X = B解算器支持長時間雙倍？

我希望反覆錘擊X = B（A^-1）的問題。也就是說，解決一個線性系統。對於C++，哪些數值求解器支持128位長雙精度（四邊形）？使用C風格的數組是一個主要優點，因爲我所有的2D數據都存儲爲一個std :: vector。我一直希望能用GCC或ICC編譯Linux上的代碼。

0熱度

1回答

如何設置UMFPACK的容差

我在C++中使用umfpack來解決稀疏矩陣。我無法找到[UMFPACK PIVOT TOLERANCE]參數設置的位置。任何人都會知道如何設置這個控制參數嗎？我當前的代碼是 umf::symbolic_type<double> Symbolic; umf::numeric_type<double> Numeric; umf::symbolic(J, Symbolic); umf::num

3熱度

2回答

R中的全奇異值分解

在大多數應用（尤其是統計學應用）中，thin SVD就足夠了。然而，有時需要全SVD才能獲得矩陣（及其共軛）的零空間的正交性。看來R中的svd（）只返回精簡版本。是否有可能生成完整版本？有替代品嗎？

6熱度

1回答

如何使用DotNumerics解決線性規劃問題？

我對數值分析非常感興趣。我一直在使用DotNumerics開源應用程序。我的線性系統如下： 1 * x + 3 * y <= 150 2 * x + 1 * y <= 100 其中x >= 0, y >= 0 z = 10 * x + 15 * y 我試圖解決z（優化...） I can use Simplex method to solve above problem as found

2熱度

2回答

用於快速計算Numpy中大數集標準偏差的方法

執行此操作的最佳（最快）方法是什麼？這會產生什麼，我相信這是正確的答案，但顯然在N = 10E6是痛苦的緩慢。我認爲我需要保持Xi值，以便我可以正確計算標準偏差，但是有沒有什麼技術可以讓這個運行更快？ def randomInterval(a,b): r = ((b-a)*float(random.random(1)) + a) return r N = 10e6 Su