numerical-integration

    2熱度

    1回答

    我很清楚如何以最簡單的形式使用scipy.integrate.ode.integrate(t)函數,但API讀取它也需要兩個可選參數,即step和relax。目前的documentation沒有關於這些參數的信息,也沒有在example中使用。我想知道,他們做了什麼以及哪些情況下他們有用?

    1熱度

    1回答

    任務是通過用數字方法創建可展曲面,在Matlab中創建一個錐頂。有3個部分我已經完成2.我的問題是關於第3部分,我需要計算可以包含帽子的最小長方形紙張表面。我需要計算紙張的材料浪費。 你也許可以跳過LONG背景並轉到最後一段 背景: 圓錐帽可以與用其尖端位於傾斜圓錐體被創建(一個; 0,B )和一個圓形的基地。 x = Rcos u, y = Rsin u z = 0 0<_ u >_2pi

    1熱度

    1回答

    我想用Python和Numpy來解決使用有限差分方法的基本方程。代碼給了我正確的第一個值,即它給了我[1];然而,之後的其他每個值都只是零? 我不知道我做錯了什麼,因爲它顯然適用於第一個值,所以我該如何解決這個問題? 任何想法都會非常有幫助。 from numpy import * import numpy as np import matplotlib.pyplot as plt impo

    -2熱度

    3回答

    我想做一些非常簡單的事情,但我無法在numpy的深處找到它。我想用數值和不斷整合由它的值給出的函數。這意味着我只需要一個包含輸入數組開頭總和的數組。例如: 輸入: [ 4, 3, 5, 8 ] 輸出: [ 4, 7, 12, 20 ] # [ sum(i[0:1]), sum(i[0:2]), sum(i[0:3]), sum(i[0:4]) ] 聽起來很直截了當,所以我希望這一定是很容易

    1熱度

    1回答

    我試圖用一個公式來擬合我的數據,該公式包含具有無限積分限值的定積分的數值計算。對於擬合,我使用要求矢量化模型函數的八度函數leasqr。以下代碼會生成錯誤,在調用數值積分時會發生這種錯誤。 不符合標準的參數(OP1是1x387,OP2是10X2) function [fGsAb] = GsAbs (x, p) Hw = 3108.0 ; fGsAb = Hw ./ (2.4 .* p(1) .

    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

    1熱度

    2回答

    我目前正在計算python中的數值積分,它實際上是四個嵌套數值積分。每個積分是不同變量的函數,但我遇到的問題是限制取決於其嵌套內積分(最外面的積分具有適當的浮點範圍)。例如,假設我有一個具有限制amin,amax,bmin,bmax,cmin,cmax,dmin,dmax的函數f(a,b,c,d)。 amin和amax是浮動的。 bmin和max是a,cmin和cmax的函數是b的函數,dmin和

    0熱度

    1回答

    我有一個程序運行ode15s幾千次,以找到一個特定的解決方案。但是,我越來越多的集成容忍錯誤,如下列: "Warning: Failure at t=5.144337e+02. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (1.818989

    1熱度

    2回答

    是否有可能meijerG函數包含負值(即{-1,0,0})?我嘗試了Mathematica和Matlab來計算這個meijerG函數,但是他們產生了一個錯誤,這個錯誤代碼是meijerG is not defined for the given parameters。任何幫助都會有所幫助。 這裏是我的代碼: D = (0.6); lg1 = lg2 = 1; G = evalin(symeng

    1熱度

    1回答

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