remove inf in quiver
이전 댓글 표시
I want quiver(X,Y,u,v), but there's inf entries in u and v at positions x=y. I'm looking for the smartest way to skip these positions with inf u and v and finish the quiver.
x=-5:0.1:5;
y=-5:0.1:5;
[X,Y]=meshgrid(x,y);
댓글 수: 1
Dyuman Joshi
2024년 4월 28일
u and v aren't defined in the code.
채택된 답변
Without having ‘u’ and ‘v’ to work with, perhaps something like this using fillmissing (or fillmissing2) —
x=-5:0.1:5;
y=-5:0.1:5;
[X,Y]=meshgrid(x,y);
u = randn(size(X)); % Create 'u'
ix = sub2ind(size(u), 1:size(u,1), 1:size(u,2)); % Linear Inmdex To Create 'u' With Diagnonal 'Inf'
u(ix) = Inf
u(isinf(u)) = NaN % Change 'Inf' To 'NaN'
u = fillmissing(u,'linear') % Use 'fillmissing'
Use whatever interpolation method you want with fillmissing. There are several options.
.
댓글 수: 4
feynman feynman
2024년 4월 28일
Thank you!
The problem with leaving them blank is that reduces the matrix by 1 in the row size. That makes plotting them with the original (X,Y) matrices inpossible, unless you also delete the corresponding diagonal elements in the (X,Y) matrices.
One way to do that is to use a version of my original code for all the matrices —
x=-5:0.1:5;
y=-5:0.1:5;
[X,Y]=meshgrid(x,y)
u = randn(size(X)) % Create 'u'
ix = sub2ind(size(u), 1:size(u,1), 1:size(u,2)); % Linear Inmdex To Create 'u' With Diagnonal 'Inf'
X(ix) = [];
X = reshape(X,100,[])
Y(ix) = [];
Y = reshape(Y,100,[])
u(ix) = [];
u = reshape(u, 100,[])
Then do the same sort of operation with ‘v’.
Another (perhaps preferable) option is to leave them as NaN values. The NaN values will not plot, and any calculations involving them will also be NaN, however if there are any recursive operations involving the matrices, that could leave many more elements a NaN values.
As I mentioned earlier, fillmissing has other options to fill the NaN values if you want to use them, for example a constant value. Interpolating them using linear or other methods is not absolutely necessary.
.
feynman feynman
2024년 4월 29일
Star Strider
2024년 4월 29일
As always, my pleasure!
추가 답변 (0개)
카테고리
도움말 센터 및 File Exchange에서 Creating and Concatenating Matrices에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
