matlab dsort,matlab中sort函数的用法

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

help文档中第一句就说明，可以用sort将参量中的元素按升序或降序排列~

格式：

B = sort(A) %沿着输入参量

A的不同维的方向、从小到大重新排列 A中的元素。A 可以是字符串的、实数的、复数的单元数组。对于 A 中完全相同的元素，则按它们在 A

中的先后位置排列在一块；若 A 为复数的，则按元素幅值的从小到大排列，若有幅值相同的复数元素，则再按它们在区间[-π ,π

]的幅角从小到大排列；若 A 中有元素为NaN，则将它们排到最后。若

A为向量，则返回从小到大的向量，若A为二维矩阵，则按列的方向进行排列；若A为多维数组，sort(A)把沿着第一非单元集的元素像向量一样进行处理。

B=sort(A,dim) %沿着矩阵A中指定维数dim方向重新排列A中的元素。

[B,IX]=sort(A);

%其中IX是一个大小等于size(A)的数组，其每一列是A中列向量的元素相对应的置换位置记号。

举例：A = [ 3 3 5

0 4 2 ];

sort(A,1) %纵向排列

ans =

0 3 2

3 4 5

sort(A,2) %横向排列

ans =

3 3 5

0 2 4

>>A = [-1.9, -0.2, 3.1415926, 5.6, 7.0, 2.4+3.6i];

>>[B1,INDEX] = sort(A)

计算结果为：

B1 =

Columns 1 through 4

-0.2000 -1.9000 3.1416 2.4000 + 3.6000i

Columns 5 through 6

5.6000 7.0000

INDEX

=

2 1 3 6 4 5

由此可知，如果A与另一个矩阵C(char或cell)一一对应，如果用[B1,INDEX] =

sort(A)对A进行排序，则可以取得索引INDEX，进而可以查询B1中每个元素与C中哪一个对应。

-----------------------------------------------------------------------------

Matlab中给一维向量排序是使用sort函数：sort(A)，排序是按升序进行的，其中A为待排序的向量；若欲保留排列前的索引，则可用

[sA,index] = sort(A) ，排序后，sA是排序好的向量，index 是 向量sA 中对 A 的索引。

索引使排列逆运算成为可能。

有久没有用Matlab都忘了。。。。。

事实上，这里A ≡sA(index) ,

[A恒等于sA(index)]，这个结论确实很奇妙，而且很有用。不信你排序之后试下键入命令sA(index)

，看看得到的是不是就是排列前的A呢。[因为鸿1018问到这个问题，特此补写这一段] 2008-10-20 13:18 +

对于一般矩阵An*m(n，m均大于等于2)，只要修正一下index的值，结论sA(index) ≡

A仍然成立，修正代码如下：

s=size(A);

for k=1:s(2)

index(:,k)=index(:,k)+(k-1)*s(1);

end

-----------------------------------------------------------------------------------

在Matlab中排序某个向量(一维)时，可以使用sort(A)，其中A为待排序的向量，如果仅是用来排序A，那么直接使用sort(A)即可，

如果排序后还需要保留原来的索引可以用返回值，即[B,ind]=sort(A)，计算后，B是A排序后的向量，A保持不变，ind是B中每一项对应于A

中项的索引。排序是安升序进行的。

在Matlab中，访问矩阵中的元素，一维用A(1)访问向量A的第一个元素；(下标从1开始)；二维用A(1，2)访问A中第一行，第二列的元素。

由于在sort函数的结果中，是安升序排序的，要转换成降序，先用X＝eye(n)生成一个n维的单位阵，然后用X＝rot90(X)将其旋转为次对角线的单位阵，再用原来矩阵乘以X即可，如要讲A逆序排列采用如下步骤：

X=eye(size(A));

X=rot90(X);

A=A*X;

假如a是一个2*n的矩阵,即两行.

b=a(1,:);

[c,pos]=sort(b);%pos为排序后的下标,c为第一行的排序结果;

a(2,:)=a(2,pos);%第二行按照第一行排序的下标对应

a(1,:)=c;%第一行结果重新赋给a 的第一行

以下适用于m*n的矩阵按第一行排序

[ b, pos ] = sort( a( 1, : ) );

a = a( :, pos );

------------------------------------------------

主要是前一篇在实现几种编码的时候遇到了好些问题。

1.元胞数组

因为要生成各个符号的码字，他们的长度都不一样，

没法用一般的矩阵来存。

c = cell(5,2)

c{1,1} = 1

反正用cell创建，用{}引用就是了。

另外，把若干字符串连接成一个元胞数组：

strvcat('Hello','Yes')

(特别鸣谢yhsa678同学的指点)

2.for循环

matlab的for让人觉得很不爽。

c:

for ( i = 1 ; i <= n ; i ++ ){...}

vb:

for i = 1 to n

...

next i

pascal:

for i := 1 to n begin

...

end

以上的东西执行完后，i的值都是n+1

但是，matlab的for i = 1:n完了之后！

i == n

。。。这在平时真的没啥。。但在有的时候，就要增加点编码了，很烦。

比如这个例子：

我们要寻找数组中有没有0。有0就干什么事，没有就干什么事。

x = [1 3 2] ;

for i = 1:3

if x(i) ==

0

break ;

end

end

这个完了之后，如果i==3，

那么你就无法确定是因为3好元素是0，break掉的。

还是循环执行完了，i停在3的。

所以相比前几种语言来说，这个设计实在是不爽。。

只有多加点东西：

x = [1 3 2] ;

has0 = 0 ;

for i = 1:3

if x(i) ==

0

has0 = 1

break ;

end

end

if has0

..

else

..

end

3.sort

第一个就是，sort是从小到大，那么从大到小就有点囧了。

可以用逆置数组的办法来完成：

p = sort(p)

p = fliplr(p)

这下就是从大到小了。

fliplr是左右逆置。

还有个flipud是上下逆置。

第二个就是，sort的返回值确实很有用处，当年还不会matlab的时候在这里看到过。

今天做编码的时候，有两个要先从大到小排序再处理。

但是处理完之后交回调用者的码字应该是按照原顺序排的。

这个时候，就好利用sort的返回值了。

>> x = [ 3 1 2]

x =

3 1 2

>> [y i] = sort(x)

y =

1 2 3

i =

2 3 1

执行完后，有这个关系：y = x(i)

比如，调用这要求得到的c是对应于以前的x的。

但我们处理的时候，是基于排序后的y的。那么这时得到的c。应该做个处理：

c(i) = c ;

这下就OK了

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

给你举了例子：

X=magic(5)

X =

17 24 1 8 15

23 5 7 14 16

4 6 13 20 22

10 12 19 21 3

11 18 25 2 9

>> [a,b]=sort(X,2)

a =

1 8 15 17 24

5 7 14 16 23

4 6 13 20 22

3 10 12 19 21

2 9 11 18 25

b =

3 4 5 1 2

2 3 4 5 1

1 2 3 4 5

5 1 2 3 4

4 5 1 2 3

结果解释：

a是原来的矩阵x按照行，每行从小到大重新排列得到的新矩阵。

b告诉你重排的详细信息，也就是做了什么样的变动。

例如b的第一行显示3 4 5 1 2，那么将原矩阵X的第一行的第3 4 5 1

2个元素取出来，顺次排列，就变成a矩阵的第一行。

@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@