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回答

用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矩阵旋转图像有没有其他的操作，我

5热度

4回答

FormPanel中的ExtJS FormPanel失败，并显示“this.body is null”

我有一个相当复杂的设置，我已经归结为下面的代码。我有一个外部FormPanel，其中我试图包括一个子类FormPanel的组件。在FF中，它导致“this.body is null”错误。这是否发生在其他人身上？有没有可能让这个工作？如果我不需要，我非常不想去触及这个子类。 var test = new Ext.Window({ title: 'test', items: [

-1热度

1回答

数值方法runge kutta

求解微分方程。 dy/dx = xz + 1，对于x = 0.3，dz/dx = -xy，假设y = 0且z = 1，当x = 0时使用四阶runge kutta方法。我无法启动此问题。请帮忙。我需要使用数值分析而不是编程的数学解决方案。