PA(find(PS==i))=S2(i+1); find(PS==i)返回的是PS矩阵中数值为i的坐标。

设为m 所以就是把S2(i+1)赋给PA中第m个数。