Fortran 95 中的数值精度:
我有以下 Fortran 代码: 奇怪的程序真实(种类=8)::Pi1=3.1415926535897932384626433832795028841971693993751058209;真实(种类=8)::Pi2=3.1415926535897932384626433832795028841971693993751058209_8;打印*, "Pi1=", Pi1;打印*, "Pi2=", P
..
numerical-analysis相关内容
为什么 Matlab 的 inv 慢且不准确?
我在几个地方读过(在文档和这篇博文中:http://blogs.mathworks.com/loren/2007/05/16/purpose-of-inv/ ) 不推荐在 Matlab 中使用 inv，因为它速度慢且不准确. 我正在尝试找出这种不准确的原因.到目前为止，谷歌没有给出有趣的结果，所以我认为这里有人可以指导我. 谢谢！ 解决方案 我提到的不准确在于方法 INV，而不
..
求和除以给定数的数组的子数组数
我被一个算法问题困住了.请为我提出一些针对以下问题的有效算法. 问题是 求和能被给定数整除的子数组的个数. 我的工作 我做了一个算法，它的复杂度是 O(N^2)，这里，N = 数组的大小. 我的代码 #include使用命名空间标准；主要的() {整数 N;内部 P;国际T;内部值；long long int count = 0;long long int answ
..
正确舍入的双精度除法
我正在使用以下算法进行双精度除法，并尝试使其在浮点软件仿真中正确取整.假设 a 为被除数，而 b 为除数. 所有操作均在Q2.62中执行. 倒数的初始近似值为 近似为 0x17504f333f9de6 (在Q2.62中为 0x5D413CCCFE779800 ). 之后，用牛顿-拉夫森迭代来近似倒数: 最后的舍入是: 如果a
..
Python中多元5度多项式回归的曲面图
我正在用Python实现一篇论文，该论文最初是在MATLAB中实现的.该论文说，使用曲线拟合从一组采样数据点中找到了一个五次多项式.我不想使用它们的多项式，因此我开始使用样本数据点(在纸上给出)，并尝试使用sklearn多项式特征和linear_model查找5度多项式.因为它是一个多元方程f(x，y)，其中x和y是某个池塘的长度和宽度，而f是污染物的初始浓度. 所以我的问题是sklearn
..
如何使用嵌套的for循环进行强行强制组合-Python
我希望执行以下操作: 1)尝试各种输入组合以搜索最佳结果 2)重置所有数组，使其与每次代码循环之前一样 我正在使用的每个变量都在一个数组中，例如f [0,1,2,3，...] 问题可能在于每个遍部分之后的重置变量，因为第一遍工作正常，但是第一遍的残差导致随后的迭代提前中断. 这是我的方法的伪.非常简单，可能是Python如何处理数据(面向对象)的问题. index_save
..
特征多项式的Souriau方法
有人知道用于发现n×n矩阵的特征多项式的Souriau方法吗?我发现第一个系数很明显，但是如何找出其他系数呢?在我需要对矩阵求逆之后，我知道如何. #include #include using namespace std; double trace(double a[5][5],int n){ int i; d
..
对于方程组，fsolve，brentq和root的用法和精度有何不同?
我已经问过这个问题对任何方程组都很好?，从中我得到了满意的答案.我在那里展示的系统 x = A * exp(x + y) y = 4 * exp(x + y) 只是一个玩具模型，与我的真实案例问题类似，fsolve进行了处理(以下答案中的代码): from scipy.optimize import fsolve import matplotlib.pyplot as pl
..
Modelica评估单
在Modelica规范中我真的找不到任何答案，请病态的你们问.规范指出 A tool is free to solve equations, reorder expressions and to not evaluate expressions if their values do not influence the result (e.g. short-circuit evaluation of
..
在傅里叶插值的情况下返回函数的“传统"符号
在数值分析中，我们必须让学生在R中实现代码，即在给定函数f(x)找到其傅立叶插值tN(x)并计算出插值误差的情况下 $||f(x)-t_{N}(x)||=\int_{0}^{2\pi}$ $|f(x)-t_{N}(x)|^2$ 或各种不同的$ N $ 我首先尝试根据以下公式计算d系数: $d = \frac 1N M y$ 其中M表示DFT矩阵，y表示一系列等距函数值，
..
如何找到特定机器上的尾数长度?
我想找到特定计算机上的尾数位数和四舍五入单位.我了解它们是什么，但不知道如何找到它们-尽管我知道它们可能因计算机而异. 我需要这个数字才能执行数值分析的某些方面，例如分析错误. 我目前在想的是，我可以编写一个小型的c ++程序来逐渐增加一个数字，直到出现溢出为止，但是我不确定要使用哪种数字. 我在正确的轨道上吗?究竟该如何计算呢? 解决方案 我认为您使用的任何语言都将指定
..
如何通过Runge-Kutta 4传递硬编码微分方程
我正在尝试实现Runge-Kutta的示例问题 dy/dt = y-t ^ 2 +1和dy/dt = t * y + t ^ 3在C#中，我似乎无法获得期望的输出.我已将程序分为几个类，以尝试分别查看工作.我认为我的主要错误来自尝试使用委托将方法作为变量传递给Runge-Kutta进程. 方程类: namespace RK4 { public class Eqn {
..
用定点迭代法求解该方程
我怎么解决这个方程式 x 3 + x-1 = 0 使用定点迭代吗? 我可以在网上找到任何定点迭代代码(尤其是在Python中)吗? 解决方案 使用import scipy.optimize as optimize def func(x): return -x**3+1 # This finds the value of x such that func(x) = x, t
..
大n的java.lang.Math.pow(x，n)有多精确?
我想计算(1.0-p)^n，其中p是0到1之间的双精度数(通常非常接近0)，n是一个正整数，可能是几百或几千(也许更大；我是还不确定).如果可能的话，我想只使用Java内置的java.lang.Math.pow(1.0-p, n)来实现，但是我有点担心，这样做可能会导致我对所关注的值范围产生巨大的准确性/精度损失.有人对使用Java的实现会有什么样的错误有一个大概的认识吗?我不确定它们的实现到底是
..
C ++的稀疏矩阵库
是否存在可以执行以下操作的稀疏矩阵库: 求解线性代数方程 支持诸如矩阵矩阵/数乘法/加法/减法，矩阵转置，获取矩阵的行/列等操作 矩阵大小可以为40k * 40k或更大，例如250k * 250k 快速 可以在Windows中使用 有人可以为我推荐一些图书馆吗? 如果您推荐，请告诉我它的优缺点，以及推荐它的原因. 顺便说一句，我已经在互联网上搜索了许多稀疏矩阵库并对其进行
..
为什么在Matlab和Octave中inv()和pinv()的输出不相等?
我注意到，如果A是一个NxN矩阵，并且它具有逆矩阵.但是inv()和pinv()函数的输出不同. -我的环境是Win7x64 SP1，Matlab R2012a，Cygwin Octave 3.6.4，FreeMat 4.2 看看Octave中的示例: A = rand(3,3) A = 0.185987 0.192125 0.046346 0.140710 0.35100
..
Jacobi迭代不会结束
我正在尝试在MATLAB中实现Jacobi迭代，但是无法使其收敛.我已经在网上和其他地方查看了可以比较的工作代码，但是找不到与我的代码相似的东西，并且仍然可以使用.这是我所拥有的: function x = Jacobi(A,b,tol,maxiter) n = size(A,1); xp = zeros(n,1); x = zeros(n,1); k=0; % number of s
..
为什么Matlab的inv缓慢而又不准确?
我在一些地方(在文档和此博客文章中阅读过: http://blogs.mathworks.com/loren/2007/05/16/purpose-of-inv/)，不建议在Matlab中使用inv，因为它速度慢且不准确 我正在尝试找出这种不准确的原因.截至目前，Google尚未给出有趣的结果，所以我认为这里有人可以指导我. 谢谢！ 解决方案 我提到的不准确之处在于方法INV，
..
绑定此程序以确定不包含零的倒数整数之和
让 A 表示正整数的集合，其十进制表示不包含数字0. 中元素的倒数之和已知 为23.10345。 Ex。 1,2,3,4,5,6,7,8,9,11-19,21-29,31-39,41-49,51-59,61-69,71-79,81-89， 91-99,111-119，... 然后取每个数字的倒数，并计算总数。 如何通过数字验证？ 编写一个计算机程序来验证这个数字。 这
..
Java中的信号处理库？
我想计算时间序列的功率谱密度;做一些带通，低通和高通滤波;也许还有一些基本的东西。 有没有一个很好的开源Java库来做这个？ I有点没有成功地搜索（例如，谷歌搜索“功率谱密度java”或“信号处理java”和点击链接，查看Apache Commons，Sourceforge，java.net等）。 有许多小程序，书籍，教程，商业产品等，不符合我的需求。 更新 ：我找到了 o
..