0

我有这样的代码：

[X, Y, Z]= meshgrid(0:5:50, 0:5:50, 0:3:3);
[t1, t2, t3]= meshgrid(35, 25, 3);

A=sqrt((X - t1).^2 + (Y - t2).^2 + (Z - t3).^2);

[rows1,cols1,pages1] = ind2sub(size(A), find(A<=5));

for i=1:1:length(rows1);
    [R1, R2, R3]= meshgrid(cols1(i)*5-5, rows1(i)*5-5, pages1(i)*3-3);
    [reader]=[R1, R2, R3]

    B=sqrt((X1 - R1).^2 + (Y1 - R2).^2 + (Z1 - R3).^2);
    [rowsB,colsB,pagesB] = ind2sub(size(B), find(B<=5));

    for j=1:1:length(rowsB);
        [T1, T2, T3]= meshgrid(colsB(j)*0.2-0.2+0.15, rowsB(j)*0.2-0.2+0.15, pagesB(j)*3-3);
        [tag]=[T1, T2, T3];
    end
end

我需要从“for”循环中获得所有结果，因为需要对这些结果应用“intersec”。但是我遇到了问题，因为“for”会覆盖结果，而我只取最后一个。如何使某事喜欢

[tag(i)]=[T1, T2, T3] where i is i=1:1:length(rows1).

我已经通过预分配应用了 smth，但没有工作。

4

2 回答 2

2

一般来说，您可以使用单元格数组：

tag = cell(length(rows1),1);
for i=1:length(rows1)
    tag{i} = [T1, T2, T3];
end

虽然看起来像一个简单的矩阵就可以了..

于 2012-07-31T15:04:25.570 回答
0

为了将向量分配给矩阵的行，您需要使用冒号作为列索引对列进行切片。像这样的东西：

tag = zeros(length(rows1), 3);
for i = 1:length(rows1)
    ...
    tag(i, :) = [T1, T2, T3];
end

编辑：修正了我的例子。我误读了你的代码。以上较为笼统。

于 2012-07-31T15:15:13.913 回答