for ii = 1:2:2*de.nP
G=[one, aux3(:,ii), aux3(:,ii) - aux2(:,ii),aux3(:,ii+1) - aux2(:,ii+1)];
try
betasPu(:,(ii+1)/2) =G\yM;
catch ME
betasPu(:,(ii+1)/2)=[Inf,Inf,Inf,Inf];
end
end
我正在使用上面的代码尝试在 G\yM 中使用 G 时捕获实例;条件很差/单数。
我已经使用了上面的代码,但现在尽管产生了几个条件不好/奇异的矩阵,但现在似乎没有触发捕获。