【气动学】基于matlab GPS卫星仿真【含Matlab源码 1094期】

一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【气动学】基于matlab GPS卫星仿真【含Matlab源码 1094期】

获取代码方式2：
通过订阅紫极神光博客付费专栏，凭支付凭证，私信博主，可获得此代码。

备注：
订阅紫极神光博客付费专栏，可免费获得1份代码（有效期为订阅日起，三天内有效）；

二、部分源代码

clear;clc;close all;
   a=26560;
    e=0.02;
    E=[0:0.1:2*pi];
    x=a*(cos(E)-e);
    y=a*sqrt((1-e^2))*sin(E);
    z=0*E;
    TimeNow=0;
    global SatellitePosition
    SatellitePosition=ones(24,4);
    figure(1);
    DrawEarth(0);
    hold on;
    
DtoR=2*pi/360;
 A1=[32.8 92.8 152.8 212.8 272.8 332.8 ];
 DrawSatelliteOrbit;
 Ctable=[10 50 160 260 ;
    80 180 220 320 ;
    10 130 250 340;
    50 150 170 300;
    100 210 310 340 ;
    140 150 240 350;];
   Wx=ones(1,1);Wy=ones(1,1);Wz=ones(1,1);
   simple=1;
   for k=1:6
       A=A1(k)*DtoR; 
       B=55*DtoR;
     for m=1:4
           C=Ctable(k,m)*DtoR+TimeNow*2*pi/24;
           x=a*(cos(C)-e);
           y=a*sqrt((1-e^2))*sin( C);
           z=0*C;
           R3=[cos(A)  -sin(A)  0;
                sin(A)  cos(A) 0;
                 0        0     1;];
    R1=[1         0    0;
        0       cos(B)  -sin(B);
        0       sin(B) cos(B);];
    R2=[cos(C) -sin(C) 0;
        sin(C) cos(C) 0; 
        0          0  1;];
    L1=length(E);
    R312=R3*R1*R2;
    Ans=R312*[x;y;z;];
    Wx=[Wx Ans(1,:)];
    Wy=[Wy Ans(2,:)];
    Wz=[Wz Ans(3,:)];
    x1=Ans(1,:);
    y1=Ans(2,:);
    z1=Ans(3,:);
    DrawSatellite(x1,y1,z1,k);
    temp=k*m;
    SatellitePosition(temp,:)=[x1 y1 z1 1]
   hold on;
   
end
end

EarthCenterPos=[0 0 0];
UserPosition=[6400 0 0];
for k=1:24
Temp=SatellitePosition(k,1:3)-UserPosition;
Dist1=Temp*Temp';
Temp=UserPosition-EarthCenterPos;
Dist2=Temp*Temp';
Temp=SatellitePosition(k,1:3)-EarthCenterPos;
Dist3=Temp*Temp';
JiaoJiao=acos((Dist1+Dist2-Dist3)/2/sqrt(Dist1)/sqrt(Dist2));
if(JiaoJiao<=pi/2)
    SatellitePosition(k,4)=0;
end
end

  
 

  
 
 

  
1
  
2
  
3
  
4
  
5
  
6
  
7
  
8
  
9
  
10
  
11
  
12
  
13
  
14
  
15
  
16
  
17
  
18
  
19
  
20
  
21
  
22
  
23
  
24
  
25
  
26
  
27
  
28
  
29
  
30
  
31
  
32
  
33
  
34
  
35
  
36
  
37
  
38
  
39
  
40
  
41
  
42
  
43
  
44
  
45
  
46
  
47
  
48
  
49
  
50
  
51
  
52
  
53
  
54
  
55
  
56
  
57
  
58
  
59
  
60
  
61
  
62
  
63
  
64
  
65
  
66
  
67
  
68
  
69
  
70
  
71
  
72
  
73
 

三、运行结果

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 门云阁.MATLAB物理计算与可视化[M].清华大学出版社，2013.

文章来源: qq912100926.blog.csdn.net，作者：海神之光，版权归原作者所有，如需转载，请联系作者。

原文链接：qq912100926.blog.csdn.net/article/details/118463909