首页 文章
  • 1 votes
     answers
     views

    使用Sympy评估一系列公式时减少计算负担

    首先,我不完全确定这是正确的stackexchange社区 . 我也在数学论坛上发帖,我认为它应该没问题,因为看问题的人应该对溢出与数学有不同的观点,但只是告诉我这是不是很糟糕我会删除其中一个帖子: https://math.stackexchange.com/questions/2330004/reducing-the-computational-burden-in-evaluating-a-s...
  • 2 votes
     answers
     views

    使用lambdify和函数评估优化运行时

    我目前正在优化我的代码的运行时间,它仍然不在我想要的消费范围内 . 在执行高斯求积时,我已经达到了 80% of the time is spent on running lambdify() on my sympy Matrix expressions and evaluating the resulting lambda functions 的程度 . 代码的所有其他方面都得到了充分的优化,因...
  • 0 votes
     answers
     views

    Scipy的优化与Multiprocessing不兼容?

    在尝试使用Scipy的优化算法来最小化在子过程中计算其值的函数时,我发现基于梯度的算法(到目前为止流域购物和L-BFGS-B)在optimize.py的第562行遇到以下错误: grad [k] =(f(*((xk d,)args)) - f0)/ d [k]TypeError:不支持的操作数类型 - :'NoneType'和'NoneType' 以下是生成此错误的代码的简单示例: impo...
  • 0 votes
     answers
     views

    scipy.optimize.minimize调用scipy.sparse.csr_matrix时出错

    我试图在filowing方式中使用scipy.optimize.minimize: fi =最小化(fun = quadratic_function,x0 = fi,args =(alfa,Ks,Kf),method ='Nelder-Mead') . x 其中fi是要优化的参数 . 并且quadratic_function具有这样的性质,即在每次调用时,都应该构造一个scipy.sparse...
  • 0 votes
     answers
     views

    Scipy.optimize.minimize目标函数ValueError

    我使用scipy.optimize.minimize来解决9个自由变量的小优化问题 . 我的目标函数基本上是另一个函数的包装器,如果我评估我的目标函数,返回类型是'numpy.float32'...这是一个标量?但是,我在尝试使用最小化函数时收到以下错误: raise ValueError("Objective function must return a scalar") V...
  • 1 votes
     answers
     views

    MySQL优化 - 使用索引

    我在表a上有以下INDEX: SHOW INDEX FROM a Table,Non_unique,Key_name,Seq_in_index,Column_name,Collation,Cardinality,Sub_part,Packed,Null,Index_type,Comment,Index_comment a,'1','oc_id','1','OC_ID','A','19 ',NU...
  • 0 votes
     answers
     views

    宏用于比较和复制数据从一张纸到另一张需要很长时间

    我使用此宏将内容从一个Excel工作表复制到另一个工作表,方法是比较两列并找到匹配的单元格 . 问题是这个宏需要很长时间(接近三天)才能完成 . 两张纸都有近4,00,000条记录可供比较 . 有人可以帮助我加快速度吗? Option Explicit Sub MatchAndCopy() Dim sheet01 As Worksheet, sheet02 As Works...
  • 8 votes
     answers
     views

    关于开销的HTTP与Websockets

    我正在使用两端的node.js构建一个文件同步程序(与Dropbox不同) . 我需要有数千个客户端同时请求数据 . 这是我目前的系统: Server通过websocket将通知推送到客户端(文件已更新) 客户端队列在空闲时下载并发出HTTP请求 我将以压缩块的形式提供数据,例如每个50 MB,因此HTTP请求开销(标头)可以忽略不计 . 如果我使用websockets进行请求 and...
  • 0 votes
     answers
     views

    优化Innodb表以释放可用空间

    我无法在innodb表中释放空间 . 我还设置了 innodb_files_per_table=1 所以每个表都有单独的ibd文件 . 当我查询information_Schema以检查数据大小,索引大小和data_free时,它总是在某些表的data_free列中显示一些值,如7Mb . 查询: SELECT table_schema“数据库名称”,table_name,(data_lengt...
  • 1 votes
     answers
     views

    计算Excel中一堆数字的平均百分比?

    基本上,我正在尝试计算我正在使用的工作表的总平均停机时间,但是我试图让它进一步细分 . 这是Excel表格的图片: Excel Calculations (不允许添加嵌入式图片:( - 需要10个声望,抱歉!) 差不多,我试图通过计算平均停机时间然后从100.00%减去它来计算当月的平均正常运行时间 . 我有什么工作,但我正在努力研究停电时间列是否可以报废,总计可以用一个更大的公式来计算 . ...
  • 10 votes
     answers
     views

    如何在Linux核心i3 / i7中击败硬件预取器

    我试图找到一种方法来打败H / W预取器来检测流模式并以随机顺序访问4KB数据,这样就不会被H / W预取器检测和预取 . 最初我想以随机模式访问所有偶数索引数据,因为H / W预取器总是预取下一个缓存行(所以当我访问偶数索引时,下一个奇数索引数据已被预取) . 我编写代码以随机模式访问所有偶数索引数据,但结果表明预取器检测到模式(不知道如何?没有固定步幅,都是随机步幅) 我正在调查原因 - 为什...
  • 4 votes
     answers
     views

    MySQL是否自动最小化重复的VARCHAR存储?

    问题是针对运行在Ubuntu 10.04 LTS服务器上的MySQL 5.5,默认的InnoDB表类型... 假设我有一个表“地址”的房屋地址,其中包含“数字”,“街道”,“地区”,“城镇”,“县”和“邮政编码”栏 . 我将在这些列中有许多具有相同值的行,并且我将单独索引它们以进行搜索 . 假设我将每列实现为VARCHAR(127)并使用town ='London'创建1000行 . 这是否意味着...
  • 0 votes
     answers
     views

    OpenCV - 如何优化色彩跟踪程序?

    我想优化我的程序,其中我使用的颜色对象跟踪算法描述为here . 唯一的区别是我使用的是cvBlob库,而不是cv :: moments(cvBlob更快更准确) . 使用profiler(valgrind kcachegrind)我发现~29%的时间是由颜色空间转换方法(cv :: cvtColor;我跟踪三种颜色的对象) . 我正在从BGR转换为HSV . 我在一些论文中读过,使用YCbCr色...
  • 0 votes
     answers
     views

    用PYOMO / PYSP解决MINLP问题

    球队, 目前我正在研究非线性随机优化问题 . 到目前为止,工具箱真的很有帮助,谢谢!但是,添加非线性约束会导致错误 . 我使用gurobi求解器 . 该问题源于以下约束 . def max_pcr_power_rule(model, t): if t == 0: return 0 <= battery.P_bat_max-model.P_sc_max[t+1]-mod...
  • 0 votes
     answers
     views

    二次程序公式:Matlab

    鉴于此目标函数: 最小化: f = (Ax + By)' * G * (Ax + By) 受到一些平等和不平等的影响 . 其中 x 和 y 分别是具有 p 和 q 元素的实值向量(决策变量) . A 大小 m * p , B 大小 m * q ,G是大小为 m * m 的对称矩阵 . 我的问题是如何以 v' * G * v 的形式编写 f ,以便它可以很容易地在quadprog中使用 . 换...
  • 0 votes
     answers
     views

    Matlab中具有线性等式约束的二次规划

    我必须在一些线性约束下识别ARX,这意味着我有一个带线性等式约束问题的二次规划 . 一种方法是在红色框中使用以下等式 . 在这种情况下可能的缺点是矩阵求逆的计算(有时Matlab给我警告:矩阵接近单数或严重缩放) 另一种方法是在Matlab中使用命令:quadprog() 另一种方法是在Matlab中使用命令:lsqlin() 这三种方法中哪一种最好? 哪个数字最强大?
  • 0 votes
     answers
     views

    为什么SCIP服用这么长时间并且记忆如此之多?

    我在matlab的OPTI工具箱中使用SCIP求解器来解决具有整数约束的二次优化问题 . 我使用以下规格运行它并且它已经运行了一天并且已经在我的系统中占用了55GB的ram并且仍在计数中 . 我是matlab中的优化新手,我做错了什么或者这是常见的吗?我尝试使用较少的maxnodes和maxtime,但程序在这些情况下停止并显示“已达到节点限制”错误 . 这是代码(H,Aeq等已经在代码中定义) ...
  • 1 votes
     answers
     views

    当原始模型可行时,为什么SCIPcopy模型不可行?

    我是SCIP的新手,所以我不确定这是一个错误还是我做错了什么 . 我有一个完全使用SCIP解决的MIP实例,但是当我试图解决模型的副本时,SCIP说这是不可行的 . 当关闭presolve时,它似乎更明显 . 我正在使用预先构建的SCIP v3.2.0的Windows . 该模型只有二进制和整数变量 . 以下代码概述了我的尝试: SCIP* _scip, subscip; SCIPcreate(&...
  • 0 votes
     answers
     views

    对于非线性凸模型,SCIP求解器中“第一LP值”的确切含义是什么

    我用二元变量,线性目标函数和一组用单个 non-linear 约束修正的线性约束解决了数学模型,使得模型成为非线性二元问题 . 在SCIP提供的输出文件中,有一个名为: First LP value 的术语,并且已赋值 . 对于我的特定非线性问题,我无法弄清 First LP value 究竟是什么意思?我非常感谢任何解释 .
  • 0 votes
     answers
     views

    solve.QP不能最大化正定二次函数

    我有一个非常简单的二次问题,原则上应该可以通过solve.QP(quadprog)来解决: max x^2+3*y^2 subject to x>=0, y>=0, x+y <=1 问题是将最大化转化为所需的最小化,矩阵(显然)变为负定 . 这不是与小特征值或舍入有关的问题 . 我认为解决问题.QP是一个普通求解器,但是尽管我努力阅读在线材料,看起来好像你不能使用solve.QP...
  • 0 votes
     answers
     views

    希望快速并行优化Lapack常规(gfortran)问题

    直到两天前我才认为我对此有一个非常清楚的认识,但现在我可能会过度思考它并使自己感到困惑 . 我将解释我正在做什么,然后提出一些可能过于简单的问题,但到目前为止我已经搜索并找到了相互矛盾的答案 . 当然有人可以帮我 . 我编写了一个使用LAPACK例程来解决特征值问题的fortran代码 . 我的问题设置是(A-LB)x = 0,其中L是我的特征值,x是我的特征向量,A和B是方形,复数,非对称,非埃...
  • 1 votes
     answers
     views

    为了对数似然估计,保持矩阵正定的约束

    我试图通过R中的GARCH模型 (VECH(1,1)) 来优化对数似然函数并计算我的协方差矩阵 . 我的问题是,如何编制使我的协方差矩阵保持正确的约束? 我知道我的矩阵alpha beta = Q其中Q的特征值必须在单位圆中 . 但是我怎么编程呢?输入是 getSymbols('^GDAXI', src='yahoo', return.class='ts',from="2005-01-0...
  • 0 votes
     answers
     views

    Cplex可以将两个稀疏矩阵作为Q的输入

    我试图使用CPLEX MATLAB API最小化具有线性约束的非常大的二进制二次规划问题 . 然而,二次函数f = x'Qx具有非常密集的矩阵Q.我可以用非常稀疏的矩阵U和A重写Q:Q = U A. 因此,有没有办法将这些矩阵传递给CPLEX而不是完整的矩阵Q(它太大而无法存储)?
  • 0 votes
     answers
     views

    在Gurobi Python Shell中访问Objective Q矩阵

    我正在尝试在Gurobi和Python解决的优化问题(具有二次约束和二次目标函数)中重新构造一个Objective Q矩阵 . Gurobi可以选择将约束和目标函数作为线性表达式添加而不是摆弄矩阵,因此我没有原始矩阵,Gurobi通过我的客观方程和系数为我创建它 . 为了对Objective Q矩阵的凸性psd性质进行分析,我需要有Q(目标矩阵)和A(约束矩阵) . 有谁知道gurobi.py ...
  • 1 votes
     answers
     views

    编译SCIP:找不到'spxsolver.h'文件

    我正在尝试使用以下命令编译SCIP: make ZLIB=false READLINE=false GMP=false 任何人都知道如何处理这个问题? -> compiling obj/O.darwin.x86_64.gnu.opt/lib/objscip/objsepa.o -> compiling obj/O.darwin.x86_64.gnu.opt/lib/objscip/o...
  • 0 votes
     answers
     views

    当达到迭代限制时,双向单形中的目标约束

    我完成了,显然既没有设置模型的 ObjBound 属性,也没有设置约束的 pi 属性 . 有没有办法检索这些属性,即使双单纯形方法没有运行完成? 谢谢!
  • 0 votes
     answers
     views

    并行优化pyomo-gurobi

    我对Pyomo和优化都很新 . 我试图解决一个迭代的局部搜索算法,我必须多次解决“相对较小”的优化问题 . 目前,我已经实现了一个循环,每次,我制定优化问题并用Gurobi解决它 . 但是,我想知道是否有一种“并行化”的方式来解决这个问题,因为本地搜索的结果并没有相互影响 . 我查看了ThreadPoolExecutor和ProcessPoolExecutor . 第一个没有提高速度,而第二个给我...
  • 1 votes
     answers
     views

    在AMPL中实现CDF(或CDF的反转)

    我的问题涉及一些数学,所以我输入了Latex格式并拍了一下屏幕截图 . 我想知道如何在AMPL中实现这一点 . 在Matlab fmincon中,您可以轻松传递CDF(或CDF的反转),但我不知道在AMPL中我们如何做到这一点 . 有没有什么方法可以在AMPL中实现CDF(或CDF的反转)的功能 .
  • 2 votes
     answers
     views

    用Matlab优化

    在matlab中编写目标函数(包括梯度向量)的常用方法如下: [L,G] = objfun(x) 其中 L 是目标函数的值, G 是梯度向量, x 是我想要优化的系数向量 . 但是,当我包含另一个输入(即 [L,G]=objfun(x,M) ,其中 M 是一个矩阵)或当我在函数 objfun 中调用另一个函数时,代码没有运行 . 如何通过保持这种优化格式来包含任何输入并调用 objfun 中的任...
  • 0 votes
     answers
     views

    在Matlab中使用Opti-toolbox优化运动规律

    我有一个关于使用Matlab中的Opti-toolbox进行理想运动定律的数值推导的问题 . 要遵循许多约束,但目标函数的主要原理如下所示: 此函数是周期性的,但仅使用1个周期 . 所有变量都是无量纲的,点数为n(奇数) . 在0和1处,函数产生0 . 在中间,在(n 1)/ 2处,该函数产生例如0.3 . 该功能需要在前半部分单调增加,在下半部分减少 . 速度需要介于-...

热门问题