Question

关于MATLAB的一个编程 请高手解释一下

r=1;R=6;v=18;n=10;
b1=sparse(1,1,v,n,1);
b=full(b1);
a1=triu(r*ones(n.n));
a2=diag(R*ones(1,n));
a3=-tril(R*ones(n,n),-1)+tril(R*ones(n,n)),-2);
a=a1+a2+a3;
I=a\b;
I0=sum(I)
画出流程图也可以 隔行如隔山啊

Answer 1

r=1;R=6;v=18;n=10;%对变量进行赋初值
b1=sparse(1,1,v,n,1);%利用1,1,v,产生一个n*1的稀疏矩阵
b=full(b1);%把系数矩阵b1转换为满阵
a1=triu(r*ones(n,n));%a1是10*10的上三角阵，上三角元素都为1
a2=diag(R*ones(1,n));%a2是对角阵，主对角元素是6
a3=-tril(R*ones(n,n),-1)+tril(R*ones(n,n),-2);%a3矩阵的主对角线下方第一个对角线的元素是-6
a=a1+a2+a3;%三个矩阵相加
I=a\b;%a阵除以b阵
I0=sum(I) %对I求和

这段程序主要的目的是求I0