37
4
I have searched on stackoverflow for people facing similar issues and this topic Replicating MATLAB's `randperm` in NumPy is the most similar.
我在stackoverflow上搜索过面临类似问题的人，这个主题在NumPy中复制MATLAB的“randperm”是最相似的。
However, although it is possible to recreate the behavior of randperm function from Matlab in Python using numpy random permutation, the numbers generated are not the same, even though I choose the same seed generator for both languages. I am a bit confused since my tests were relevant for other random functions between Matlab and Python.
然而，尽管可以使用numpy随机排列在Python中从Matlab中重新创建randperm函数的行为，但生成的数字并不相同，即使我为两种语言选择了相同的种子生成器。我有点困惑，因为我的测试与Matlab和Python之间的其他随机函数有关。
Here is what I have tried:
以下是我尝试过的：
Matlab
Matlab
个人中心
文章发布
退出rng(42);
randperm(15)
which returns
返回
ans =
11 7 6 5 15 14 1 4 9 10 3 13 8 2 12
Python
蟒蛇
np.random.seed(42)
print(np.random.permutation(range(1,16)))
which returns
返回
[10 12 1 14 6 9 3 2 15 5 8 11 13 4 7]
How can I change my Python code so it can reproduce the same order of random numbers than Matlab ?
我如何更改我的Python代码，使其能够重现与Matlab相同的随机数顺序？
if you want to guarantee the equality of outcomes, call Matlab from Python using its API. Otherwise, you either know that their random generation algorithm is equal or you have no guarntee.
如果您想保证结果的平等性，请使用其API从Python调用Matlab。否则，你要么知道他们的随机生成算法是相等的，要么就没有保证。
It seems that Matlab and Numpy use the same random number generators by default, and the discrepancy is caused by the inner workings of randperm being different in the two languages.
默认情况下，Matlab和Numpy似乎使用相同的随机数生成器，而这种差异是由于两种语言中randperm的内部工作方式不同造成的。
In old Matlab versions, randperm worked by generating a random array and outputting the indices that would make the array sorted (using the second output of sort). In more modern Matlab versions (I'm using R2017b), randperm is a built-in function, so the source code cannot be seen, but it seems to use the same method:
在旧的Matlab版本中，randperm通过生成一个随机数组并输出索引来进行排序（使用排序的第二个输出）。在更现代的Matlab版本中（我使用的是R2017b），randperm是一个内置函数，因此看不到源代码，但它似乎使用了相同的方法：
>> rng('default')
>> rng(42)
>> randperm(15)
ans =
11 7 6 5 15 14 1 4 9 10 3 13 8 2 12
>> rng(42)
>> [~, ind] = sort(rand(1,15))
ind =
11 7 6 5 15 14 1 4 9 10 3 13 8 2 12
So, if the random number generators are actually the same in the two languages, which seems to be the case, you can replicate that behaviour in Numpy by defining your own version of randperm using argsort:
因此，如果这两种语言中的随机数生成器实际上是相同的（似乎是这样），那么您可以通过使用argsort定义自己版本的randperm来在Numpy中复制这种行为：
>>> import numpy as np
np.random.seed(42)
ind = np.argsort(np.random.random((1,16)))+1
print(ind)
[[11 7 6 5 15 16 14 1 4 9 10 3 13 8 2 12]]
Note, however, that relying on the random number generators being the same in the two languages is risky, and probably version-dependent.
然而，请注意，依赖两种语言中相同的随机数生成器是有风险的，并且可能依赖于版本。
I sometimes am deeply surprised learning that MATLAB implements a weird, sub-optimal algorithm for a standard task with a standard, efficient solution.
当我得知MATLAB用标准高效的解决方案为标准任务实现了一种奇怪的次优算法时，我有时会深感惊讶。
@CrisLuengo Do you mean this? That seems to be linear, an sorting is about n*log(n). so not much more expensive. Also, the linear algorithm requires looping, which was an issue in old Matlab versions. Note that randperm was an m-file. And you have to admit that [~, out] = sort(in); looks elegant :-)
@CrisLuengo你是这个意思吗？这似乎是线性的，排序大约是n*log（n）。所以不会贵多少。此外，线性算法需要循环，这在旧的Matlab版本中是一个问题。请注意，randperm是一个m文件。你必须承认[~，out]=排序（in）；看起来很优雅：-）
I guess the sorting method is more efficient if you don’t have a JIT and can rely on fast, vectorized sort and rand. But in a compiled language it’s much less efficient, it uses more memory and sorting is a lot more expensive. O notation talks about how it scales with array size, not about how much it costs. Fisher-Yates doesn’t use any comparisons at all, sorting needs to do O(n log n) comparisons. Sorting also needs to move each array element multiple times, Fisher-Yates does exactly n-1 swaps.
我想，如果你没有JIT，并且可以依靠快速、矢量化的排序和rand，排序方法会更有效。但在编译语言中，它的效率要低得多，占用的内存更多，排序也要贵得多。O表示法谈论的是它如何随数组大小而扩展，而不是它的成本。Fisher Yates根本不使用任何比较，排序需要进行O（n-logn）比较。排序还需要多次移动每个数组元素，Fisher Yates正好进行n-1次交换。
37
4
0
前言： 有时候，一个数据库有多个帐号，包括数据库管理员，开发人员，运维支撑人员等，可能有很多帐号都有比较大的权限，例如DDL操作权限(创建，修改，删除存储过程，创建，修改，删除表等），账户多了，管理
所以我用 Create React App 创建并设置了一个大型 React 应用程序。最近我们开始使用 Storybook 来处理和创建组件。它很棒。但是，当我们尝试运行或构建应用程序时，我们不断遇
遵循我正在创建的控件的代码片段。这个控件用在不同的地方，变量也不同。 我正在尝试编写指令来清理代码，但在 {{}} 附近插入值时出现解析错误。 刚接触 Angular ，无法确定我错过了什么。请帮忙。
我正在尝试创建一个 image/jpeg jax-rs 提供程序类，它为我的基于 post rest 的 Web 服务创建一个图像。我无法制定请求来测试以下内容，最简单的测试方法是什么？ @POST
我一直在 Windows 10 的模拟器中练习 c。后来我改用dev C++ IDE。当我在 C 中使用 FILE 时。创建的文件的名称为 test.txt ，而我给出了其他名称。请帮助解决它。 下面
当我们创建自定义 View 时，我们将 View 文件的所有者设置为自定义类，并使用 initWithFrame 或 initWithCode 对其进行实例化。 当我们创建 customUITable
我正在尝试为函数 * Producer 创建一个线程，但用于创建线程的行显示错误。我为这句话加了星标，但我无法弄清楚它出了什么问题...... #include #include #include
今天在做项目时，遇到了需要创建JavaScript对象的情况。所以Bing了一篇老外写的关于3种创建JavaScript对象的文章，看后跟着打了一遍代码。感觉方法挺好的，在这里与大家分享一下。 &
我正在阅读将查询字符串传递给 Amazon 的 S3 以进行身份验证的文档，但似乎无法理解 StringToSign 的创建和使用方式。我正在寻找一个具体示例来说明 (1) 如何构造 String
前言:我对 C# 中任务的底层实现不太了解，只了解它们的用法。为我在下面屠宰的任何东西道歉: 对于“我怎样才能开始一项任务但不等待它？”这个问题，我找不到一个好的答案。在 C# 中。更具体地说，即使任
我有一个由一些复杂的表达式生成的 ILookup。假设这是按姓氏查找人。 (在我们简单的世界模型中，姓氏在家庭中是唯一的) ILookup families; 现在我有两个对如何构建感兴趣的查询。 首
我试图创建一个 MSI，其中包含 和 exe。在 WIX 中使用了捆绑选项。这样做时出错。有人可以帮我解决这个问题。下面是代码: 错误 error LGH
在 Yii 中，Create 和 Update 通常使用相同的形式。因此，如果我在创建期间有电子邮件、密码、...other_fields...等字段，但我不想在更新期间专门显示电子邮件和密码字段，但
上周我一直在努力创建一个给定一行和一列的 QModelIndex。 或者，我会满足于在已经存在的 QModelIndex 中更改 row() 的值。 任何帮助，将不胜感激。 编辑: QModelInd
出于某种原因，这不起作用: const char * str_reset_command = "\r\nReset"; const char * str_config_command = "\r\nC
现在，我有以下由 original.df %.% group_by(Category) %.% tally() %.% arrange(desc(n)) 创建的 data.frame。 DF 5),
在今天之前，我使用/etc/vim/vimrc来配置我的vim设置。今天，我想到了创建.vimrc文件。所以，我用 touch .vimrc cat /etc/vim/vimrc > .vimrc 所
我可以创建一个 MKAnnotation，还是只读的？我有坐标，但我发现使用 setCooperative 手动创建 MKAnnotation 并不容易。 想法？ 最佳答案 MKAnnotation
在以下代码中，第一个日志语句按预期显示小数，但第二个日志语句记录 NULL。我做错了什么？ NSDictionary *entry = [[NSDictionary alloc] initWithOb
我正在使用与此类似的代码动态添加到数组； $arrayF[$f+1][$y][$x+1] = $value+1; 但是我在错误报告中收到了这个: undefined offset :1 问题:尝试创
我是一名优秀的程序员,十分优秀！