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

基于matlab GPS卫星仿真

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
基于matlab GPS卫星仿真5 i. _& U  {% }& N
- c, h: o2 v2 N. L/ w$ o4 J! H8 e  l
一、部分源代码/ T6 B% Q8 t/ w, u3 Q+ q, ^  h/ q
  • 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
    # G0 X) d3 g2 c- u# O* I, P

& E! y% r" e4 S* ~7 u/ N
+ O1 I8 X4 y: h, v
5 b& f# U( d3 P4 Q6 `$ S二、运行结果
8 e" R! F7 X( u& h' A+ I
: t4 G" Y' a! w- W% q $ h* T. h" x+ C; O
8 v# M! ?9 X9 J6 Y# |, P

: ?# m, g9 f* h$ X: k6 a
) D1 ~* W& X' M0 r0 o
3 h( K! h$ z* A4 Z, T' T
0 w. \- @8 _' s0 b1 Q
9 C  f+ f! d* A3 J

该用户从未签到

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-28 20:20 , Processed in 0.156250 second(s), 26 queries , Gzip On.

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

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

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