周跳探测程序Matlab

本资源使用matlab读取卫星星历文件进行周跳探测，计算出周跳位置和大概周跳的周数。具有很强的实用性，读者也可以在此基础上进行改进，写出更好的程序来。

function zttc
% 探测周跳
% 调用格式：zttc
%
fid = fopen（‘zt.txt‘‘r‘）;
for i=1:97
    status = fseek（fid720）;     
    Oo（i） = fscanf（fid‘%f‘1）;    
    status = fseek（fid1980）;    
end
fclose（fid）;
%fid = fopen（‘周跳数据.dat‘‘w‘）;
%count = fprintf（fid‘%.5f  ‘O）
%closestatus=fclose（fid）
for i=1:48
    O（i） = Oo（i）;
end
for i=49:97
    O（i） = Oo（i）+100;    % 加周跳  （经测试，此组数据最小只能准确探测到 + -3周）
end
to = 1:97;
t1=1:96; t2=1:95; t3=1:94; t4=1:93;
%*********************************
for i=1:96
    A（i） = O（i+1）-O（i）;  % 一次差
end
for i=1:95
    B（i） = A（i+1）-A（i）;  % 二次差
end
for i=1:94
    C（i） = B（i+1）-B（i）;  % 三次差
end
for i=1:93
    D（i） = C（i+1）-C（i）;  % 四次差
end
%*********************************
flag=0;
for i=1:92                % 正序查找
    if abs（D（i））<2 & abs（D（i+

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       27742  2010-09-28 19:53  周跳探测程序（Matlab）\zt.txt
     文件        3689  2010-09-28 22:46  周跳探测程序（Matlab）\zttc.m
     文件       68286  2010-09-28 22:49  周跳探测程序（Matlab）\周跳运行结果.JPG
     文件         383  2014-04-01 08:50  周跳探测程序（Matlab）\程序说明.txt
     目录           0  2014-04-01 08:48  周跳探测程序（Matlab）\