|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
5 Z6 W2 U4 M9 g P' e) G3 E; H( _
一、稀疏矩阵
8 }8 n! J, g$ n( T0 Z! E
5 ]+ V s2 i& ^对于一个 n 阶矩阵,通常需要 n2 的存储空间,当 n 很大时,进行矩阵运算时会占用大量的内存空间和运算时间。在许多实际问题中遇到的大规模矩阵中通常含有大量0元素,这样的矩阵称为稀疏矩阵。Matlab支持稀疏矩阵,只存储矩阵的非零元素。由于不存储那些”0″元素,也不对它们进行操作,从而节省内存空间和计算时间,其计算的复杂性和代价仅仅取决于稀疏矩阵的非零元素的个数,这在矩阵的存储空间和计算时间上都有很大的优点。1 r' _$ X' D+ P% B6 q% B
矩阵的密度定义为矩阵中非零元素的个数除以矩阵中总的元素个数。对于低密度的矩阵,采用稀疏方式存储是一种很好的选择。 I7 L) @+ I+ x5 B* n6 Y
' {$ u- f* i. y1、稀疏矩阵的创建
# `( u: ~9 `7 p/ m(1) 将完全存储方式转化为稀疏存储方式 函数A=sparse(S)将矩阵S转化为稀疏存储方式的矩阵A。当矩阵S是稀疏存储方式时,则函数调用相当于A=S。 sparse函数还有其他一些调用格式: sparse(m,n):生成一个m*n的所有元素都是0的稀疏矩阵。 sparse(u,v,S)--:u,v,S是3个等长的向量。S是要建立的稀疏矩阵的非0元素,u(i)、v(i)分别是S(i)的行和列下标,该函数建立一个max(u)行、max(v)列并以S为稀疏元素的稀疏矩阵。 此外,还有一些和稀疏矩阵操作有关的函数。full(A):返回和稀疏存储矩阵A对应的完全存储方式矩阵。: G* X: e- j2 p* o6 c
(2) 直接创建稀疏矩阵 S=sparse(i,j,s,m,n),其中i 和j 分别是矩阵非零元素的行和列指标向量,s 是非零元素值向量,m,n 分别是矩阵的行数和列数。 {4 M- `- C: U5 m( n
(3) 从文件中创建稀疏矩阵 利用load和spconvert函数可以从包含一系列下标和非零元素的文本文件中输入稀疏矩阵。例:设文本文件 T.txt 中有三列内容\begin{bmatrix}1\; 3 \; 5\\ 2 \; 4 \; 6\\ 2 \; 5 \; 8\\ 3 \; 6 \; 9\end{bmatrix},第一列是一些行下标,第二列是列下标,第三列是非零元素值。load T.txt S=spconvert(T)。
8 V8 R2 T4 ~) ^# G7 X% G" A(4) 稀疏带状矩阵的创建 S=spdiags(B,d,m,n) 其中m 和n 分别是矩阵的行数和列数;d是长度为p的整数向量,它指定矩阵S的对角线位置;B是全元素矩阵,用来给定S对角线位置上的元素,行数为min(m,n),列数为p 。( q) E" A5 X/ s9 x
(5) 其它稀疏矩阵创建函数 @6 h. H5 p9 L( {, ?/ \
S=speye(m,n)
9 r( s0 s) {6 p1 F3 Y A1 t7 BS=speye(size(A)) % has the same size as A
3 _. J' c: V( US=buchy % 一个内置的稀疏矩阵(邻接矩阵)
4 T6 B6 x9 g2 ~等等
8 \$ s7 v& S/ g- u' Z6 x0 {1 a0 l8 ^; z* d. u' j3 V2 {
2、稀疏矩阵的运算
6 W. e9 H) J0 M8 N$ v, I2 h# C; Q6 A+ r2 U0 K; e
稀疏存储矩阵只是矩阵的存储方式不同,它的运算规则与普通矩阵是一样的,可以直接参与运算。所以,Matlab中对满矩阵的运算和函数同样可用在稀疏矩阵中。结果是稀疏矩阵还是满矩阵,取决于运算符或者函数。当参与运算的对象不全是稀疏存储矩阵时,所得结果一般是完全存储形式。
1 D+ C1 o/ K! x0 z! W' S9 \$ v2 ?
% a( N( R- K8 z5 i1 N! h8 \3、其他; H8 N! f. G8 F" d8 E: {. D
0 c5 [: M! h) R
(1) 非零元素信息+ S7 k" W9 V( \! [. U+ U# a, ^* K
nnz(S) % 返回非零元素的个数
! S' z1 f! L. ]% c4 e" }nonzeros(S) % 返回列向量,包含所有的非零元素, ]& w' }. u7 N4 I7 |: h0 n* c, d$ r
nzmax(S) % 返回分配给稀疏矩阵中非零项的总的存储空间
+ C6 B& f/ w( @3 v( V! t(2) 查看稀疏矩阵的形状 spy(S)9 e* Z; U \* H
(3) find函数与稀疏矩阵( P# F$ Z/ F* F5 e2 j
[i,j,s]=find(S)/ V& q7 s3 s1 g# F" Z
[i,j]=find(S)# z6 I& z* B3 g6 Z6 {
返回 S 中所有非零元素的下标和数值,S 可以是稀疏矩阵或满矩阵。
" K8 l, d+ L* p
9 T+ Z) [+ s, D$ ~- ^, O% N二、有限域中的矩阵
2 @# K: m3 Q- E, }
) P0 i) X$ x+ a; ~8 m. j信道编码中的矩阵运算一般都是基于有限域的,因此需要将普通矩阵转换为有限域中的矩阵,使其运算在有限域GF(m)中。可以通过命令gf(data,m)将数据限制在有限域中,这样如矩阵求逆、相加、相乘等运算就均是基于有限域GF(m)的运算了。
6 Q' T& U9 l& d1 c8 }
% v5 [4 f4 q3 x% I- o+ W9 n1 U4 B2 ^那么如何将有限域元素转换为double型的呢?可以利用命令 double(data.x) 其中x是后缀。关于有限域的详细情况请参考 这里。# ^0 B# f+ ^$ }& s U! \: d5 w- o
+ a% b3 y* x) M9 H/ b$ E" ]
1 [8 R* _* N- P Y7 ~6 n/ R7 ~
解决方法:用\;代替&。估计这个问题是Latex Math插件的bug。呵呵,不知道有没有更好的解决办法。 |
|