GPS卫星运动及定位matlab仿真
1.2 本課題研究的意義和方法
GPS系統(tǒng)是一個(gè)很龐大的系統(tǒng)，包含了天文，地理，計(jì)算機(jī)，電磁學(xué)，通信學(xué)，信息學(xué)等等。通過(guò)本文對(duì)GPS的學(xué)習(xí)研究，最重要的還是要學(xué)習(xí)其原理：衛(wèi)星運(yùn)動(dòng)原理；衛(wèi)星定位原理；衛(wèi)星跟蹤原理等等。通過(guò)基礎(chǔ)原理的學(xué)習(xí)，一方面，可以使我們更進(jìn)一步的理解衛(wèi)星運(yùn)動(dòng)，定位的實(shí)現(xiàn)方法；通過(guò)仿真，進(jìn)一步了解簡(jiǎn)單定位的方法及其在仿真平臺(tái)上的實(shí)現(xiàn)途徑；另一方面，也可以培養(yǎng)我們自學(xué)的能力，訓(xùn)練仿真模擬的技巧和方法。
至今，基本上完成了課題的要求，通過(guò)不斷的注入既定參數(shù)，可以更加詳細(xì)，直觀的理解基本的定位原理和實(shí)現(xiàn)方法！
1.3 GPS前景
GPS導(dǎo)航定位以其定位精度高、觀測(cè)時(shí)間短、測(cè)站間無(wú)需通視、可提供三維坐標(biāo)、操作簡(jiǎn)便、全天候作業(yè)、功能多、應(yīng)用廣泛等特點(diǎn)著稱。
用GPS信號(hào)可以進(jìn)行海、空和陸地的導(dǎo)航、導(dǎo)彈的制導(dǎo)、大地測(cè)量和工程測(cè)量的精密定位、時(shí)間的傳遞和速度的測(cè)量等。對(duì)于測(cè)繪領(lǐng)域，GPS衛(wèi)星定位技術(shù)已經(jīng)用于建立高精度的全國(guó)性的大地測(cè)量控制網(wǎng)，測(cè)定全球性的地球動(dòng)態(tài)參數(shù)；用于建立陸地海洋大地測(cè)量基準(zhǔn)，進(jìn)行高精度的海島陸地聯(lián)測(cè)以及海洋測(cè)繪；用于檢測(cè)地球板塊運(yùn)動(dòng)狀態(tài)和地殼形變；用于工程測(cè)量，成為建立城市與工程控制網(wǎng)的主要手段。用于測(cè)定航空航天攝影瞬間相機(jī)位置，實(shí)現(xiàn)僅有少量的地面控制或無(wú)地面控制的航測(cè)快速成圖，導(dǎo)致地理信息系統(tǒng)、全球環(huán)境遙感監(jiān)測(cè)的技術(shù)革命[4]。
目前，GPS、GLONASS、INMARSAT等系統(tǒng)都具備了導(dǎo)航定位功能，形成了多元化的空間資源環(huán)境。這一多元化的空間資源環(huán)境，促使國(guó)際民間形成了一個(gè)共同的策略，即一方面對(duì)現(xiàn)有系統(tǒng)充分利用，一方面積極籌建民間GNSS系統(tǒng)，待2011年左右，GNSS純民間系統(tǒng)建成，全球?將形成GPS/GLONASS/GNSS三足鼎立之勢(shì)，才能從根本上擺脫對(duì)單一系統(tǒng)的依賴，形成國(guó)際共有、國(guó)際共享的安全資源環(huán)境。世界才可以將衛(wèi)星導(dǎo)航作為單一導(dǎo)航手段的最高應(yīng)用境界。國(guó)際民間的這一策略，反過(guò)來(lái)又影響和迫使美國(guó)對(duì)其GPS使用政策作出更開放的調(diào)整。多元化的空間資源環(huán)境的確立，給GPS的發(fā)發(fā)展應(yīng)用創(chuàng)造了一個(gè)前所未有的良好的國(guó)際環(huán)境。
?
?
第二章 GPS測(cè)量原理
GPS導(dǎo)航系統(tǒng)的基本原理是測(cè)量出已知位置的衛(wèi)星到用戶接收機(jī)之間的距離 ,然后綜合多顆衛(wèi)星的數(shù)據(jù)就可知道接收機(jī)的具體位置。要達(dá)到這一目的 ,衛(wèi)星的位置可以根據(jù)星載時(shí)鐘所記錄的時(shí)間在衛(wèi)星星歷中查出。而用戶到衛(wèi)星的距離則通過(guò)紀(jì)錄衛(wèi)星信號(hào)傳播到用戶所經(jīng)歷的時(shí)間 ,再將其乘以光速得到(由于大氣層電離層的干擾 ,這一距離并不是用戶與衛(wèi)星之間的真實(shí)距離 ,而是偽距( PR) :當(dāng)GPS衛(wèi)星正常工作時(shí) ,會(huì)不斷地用 1 和 0 二進(jìn)制碼元組成的偽隨機(jī)碼(簡(jiǎn)稱偽碼)發(fā)射導(dǎo)航電文。GPS系統(tǒng)使用的偽碼一共有兩種 ,分別是民用的 C/ A 碼和軍用的 P( Y)碼。C/ A 碼頻率 1. 023MHz ,重復(fù)周期一毫秒 ,碼間距 1 微秒 ,相當(dāng)于 300m; P 碼頻率10. 23MHz ,重復(fù)周期266. 4 天 ,碼間距0. 1 微秒 ,相當(dāng)于 30m。而 Y碼是在 P碼的基礎(chǔ)上形成的 ,保密性能更佳。
GPS導(dǎo)航系統(tǒng)衛(wèi)星部分的作用就是不斷地發(fā)射導(dǎo)航電文。然而 ,由于用戶接受機(jī)使用的時(shí)鐘與衛(wèi)星星載時(shí)鐘不可能總是同步 ,所以除了用戶的三維坐標(biāo) x、 y、 z外 ,還要引進(jìn)一個(gè)Δt 即衛(wèi)星與接收機(jī)之間的時(shí)間差作為未知數(shù) ,然后用 4 個(gè)方程將這 4個(gè)未知數(shù)解出來(lái)。所以如果想知道接收機(jī)所處的位置 ,至少要能接收到 4 個(gè)衛(wèi)星的信號(hào)。
2.1 偽距測(cè)量的原理
GPS定位采用的是被動(dòng)式單程測(cè)距。它的信號(hào)發(fā)射書機(jī)由衛(wèi)星鐘確定，收到時(shí)刻是由接收機(jī)鐘確定，這就在測(cè)定的衛(wèi)星至接收機(jī)的距離中，不可避免地包含著兩臺(tái)鐘不同步的誤差和電離層、對(duì)流層延遲誤差影響，它并不是衛(wèi)星與接受機(jī)之間的實(shí)際距離，所以稱之為偽距。
偽距定位法是利用全球衛(wèi)星定位系統(tǒng)進(jìn)行導(dǎo)航定位的最基本的方法，其基本原理是：在某一瞬間利用GPS接收機(jī)同時(shí)測(cè)定至少四顆衛(wèi)星的偽距，根據(jù)已知的衛(wèi)星位置和偽距觀測(cè)值，采用距離交會(huì)法求出接收機(jī)的三維坐標(biāo)和時(shí)鐘改正數(shù)。它的優(yōu)點(diǎn)是速度快、無(wú)多值性問(wèn)題，利用增加觀測(cè)時(shí)間可以提高定位精度；缺點(diǎn)是測(cè)量定位精度低，但足以滿足部分用戶的需要。
?
?
|
|
?
圖5-5 用戶可見的衛(wèi)星分布
SatellitePosition =?????????????????????
??1.0e+004 *
??? 1.7746??? 1.7572??? 0.7365??? 0.0001
?? -1.2161??? 0.9732??? 2.1091??? 0.0001
??? 2.2883??? 0.2975?? -1.4132??? 0.0001
?? -2.3882?? -0.8836??? 0.7869???????? 0
?? -0.3681?? -2.5255??? 0.7012???????? 0
?? -0.1323??? 2.7059?? -0.0000??? 0.0001
?? -1.5424??? 0.4273??? 2.1704??? 0.0001
??? 1.4582??? 0.4926?? -2.1144???????? 0
?? -2.4090??? 0.6583??? 0.7365???????? 0
??? 1.0759??? 1.1599?? -2.1757?????? ??0
??? 1.3324?? -1.8178??? 1.4392??? 0.0001
?? -1.3225??? 1.7688?? -1.3835???????? 0
??? 1.2127?? -0.9774??? 2.1091??? 0.0001
?? -1.8766??? 0.3838?? -1.9058???????? 0
?? -2.4302?? -0.9295?? -0.7516???????? 0
??? 0.4443??? 1.8187?? -1.8463???????? 0
?? -0.6718??? 2.4562?? -0.7869???????? 0
??? 1.4000?? -1.3073??? 1.9058??? 0.0001
?? -1.4992??? 0.4223?? -2.1091???????? 0
?? -0.8706?? -2.0297?? -1.3835???????? 0
?? -1.7722?? -0.6021?? -1.9217???????? 0
?? -0.2481?? -1.5812?? -2.1704???????? 0
?? -0.5421??? 1.7915??? 1.9217??? 0.0001
??? 1.9376?? -1.5756?? -0.7365??? 0.0001
仿真程序四：用可見衛(wèi)星計(jì)算用戶位置
程序見附錄
多于四顆的衛(wèi)星用最小二乘法逼近計(jì)算。對(duì)于前面的假設(shè)用戶位置（6400，3352，5410）進(jìn)行計(jì)算，結(jié)果如下：
?
圖5-6 用戶位置的計(jì)算
由圖13可見，根據(jù)最小二乘法的原理，計(jì)算出來(lái)了用戶在一個(gè)時(shí)刻的位置坐標(biāo)，并把它表示在了這個(gè)天球坐標(biāo)系當(dāng)中。
calculaterecord =
? 1.0e+003 *
?????? 0? ???????0???????? 0
??? 7.3399??? 3.8857??? 6.1388
??? 6.3034??? 3.2976??? 5.3346
??? 6.3914??? 3.3470??? 5.4034
??? 6.4009??? 3.3525??? 5.4107
??? 6.4001??? 3.3520??? 5.4101
先后經(jīng)過(guò)了六次的迭代算法，吧用戶的計(jì)算位置一步一步的逼近了用戶的實(shí)際位置，根據(jù)部同的精度要求，我們運(yùn)算的量的大小也有不同。這點(diǎn)，可以根據(jù)程序當(dāng)中的參數(shù)Error的設(shè)定而有所出入。
?
圖5-7 用戶位置的計(jì)算
通過(guò)對(duì)用戶仿真計(jì)算的進(jìn)一步放大，可以很清楚地看到：在每一次迭代的過(guò)程當(dāng)中，我們都實(shí)時(shí)的把每一步迭代的值也表示在了同一個(gè)坐標(biāo)軸里面：目的就是很直觀的反應(yīng)出我們基于最小二乘法的偽距算法的主要思路。如下圖所示：
?
圖5-8 用戶位置的計(jì)算
?? 上圖15當(dāng)中的黑色箭頭從上到下一次表示的是第一次到第六次計(jì)算出來(lái)的用戶的位置，第七箭頭表示的用戶的實(shí)際位置（用白色的柱體表示）。
?
?
程序一（繪制衛(wèi)星的軌道平面）代碼
function drawsatelliteorbit
a=26560;% 衛(wèi)星軌道的長(zhǎng)半軸.
e=0.02;%e是橢圓的偏心率.
E=[0:0.1:2*pi];
x=a*(cos(E)-e);
y=a*sqrt((1-e^2))*sin(E);
z=0*E;
DtoR=2*pi/360;
A1=[32.8? 92.8? 152.8? 212.6? 272.8? 332.8];%衛(wèi)星星座數(shù)據(jù).
for k=1:6
??? A=A1(k)*DtoR;%升交點(diǎn)的經(jīng)度
??? B=55*DtoR;%軌道的傾角
??? C=pi/100;%近地點(diǎn)幅角
??? %總共有6個(gè)衛(wèi)星軌道平面
??? R1=[cos(A)? -sin(A)? 0;
??????? sin(A)? cos(A)?? 0;
??????? 0??????? 0?????? 1;];
??? R2=[1???????? 0?????? 0;
??????? 0??????? cos(B)? -sin(B);
??????? 0??????? sin(B)?? cos(B);];
??? R3=[cos(C)??? -sin(C)??? 0;
??????? sin(C)??? cos(C)???? 0;
??????? 0????????? 0??????? 1;];
??? L1=length(E);
??? R312=R1*R2*R3;
??? Ans=R312*[x;y;z;];
??? x1=Ans(1,:);
??? y1=Ans(2,:);
??? z1=Ans(3,:);
??? plot3c(x1,y1,z1,k);
??? boxplot3(0,0,0,200,200,200,7);
??? hold on;
? axis equal;
? axis off;
end
程序二（單顆衛(wèi)星不同時(shí)刻的動(dòng)態(tài)仿真）代碼
clear; clc;close all;
DtoR=2*pi/360;
jiaostep=360/24*DtoR;
?j=0;
??? a=26560;% 衛(wèi)星軌道的長(zhǎng)半軸.
e=0.02;%e是橢圓的偏心率.
E=[0:0.1:2*pi];
x=a*(cos(E)-e);
y=a*sqrt((1-e^2))*sin(E);
z=0*E;
??? drawearth(0);
hold on;
for time=1:12
??? axis on;
A1=32.8 ;%衛(wèi)星星座數(shù)據(jù).
??? A=A1*DtoR;%升交點(diǎn)的經(jīng)度
??? B=55*DtoR;%軌道的傾角
??? C=pi/100;%近地點(diǎn)幅角
??? R1=[cos(A)? -sin(A)? 0;
??????? sin(A)? cos(A)?? 0;
??????? 0??????? 0?????? 1;];
??? R2=[1???????? 0?????? 0;
??????? 0??????? cos(B)? -sin(B);
??????? 0??????? sin(B)?? cos(B);];
??? R3=[cos(C)??? -sin(C)??? 0;
??????? sin(C)??? cos(C)???? 0;
??????? 0????????? 0??????? 1;];
??? L1=length(E);
??? R312=R1*R2*R3;
??? Ans=R312*[x;y;z;];
??? x1=Ans(1,:);
??? y1=Ans(2,:);
??? z1=Ans(3,:);
??? plot3c(x1,y1,z1,2);
??? hold on;
? axis equal;
? axis on;
? grid on;
ctable=10;
??????? A=A1*DtoR;
??????? B=55*DtoR;
??????????? C=ctable*DtoR+time*2*pi/24;???????? %近地點(diǎn)幅角
????????????? x=a*(cos(C)-e);
y=a*sqrt((1-e^2))*sin(C);
z=0*C;
?R1=[cos(A)??? -sin(A)?? ????0;
????? sin(A)??? cos(A)?????? 0;
??????? 0?????????? 0?????? 1;];
??? R2=[1???????? 0??????? 0;
??????? 0??????? cos(B)? -sin(B);
??????? 0??????? sin(B)?? cos(B);];
??? R3=[cos(C)??? -sin(C)???? 0;
??????? sin(C)??? cos(C)????? 0;
??????? 0???????? ?0??????? 1;];
??? L1=length(E);
??? R312=R1*R2*R3;
??? Ans=R312*[x;y;z;];
??? x1=Ans(1,:);
??? y1=Ans(2,:);
??? z1=Ans(3,:);
??? drawsatellite(x1,y1,z1,6);??
??? M(time)=getframe;
??? M(time+1) = getframe;
?? end
axis on;???????????
程序三（衛(wèi)星在某個(gè)時(shí)刻的全軌道平面的分布和可見衛(wèi)星）代碼
clear;clc;close all;
a=26560;????????????????????????????????????????? % 衛(wèi)星軌道的長(zhǎng)半軸.
e=0.02;???
temp=0;%便于把衛(wèi)星的指標(biāo)點(diǎn)順序放到矩陣指標(biāo)當(dāng)中。
%e是橢圓的偏心率.
E=[0:0.1:2*pi];
x=a*(cos(E)-e);
y=a*sqrt((1-e^2))*sin(E);
z=0*E;
timenow=0;%單位是小時(shí)
global SatellitePosition
global ttum
SatellitePosition=ones(24,4);
figure(1);
drawearth(0);??
ttum=ones(24,4);
hold on;
DtoR=2*pi/360;
A1=[32.8? 92.8? 152.8? 212.6? 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;
??????? 120 140 240 350;];
??? simple=1;
???? for k=1:6
??????? A=A1(k)*DtoR;??????????????????????????????? %升交點(diǎn)經(jīng)度
???? ???B=55*DtoR;?????????????????????????????????? %軌道傾角
??????? for m=1:4
??????????? C=ctable(k,m)*DtoR+timenow*2*pi/24;???????? %近地點(diǎn)幅角
????????????? x=a*(cos(C)-e);
y=a*sqrt((1-e^2))*sin(C);
z=0*C;
?R1=[cos(A)??? -sin(A)?????? 0;
????? sin(A)??? cos(A)??? ???0;
??????? 0?????????? 0?????? 1;];
R2=[1???????? 0??????? 0;
??????? 0??????? cos(B)? -sin(B);
??????? 0??????? sin(B)?? cos(B);];
R3=[cos(C)??? -sin(C)???? 0;
??????? sin(C)??? cos(C)????? 0;
??????? 0????????? 0??????? 1;];
??? L1=length(E);
??? R312=R1*R2*R3;
??? Ans=R312*[x;y;z;];
??? x1=Ans(1,:);
??? y1=Ans(2,:);
??? z1=Ans(3,:);
?? drawsatellite(x1,y1,z1,k);?????????????????????????????
??? temp=temp+1;
??? SatellitePosition(temp,:)=[x1 y1 z1 1];
??? ttum(temp,:)=[x1 y1 z1 1];
??? hold on;
????? ??end
???? end
?%用余弦定理計(jì)算地心到用戶和用戶到衛(wèi)星證件的夾角要是兩者證件的夾角小于90°,就認(rèn)為改衛(wèi)星是不可見的.
???? earthcenterpos=[0 0 0];
???? userposition=[6400? 3352? 5410];?????????????? %假設(shè)一個(gè)用戶的位置
???? 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';
???????? jiajiao=acos((Dist2+Dist3-Dist1)/2/sqrt(Dist3)/sqrt(Dist2));
???????? if(jiajiao>=pi/2)
?????? ??????SatellitePosition(k,4)=0;
???????????? %ttum(k,:)=SatellitePosition(k,1:3);
???????? end
???? end
?? %? ttum????
???? figure(2)
??? drawearth(0);????????????????????????????????????
???? hold on;
???? drawsatelliteorbit;?
???????? j=1;
????? for k=1:6
????????? for m=1:4
???????? if(SatellitePosition(j,4)==1)
???????????? tempx=SatellitePosition(j,1);
???????????? tempy=SatellitePosition(j,2);
???????????? tempz=SatellitePosition(j,3);
???????????? drawsatellite(tempx,tempy,tempz,k);???????????
? ???????end
????????? j=j+1;
????????? end
????? end
???? tempx=userposition(1);
???? tempy=userposition(2);
???? tempz=userposition(3);
???? cube=100;
???? boxplot3(tempx,tempy,tempz,cube,cube,cube,0); %標(biāo)記出用戶的位置
程序四（用可見衛(wèi)星計(jì)算用戶的位置）代碼（mian3）
%用余弦定理計(jì)算地心到用戶和用戶到衛(wèi)星證件的夾角要是兩者證件的夾角小于90°,就認(rèn)為該衛(wèi)星是不可見的.
???? earthcenterpos=[0 0 0];
???? userposition=[6400? 3352? 5410];?????????????? %假設(shè)一個(gè)用戶的位置
???? for k=1:24
???? SatellitePosition(k,:)=ttum(k,:);
???? end
???? 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';
???????? jiajiao=acos((Dist3+Dist2-Dist1)/2/sqrt(Dist3)/sqrt(Dist2));
? ???????if(jiajiao>=pi/2)
???????????? SatellitePosition(k,4)=0;
???????? end
???? end
???? figure(1)
??? drawearth(0);??????????????????????????????????
???? hold on;
???? drawsatelliteorbit;?
????? tempx=userposition(1);
???? tempy=userposition(2);
??? ?tempz=userposition(3);
???????? cube=55;
???????? boxplot3(tempx,tempy,tempz,cube,cube,cube,7);????????????
???????? hold on; %先畫出用戶的位置，以便和后面的做比較
????? j=1;
????? for k=1:6
????????? for m=1:4
???????? if(SatellitePosition(j,4)==1)
???????????? tempx=SatellitePosition(j,1);
???????????? tempy=SatellitePosition(j,2);
???????????? tempz=SatellitePosition(j,3);
???????????? drawsatellite(tempx,tempy,tempz,k);????????????
???????? end
????????? j=j+1;
????????? end
????? end
????? tempx=userposition(1);
??? ?tempy=userposition(2);
???? tempz=userposition(3);
???? hold on;
???? SatellitePosition=[SatellitePosition;userposition 0];????
?? %用變量caluserpos來(lái)記錄迭代計(jì)算的結(jié)果.
???? caluserpos=calculateuserposition(SatellitePosition);
???? [m,n]=size(caluserpos);
???? for(k=3:m);
???????? tempx=caluserpos(k,1);
???????? tempy=caluserpos(k,2);
???????? tempz=caluserpos(k,3);
???????? cube=50;
???????? boxplot3(tempx,tempy,tempz,cube,cube,cube,0);????????????
???????? hold on;
???? end
以下是一個(gè)和上面程序配套使用的功能函數(shù)，把它命名為calculateuserposition.m。并且要和上面的程序放在一個(gè)文件夾當(dāng)中。
function caluserposition=calculateuserposition(SatellitePosition)
global SatellitePosition?????????????
?%單獨(dú)運(yùn)行此程序時(shí)候，為調(diào)用前面程序而得的satelliteposotion結(jié)果，應(yīng)該聲明satelliteposotion為全局變量，即應(yīng)該在satelliteposition前面加上global
r1=6400;?????????? ?????????????????????????????????????????%地球半徑
c=300000;
deltat=0.00001;
satelliteposnew=ones(1,3);
vissatnum=0;
calculateok=1;
for k=1:24
??? if(SatellitePosition(k,4)==1)
??????? vissatnum=vissatnum+1;
??????? satelliteposnew=[satelliteposnew;SatellitePosition(k,1:3)];
??? end%if
end%for
satelliteposnew(1,:)=[];
if(vissatnum<4)
??? calculateok=0;
??? caluserposition=[0 0 0 0];
??? return
end
prange=ones(1,vissatnum);
userpos=SatellitePosition(25,1:3);
for n=1:vissatnum??? prange(1,n)=sqrt((satelliteposnew(n,:)-userpos)*(satelliteposnew(n,:)-userpos)')+c*deltat;
end
calculaterecord=[0??? 0?? 0];
xyz0=[0 0 0];
deltat0=0;
wxyz=satelliteposnew;
Error=1000;
computertime=0;
while((Error>0.00001)&(computertime<1000))
??? computertime=computertime+1;
??? r=ones(1,vissatnum);
??? for n=1:vissatnum
??????? r(1,n)=sqrt((wxyz(n,:)-xyz0)*(wxyz(n,:)-xyz0)')+deltat0*c;
??? end
??? deltap=r-prange;
??? A=ones(vissatnum,3);
??? for n=1:vissatnum
??? A(n,:)=(wxyz(n,:)-xyz0)./r(1,n);
??? end%for
??? H=[A ones(vissatnum,1)];
??? deltax=inv(H'*H)*H'*deltap';
??? tempdeltax=deltax(1:3)
??? Error=max(abs(tempdeltax));
??? xyz0=xyz0+deltax(1:3,:)';
??? if(computertime<1000)
??????? calculaterecord=[calculaterecord;xyz0];
??? end
??? deltat0=deltax(4,1)/(-c);
end%while
calculaterecord;
if(computertime==1000)
??? caluserposition=[99999999999 99999999999 999999999999];
else
??? caluserposition=[xyz0;calculaterecord];
end
總結(jié)
以上是生活随笔為你收集整理的GPS卫星运动及定位matlab仿真的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇： 主流的6个Go语言Web框架
- 下一篇： matlab中omg什么意思,英文聊天中