www.gusucode.com > gps卫星位置预报，MATLAB编写，有动画界面源码程序 > code5/gps_final_new/DrawEarth.m

    function DrawEarth(time)
%绘制地球
%time是参数
%利用这个参数，可以绘制一个看起来是旋转的地球
r=6400;              %地球半径
j1=[0:pi/12:2*pi];
w1=[-pi/2:pi/12:pi/2];
L1=length(w1);
L2=length(j1);
for n=1:L1
    z=ones(L2,1);
    z=z*r*sin(w1(n));
    temp=r*cos(w1(n));
    x=temp*sin(j1);
    y=temp*cos(j1);
    plot3(x,y,z);
    hold on;
    grid;
end

unit=ones(1,1);
z0=ones(1,1);
x0=ones(1,1);
y0=ones(1,1);


for n=1:L2
    for m=1:L1
        temp=w1(m);
        temp2=j1(n)+time*pi/12;
        z=r*sin(temp);
        x=r*cos(temp)*sin(temp2);
        y=r*cos(temp)*cos(temp2);
        z1=unit*z;
        x1=unit*x;
        y1=unit*y;
            
        z0=[z0 z1];
        x0=[x0 x1];
        y0=[y0 y1];
    end
    
    
    z0(:,1)=[];
    x0(:,1)=[];
    y0(:,1)=[];
    plot3(x0,y0,z0);
    axis equal;
    %axis off;
    hold on;
end