我编写了一段数据处理程序,如下:
clc;
clear all;
close all;
load F:\科研立项\Data\kuai2_e14.mat
s=data1(1*10^5:4*10^5)-mean(data1(1*10^5:4*10^5));
clear data1
n=length(s); %样本长度
fs=25000; %采样频率
s=s/max(abs(s)); %归一化
pp=downsample(s,5); %降采样
nn=length(s);
fss=fs/5;
nn=fss; %fix(length(pp)/5);
s=pp(nn+1:2*nn)';
max_tau=500;
%[ a,R ] = ComputeRalpha(s',20);
%[ a,R ] = Blanchit(alpha,R);
%I = floor(0.50*length(R)):length(R);
T=1000;
alpha=1/T;%1./(1:T);%Ta(floor(0.50*length(R))+1:length(R));
%clear a;
%nm=length(alpha);
%rb=zeros(nm,max_tau*2+1,T+1);
%for k&#