numerical-analysis

0熱度

1回答

我正在用Python實現一篇文章，最初是在MATLAB中實現的。該論文指出，使用來自一組採樣數據點的曲線擬合找到五次多項式。我不想使用它們的多項式，所以我開始使用樣本數據點（在紙中給出），並嘗試使用sklearn多項式特徵和linear_model來查找5度多項式。因爲它是一個多元方程f（x，y），其中x和y是某個池塘的長度和寬度，f是污染物的初始濃度。所以我的問題是，sklearn多項式特徵將

-2熱度

1回答

計算和確定收斂速度

寫MATLAB代碼來計算和確定的收斂速度： (exp(h)-(1+h+1/2*h^2))/h ‍‍‍‍‍‍ ‍‍‍‍‍‍ with h=1/2, 1/2^2,..., 1/2^10 我的代碼是： h0=(0.5)^i; TOL=10^(-8); N=10; i=1; flag=0; table=zeros(30,1); table(1)=h0 while i < N h=

0熱度

1回答

找到根確定區間的個數

我有n度多項式系統，只是我想了解之前確定的區間之間的根數。但我不想找到root.I需要root的數量。我需要編寫python代碼。例如： x^8 + 2x^6 + 5x^4 + x^2 + 45x + 1 = 0 我們在3-5之間有多少根？強調=我不想找到根，只是我想了解我有多少根。

0熱度

1回答

用apache-commons-math計算數組的派生值

早上好，我有一個數組，包含大約3000個double值，我需要找到所有局部最小值和最大值，爲此我對一階和二階導數感興趣，用Apache Commons Math實現這個最好的方法是什麼？我的麻煩是我直接從數組開始，而不是從像sin（x）這樣的函數開始。謝謝

1熱度

2回答

MATLAB - 自適應步長Runge-Kutta

我在MATLAB中編寫了一個自適應步長RK4來解決一個ODE系統。代碼運行時沒有錯誤，但是當我嘗試對y繪製x時，它不會生成所需的曲線。而不是一個環形的形狀，我只是得到一條平坦的線。這從r正在輸出一個常數值的事實可以看出。在檢查每行的輸出後，它們不輸出常量或錯誤或inf或NaN，而是輸出實數和虛數分量（複數）。我不知道爲什麼會發生這種情況，我相信這是我麻煩的根源。 function AdaptRK4

1熱度

1回答

波動方程FDM，MATLAB

我想實現MATLAB代碼解決了波動方程，我的功能看起來是這樣的： function [x,t,w] = wave_eqn(xl,xr,yb,yt,M,N,f,l,r,p) % input: space interval [xl,xr], time interval [yb,yt] % number of space steps M, number of time steps N % outp

1熱度

1回答

二階法比一階法差嗎？

我在想數值分析中的一個基本問題。由於二階方法的截斷誤差爲O（dx^2）和O（dx）爲（），因此當離散常微分方程時，衆所周知二階方法比一階方法更精確。一階方法。 0 < dx < 1. 如果dx> 1，該怎麼辦？例如，域爲0到10000，網格大小爲1000，則dx = 10。在這種情況下，由於dx^2 = 100且dx = 10，所以二階方法不像一階方法那樣精確。在處理大規模問題時，我們可以遇到這個

0熱度

1回答

Romberg積分算法

我想解決在給定區間[a，b]上的基本積分，對於Romberg積分中的R（n，m）中的任意n和m值。我從梯形法則中推導出布爾的規則，所以我知道如何在紙上做到這一點。我甚至畫出了一張顯示所有依賴關係的流程圖。這不是幫助我編碼。我有一種感覺，這需要某種遞歸。我正在用Java編程。編輯：我沒有問任何人爲我代碼。上述所有信息，包括其他人編輯的信息，都是爲了說明我的相關經驗和對技術詞彙的理解以及我應用

-4熱度

1回答

python中的數值分析矩陣圖像

我可以用矩陣圖像做些什麼？ import matplotlib.pyplot as plt import matplotlib.image as mpimg img=mpimg.imread('image.png') imgplot = plt.imshow(img) 我已經想到：相同的圖像之間 -difference注意到一些運動 -Transpose矩陣旋轉圖像有沒有其他的操作，我

-1熱度

1回答

數值方法runge kutta

求解微分方程。 dy/dx = xz + 1，對於x = 0.3，dz/dx = -xy，假設y = 0且z = 1，當x = 0時使用四階runge kutta方法。我無法啓動此問題。請幫忙。我需要使用數值分析而不是編程的數學解決方案。