numerical-methods

0熱度

1回答

我想用C++找到一個具有簡單定點方法的根，但問題是Xr是f(x)的一個根，也是一個拐點。另外，A方程比正常的定點方法稍微複雜一些。該方程式增加了常數c以檢查多快地收斂到根xr。我打算找到一個根，然後檢查根是否是一個拐點，但它不工作，我在代碼中找不到問題。我需要你的幫助。真正的問題是考慮求根問題f(x)=0與根xr，其中f '(x)=0。將其轉換爲簡單的定點問題。 x=x+c*f(x)=

4熱度

1回答

matlab中fft，ifft和fftshift的使用

我試圖實現split-step fourier方法來求解光學中的非線性薛定諤方程。它基本上分別對待線性部分和非線性部分。它利用傅里葉變換和時域非線性部分解決了線性部分。下面的代碼從一本書複製： alpha = 0 beta_2 = 1 gamma = 1 T = linspace(-5,5,2^13); delta_T = T(2)-T(1); L = max(size(A));

0熱度

1回答

離散微分方程並在每個時間步中找到參數的值

我想通過在t上使用for循環來找到f(t)的值。管理f的公式爲： df/dt = cos(f) g(x,y,t) x和y在二維座標是。我寫了下面的代碼： dt=1; a=2; f(1,1)= a*rand(1,1); for t = 2:100 f(1,t)= cos(f(1,t-1)) g(x,y,t) dt; end 但我不確定此方法是否正確，或者如果可以寫得更好。任何

-2熱度

2回答

如何生成在Python所有點爲0 <X_1 + X_2 + X_3 <1, x_1> 0，X_2> 0，X_3> 0，說DX = 0.01

如何生成的所有點的Python爲0<x_1+x_2+x_3<1，x_1>0，x_2>0，x_3>0，說$ DX = $ 0.01。嘗試在解決約束優化問題時將其用作網格。

1熱度

2回答

不知道爲什麼我得到這個錯誤，我有一個想法，但不知道如何解決它

這是爲我正在處理的項目。我正在試圖在地殼內的一些主體或鄉村岩石中模擬岩漿冷卻的侵入堤防。我對編碼相當陌生。我盡我所能將這種代碼格式從另一種編碼語言轉換爲python。我對發生了什麼有一個基本的想法。我知道我正在嘗試索引超出範圍的內容，但我不確定在哪裏以及如何解決它。任何幫助，我可以得到將不勝感激！提前致謝。 import numpy as np import matplotlib.pyplot

0熱度

1回答

使用Python來計算m階下降階乘

下降階乘：https://en.wikipedia.org/wiki/Falling_and_rising_factorials 如何使用Python來計算呢？我嘗試 def factorial(gamma, m): prod = 1 for i in range(m): prod *= (1-gamma-i)

3熱度

2回答

將我的有限差分模型從2D更改爲3D會導致不穩定行爲

我一直在寫一個有限差分代碼，用於使用激光誘導熱成像進行模擬和裂紋檢測。裂縫由因子a和b實現，這些因子通過使用鬼點方法「阻尼」通過充氣裂縫的熱流。二維模型按預期運行，穩定條件滿足，一切正常。它甚至可以用實驗數據證明。只需複製並粘貼即可使用。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 2

1熱度

1回答

矩陣求逆方法

當一個具有與載體基質逆乘法的一個問題，因爲這樣的：一個可利用一個Cholesky分解甲和backsubstitute b找到產生的矢量x。但是，如果不像上面那樣制定問題，有時需要矩陣求逆。我的問題是處理這種情況的最佳方式是什麼。下面，我比較了各種方式（使用numpy的）反轉正定矩陣：首先，生成矩陣： >>> A = np.random.rand(5,5) >>> A array([[ 0.

7熱度

1回答

集成一個不直接操作元素的函數

早上好/下午/晚上，我正在研究Matlab腳本，涉及四階張量計算量積分。設H（r，theta，phi）是我想要整合的函數。假設H不能通過r，theta和phi的簡單操作獲得。我的問題是，在Matlab中任何其他代碼我知道： All input functions must accept arrays and operate elementwise. The function FUN(X,Y,Z

-1熱度

3回答

整數不溢出

如何計算hypot兩個整數，每個小於2^63，這樣在任何中間計算都不會溢出64位？（如傳統方法中的x^2+y^2）。鏈接的文章提到了一個浮點算法，由於整數爲0，所以不可以使用，因爲它是t = t/x;。我能找到的最接近的算法是從here但不幸的是它不夠精確： int ihypot(xd1, yd1) double xd1, yd1; { register x1 = (