一直调试的问题是：Error using ==> mtimes Inner matrix dimensions must agree.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); B=(inv(D-L*w)*((1-w)*D+w*U));[v,d]=eigs(B);w=-2:0.1:4;c=myTest(d);plot(w,c);A=[1 0.5 0;0 1

一直调试的问题是：Error using ==> mtimes Inner matrix dimensions must agree.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); B=(inv(D-L*w)*((1-w)*D+w*U));[v,d]=eigs(B);w=-2:0.1:4;c=myTest(d);plot(w,c);A=[1 0.5 0;0 1
一直调试的问题是：Error using ==> mtimes Inner matrix dimensions must agree.
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);
B=(inv(D-L*w)*((1-w)*D+w*U));
[v,d]=eigs(B);
w=-2:0.1:4;
c=myTest(d);
plot(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);
B=(inv(D-L*w)*((1-w)*D+w*U));
[v,d]=eigs(B);
w=-2:0.1:4;
c=max(d);
plot(w,c);
Error using ==> mtimes
Inner matrix dimensions must agree.

一直调试的问题是：Error using ==> mtimes Inner matrix dimensions must agree.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); B=(inv(D-L*w)*((1-w)*D+w*U));[v,d]=eigs(B);w=-2:0.1:4;c=myTest(d);plot(w,c);A=[1 0.5 0;0 1
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 输出的是个数字
end
plot(w,c);