TSA 时域同步平均法

时域同步平均法，在降噪上有很明显的帮助TSA

clcclear allclose all
%% 时域同步平均仿真 

% 产生一个无噪的正弦信号
w=10;
t=0:（2*pi/w/500）:（300*2*pi/w）;
x=sin（w*t）;
% 产生白噪声
n=length（t）;
ss=888;
randn（‘seed‘ss）;
z=randn（1n）;
z1=z-mean（z）;
z2=z1/max（abs（z1））;
z2=z2*5;
% 含噪的正弦信号 
x=x+z2;
% 绘制含噪的正弦信号
figure
plot（x）;title（‘仿真含噪信号‘）
N = 2048;
n1 = 0:N-1;
fs = 2500/pi;
f=n1*fs/N;
x1=fft（xN）;
mag1=abs（x1）*2/N;
figure
plot（f（1:N/2）mag1（1:N/2））;
title（‘仿真信号频谱‘）

%% Averaging in time domain

M=10000;  % 数据长度
xx=zeros（1M）; % 用于存储时域平均后的数据
m1=1;m2=M;

nnnn=fix（n/M）; % 确定分段数

% 将数据分成nnnn段，

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

     文件       1249  2014-11-26 14:58  TSA.m

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

                 1249                    1