作为我目前正在进行的项目的一部分,我必须解决国际空间站绕地球轨道运行的两体问题。到目前为止,我已经设法通过使用sphere / surf函数来近似此值,但是,我想知道是否有任何方法可以创建一个更真实的代表ISS的图形?不幸的是,这个项目只能通过MATLAB来完成,因此我不能使用任何其他工具来提供更好的可视化效果
答案 0 :(得分:2)
NASA具有许多对象的3D模型,包括ISS which can be found here。可以根据需要将这个文件转换为STL,我发现this random website对我有用。
在Matlab中,您可以通过以下方式读取该文件
stl = stlread('isscombined.stl');
V = stl.Points;
F = stl.ConnectivityList
然后,您可以使用
进行绘制p = patch('vertices',V,'faces',F,'FaceColor',[.8 .8 .8]);
,然后您可以在测站绕地球旋转时用新的顶点位置更新对象。显然,您也可以通过将顶点乘以一定数量来缩放对象。如果您不想绘制刻面边缘,也可以将'EdgeAlpha', 0
添加到patch
选项中。
这是一个简单的示例,显示了ISS绕球飞行
% Note: not to scale
ISS_radius = 2; % distance from center of Earth
RE = 1; % radius of earth
theta = 0:.05:2*pi;
x = ISS_radius*cos(theta);
y = ISS_radius*sin(theta);
stl = stlread('isscombined.stl');
r = .01; % scaling factor
V = stl.Points * r;
V = V - mean(V); % center at origin
F = stl.ConnectivityList;
figure; hold on;
plot3(x,y,zeros(numel(theta)),'--');
[X,Y,Z] = sphere(50);
surf(RE*X,RE*Y,RE*Z,'FaceColor',[0 0 .8],'EdgeAlpha',0);
p = patch('Vertices', V*r, 'Faces', F, 'FaceColor', [0 0 0], 'EdgeAlpha', 0);
axis equal;
set(gca,'View',[200 13])
grid on;
counter = 1;
while true
p.Vertices = V + [x(counter), y(counter), 0];
pause(0.01);
drawnow
counter = mod(counter + 1, numel(theta)) + 1;
axis([-1 1 -1 1 -1 1]*ISS_radius*1.2)
end