3

我想对不属于特定范围的数组内容执行度量。

例如,我有一个包含 1000 行和 2 列的数组。我想对一列(假设第 2 列)中不属于第 50-150、250-300、400-700 和 900-950 行的所有元素执行 mean() 计算。

因此,应根据第 1-49、151-249、301-399、701-899 和 951-1000 行计算平均值。

任何想法如何去做?

编辑:我应该指出,每次运行程序时,包含的那些项目都会改变。因此,我不能只硬编码包含的内容;它们需要根据排除项来制定。

4

2 回答 2

6

怎么样:

M = rand(1000,2);
idx = setdiff(1:size(M,1), [50:150, 250:300, 400:700, 900:950]);

MM = M(idx,:)

现在将任何函数应用于过滤后的矩阵:

mean(MM,1)
于 2012-06-15T16:31:55.857 回答
0

您可以定义排除范围,然后使用逻辑寻址:

LowerLimit1 = 1;
UpperLimit1 = 50;

LowerLimit2 = 151;
UpperLimit2 = 249;

LowerLimit3 = 301;
UpperLimit3 = 399;

LowerLimit4 = 701;
UpperLimit4 = 899;

LowerLimit5 = 951;
UpperLimit5 = 1000;

MyVector = MyMatrix(:,2);

MeanValue = mean(MyVector(~(MyVector > LowerLimit1 & MyVector < UpperLimit1) | (MyVector > LowerLimit2 & MyVector < UpperLimit2) |  (MyVector > LowerLimit3 & MyVector < UpperLimit3)  | (MyVector > LowerLimit4 & MyVector < UpperLimit4) | (MyVector > LowerLimit5 & MyVector < UpperLimit5)));
于 2012-06-15T16:43:51.780 回答