python关于矩阵重复赋值覆盖问题的解决方法

作者:古德蜡克 时间:2021-12-11 03:07:17 

本文实例讲述了python关于矩阵重复赋值覆盖问题的解决方法。分享给大家供大家参考,具体如下:


import itertools
import numpy as np
comb = list(itertools.combinations(list(range(regions)), 2))
bands_info = []
coeff = np.zeros([bands, len(comb)])
for cla in range(classes):
 class_info = data[:,cla*bands*regions:(cla+1)*bands*regions]
 for bs in range(bands):
   n = bs*regions
   for i in range(len(comb)):
     index1 = comb[i][0]+n
     index2 = comb[i][1]+n
     part1 = class_info[:, index1]
     part2 = class_info[:, index2]
     coeff[bs, i] = (np.corrcoef(part1, part2))[0, 1]
 bands_info.append(coeff.reshape([1,-1]))
coeff_info = np.vstack((bands_info[0], bands_info[1], bands_info[2],bands_info[3]))

例如这个循环赋值过程,最终得出来的结果是bands_info这个List里面每一个矩阵都是一样的,这是为什么呢?我一开始也在这里纠结了很长时间,思来想去感觉没错的呀。后来想想以前学的C语言知识,才有点明白。原来python里面有浅层copy和深层copy这一说,同是一个矩阵的话占用的是同一个地址,在里面进行重复赋值的话前面的值都会被覆盖掉。不只是当前变量被覆盖掉,就是你之后用到这个变量的也会被覆盖。比如说你a的变量被b覆盖了,那你后面用到a的变量的地方,a的值也会变成b的值。是不是很可怕。

那么应该怎么进行修改呢?


import itertools
import numpy as np
comb = list(itertools.combinations(list(range(regions)), 2))
bands_info = []
for cla in range(classes):
 coeff = np.zeros([bands, len(comb)])
 class_info = data[:,cla*bands*regions:(cla+1)*bands*regions]
 for bs in range(bands):
   n = bs*regions
   for i in range(len(comb)):
     index1 = comb[i][0]+n
     index2 = comb[i][1]+n
     part1 = class_info[:, index1]
     part2 = class_info[:, index2]
     coeff[bs, i] = (np.corrcoef(part1, part2))[0, 1]
 bands_info.append(coeff.reshape([1,-1]))
coeff_info = np.vstack((bands_info[0], bands_info[1], bands_info[2],bands_info[3]))

改成这样就可以了。就是将初始矩阵在循环里在重新设定一遍,意思相当就是给了一次新的地址,再进行复制的话就不会覆盖前面的结果了。

希望本文所述对大家Python程序设计有所帮助。

来源:https://www.jianshu.com/p/472b7f615256

标签:python,矩阵
0
投稿

猜你喜欢

  • 基于Python绘制一个摸鱼倒计时界面

    2022-09-02 03:03:26
  • document.execCommand() 解析

    2007-11-13 17:22:00
  • python中出现invalid syntax报错的几种原因分析

    2023-01-26 22:01:24
  • Python基于pygame实现的font游戏字体(附源码)

    2021-04-16 05:06:17
  • django settings.py 配置文件及介绍

    2022-01-17 12:52:44
  • 20行Python代码实现一款永久免费PDF编辑工具

    2023-11-17 23:51:45
  • 防采集,几种觉得有用的防采集方法

    2009-09-03 13:30:00
  • 关于生成目录树结构的类

    2007-09-13 12:19:00
  • python使用rsa非对称加密过程解析

    2021-06-15 00:14:51
  • 如何基于Python实现word文档重新排版

    2023-06-09 11:19:24
  • MySQL错误中文参照列表

    2010-09-30 14:41:00
  • Python 解析pymysql模块操作数据库的方法

    2021-05-07 07:46:13
  • 解决numpy矩阵相减出现的负值自动转正值的问题

    2021-10-07 19:49:16
  • js实现GIF动图分解成多帧图片上传

    2023-08-15 00:06:38
  • python3.6使用pymysql连接Mysql数据库

    2024-01-27 13:00:48
  • AJAX缓存问题的两种解决方法(IE)

    2008-05-02 20:57:00
  • 分享下GET和POST的真正区别

    2023-10-19 23:23:09
  • 解决Python3中的中文字符编码的问题

    2022-05-04 04:22:09
  • 使用c#构造date数据类型

    2024-01-15 22:19:15
  • python如何生成任意n阶的三对角矩阵

    2021-09-18 07:41:37
  • asp之家 网络编程 m.aspxhome.com