假设矩阵为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