简介
利用MATLAB来做系统模拟，以其显示功能、动画功能而突出。这里介绍一个卫星（月球）绕地球运行系统。
（注：部分代码参考网络分享）
实例
% By lyqmath
% Dalian University of Technology
% School of Mathematical Sciences
%
http://blog.csdn.net/lyqmath
clc; clear all; close all;
% 地球的半径设为100
R0 = 100;
a = 12*R0; b = 9*R0;
% 轨道周期设为T0
T0 = 2*pi; T = 5*T0;
dt = pi/100;
t = [0:dt:T]’;
% 地球与另一焦点的距离
f = sqrt(a^2-b^2);
% 卫星轨道与xoy面的倾角
th = 12.5*pi/180;
E = exp(-t/20);
x = E.*(a*cos(t)-f);
y = E.*(b*cos(th)*sin(t));
z = E.*(b*sin(th)*sin(t));
figure; box on; hold on;
axis off
plot3(x, y, z, ‘k’);
[X, Y, Z] = sphere(30);
X = R0*X; Y = R0*Y; Z = Z*R0;
surf(X, Y, Z);
colormap(‘hot’);
shading interp;
% 确定坐标范围
axis([-18 6 -12 12 -6 6]*R0)
view([117 37]);
h = plot3(x(1), y(1), z(1), ‘ro’, ‘MarkerFaceColor’, ‘g’);
set(gcf, ‘Color’, ‘w’);
title(‘http://blog.sina.com.cn/lyqmath’, ‘FontWeight’, ‘Bold’, ‘Color’, ‘r’);
for i = 1 : length(x)
set(h, ‘xdata’, x(i), ‘ydata’, y(i), ‘zdata’, z(i));
pause(0.01);
end
结果
总结
借助MATLAB来做动态系统的模拟，并利用显示功能来表达出动态效果。可以借助于此方法来应用于不同的研究领域与生活应用。