function Ipp = find_ionodelay( latpp,lonpp,gpstime )
%根据Klobuchar模型计算垂直电离层延迟
%clear all
%close all
%clc
%latpp=39.305*pi/180;
%lonpp=109.87*pi/180;
%gpstime=5000;
ALPHA=[1.583e-008 1.49e-008 -1.192e-007 -1.192e-007];
BETA=[106500 65540 -196600 -196600];
%F = 1 + 16*(0.53 - E)^3;
phim = latpp/pi + 0.064*cos((lonpp/pi-1.616)*pi);%因为地磁纬度是半圆的所以要除以pi
%phim = latpp + 0.064*cos(lonpp-1.616);
%t = 4.32e4*lonpp + gpstime;
t = 4.32e4*lonpp/pi + gpstime;
while ( (t < 0) || (t >= 86400) )
if t >= 86400,
t = t - 86400;
end
if t < 0
t = t + 86400;
end
end
A1=5e-9;
A2 = ALPHA(1) + ALPHA(2)*phim + ALPHA(3)*phim^2 + ALPHA(4)*phim^3;
if A2 < 0
A2 = 0;
end
A3=50400;
A4 = BETA(1) + BETA(2)*phim + BETA(3)*phim^2 + BETA(4)*phim^3;
if A4 < 72000
A4 = 72000;
end
x = 2*pi*(t - A3)/A4;
if abs(x) < 1.57,
Tiono = A1 + A2*cos(x);
%Tiono = A1 + A2*(1 - x^2/2 + x^4/24);
else
Tiono = A1;
end
Ipp=Tiono*299792458;
matlab_计算GPS导航系统在Klobuchar模型下的电离层延迟
版权申诉
5星 · 超过95%的资源 172 浏览量
2022-05-06
23:40:36
上传
评论 2
收藏 915B ZIP 举报
wouderw
- 粉丝: 286
- 资源: 2961
- 1
- 2
前往页