怎样构造稀疏矩阵

allenfieldin

要再次麻烦各位大侠了，呵呵
最新用matlab编有限元，遇到如附件中的稀疏矩阵，其中A,B为3*3的方阵，上标T表示转置，想构造具有M个对角分块矩阵A的稀疏阵，该怎么构造？
如：A=

7.96E+13	0	0
0	9.62E+11	0
0	0	2.65E+13

B=

-4E+13	0	-2E+13
0	4.0E+11	0
1.99E+13	0	6.63E+12

wxm27811

这个需要自己写几行程序的吧。貌似没有相关函数。

nwcwww

这个问题不复杂，解决方法也不止一种。。用blkdiag()的话挺容易的。

方法之一是先构造一个m个A+1个A/2的对角线矩阵，再构造一个m个矩阵B组成的对角线矩阵并补0保证大小一样。然后矩阵1+矩阵2+矩阵2的转置就行了。

1. %定义
   
2. m=4;
   
3. Adiag=[2 4 6];     %A的对角线元素
   
4. B=[1,0,3;0,5,0;7,0,9];
   
5. 
   
6. 
   
7. %M1是由m个 A 和一个 A/2组成的稀疏矩阵
   
8. M1=diag([repmat(Adiag,1,m),Adiag/2]);
   
9. 
   
10. %M2由矩阵B组成
    
11. BB=cell(1,m);       %共计m个B矩阵
    
12. [BB{:}]=deal(B);
    
13. M2_0=blkdiag(BB{:});
    
14. M2=[zeros(3*m,3),M2_0;zeros(3,3*m+3)];   %补0
    
15. 
    
16. %M1，M2和M2的转置矩阵相加即可
    
17. M=M1+M2+M2'

复制代码

M =

     2     0     0     1     0     3     0     0     0     0     0     0     0     0     0
     0     4     0     0     5     0     0     0     0     0     0     0     0     0     0
     0     0     6     7     0     9     0     0     0     0     0     0     0     0     0
     1     0     7     2     0     0     1     0     3     0     0     0     0     0     0
     0     5     0     0     4     0     0     5     0     0     0     0     0     0     0
     3     0     9     0     0     6     7     0     9     0     0     0     0     0     0
     0     0     0     1     0     7     2     0     0     1     0     3     0     0     0
     0     0     0     0     5     0     0     4     0     0     5     0     0     0     0
     0     0     0     3     0     9     0     0     6     7     0     9     0     0     0
     0     0     0     0     0     0     1     0     7     2     0     0     1     0     3
     0     0     0     0     0     0     0     5     0     0     4     0     0     5     0
     0     0     0     0     0     0     3     0     9     0     0     6     7     0     9
     0     0     0     0     0     0     0     0     0     1     0     7     1     0     0
     0     0     0     0     0     0     0     0     0     0     5     0     0     2     0
     0     0     0     0     0     0     0     0     0     3     0     9     0     0     3

allenfieldin

3# nwcwww
nwcwww 起的真早，我才刚刚看到，不过程序挺好的，我得看看cell的用法了，消化一下deal和bikdiag，呵呵

qibbxxt

3# nwcwww
不错，用稀疏矩阵的办法也能实现

1. 
   
2. m=4;
   
3. Adiag=[2 4 6];     %A的对角线元素
   
4. B=[1,0,3;0,5,0;7,0,9];
   
5. md = sparse(1:m+1,1:m+1,[ones(1,m),1/2],m+1,m+1);
   
6. M = kron(md,spdiags(Adiag',0,length(Adiag),length(Adiag)))+...
   
7.     kron(spdiags(ones(m+1,2),1,m+1,m+1),B)+...
   
8.     kron(spdiags(ones(m+1,2),-1,m+1,m+1),B');

复制代码

1. 
   
2. >> full(M)
   

复制代码

ans =
     2     0     0     1     0     3     0     0     0     0     0     0     0     0     0
     0     4     0     0     5     0     0     0     0     0     0     0     0     0     0
     0     0     6     7     0     9     0     0     0     0     0     0     0     0     0
     1     0     7     2     0     0     1     0     3     0     0     0     0     0     0
     0     5     0     0     4     0     0     5     0     0     0     0     0     0     0
     3     0     9     0     0     6     7     0     9     0     0     0     0     0     0
     0     0     0     1     0     7     2     0     0     1     0     3     0     0     0
     0     0     0     0     5     0     0     4     0     0     5     0     0     0     0
     0     0     0     3     0     9     0     0     6     7     0     9     0     0     0
     0     0     0     0     0     0     1     0     7     2     0     0     1     0     3
     0     0     0     0     0     0     0     5     0     0     4     0     0     5     0
     0     0     0     0     0     0     3     0     9     0     0     6     7     0     9
     0     0     0     0     0     0     0     0     0     1     0     7     1     0     0
     0     0     0     0     0     0     0     0     0     0     5     0     0     2     0
     0     0     0     0     0     0     0     0     0     3     0     9     0     0     3

allenfieldin

5# qibbxxt
版主太强大了，以前用软件计算时，总有一个Sparse组装，版主竟然也用上了，呵呵长了见识:lol

rend

5# qibbxxt
长见识了