Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
高级会员
注册日期: 2019-11-21
帖子: 3,006
声望力: 66 ![]() |
![]()
我在MATLAB中的find函数遇到一些问题。我有一个由零和一组成的矩阵(代表结构元素的几何形状),其中当矩阵元素= 1时存在材料,而当矩阵元素= 0时不存在材料。矩阵可以具有一般形式如下图所示(它会随着几何形状的变化而更新,但这并不是太重要)。
Geometry = [0 0 0 0 0 0 0 0 0 0; 0 0 1 0 1 0 1 1 0 0; 0 0 1 0 0 0 0 1 0 0; 0 0 1 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 1 0 0; 0 0 0 0 0 0 0 0 0 0; 0 0 1 0 0 0 0 1 0 0; 0 0 1 0 0 0 0 1 0 0; 0 0 1 1 1 1 0 1 0 0; 0 0 0 0 0 0 0 0 0 0;] 我试图找到不连续连接的行和列(即,行和列的外部范围之间的行和列不均等于1),然后更新它们,以使它们都连接在一起。即上面的矩阵变成: Geometry = [0 0 0 0 0 0 0 0 0 0; 0 0 1 1 1 1 1 1 0 0; 0 0 1 0 0 0 0 1 0 0; 0 0 1 0 0 0 0 1 0 0; 0 0 1 0 0 0 0 1 0 0; 0 0 1 0 0 0 0 1 0 0; 0 0 1 0 0 0 0 1 0 0; 0 0 1 0 0 0 0 1 0 0; 0 0 1 1 1 1 1 1 0 0; 0 0 0 0 0 0 0 0 0 0;] 我遇到的问题是,我希望能够找到每行(和每一列)中等于1的第一个元素和最后一个元素的索引,然后将其用于更新geomtry矩阵。 理想情况下,我想用向量表示它们,因此遍历各列,找到第一个元素的行号等于1并将其存储在称为rowfirst的向量中。 即: rowfirst = zeros(1,numcols) for i = 1:numcols % Going across the columns rowfirst(i) = find(Geometry(i,1) == 1, 1,'first') % Store values in vector called rowfirst end 然后对各列重复此操作,并找到每一行的最后一个元素。 由于某种原因,我无法将值正确存储在向量中,有人知道我要去哪儿吗? 提前致谢。如果不清楚,请告诉我,因为我可能没有很好地解释问题。 回答: 好的,我可以发现两个错误:
for i = 1:numcols tmp = find(Geometry(:, i), 1, 'first'); if(tmp) rowfirst(i) = tmp; else rowfirst(i) = -1; end; end; 更多&回答... |
![]() |
![]() |