How can I plot the streamlines in all region around the inside sphere ?

4 次查看(过去 30 天)

采纳的回答

Star Strider
Star Strider 2024-1-2
编辑:Star Strider 2024-1-2
One approach is to change the code to add a patch object in this part of it (add ‘pv’ and ‘cv’ and the patch call):
set(polar(th,r1,'-k'),'LineWidth',1.2);set(polar(th,r2,'-k'),'LineWidth',1.2);
pv = linspace(0, 2*pi*1.01, 250); % Used in 'patch'
cv = [cos(pv); sin(pv)]; % Used in 'patch'
patch(cv(1,:), cv(2,:), 'w') % Plot 'patch'
title('$\eta_1=\eta_2=0.0$','Interpreter','latex','FontSize',9,'FontName','Times New Roman','FontWeight','Normal')
First, delete (or comment-out) these lines:
if ~isempty(I);
x(I,J) = 0;
y(I,J) = 0;
end
They create the centre square, and you need to get rid of it.
This plots everything, then creates a white patch object for the inside circle to cover the contour lines in it. (An alternative would be to get all the contour coordinates from ‘DH1’ and then delete the (x,y) pairs that are inside the inner circle, defined here by the patch coordinates ‘pv’ and ‘cv’. That could be done, however it would significantly slow the code. Creating the patch object is just easier.)
(The entire code is too long to post, exceeding the maximum character limit.)
EDIT — (2 Jan 2024 at 22:20)
The result is something like this —
.
  6 个评论

请先登录,再进行评论。

更多回答(1 个)

Shreen El-Sapa
Shreen El-Sapa 2024-1-7
Dear Star,can I plot half circles only?
  6 个评论

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Specifying Target for Graphics Output 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by