注册
登录
矩阵交换机
使用切片/索引在矩阵内切换numpy向量
返回
使用切片/索引在矩阵内切换numpy向量
作者:
易米烊光
发布时间:
2025-03-03 05:29:26 (4小时前)
转自：
<p> <br/> 让我们考虑以下几点<br/> <br/> 矩阵<br/> </跨度><br/> ：</p><p>[v1 v2 v3 v4 v5 v6] =<br/><a href="./1.%EF%BC%8C0.%EF%BC%8C0.%EF%BC%8C0.%EF%BC%8C0.%EF%BC%8C0%5D%EF%BC%8C%0A%5B0.%EF%BC%8C0.%EF%BC%8C1.%EF%BC%8C0.%EF%BC%8C0.%EF%BC%8C0%E3%80%82%5D%EF%BC%8C%0A%5B0.%EF%BC%8C1.%EF%BC%8C0.%EF%BC%8C0.%EF%BC%8C0.%EF%BC%8C0%E3%80%82%5D%EF%BC%8C%0A%5B0.%EF%BC%8C0.%EF%BC%8C0.%EF%BC%8C1.%EF%BC%8C0.%EF%BC%8C0%E3%80%82%5D%EF%BC%8C%0A%5B0%2C0%2C0%2C0%2C0%2C0%2C1%5D%EF%BC%8C%0A%5B0%2C0%E3%80%82%EF%BC%8C0.%EF%BC%8C0.%EF%BC%8C1.%EF%BC%8C0%E3%80%82.html">1.，0.，0.，0.，0.，0]， [0.，0.，1.，0.，0.，0。]， [0.，1.，0.，0.，0.，0。]， [0.，0.，0.，1.，0.，0。]， [0,0,0,0,0,0,1]， [0,0。，0.，0.，1.，0。</a>，</p><p>其中向量[v2 v3]和[v5 v6]相对于身份进行切换<br/> <br/> 矩阵<br/> </跨度><br/> 。</p><p>从身份开始<br/> <br/> 矩阵<br/> </跨度><br/> ，我们怎样才能建立前一个<br/> <br/> 矩阵<br/> </跨度><br/> 使用切片鈥<br/> <br/></p>
收藏
举报
4 条回复
0#
回复此人
撩心
|
2019-08-31 10-32
<div class =“post-text”itemprop =“text”> <P> 没有一个函数（我知道）将交换列的位置，因为您可以在发布时使用基本切片来执行此操作。您始终可以构建该类型的函数： </p> <pre> <code> def swap(arr, *tups, copy=False): ix = np.array(tups) if copy: arr = arr.copy() arr[ix.ravel()] = arr[ix[:,::-1].ravel()] return arr swap(np.eye(6), (0,2), (1,3)) >> array([[0., 0., 1., 0., 0., 0.], [0., 0., 0., 1., 0., 0.], [1., 0., 0., 0., 0., 0.], [0., 1., 0., 0., 0., 0.], [0., 0., 0., 0., 1., 0.], [0., 0., 0., 0., 0., 1.]]) </code> </pre> <P> 如果你想要性能，我会建议使用稀疏矩阵。该 <code> csr_matrix </code> 存储非零列 <code> indices </code> ，使切换列容易： </p> <pre> <code> import scipy.sparse as sparse I = sparse.eye(6, format='csr') I.indices[[0,2,1,3]] = I.indices[[2,0,3,1]] print(I.todense()) >> array([[0., 0., 1., 0., 0., 0.], [0., 0., 0., 1., 0., 0.], [1., 0., 0., 0., 0., 0.], [0., 1., 0., 0., 0., 0.], [0., 0., 0., 0., 1., 0.], [0., 0., 0., 0., 0., 1.]]) </code> </pre> </DIV>
编辑
1#
回复此人
圈圈红
|
2019-08-31 10-32
<div class =“post-text”itemprop =“text”> <P> 矢量解决方案可扩展到任何 <code> N </code> ： </p> <pre> <code> N=2 s = np.zeros((3*N,3*N)) pattern = [[1,0,0],[0,0,1],[0,1,0]] s.reshape(N,3,N,3)[range(N),:,range(N),:] = pattern </code> </pre> <P> 对于 </p> <pre> <code> array([[1., 0., 0., 0., 0., 0.], [0., 0., 1., 0., 0., 0.], [0., 1., 0., 0., 0., 0.], [0., 0., 0., 1., 0., 0.], [0., 0., 0., 0., 0., 1.], [0., 0., 0., 0., 1., 0.]]) </code> </pre> </DIV>
编辑
2#
回复此人
不想吃东西
|
2019-08-31 10-32
<div class =“post-text”itemprop =“text”> <P> 像这样的东西？ </p> <pre class="lang-py prettyprint-override"> <code> # helper function, swaps two rows def swaprows(arr, row1, row2): temp = arr[row1].copy() arr[row1] = arr[row2] arr[row2] = temp return arr a = np.eye(6) a = swaprows(a, 1, 2) a = swaprows(a, 4, 5) print(a) # => # [[1. 0. 0. 0. 0. 0.] # [0. 0. 1. 0. 0. 0.] # [0. 1. 0. 0. 0. 0.] # [0. 0. 0. 1. 0. 0.] # [0. 0. 0. 0. 0. 1.] # [0. 0. 0. 0. 1. 0.]] </code> </pre> <P> 请注意，python使用0索引，我猜你来自R或Matlab，你的最后一个代码块将a）不使用双括号和b）你会因为尝试索引大小数组而得到索引超出范围的错误6位于第6位。 </p> </DIV>
编辑
登录
后才能参与评论
