[绘图] 想请教有关画四维图方法(更新问题)

楼主: zx168zx (zx168zx)   2015-04-02 15:11:22
小弟我想画一个f(x,y,z)的图
x,y,z为3个矩阵
大小分别为
x => 64,1
y => 64,1
z => 57,1
f => 64x64x57,1
f的值就是 当x=1,y=1,z=1 带入以后得出的值
我想要在三维座标上 将f的大小用颜色(colorbar)表示
我试过scatter3(x(:,1),y(:,1),z(:,1),f,f,'filled')
本来想说用surf画 但meshgrid的矩阵也是要同样大小
但他显示我的x,y,z不一样大 而且scatter3画出来是点
我想问看看 是不是有什么方法可以画出f的曲面 然后颜色表示大小?
以下是P大提供的方法 给大家参考
[xx yy]=meshgrid( 64 x 64 ) %将资料做成meshgrid
for n = 1:57
f =f2(:,:,n); %用循环一页一页画
mesh(xx,yy,f);
hold on %叠图
colorbar
drawnow
frame = getframe(1);
im = frame2im(frame);
[imind,cm] = rgb2ind(im,256);
if n == 1;
imwrite(imind,cm,filename,'gif', 'Loopcount',inf);
else
imwrite(imind,cm,filename,'gif','WriteMode','append');
end
end
将z的变化改用叠图的方式 处理出来效果蛮不错的
另外我想问的是 如果我现在x,y,z很大的话 要如何处理
比如 f(x,y,z) 里面的 x y z各是10000的话
这样meshgrid会因为内存问题无法建立
请问有什么方法可以处理这么大的资料量吗?
作者: profyang (prof)   2015-04-02 20:55:00
感觉你需要的是3d的contour(一般contour是z=f(x,y), 你需要的是t=f(x,y,z), 但是这样画真的会清楚吗? 感觉例如f=1的曲面跟f=2的曲面就有可能互相挡住之类的如果是我宁可做成gif动画 让x,y,z其中一自变量随时间变仍然画成f(x,y)然后z随时间变的2d contour
作者: sin55688 (单手挑藏獒)   2015-04-02 21:33:00
依你叙述,结果会是充满颜色的方块耶?这是你要的吗?
作者: profyang (prof)   2015-04-04 15:27:00
你只是要叠图就不用后面写成gif的东西了阿XD 我本意是怕叠出来曲面之间可能互相挡到才建议你用动画 既然你觉得效果不错代表不会挡到 那就for+hold on去叠就好了 没必要写gif另外资料量太大感觉很难...可能只能用for把x,y也分段画去叠对阿 把x y分小一点

Links booklink

Contact Us: admin [ a t ] ucptt.com