getCOT 阶比分析文件。特别适合对数据进行处理

阶比分析文件。特别适合对数据进行处理，特别别是信号处理。

%%－－－－－－阶比跟踪算法－－－－－
function [Tnxtn] = getCOT（xtfsDmaxpforderwu）
    %xtn为输出：等角度采样的信号序列
    %输入：x为等时间间隔采样信号序列，t为时间，fs采样频率，Dmax为最大阶次，pf为频率曲线序列order:拟合频率曲线的阶次，wu：舍弃的点数
    
    t = t - min（t）;
     
    dw = pi/Dmax; %重采样角度间隔
    dt = 1/fs ;%采样时间间隔
    
    a = polyfit（tpforder）; %3阶拟合：ft = a（1）*t.^3 + a（2）*t.^2 + a（3）*t+a（4）;
    ft = polyval（at）;  %得到拟合后的频率曲线序列
    Na = length（a）;
    for j = 1:Na
        a（j） = a（j）/（Na-j+1）;
    end
    
    lenXtn = fix（2*pi*sum（ft*dt）/dw）; %计算重采样后的数据长度
    lenXtn = lenXtn -wu;
    
    Tn = zeros（1lenXtn）; % 计算键相时标 
    for ii = 1 : lenXtn      %   求解方程
        temp = ii/（2*Dmax）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1455  2011-07-01 11:30  getCOT.m

----------- ---------  ---------- -----  ----

                 1455                    1