时域平均程序

% 根据刘红星文章，实现无时标信号的直接时域平均
function out = Liu_TDA（xfsf0）

% 获取x的数据长度
len = length（x）;

x = x（:）‘;

% 粗略计算一个周期有多少个数据点
N = round（fs/f0）;

% 精确计算一个周期有多少个数据点
Nfr = fs/f0;

% 计算原数据点有包含多少个整周期
M = floor（len/Nfr）;

% 结果存储数组
out = zeros（1N）;
 
% 计算每个周期的起始点
ind = zeros（1M）;

for i=0:M-1
    ind（i+1） = roun