numerical-methods

    0熱度

    1回答

    多維拉普拉斯變換實現的數值求逆

    在Mathematica，Matlab和Python等語言中，在一維（例如Talbot或歐拉算法）中對拉普拉斯變換進行數值逆運算的算法有許多實現。 我正在尋找一個實現（任何語言）在多維的情況下，或者我可以自己編碼的東西。如果您知道這一點，請告訴我。
    numerical-methods 2017-03-09

    0熱度

    1回答

    沒有SciPy的Python數值整合？

    我需要整合功能：Y（X）= E -ax COS（x）的不使用SciPy的/ NumPy的等我掙扎了不少。我知道這一點： def integrand(x): return 對於返回位我想用這個公式：從i=1到Y（XI）（XI - X（I-1））N總和 我該怎麼寫這個？
    python math numerical-methods numerical-integration 2017-03-10

    -1熱度

    2回答

    用numpy精確求和

    我正在編寫一個程序，它將數字x和y數據作爲一個numpy數組，然後將其整合（數值）在x的一個區域上。我這樣做了，現在我開始懷疑精度了。 我注意到math.fsum()以某種方式糾正了某些錯誤sum()沒有 - numpy.sum也做出這些更正？ 如果我寫sum = numpy.sum(array)或類似的東西，其結果（總和）是否儘可能精確？我總結了數千個數據點，這些數據點本身就是幾個操作的產物，所
    python numpy numerical-methods 2017-03-15

    4熱度

    3回答

    Runge-Kutta算法C++

    下面是我的4階Runge-Kutta算法來求解一階ODE。我檢查它反對維基百科例如發現here解決： \frac{dx}{dt} = tan(x) + 1 遺憾的是，通過了一點點。我玩了很長一段時間，但我找不到錯誤。答案應該是t = 1.1和x = 1.33786352224364362。下面的代碼給出t = 1.1和x = 1.42223。 /* This code is a 1D cl
    c++ math numerical-methods runge-kutta 2017-03-17

    1熱度

    1回答

    數值積分後顯示的數值範圍不正確

    我想使用MATLAB來繪製下面的公式的誤差概率，我想使用命令trapz進行數值積分，問題是我得到了一個很好的形狀但是y軸的值是錯誤的，整個曲線應該在0到1.2之間，但是在0.492到0.5之間！任何人都可以告訴我我的代碼有什麼問題，或者只是給我一個提示？我真的需要幫助。這裏是我的公式，我需要繪製（使用Maketex書面）： 這是我的代碼： close all; clear;clc; Nr=2;N
    matlab plot numerical-methods axis-labels numerical-integration 2017-03-19

    1熱度

    1回答

    數值的Hessian使用梯度函數R

    我需要計算使用我的梯度函數（由式，不是數字編程）我的函數數值的的Hessian。像numDeriv或rootSolve這樣的軟件包使用數字計算黑森州不符合我的需求的梯度。我需要執行的內容內部（我可以調用沒有單獨的方法）在optim包，但唯一的方法來處理我的優化任務良好實施在nlopt包和傳遞其最佳值爲optim爲了得到粗麻布太昂貴我的程序。 所以我需要一些功能，使用非數字梯度（見例如這些公式htt
    r optimization numeric numerical-methods hessian 2017-02-05

    0熱度

    1回答

    在Python中獲取勒讓德多項式的導數

    我有一個引力勢的表達式（來自here的eq.15），並計算我需要的軌道以評估作爲局部梯度的引力，對於我來說意味着評估Legendre polynomialsP2，P4和P6的衍生物，其單值爲成千上萬次。 我可以使用this question表達式計算的話，但我不知道是否有辦法讓蟒蛇對於沒有明確涉及我評估衍生物作爲有限差的衍生物。 我無法在SciPy中找到任何要自動執行此操作的內容。在numpy.p
    python numpy scipy numerical-methods 2017-02-10

    0熱度

    1回答

    R中的digamma函數的根

    我正在研究最大似然估計器，其中一個參數是使用digamma函數估計的。我試圖用uniroot解決方程，但無法做到這一點。這裏是我的代碼： dig = function(alpha){ digamma(2 + alpha) - digamma(alpha) - (1/(2+alpha)) + (2/(2+alpha)) } curve(dig, from = 0, to = 10)
    r statistics numerical-methods mle gamma-function 2017-02-12

    0熱度

    2回答

    使用來自scipy.integrate的四元數值積分結果中的虛假陷波

    我試圖從數字上解決使用NumPy和quad從scipy.integrate以下積分。該代碼是有點兒工作，但我得到虛假的缺口中得到的數據： 任何人有他們爲什麼發生，以及如何得到正確的結果，平穩的想法？ 這裏是在Python原代碼： #!/usr/bin/env python import numpy as np from scipy.integrate import quad import
    python numpy scipy numerical-methods numerical-integration 2017-02-17

    1熱度

    1回答

    Lapack例程中的分割錯誤dstev

    我目前正在編寫一個小程序，用C++中的有限差分方法求解微分方程。問題是我正在使用lapack來解決我的三角對角矩陣，但是在例程被調用後我得到了一個Segmentation Fault。該例程仍然退出info = 0，所以我不知道問題是什麼。這是該計劃。 #include <iostream> #include <cmath> #include <vector> #define pi 3.1
    c++ vector numerical-methods lapack 2017-03-05
最新問題