问题描述
分别实现矩阵相乘的3种算法,比较三种算法在矩阵大小分别为22"htmlcode">
#-*- coding: utf-8 -*- from matplotlib.font_manager import FontProperties import numpy as np import time import random import math import copy import matplotlib.pyplot as plt #n = [2**2, 2**3, 2**4, 2**5, 2**6, 2**7, 2**8, 2**9, 2**10, 2**11, 2**12] n = [2**2, 2**3, 2**4, 2**5, 2**6, 2**7, 2**8, 2**9, 2**10, 2**11] Sum_time1 = [] Sum_time2 = [] Sum_time3 = [] Sum_time4 = [] for m in n: A = np.random.randint(0, 2, [m, m]) B = np.random.randint(0, 2, [m, m]) A1 = np.mat(A) B1 = np.mat(B) time_start = time.time() C1 = A1*B1 time_end = time.time() Sum_time1.append(time_end - time_start) C2 = np.zeros([m, m], dtype = np.int) time_start = time.time() for i in range(m): for k in range(m): for j in range(m): C2[i, j] = C2[i, j] + A[i, k] * B[k, j] time_end = time.time() Sum_time2.append(time_end - time_start) A11 = np.mat(A[0:m//2, 0:m//2]) A12 = np.mat(A[0:m//2, m//2:m]) A21 = np.mat(A[m//2:m, 0:m//2]) A22 = np.mat(A[m//2:m, m//2:m]) B11 = np.mat(B[0:m//2, 0:m//2]) B12 = np.mat(B[0:m//2, m//2:m]) B21 = np.mat(B[m//2:m, 0:m//2]) B22 = np.mat(B[m//2:m, m//2:m]) time_start = time.time() C11 = A11 * B11 + A12 * B21 C12 = A11 * B12 + A12 * B22 C21 = A21 * B11 + A22 * B21 C22 = A21 * B12 + A22 * B22 C3 = np.vstack((np.hstack((C11, C12)), np.hstack((C21, C22)))) time_end = time.time() Sum_time3.append(time_end - time_start) time_start = time.time() M1 = A11 * (B12 - B22) M2 = (A11 + A12) * B22 M3 = (A21 + A22) * B11 M4 = A22 * (B21 - B11) M5 = (A11 + A22) * (B11 + B22) M6 = (A12 - A22) * (B21 + B22) M7 = (A11 - A21) * (B11 + B12) C11 = M5 + M4 - M2 + M6 C12 = M1 + M2 C21 = M3 + M4 C22 = M5 + M1 - M3 - M7 C4 = np.vstack((np.hstack((C11, C12)), np.hstack((C21, C22)))) time_end = time.time() Sum_time4.append(time_end - time_start) f1 = open('python_time1.txt', 'w') for ele in Sum_time1: f1.writelines(str(ele) + '\n') f1.close() f2 = open('python_time2.txt', 'w') for ele in Sum_time2: f2.writelines(str(ele) + '\n') f2.close() f3 = open('python_time3.txt', 'w') for ele in Sum_time3: f3.writelines(str(ele) + '\n') f3.close() f4 = open('python_time4.txt', 'w') for ele in Sum_time4: f4.writelines(str(ele) + '\n') f4.close() font = FontProperties(fname=r"c:\windows\fonts\simsun.ttc", size=8) plt.figure(1) plt.subplot(221) plt.semilogx(n, Sum_time1, 'r-*') plt.ylabel(u"时间(s)", fontproperties=font) plt.xlabel(u"矩阵的维度n", fontproperties=font) plt.title(u'python自带的方法', fontproperties=font) plt.subplot(222) plt.semilogx(n, Sum_time2, 'b-*') plt.ylabel(u"时间(s)", fontproperties=font) plt.xlabel(u"矩阵的维度n", fontproperties=font) plt.title(u'定义法', fontproperties=font) plt.subplot(223) plt.semilogx(n, Sum_time3, 'y-*') plt.ylabel(u"时间(s)", fontproperties=font) plt.xlabel(u"矩阵的维度n", fontproperties=font) plt.title( u'分治法', fontproperties=font) plt.subplot(224) plt.semilogx(n, Sum_time4, 'g-*') plt.ylabel(u"时间(s)", fontproperties=font) plt.xlabel(u"矩阵的维度n", fontproperties=font) plt.title( u'Strasses法', fontproperties=font) plt.figure(2) plt.semilogx(n, Sum_time1, 'r-*', n, Sum_time2, 'b-+', n, Sum_time3, 'y-o', n, Sum_time4, 'g-^') #plt.legend(u'python自带的方法', u'定义法', u'分治法', u'Strasses法', fontproperties=font) plt.show()
以上这篇Python实现矩阵相乘的三种方法小结就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
标签:
Python,矩阵,相乘
免责声明:本站文章均来自网站采集或用户投稿,网站不提供任何软件下载或自行开发的软件!
如有用户或公司发现本站内容信息存在侵权行为,请邮件告知! 858582#qq.com
白云城资源网 Copyright www.dyhadc.com
暂无“Python实现矩阵相乘的三种方法小结”评论...
《魔兽世界》大逃杀!60人新游玩模式《强袭风暴》3月21日上线
暴雪近日发布了《魔兽世界》10.2.6 更新内容,新游玩模式《强袭风暴》即将于3月21 日在亚服上线,届时玩家将前往阿拉希高地展开一场 60 人大逃杀对战。
艾泽拉斯的冒险者已经征服了艾泽拉斯的大地及遥远的彼岸。他们在对抗世界上最致命的敌人时展现出过人的手腕,并且成功阻止终结宇宙等级的威胁。当他们在为即将于《魔兽世界》资料片《地心之战》中来袭的萨拉塔斯势力做战斗准备时,他们还需要在熟悉的阿拉希高地面对一个全新的敌人──那就是彼此。在《巨龙崛起》10.2.6 更新的《强袭风暴》中,玩家将会进入一个全新的海盗主题大逃杀式限时活动,其中包含极高的风险和史诗级的奖励。
《强袭风暴》不是普通的战场,作为一个独立于主游戏之外的活动,玩家可以用大逃杀的风格来体验《魔兽世界》,不分职业、不分装备(除了你在赛局中捡到的),光是技巧和战略的强弱之分就能决定出谁才是能坚持到最后的赢家。本次活动将会开放单人和双人模式,玩家在加入海盗主题的预赛大厅区域前,可以从强袭风暴角色画面新增好友。游玩游戏将可以累计名望轨迹,《巨龙崛起》和《魔兽世界:巫妖王之怒 经典版》的玩家都可以获得奖励。
更新日志
2025年01月12日
2025年01月12日
- 小骆驼-《草原狼2(蓝光CD)》[原抓WAV+CUE]
- 群星《欢迎来到我身边 电影原声专辑》[320K/MP3][105.02MB]
- 群星《欢迎来到我身边 电影原声专辑》[FLAC/分轨][480.9MB]
- 雷婷《梦里蓝天HQⅡ》 2023头版限量编号低速原抓[WAV+CUE][463M]
- 群星《2024好听新歌42》AI调整音效【WAV分轨】
- 王思雨-《思念陪着鸿雁飞》WAV
- 王思雨《喜马拉雅HQ》头版限量编号[WAV+CUE]
- 李健《无时无刻》[WAV+CUE][590M]
- 陈奕迅《酝酿》[WAV分轨][502M]
- 卓依婷《化蝶》2CD[WAV+CUE][1.1G]
- 群星《吉他王(黑胶CD)》[WAV+CUE]
- 齐秦《穿乐(穿越)》[WAV+CUE]
- 发烧珍品《数位CD音响测试-动向效果(九)》【WAV+CUE】
- 邝美云《邝美云精装歌集》[DSF][1.6G]
- 吕方《爱一回伤一回》[WAV+CUE][454M]