numerical-methods

0热度

1回答

我正在尝试使用高斯 - 牛顿方法在Python中找到近似交点n。这里的目标是当给定点（x1，y1），（x2，y2），...，（xn，yn）以及半径R1，R2，...，Rn，高斯 - 牛顿法是用于找到最小化到圆的距离的平方和的点的最小化。这是一种迭代方法，我们从一个初始向量v0 = [0,0]开始。第一次迭代完全正确，但后续迭代不正确。我不能在代码发现错误： solutions = [] rms

1热度

1回答

在包含未知变量的Matlab中求解积分？

任务是通过用数字方法创建可展曲面，在Matlab中创建一个锥顶。有3个部分我已经完成2.我的问题是关于第3部分，我需要计算可以包含帽子的最小长方形纸张表面。我需要计算纸张的材料浪费。你也许可以跳过LONG背景并转到最后一段背景：圆锥帽可以与用其尖端位于倾斜圆锥体被创建（一个; 0，B ）和一个圆形的基地。 x = Rcos u, y = Rsin u z = 0 0<_ u >_2pi

0热度

1回答

用python解决二阶ODE

我在写一些代码来求解二阶微分方程，但是它给出了一个完全错误的结果。我认为问题在于以正确的方式表达欧拉方法。我也尝试了另一个二阶ODE，但我也未能近似y（x）。你能指出哪里可能出错吗？请参阅图形和代码：求解ODE： y"(x)=-1001y'(x)-1000y(x) y(0)=1, y'(0)=0 解析解： y(x)=(1000exp(-x)-exp(-1000x))/999 重写为2

0热度

1回答

如何在CSR格式中表示对称矩阵？

尽管我使用C语言，但我相信编程语言和代码是无关紧要的。我想表示CSR格式的对称矩阵，用于快速矩阵向量乘积。我读取的输入流只给出矩阵的上三角部分，按数字编号，从左到右，按行排列。实施例： 0 2 3 5 2 0 1 4 3 1 0 9 5 4 9 0 输入流： 0 2 3 5 0 1 4 0 9 0 如果我建立，比如说，从二维阵列的规则的矩阵，这将是容易的。我需要的是我得到的数据

0热度

1回答

Swift 3中的Lotka-Volterra模型

我正在尝试实现一个rk4函数来求解2个微分方程。我有这样的代码，实现了龙格库塔4方法： //RK4 method func rk4_func(y_array: [Double], f_array: [(([Double], Double) -> Double)], t_val: Double, h_val: Double) -> [Double] { let length

2热度

2回答

定义一个函数用于从数据文件求解方程

我对Python和实际上任何基本的编程语言都是全新的，我将Mathematica用于我所有的符号和数字计算。我正在学习使用python并发现它非常棒！这是我正在努力解决的一个问题，但没有任何线索而陷入困境！我有一个数据文件例如 0. 1. 0.01 0.9998000066665778 0.02 0.9992001066609779 ... .. 哪些只是{吨，COS [2吨

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

0热度

1回答

最小化过程快速减少然后增加

我运行MATLAB的FMINCON的“内点”算法。该算法使用BFGS，因此考虑到数值Hessian中的信息。我正在优化一个参数空间，这个参数空间就一个独特的解决方案而言是非凸的。我惩罚那些不是arctan函数唯一的参数组合。我发现优化过程开始时高（约10），在几次迭代中减小到4（接近最小值），然后在几百次迭代中缓慢增加回10。我有选项设置重新调整问题和约束。我有点卡在这一点。你怎么看？

0热度

2回答

'pdepe'（MATLAB函数）使用什么样的数值方法？

我正在使用MATLAB的函数'pdepe'来解决一些偏微分方程（抛物线方程）的问题。我需要知道函数使用的数值方法的种类，因为我必须在报告中通知它。 MathWorks中对函数的描述是“在一个空间变量和时间内解决抛物型和椭圆型偏微分方程系统的初始边值问题”。它是一种有限差分方法吗？感谢您的帮助。

0热度

1回答

在MATLAB中的3D拉普拉斯松弛

我能写出一个MATLAB程序，绘制两个金属板之间的1D拉普拉斯松弛，以使用Jacobi方法找到平衡电位。我会说实话，我不能完全肯定，我明白我做了什么，但这里是代码的相关部分： N=100; V = zeros(N,N); V(1,:)=-1; V(N,:)=1; for n=1:400 for i=2:99 for j=2:99 V(i,j)=(V(