一直调试的问题是：Error using ==> mtimes Inner matrix dimensions must

A=[1 0。5 0;0 1 0。5;0。5 0 1];D=diag(diag(A)); L=-tril(A,-1); U=-triu(A,1);for i=1:61 w(i)=(i-1)*0。1-2; B=(inv(D-L*w)*((1-w)*D+w*U)); [v,d]=eigs(B); c(i)=mytest(d); % 确定mytest 输出的是个数字endplot(w,c); 再问： 问题是出在B=(inv(D-L*w)*((1-w)*D+w*U));这一行，这一行只要出现w，就会报错，是不是因为w是个未知数，不能出现在里面？你能帮我调试一下吗？谢谢了。肯定加分 再答： 对的，你又没定义函数，没有定义w，怎么能用呢 你的c想表示什么呢？c是什么含义 看着是超松弛和欠松弛 我班你改的已经改掉了w的问题再问： 但是我还是报错额，你们吧得到的结果截图发给我嘛？下午要交作业额 再答： 我问你，你的c是啥？再问： A=[1 0。5 0;0 1 0。5;0。5 0 1]; D=diag(diag(A)); L=-tril(A,-1); U=-triu(A,1); for i=1:61 w(i)=(i-1)*0。1-2; B=(inv(D-L*w)*((1-w)*D+w*U)); d=eigs(B); c=max(d); % 确定mytest 输出的是个数字 end plot(w,c); 我是要得到c关于w的图像。 再答： d是个矩阵 c=max(d)是个向量 不对的，我就问你，你的c表示的是什么意思，不要给我code。。。再问： c表示的是B的特征值的绝对值的最大值 再答： OK，搞定咯 A=[1 0。5 0;0 1 0。5;0。5 0 1]; D=diag(diag(A)) L=-tril(A,-1) U=-triu(A,1) for i=1:61 w(i)=(i-1)*0。1-2; B=inv(D-L*w(i))*((1-w(i))*D+w(i)*U); [v,d]=eigs(B); c(i)=max(max(abs(d))); end plot(w,c); xlabel('w') ylabel('max eig val')