level 4
任道正
楼主
我自己写了一个函数，用于交换矩阵的行或者列，该函数首先将矩阵分为上下两部分并交换位置(交换行)；随后又将矩阵分为左右两部分并交换位置(交换列)；然后返回交换后的矩阵即可！
代码如下：
FFTShift[mtx_] := Block[{dim = Dimensions@mtx},
midx = Round[dim[[1]]/2];
midy = Round[dim[[2]]/2];
mtxtmp = ConstantArray[0, dim];
mtxtmp[[All, 1 ;; midy]] = mtx[[All, -midy ;; -1]];
mtxtmp[[All, midy + 1 ;;]] =
mtx[[All, 1 ;; dim[[2]] - midy]];(*交换行*)
mtxout = ConstantArray[0, dim];
mtxout[[1 ;; midx, All]] = mtxtmp[[-midx ;; -1, All]];
mtxout[[midx + 1 ;;, All]] =
mtxtmp[[1 ;; dim[[1]] - midx, All]];(*交换列*)
mtxout];
然而可以看到以上函数其实就得定义mtxtmp和mtxout俩中间变量，这将占用大量额外内存，有没有什么办法达到同样的效果而不占用更多内存的？或者Mathematica有没有类似MataLab中FFTShift的内置函数？
2018年10月29日 02点10分
1
代码如下：
FFTShift[mtx_] := Block[{dim = Dimensions@mtx},
midx = Round[dim[[1]]/2];
midy = Round[dim[[2]]/2];
mtxtmp = ConstantArray[0, dim];
mtxtmp[[All, 1 ;; midy]] = mtx[[All, -midy ;; -1]];
mtxtmp[[All, midy + 1 ;;]] =
mtx[[All, 1 ;; dim[[2]] - midy]];(*交换行*)
mtxout = ConstantArray[0, dim];
mtxout[[1 ;; midx, All]] = mtxtmp[[-midx ;; -1, All]];
mtxout[[midx + 1 ;;, All]] =
mtxtmp[[1 ;; dim[[1]] - midx, All]];(*交换列*)
mtxout];
然而可以看到以上函数其实就得定义mtxtmp和mtxout俩中间变量，这将占用大量额外内存，有没有什么办法达到同样的效果而不占用更多内存的？或者Mathematica有没有类似MataLab中FFTShift的内置函数？
