matlab代码
close all;
clear;
oR=[cos(pi),-sin(pi),0;
sin(pi),cos(pi),0;
0,0,1];
stp=1;
nsp=20;
ox=-11;
oy=0;
oz=0;
%[u,v,w]=sphere;
x=0;
%[u,v,w] = sphere(20);
for i=0:1:nsp
%生成一个地球
R=3; %地球半径
[k,j,l] = sphere(200); %获得三维球面得坐标,生成三个n乘n得矩阵,n默认为20,三个矩阵分别是x,y,z单位球得坐标,然后是一列一列来的
k=reshape(k,1,[]); %B = reshape(A,m,n) 将矩阵A的元素返回到一个m×n的矩阵B。[]代表让他自己生成
j=reshape(j,1,[]); %这里相当于把矩阵展平,变成一向量
l=reshape(l,1,[]);
Tz=[cos(x),-sin(x),0; %地球自转速度
sin(x),cos(x),0;
0,0,1];
p=Tz*[k;j;l];
% k=k'*Tz(1,:);
% j=j'*Tz(2,:);
% l=i'*Tz(3,:);
plot3(p(1,:),p(2,:),p(3,:)); %分别取第一列第二列三列,也就是x