main.m

大学时候跟导师做的一个暂态电磁场数据分析的程序

filename='A019.FLD';
[wav,length]=waveread(filename);
lev=5;
[c,l] = wavedec(wav,lev,'sym8');
 % threshold the decomposition structure [c,l].
xd = wden(c,l,'minimaxi','s','sln',lev,'sym8');
figure(1);plot(xd);xlabel('ns');ylabel('kv');grid on;
[c1,l1] = wavedec(xd,lev,'sym8');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
segoffset=1;    % 三个参数分别为：分段起始点-－结束点-－信号限
len=length;     % 
signalval=0.02;   % 默认值结束点等于数据总长，起始为总开始点
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
                   % 可选。
pulse=waveout(xd,segoffset,len,signalval);
wavsave(filename,signalval,pulse);