专业今年开了MATLAB的课程(土木工程学MATLAB?),老师给出了保过题,为了不挂科也是苦逼百度+谷歌了半天,加上自己的琢磨,总算是明白些surf颜色控制方法了,在这里分享记录一下~初学matlab,理解的不是太深,望海涵。
题目
从excel中读取粒子,这些粒子具有三维的坐标,xyz,并且具有三维的速度,excel中的每一页内存储一个时间点的例子信息,x,y,z,vx,vy,vz,页的命名为sheet+i,希望能够先读取粒子信息到向量当中,然后粒子的位置进行图像显示,粒子采用球体进行显示,球体的颜色和合速度的大小有关,采用从绿色到红色的渐变关系,然后粒子在三维空间中,根据时间的变化运动起来。
说明:我希望时间步数是作为函数的输入参数的,也就是显示多少个时间点是可控的,起始的时间点是可控的
从我多年编程的角度,嗯,逻辑非常简单,但是动画这块就比较蛋疼了,毕竟完全没接触过matlab动画~于是开始在网上搜相关的教程,这些东西还是挺多的。然后就做成了生成一个球并使其运动的代码,然后hold on绘制多个球。
然后就是颜色控制,第一想法就是使用colormap,用了之后就发现后面设置的会覆盖之前的,就是所有球都是同一颜色,,,。
之后继续搜,发现surf(x,y,z,c)的第四个参数就是颜色控制,,,但是搜了好多都说是颜色控制,就是不说到底怎么控制,,这就比较尴尬了。
于是在调试中和搜索中,发现错误提示 cDATA必须是M*N矩阵或者M*N*3数组,而c的默认值为z,,,那么不得不去联想,surf函数是使用网格绘制,c对应每一个网格的颜色,那也就是每个网格可以使用rgb进行颜色设置。那么设置一个和z同大小的矩阵,然后升级为三维数组,第三维(列)对应rgb三色,传给surf函数,那么正好实现了颜色控制。于是测试(以下为代码片段,无上下文)
[x,y,z] = sphere(n);
c=zeros(size(x));%获得o阵大小和x相同
for i=1:1:length(c(1,:))
for j=1:1:length(c(:,1))
c(i,j,1)=1;
c(i,j,2)=0;
c(i,j,3)=0;%红色
end
end
没问题,所有球都是红的
那么题目就迎刃而解了
function anim=ani2(numb,start)
dian=[];
v=[];
i=1;
vmax=0;
vmin=100;
while 1
lsp=xlsread(‘dian.xlsx‘,i);
if length(lsp)==0
break;
end
lsp(1)=lsp(1)+lsp(4)*start;
lsp(2)=lsp(2)+lsp(5)*start;
lsp(3)=lsp(3)+lsp(6)*start;
dian=[dian;lsp(1),lsp(2),lsp(3)];
v=[v;lsp(4),lsp(5),lsp(6)];
lsv=sqrt(lsp(4)*lsp(4)+lsp(5)*lsp(5)+lsp(6)*lsp(6));
if lsv>vmax
vmax=lsv;
end
if lsv<vmin
vmin=lsv;
end
i=i+1;
end
k = 5;
n = 2^k-1;
[x,y,z] = sphere(n);
for i=1:1:length(dian(:,1))
surf(x+dian(i,1),y+dian(i,2),z+dian(i,3));
hold on;
end
axis equal
axis tight
set(gcf,‘renderer‘,‘zbuffer‘);
set(gca,‘nextplot‘,‘replacechildren‘);
title(‘球体动画‘);
j_color=vmax-vmin;
for j=1:1:numb
k = 5;
n = 2^k-1;
[x,y,z] = sphere(n);
set(gca,‘nextplot‘,‘replacechildren‘);
for i=1:1:length(dian(:,1))
dian(i,1)=dian(i,1)+v(i,1);
dian(i,2)=dian(i,2)+v(i,2);
dian(i,3)=dian(i,3)+v(i,3);
g=1-(sqrt(v(i,1)*v(i,1)+v(i,2)*v(i,2)+v(i,3)*v(i,3))-vmin)/j_color;
r=1-g;
c=zeros(size(z));
for ll=1:1:length(c(1,:))
for j=1:1:length(c(:,1))
c(ll,j,1)=r;
c(ll,j,2)=g;
c(ll,j,3)=0;
end
end
surf(x+dian(i,1),y+dian(i,2),z+dian(i,3),c);
hold on;
end
axis equal
F(j)=getframe;
end
希望对大家使用surf时有所帮助~