- iOS/Objective-C 元类和类别
- objective-c - -1001 错误,当 NSURLSession 通过 httpproxy 和/etc/hosts
- java - 使用网络类获取 url 地址
- ios - 推送通知中不播放声音
我正在使用具有 4 个适应度函数的遗传算法。我有一个文件,其中每一行都列出了单个解决方案的适应度值。遗传算法的目的是最大化这些值。
我需要从文件中列出的解决方案中找出非支配解决方案。我有一个名为“identify_pareto_observations”的函数,用于识别遗传算法输出的非支配解。但是,我所得到的结果不正确。我在下面发布了一些代码,代表我的程序中的内容。我希望有人可以帮助我确定我要去哪里错了。如果有人能简单地解释在 4 维解空间中识别非支配解的过程,那就太好了!
请注意,我很难理解我找到的解释如何识别非支配解决方案的伪代码,所以我只是使用了我找到的代码 here .
感谢您的帮助/时间!
import os
import sys
import numpy as np
fitness_observations_matrix = []
fitness_observations_matrix.append([-447.117,928641,-0,4131.33])
fitness_observations_matrix.append([-838.066,977440,-0,3859.25])
fitness_observations_matrix.append([-803.34,385226,-0,4799.27])
fitness_observations_matrix.append([-790.052,4919.54,-0,6672])
fitness_observations_matrix.append([-403.629,1.9091e+06,-0.775081,3011.89])
fitness_observations_matrix.append([-99.4555,413201,-0,7615.26])
fitness_observations_matrix.append([-78.8706,9472.03,-0,7071.58])
fitness_observations_matrix.append([-1096.82,127109,-0,7227.28])
fitness_observations_matrix.append([-1058.26,1.97533e+06,-0,4959.18])
fitness_observations_matrix.append([-1324.52,2.56381e+06,-2.66017,2669.81])
fitness_observations_matrix.append([-815.861,247594,-0,4674.39])
fitness_observations_matrix.append([-370.815,1.41723e+06,-0,6036.94])
fitness_observations_matrix.append([-202.34,1.21506e+06,-0,3587.2])
fitness_observations_matrix.append([-718.299,3.0879e+06,-0,3903.51])
fitness_observations_matrix.append([-1514.39,2.62657e+06,-0,6960.03])
fitness_observations_matrix.append([-815.018,414420,-0,-0])
fitness_observations_matrix.append([-747.015,2.5201e+06,-0.913635,7153])
fitness_observations_matrix.append([-246.017,683330,-0,10110.5])
fitness_observations_matrix.append([-5884.85,1.40051e+07,-8.74525,3037.08])
fitness_observations_matrix.append([-4045.09,5.83051e+06,-1.99262,4397.97])
def identify_pareto_observations(fitness_observations):
#Convert to np array
np_array_fitness_observations = np.array(fitness_observations)
#Sort on first dimension
np_array_fitness_observations = np_array_fitness_observations[np_array_fitness_observations[:,0].argsort()]
# Add first row to pareto_frontier
np_pareto_frontier = np_array_fitness_observations[0:1,:]
# Test next row against the last row in pareto_frontier
for row in np_array_fitness_observations[1:,:]:
if sum([row[x] >= np_pareto_frontier[-1][x] for x in range(len(row))]) == len(row):
# If it is better on all features add the row to pareto_frontier
np_pareto_frontier = np.concatenate((np_pareto_frontier, [row]))
#Return
print(np_pareto_frontier)
identify_pareto_observations(fitness_observations_matrix)
最佳答案
如果存在一个解决方案 y
,其中 y
不比 x
差,则解决方案 x
是帕累托支配的在任何方面,并且在某些方面严格来说更好。我假设您正在尝试在所有方面最小化您的问题。然后你可以采取以下方法(你发布的代码有类似的风格,但有很多问题):
numpy.lexsort
函数高效地执行。请注意,在此排序列表中,任何元素都不能被列表中后面的元素支配。np_pareto_frontier
初始化为第一个排序元素。由于步骤 1 中的排序,我们知道它是帕累托非支配的。np_pareto_frontier
中没有任何解决方案支配,则将一行添加到 np_pareto_frontier
-- 因为这些元素不受 np_pareto_frontier
中的任何支配np_pareto_frontier
并且不被排序列表中后面的任何东西支配,它们必须是 pareto 非支配的。下面是这种方法的实现:
import numpy as np
def identify_pareto_observations(fitness_observations):
np_array_fitness_observations = np.array(fitness_observations)
np_array_fitness_observations = np_array_fitness_observations[np.lexnp_array_fitness_observations[:,0].argsort()]
np_pareto_frontier = np_array_fitness_observations[[0],:]
for row in np_array_fitness_observations[1:,:]:
dominated = False
for nondom in np_pareto_frontier:
if all([row[x] > nondom[x] for x in range(len(row))]):
dominated = True
break
if not dominated:
np_pareto_frontier = np.concatenate((np_pareto_frontier, [row]))
return np_pareto_frontier
我们可以在您的数据集上对其进行测试:
print identify_pareto_observations(fitness_observations_matrix)
# [[ -8.15018000e+02 4.14420000e+05 0.00000000e+00 0.00000000e+00]
# [ -1.32452000e+03 2.56381000e+06 -2.66017000e+00 2.66981000e+03]
# [ -4.03629000e+02 1.90910000e+06 -7.75081000e-01 3.01189000e+03]
# [ -5.88485000e+03 1.40051000e+07 -8.74525000e+00 3.03708000e+03]
# [ -2.02340000e+02 1.21506000e+06 0.00000000e+00 3.58720000e+03]
# [ -8.38066000e+02 9.77440000e+05 0.00000000e+00 3.85925000e+03]
# [ -4.47117000e+02 9.28641000e+05 0.00000000e+00 4.13133000e+03]
# [ -4.04509000e+03 5.83051000e+06 -1.99262000e+00 4.39797000e+03]
# [ -8.15861000e+02 2.47594000e+05 0.00000000e+00 4.67439000e+03]
# [ -8.03340000e+02 3.85226000e+05 0.00000000e+00 4.79927000e+03]
# [ -1.05826000e+03 1.97533000e+06 0.00000000e+00 4.95918000e+03]
# [ -3.70815000e+02 1.41723000e+06 0.00000000e+00 6.03694000e+03]
# [ -7.90052000e+02 4.91954000e+03 0.00000000e+00 6.67200000e+03]
# [ -1.51439000e+03 2.62657000e+06 0.00000000e+00 6.96003000e+03]
# [ -7.88706000e+01 9.47203000e+03 0.00000000e+00 7.07158000e+03]
# [ -7.47015000e+02 2.52010000e+06 -9.13635000e-01 7.15300000e+03]
# [ -1.09682000e+03 1.27109000e+05 0.00000000e+00 7.22728000e+03]
# [ -9.94555000e+01 4.13201000e+05 0.00000000e+00 7.61526000e+03]
# [ -2.46017000e+02 6.83330000e+05 0.00000000e+00 1.01105000e+04]]
正如我们所见,您只有一个帕累托支配的解决方案,行 [ -7.18299000e+02 3.08790000e+06 0.00000000e+00 3.90351000e+03]
。它由 [ -1.32452000e+03 2.56381000e+06 -2.66017000e+00 2.66981000e+03] 行控制
。
关于python - 无法正确识别 4 维解空间中的 Pareto 观测值,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/23624343/
运行 Tomcat 失败并出现 java.lang.OutOfMemoryError - 与缺少 PermGen 空间相关的错误。 我最近将 Tomcat 更改为以自己的用户(而非 root)运行。
我们有一个表,其中包含数百万行,其中包含 PostGIS 几何图形。我们要执行的查询是:落在边界几何内的最新条目是什么?这个查询的问题是我们经常会有大量的项目匹配边界框(半径大约为 5 公里),然后
我有一个Elasticsearch设置,它将允许用户搜索通配符作为索引。 array:3 [ "index" => "users" "type" => "user" "body" => arra
我创建了一个表,其中每行包含两个按钮,并且两个按钮连接在一起,我想将两个按钮分开。我用过 不起作用,css 也是,这是他们的另一种方式。 我有另一个问题,因为我不想在表格边框内显示操作按钮,而是在靠近
我试图在 jQuery Mobile 中的两个按钮之间留出空白。现实中的布局是这样的: Button 1 Button 2 (Hidden w/ display: none)
按照目前的情况,这个问题不适合我们的问答形式。我们希望答案得到事实、引用或专业知识的支持,但这个问题可能会引发辩论、争论、投票或扩展讨论。如果您觉得这个问题可以改进并可能重新打开,visit the
您好,我对图表应用程序还很陌生。现在我为我的应用程序创建了条形图。当我运行 create bar chart as separate project 时,输出如下所示。 然后当我将条形图与我的应用程序
我在使用 H2 和 GeoDB(内存中,junit)时遇到问题。 另外,使用 Hibernate 5(每个包的最新版本,包括 hibernate-spatial)和 Spring 4。 通过 id 实
我想画一张澳大利亚的 map ,并将每个城市表示为一个点。 然后突出显示人口众多(> 1M)的城市 library(sp) library(maps) data(canada.cities) head
关闭。这个问题是opinion-based .它目前不接受答案。 想改进这个问题?更新问题,以便 editing this post 提供事实和引用来回答它. 6年前关闭。 Improve this
如何保持.txt文件中存在的空格?在.txt文件中,它表示: text :text text1 :text1 text23 :text2 text345 :text3 如果我写这段
以下哪个键最大? 选项 1:16 个数字 [0,9] 选项 2:30 个元音 选项 3:字母表中的 16 个字母 选项 4:32 位 有人可以帮助我,告诉我哪一个是正确的答案以及我们如何计算它吗?我知
在 Unity 3d 中使用 Azure 空间 anchor 来实现在 iOS 和 Android 上部署的室内和室外增强现实体验是否有益? 最佳答案 是的,对于 Azure Spatial Anch
我有一个绝对定位的圆形图像。图像只需占据屏幕宽度的 17%,并且距离顶部 5 个像素。 问题是,当我调整图像大小以占据屏幕宽度的 17% 时,它会这样做,但同时容器会变长。图像本身不会拉伸(stret
我在 Ubuntu 14.04 上使用 Cassandra。从文档中,我可以看到运行命令: nodetool snapshot 创建我的 key 空间的快照。 命令的输出是: nodetool sn
Heroku引入了“私有(private)空间”,是否可以将现有应用迁移到私有(private)空间? https://blog.heroku.com/archives/2015/9/10/herok
是否允许在语义记录中使用非绑定(bind)空格 或其他 HTML 编码字符?我遇到的问题是 ; 字符被软件视为记录的结尾。 例如:假设我有一份婚姻记录,其中包含 2 个结婚者的姓氏、结婚年份以及结
我正在研究“智能 parking ”项目,偶然发现了包含我们真正需要的YouTube视频。我们已经实现了第一部分,即从视频源进行实时透视变换,下一步是将其定义为一组矩形 我基本上需要知道他是如何做到的
我有两个类:Engine 和 Trainset(多个单元),这两个类共享其 ID 空间,其中包含名称和系列 id=- . 这是我的Engine类(它是抽象的,因为有引擎的子类型(DieselEngin
如果有人能帮助我,那就太好了。 我正在尝试使用Java的Split命令,使用空格分割字符串,但问题是,字符串可能没有空格,这意味着它将只是一个简单的顺序(而不是“输入2”将是“退出”) Scanner
我是一名优秀的程序员,十分优秀!