Index of points which satisfy a given equation

4 次查看（过去 30 天）
Ali Baig
Ali Baig 2019-8-3
评论： Ali Baig 2019-8-3
Hellow everyone,
I am working on a problem in which I have defined following variables:
NodesAlongOneDirection = 10
x1 = 0; x2 = 1; Nx = NodesAlongOneDirection; x = linspace(x1, x2, Nx); dx = x(2) - x(1);
y1 = 0; y2 = 1; Ny = NodesAlongOneDirection; y = linspace(y1, y2, Ny); dy = y(2) - y(1);
[X, Y] = meshgrid(x, y); X = X(:); Y = Y(:); N = length(X);
Now I want to find index of points (index of X, Y) which satisfy any given equation, say, Y = -X +1, which is equation of "secondary diagonal" of square domain.
I have tried following lines:
A = find(Y == -X +1);
But it doesn't contain all the points on secondary diagonal.
Regards,
Ali Baig
  2 个评论
Walter Roberson
Walter Roberson 2019-8-3
You should not be using floating point equality as it demands bit for bit exact comparisons. You should compare with a tolerance.
Ali Baig
Ali Baig 2019-8-3
Dear Walter Roberson
Thank you for reply. Keeping what you have suggested, I have written the following lines:
Y_Low = -X + 1 - 0.25 * sqrt(dx^2 + dy^2);
Y_High = -X + 1 + 0.25 * sqrt(dx^2 + dy^2);
A = find(Y > Y_Low & Y < Y_High);
figure()
plot(X(A), Y(A), '*')
This works fine at least for points on secondary diagonal. The "0.25 * sqrt(dx^2 + dy^2)" is used as a tolerance value.
Once again, thank you so much for your help.
Regards,
Ali Baig

请先登录，再进行评论。

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

产品


版本

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by