我正在尝试编写一个模拟,将光源移动到物体表面上,每一步都会测量反射光强度回到光源。
然后绘制表面每个点的强度。随着每一步添加到第一个图。从而建立了表面地形图。
如果你能提出一种整理的方法,那就太麻烦了。
x = -10:0.25:10;
y = -10:0.25:10;
xlength=length(x);
ylength=length(y);
ymin=min(y);
ymax=max(y);
xmin=min(x);
xmax=max(x);
D=zeros(xlength);
[X,Y] = meshgrid(x,y);
z = 10-(X.^2)-(Y.^2)+5*sin(X);
rnd=rand(xlength);
c=randi(xlength);
d=randi(xlength);
s=1;
t=0;
for i=1:ylength
t=t+1;
for j=1:xlength
if s==c
if t==d
D(s,t)=z(s,t);
elseif t==d+1
D(s,t)=z(s,t);
elseif t==d+2
D(s,t)=z(s,t);
elseif t==d+3
D(s,t)=z(s,t);
elseif t==d+4
D(s,t)=z(s,t);
elseif t==d+5
D(s,t)=z(s,t);
else
D(s,t)=0;
end
elseif s==c+1
if t==d
D(s,t)=z(s,t);
elseif t==d+1
D(s,t)=z(s,t);
elseif t==d+2
D(s,t)=z(s,t);
elseif t==d+3
D(s,t)=z(s,t);
elseif t==d+4
D(s,t)=z(s,t);
elseif t==d+5
D(s,t)=z(s,t);
else
D(s,t)=0;
end
elseif s==c+2
if t==d
D(s,t)=z(s,t);
elseif t==d+1
D(s,t)=z(s,t);
elseif t==d+2
D(s,t)=z(s,t);
elseif t==d+3
D(s,t)=z(s,t);
elseif t==d+4
D(s,t)=z(s,t);
elseif t==d+5
D(s,t)=z(s,t);
else
D(s,t)=0;
end
elseif s==c+3
if t==d
D(s,t)=z(s,t);
elseif t==d+1
D(s,t)=z(s,t);
elseif t==d+2
D(s,t)=z(s,t);
elseif t==d+3
D(s,t)=z(s,t);
elseif t==d+4
D(s,t)=z(s,t);
elseif t==d+5
D(s,t)=z(s,t);
else
D(s,t)=0;
end
elseif s==(c+4)
if t==d
D(s,t)=z(s,t);
elseif t==d+1
D(s,t)=z(s,t);
elseif t==d+2
D(s,t)=z(s,t);
elseif t==d+3
D(s,t)=z(s,t);
elseif t==d+4
D(s,t)=z(s,t);
elseif t==d+5
D(s,t)=z(s,t);
else
D(s,t)=0;
end
else
D(s,t)=0;
end
s=s+1;
end
s=1;
end
z1=z -(D/20);
z2=z1-z;
s=0;
figure
surf(X,Y,z)
axis([xmin xmax ymin ymax])
xlabel('X')
ylabel('Y')
zlabel('Z')
for i=-10:2.5:10
hold on
light('position',[i,0,50])
surf(X,Y,z,'EdgeColor', 'none')
axis([xmin xmax ymin ymax])
drawnow
pause (1)
delete(findall(gcf,'Type','light'))
hold off
end
这是我所拥有的。
答案 0 :(得分:0)
完全符合您要求的代码:(感谢@Hoki帮助您更简化代码)
clear;clc
x = -10:0.25:10;
y = -10:0.25:10;
xlength=length(x);
ylength=length(y);
ymin=min(y);
ymax=max(y);
xmin=min(x);
xmax=max(x);
[X,Y] = meshgrid(x,y);
z = 10-(X.^2)-(Y.^2)+5*sin(X);
%% plot!
figure
surf(X,Y,z,'EdgeColor', 'none','lineStyle','none') % taken out the linestyle because it looks cooler
axis([xmin xmax ymin ymax])
xlabel('X')
ylabel('Y')
zlabel('Z')
% Added repetitions and decreased step size, plus making it go forward and backward.
repetitions=4;
for jj=1:repetitions
hl = light('position',[-10,0,50]) ; %// create the light source
for i=-10:2.5:10
set(hl,'position',[i,0,50]) ; %// set the new light position
drawnow %// flush graphic pipeline
pause (0.1) %// let human user see something
end
delete(hl) %// delete the light source
end
这给出了以下gif:
http://i.imgur.com/42ffYll.gifv
关于代码的其余部分......
请记住
if t==d
D(s,t)=z(s,t);
elseif t==d+1
D(s,t)=z(s,t);
elseif t==d+2
D(s,t)=z(s,t);
elseif t==d+3
D(s,t)=z(s,t);
elseif t==d+4
D(s,t)=z(s,t);
elseif t==d+5
D(s,t)=z(s,t);
else
D(s,t)=0;
end
与字面意思相同:
if (t>=d && t<=d+5)
D(s,t)=z(s,t);
else
D(s,t)=0;
end
这可以帮助您减少代码。 还有矢量化。
你可以做到
D(s,(t>=d && t<=d+5))=z(s,(t>=d && t<=d+5)));
对于s
等所有值等等。希望这能让你前进;)