numerical-analysis相关内容
我有以下 Fortran 代码: 奇怪的程序真实(种类=8)::Pi1=3.1415926535897932384626433832795028841971693993751058209;真实(种类=8)::Pi2=3.1415926535897932384626433832795028841971693993751058209_8;打印*, "Pi1=", Pi1;打印*, "Pi2=", P
..
我在几个地方读过(在文档和这篇博文中: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实现一篇论文,该论文最初是在MATLAB中实现的.该论文说,使用曲线拟合从一组采样数据点中找到了一个五次多项式.我不想使用它们的多项式,因此我开始使用样本数据点(在纸上给出),并尝试使用sklearn多项式特征和linear_model查找5度多项式.因为它是一个多元方程f(x,y),其中x和y是某个池塘的长度和宽度,而f是污染物的初始浓度. 所以我的问题是sklearn
..
我希望执行以下操作: 1)尝试各种输入组合以搜索最佳结果 2)重置所有数组,使其与每次代码循环之前一样 我正在使用的每个变量都在一个数组中,例如f [0,1,2,3,...] 问题可能在于每个遍部分之后的重置变量,因为第一遍工作正常,但是第一遍的残差导致随后的迭代提前中断. 这是我的方法的伪.非常简单,可能是Python如何处理数据(面向对象)的问题. index_save
..
有人知道用于发现n×n矩阵的特征多项式的Souriau方法吗?我发现第一个系数很明显,但是如何找出其他系数呢?在我需要对矩阵求逆之后,我知道如何. #include #include using namespace std; double trace(double a[5][5],int n){ int i; d
..
我已经问过这个问题对任何方程组都很好?,从中我得到了满意的答案.我在那里展示的系统 x = A * exp(x + y) y = 4 * exp(x + y) 只是一个玩具模型,与我的真实案例问题类似,fsolve进行了处理(以下答案中的代码): from scipy.optimize import fsolve import matplotlib.pyplot as pl
..
在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的示例问题 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
..
我想计算(1.0-p)^n,其中p是0到1之间的双精度数(通常非常接近0),n是一个正整数,可能是几百或几千(也许更大;我是还不确定).如果可能的话,我想只使用Java内置的java.lang.Math.pow(1.0-p, n)来实现,但是我有点担心,这样做可能会导致我对所关注的值范围产生巨大的准确性/精度损失.有人对使用Java的实现会有什么样的错误有一个大概的认识吗?我不确定它们的实现到底是
..
是否存在可以执行以下操作的稀疏矩阵库: 求解线性代数方程 支持诸如矩阵矩阵/数乘法/加法/减法,矩阵转置,获取矩阵的行/列等操作 矩阵大小可以为40k * 40k或更大,例如250k * 250k 快速 可以在Windows中使用 有人可以为我推荐一些图书馆吗? 如果您推荐,请告诉我它的优缺点,以及推荐它的原因. 顺便说一句,我已经在互联网上搜索了许多稀疏矩阵库并对其进行
..
我注意到,如果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
..
我正在尝试在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
..
我在一些地方(在文档和此博客文章中阅读过: 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库来做这个? I有点没有成功地搜索(例如,谷歌搜索“功率谱密度java”或“信号处理java”和点击链接,查看Apache Commons,Sourceforge,java.net等)。 有许多小程序,书籍,教程,商业产品等,不符合我的需求。 更新 :我找到了 o
..