3-7.m

来自「《MATLAB 7.0编程基础》第3章 （源码实例）主要讲解MATLAB的数值计」· M 代码 · 共 19 行

A = [11    0   13    0
      0   22    0   24
      0    0   33    0
     41    0    0   44
      0   52    0    0
      0    0   63    0
      0    0    0   74]; 
% 创建稀疏矩阵A
[B,d] = spdiags(A)
% 提取对角线全元素矩阵B和位置向量d 

B1 = magic(4)
% 生成新的全元素矩阵B1

d1 = [-5,-2,0,2]; 
% 生成新的位置向量d1
A = spdiags(B1,d1,A);
% 将原稀疏矩阵相应位置的元素替换为新的元素，未替换的原元素不变
full(A)