Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 其它 > 资料存档
资料存档 资料存档
回复
 
主题工具 显示模式
旧 2019-12-10, 20:48   #1
poster
高级会员
 
注册日期: 2019-11-21
帖子: 3,006
声望力: 66
poster 正向着好的方向发展
帖子 Matlab如何进行排序?

在Matlab中sort()如何工作?
纯Matlab中的代码:
q是一个数组:

q = -0.2461 2.9531 -15.8867 49.8750 -99.1172 125.8438 -99.1172 49.8750 -15.8867 2.9531 -0.2461 在q = sort(roots(q)) ,我得到了:
q = 0.3525 0.3371 - 0.1564i 0.3371 + 0.1564i 0.2694 - 0.3547i 0.2694 + 0.3547i 1.3579 - 1.7880i 1.3579 + 1.7880i 2.4410 - 1.1324i 2.4410 + 1.1324i 2.8365

好吧,似乎工作正常!然后在python中,我使用(q与上面相同,它是一个np.array ):

import numpy as np q = np.sort(np.roots(q))

我得到:

[ 0.26937874-0.35469815j 0.26937874+0.35469815j 0.33711562-0.15638427j 0.33711562+0.15638427j 0.35254298+0.j 1.35792218-1.78801226j 1.35792218+1.78801226j 2.44104520-1.13237431j 2.44104520+1.13237431j 2.83653354+0.j ]

好吧...这两个结果似乎有所不同,因为它们的排序方式不同,那么原因是什么呢?我做错什么了吗?先感谢您!

我的答案:



def sortComplex(complexList): complexList.sort(key=abs) # then sort by the angles, swap those in descending orders return complexList 然后在python代码中调用它,效果很好:p



回答:

SORT的MATLAB文档中:
如果A具有复杂的条目r和s ,则按以下规则sort它们进行sort(A)如果满足以下任一条件,则r在sort(A) s之前出现:
  • abs(r) < abs(s)
  • abs(r) = abs(s)和angle(r) < angle(s)
换句话说,具有复杂条目的阵列首先基于那些条目的绝对值 (即,复数幅值)进行排序,并且具有相同绝对值的任何条目都基于其相角进行排序。

Python(即numpy)对事物的排序不同。 在他的评论中链接到的Amro文档
复数的排序顺序为字典顺序。如果实部和虚部都不都是南,则顺序由实部确定,除非它们相等,在这种情况下,顺序由虚部确定。

换句话说,首先根据条目的实分量对具有复杂条目的数组进行排序,然后根据虚分量对具有相等实数的任何条目进行排序。

编辑:

如果要在MATLAB中重现numpy行为,一种方法是使用函数SORTROWS根据数组条目的实部虚部创建排序索引,然后将该排序索引应用于复杂数组值:

>> r = roots(q); %# Compute your roots >> [junk,index] = sortrows([real(r) imag(r)],[1 2]); %# Sort based on real, %# then imaginary parts >> r = r(index) %# Apply the sort index to r r = 0.2694 - 0.3547i 0.2694 + 0.3547i 0.3369 - 0.1564i 0.3369 + 0.1564i 0.3528 1.3579 - 1.7879i 1.3579 + 1.7879i 2.4419 - 1.1332i 2.4419 + 1.1332i 2.8344

更多&回答...
poster 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛禁用 表情符号
论坛启用 [IMG] 代码
论坛启用 HTML 代码



所有时间均为北京时间。现在的时间是 19:43


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.