假设矩阵为a
判断矩阵中是否有nan
display([num2str(sum(sum(1*(isnan(a)))))]);
找出nan的位置
ind_nan = find(isnan(a)); % 输出nan的线性索引
sub_nan = zeros(length(ind_nan), 2); % 行列式矩阵位置索引
for n=1:length(ind_nan)
[sub_nan(n,1), sub_nan(n,2)] = ind2sub(size(a),ind_nan(n));
end