找回密码
 注册
关于网站域名变更的通知
查看: 686|回复: 3
打印 上一主题 下一主题

基于matlab GPS卫星仿真

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-9-29 14:34 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
基于matlab GPS卫星仿真7 z% n; i  [, z' r# a# v

5 N0 ^/ r& E4 Q  T6 `一、部分源代码: x% w8 `6 W: |4 A) t4 y7 C, G
  • 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! r& [2 _# I) O

$ d  o, p/ R: d: C" Q- `! j2 @
" g) T: }6 f1 N0 m2 |. ?. E. o) m/ f9 y0 o. l
二、运行结果% M! K* ]5 u1 v& h. l+ V
- {( O7 W2 n  g! T2 {4 c  _5 _
! w" U* m: E8 a; }/ @$ F
7 N9 Y5 q- }6 Q2 |5 Z( ^+ w( J9 q
' ~- Q& R& G! s: W/ w
: ~0 v0 U' H. d7 f! w! B
2 H% o$ u1 Y' b

: V# q- s# i7 t+ V' \& ?. T1 h+ z+ W7 I1 x  |, e

该用户从未签到

2#
发表于 2021-9-29 15:28 | 只看该作者
基于matlab GPS卫星仿真

该用户从未签到

3#
发表于 2021-9-29 16:22 | 只看该作者
基于matlab GPS卫星仿真
  • TA的每日心情

    2019-11-19 15:29
  • 签到天数: 1 天

    [LV.1]初来乍到

    4#
    发表于 2021-9-29 16:23 | 只看该作者
    基于matlab GPS卫星仿真
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-10-5 01:07 , Processed in 0.125000 second(s), 27 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表