D=diag(A);
U=triu(A,1);
L=tril(A,-1);
G=-(D+L)\U; f=(D+L)\b;
运行到红色那一句的时候就出现Matrix dimensions must agree.错误,哪位高手能告诉我为什么呢?
------解决方案--------------------------------------------------------
A=ones(3,3);
D=diag(A,3)
U=triu(A,1);
L=tril(A,-1)
G=-(D+L)\U;
f=(D+L)\b;
你把D和L打印出来就知道了,不一致!
D =
1
1
1
L =
0 0 0
1 0 0
1 1 0
没法做运算!