码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • Python实现Matlab中的circshift函数


    circshift是Matlab中矩阵循环移位函数,具体使用参照该链接。
    但是python中并没有封装好的该函数,因此需要自己实现。

    思路:将矩阵分为两部分,然后按照自己的需要堆叠在一起就可以了。
    np.vstack用在行的分割;
    np.hstack用在列的分割。

    针对向量
    ##  u代表原向量,shiftnum2代表移动的单位。
    def circshift(u,shiftnum2):
        w=len(u)
        if shiftnum2 > 0:
            u = np.hstack((u[(w - shiftnum2):], u[:(w - shiftnum2)]))
        else:
            u = np.hstack((u[-shiftnum2:],u[:-shiftnum2]))
        return u
    X=np.zeros(10)
    X[0]=1
    X[1]=2
    X[2]=1
    p=X/2
    print(p)
    print(circshift(p,2))
    
    [0.5 1.  0.5 0.  0.  0.  0.  0.  0.  0. ]
    [0.  0.  0.5 1.  0.5 0.  0.  0.  0.  0. ]
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9
    • 10
    • 11
    • 12
    • 13
    • 14
    • 15
    • 16
    • 17
    • 18
    针对矩阵
    ##   u代表原矩阵,shiftnum1代表行,shiftnum2代表列。
    def circshift(u,shiftnum1,shiftnum2):
        h,w = u.shape
        if shiftnum1 < 0:
            u = np.vstack((u[-shiftnum1:,:],u[:-shiftnum1,:]))
        else:
            u = np.vstack((u[(h-shiftnum1):,:],u[:(h-shiftnum1),:]))
        if shiftnum2 > 0:
            u = np.hstack((u[:, (w - shiftnum2):], u[:, :(w - shiftnum2)]))
        else:
            u = np.hstack((u[:,-shiftnum2:],u[:,:-shiftnum2]))
        return u
    
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9
    • 10
    • 11
    • 12
    • 13

    若 shiftnum1大于0,则代表从上往下移动,否则就是从下往上移动;
    若 shiftnum2大于0,则代表从左往右移动,否则就是从右往左移动;

    示例1:

    a = np.matrix([[1,2,3],[4,5,6],[7,8,9]])
    b = circshift(a,-1,0)
    print(b)
    
    输出:
    [[4 5 6]
     [7 8 9]
     [1 2 3]]
    
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9

    示例2:

    a = np.matrix([[1,2,3],[4,5,6],[7,8,9]])
    b = circshift(a,1,0)
    print(b)
    
    输出:
    [[7 8 9]
     [1 2 3]
     [4 5 6]]
    
    
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9
    • 10

    示例3:

    a = np.matrix([[1,2,3],[4,5,6],[7,8,9]])
    b = circshift(a,0,1)
    print(b)
    
    输出:
    [[3 1 2]
     [6 4 5]
     [9 7 8]]
    
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9

    示例4:

    a = np.matrix([[1,2,3],[4,5,6],[7,8,9]])
    b = circshift(a,0,-1)
    print(b)
    
    输出:
    [[2 3 1]
     [5 6 4]
     [8 9 7]]
    
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7
    • 8
    • 9
  • 相关阅读:
    13个培训师游戏
    STM8的C语言编程(12)--+AD转换
    C#,数值计算——函数计算,切比雪夫近似算法(Chebyshev approximation)的计算方法与源程序
    【内网安全】——Linux信息收集
    小样本规模船型优化策略的选择研究
    【RocketMQ集群】Linux搭建RocketMQ双主双从集群
    淘宝API接口参考
    线性代数 | 矩阵运算 加减 数乘 矩阵的幂运算
    含镧系配合物荧光/有机硅烷/二氧化硅/二氧化硅气凝胶聚苯乙烯微球改性与制备
    【学习笔记41】DOM操作的练习
  • 原文地址:https://blog.csdn.net/y15520833229/article/details/136235240
  • 最新文章
  • 攻防演习之三天拿下官网站群
    数据安全治理学习——前期安全规划和安全管理体系建设
    企业安全 | 企业内一次钓鱼演练准备过程
    内网渗透测试 | Kerberos协议及其部分攻击手法
    0day的产生 | 不懂代码的"代码审计"
    安装scrcpy-client模块av模块异常,环境问题解决方案
    leetcode hot100【LeetCode 279. 完全平方数】java实现
    OpenWrt下安装Mosquitto
    AnatoMask论文汇总
    【AI日记】24.11.01 LangChain、openai api和github copilot
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1
正则表达式工具 cron表达式工具 密码生成工具

京公网安备 11010502049817号