-
0 votesanswersviews
用python中的积分方程拟合数据
我有一些数据,我试图适应包含和定积分方程的模型 . 我的策略是使用optimize.leastsq和integrate.quad,我不断收到类型错误:“只有length-1数组可以转换为Python标量” 任何帮助将不胜感激 . 这是我的代码的相关部分(请记住self.vvals和self.bvals是一维数组,self.L是一个浮点数): def NLFit(self): '''fits... -
0 votesanswersviews
在Scipy LeastSq中 - 如何添加惩罚词
如果对象函数是 如何在python中编码?我已编码正常的: import numpy as np import scipy as sp from scipy.optimize import leastsq import pylab as pl m = 9 #the degree of the polynomial def real_func(... -
0 votesanswersviews
如何用约束求解这个线性方程?
我有a1,a2,a3,b1,b2,b3作为已知值,所以我可以列出如下的两个方程式 . 但是,这还不足以让H解决 . 所以,还有额外的约束|| H || = 1 . matrix H = [x y z]^T [a1, a2, a3] H = 0 [b1, b2, b3] H = 0 我想问一下如何解决这个问题? -
4 votesanswersviews
最小的置信区间适合scipy python
如何在python中计算最小二乘拟合(scipy.optimize.leastsq)的置信区间? -
1 votesanswersviews
Scipy LeastSq错误栏
我使用来自SciPy的LeastSq将实验光谱拟合到理论预期 . 当然有与实验值相关的错误 . 如何将这些提供给LeastSq或者我需要一个不同的例程?我在文档中找不到任何内容 . -
1 votesanswersviews
python scipy leastsq适合复杂的数字
我有一个复数的数据集,我希望能够找到最适合数据的参数 . 您可以使用pyp中的scipy实现的leastsq来复杂数字的数据吗? 例如,我的代码是这样的: import cmath from scipy.optimize import leastsq def residuals(p,y,x): L,Rs,R1,C=p denominator=1+(x**2)*(C**2... -
4 votesanswersviews
Scipy leastsq()函数开销
我正在研究一个图像分析程序,我已经缩小了我的瓶颈,试图将2D高斯镜片多次安装到一个小窗口(20x20)像素上 . 在此代码中花费了90%的执行时间 . 我正在使用scipy cookbook中给出的代码来解决这个问题: def gaussian(height, center_x, center_y, width_x, width_y): """R... -
0 votesanswersviews
具有权重的scipy(leastsq)中的NonLinearModelFit
我非常喜欢scipy和我一起请:-) 我最近一直在使用mathematica来搞乱我的数据 . 我有一种计算来自静态接收器(也是x,y坐标)的4个或更多距离测量的x,y坐标的方法 . 我用来最有效地使用我的数据执行此操作的函数是mathematica函数: NonlinearModelFit[data, Norm[{x, y} - {x0, y0}], {x0, y0}, {x, y}, ... -
0 votesanswersviews
SciPy LeastSq未能提供最佳解决方案
我试图使用Leastsq来拟合一条非常简单的曲线 . 但是,其解决方案尚未优化 . 有人能给我一些建议吗?以下是我的代码: from scipy import optimize import numpy as np hl_obs = np.array([10.0, 23.0, 20.0]) ph=np.array([5.0,7.0,9.0]) tp=60 def residuals(k_ab... -
0 votesanswersviews
SciPy Leastsq实施问题
我使用下面的代码来实现我的模型与一些实验数据的最小平方拟合 . 该程序只是继续运行,似乎我的实现有问题,但我无法搞清楚 . 两个函数中的打印语句(simulateFunc和fit)表示这些语句被调用了几次(准确地说是9),但之后没有响应 . 谁能告诉我我的实施有什么问题?程序运行正常,没有使用minimalsq函数,我可以通过调用带有一些参数的fit函数来获取残差 . 示例数据点:(我的文件中共有... -
0 votesanswersviews
scipy leastsq fit - 惩罚某些解决方案
我已经实现了一种能够同时适应多个数据集的算法 . 它基于这个解决方案:multi fit 目标函数太复杂而无法显示(LaFortune散点模型),因此我将使用解决方案中的目标函数进行解释: def lor_func(x,c,par): a,b,d=par return a/((x-c)**2+b**2) 如果选择导致 lor_func < 0 的参数集 par ,如何惩罚拟... -
6 votesanswersviews
python中稀疏矩阵的伪逆
我正在使用来自神经成像的数据,并且由于数据量很大,我想为我的代码使用稀疏矩阵(scipy.sparse.lil_matrix或csr_matrix) . 特别是,我需要计算矩阵的伪逆以解决最小二乘问题 . 我找到了sparse.lsqr方法,但效率不高 . 有没有一种方法来计算Moore-Penrose的伪逆(对应于正常矩阵的pinv) . 我的矩阵A的大小约为600'000x2000,并且在矩阵... -
0 votesanswersviews
手动计算多元线性回归(OLS)的置信区间
我试图了解如何手动计算多元线性回归(OLS)的置信区间 . 我的问题是我不知道如何计算所有单个系数的标准误差 . 对于只有一个自变量的回归,我遵循以下教程:http://stattrek.com/regression/slope-confidence-interval.aspx . 本教程提供以下公式: 事实证明,该公式有效 . 但是,我并没有完全理解这个公式 . 例如,为什么(-2)位于公式的... -
0 votesanswersviews
内存泄漏使用numpy linalg并行
我有以下代码在矩阵U的每一行上运行最小二乘法 . import numpy as np from numpy.linalg import norm,lstsq from sklearn.externals.joblib import Parallel,delayed k = 25 max_iter = 50 lmbda = 1e-3 R = np.random.randint(0,2,(1500,... -
2 votesanswersviews
最小二乘拟合python中具有错误的复数
这里抱歉是代码 . 我试图将贝塞尔函数拟合到某些2D数据,贝塞尔函数很复杂 . 但是,当我适应数据时,我总是最终得到一个数组错误,我无法排序(我是python的新手) . 最后,代码理想情况下会给出错误来解释导入数据的分散 . from numpy import * from pylab import * from scitools.std import * import csv import m... -
0 votesanswersviews
解决非负最小二乘问题“xA = b”[重复]
这个问题在这里已有答案: Constrained linear least-squares for xA=b in matlab 1回答 我想找到“xA = b”的非负最小二乘解 . 我很高兴获得Python,Matlab或R的答案 . A 是6 * 10矩阵, b 是8192 * 10矩阵 . 我发现了一些函数:Python中的 least_squares 和 nnls ,以及Matlab中... -
1 votesanswersviews
非线性拟合
我有一些实验数据,我想使用最小二乘法(Levenberg-Marquardt)来拟合我的参数 . 我正在使用两个非线性方程式,我正在使用一些计算程序(Origin和Matlab) . 首先是: y=A+B*(((2*pi*x)^2+Alfa4^2)*((2*pi*x)^2+Alfa5^2))/(((2*pi*x)^2+Alfa1^2)*((2*pi*x)^2+Alfa2^2)*((2*pi*x)^... -
2 votesanswersviews
稀疏线性约束非线性最小二乘的求解器
有没有任何算法或求解器可以解决非线性最小二乘问题,其中已知jacobian总是稀疏的,并且解决方案受到以下任一方面的约束: 线性等式约束和框约束 一般线性不等式约束 -
1 votesanswersviews
在python中使用复杂函数拟合数据
我有一些数据,我试图适应模型 这是我的代码的相关部分 path='D:/ParPhy/2-BESIII15_new.dat' data = pd.read_table(path,header=None) y=np.array(data[1]) x=np.array(data[0]**(1/2)) s=x**2 def F_w(s,alpha,m_p,gamma_p): P_s=1+alph...