基于MATLAB的卫星绕地球模拟系统

夜神月

以下内容摘自好友lyqmath新浪博客，在此借花献佛。

简介

线性代数相关领域中，矩阵操作一般是核心内容。对矩阵的特征值和特征向量的计算一直在数值计算占有重要位置。这里介绍一种根据指定个数来返回矩阵特征值及特征向量的方式。

代码

1. % By lyqmath
   
2. % Dalian University of Technology
   
3. % School of Mathematical Sciences
   
4. % http://blog.sina.com.cn/lyqmath
   
5. clc; clear all; close all;
   
6. % 地球的半径设为100
   
7. R0 = 100;
   
8. a = 12*R0; b = 9*R0;
   
9. % 轨道周期设为T0
   
10. T0 = 2*pi; T = 5*T0;
    
11. dt = pi/100;
    
12. t = [0:dt:T]';
    
13. % 地球与另一焦点的距离
    
14. f = sqrt(a^2-b^2);
    
15. % 卫星轨道与xoy面的倾角
    
16. th = 12.5*pi/180;
    
17. E = exp(-t/20);
    
18. x = E.*(a*cos(t)-f);
    
19. y = E.*(b*cos(th)*sin(t));
    
20. z = E.*(b*sin(th)*sin(t));
    
21. figure; box on; hold on;
    
22. axis off
    
23. plot3(x, y, z, 'k');
    
24. [X, Y, Z] = sphere(30);
    
25. X = R0*X; Y = R0*Y; Z = Z*R0;
    
26. surf(X, Y, Z);
    
27. colormap('hot');
    
28. shading interp;
    
29. % 确定坐标范围
    
30. axis([-18 6 -12 12 -6 6]*R0)
    
31. view([117 37]);
    
32. h = plot3(x(1), y(1), z(1), 'ro', 'MarkerFaceColor', 'g');
    
33. set(gcf, 'Color', 'w');
    
34. title('http://blog.sina.com.cn/lyqmath', 'FontWeight', 'Bold', 'Color', 'r');
    
35. for i = 1 : length(x)
    
36.     set(h, 'xdata', x(i), 'ydata', y(i), 'zdata', z(i));
    
37.     pause(0.01);
    
38. end

复制代码

结果

结果

借助MATLAB来做动态系统的模拟，并利用显示功能来表达出动态效果。可以借助于此方法来应用于不同的研究领域与生活应用。