-
7 votesanswersviews
如何在Python中进行非线性复杂根发现
我想对以下非线性方程进行根搜索,我在Python中进行,但它不起作用 . 我的代码如下 from pylab import * import scipy import scipy.optimize def z1(x,y): temp=1+1j+x+2*y; return temp def z2(x,y): temp=-1j-2*x+sqrt(3)*y; retu... -
0 votesanswersviews
模数表达式算术方程
我如何解决:(2a ^ 2-195)mod26 = 1 我尝试了下一个方法:x =(2a ^ 2-195)并且如果x mod26 = 1则x = 27,53,79,105 ....但是找不到答案,我怎样才能在数学上解决这个问题? 谢谢! -
2 votesanswersviews
实时神经网络
In short: 是否可以让神经网络实时响应用户输入? 想象一下这个世界的小游戏,其中包含接收周围环境输入的实体,并使用神经网络生成输出以实现生存 . 这些实体应该有某种能力来杀死和消耗另一个实体,以延长饥饿倒计时,最终导致死亡,从而尽可能长时间地生存 . 一个简单的解决方案是使用遗传算法来改进神经网络并找到一组具有更高生存能力的实体(适用于这种情况) . 现在,如果用户应该能够控制一个这样的实... -
0 votesanswersviews
如何用数学方程绘制图片
最近我看到一些图片可以通过一些数学方程式绘制,如Batman Logo和Heart . 有没有特定的方法来找到绘制所需图像的方程? (例如,我想用一些数学方程式绘制字母S) . 谢谢 . 附:我想这是一个优化问题 . 首先从所需图片的边界获取一些样本,然后找到与这些样本具有最小差异的函数 . -
363 votesanswersviews
什么是计算机科学的NP-complete?
什么是NP完全问题?为什么它是计算机科学中如此重要的话题? -
0 votesanswersviews
如何解决此CPLEX异常?
我已经制作了一个程序来解决CPLEX的(非平凡)优化问题 . 我已经多次审查了代码和模型,我没有看到问题 . 但是,当我启动该过程时,会出现以下错误消息: Concert exception'ilog.cplex.CpxException:CPLEX Error 5002:'q1'中的Q不是正半正定 这是否意味着我的东西中存在错误或是否有诀窍可以避免这种情况? 提前致谢 ! -
3 votesanswersviews
BFGS无法融合
我正在研究的模型是一个多项logit选择模型 . 这是一个非常具体的数据集,因此其他现有的MNLogit库不适合我的数据 . 所以基本上,它是一个非常复杂的函数,它接受11个参数并返回对数似然值 . 然后我需要找到可以使用scipy.optimize.minimize最小化对数似然的最佳参数值 . 以下是我遇到的不同方法的问题: 'Nelder-Mead':效果很好,总能给我正确答案 . 但是,它... -
0 votesanswersviews
遗传算法立即在非常低的最佳值上收敛并且非常缓慢
我已经开始高尚的努力,让我工作的公司(有些)更高层,至少考虑ML . 我在汽车行业工作,所以我们的“问题”主要与找到汽车重心与后轮和前轮之间的完美距离,理想重量等等有关 . 因此,整个交易可以被制定为适合优化的适应度函数 . 实际功能的工作方式如下:我们得到了许多测试用例,比如10-15,代表100转/小时转弯的测试用例,150公里/小时的车道切换等等 . 我们有一个通用的测试功能,给定一组或多个... -
3 votesanswersviews
使用给定(错误的)节点距离优化图形的布局
我有一个松散连接的图表 . 对于该图中的每个边,我知道位置p(v)和p(w)处的节点v和w之间的近似距离d(v,w)为 a vector in R3 ,而不仅仅是欧氏距离 . 错误应该很小(比如说<3%),第一个节点是<0,0,0> . 如果根本没有错误,我可以这样计算节点位置: set p(first_node) = <0,0,0> calculate_positi... -
2 votesanswersviews
GA包中的“缺少值需要TRUE / FALSE”错误
当使用GA包运行遗传算法时,我收到"missing value where TRUE/FALSE needed"错误 . 在下面的代码中, like 是适应度函数, X 是长度为6的数值向量 . library(GA) y<-c(46,38,49,55) W<-matrix(c(0,1,0,0,1,0,0,0,0,0,0,1,0,0,1,0),ncol=4) ei&... -
16 votesanswersviews
Python中的约束线性回归
我有一个形式的classic linear回归问题: y = X b 其中 y 是响应向量 X 是输入变量的矩阵, b 是我正在搜索的拟合参数的向量 . Python提供 b = numpy.linalg.lstsq( X , y ) 来解决此形式的问题 . 但是,当我使用它时,我倾向于为 b 的组件获得极大或极小的值 . 我想执行相同的拟合,但将 b 的值限制在0到255之间 . 看起来 sci... -
2 votesanswersviews
尽管数学上是不可能的,但Gurobi报告了无界模型
我正在使用Julia的精彩JuMP包来解决一个使用Gurobi 6.0.4作为求解器的线性程序 . 目标函数是决策变量的总和,明确定义为非负,并且问题需要将其最小化 . 出于某种原因,Gurobi认为该模型是无限的 . 以下是变量和目标的定义: @defVar(model, delta2[i=irange,j=pair[i]] >= 0) @setObjective(model, Min, ... -
1 votesanswersviews
我应该使用gurobi线性化或尝试用python解决MINLP,还是尝试完全不同的方法?
我对此很陌生,所以我只是想拍摄并希望我尽可能准确,你会认为它值得回答 . 我正在尝试优化(最小化)成本/数量模型,其中两者都是连续变量 . 应尽量减少全球成本,但这取决于总量,这取决于具体成本 . 到目前为止我的代码看起来像这样: # create model m = Model('Szenario1') # create variables X_WP = {} X_IWP = {} P_W... -
0 votesanswersviews
数学问题:三角形,金字塔,旋转,平移,缩放
另一个棘手的问题 . 你在这里可以看到的是我的物理金字塔,它由3个LED组成,在1个平面上形成一个三角形,另一个在中间的中心,在另一个上方约18mm处 . 第四个使三角形成为金字塔 . (如果你看一下正确的三角形,你可能会看到它更好 . 这一个绕水平方向旋转,你可以很好地看到一个二极管上的二极管) . 第二张图显示了我正在运行的程序 . 左侧框显示LED的原始图片(带有红外滤镜的照片) . 中... -
-2 votesanswersviews
求解scipy中的非线性方程组,为什么它说“输入和输出形状之间不匹配”
我想解决一个非线性方程组(我有17个变量和17个方程) . 我使用了scipy.optimize.fsolve,但错误信息显示为“fsolve: there is a mismatch between the input and output shape of the 'func' argument '' . ”,有人知道为什么吗?我在这里附上我的代码(它可能看起来很复杂,但实际上只是一个非线性... -
1 votesanswersviews
将Hessian提供给optim()
我有一个目标函数(对数似然),我希望使用R作为输入向量(参数)最大化 . 我有函数的渐变(即得分向量),我也碰巧知道函数的Hessian . 在Matlab中,我可以通过使用optimset('GradObj','on')和optimset('Hessian','on')将渐变和Hessian包含在最小化中,轻松地最大化函数并显着提高性能 . 特别是,后者在这种情况下产生了巨大的差异 . 但是,我... -
1 votesanswersviews
如何为nlm提供渐变(或粗麻布)?
在非线性最小化函数 stats::nlm 的帮助文件中,它声明其中一个参数 f 是(强调我的): 要最小化的函数,返回单个数值 . 这应该是一个函数,第一个参数是长度为p的向量,后跟由...参数指定的任何其他参数 . 如果函数值具有称为渐变的属性或梯度和粗麻布属性,则这些属性将用于计算更新的参数值 . 这是否意味着渐变和hessians由以下命令提供: attr(f, 'gradient') ... -
1 votesanswersviews
R优化:每次迭代都将值从函数传递到梯度
我有一个函数,我正在使用R中的 optimx 函数进行优化(我也打开使用 optim ,因为我试图这样做) . 我有一个渐变,我传递给 optimx (希望)与不使用渐变相比更快收敛 . 函数和梯度都使用从每个新参数集计算的许多相同量 . 特别是这些量中的一个在计算上非常昂贵,并且它试图找到一种计算该量的方法,然后将其传递给函数和梯度 . 所以这就是我在做的事情 . 到目前为止,这是有效的,但效率... -
1 votesanswersviews
MIP与CP用于调度问题
众所周知,精确的数学方法对于解决柔性作业车间问题的大型实例并不高效 . 在文献中常见的是关于FJS问题的MILP公式 . 我读到有趣的是将MILP模型用于涉及非精确方法的实验作为元启发式(GA,FA,TS等),因为它提供了下限和上限 . 我还读到,当找到可行的解决方案比最优解决方案更重要时,应该选择CP . 这是真实的陈述吗? -
1 votesanswersviews
隐马尔可夫模型初始猜测
我正在训练大型英文文本(布朗语料库的前50,000个字符,包括字母和空格)的2状态HMM,我的算法遵循Mark Stamp的教程(www.cs.sjsu.edu/~stamp/RUA/) HMM.pdf) . 由于观察仅包括26个字母和空格,所以最初我给每个观察(在一个状态内)的概率为1/27,然后将每个修改为0.0001,同时保持行随机 . 运行训练器50次迭代使我在log [P(O | lam... -
0 votesanswersviews
CPLEX的限制
我们在目标函数和1500条件中涉及100K变量的问题 . CPLEX是否支持这样的数字?如果是这样,CPLEX的限制和性能基准是什么? 问候, 佳日 -
0 votesanswersviews
Cplex Java用于解决小型模型
我试图用额外的约束(权重)来解决类似于TSP的问题,我已经用Java编写了我的代码,但是我在编译时遇到了问题:这是我的代码: package example01; import ilog.concert.*; import ilog.cplex.*; import java.io.*; import java.util.Arrays; public class modele ... -
1 votesanswersviews
凹面船体在边界上取多边形的所有点
在我的工作中,我必须在边界中包含一些随机的点组 . 凸形船体占用了额外的空间而且没有最紧凑的形状,所以我修改它以通过以下方式放松边缘: i)为给定的点数绘制凸包 . ii)现在对于不在凸包边界上的每个点,检查是否可以将其添加到边界(当然这会改变边界整形),同时确保没有给定的点位于新的多边形形状之外 . (点多边形算法) iii)如果所有点都位于多边形内部,则重复步骤2以获得其他点 . iv)如果... -
5 votesanswersviews
对于fn和gr,R optim具有相同的功能
我想使用optim()来优化成本函数(fn参数),我将提供一个渐变(gr参数) . 我可以为fn和gr编写单独的函数 . 但是,它们有很多共同的代码,我不希望优化器浪费时间重复这些计算 . 那么有可能提供一个计算成本和梯度的函数吗?如果是这样,那么optim()的调用语法是什么? 举个例子,假设我想要最小化的函数是 cost <- function(x) { x*exp(x) } ... -
1 votesanswersviews
即使代码没有错误,我也无法模拟CPLEX模型
我对我正在尝试编写的CPLEX代码有一些疑问 . 代码本身(模型)似乎编写得很好,但是当填写数据时我有错误 . NOTE: there are no constraints in the model because i'm trying to make a trial of the model just to see that it works. 这是代码: using CP; // NET... -
3 votesanswersviews
将分布拟合到数据:如何惩罚“坏”参数估计?
我正在使用scipy's least-squares optimization来调整exponentially-modified gaussian distribution到一组反应时间测量值 . 一般情况下,它运行良好,但有时候,优化会偏离轨道并为参数选择一个疯狂的值 - 结果图显然不能很好地拟合数据 . 一般来说,看起来问题是浮点精度误差引起的 - 我们进入0或inf或nan-land . 我... -
0 votesanswersviews
如何为雅可比提供Fsolve?
pow=fsolve(@eqns,pop); 这是我用来解决2x2非线性方程组的代码,在函数eqns.m中定义 . pop是一个非常接近解决方案的2x1初始化向量 . 当我运行它时,输出说 No solution found.fsolve stopped because the relative size of the current step is less than the default ... -
1 votesanswersviews
数据拟合3D空间中的椭圆
论坛 我有一组数据显然在3D空间中形成一个椭圆(不是椭球,而是3D中的曲线) . 受到跟随线程http://au.mathworks.com/matlabcentral/newsreader/view_thread/65773的启发并在某人的帮助下,我设法运行优化代码并输出一组最佳参数x(向量) . 但是,当我尝试使用此x来复制椭圆时,结果是空间中的奇怪直线 . 我已经被困在这里好几天了,仍然不知... -
2 votesanswersviews
MATLAB中使用Opti Toolbox的混合整数二次规划
我希望使用MATLAB中的OPTI工具箱解决具有线性约束的混合整数二次规划 . 我希望我的一些决策变量是连续的,一些决策变量是二进制的 . 我该如何指定? -
0 votesanswersviews
在分支或自定义分支规则之前重置优先级
给定MIP求解器即将选择一个变量进行分支的节点,我想建议一小部分变量可供选择,但与解算器的启发式断开关系 . 我有充分的理由相信这可以显着减少解决整数编程问题所需的时间 . I prefer Gurobi (Python API) 但如果有必要,我愿意切换到另一个求解器(SCIP,CPLEX) . Issues: I failed to figure out which Gurobi c...