求用matlab编写turboedit方法探测和修复gps周跳的程序

function [t_cy t_ou L_r]=turboedit(P,L_l)

% 本程序通过turboedit方法探测周跳与野值，并修复。

% 输入稿孙并time历元信息,P伪距观测值，L_l相位观测值

% 输出time_cy为发生周跳的历元，time_ou为存在野值的历元。

% L_2为去了周跳后凯清的伪距观测值；

%全局变量

c=299792458

f=[157542000 122760000]

t_cy=[]t_ou=[]L_r=[]

%MW组合探测周跳

%载波相位转化为距离

for i=1:size(L_l,2)

L(:,i)=(L_l(:,i)*c)/f(i)

end

L_mw=(f(1)*L(:,1)-f(2)*L(:,2))/(f(1)-f(2))-(f(1)*P(:,1)+f(2)*P(:,2))/(f(1)+f(2))

b_w=(L_mw*(f(1)-f(2)))/c

b_wp(1)=b_w(1)var(1)=0j=1k=1l=2

for i=2:length(b_w)-1

var(i)=sqrt((((l-1)/l)*(var(i-1)^2))+((b_w(i)-b_wp(i-1))^2)/l)

%判断周跳条件

if abs(b_wp(i-1)-b_w(i))>=4*var(i) &&abs(b_w(i+1)-b_w(i))<=1

b_wp(i)=b_w(i)var(i)=0l=2

t_cy(j,1)=(i-1)*30 t_cy(j,2)=i

j=j+1

else

b_wp(i)=(l-1)*b_wp(i-1)/l+b_w(i)/l

end

%判断野值条件

if abs(b_wp(i-1)-b_w(i))>=4*var(i) &&abs(b_w(i+1)-b_w(i))>1

t_ou(k,1)=(i-1)*30 t_ou(k,2)=i

k=k+1

end

l=l+1

end

%GF组合探测周跳

L_I = L(:,1)-L(:,2)

P_I = P(:,2)-P(:,1)

c_fI = c/f(2)-c/f(1)

yo0 = polyfit(1:size(P,1),P_I',6)

Q_I = polyval(yo0,1:size(P,1))

d_LI=(L_I-Q_I')'

for i=2:length(L_I)-1

if (L_I(i)-Q_I(i))-(L_I(i-1)-Q_I(i-1))>6*c_fI &&(L_I(i+1)-Q_I(i+1))-(L_I(i)-Q_I(i))<c_fI

t_cy(j,1)=(i-1)*30 t_cy(j,2)=i

j=j+1

end

end

t_cy=sort(t_cy,1)

% 修复周跳

%注！！！这里选择30个点平均，键迹根据实际情况调整

if ~isempty(t_cy) &&t_cy(1,2)>=32 &&t_cy(end,2)<=length(b_w)-30

for i=1:size(t_cy,1)

d_bw=(sum(b_w(t_cy(i,2):t_cy(i,2)+30))-sum(b_w(t_cy(i,2)-31:t_cy(i,2)-1)))/30

yyy_y=polyfit(t_cy(i,2)-31:t_cy(i,2)-1,d_LI(t_cy(i,2)-31:t_cy(i,2)-1),2)

dd_LI=d_LI(t_cy(i,2))-polyval(yyy_y,t_cy(i,2))

b2=(c/f(1)*d_bw-dd_LI)/c_fI

b1=d_bw+b2

b(i,1)=round(b1)

b(i,2)=round(b2)

L_r(:,1)=[L_l(1:(t_cy(i,2)-1),1)L_l(t_cy(i,2):end,1)-b1]

L_r(:,2)=[L_l(1:(t_cy(i,2)-1),2)L_l(t_cy(i,2):end,2)-b2]

end

else

L_r=L_l

end

t_cy=sort(t_cy,1)

周跳探测的电离层残差法是一种用于评估电离层对GPS信号的影响的方法，它通过计算GPS卫星发射的信号与地面接收机测量到的信号之间的时间差来评估电离郑老层对GPS信号的影响。这个时间差被称为“周跳”，它是由裂悉于地球上肆丛乎不同位置处的电离层厚度不同而引起的。

---