pcolor

2011 6 7
2 个回答
回答已采纳
5 次查看（30 天）

0 个投票

Hi dear MATLAB user,
Im making a FDTD (acoustic preassure) model and at the end of each calculation representing the preassure. The problem is by wieving the contours, it only represents the actual preassure but if the wave hasnt arrived yet you wont see nothing from the contours.
for ...
for ...
end
...
pcolor((1:nx)*dh,(1:ny)*dh,10*log10(min(max(pp'.^2,1d-6),1)) );
shading flat,axis equal; colorbar; title(['time=' num2str(round((tt+maxtt*(ttt-1))*1000*dt)) ' ms']);drawnow
hold on
[x,y]=find(mapaz == 0);
contorno=plot(x,y,'+k');
set(contorno,'MarkerSize',10);
hold off
drawnow
end
Im trying to draw the contours stored in matrix 'mapaz': 0 = contour; 1 nothing.
nx and ny are the sizes of the matrix (1600x1200)
I really don't know to mix this together. Could you help me?
Thanks, Daniel Serra

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

 采纳的回答

Walter Roberson
Walter Roberson 2011-6-7

1 个投票

Perhaps
contour((1:nx)*dh, (1:ny)*dh, 0 + (mapaz ~= 0), [1 1])

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Walter Roberson
Walter Roberson 2011-6-7
I just realized that I should probably have used == 0 instead of ~= 0

请先登录，再进行评论。

更多回答（1 个）

Daniel
Daniel 2011-6-7

0 个投票

Thankyou very much! It worked fine :)
I had to transpose 0 + (mapaz ~...) and it worked fine.
Here is the final code if someone needs it.
Where:
nx,ny = 1, 2, 3...
dh = 0.01 % 1 cm
mapaz = matrix with 0 and 1
CODE:
contour((1:nx)*dh, (1:ny)*dh, (0 + (mapaz ~= 0))', [1 1])
hold on;
pcolor((1:nx)*dh,(1:ny)*dh,10*log10(min(max(pp'.^2,1d-6),1)) );
shading flat,axis equal; colorbar; title(['time=' num2str(round((tt+maxtt*(ttt-1))*1000*dt)) ' ms']);drawnow
Thanks!

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Walter Roberson
Walter Roberson 2011-6-7
You might not have to add the 0, but I put it in for safety.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Contour Plots 的更多信息

产品

标签

Community Treasure Hunt

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

Start Hunting!

Translated by