稀疏矩阵数值分析及MATLAB实现示例.pdfVIP

clearall;clc

n=8;

p=ones(n,1);%产生n维列向量，每项都为1

sp2=spdiags([-p/4,-p/2,p*3,-p/2,-p/4],-2:2,n,n);

%其中-2：2表示参数为-2，-1,0,1,2，其意思代表下次二对角线，下次对角线，主对角线，上次对角线，上次二对角线

sf=sym(full(sp2))

%sym(diag(linspace(-1/4,-1/4,6),2)+diag(linspace(-1/2,-1/2,7),1)+diag(linspace(3,3,8))+diag(linspace(-1/4,-

1/4,6),-2)+diag(linspace(-1/2,-1/2,7),-1))

%clearall;clc;n=8;p=ones(n,1);sp2=spdiags([-p/4,-p/2,p*3,-p/2,-p/4],-2:2,n,n);sf=sym(full(sp2))

%sym(diag(diag(-1/4*ones(6)),2)+diag(diag(-1/2*ones(7)),1)+diag(linspace(3,3,8))+diag(linspace(-1/4,-

1/4,6),-2)+diag(linspace(-1/2,-1/2,7),-1))